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 modeled 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 modeled niche optima obtained from European-scale SDMs of 1476 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 modeled 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 modeled 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.     

The effect of positional error on fine scale species distribution models increases for specialist species

Species occurrences inherently include positional error. Such error can be problematic for species distribution models (SDMs), especially those based on fine-resolution environmental data. It has been suggested that there could be a link between the influence of positional error and the width of the species ecological niche. Although positional errors in species occurrence data may imply serious limitations, especially for modelling species with narrow ecological niche, it has never been thoroughly explored. We used a virtual species approach to assess the effects of the positional error on fine-scale SDMs for species with environmental niches of different widths. We simulated three virtual species with varying niche breadth, from specialist to generalist. The true distribution of these virtual species was then altered by introducing different levels of positional error (from 5 to 500 m). We built generalized linear models and MaxEnt models using the distribution of the three virtual species (unaltered and altered) and a combination of environmental data at 5 m resolution. The models’ performance and niche overlap were compared to assess the effect of positional error with varying niche breadth in the geographical and environmental space. The positional error negatively impacted performance and niche overlap metrics. The amplitude of the influence of positional error depended on the species niche, with models for specialist species being more affected than those for generalist species. The positional error had the same effect on both modelling techniques. Finally, increasing sample size did not mitigate the negative influence of positional error. We showed that fine-scale SDMs are considerably affected by positional error, even when such error is low. Therefore, where new surveys are undertaken, we recommend paying attention to data collection techniques to minimize the positional error in occurrence data and thus to avoid its negative effect on SDMs, especially when studying specialist species.     

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.     

Equilibrium or not? Modelling potential distribution of invasive species in different stages of invasion

Aim The assumption of equilibrium between organisms and their environment is a standard working postulate in species distribution models (SDMs). However, this assumption is typically violated in models of biological invasions where range expansions are highly constrained by dispersal and colonization processes. Here, we examined how stage of invasion affects the extent to which occurrence data represent the ecological niche of organisms and, in turn, influences spatial prediction of species’ potential distributions. Location Six ecoregions in western Oregon, USA. Methods We compiled occurrence data from 697 field plots collected over a 9‐year period (2001–09) of monitoring the spread of invasive forest pathogen Phytophthora ramorum. Using these data, we applied ecological‐niche factor analysis to calibrate models of potential distribution across different years of colonization. We accounted for natural variation and uncertainties in model evaluation by further investigating three hypothetical scenarios of varying equilibrium in a simulated virtual species, for which the ‘true’ potential distribution was known. Results We confirm our hypothesis that SDMs calibrated in early stages of invasion are less accurate than models calibrated under scenarios closer to equilibrium. SDMs that are developed in early stages of invasion tend to underpredict the potential range compared to models that are built in later stages of invasion. Main conclusions A full environmental niche of invasive species cannot be effectively captured with data from a realized distribution that is restricted by processes preventing full occupancy of suitable habitats. If SDMs are to be used effectively in conservation and management, stage of invasion needs to be considered to avoid underestimation of habitats at risk of invasion.     

Predicting the occurrence of rare Brazilian birds with species distribution models

Species distribution models (SDMs) yield reliable and needed predictions to identify regions that have similar environmental conditions and were used here to predict potential ranges of rare species to identify new localities were they might occur based on their occurrence probability (i.e. niche suitability). We modeled the potential distribution ranges of ten endangered or rare birds from the South American Cerrado biome, using four temperature- and four precipitation-related bioclimatic variables, three topographical variables, and nine different niche modeling methods for each species. We used an ensemble-forecasting approach to reach a consensus scenario to obtain the average distribution for each species based on the five best models generating a distribution map of each species. Model efficiency was related to sample size and not appropriate below ten independent spatial occurrences. The potential distributions of seven species revealed that their occurrence ranges might go beyond their known ranges, but that most of them seem to occur near the regions where they have already been reported. The models of only three species were considered unsatisfactory in helping identify their potential distribution. Models created maps with higher occurrence probability regions where rare Cerrado birds might occur. These range projections can potentially decrease the costs and improve the efficiency of future field searches. On methodological terms, the application of SDMs to predict species ranges should compare different modeling methods and evaluate the effect of sample size on their performance.     

Alien Invasive Slider Turtle in Unpredicted Habitat: A Matter of Niche Shift or of Predictors Studied?

Background

Species Distribution Models (SDMs) aim on the characterization of a species'' ecological niche and project it into geographic space. The result is a map of the species'' potential distribution, which is, for instance, helpful to predict the capability of alien invasive species. With regard to alien invasive species, recently several authors observed a mismatch between potential distributions of native and invasive ranges derived from SDMs and, as an explanation, ecological niche shift during biological invasion has been suggested. We studied the physiologically well known Slider turtle from North America which today is widely distributed over the globe and address the issue of ecological niche shift versus choice of ecological predictors used for model building, i.e., by deriving SDMs using multiple sets of climatic predictor.

Principal Findings

In one SDM, predictors were used aiming to mirror the physiological limits of the Slider turtle. It was compared to numerous other models based on various sets of ecological predictors or predictors aiming at comprehensiveness. The SDM focusing on the study species'' physiological limits depicts the target species'' worldwide potential distribution better than any of the other approaches.

Conclusion

These results suggest that a natural history-driven understanding is crucial in developing statistical models of ecological niches (as SDMs) while “comprehensive” or “standard” sets of ecological predictors may be of limited use.     

Spatial modelling of ecological indicator values improves predictions of plant distributions in complex landscapes

Ecologically meaningful predictors are often neglected in plant distribution studies, resulting in incomplete niche quantification and low predictive power of species distribution models (SDMs). Because environmental data are rare and expensive to collect, and because their relationship with local climatic and topographic conditions are complex, mapping them over large geographic extents and at high spatial resolution remains a major challenge. Here, we propose to derive environmental data layers by mapping ecological indicator values in space. We combined ~6 million plant occurrences with expert-based plant ecological indicator values (EIVs) of 3600 species in Switzerland. EIVs representing local soil properties (pH, moisture, moisture variability, aeration, humus and nutrients) and climatic conditions (continentality, light) were modelled at 93 m spatial resolution with the Random Forest algorithm and 16 predictors representing meso-climate, land use, topography and geology. Models were evaluated and predictions of EIVs were compared with soil inventory data. We mapped each EIV separately and evaluated EIV importance in explaining the distribution of 500 plant species using SDMs with a set of 30 environmental predictors. Finally, we tested how they improve an ensemble of SDMs compared to a standard set of predictors for ca 60 plant species. All EIV models showed excellent performance (|r| > 0.9) and predictions were correlated reasonably (|r| > 0.4) to soil properties measured in the field. Resulting EIV maps were among the most important predictors in SDMs. Also, in ensemble SDMs overall predictive performance increased, mainly through improved model specificity reducing species range overestimation. Combining large citizen science databases to expert-based EIVs is a powerful and cost–effective approach for generalizing local edaphic and climatic conditions over large areas. Producing ecologically meaningful predictors is a first step for generating better predictions of species distribution which is of main importance for decision makers in conservation and environmental management projects.     

Improving species distribution models for climate change studies: variable selection and scale

Statistical species distribution models (SDMs) are widely used to predict the potential changes in species distributions under climate change scenarios. We suggest that we need to revisit the conceptual framework and ecological assumptions on which the relationship between species distributions and environment is based. We present a simple conceptual framework to examine the selection of environmental predictors and data resolution scales. These vary widely in recent papers, with light inconsistently included in the models. Focusing on light as a necessary component of plant SDMs, we briefly review its dependence on aspect and slope and existing knowledge of its influence on plant distribution. Differences in light regimes between north‐ and south‐facing aspects in temperate latitudes can produce differences in temperature equivalent to moves 200 km polewards. Local topography may create refugia that are not recognized in many climate change SDMs using coarse‐scale data. We argue that current assumptions about the selection of predictors and data resolution need further testing. Application of these ideas can clarify many issues of scale, extent and choice of predictors, and potentially improve the use of SDMs for climate change modelling of biodiversity.     

Weak climatic associations among British plant distributions

Aim Species distribution models (SDMs) are used to infer niche responses and predict climate change‐induced range shifts. However, their power to distinguish real and chance associations between spatially autocorrelated distribution and environmental data at continental scales has been questioned. Here this is investigated at a regional (10 km) scale by modelling the distributions of 100 plant species native to the UK. Location UK. Methods SDMs fitted using real climate data were compared with those utilizing simulated climate gradients. The simulated gradients preserve the exact values and spatial structure of the real ones, but have no causal relationships with any species and so represent an appropriate null model. SDMs were fitted as generalized linear models (GLMs) or by the Random Forest machine‐learning algorithm and were either non‐spatial or included spatially explicit trend surfaces or autocovariates as predictors. Results Species distributions were significantly but erroneously related to the simulated gradients in 86% of cases (P < 0.05 in likelihood‐ratio tests of GLMs), with the highest error for strongly autocorrelated species and gradients and when species occupied 50% of sites. Even more false effects were found when curvilinear responses were modelled, and this was not adequately mitigated in the spatially explicit models. Non‐spatial SDMs based on simulated climate data suggested that 70–80% of the apparent explanatory power of the real data could be attributable to its spatial structure. Furthermore, the niche component of spatially explicit SDMs did not significantly contribute to model fit in most species. Main conclusions Spatial structure in the climate, rather than functional relationships with species distributions, may account for much of the apparent fit and predictive power of SDMs. Failure to account for this means that the evidence for climatic limitation of species distributions may have been overstated. As such, predicted regional‐ and national‐scale impacts of climate change based on the analysis of static distribution snapshots will require re‐evaluation.     

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.     

Mechanistic niche modelling: combining physiological and spatial data to predict species' ranges

Species distribution models (SDMs) use spatial environmental data to make inferences on species' range limits and habitat suitability. Conceptually, these models aim to determine and map components of a species' ecological niche through space and time, and they have become important tools in pure and applied ecology and evolutionary biology. Most approaches are correlative in that they statistically link spatial data to species distribution records. An alternative strategy is to explicitly incorporate the mechanistic links between the functional traits of organisms and their environments into SDMs. Here, we review how the principles of biophysical ecology can be used to link spatial data to the physiological responses and constraints of organisms. This provides a mechanistic view of the fundamental niche which can then be mapped to the landscape to infer range constraints. We show how physiologically based SDMs can be developed for different organisms in different environmental contexts. Mechanistic SDMs have different strengths and weaknesses to correlative approaches, and there are many exciting and unexplored prospects for integrating the two approaches. As physiological knowledge becomes better integrated into SDMs, we will make more robust predictions of range shifts in novel or non-equilibrium contexts such as invasions, translocations, climate change and evolutionary shifts.     

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.     

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

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

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.     

Average niche breadths of species in lake macrophyte communities respond to ecological gradients variably in four regions on two continents

Different species’ niche breadths in relation to ecological gradients are infrequently examined within the same study and, moreover, species niche breadths have rarely been averaged to account for variation in entire ecological communities. We investigated how average environmental niche breadths (climate, water quality and climate–water quality niches) in aquatic macrophyte communities are related to ecological gradients (latitude, longitude, altitude, species richness and lake area) among four distinct regions (Finland, Sweden and US states of Minnesota and Wisconsin) on two continents. We found that correlations between the three different measures of average niche breadths and ecological gradients varied considerably among the study regions, with average climate and average water quality niche breadth models often showing opposite trends. However, consistent patterns were also found, such as widening of average climate niche breadths and narrowing of average water quality niche breadths of aquatic macrophytes along increasing latitudinal and altitudinal gradients. This result suggests that macrophyte species are generalists in relation to temperature variations at higher latitudes and altitudes, whereas species in southern, lowland lakes are more specialised. In contrast, aquatic macrophytes growing in more southern nutrient-rich lakes were generalists in relation to water quality, while specialist species are adapted to low-productivity conditions and are found in highland lakes. Our results emphasise that species niche breadths should not be studied using only coarse-scale data of species distributions and corresponding environmental conditions, but that investigations on different kinds of niche breadths (e.g., climate vs. local niches) also require finer resolution data at broad spatial extents.     

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.     

Niche variability and its consequences for species distribution modeling

When species distribution models (SDMs) are used to predict how a species will respond to environmental change, an important assumption is that the environmental niche of the species is conserved over evolutionary time-scales. Empirical studies conducted at ecological time-scales, however, demonstrate that the niche of some species can vary in response to environmental change. We use habitat and locality data of five species of stream fishes collected across seasons to examine the effects of niche variability on the accuracy of projections from Maxent, a popular SDM. We then compare these predictions to those from an alternate method of creating SDM projections in which a transformation of the environmental data to similar scales is applied. The niche of each species varied to some degree in response to seasonal variation in environmental variables, with most species shifting habitat use in response to changes in canopy cover or flow rate. SDMs constructed from the original environmental data accurately predicted the occurrences of one species across all seasons and a subset of seasons for two other species. A similar result was found for SDMs constructed from the transformed environmental data. However, the transformed SDMs produced better models in ten of the 14 total SDMs, as judged by ratios of mean probability values at known presences to mean probability values at all other locations. Niche variability should be an important consideration when using SDMs to predict future distributions of species because of its prevalence among natural populations. The framework we present here may potentially improve these predictions by accounting for such variability.     

No species is an island: testing the effects of biotic interactions on models of avian niche occupation

Traditionally, the niche of a species is described as a hypothetical 3D space, constituted by well‐known biotic interactions (e.g. predation, competition, trophic relationships, resource–consumer interactions, etc.) and various abiotic environmental factors. Species distribution models (SDMs), also called “niche models” and often used to predict wildlife distribution at landscape scale, are typically constructed using abiotic factors with biotic interactions generally been ignored. Here, we compared the goodness of fit of SDMs for red‐backed shrike Lanius collurio in farmlands of Western Poland, using both the classical approach (modeled only on environmental variables) and the approach which included also other potentially associated bird species. The potential associations among species were derived from the relevant ecological literature and by a correlation matrix of occurrences. Our findings highlight the importance of including heterospecific interactions in improving our understanding of niche occupation for bird species. We suggest that suite of measures currently used to quantify realized species niches could be improved by also considering the occurrence of certain associated species. Then, an hypothetical “species 1” can use the occurrence of a successfully established individual of “species 2” as indicator or “trace” of the location of available suitable habitat to breed. We hypothesize this kind of biotic interaction as the “heterospecific trace effect” (HTE): an interaction based on the availability and use of “public information” provided by individuals from different species. Finally, we discuss about the incomes of biotic interactions for enhancing the predictive capacities on species distribution models.     

A stochastic biodiversity model with overlapping niche structure

The niche is a fundamental ecological concept that underpins many explanations of patterns of biodiversity. The complexity of niche processes in ecological systems, however, means that it is difficult to capture them accurately in theoretical models of community assembly. In this study, we build upon simple neutral biodiversity models by adding the important ingredient of overlapping niche structure. Our model is spatially implicit and contains a fixed number of equal-sized habitats. Each species in the metacommunity arises through a speciation event; at which time, it is randomly assigned a fundamental niche or set of environments/habitats in which it can persist. Within each habitat, species compete with other species that have different but overlapping fundamental niches. Species abundances then change through ecological drift; each, however, is constrained by its maximum niche breadth and by the presence of other species in its habitats. Using our model, we derive analytical expressions for steady-state species abundance distributions, steady-state distributions of niche breadth across individuals and across species, and dynamic distributions of niche breadth across species. With this framework, we identify the conditions that produce the log-series species abundance distribution familiar from neutral models. We then identify how overlapping niche structure can lead to other species abundance distributions and, in particular, ask whether these new distributions differ significantly from species abundance distributions predicted by non-overlapping niche models. Finally, we extend our analysis to consider additional distributions associated with realized niche breadths. Overall, our results show that models with overlapping niches can exhibit behavior similar to neutral models, with the caveat that species with narrow fundamental niche breadths will be very rare. If narrow-niche species are common, it must be because they are in a non-overlapping niche or have countervailing advantages over broad-niche species. This result highlights the role that niches can play in establishing demographic neutrality.