物种分布模型在海洋潜在生境预测的应用研究进展

海洋生物的栖息分布与环境要素的关联性一直是海洋生态学研究的热点之一.近年来,物种分布模型被广泛应用于预测海洋物种分布、潜在适宜性生境评价等研究,为保护海洋生物多样性、防治外来物种入侵及制定渔业管理措施等提供了一条有效途径.物种分布模型主要包括生境适宜性指数模型、机理模型和统计模型.本文对物种分布模型的理论基础进行了归纳和总结,回顾了物种分布模型在预测海洋物种潜在地理分布研究中的开发与应用,重点介绍了不同类型统计模型在海洋物种潜在分布预测中的研究实例.比较各种选取变量和模型验证方法,认为赤池信息准则对于选取模型变量具有优势,Kappa系数和受试者操作特征曲线下面积在验证模型精度中应用最广泛.阐述了物种分布模型存在的问题及未来发展趋势,随着海洋生物生理机制研究的进一步深入,机理模型将是今后物种分布模型发展的重点.     

Testing projected wild bee distributions in agricultural habitats: predictive power depends on species traits and habitat type

Species distribution models (SDM) are increasingly used to understand the factors that regulate variation in biodiversity patterns and to help plan conservation strategies. However, these models are rarely validated with independently collected data and it is unclear whether SDM performance is maintained across distinct habitats and for species with different functional traits. Highly mobile species, such as bees, can be particularly challenging to model. Here, we use independent sets of occurrence data collected systematically in several agricultural habitats to test how the predictive performance of SDMs for wild bee species depends on species traits, habitat type, and sampling technique. We used a species distribution modeling approach parametrized for the Netherlands, with presence records from 1990 to 2010 for 193 Dutch wild bees. For each species, we built a Maxent model based on 13 climate and landscape variables. We tested the predictive performance of the SDMs with independent datasets collected from orchards and arable fields across the Netherlands from 2010 to 2013, using transect surveys or pan traps. Model predictive performance depended on species traits and habitat type. Occurrence of bee species specialized in habitat and diet was better predicted than generalist bees. Predictions of habitat suitability were also more precise for habitats that are temporally more stable (orchards) than for habitats that suffer regular alterations (arable), particularly for small, solitary bees. As a conservation tool, SDMs are best suited to modeling rarer, specialist species than more generalist and will work best in long‐term stable habitats. The variability of complex, short‐term habitats is difficult to capture in such models and historical land use generally has low thematic resolution. To improve SDMs’ usefulness, models require explanatory variables and collection data that include detailed landscape characteristics, for example, variability of crops and flower availability. Additionally, testing SDMs with field surveys should involve multiple collection techniques.     

Using a new framework of two-phase generalized additive models to incorporate prey abundance in spatial distribution models of juvenile slender lizardfish in Haizhou Bay,China

The predictive skill of species distribution models depends on the quality and quantity of input information. In addition to the physical environmental variables, prey availability is also one of the main drivers regulating spatial distribution of marine species. However, prey distribution data have rarely been considered in habitat models due to the lack of information on non-commercial prey species. This may lead to an incomplete view of species distributions and biased model predictions. In this study, we developed a new framework of two-phase generalized additive models (GAMs) based on the Tweedie distribution to incorporate the predicted prey abundance as covariates in habitat models, and applied this framework to juvenile slender lizardfish Saurida elongata in Haizhou Bay, China. This study demonstrated that the predictive skill of habitat models could be greatly improved through incorporating prey abundance as explanatory variables. The importance of prey distribution data in the habitat model confirms the essentiality of including prey data while modelling species distribution. Spatial overlap and GAM analysis demonstrated that not all dominant prey can be selected as potential explanatory variables and only those prey species showing high spatiotemporal occurrences with predators should be incorporated. The framework derived in this study could be extended to other marine organisms to improve the predictive skill of habitat models and enhance our understanding of the ecological mechanisms underlying the distribution of marine species.     

Complementary strengths of spatially-explicit and multi-species distribution models

Species distribution models (SDMs) project the outcome of community assembly processes – dispersal, the abiotic environment and biotic interactions – onto geographic space. Recent advances in SDMs account for these processes by simultaneously modeling the species that comprise a community in a multivariate statistical framework or by incorporating residual spatial autocorrelation in SDMs. However, the effects of combining both multivariate and spatially-explicit model structures on the ecological inferences and the predictive abilities of a model are largely unknown. We used data on eastern hemlock Tsuga canadensis and five additional co-occurring overstory tree species in 35 569 forest stands across Michigan, USA to evaluate how the choice of model structure, including spatial and non-spatial forms of univariate and multivariate models, affects ecological inference about the processes that shape community composition as well as model predictive ability. Incorporating residual spatial autocorrelation via spatial random effects did not improve out-of-sample prediction for the six tree species, although in-sample model fit was higher in the spatial models. Spatial models attributed less variation in occurrence probability to environmental covariates than the non-spatial models for all six tree species, and estimated higher (more positive) residual co-occurrence values for most species pairs. The non-spatial multivariate model was better suited for evaluating habitat suitability and hypotheses about the processes that shape community composition. Environmental correlations and residual correlations among species pairs were positively related, perhaps indicating that residual correlations were due to shared responses to unmeasured environmental covariates. This work highlights the importance of choosing a non-spatial model formulation to address research questions about the species–environment relationship or residual co-occurrence patterns, and a spatial model formulation when within-sample prediction accuracy is the main goal.     

Microclimate species distribution models estimate lower levels of climate-related habitat loss for salamanders

Species distribution models (SDMs) largely rely on free-air temperatures at coarse spatial resolutions to predict habitat suitability, potentially overlooking important microhabitat. Integrating microclimate data into SDMs may improve predictions of organismal responses to climate change and support targeting of conservation assets at biologically relevant scales, especially for small, dispersal-limited species vulnerable to climate-change-induced range loss. We integrated microclimate data that account for the buffering effects of forest vegetation into SDMs at a very high spatial resolution (3 m²) for three plethodontid salamander species in Great Smoky Mountains National Park (North Carolina and Tennessee). Microclimate SDMs were used to characterize potential changes to future plethodontid habitat, including habitat suitability and habitat spatial patterns. Additionally, we evaluated spatial discrepancies between predictions of habitat suitability developed with microclimate and coarse-resolution, free-air climate data. Microclimate SDMs indicated substantial losses to plethodontid ranges and highly suitable habitat by mid-century, but at much more conservative levels than coarse-resolution models. Coarse-resolution SDMs generally estimated higher mid-century losses to plethodontid habitat compared to microclimate models and consistently undervalued areas containing highly suitable microhabitat. Furthermore, microclimate SDMs revealed potential areas of future gain in highly suitable habitat within current species’ ranges, which may serve as climatic microrefugia. Taken together, this study highlights the need to develop microclimate SDMs that account for vegetation and its biophysical effects on near-surface temperatures. As microclimate datasets become increasingly available across the world, their integration into correlative and mechanistic SDMs will be imperative for accurately estimating organismal responses to climate change and helping environmental managers tasked with spatially prioritizing conservation assets.     

Testing whether ensemble modelling is advantageous for maximising predictive performance of species distribution models

Predictive performance is important to many applications of species distribution models (SDMs). The SDM ‘ensemble’ approach, which combines predictions across different modelling methods, is believed to improve predictive performance, and is used in many recent SDM studies. Here, we aim to compare the predictive performance of ensemble species distribution models to that of individual models, using a large presence–absence dataset of eucalypt tree species. To test model performance, we divided our dataset into calibration and evaluation folds using two spatial blocking strategies (checkerboard-pattern and latitudinal slicing). We calibrated and cross-validated all models within the calibration folds, using both repeated random division of data (a common approach) and spatial blocking. Ensembles were built using the software package ‘biomod2’, with standard (‘untuned’) settings. Boosted regression tree (BRT) models were also fitted to the same data, tuned according to published procedures. We then used evaluation folds to compare ensembles against both their component untuned individual models, and against the BRTs. We used area under the receiver-operating characteristic curve (AUC) and log-likelihood for assessing model performance. In all our tests, ensemble models performed well, but not consistently better than their component untuned individual models or tuned BRTs across all tests. Moreover, choosing untuned individual models with best cross-validation performance also yielded good external performance, with blocked cross-validation proving better suited for this choice, in this study, than repeated random cross-validation. The latitudinal slice test was only possible for four species; this showed some individual models, and particularly the tuned one, performing better than ensembles. This study shows no particular benefit to using ensembles over individual tuned models. It also suggests that further robust testing of performance is required for situations where models are used to predict to distant places or environments.     

Species distribution models as a tool to estimate reproductive parameters: a case study with a passerine bird species

1.?Correlative species distribution models (SDMs) assess relationships between species distribution data and environmental features, to evaluate the environmental suitability (ES) of a given area for a species, by providing a measure of the probability of presence. If the output of SDMs represents the relationships between habitat features and species performance well, SDM results can be related also to other key parameters of populations, including reproductive parameters. To test this hypothesis, we evaluated whether SDM results can be used as a proxy of reproductive parameters (breeding output, territory size) in red-backed shrikes (Lanius collurio). 2.?The distribution of 726 shrike territories in Northern Italy was obtained through multiple focused surveys; for a subset of pairs, we also measured territory area and number of fledged juveniles. We used Maximum Entropy modelling to build a SDM on the basis of territory distribution. We used generalized least squares and spatial generalized mixed models to relate territory size and number of fledged juveniles to SDM suitability, while controlling for spatial autocorrelation. 3.?Species distribution models predicted shrike distribution very well. Territory size was negatively related to suitability estimated through SDM, while the number of fledglings significantly increased with the suitability of the territory. This was true also when SDM was built using only spatially and temporally independent data. 4.?Results show a clear relationship between ES estimated through presence-only SDMs and two key parameters related to species' reproduction, suggesting that suitability estimated by SDM, and habitat quality determining reproduction parameters in our model system, are correlated. Our study shows the potential use of SDMs to infer important fitness parameters; this information can have great importance in management and conservation.     

Livestock-driven land use change to model species distributions: Egyptian vulture as a case study

Species distribution models (SDMs) are increasingly used to predict species ranges and their shifts under future scenarios of global environmental change (GEC). SDMs are thus incorporating key drivers of GEC (e.g. climate, land use) to improve predictions of species’ habitat suitability (i.e. as an indicator of species occurrence). Yet, most SDMs incorporating land use only consider dominant land cover types, largely ignoring other key aspects of land use such as land management intensity and livestock. We developed SDMs including main land use components (i.e. land cover, livestock and its management intensity) to assess their relative importance in shaping habitat suitability for the Egyptian vulture, an endangered raptor linked to livestock presence. We modelled current and future (2020 and 2050) habitat suitability for this vulture using an organism-centred approach. This allowed us to account for basic species’ habitat needs (i.e. nesting cliff) while gaining insight into our variables of interest (i.e. livestock and land cover). Once nest-site requirements were fulfilled, land use variables (i.e. openland and sheep and goat density) were the main factors determining species’ habitat suitability. Current suitable area could decrease by up to 6.81% by 2050 under scenarios with rapid economic growth but no focus on environmental conservation and rural development. Local solutions to environmental sustainability and rural development could double current habitat suitability by 2050. Land use is expected to play a key role in determining Egyptian vulture's distribution through land cover change but also through changes in livestock management (i.e. species and stocking density). Change in stocking densities (sheep and goats/km²) becomes thus an indicator of habitat suitability for this vulture in our study area. Abandonment of agro-pastoral practises (i.e. below ∼15–20 sheep and goats/km²) will negatively influence the species distribution. Nonetheless, livestock densities above these values will not further increase habitat suitability. Given the widespread impacts of livestock on ecosystems, the role of livestock and its management intensity in SDMs for other (non-livestock-related) species should be further explored.     

Modeling bird species distribution change in fire prone Mediterranean landscapes: incorporating species dispersal and landscape dynamics

Species distribution models (SDMs) have traditionally been founded on the assumption that species distributions are in equilibrium with environmental conditions and that these species–environment relationships can be used to estimate species responses to environmental changes. Insight into the validity of this assumption can be obtained from comparing the performance of correlative species distribution models with more complex hybrid approaches, i.e. correlative and process‐based models that explicitly include ecological processes, thereby accounting for mismatches between habitat suitability and species occupancy patterns. Here we compared the ability of correlative SDMs and hybrid models, which can accommodate non‐equilibrium situations arising from dispersal constraints, to reproduce the distribution dynamics of the ortolan bunting Emberiza hortulana in highly dynamic, early successional, fire driven Mediterranean landscapes. Whereas, habitat availability was derived from a correlative statistical SDM, occupancy was modeled using a hybrid approach combining a grid‐based, spatially‐explicit population model that explicitly included bird dispersal with the correlative model. We compared species occupancy patterns under the equilibrium assumption and different scenarios of species dispersal capabilities. To evaluate the predictive capability of the different models, we used independent species data collected in areas affected to different degree by fires. In accordance with the view that disturbance leads to a disparity between the suitable habitat and the occupancy patterns of the ortolan bunting, our results indicated that hybrid modeling approaches were superior to correlative models in predicting species spatial dynamics. Furthermore, hybrid models that incorporated short dispersal distances were more likely to reproduce the observed changes in ortolan bunting distribution patterns, suggesting that dispersal plays a key role in limiting the colonization of recently burnt areas. We conclude that SDMs used in a dynamic context can be significantly improved by using combined hybrid modeling approaches that explicitly account for interactions between key ecological constraints such as dispersal and habitat suitability that drive species response to environmental changes.     

Testing the ability of species distribution models to infer variable importance

Models of species’ distributions and niches are frequently used to infer the importance of range- and niche-defining variables. However, the degree to which these models can reliably identify important variables and quantify their influence remains unknown. Here we use a series of simulations to explore how well models can 1) discriminate between variables with different influence and 2) calibrate the magnitude of influence relative to an ‘omniscient’ model. To quantify variable importance, we trained generalized additive models (GAMs), Maxent and boosted regression trees (BRTs) on simulated data and tested their sensitivity to permutations in each predictor. Importance was inferred by calculating the correlation between permuted and unpermuted predictions, and by comparing predictive accuracy of permuted and unpermuted predictions using AUC and the continuous Boyce index. In scenarios with one influential and one uninfluential variable, models failed to discriminate reliably between variables when training occurrences were < 8–64, prevalence was > 0.5, spatial extent was small, environmental data had coarse resolution and spatial autocorrelation was low, or when pairwise correlation between environmental variables was |r| > 0.7. When two variables influenced the distribution equally, importance was underestimated when species had narrow or intermediate niche breadth. Interactions between variables in how they shaped the niche did not affect inferences about their importance. When variables acted unequally, the effect of the stronger variable was overestimated. GAMs and Maxent discriminated between variables more reliably than BRTs, but no algorithm was consistently well-calibrated vis-à-vis the omniscient model. Algorithm-specific measures of importance like Maxent's change-in-gain metric were less robust than the permutation test. Overall, high predictive accuracy did not connote robust inferential capacity. As a result, requirements for reliably measuring variable importance are likely more stringent than for creating models with high predictive accuracy.     

Testing the utility of species distribution modelling using Random Forests for a species in decline

Habitat suitability estimates derived from species distribution models (SDMs) are increasingly used to guide management of threatened species. Poorly estimating species’ ranges can lead to underestimation of threatened status, undervaluing of remaining habitat and misdirection of conservation funding. We aimed to evaluate the utility of a SDM, similar to the models used to inform government regulation of habitat in our study region, in estimating the contemporary distribution of a threatened and declining species. We developed a presence‐only SDM for the endangered New Holland Mouse (Pseudomys novaehollandiae) across Victoria, Australia. We conducted extensive camera trap surveys across model‐predicted and expert‐selected areas to generate an independent data set for use in evaluating the model, determining confidence in absence data from non‐detection sites with occupancy and detectability modelling. We assessed the predictive capacity of the model at thresholds based on (1) sum of sensitivity and specificity (SSS), and (2) the lowest presence threshold (LPT; i.e. the lowest non‐zero model‐predicted habitat suitability value at which we detected the species). We detected P. novaehollandiae at 40 of 472 surveyed sites, with strong support for the species’ probable absence from non‐detection sites. Based on our post hoc optimised SSS threshold of the SDM, 25% of our detection sites were falsely predicted as non‐suitable habitat and 75% of sites predicted as suitable habitat did not contain the species at the time of our survey. One occupied site had a model‐predicted suitability value of zero, and at the LPT, 88% of sites predicted as suitable habitat did not contain the species at the time of our survey. Our findings demonstrate that application of generic SDMs in both regulatory and investment contexts should be tempered by considering their limitations and currency. Further, we recommend engaging species experts in the extrapolation and application of SDM outputs.     

Terrestrial or marine species distribution model: Why not both? A case study with seabirds

Species reliant on both the terrestrial and marine realms present a challenge for conventional species distribution models (SDMs). For such species, standard single‐realm SDMs may omit key information that could result in decreased model accuracy and performance. Existing approaches to habitat suitability modeling typically do not effectively combine information from multiple realms; this methodological gap can ultimately hamper management efforts for groups such as seabirds, seals, and turtles. This study, for the first time, jointly incorporates both terrestrial information and marine information into a single species distribution model framework. We do this by sampling nearby marine conditions for a given terrestrial point and vice versa using parameters set by each species’ mean maximum foraging distance and then use standard SDM methods to generate habitat suitability predictions; therefore, our method does not rely on post hoc combination of several different models. Using three seabird species with very different ecologies, we investigate whether this new multi‐realm approach can improve our ability to identify suitable habitats for these species. Results show that incorporating terrestrial information into marine SDMs, or vice versa, generally improves model performance, sometimes drastically. However, there is considerable variability between species in the level of improvement as well as in the particular method that produces the most improvement. Our approach provides a repeatable and transparent method to combine information from multiple ecological realms in a single SDM framework. Important advantages over existing solutions include the opportunity to, firstly, easily combine terrestrial and marine information for species that forage large distances inland or out to sea and, secondly, consider interactions between terrestrial and marine variables.     

Relationships between survival and habitat suitability of semi‐aquatic mammals

Spatial distribution and habitat selection are integral to the study of animal ecology. Habitat selection may optimize the fitness of individuals. Hutchinsonian niche theory posits the fundamental niche of species would support the persistence or growth of populations. Although niche‐based species distribution models (SDMs) and habitat suitability models (HSMs) such as maximum entropy (Maxent) have demonstrated fair to excellent predictive power, few studies have linked the prediction of HSMs to demographic rates. We aimed to test the prediction of Hutchinsonian niche theory that habitat suitability (i.e., likelihood of occurrence) would be positively related to survival of American beaver (Castor canadensis), a North American semi‐aquatic, herbivorous, habitat generalist. We also tested the prediction of ideal free distribution that animal fitness, or its surrogate, is independent of habitat suitability at the equilibrium. We estimated beaver monthly survival probability using the Barker model and radio telemetry data collected in northern Alabama, United States from January 2011 to April 2012. A habitat suitability map was generated with Maxent for the entire study site using landscape variables derived from the 2011 National Land Cover Database (30‐m resolution). We found an inverse relationship between habitat suitability index and beaver survival, contradicting the predictions of niche theory and ideal free distribution. Furthermore, four landscape variables selected by American beaver did not predict survival. The beaver population on our study site has been established for 20 or more years and, subsequently, may be approaching or have reached the carrying capacity. Maxent‐predicted increases in habitat use and subsequent intraspecific competition may have reduced beaver survival. Habitat suitability‐fitness relationships may be complex and, in part, contingent upon local animal abundance. Future studies of mechanistic SDMs incorporating local abundance and demographic rates are needed.     

From topography to hydrology—The modifiable area unit problem impacts freshwater species distribution models

Species distribution models (SDMs) are statistical tools to identify potentially suitable habitats for species. For SDMs in river ecosystems, species occurrences and predictor data are often aggregated across subcatchments that serve as modeling units. The level of aggregation (i.e., model resolution) influences the statistical relationships between species occurrences and environmental predictors—a phenomenon known as the modifiable area unit problem (MAUP), making model outputs directly contingent on the model resolution. Here, we test how model performance, predictor importance, and the spatial congruence of species predictions depend on the model resolution (i.e., average subcatchment size) of SDMs. We modeled the potential habitat suitability of 50 native fish species in the upper Danube catchment at 10 different model resolutions. Model resolutions were derived using a 90‐m digital‐elevation model by using the GRASS‐GIS module r.watershed. Here, we decreased the average subcatchment size gradually from 632 to 2 km². We then ran ensemble SDMs based on five algorithms using topographical, climatic, hydrological, and land‐use predictors for each species and resolution. Model evaluation scores were consistently high, as sensitivity and True Skill Statistic values ranged from 86.1–93.2 and 0.61–0.73, respectively. The most contributing predictor changed from topography at coarse, to hydrology at fine resolutions. Climate predictors played an intermediate role for all resolutions, while land use was of little importance. Regarding the predicted habitat suitability, we identified a spatial filtering from coarse to intermediate resolutions. The predicted habitat suitability within a coarse resolution was not ported to all smaller, nested subcatchments, but only to a fraction that held the suitable environmental conditions. Across finer resolutions, the mapped predictions were spatially congruent without such filter effect. We show that freshwater SDM predictions can have consistently high evaluation scores while mapped predictions differ significantly and are highly contingent on the underlying subcatchment size. We encourage building freshwater SDMs across multiple catchment sizes, to assess model variability and uncertainties in model outcomes emerging from the MAUP.     

Utilizing the density of inventory samples to define a hybrid lattice for species distribution models: DISTRIB‐II for 135 eastern U.S. trees

相似文献   

Static species distribution models in dynamically changing systems: how good can predictions really be?

It is widely acknowledged that species respond to climate change by range shifts. Robust predictions of such changes in species’ distributions are pivotal for conservation planning and policy making, and are thus major challenges in ecological research. Statistical species distribution models (SDMs) have been widely applied in this context, though they remain subject to criticism as they implicitly assume equilibrium, and incorporate neither dispersal, demographic processes nor biotic interactions explicitly. In this study, the effects of transient dynamics and ecological properties and processes on the prediction accuracy of SDMs for climate change projections were tested. A spatially explicit multi‐species dynamic population model was built, incorporating species‐specific and interspecific ecological processes, environmental stochasticity and climate change. Species distributions were sampled in different scenarios, and SDMs were estimated by applying generalised linear models (GLMs) and boosted regression trees (BRTs). Resulting model performances were related to prevailing ecological processes and temporal dynamics. SDM performance varied for different range dynamics. Prediction accuracies decreased when abrupt range shifts occurred as species were outpaced by the rate of climate change, and increased again when a new equilibrium situation was realised. When ranges contracted, prediction accuracies increased as the absences were predicted well. Far‐dispersing species were faster in tracking climate change, and were predicted more accurately by SDMs than short‐dispersing species. BRTs mostly outperformed GLMs. The presence of a predator, and the inclusion of its incidence as an environmental predictor, made BRTs and GLMs perform similarly. Results are discussed in light of other studies dealing with effects of ecological traits and processes on SDM performance. Perspectives are given on further advancements of SDMs and for possible interfaces with more mechanistic approaches in order to improve predictions under environmental change.     

Exploring the spatiotemporal dynamics of habitat suitability to improve conservation management of a vulnerable plant species

Ongoing declines in biodiversity caused by global environmental changes call for adaptive conservation management, including the assessment of habitat suitability spatiotemporal dynamics potentially affecting species persistence. Remote sensing (RS) provides a wide-range of satellite-based environmental variables that can be fed into species distribution models (SDMs) to investigate species-environment relations and forecast responses to change. We address the spatiotemporal dynamics of species’ habitat suitability at the landscape level by combining multi-temporal RS data with SDMs for analysing inter-annual habitat suitability dynamics. We implemented this framework with a vulnerable plant species (Veronica micrantha), by combining SDMs with a time-series of RS-based metrics of vegetation functioning related to primary productivity, seasonality, phenology and actual evapotranspiration. Besides RS variables, predictors related to landscape structure, soils and wildfires were ranked and combined through multi-model inference (MMI). To assess recent dynamics, a habitat suitability time-series was generated through model hindcasting. MMI highlighted the strong predictive ability of RS variables related to primary productivity and water availability for explaining the test-species distribution, along with soil, wildfire regime and landscape composition. The habitat suitability time-series revealed the effects of short-term land cover changes and inter-annual variability in climatic conditions. Multi-temporal SDMs further improved predictions, benefiting from RS time-series. Overall, results emphasize the integration of landscape attributes related to function, composition and spatial configuration for improving the explanation of ecological patterns. Moreover, coupling SDMs with RS functional metrics may provide early-warnings of future environmental changes potentially impacting habitat suitability. Applications discussed include the improvement of biodiversity monitoring and conservation strategies.     

Ecological principles of species distribution models: the habitat matching rule

Most species distribution models (SDMs) assume that habitats are closed, stable and without competition. In that environmental context, it is ecologically correct to assume that members of a species will be distributed in direct relation to the suitability of the habitat, that is, according to the so‐called habitat matching rule. This paper examines whether it is possible to maintain the assumption of the habitat matching rule in the following circumstances: (1) when habitats are connected and organisms can move between them, (2) when there are disturbances and seasonal cycles that generate instability, and (3) when there is inter‐specific and intra‐specific competition. Here I argue that it is possible as long as the following aspects are taken into account. In open habitats at equilibrium, in which habitat selection and competition operate, the habitat matching rule can be applied in some conditions, while competition tends to homogenize the species distribution in other environmental contexts. In the latter case, two methods can be used to incorporate these effects into SDMs: new parameters can be incorporated into the response functions, or the occurrence of proportions of categories of individuals (adult/young, male/female, or dominant/subordinate species in guilds) can be used instead of the occurrence of organisms. The habitat matching rule is not fulfilled in non‐equilibrium environments. The solution to this problem lies in the design of SDMs with two strategies that depend on scale. Locally, the disequilibrium can be encapsulated using average environmental conditions, with sufficiently large cells (in the case of metapopulations) and/or long enough sampling periods (in the case of seasonal cycles). At coarse scales, the use of presence‐only models can in some cases avoid the destabilizing effect of catastrophic historical processes. The matching law is a strong assumption of SDMs because it is based on population ecology theory and the principle of evolution by natural selection.     

Transferable species distribution modelling: Comparative performance of Generalised Functional Response models

Predictive species distribution models (SDMs) are becoming increasingly important in ecology, in the light of rapid environmental change. However, the predictions of most current SDMs are specific to the habitat composition of the environments in which they were fitted. This may limit SDM predictive power because species may respond differently to a given habitat depending on the availability of all habitats in their environment, a phenomenon known as a functional response in resource selection. The Generalised Functional Response (GFR) framework captures this dependence by formulating the SDM coefficients as functions of habitat availability. The original GFR implementation used global polynomial functions of habitat availability to describe the functional responses. In this study, we develop several refinements of this approach and compare their predictive performance using two simulated and two real datasets. We first use local radial basis functions (RBF), a more flexible approach than global polynomials, to represent the habitat selection coefficients, and balance bias with precision via regularization to prevent overfitting. Second, we use the RBF-GFR and GFR models in combination with the classification and regression tree CART, which has more flexibility and better predictive powers for non-linear modelling. As further extensions, we use random forests (RFs) and extreme gradient boosting (XGBoost), ensemble approaches that consistently lead to variance reduction in generalization error. We find that the different methods are ranked consistently across the datasets for out-of-data prediction. The traditional stationary approach to SDMs and the GFR model consistently perform at the bottom of the ranking (simple SDMs underfit, and polynomial GFRs overfit the data). The best methods in our list provide non-negligible improvements in predictive performance, in some cases taking the out-of-sample R2 from 0.3 up to 0.7 across datasets. At times of rapid environmental change and spatial non-stationarity ignoring the effects of functional responses on SDMs, results in two different types of prediction bias (under-prediction or mis-positioning of distribution hotspots). However, not all functional response models perform equally well. The more volatile polynomial GFR models can generate biases through over-prediction. Our results indicate that there are consistently robust GFR approaches that achieve impressive gains in transferability across very different datasets.     

Calibration of hybrid species distribution models: the value of general‐purpose vs. targeted monitoring data

Aim Temporally replicated observations are essential for the calibration and validation of species distribution models (SDMs) aiming at making temporal extrapolations. We study here the usefulness of a general‐purpose monitoring programme for the calibration of hybrid SDMs. As a benchmark case, we take the calibration with data from a monitoring programme that specifically surveys those areas where environmental changes expected to be relevant occur. Location Catalonia, north‐east of Spain. Methods We modelled the distribution changes of twelve open‐habitat bird species in landscapes whose dynamics are driven by fire and forest regeneration. We developed hybrid SDMs combining correlative habitat suitability with mechanistic occupancy models. We used observations from two monitoring programmes to provide maximum‐likelihood estimates for spread parameters: a common breeding bird survey (CBS) and a programme specifically designed to monitor bird communities within areas affected by wildfires (DINDIS). Results Both calibration with CBS and DINDIS data yielded sound spread parameter estimates and range dynamics that suggested dispersal limitations. However, compared to calibration with DINDIS data, calibration with CBS data leads to biased estimates of spread distance for seven species and to a higher degree of uncertainty in predicted range dynamics for six species. Main conclusions We have shown that available monitoring data can be used in the calibration of the mechanistic component of hybrid SDMs. However, if the dynamics of the target species occur within areas not well covered, general‐purpose monitoring data can lead to biased and inaccurate parameter estimates. To determine the potential usefulness of a given monitoring data set for the calibration of the mechanistic component of a hybrid SDM, we recommend quantifying the number of surveyed sites that are predicted to undergo habitat suitability changes.