Connecting species’ geographical distributions to environmental variables: range maps versus observed points of occurrence

Connecting the geographical occurrence of a species with underlying environmental variables is fundamental for many analyses of life history evolution and for modeling species distributions for both basic and practical ends. However, raw distributional information comes principally in two forms: points of occurrence (specific geographical coordinates where a species has been observed), and expert-prepared range maps. Each form has potential short-comings: range maps tend to overestimate the true occurrence of a species, whereas occurrence points (because of their frequent non-random spatial distribution) tend to underestimate it. Whereas previous comparisons of the two forms have focused on how they may differ when estimating species richness, less attention has been paid to the extent to which the two forms actually differ in their representation of a species’ environmental associations. We assess such differences using the globally distributed avian order Galliformes (294 species). For each species we overlaid range maps obtained from IUCN and point-of-occurrence data obtained from GBIF on global maps of four climate variables and elevation. Over all species, the median difference in distribution centroids was 234 km, and median values of all five environmental variables were highly correlated, although there were a few species outliers for each variable. We also acquired species’ elevational distribution mid-points (mid-point between minimum and maximum elevational extent) from the literature; median elevations from point occurrences and ranges were consistently lower (median −420 m) than mid-points. We concluded that in most cases occurrence points were likely to produce better estimates of underlying environmental variables than range maps, although differences were often slight. We also concluded that elevational range mid-points were biased high, and that elevation distributions based on either points or range maps provided better estimates.     

Atlas versus range maps: robustness of chorological relationships to distribution data types in European mammals

Aim Chorological relationships describe the patterns of distributional overlap among species. In addition to revealing biogeographical structure, the resulting clusters of species with similar geographical distributions can serve as natural units in conservation planning. Here, we assess the extent to which temporal, methodological and taxonomical differences in the source of species’ distribution data can affect the relationships that are found. Location Western Europe. Methods We used two data sets – the Atlas of European mammals and polygon range maps from the IUCN Global Mammal Assessment – both as presence–absence data for UTM 50 km × 50 km squares. We performed pairwise comparisons among 156 species for each data set to build matrices of the similarity in distribution across species, using both Jaccard’s and Baroni‐Urbani & Buser’s indices. We then compared these similarity matrices (chorological relationships), as well as the species richness and occurrence patterns from the two data sets. Results As expected, range maps increased both the mean prevalence per species and mean species richness per grid cell in comparison to atlas data, reflecting the general view that these data types respectively over‐ and underestimate species occurrence. However, species richness and occurrence patterns in atlas and range map data were positively associated and, most importantly, the chorological relationships underlying the two data sets were highly similar. Main conclusions Despite many methodological, temporal and taxonomical differences between atlas data and range maps, the chorological relationships encountered between species were similar for both data sets. Chorological analyses can thus be robust to the data source used and provide a solid basis for analytical biogeographical studies, even over broad spatial scales.     

Geographic range estimates and environmental requirements for the harpy eagle derived from spatial models of current and past distribution

Understanding species–environment relationships is key to defining the spatial structure of species distributions and develop effective conservation plans. However, for many species, this baseline information does not exist. With reliable presence data, spatial models that predict geographic ranges and identify environmental processes regulating distribution are a cost‐effective and rapid method to achieve this. Yet these spatial models are lacking for many rare and threatened species, particularly in tropical regions. The harpy eagle (Harpia harpyja) is a Neotropical forest raptor of conservation concern with a continental distribution across lowland tropical forests in Central and South America. Currently, the harpy eagle faces threats from habitat loss and persecution and is categorized as Near‐Threatened by the International Union for the Conservation of Nature (IUCN). Within a point process modeling (PPM) framework, we use presence‐only occurrences with climatic and topographical predictors to estimate current and past distributions and define environmental requirements using Ecological Niche Factor Analysis. The current PPM prediction had high calibration accuracy (Continuous Boyce Index = 0.838) and was robust to null expectations (pROC ratio = 1.407). Three predictors contributed 96% to the PPM prediction, with Climatic Moisture Index the most important (72.1%), followed by minimum temperature of the warmest month (15.6%) and Terrain Roughness Index (8.3%). Assessing distribution in environmental space confirmed the same predictors explaining distribution, along with precipitation in the wettest month. Our reclassified binary model estimated a current range size 11% smaller than the current IUCN range polygon. Paleoclimatic projections combined with the current model predicted stable climatic refugia in the central Amazon, Guyana, eastern Colombia, and Panama. We propose a data‐driven geographic range to complement the current IUCN range estimate and that despite its continental distribution, this tropical forest raptor is highly specialized to specific environmental requirements.     

Reconstructing range dynamics and range fragmentation of European bison for the last 8000 years

Aim Understanding what constituted species’ ranges prior to large‐scale human influence, and how past climate and land use change have affected range dynamics, provides conservation planners with important insights into how species may respond to future environmental change. Our aim here was to reconstruct the Holocene range of European bison (Bison bonasus) by combining a time‐calibrated species distribution models (SDM) with a dynamic vegetation model. Location Europe. Method We used European bison occurrences from the Holocene in a maximum entropy model to assess bison range dynamics during the last 8000 years. As predictors, we used bioclimatic variables and vegetation reconstructions from the generalized dynamic vegetation model LPJ‐GUESS. We compared our range maps with maps of farmland and human population expansion to identify the main species range constraints. Results The Holocene distribution of European bison was mainly determined by vegetation patterns, with bison thriving in both broadleaved and coniferous forests, as well as by mean winter temperature. The heartland of European bison was in Central and Eastern Europe, whereas suitable habitat in Western Europe was scarce. While environmentally suitable regions were overall stable, the expansion of settlements and farming severely diminished available habitat. Main conclusions European bison habitat preferences may be wider than previously assumed, and our results suggest that the species had a more eastern and northern distribution than previously reported. Vegetation and climate transformation during the Holocene did not affect the bison’s range substantially. Conversely, human population growth and the spread of farming resulted in drastic bison habitat loss and fragmentation, likely reaching a tipping point during the last 1000 years. Combining SDM and dynamic vegetation models can improve range reconstructions and projections, and thus help to identify resilient conservation strategies for endangered species.     

Determinants of palm species distributions across Africa: the relative roles of climate,non‐climatic environmental factors,and spatial constraints

Most of the Earth's biodiversity resides in the tropics. However, a comprehensive understanding of which factors control range limits of tropical species is still lacking. Climate is often thought to be the predominant range‐determining mechanism at large spatial scales. Alternatively, species’ ranges may be controlled by soil or other environmental factors, or by non‐environmental factors such as biotic interactions, dispersal barriers, intrinsic population dynamics, or time‐limited expansion from place of origin or past refugia. How species ranges are controlled is of key importance for predicting their responses to future global change. Here, we use a novel implementation of species distribution modelling (SDM) to assess the degree to which African continental‐scale species distributions in a keystone tropical group, the palms (Arecaceae), are controlled by climate, non‐climatic environmental factors, or non‐environmental spatial constraints. A comprehensive data set on African palm species occurrences was assembled and analysed using the SDM algorithm Maxent in combination with climatic and non‐climatic environmental predictors (habitat, human impact), as well as spatial eigenvector mapping (spatial filters). The best performing models always included spatial filters, suggesting that palm species distributions are always to some extent limited by non‐environmental constraints. Models which included climate provided significantly better predictions than models that included only non‐climatic environmental predictors, the latter having no discernible effect beyond the climatic control. Hence, at the continental scale, climate constitutes the only strong environmental control of palm species distributions in Africa. With regard to the most important climatic predictors of African palm distributions, water‐related factors were most important for 25 of the 29 species analysed. The strong response of palm distributions to climate in combination with the importance of non‐environmental spatial constraints suggests that African palms will be sensitive to future climate changes, but that their ability to track suitable climatic conditions will be spatially constrained.     

Model-based uncertainty in species range prediction

Aim Many attempts to predict the potential range of species rely on environmental niche (or ‘bioclimate envelope’) modelling, yet the effects of using different niche‐based methodologies require further investigation. Here we investigate the impact that the choice of model can have on predictions, identify key reasons why model output may differ and discuss the implications that model uncertainty has for policy‐guiding applications. Location The Western Cape of South Africa. Methods We applied nine of the most widely used modelling techniques to model potential distributions under current and predicted future climate for four species (including two subspecies) of Proteaceae. Each model was built using an identical set of five input variables and distribution data for 3996 sampled sites. We compare model predictions by testing agreement between observed and simulated distributions for the present day (using the area under the receiver operating characteristic curve (AUC) and kappa statistics) and by assessing consistency in predictions of range size changes under future climate (using cluster analysis). Results Our analyses show significant differences between predictions from different models, with predicted changes in range size by 2030 differing in both magnitude and direction (e.g. from 92% loss to 322% gain). We explain differences with reference to two characteristics of the modelling techniques: data input requirements (presence/absence vs. presence‐only approaches) and assumptions made by each algorithm when extrapolating beyond the range of data used to build the model. The effects of these factors should be carefully considered when using this modelling approach to predict species ranges. Main conclusions We highlight an important source of uncertainty in assessments of the impacts of climate change on biodiversity and emphasize that model predictions should be interpreted in policy‐guiding applications along with a full appreciation of uncertainty.     

Species distribution modelling as a macroecological tool: a case study using New World amphibians

Although species distribution modelling (SDM) is widely accepted among the scientific community and is increasingly used in ecology, conservation biology and biogeography, methodological limitations generate potential problems for its application in macroecology. Using amphibian species richness in North and South America, we compare species richness patterns derived from SDM maps and ‘expert’ maps to evaluate if: 1) richness patterns derived from SDM are biased toward climate‐based explanations for diversity when compared to expert maps, since SDM methods are typically based on climatic variables; and 2) SDM is a reliable tool for generating richness maps in hyperrich regions where point occurrence data are limited for many species. We found that although three widely used SDM methods overestimated amphibian species richness in grid cells when compared to expert richness maps in both North and South America due to systematic overestimation of range sizes, diversity gradients were reasonably robust at broad scales. Further, climatic variables statistically explained patterns of richness at similar levels among the different richness sources, although climatic relationships were stronger in the much better known North America than in South America. We conclude that in the face of the high deforestation rates coupled with incomplete data on species distributions, especially in the tropics, SDM represents a useful macroecological tool for investigating broad‐scale richness patterns and the dynamics between species richness and climate.     

Climate limitation at the cold edge: contrasting perspectives from species distribution modelling and a transplant experiment

The role of climate in determining range margins is often studied using species distribution models (SDMs), which are easily applied but have well-known limitations, e.g. due to their correlative nature and colonization and extinction time lags. Transplant experiments can give more direct information on environmental effects, but often cover small spatial and temporal scales. We simultaneously applied a SDM using high-resolution spatial predictors and an integral projection (demographic) model based on a transplant experiment at 58 sites to examine the effects of microclimate, light and soil conditions on the distribution and performance of a forest herb, Lathyrus vernus, at its cold range margin in central Sweden. In the SDM, occurrences were strongly associated with warmer climates. In contrast, only weak effects of climate were detected in the transplant experiment, whereas effects of soil conditions and light dominated. The higher contribution of climate in the SDM is likely a result from its correlation with soil quality, forest type and potentially historic land use, which were unaccounted for in the model. Predicted habitat suitability and population growth rate, yielded by the two approaches, were not correlated across the transplant sites. We argue that the ranking of site habitat suitability is probably more reliable in the transplant experiment than in the SDM because predictors in the former better describe understory conditions, but that ranking might vary among years, e.g. due to differences in climate. Our results suggest that L. vernus is limited by soil and light rather than directly by climate at its northern range edge, where conifers dominate forests and create suboptimal conditions of soil and canopy-penetrating light. A general implication of our study is that to better understand how climate change influences range dynamics, we should not only strive to improve existing approaches but also to use multiple approaches in concert.     

Performance tradeoffs in target‐group bias correction for species distribution models

Species distribution models (SDMs) are often calibrated using presence‐only datasets plagued with environmental sampling bias, which leads to a decrease of model accuracy. In order to compensate for this bias, it has been suggested that background data (or pseudoabsences) should represent the area that has been sampled. However, spatially‐explicit knowledge of sampling effort is rarely available. In multi‐species studies, sampling effort has been inferred following the target‐group (TG) approach, where aggregated occurrence of TG species informs the selection of background data. However, little is known about the species‐ specific response to this type of bias correction. The present study aims at evaluating the impacts of sampling bias and bias correction on SDM performance. To this end, we designed a realistic system of sampling bias and virtual species based on 92 terrestrial mammal species occurring in the Mediterranean basin. We manipulated presence and background data selection to calibrate four SDM types. Unbiased (unbiased presence data) and biased (biased presence data) SDMs were calibrated using randomly distributed background data. We used real and TG‐estimated sampling efforts in background selection to correct for sampling bias in presence data. Overall, environmental sampling bias had a deleterious effect on SDM performance. In addition, bias correction improved model accuracy, and especially when based on spatially‐explicit knowledge of sampling effort. However, our results highlight important species‐specific variations in susceptibility to sampling bias, which were largely explained by range size: widely‐distributed species were most vulnerable to sampling bias and bias correction was even detrimental for narrow‐ranging species. Furthermore, spatial discrepancies in SDM predictions suggest that bias correction effectively replaces an underestimation bias with an overestimation bias, particularly in areas of low sampling intensity. Thus, our results call for a better estimation of sampling effort in multispecies system, and cautions the uninformed and automatic application of TG bias correction.     

Consensus Forecasting of Species Distributions: The Effects of Niche Model Performance and Niche Properties

Ensemble forecasting is advocated as a way of reducing uncertainty in species distribution modeling (SDM). This is because it is expected to balance accuracy and robustness of SDM models. However, there are little available data regarding the spatial similarity of the combined distribution maps generated by different consensus approaches. Here, using eight niche-based models, nine split-sample calibration bouts (or nine random model-training subsets), and nine climate change scenarios, the distributions of 32 forest tree species in China were simulated under current and future climate conditions. The forecasting ensembles were combined to determine final consensual prediction maps for target species using three simple consensus approaches (average, frequency, and median [PCA]). Species’ geographic ranges changed (area change and shifting distance) in response to climate change, but the three consensual projections did not differ significantly with respect to how much or in which direction, but they did differ with respect to the spatial similarity of the three consensual predictions. Incongruent areas were observed primarily at the edges of species’ ranges. Multiple stepwise regression models showed the three factors (niche marginality and specialization, and niche model accuracy) to be related to the observed variations in consensual prediction maps among consensus approaches. Spatial correspondence among prediction maps was the highest when niche model accuracy was high and marginality and specialization were low. The difference in spatial predictions suggested that more attention should be paid to the range of spatial uncertainty before any decisions regarding specialist species can be made based on map outputs. The niche properties and single-model predictive performance provide promising insights that may further understanding of uncertainties in SDM.     

Modelling the occurrence of threatened plant species in taiga landscapes: methodological and ecological perspectives

Aim Understanding the spatial patterns of species distribution and predicting the occurrence of high biological diversity and rare species are central themes in biogeography and environmental conservation. The aim of this study was to model and scrutinize the relative contributions of climate, topography, geology and land‐cover factors to the distributions of threatened vascular plant species in taiga landscapes in northern Finland. Location North‐east Finland, northern Europe. Methods The study was performed using a data set of 28 plant species and environmental variables at a 25‐ha resolution. Four different stepwise selection algorithms [Akaike information criterion (AIC), Bayesian information criterion (BIC), adaptive backfitting, cross selection] with generalized additive models (GAMs) were fitted to identify the main environmental correlates for species occurrences. The accuracies of the distribution models were evaluated using fourfold cross‐validation based on the area under the curve (AUC) derived from receiver operating characteristic plots. The GAMs were tentatively extrapolated to the whole study area and species occurrence probability maps were produced using GIS techniques. The effect of spatial autocorrelation on the modelling results was also tested by including autocovariate terms in the GAMs. Results According to the AUC values, the model performance varied from fair to excellent. The AIC algorithm provided the highest mean performance (mean AUC = 0.889), whereas the lowest mean AUC (0.851) was obtained from BIC. Most of the variation in the distribution of threatened plant species was related to growing degree days, temperature of the coldest month, water balance, cover of mire and mean elevation. In general, climate was the most powerful explanatory variable group, followed by land cover, topography and geology. Inclusion of the autocovariate only slightly improved the performance of the models and had a minor effect on the importance of the environmental variables. Main conclusions The results confirm that the landscape‐scale distribution patterns of plant species can be modelled well on the basis of environmental parameters. A spatial grid system with several environmental variables derived from remote sensing and GIS data was found to produce useful data sets, which can be employed when predicting species distribution patterns over extensive areas. Landscape‐scale maps showing the predicted occurrences of individual or multiple threatened plant species may provide a useful basis for focusing field surveys and allocating conservation efforts.     

Disentangling biotic interactions,environmental filters,and dispersal limitation as drivers of species co‐occurrence

A key focus in ecology is to search for community assembly rules. Here we compare two community modelling frameworks that integrate a combination of environmental and spatial data to identify positive and negative species associations from presence–absence matrices, and incorporate an additional comparison using joint species distribution models (JSDM). The frameworks use a dichotomous logic tree that distinguishes dispersal limitation, environmental requirements, and interspecific interactions as causes of segregated or aggregated species pairs. The first framework is based on a classical null model analysis complemented by tests of spatial arrangement and environmental characteristics of the sites occupied by the members of each species pair (Classic framework). The second framework, (SDM framework) implemented here for the first time, builds on the application of environmentally‐constrained null models (or JSDMs) to partial out the influence of the environment, and includes an analysis of the geographical configuration of species ranges to account for dispersal effects. We applied these approaches to examine plot‐level species co‐occurrence in plant communities sampled along a wide elevation gradient in the Swiss Alps. According to the frameworks, the majority of species pairs were randomly associated, and most of the non‐random positive and negative species associations could be attributed to environmental filtering and/or dispersal limitation. These patterns were partly detected also with JSDM. Biotic interactions were detected more frequently in the SDM framework, and by JSDM, than in the Classic framework. All approaches detected species aggregation more often than segregation, perhaps reflecting the important role of facilitation in stressful high‐elevation environments. Differences between the frameworks may reflect the explicit incorporation of elevational segregation in the SDM framework and the sensitivity of JSDM to the environmental data. Nevertheless, all methods have the potential to reveal general patterns of species co‐occurrence for different taxa, spatial scales, and environmental conditions.     

Beyond climate control on species range: The importance of soil data to predict distribution of Amazonian plant species

Aim

To evaluate the relative importance of climatic versus soil data when predicting species distributions for Amazonian plants and to gain understanding of potential range shifts under climate change.

Location

Amazon rain forest.

Methods

We produced species distribution models (SDM) at 5‐km spatial resolution for 42 plant species (trees, palms, lianas, monocot herbs and ferns) using species occurrence data from herbarium records and plot‐based inventories. We modelled species distribution with Bayesian logistic regression using either climate data only, soil data only or climate and soil data together to estimate their relative predictive powers. For areas defined as unsuitable to species occurrence, we mapped the difference between the suitability predictions obtained with climate‐only versus soil‐only models to identify regions where climate and soil might restrict species ranges independently or jointly.

Results

For 40 out of the 42 species, the best models included both climate and soil predictors. The models including only soil predictors performed better than the models including only climate predictors, but we still detected a drought‐sensitive response for most of the species. Edaphic conditions were predicted to restrict species occurrence in the centre, the north‐west and in the north‐east of Amazonia, while the climatic conditions were identified as the restricting factor in the eastern Amazonia, at the border of Roraima and Venezuela and in the Andean foothills.

Main conclusions

Our results revealed that soil data are a more important predictor than climate of plant species range in Amazonia. The strong control of species ranges by edaphic features might reduce species’ abilities to track suitable climate conditions under a drought‐increase scenario. Future challenges are to improve the quality of soil data and couple them with process‐based models to better predict species range dynamics under climate change.     

Where are threatened ferns found? Global conservation priorities for pteridophytes

Recent efforts to improve the representation of plant species included on the IUCN Red List of Threatened Species through the IUCN Sampled Red List Index (SRLI) for Plants have led to the assessment of almost 1000 additional species of pteridophytes and lycophytes under IUCN Red List criteria. Species were selected at random from all lineages of pteridophytes and lycophytes and are taxonomically as well as ecologically representative of pteridophyte and lycophyte diversity. 16% of pteridophyte and lycophyte species are globally threatened with extinction and 22% are of elevated conservation concern (threatened or Near Threatened); of species of pteridophytes and lycophytes previously included on the Red List, 54% were considered threatened. Over half of pteridophyte and lycophyte species assessed for the SRLI use estimates of range size; therefore the method used to measure range may affect the Red List category assigned. We evaluated this using two alternative metrics for estimating range, species distribution modelling (SDM) and ecologically suitable habitat (ESH), for 227 species endemic to the Neotropical biogeographic realm. Differences between range estimates were small when ranges were small but increased with increasing range size. For 58 (25.6%) species alternative modelling techniques result in the species meeting the threshold for a different IUCN Red List category from using extent of occurrence. Modelling threatened species distributions also highlights priority areas for conservation in tropical and subtropical montane forests that are the most species-rich habitat for small-range pteridophyte and lycophyte species, but which are now increasingly subject to rapid conversion to agriculture.     

Predicting the Fate of Biodiversity Using Species�� Distribution Models: Enhancing Model Comparability and Repeatability

Species distribution modeling (SDM) is an increasingly important tool to predict the geographic distribution of species. Even though many problems associated with this method have been highlighted and solutions have been proposed, little has been done to increase comparability among studies. We reviewed recent publications applying SDMs and found that seventy nine percent failed to report methods that ensure comparability among studies, such as disclosing the maximum probability range produced by the models and reporting on the number of species occurrences used. We modeled six species of Falco from northern Europe and demonstrate that model results are altered by (1) spatial bias in species’ occurrence data, (2) differences in the geographic extent of the environmental data, and (3) the effects of transformation of model output to presence/absence data when applying thresholds. Depending on the modeling decisions, forecasts of the future geographic distribution of Falco ranged from range contraction in 80% of the species to no net loss in any species, with the best model predicting no net loss of habitat in Northern Europe. The fact that predictions of range changes in response to climate change in published studies may be influenced by decisions in the modeling process seriously hampers the possibility of making sound management recommendations. Thus, each of the decisions made in generating SDMs should be reported and evaluated to ensure conclusions and policies are based on the biology and ecology of the species being modeled.     

Integrating anthropogenic factors into regional‐scale species distribution models—A novel application in the imperiled sagebrush biome

Species distribution models (SDMs) that rely on regional‐scale environmental variables will play a key role in forecasting species occurrence in the face of climate change. However, in the Anthropocene, a number of local‐scale anthropogenic variables, including wildfire history, land‐use change, invasive species, and ecological restoration practices can override regional‐scale variables to drive patterns of species distribution. Incorporating these human‐induced factors into SDMs remains a major research challenge, in part because spatial variability in these factors occurs at fine scales, rendering prediction over regional extents problematic. Here, we used big sagebrush (Artemisia tridentata Nutt.) as a model species to explore whether including human‐induced factors improves the fit of the SDM. We applied a Bayesian hurdle spatial approach using 21,753 data points of field‐sampled vegetation obtained from the LANDFIRE program to model sagebrush occurrence and cover by incorporating fire history metrics and restoration treatments from 1980 to 2015 throughout the Great Basin of North America. Models including fire attributes and restoration treatments performed better than those including only climate and topographic variables. Number of fires and fire occurrence had the strongest relative effects on big sagebrush occurrence and cover, respectively. The models predicted that the probability of big sagebrush occurrence decreases by 1.2% (95% CI: ?6.9%, 0.6%) when one fire occurs and cover decreases by 44.7% (95% CI: ?47.9%, ?41.3%) if at least one fire occurred over the 36 year period of record. Restoration practices increased the probability of big sagebrush occurrence but had minimal effect on cover. Our results demonstrate the potential value of including disturbance and land management along with climate in models to predict species distributions. As an increasing number of datasets representing land‐use history become available, we anticipate that our modeling framework will have broad relevance across a range of biomes and species.     

An evaluation of four threshold selection methods in species occurrence modelling with random forest: Case studies with Davidia involucrata and Cunninghamia lanceolata

Aims Predictive species distribution models (SDMs) are increasingly applied in resource assessment, environmental conservation and biodiversity management. However, most SDM models often yield a predicted probability (suitability) surface map. In conservation and environmental management practices, the information presented as species presence/absence (binary) may be more practical than presented as probability or suitability. Therefore, a threshold is needed to transform the probability or suitability data to presence/absence data. However, little is known about the effects of different threshold-selection methods on model performance and species range changes induced by future climate. Of the numerous SDM models, random forest (RF) can produce probabilistic and binary species distribution maps based on its regression and classification algorisms, respectively. Studies dealing with the comparative test of the performances of RF regression and classification algorisms have not been reported.
Methods Here, the RF was used to simulate the current and project the future potential distributions of Davidia involucrata and Cunninghamia lanceolata. Then, four threshold-setting methods (Default 0.5, MaxKappa, MaxTSS and MaxACC) were selected and used to transform modelled probabilities of occurrence into binary predictions of species presence and absence. Lastly, we investigated the difference in model performance among the threshold selection methods by using five model accuracy measures (Kappa, TSS, Overall accuracy, Sensitivity and Specificity). We also used the map similarity measure, Kappa, for a cell-by-cell comparison of similarities and differences of distribution map under current and future climates.
Important findings We found that the choice of threshold method altered estimates of model performance, species habitat suitable area and species range shifts under future climate. The difference in selected threshold cut-offs among the four threshold methods was significant for D. involucrata, but was not significant for C. lanceolata. Species’ geographic ranges changed (area change and shifting distance) in response to climate change, but the projections of the four threshold methods did not differ significantly with respect to how much or in which direction, but they did differ against RF classification predictions. The pairwise similarity analysis of binary maps indicated that spatial correspondence among prediction maps was the highest between the MaxKappa and the MaxTSS, and lowest between RF classification algorism and the four threshold-setting methods. We argue that the MaxTSS and the MaxKappa are promising methods for threshold selection when RF regression algorism is used for the distribution modeling of species. This study also provides promising insights to our understanding of the uncertainty of threshold selection in species distribution modeling.     

Quantification of climatic niches in birds: adding the temporal dimension

Quantification of the climatic niche from geographic occurrences is an increasingly important tool for studying species’ relationships to their environment, for example to predict responses to climate change. However, as the geographic distributions of birds are seasonally dynamic, they pose a challenge to carrying out comparable and appropriate quantification of climatic niches. In this review, we first assess how relevant seasonal dynamics are across birds as a whole by compiling a database of migratory behaviour for 10 443 bird species. Second, we examine how studies have quantified climatic niches of birds. Finally, using Australia as a case study, we investigate how well existing distribution datasets represent temporal dynamics by comparing seasonal patterns of species richness obtained from point‐occurrence data with those from range maps and assess the consequences for niche quantification. We provide a consistent classification of migratory behaviour across all birds, and find that a huge variety exists between and within species that should be considered when quantifying climatic niches. Despite this, our review of the literature revealed that seasonal dynamics have often not been accounted for. For future studies, we provide a framework for selecting appropriate occurrence data depending on migratory behaviour and data availability. Our comparison of seasonal species richness patterns obtained from extent‐of‐occurrence range maps and point‐occurrence data suggests that range maps are less able to detect temporal dynamics of bird distributions than point‐occurrence data. We conclude that seasonally explicit range maps combined with climatic data for the corresponding time period can be used to adequately quantify climatic niches for resident birds, but are not adequate to quantify the climatic niches of migratory and nomadic species. Therefore, consistent quantification of climatic niches across all birds requires temporally explicit occurrence points. As such, increasing the availability of occurrence data and methods correcting biases should be a priority.