Mountain landscapes offer few opportunities for high‐elevation tree species migration

Climate change is anticipated to alter plant species distributions. Regional context, notably the spatial complexity of climatic gradients, may influence species migration potential. While high‐elevation species may benefit from steep climate gradients in mountain regions, their persistence may be threatened by limited suitable habitat as land area decreases with elevation. To untangle these apparently contradictory predictions for mountainous regions, we evaluated the climatic suitability of four coniferous forest tree species of the western United States based on species distribution modeling (SDM) and examined changes in climatically suitable areas under predicted climate change. We used forest structural information relating to tree species dominance, productivity, and demography from an extensive forest inventory system to assess the strength of inferences made with a SDM approach. We found that tree species dominance, productivity, and recruitment were highest where climatic suitability (i.e., probability of species occurrence under certain climate conditions) was high, supporting the use of predicted climatic suitability in examining species risk to climate change. By predicting changes in climatic suitability over the next century, we found that climatic suitability will likely decline, both in areas currently occupied by each tree species and in nearby unoccupied areas to which species might migrate in the future. These trends were most dramatic for high elevation species. Climatic changes predicted over the next century will dramatically reduce climatically suitable areas for high‐elevation tree species while a lower elevation species, Pinus ponderosa, will be well positioned to shift upslope across the region. Reductions in suitable area for high‐elevation species imply that even unlimited migration would be insufficient to offset predicted habitat loss, underscoring the vulnerability of these high‐elevation species to climatic changes.     

Less favourable climates constrain demographic strategies in plants

Correlative species distribution models are based on the observed relationship between species’ occurrence and macroclimate or other environmental variables. In climates predicted less favourable populations are expected to decline, and in favourable climates they are expected to persist. However, little comparative empirical support exists for a relationship between predicted climate suitability and population performance. We found that the performance of 93 populations of 34 plant species worldwide – as measured by in situ population growth rate, its temporal variation and extinction risk – was not correlated with climate suitability. However, correlations of demographic processes underpinning population performance with climate suitability indicated both resistance and vulnerability pathways of population responses to climate: in less suitable climates, plants experienced greater retrogression (resistance pathway) and greater variability in some demographic rates (vulnerability pathway). While a range of demographic strategies occur within species’ climatic niches, demographic strategies are more constrained in climates predicted to be less suitable.     

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.     

Ecological niche models display nonlinear relationships with abundance and demographic performance across the latitudinal distribution of Astragalus utahensis (Fabaceae)

The potential for ecological niche models (ENMs) to accurately predict species' abundance and demographic performance throughout their geographic distributions remains a topic of substantial debate in ecology and biogeography. Few studies simultaneously examine the relationship between ENM predictions of environmental suitability and both a species' abundance and its demographic performance, particularly across its entire geographic distribution. Yet, studies of this type are essential for understanding the extent to which ENMs are a viable tool for identifying areas that may promote high abundance or performance of a species or how species might respond to future climate conditions. In this study, we used an ensemble ecological niche model to predict climatic suitability for the perennial forb Astragalus utahensis across its geographic distribution. We then examined relationships between projected climatic suitability and field‐based measures of abundance, demographic performance, and forecasted stochastic population growth (λ_s). Predicted climatic suitability showed a J‐shaped relationship with A. utahensis abundance, where low‐abundance populations were associated with low‐to‐intermediate suitability scores and abundance increased sharply in areas of high predicted climatic suitability. A similar relationship existed between climatic suitability and λ_s from the center to the northern edge of the latitudinal distribution. Patterns such as these, where density or demographic performance only increases appreciably beyond some threshold of climatic suitability, support the contention that ENM‐predicted climatic suitability does not necessarily represent a reliable predictor of abundance or performance across large geographic regions.     

Topographic,latitudinal and climatic distribution of Pinus coulteri: geographic range limits are not at the edge of the climate envelope

With changing climate, many species are projected to move poleward or to higher elevations to track suitable climates. The prediction that species will move poleward assumes that geographically marginal populations are at the edge of the species' climatic range. We studied Pinus coulteri from the center to the northern (poleward) edge of its range, and examined three scenarios regarding the relationship between the geographic and climatic margins of a species' range. We used herbarium and iNaturalist.org records to identify P. coulteri sites, generated a species distribution model based on temperature, precipitation, climatic water deficit, and actual evapotranspiration, and projected suitability under future climate scenarios. In fourteen populations from the central to northern portions of the range, we conducted field studies and recorded elevation, slope and aspect (to estimate solar insolation) to examine relationships between local and regional distributions. We found that northern populations of P. coulteri do not occupy the cold or wet edge of the species' climatic range; mid‐latitude, high elevation populations occupy the cold margin. Aspect and insolation of P. coulteri populations changed significantly across latitudes and elevations. Unexpectedly, northern, low‐elevation stands occupy north‐facing aspects and receive low insolation, while central, high‐elevation stands grow on more south‐facing aspects that receive higher insolation. Modeled future climate suitability is projected to be highest in the central, high elevation portion of the species range, and in low‐lying coastal regions under some scenarios, with declining suitability in northern areas under most future scenarios. For P. coulteri, the lack of high elevation habitat combined with a major dispersal barrier may limit northward movement in response to a warming climate. Our analyses demonstrate the importance of distinguishing geographically vs. climatically marginal populations, and the importance of quantitative analysis of the realized climate space to understand species range limits.     

Plant performance in central and northern peripheral populations of the widespread Plantago coronopus

Peripheral populations have long been predicted to show lower vital rates, higher demographic fluctuations, and lower densities than central populations. However, recent research has questioned the existence of clear patterns across species’ ranges. To test these hypotheses, we monitored five central and six northern peripheral populations of the widespread herb Plantago coronopus along the European Atlantic coast during 5 yr. We estimated population density, and calculated mean values and temporal variability of four vital rates (survival, individual growth, fecundity and recruitment) in hundreds of plants in permanent plots. Central populations showed higher fecundity, whereas peripheral populations had higher recruitment per reproductive plant, indicating a higher overall reproductive success in the periphery. Central populations showed a marginally significant tendency for higher growth, and there were no differences between range positions in survival. Fecundity and growth were affected by intraspecific competition, and recruitment was affected by precipitation, highlighting the importance of local environmental conditions for population performance. Central and peripheral populations showed no significant differences in temporal variability of vital rates. Finally, density was significantly higher in peripheral than in central populations, in discrepancy with the abundant‐centre model. Density was correlated to seedling recruitment, which would counterbalance in peripheral populations the lower fecundity and the tendency for lower growth of established plants. Such compensations among vital rates might be particularly common in widespread plants, and advise against simplistic assumptions of population performance across ranges. The whole species’ life cycle should be considered, since different arrangements of vital rates are expected to maximize fitness in local environments. Our results show also the importance of discerning between geographical periphery and ecological marginality. In a context of climate‐induced range shifts, these considerations are crucial for the reliability of niche‐models and the management of plant peripheral populations.     

Incorporating intraspecific variation into species distribution models improves distribution predictions,but cannot predict species traits for a wide-spread plant species

The most common approach to predicting how species ranges and ecological functions will shift with climate change is to construct correlative species distribution models (SDMs). These models use a species’ climatic distribution to determine currently suitable areas for the species and project its potential distribution under future climate scenarios. A core, rarely tested, assumption of SDMs is that all populations will respond equivalently to climate. Few studies have examined this assumption, and those that have rarely dissect the reasons for intraspecific differences. Focusing on the arctic-alpine cushion plant Silene acaulis, we compared predictive accuracy from SDMs constructed using the species’ full global distribution with composite predictions from separate SDMs constructed using subpopulations defined either by genetic or habitat differences. This is one of the first studies to compare multiple ways of constructing intraspecific-level SDMs with a species-level SDM. We also examine the contested relationship between relative probability of occurrence and species performance or ecological function, testing if SDM output can predict individual performance (plant size) and biotic interactions (facilitation). We found that both genetic- and habitat-informed SDMs are considerably more accurate than a species-level SDM, and that the genetic model substantially differs from and outperforms the habitat model. While SDMs have been used to infer population performance and possibly even biotic interactions, in our system these relationships were extremely weak. Our results indicate that individual subpopulations may respond differently to climate, although we discuss and explore several alternative explanations for the superior performance of intraspecific-level SDMs. We emphasize the need to carefully examine how to best define intraspecific-level SDMs as well as how potential genetic, environmental, or sampling variation within species ranges can critically affect SDM predictions. We urge caution in inferring population performance or biotic interactions from SDM predictions, as these often-assumed relationships are not supported in our study.     

Uncertainty in predictions of range dynamics: black grouse climbing the Swiss Alps

Empirical species distribution models (SDMs) constitute often the tool of choice for the assessment of rapid climate change effects on species’ vulnerability. Conclusions regarding extinction risks might be misleading, however, because SDMs do not explicitly incorporate dispersal or other demographic processes. Here, we supplement SDMs with a dynamic population model 1) to predict climate‐induced range dynamics for black grouse in Switzerland, 2) to compare direct and indirect measures of extinction risks, and 3) to quantify uncertainty in predictions as well as the sources of that uncertainty. To this end, we linked models of habitat suitability to a spatially explicit, individual‐based model. In an extensive sensitivity analysis, we quantified uncertainty in various model outputs introduced by different SDM algorithms, by different climate scenarios and by demographic model parameters. Potentially suitable habitats were predicted to shift uphill and eastwards. By the end of the 21st century, abrupt habitat losses were predicted in the western Prealps for some climate scenarios. In contrast, population size and occupied area were primarily controlled by currently negative population growth and gradually declined from the beginning of the century across all climate scenarios and SDM algorithms. However, predictions of population dynamic features were highly variable across simulations. Results indicate that inferring extinction probabilities simply from the quantity of suitable habitat may underestimate extinction risks because this may ignore important interactions between life history traits and available habitat. Also, in dynamic range predictions uncertainty in SDM algorithms and climate scenarios can become secondary to uncertainty in dynamic model components. Our study emphasises the need for principal evaluation tools like sensitivity analysis in order to assess uncertainty and robustness in dynamic range predictions. A more direct benefit of such robustness analysis is an improved mechanistic understanding of dynamic species’ responses to climate change.     

Applying landscape metrics to species distribution model predictions to characterize internal range structure and associated changes

Distributional shifts in species ranges provide critical evidence of ecological responses to climate change. Assessments of climate-driven changes typically focus on broad-scale range shifts (e.g. poleward or upward), with ecological consequences at regional and local scales commonly overlooked. While these changes are informative for species presenting continuous geographic ranges, many species have discontinuous distributions—both natural (e.g. mountain or coastal species) or human-induced (e.g. species inhabiting fragmented landscapes)—where within-range changes can be significant. Here, we use an ecosystem engineer species (Sabellaria alveolata) with a naturally fragmented distribution as a case study to assess climate-driven changes in within-range occupancy across its entire global distribution. To this end, we applied landscape ecology metrics to outputs from species distribution modelling (SDM) in a novel unified framework. SDM predicted a 27.5% overall increase in the area of potentially suitable habitat under RCP 4.5 by 2050, which taken in isolation would have led to the classification of the species as a climate change winner. SDM further revealed that the latitudinal range is predicted to shrink because of decreased habitat suitability in the equatorward part of the range, not compensated by a poleward expansion. The use of landscape ecology metrics provided additional insights by identifying regions that are predicted to become increasingly fragmented in the future, potentially increasing extirpation risk by jeopardising metapopulation dynamics. This increased range fragmentation could have dramatic consequences for ecosystem structure and functioning. Importantly, the proposed framework—which brings together SDM and landscape metrics—can be widely used to study currently overlooked climate-driven changes in species internal range structure, without requiring detailed empirical knowledge of the modelled species. This approach represents an important advancement beyond predictive envelope approaches and could reveal itself as paramount for managers whose spatial scale of action usually ranges from local to regional.     

Invasion hotspots for non‐native plants in Australia under current and future climates

We apply the concept of biodiversity hotspot analysis (the identification of biogeographical regions of high species diversity) to identify invasion hotspots – areas of potentially suitable climate for multiple non‐native plant species – in Australia under current and future climates. We used the species distribution model Maxent to model climate suitability surfaces for 72 taxa, recognized as ‘Weeds of National Significance’ (WoNS) in Australia, under current and projected climate for 2020 and 2050. Current climate suitability layers were summed across all 72 species, and we observed two regions of high climatic suitability corresponding to the top 25th percentile of combined climatic suitability values across Australia. We defined these as potential invasion hotspots. Areas of climatic suitability equivalent to the hotspot regions were identified in the composite maps for 2020 and 2050, to track spatial changes in the hotspots over the two time steps. Two potential invasion hotspot regions were identified under current and projected climates: the south west corner of Western Australia (SW), and south eastern Australia (SE). Herbarium data confirmed the presence of 73% and 99% of those species predicted to be in each hotspot respectively, suggesting that the SE has greater invasion potential. The area of both hotspots was predicted to retract southward and towards the coast under future climate scenarios, reducing in size by 81% (SW) and 71% (SE) by 2050. This reduction was driven by the dominance of southern temperate invasive plant species in the WoNS list (47 of the 72), of which 44 were predicted to experience reductions in their bioclimatic range by 2050. While climate is likely to become less suitable for the majority of WoNS in the future, potential invasion hotspots based on climate suitability are likely to remain in the far south of eastern Australia, and in the far south west of Western Australia by 2050.     

Ecophysiological variation across a forest‐ecotone gradient produces divergent climate change vulnerability within species

Climate change related risks and impacts on ectotherms will be mediated by habitats and their influence on local thermal environments. While many studies have documented morphological and genetic aspects of niche divergence across habitats, few have examined thermal performance across such gradients and directly linked this variation to contemporary climate change impacts. In this study, we quantified variation in thermal performance across a gradient from forest to gallery forest‐savanna mosaic in Cameroon for a skink species (Trachylepis affinis) known to be diverging genetically and morphologically across that habitat gradient. Based on these results, we then applied a mechanistic modelling approach (NicheMapR) to project changes in potential activity, as constrained by thermal performance, in response to climate change. As a complimentary approach, we also compared mechanistic projections with climate‐driven changes in habitat suitability based on species distribution models of forest and ecotone skinks. We found that ecotone skinks may benefit from warming and experience increased activity while forest skinks will likely face a drastic decrease in thermal suitability across the forest zone. Species distribution models projected that thermal suitability for forest skinks in coastal forests would decline but in other parts of the forest zone skinks are projected to experience increased thermal suitability. The results here highlight the utility of mechanistic approaches in revealing and understanding patterns of climate change vulnerability which may not be detected with species distribution models alone. This study also emphasizes the importance of intra‐specific physiological variation, and habitat‐specific thermal performance relationships in particular, in determining warming responses.     

Species Distribution Modelling predicts habitat suitability and reduction of suitable habitat under future climatic scenario for Sclerophrys perreti: A critically endangered Nigerian endemic toad

Sclerophrys perreti is a critically endangered Nigerian native frog currently imperilled by human activities. A better understanding of its potential distribution and habitat suitability will aid in conservation; however, such knowledge is limited for S. perreti. Herein, we used a species distribution model (SDM) approach with all known occurrence data (n = 22) from our field surveys and primary literature, and environmental variable predictors (19 bioclimatic variables, elevation and land cover) to elucidate habitat suitability and impact of climate change on this species. The SDM showed that temperature and precipitation were the predictors of habitat suitability for S. perreti with precipitation seasonality as the strongest predictor of habitat suitability. The following variable also had a significant effect on habitat suitability: temperature seasonality, temperature annual range, precipitation of driest month, mean temperature of wettest quarter and isothermality. The model predicted current suitable habitat for S. perreti covering an area of 1,115 km². However, this habitat is predicted to experience 60% reduction by 2050 owing to changes in temperature and precipitation. SDM also showed that suitable habitat exists in south-eastern range of the inselberg with predicted low impact of climate change compared to other ranges. Therefore, this study recommends improved conservation measures through collaborations and stakeholder's meeting with local farmers for the management and protection of S. perreti.     

Climate warming alters effects of management on population viability of threatened species: results from a 30‐year experimental study on a rare orchid

Climate change is expected to influence the viability of populations both directly and indirectly, via species interactions. The effects of large‐scale climate change are also likely to interact with local habitat conditions. Management actions designed to preserve threatened species therefore need to adapt both to the prevailing climate and local conditions. Yet, few studies have separated the direct and indirect effects of climatic variables on the viability of local populations and discussed the implications for optimal management. We used 30 years of demographic data to estimate the simultaneous effects of management practice and among‐year variation in four climatic variables on individual survival, growth and fecundity in one coastal and one inland population of the perennial orchid Dactylorhiza lapponica in Norway. Current management, mowing, is expected to reduce competitive interactions. Statistical models of how climate and management practice influenced vital rates were incorporated into matrix population models to quantify effects on population growth rate. Effects of climate differed between mown and control plots in both populations. In particular, population growth rate increased more strongly with summer temperature in mown plots than in control plots. Population growth rate declined with spring temperature in the inland population, and with precipitation in the coastal population, and the decline was stronger in control plots in both populations. These results illustrate that both direct and indirect effects of climate change are important for population viability and that net effects depend both on local abiotic conditions and on biotic conditions in terms of management practice and intensity of competition. The results also show that effects of management practices influencing competitive interactions can strongly depend on climatic factors. We conclude that interactions between climate and management should be considered to reliably predict future population viability and optimize conservation actions.     

Ancient woodland indicators signal the climate change risk for dispersal-limited species

The climate change risk to biodiversity operates alongside a range of anthropogenic pressures. These include habitat loss and fragmentation, which may prevent species from migrating between isolated habitat patches in order to track their suitable climate space. Predictive modelling has advanced in scope and complexity to integrate: (i) projected shifts in climate suitability, with (ii) spatial patterns of landscape habitat quality and rates of dispersal. This improved ecological realism is suited to data-rich model species, though its broader generalisation comes with accumulated uncertainties, e.g. incomplete knowledge of species response to variable habitat quality, parameterisation of dispersal kernels etc. This study adopts ancient woodland indicator species (lichen epiphytes) as a guild that couples relative simplicity with biological rigour. Subjectively-assigned indicator species were statistically tested against a binary habitat map of woodlands of known continuity (>250 yr), and bioclimatic models were used to demonstrate trends in their increased/decreased environmental suitability under conditions of ‘no dispersal’. Given the expectation of rapid climate change on ecological time-scales, no dispersal for ancient woodland indicators becomes a plausible assumption. The risk to ancient woodland indicators is spatially structured (greater in a relative continental compared to an oceanic climatic zone), though regional differences are weakened by significant variation (within regions) in woodland extent. As a corollary, ancient woodland indicators that are sensitive to projected climate change scenarios may be excellent targets for monitoring climate change impacts for biodiversity at a site-scale, including the outcome of strategic habitat management (climate change adaptation) designed to offset risk for dispersal-limited species.     

Range position and climate sensitivity: The structure of among‐population demographic responses to climatic variation

Species’ distributions will respond to climate change based on the relationship between local demographic processes and climate and how this relationship varies based on range position. A rarely tested demographic prediction is that populations at the extremes of a species’ climate envelope (e.g., populations in areas with the highest mean annual temperature) will be most sensitive to local shifts in climate (i.e., warming). We tested this prediction using a dynamic species distribution model linking demographic rates to variation in temperature and precipitation for wood frogs (Lithobates sylvaticus) in North America. Using long‐term monitoring data from 746 populations in 27 study areas, we determined how climatic variation affected population growth rates and how these relationships varied with respect to long‐term climate. Some models supported the predicted pattern, with negative effects of extreme summer temperatures in hotter areas and positive effects on recruitment for summer water availability in drier areas. We also found evidence of interacting temperature and precipitation influencing population size, such as extreme heat having less of a negative effect in wetter areas. Other results were contrary to predictions, such as positive effects of summer water availability in wetter parts of the range and positive responses to winter warming especially in milder areas. In general, we found wood frogs were more sensitive to changes in temperature or temperature interacting with precipitation than to changes in precipitation alone. Our results suggest that sensitivity to changes in climate cannot be predicted simply by knowing locations within the species’ climate envelope. Many climate processes did not affect population growth rates in the predicted direction based on range position. Processes such as species‐interactions, local adaptation, and interactions with the physical landscape likely affect the responses we observed. Our work highlights the need to measure demographic responses to changing climate.     

Habitat suitability modelling for species at risk is sensitive to algorithm and scale: A case study of Blanding's turtle,Emydoidea blandingii,in Ontario,Canada

Species distribution modelling (SDM) can help conservation by providing information on the ecological requirements of species at risk. We developed habitat suitability models at multiple spatial scales for a threatened freshwater turtle, Emydoidea blandingii, in Ontario as a case study. We also explored the effect of background data selection and modelling algorithm selection on habitat suitability predictions. We used sighting records, high-resolution land cover data (25 m), and two SDM techniques: boosted regression trees; and maximum entropy modelling. The area under the receiver characteristic operating curve (AUC) for habitat suitability models tested on independent data ranged from 0.878 to 0.912 when using random background and from 0.727 to 0.741 with target-group background. E. blandingii habitat suitability was best predicted by air temperature, wetland area, open water area, road density, and cropland area. Habitat suitability increased with increasing air temperature and wetland area, and decreased with increasing cropland area. Low road density and open water increased habitat suitability, while high levels of either variable decreased habitat suitability. Robust habitat suitability maps for species at risk require using a multi-scale and multi-algorithm approach. If well used, SDM can offer insight on the habitat requirements of species at risk and help guide the development of management plans. Our results suggest that E. blandingii management plans should promote the protection of terrestrial habitat surrounding residential wetlands, halt the building of roads within and adjacent to currently occupied habitat, and identify movement corridors for isolated populations.     

On the use of climate covariates in aquatic species distribution models: are we at risk of throwing out the baby with the bath water?

Species distribution models (SDMs) in river ecosystems can incorporate climate information by using air temperature and precipitation as surrogate measures of instream conditions or by using independent models of water temperature and hydrology to link climate to instream habitat. The latter approach is preferable but constrained by the logistical burden of developing water temperature and hydrology models. We therefore assessed whether regional scale, freshwater SDM predictions are fundamentally different when climate data versus instream temperature and hydrology are used as covariates. Maximum entropy (MaxEnt) SDMs were built for 15 freshwater fishes using one of two covariate sets: 1) air temperature and precipitation (climate variables) in combination with physical habitat variables; or 2) water temperature, hydrology (instream variables) and physical habitat. Three procedures were then used to compare results from climate vs instream models. First, equivalence tests assessed average pairwise differences (site‐specific comparisons throughout each species’ range) among climate and instream models. Second, ‘congruence’ tests determined how often the same stream segments were assigned high habitat suitability by climate and instream models. Third, Schoener's D and Warren's I niche overlap statistics quantified range‐wide similarity in predicted habitat suitability from climate vs instream models. Equivalence tests revealed small, pairwise differences in habitat suitability between climate and instream models (mean pairwise differences in MaxEnt raw scores for all species < 3 × 10^–4). Congruence tests showed a strong tendency for climate and instream models to predict high habitat suitability at the same stream segments (median congruence = 68%). D and I statistics reflected a high margin of overlap among climate and instream models (median D = 0.78, median I = 0.96). Overall, we found little support for the hypothesis that SDM predictions are fundamentally different when climate versus instream covariates are used to model fish species’ distributions at the scale of the Columbia Basin.     

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.     

Using a species distribution model to guide NSW surveys of the long‐footed potoroo (Potorous longipes)

Knowledge of threatened species’ distributions is essential for effective conservation decision‐making. Species distribution models (SDMs) are widely used to map species’ geographic ranges, identify new areas of suitable habitat and guide field surveys. In New South Wales (NSW), Australia, there are grave doubts about whether populations of the critically endangered long‐footed potoroo (Potorous longipes) remain extant, and identification of occupied sites is a high priority for its conservation. We used an SDM (Maxent) to identify regions in NSW that may have suitable habitat for the potoroo. The SDM was built with seven climate layers and had strong predictive performance (cross‐validated AUC = 0.94). We then combined this information on habitat suitability with vegetation and topography, to identify 58 survey sites across NSW. From April 2016 to May 2017, we undertook six field trips deploying six to eight cameras at each site for 52–63 days, resulting in 25 120 camera trap nights. A total of 215 759 images captured 43 native and feral animal species, but no long‐footed potoroos. Following the survey, newly available, independent presence and absence data were used to validate our model. A Kruskal–Wallis H test indicated that habitat suitability values were significantly higher at presence locations than absence locations (H = 58.66, d.f. = 1, P < 0.001). Finally, we refitted the Maxent model with the new data and identified additional regions that future surveys could explore. We conclude, however, that if the long‐footed potoroo remains extant in NSW, it is extremely rare.     

Intraspecific variation buffers projected climate change impacts on Pinus contorta

Species distribution modeling (SDM) is an important tool to assess the impact of global environmental change. Many species exhibit ecologically relevant intraspecific variation, and few studies have analyzed its relevance for SDM. Here, we compared three SDM techniques for the highly variable species Pinus contorta. First, applying a conventional SDM approach, we used MaxEnt to model the subject as a single species (species model), based on presence–absence observations. Second, we used MaxEnt to model each of the three most prevalent subspecies independently and combined their projected distributions (subspecies model). Finally, we used a universal growth transfer function (UTF), an approach to incorporate intraspecific variation utilizing provenance trial tree growth data. Different model approaches performed similarly when predicting current distributions. MaxEnt model discrimination was greater (AUC – species model: 0.94, subspecies model: 0.95, UTF: 0.89), but the UTF was better calibrated (slope and bias – species model: 1.31 and −0.58, subspecies model: 1.44 and −0.43, UTF: 1.01 and 0.04, respectively). Contrastingly, for future climatic conditions, projections of lodgepole pine habitat suitability diverged. In particular, when the species'' intraspecific variability was acknowledged, the species was projected to better tolerate climatic change as related to suitable habitat without migration (subspecies model: 26% habitat loss or UTF: 24% habitat loss vs. species model: 60% habitat loss), and given unlimited migration may increase amount of suitable habitat (subspecies model: 8% habitat gain or UTF: 12% habitat gain vs. species model: 51% habitat loss) in the climatic period 2070–2100 (SRES A2 scenario, HADCM3). We conclude that models derived from within-species data produce different and better projections, and coincide with ecological theory. Furthermore, we conclude that intraspecific variation may buffer against adverse effects of climate change. A key future research challenge lies in assessing the extent to which species can utilize intraspecific variation under rapid environmental change.