Intra‐specific variability and plasticity influence potential tree species distributions under climate change

Aim To assess the effect of local adaptation and phenotypic plasticity on the potential distribution of species under future climate changes. Trees may be adapted to specific climatic conditions; however, species range predictions have classically been assessed by species distribution models (SDMs) that do not account for intra‐specific genetic variability and phenotypic plasticity, because SDMs rely on the assumption that species respond homogeneously to climate change across their range, i.e. a species is equally adapted throughout its range, and all species are equally plastic. These assumptions could cause SDMs to exaggerate or underestimate species at risk under future climate change. Location The Iberian Peninsula. Methods Species distributions are predicted by integrating experimental data and modelling techniques. We incorporate plasticity and local adaptation into a SDM by calibrating models of tree survivorship with adaptive traits in provenance trials. Phenotypic plasticity was incorporated by calibrating our model with a climatic index that provides a measure of the differences between sites and provenances. Results We present a new modelling approach that is easy to implement and makes use of existing tree provenance trials to predict species distribution models under global warming. Our results indicate that the incorporation of intra‐population genetic diversity and phenotypic plasticity in SDMs significantly altered their outcome. In comparing species range predictions, the decrease in area occupancy under global warming conditions is smaller when considering our survival–adaptation model than that predicted by a ‘classical SDM’ calibrated with presence–absence data. These differences in survivorship are due to both local adaptation and plasticity. Differences due to the use of experimental data in the model calibration are also expressed in our results: we incorporate a null model that uses survival data from all provenances together. This model always predicts less reduction in area occupancy for both species than the SDM calibrated with presence–absence. Main conclusions We reaffirm the importance of considering adaptive traits when predicting species distributions and avoiding the use of occurrence data as a predictive variable. In light of these recommendations, we advise that existing predictions of future species distributions and their component populations must be reconsidered.     

Mutualism influences species distribution predictions for a bromeliad‐breeding anuran under climate change

Ecological niche models, or species distribution models, have been widely used to identify potentially suitable areas for species in future climate change scenarios. However, there are inherent errors to these models due to their inability to evaluate species occurrence influenced by non‐climatic factors. With the intuit to improve the modelling predictions for a bromeliad‐breeding treefrog (Phyllodytes melanomystax, Hylidae), we investigate how the climatic suitability of bromeliads influences the distribution model for the treefrog in the context of baseline and 2050 climate change scenarios. We used point occurrence data on the frog and the bromeliad (Vriesea procera, Bromeliaceae) to generate their predicted distributions based on baseline and 2050 climates. Using a consensus of five algorithms, we compared the accuracy of the models and the geographic predictions for the frog generated from two modelling procedures: (i) a climate‐only model for P. melanomystax and V. procera; and (ii) a climate‐biotic model for P. melanomystax, in which the climatic suitability of the bromeliad was jointly considered with the climatic variables. Both modelling approaches generated strong and similar predictive power for P. melanomystax, yet climate‐biotic modelling generated more concise predictions, particularly for the year 2050. Specifically, because the predicted area of the bromeliad overlaps with the predictions for the treefrog in the baseline climate, both modelling approaches produce reasonable similar predicted areas for the anuran. Alternatively, due to the predicted loss of northern climatically suitable areas for the bromeliad by 2050, only the climate‐biotic models provide evidence that northern populations of P. melanomystax will likely be negatively affected by 2050.     

Unpacking the mechanisms captured by a correlative species distribution model to improve predictions of climate refugia

Climate refugia are regions that animals can retreat to, persist in and potentially then expand from under changing environmental conditions. Most forecasts of climate change refugia for species are based on correlative species distribution models (SDMs) using long‐term climate averages, projected to future climate scenarios. Limitations of such methods include the need to extrapolate into novel environments and uncertainty regarding the extent to which proximate variables included in the model capture processes driving distribution limits (and thus can be assumed to provide reliable predictions under new conditions). These limitations are well documented; however, their impact on the quality of climate refugia predictions is difficult to quantify. Here, we develop a detailed bioenergetics model for the koala. It indicates that range limits are driven by heat‐induced water stress, with the timing of rainfall and heat waves limiting the koala in the warmer parts of its range. We compare refugia predictions from the bioenergetics model with predictions from a suite of competing correlative SDMs under a range of future climate scenarios. SDMs were fitted using combinations of long‐term climate and weather extremes variables, to test how well each set of predictions captures the knowledge embedded in the bioenergetics model. Correlative models produced broadly similar predictions to the bioenergetics model across much of the species' current range – with SDMs that included weather extremes showing highest congruence. However, predictions in some regions diverged significantly when projecting to future climates due to the breakdown in correlation between climate variables. We provide unique insight into the mechanisms driving koala distribution and illustrate the importance of subtle relationships between the timing of weather events, particularly rain relative to hot‐spells, in driving species–climate relationships and distributions. By unpacking the mechanisms captured by correlative SDMs, we can increase our certainty in forecasts of climate change impacts on species.     

Choice of predictor variables as a source of uncertainty in continental‐scale species distribution modelling under climate change

Aim Species distribution modelling is commonly used to guide future conservation policies in the light of potential climate change. However, arbitrary decisions during the model‐building process can affect predictions and contribute to uncertainty about where suitable climate space will exist. For many species, the key climatic factors limiting distributions are unknown. This paper assesses the uncertainty generated by using different climate predictor variable sets for modelling the impacts of climate change. Location Europe, 10° W to 50° E and 30° N to 60° N. Methods Using 1453 presence pixels at 30 arcsec resolution for the great bustard (Otis tarda), predictions of future distribution were made based on two emissions scenarios, three general climate models and 26 sets of predictor variables. Twenty‐six current models were created, and 156 for both 2050 and 2080. Map comparison techniques were used to compare predictions in terms of the quantity and the location of presences (map comparison kappa, MCK) and using a range change index (RCI). Generalized linear models (GLMs) were used to partition explained deviance in MCK and RCI among sources of uncertainty. Results The 26 different variable sets achieved high values of AUC (area under the receiver operating characteristic curve) and yet introduced substantial variation into maps of current distribution. Differences between maps were even greater when distributions were projected into the future. Some 64–78% of the variation between future maps was attributable to choice of predictor variable set alone. Choice of general climate model and emissions scenario contributed a maximum of 15% variation and their order of importance differed for MCK and RCI. Main conclusions Generalized variable sets produce an unmanageable level of uncertainty in species distribution models which cannot be ignored. The use of sound ecological theory and statistical methods to check predictor variables can reduce this uncertainty, but our knowledge of species may be too limited to make more than arbitrary choices. When all sources of modelling uncertainty are considered together, it is doubtful whether ensemble methods offer an adequate solution. Future studies should explicitly acknowledge uncertainty due to arbitrary choices in the model‐building process and develop ways to convey the results to decision‐makers.     

From climate change predictions to actions – conserving vulnerable animal groups in hotspots at a regional scale

Current climate change is a major threat to biodiversity. Species unable to adapt or move will face local or global extinction and this is more likely to happen to species with narrow climatic and habitat requirements and limited dispersal abilities, such as amphibians and reptiles. Biodiversity losses are likely to be greatest in global biodiversity hotspots where climate change is fast, such as the Iberian Peninsula. Here we assess the impact of climate change on 37 endemic and nearly endemic herptiles of the Iberian Peninsula by predicting species distributions for three different times into the future (2020, 2050 and 2080) using an ensemble of bioclimatic models and different combinations of species dispersal ability, emission levels and global circulation models. Our results show that species with Atlantic affinities that occur mainly in the North‐western Iberian Peninsula have severely reduced future distributions. Up to 13 species may lose their entire potential distribution by 2080. Furthermore, our analysis indicates that the most critical period for the majority of these species will be the next decade. While there is considerable variability between the scenarios, we believe that our results provide a robust relative evaluation of climate change impacts among different species. Future evaluation of the vulnerability of individual species to climate change should account for their adaptive capacity to climate change, including factors such as physiological climate tolerance, geographical range size, local abundance, life cycle, behavioural and phenological adaptability, evolutionary potential and dispersal ability.     

Modeling nonbreeding distributions of shorebirds and waterfowl in response to climate change

To identify areas on the landscape that may contribute to a robust network of conservation areas, we modeled the probabilities of occurrence of several en route migratory shorebirds and wintering waterfowl in the southern Great Plains of North America, including responses to changing climate. We predominantly used data from the eB ird citizen‐science project to model probabilities of occurrence relative to land‐use patterns, spatial distribution of wetlands, and climate. We projected models to potential future climate conditions using five representative general circulation models of the Coupled Model Intercomparison Project 5 (CMIP 5). We used Random Forests to model probabilities of occurrence and compared the time periods 1981–2010 (hindcast) and 2041–2070 (forecast) in “model space.” Projected changes in shorebird probabilities of occurrence varied with species‐specific general distribution pattern, migration distance, and spatial extent. Species using the western and northern portion of the study area exhibited the greatest likelihoods of decline, whereas species with more easterly occurrences, mostly long‐distance migrants, had the greatest projected increases in probability of occurrence. At an ecoregional extent, differences in probabilities of shorebird occurrence ranged from ?0.015 to 0.045 when averaged across climate models, with the largest increases occurring early in migration. Spatial shifts are predicted for several shorebird species. Probabilities of occurrence of wintering Mallards and Northern Pintail are predicted to increase by 0.046 and 0.061, respectively, with northward shifts projected for both species. When incorporated into partner land management decision tools, results at ecoregional extents can be used to identify wetland complexes with the greatest potential to support birds in the nonbreeding season under a wide range of future climate scenarios.     

Biotic interactions influence the projected distribution of a specialist mammal under climate change

Aim

To measure the effects of including biotic interactions on climate‐based species distribution models (SDMs) used to predict distribution shifts under climate change. We evaluated the performance of distribution models for an endangered marsupial, the northern bettong (Bettongia tropica), comparing models that used only climate variables with models that also took into account biotic interactions.

Location

North‐east Queensland, Australia.

Methods

We developed separate climate‐based distribution models for the northern bettong, its two main resources and a competitor species. We then constructed models for the northern bettong by including climate suitability estimates for the resources and competitor as additional predictor variables to make climate + resource and climate + resource + competition models. We projected these models onto seven future climate scenarios and compared predictions of northern bettong distribution made by these differently structured models, using a ‘global’ metric, the I similarity statistic, to measure overlap in distribution and a ‘local’ metric to identify where predictions differed significantly.

Results

Inclusion of food resource biotic interactions improved model performance. Over moderate climate changes, up to 3.0 °C of warming, the climate‐only model for the northern bettong gave similar predictions of distribution to the more complex models including interactions, with differences only at the margins of predicted distributions. For climate changes beyond 3.0 °C, model predictions diverged significantly. The interactive model predicted less contraction of distribution than the simpler climate‐only model.

Main conclusions

Distribution models that account for interactions with other species, in particular direct resources, improve model predictions in the present‐day climate. For larger climate changes, shifts in distribution of interacting species cause predictions of interactive models to diverge from climate‐only models. Incorporating interactions with other species in SDMs may be needed for long‐term prediction of changes in distribution of species under climate change, particularly for specialized species strongly dependent on a small number of biotic interactions.     

Detecting nonlinear response of spring phenology to climate change by Bayesian analysis

The impact of climate change on the advancement of plant phenological events has been heavily studied in the last decade. Although the majority of spring plant phenological events have been trending earlier, this is not universally true. Recent work has suggested that species that are not advancing in their spring phenological behavior are responding more to lack of winter chill than increased spring heat. One way to test this hypothesis is by evaluating the behavior of a species known to have a moderate to high chilling requirement and examining how it is responding to increased warming. This study used a 60‐year data set for timing of leaf‐out and male flowering of walnut (Juglans regia) cultivar ‘Payne’ to examine this issue. The spring phenological behavior of ‘Payne’ walnut differed depending on bud type. The vegetative buds, which have a higher chilling requirement, trended toward earlier leaf‐out until about 1994, when they shifted to later leaf‐out. The date of male bud pollen shedding advanced over the course of the whole record. Our findings suggest that many species which have exhibited earlier bud break are responding to warmer spring temperatures, but may shift into responding more to winter temperatures (lack of adequate chilling) as warming continues.     

Predicting shifts in parasite distribution with climate change: a multitrophic level approach

Climate change likely will lead to increasingly favourable environmental conditions for many parasites. However, predictions regarding parasitism's impacts often fail to account for the likely variability in host distribution and how this may alter parasite occurrence. Here, we investigate potential distributional shifts in the meningeal worm, Parelaphostrongylosis tenuis, a protostrongylid nematode commonly found in white‐tailed deer in North America, whose life cycle also involves a free‐living stage and a gastropod intermediate host. We modelled the distribution of the hosts and free‐living larva as a complete assemblage to assess whether a complex trophic system will lead to an overall increase in parasite distribution with climate change, or whether divergent environmental niches may promote an ecological mismatch. Using an ensemble approach to climate modelling under two different carbon emission scenarios, we show that whereas the overall trend is for an increase in niche breadth for each species, mismatches arise in habitat suitability of the free‐living larva vs. the definitive and intermediate hosts. By incorporating these projected mismatches into a combined model, we project a shift in parasite distribution accounting for all steps in the transmission cycle, and identify that overall habitat suitability of the parasite will decline in the Great Plains and southeastern USA, but will increase in the Boreal Forest ecoregion, particularly in Alberta. These results have important implications for wildlife conservation and management due to the known pathogenicity of parelaphostrongylosis to alternate hosts including moose, caribou and elk. Our results suggest that disease risk forecasts which fail to consider biotic interactions may be overly simplistic, and that accounting for each of the parasite's life stages is key to refining predicted responses to climate change.     

Choice of baseline climate data impacts projected species' responses to climate change

Climate data created from historic climate observations are integral to most assessments of potential climate change impacts, and frequently comprise the baseline period used to infer species‐climate relationships. They are often also central to downscaling coarse resolution climate simulations from General Circulation Models (GCMs) to project future climate scenarios at ecologically relevant spatial scales. Uncertainty in these baseline data can be large, particularly where weather observations are sparse and climate dynamics are complex (e.g. over mountainous or coastal regions). Yet, importantly, this uncertainty is almost universally overlooked when assessing potential responses of species to climate change. Here, we assessed the importance of historic baseline climate uncertainty for projections of species' responses to future climate change. We built species distribution models (SDMs) for 895 African bird species of conservation concern, using six different climate baselines. We projected these models to two future periods (2040–2069, 2070–2099), using downscaled climate projections, and calculated species turnover and changes in species‐specific climate suitability. We found that the choice of baseline climate data constituted an important source of uncertainty in projections of both species turnover and species‐specific climate suitability, often comparable with, or more important than, uncertainty arising from the choice of GCM. Importantly, the relative contribution of these factors to projection uncertainty varied spatially. Moreover, when projecting SDMs to sites of biodiversity importance (Important Bird and Biodiversity Areas), these uncertainties altered site‐level impacts, which could affect conservation prioritization. Our results highlight that projections of species' responses to climate change are sensitive to uncertainty in the baseline climatology. We recommend that this should be considered routinely in such analyses.     

Probabilistic accounting of uncertainty in forecasts of species distributions under climate change

Forecasts of species distributions under future climates are inherently uncertain, but there have been few attempts to describe this uncertainty comprehensively in a probabilistic manner. We developed a Monte Carlo approach that accounts for uncertainty within generalized linear regression models (parameter uncertainty and residual error), uncertainty among competing models (model uncertainty), and uncertainty in future climate conditions (climate uncertainty) to produce site‐specific frequency distributions of occurrence probabilities across a species' range. We illustrated the method by forecasting suitable habitat for bull trout (Salvelinus confluentus) in the Interior Columbia River Basin, USA, under recent and projected 2040s and 2080s climate conditions. The 95% interval of total suitable habitat under recent conditions was estimated at 30.1–42.5 thousand km; this was predicted to decline to 0.5–7.9 thousand km by the 2080s. Projections for the 2080s showed that the great majority of stream segments would be unsuitable with high certainty, regardless of the climate data set or bull trout model employed. The largest contributor to uncertainty in total suitable habitat was climate uncertainty, followed by parameter uncertainty and model uncertainty. Our approach makes it possible to calculate a full distribution of possible outcomes for a species, and permits ready graphical display of uncertainty for individual locations and of total habitat.     

Morphological and dietary responses of chipmunks to a century of climate change

Predicting how individual taxa will respond to climatic change is challenging, in part because the impacts of environmental conditions can vary markedly, even among closely related species. Studies of chipmunks (Tamias spp.) in Yosemite National Park provide an important opportunity to explore the reasons for this variation in response. While the alpine chipmunk (T. alpinus) has undergone a significant elevational range contraction over the past century, the congeneric and partially sympatric lodgepole chipmunk (T. speciosus) has not experienced an elevational range shift during this period. As a first step toward identifying the factors underlying this difference in response, we examined evidence for dietary changes and changes in cranial morphology in these species over the past century. Stable isotope analyses of fur samples from modern and historical museum specimens of these species collected at the same localities indicated that signatures of dietary change were more pronounced in T. alpinus, although diet breadth did not differ consistently between the study species. Morphometric analyses of crania from these specimens revealed significant changes in cranial shape for T. alpinus, with less pronounced changes in shape for T. speciosus; evidence of selection on skull morphology was detected for T. alpinus, but not for T. speciosus. These results are consistent with growing evidence that T. alpinus is generally more responsive to environmental change than T. speciosus, but emphasize the complex and often geographically variable nature of such responses. Accordingly, future studies that make use of the taxonomically and spatially integrative approach employed here may prove particularly informative regarding relationships between environmental conditions, range changes, and patterns of phenotypic variation.     

Evaluating within‐population variability in behavior and demography for the adaptive potential of a dispersal‐limited species to climate change

Multiple pathways exist for species to respond to changing climates. However, responses of dispersal‐limited species will be more strongly tied to ability to adapt within existing populations as rates of environmental change will likely exceed movement rates. Here, we assess adaptive capacity in Plethodon cinereus, a dispersal‐limited woodland salamander. We quantify plasticity in behavior and variation in demography to observed variation in environmental variables over a 5‐year period. We found strong evidence that temperature and rainfall influence P. cinereus surface presence, indicating changes in climate are likely to affect seasonal activity patterns. We also found that warmer summer temperatures reduced individual growth rates into the autumn, which is likely to have negative demographic consequences. Reduced growth rates may delay reproductive maturity and lead to reductions in size‐specific fecundity, potentially reducing population‐level persistence. To better understand within‐population variability in responses, we examined differences between two common color morphs. Previous evidence suggests that the color polymorphism may be linked to physiological differences in heat and moisture tolerance. We found only moderate support for morph‐specific differences for the relationship between individual growth and temperature. Measuring environmental sensitivity to climatic variability is the first step in predicting species' responses to climate change. Our results suggest phenological shifts and changes in growth rates are likely responses under scenarios where further warming occurs, and we discuss possible adaptive strategies for resulting selective pressures.     

Early response of the platypus to climate warming

Combining a climatic envelope modelling technique with more than two centuries (1800–2009) of distribution records has revealed the effects of a changing climate on the egg‐laying monotreme, the platypus, Ornithorhynchus anatinus. We show that the main factor associated with platypus occurrence switched from aquatic habitat availability (estimated by rainfall) to thermal tolerances (estimated by annual maximum temperature) in the 1960s. This correlates directly with the change in the annual maximum temperature anomaly from cooler to warmer conditions in southeastern Australia. Modelling of platypus habitat under emission scenarios (A1B, A2, B1 and B2) revealed large decreases (>30%) in thermally suitable habitat by 2070. This reduction, compounded by increasing demands for water for agriculture and potable use, suggests that there is real cause for concern over the future status of this species, and highlights the need for restoration of thermal refugia within the platypus’ modelled range.     

Rapid evolution and behavioral plasticity following introduction to an environment with reduced predation risk

Adaptive behavioral plasticity can play a beneficial role when a population becomes established in a novel environment if environmental cues allow the expression of appropriate behavior. Further, plasticity itself can evolve over time in a new environment causing changes in the way or degree to which animals respond to environmental cues. Colonization events provide an opportunity to investigate such relationships between behavioral plasticity and adaptation to new environments. Here, we investigated the evolution of behavior and behavioral plasticity during colonization of a new environment, by testing if female mate‐choice behavior diverged in Trinidadian guppies 2–3 years (~6–9 generations) after being introduced to four locations with reduced predation risk. We collected wild‐caught fish from the source and introduced populations, and we reared out second‐generation females in the laboratory with and without predator cues to examine their plastic responses to a bright and dull male. We found introduced females were less responsive to males when reared without predator cues, but both introduced and source females were similarly responsive when reared with predator cues. Thus, the parallel evolution of behavior across multiple populations in the low‐predation environment was only observed in the treatment mimicking the introduction environment. Such results are consistent with theory predicting that the evolution of plasticity is a by‐product of differential selection across environments.