Selecting from correlated climate variables: a major source of uncertainty for predicting species distributions under climate change

Correlative species distribution models are frequently used to predict species’ range shifts under climate change. However, climate variables often show high collinearity and most statistical approaches require the selection of one among strongly correlated variables. When causal relationships between species presence and climate parameters are unknown, variable selection is often arbitrary, or based on predictive performance under current conditions. While this should only marginally affect current range predictions, future distributions may vary considerably when climate parameters do not change in concert. We investigated this source of uncertainty using four highly correlated climate variables together with a constant set of landscape variables in order to predict current (2010) and future (2050) distributions of four mountain bird species in central Europe. Simulating different parameterization decisions, we generated a) four models including each of the climate variables singly, b) a model taking advantage of all variables simultaneously and c) an un‐weighted average of the predictions of a). We compared model accuracy under current conditions, predicted distributions under four scenarios of climate change, and – for one species – evaluated back‐projections using historical occurrence data. Although current and future variable‐correlations remained constant, and the models’ accuracy under contemporary conditions did not differ, future range predictions varied considerably in all climate change scenarios. Averaged models and models containing all climate variables simultaneously produced intermediate predictions; the latter, however, performed best in back‐projections. This pattern, consistent across different modelling methods, indicates a benefit from including multiple climate predictors in ambiguous situations. Variable selection proved to be an important source of uncertainty for future range predictions, difficult to control using contemporary information. Small, but diverging changes of climate variables, masked by constant overall correlation patterns, can cause substantial differences between future range predictions which need to be accounted for, particularly when outcomes are intended for conservation decisions.     

Uncertainty in assessing the impacts of global change with coupled dynamic species distribution and population models

Concern over rapid global changes and the potential for interactions among multiple threats are prompting scientists to combine multiple modelling approaches to understand impacts on biodiversity. A relatively recent development is the combination of species distribution models, land‐use change predictions, and dynamic population models to predict the relative and combined impacts of climate change, land‐use change, and altered disturbance regimes on species' extinction risk. Each modelling component introduces its own source of uncertainty through different parameters and assumptions, which, when combined, can result in compounded uncertainty that can have major implications for management. Although some uncertainty analyses have been conducted separately on various model components – such as climate predictions, species distribution models, land‐use change predictions, and population models – a unified sensitivity analysis comparing various sources of uncertainty in combined modelling approaches is needed to identify the most influential and problematic assumptions. We estimated the sensitivities of long‐run population predictions to different ecological assumptions and parameter settings for a rare and endangered annual plant species (Acanthomintha ilicifolia, or San Diego thornmint). Uncertainty about habitat suitability predictions, due to the choice of species distribution model, contributed most to variation in predictions about long‐run populations.     

A spatio‐temporal framework for efficient inventories of natural resources: A case study with submersed macrophytes

Questions: Does a reduced nutrient load in open water increase species richness and the importance of regional and local site characteristics for species abundance and spatial distribution? Can we build lake‐specific models of macrophyte abundance and distribution based on site characteristics in order to prepare a cost‐efficient framework for future surveys? Location: Lake Constance, 47°39′N, 9°18′E. Methods: Generalized additive models (GAMs) were used to predict the potential distributions of eight species and overall species richness. Submersed macrophyte distribution in 1993 was compared with corresponding data from 1978, when eutrophication was at its maximum. Results: Spatial predictions for eight species and overall species richness were relatively accurate and independent of water chemistry. Depth was confirmed as a main predictor of species distribution, while effective fetch distance was retained in many models. Mineralogical variables of sediment composition represent allogenic and autogenic sediment sources and their east‐west gradient in Lake Constance corresponded to east‐west gradients of species distribution and richness. GAMs appeared more efficient than generalized linear models (GLMs) for modelling species responses to environmental gradients. Conclusions: Reduced trophic status increases species richness and the importance of regional and local site characteristics for species abundance and distribution. Our models represent a spatio‐temporal framework for future lake monitoring purposes and allow the development of effective monitoring; this could be generalized for many ecosystem types and would be particularly efficient for large lakes such as Lake Constance.     

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.     

Climate versus weather extremes: Temporal predictor resolution matters for future rather than current regional species distribution models

Aim

Climate is considered a major driver of species distributions. Long‐term climatic means are commonly used as predictors in correlative species distribution models (SDMs). However, this coarse temporal resolution does not reflect local conditions that populations experience, such as short‐term weather extremes, which may have a strong impact on population dynamics and local distributions. We here compare the performance of climate‐ and weather‐based predictors in regional SDMs and their influence on future predictions, which are increasingly used in conservation planning.

Location

South‐western Germany.

Methods

We built different SDMs for 20 Orthoptera species based on three predictor sets at a regional scale for current and future climate scenarios. We calculated standard bioclimatic variables and yearly and seasonal sets of climate change indicating variables of weather extremes. As the impact of extreme events may be stronger for habitat specialists than for generalists, we distinguished species’ degrees of specialization. We computed linear mixed‐effects models to identify significant effects of algorithm, predictor set and specialization on model performance and calculated correlations and geographical niche overlap between spatial predictions.

Results

Current predictions were rather similar among all predictor sets, but highly variable for future climate scenarios. Bioclimatic and seasonal weather predictors performed slightly better than yearly weather predictors, though performance differences were minor. We found no evidence that specialists are more sensitive to weather extremes than generalists.

Main conclusions

For future projections of species distributions, SDM predictor selection should not solely be based on current performances and predictions. As long‐term climate and short‐term weather predictors represent different environmental drivers of a species’ distribution, we argue to interpret diverging future projections as complements. Even if similar current performances and predictions might imply their equivalency, favouring one predictor set neglects important aspects of future distributions and might mislead conservation decisions based on them.

     

Modelling and predicting fungal distribution patterns using herbarium data

Aim The main aims of this study are: (1) to test if temperature and related parameters are the primary determinants of the regional distribution of macrofungi (as is commonly recognized for plants); (2) to test if the success of modelling fungal distribution patterns depends on species and distribution characteristics; and (3) to explore the potential of using herbarium data for modelling and predicting fungal species’ distributions. Location The study area, Norway, spans 58–71° N latitude and 4–32° E longitude, and embraces extensive ecological gradients in a small area. Methods The study is based on 1020 herbarium collections of nine selected species of macrofungi and a set of 75 environmental predictor variables, all recorded in a 5 × 5‐km grid covering Norway. Primarily, generalized linear model (GLM; logistic regression) analyses were used to identify the environmental variables that best accounted for the species’ recorded distributions in Norway. Second, Maxent analyses (using variables identified by GLM) were used to produce predictive potential distribution maps for these species. Results Variables relating to temperature and radiation were most frequently included in the GLMs, and between 24.8% and 59.8% of the variation in single‐species occurrence was accounted for. The fraction of variation explained by the GLMs ranged from 41.6% to 59.8% for species with restricted distributions, and from 24.8% to 39.3% for species with widespread/scattered and intermediate distributions. The two‐step procedure of GLM followed by Maxent gave predictions with very high values for the area under the curve (0.927–0.997), and maps of potential distribution were generally credible. Main conclusions We show that temperature is a key factor governing the distribution of macrofungi in Norway, indicating that fungi may respond strongly to global warming. We confirm that modelling success depends partly on species and distribution characteristics, notably on how the distribution relates to the extent of the study area. Our study demonstrates that the combination of GLM and Maxent may be a fruitful approach for biogeography. We conclude that herbarium data improve insight into factors that control the distributions of fungi, of particular value for research on fleshy fungi (mushrooms), which have largely cryptic life cycles.     

To Be Or Not to Be? Variable selection can change the projected fate of a threatened species under future climate

Species distribution models (SDMs) are commonly used to project future changes in the geographic ranges of species, to estimate extinction rates and to plan biodiversity conservation. However, these models can produce a range of results depending on how they are parameterized, and over‐reliance on a single model may lead to overconfidence in maps of future distributions. The choice of predictor variable can have a greater influence on projected future habitat than the range of climate models used. We demonstrate this in the case of the Ptunarra Brown Butterfly, a species listed as vulnerable in Tasmania, Australia. We use the Maxent model to develop future projections for this species based on three variable sets; all 35 commonly used so‐called ‘bioclimatic’ variables, a subset of these based on expert knowledge, and a set of monthly climate variables relevant to the species’ primary activity period. We used a dynamically downscaled regional climate model based on three global climate models. Depending on the choice of variable set, the species is projected either to experience very little contraction of habitat or to come close to extinction by the end of the century due to lack of suitable climate. The different conclusions could have important consequences for conservation planning and management, including the perceived viability of habitat restoration. The output of SDMs should therefore be used to define the range of possible trajectories a species may be on, and ongoing monitoring used to inform management as changes occur.     

Biotic interactions improve prediction of boreal bird distributions at macro-scales

Aim The role of biotic interactions in influencing species distributions at macro‐scales remains poorly understood. Here we test whether predictions of distributions for four boreal owl species at two macro‐scales (10 × 10 km and 40 × 40 km grid resolutions) are improved by incorporating interactions with woodpeckers into climate envelope models. Location Finland, northern Europe. Methods Distribution data for four owl and six woodpecker species, along with data for six land cover and three climatic variables, were collated from 2861 10 × 10 km grid cells. Generalized additive models were calibrated using a 50% random sample of the species data from western Finland, and by repeating this procedure 20 times for each of the four owl species. Models were fitted using three sets of explanatory variables: (1) climate only; (2) climate and land cover; and (3) climate, land cover and two woodpecker interaction variables. Models were evaluated using three approaches: (1) examination of explained deviance; (2) four‐fold cross‐validation using the model calibration data; and (3) comparison of predicted and observed values for independent grid cells in eastern Finland. The model accuracy for approaches (2) and (3) was measured using the area under the curve of a receiver operating characteristic plot. Results At 10‐km resolution, inclusion of the distribution of woodpeckers as a predictor variable significantly improved the explanatory power, cross‐validation statistics and the predictive accuracy of the models. Inclusion of land cover led to similar improvements at 10‐km resolution, although these improvements were less apparent at 40‐km resolution for both land cover and biotic interactions. Main conclusions Predictions of species distributions at macro‐scales may be significantly improved by incorporating biotic interactions and land cover variables into models. Our results are important for models used to predict the impacts of climate change, and emphasize the need for comprehensive evaluation of the reliability of species–climate impact models.     

Model-based uncertainty in species range prediction

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

Mechanistic variables can enhance predictive models of endotherm distributions: the American pika under current,past, and future climates

How climate constrains species’ distributions through time and space is an important question in the context of conservation planning for climate change. Despite increasing awareness of the need to incorporate mechanism into species distribution models (SDMs), mechanistic modeling of endotherm distributions remains limited in this literature. Using the American pika (Ochotona princeps) as an example, we present a framework whereby mechanism can be incorporated into endotherm SDMs. Pika distribution has repeatedly been found to be constrained by warm temperatures, so we used Niche Mapper, a mechanistic heat‐balance model, to convert macroclimate data to pika‐specific surface activity time in summer across the western United States. We then explored the difference between using a macroclimate predictor (summer temperature) and using a mechanistic predictor (predicted surface activity time) in SDMs. Both approaches accurately predicted pika presences in current and past climate regimes. However, the activity models predicted 8–19% less habitat loss in response to annual temperature increases of ~3–5 °C predicted in the region by 2070, suggesting that pikas may be able to buffer some climate change effects through behavioral thermoregulation that can be captured by mechanistic modeling. Incorporating mechanism added value to the modeling by providing increased confidence in areas where different modeling approaches agreed and providing a range of outcomes in areas of disagreement. It also provided a more proximate variable relating animal distribution to climate, allowing investigations into how unique habitat characteristics and intraspecific phenotypic variation may allow pikas to exist in areas outside those predicted by generic SDMs. Only a small number of easily obtainable data are required to parameterize this mechanistic model for any endotherm, and its use can improve SDM predictions by explicitly modeling a widely applicable direct physiological effect: climate‐imposed restrictions on activity. This more complete understanding is necessary to inform climate adaptation actions, management strategies, and conservation plans.     

Simulating tick distributions over sub-Saharan Africa: the use of observed and simulated climate surfaces

Aim A broad suit of climate data sets is becoming available for use in predictive species modelling. We compare the efficacy of using interpolated climate surfaces [Center for Resource and Environmental Studies (CRES) and Climate Research Unit (CRU)] or high‐resolution model‐derived climate data [Division of Atmospheric Research limited‐area model (DARLAM)] for predictive species modelling, using tick distributions from sub‐Saharan Africa. Location The analysis is restricted to sub‐Saharan Africa. The study area was subdivided into 3000 grids cells with a resolution of 60 × 60 km. Methods Species distributions were predicted using an established multivariate climate envelope modelling approach and three very different climate data sets. The recorded variance in the climate data sets was quantified by employing omnidirectional variograms. To further compare the interpolated tick distributions that flowed from using three climate data sets, we calculated true positive (TP) predictions, false negative (FN) predictions as well as the proportional overlaps between observed and modelled tick distributions. In addition, the effect of tick data set size on the performance of the climate data sets was evaluated by performing random draws of known tick distribution records without replacement. Results The predicted distributions were consistently wider ranging than the known records when using any of the three climate data sets. However, the proportional overlap between predicted and known distributions varied as follows: for Rhipicephalus appendiculatus Neumann (Acari: Ixodidae), these were 60%, 60% and 70%; for Rhipicephalus longus Neumann (Acari: Ixodidae) 60%, 57% and 75%; for Rhipicephalus zambeziensis Walker, Norval & Corwin (Acari: Ixodidae) 57%, 51% and 62%, and for Rhipicephalus capensis Koch (Acari: Ixodidae) 70%, 60% and 60% using the CRES, CRU and DARLAM climate data sets, respectively. All data sets were sensitive to data size but DARLAM performed better when using smaller species data sets. At a 20% data subsample level, DARLAM was able to capture more than 50% of the known records and captured more than 60% of known records at higher subsample levels. Main conclusions The use of data derived from high‐resolution nested climate models (e.g. DARLAM) provided equal or even better species distribution modelling performance. As the model is dynamic and process based, the output data are available at the modelled resolution, and are not hamstrung by the sampling intensity of observed climate data sets (c. one sample per 30,000 km² for Africa). In addition, when exploring the biodiversity consequences of climate change, these modelled outputs form a more useful basis for comparison with modelled future climate scenarios.     

Validation of species–climate impact models under climate change

Increasing concern over the implications of climate change for biodiversity has led to the use of species–climate envelope models to project species extinction risk under climate‐change scenarios. However, recent studies have demonstrated significant variability in model predictions and there remains a pressing need to validate models and to reduce uncertainties. Model validation is problematic as predictions are made for events that have not yet occurred. Resubstituition and data partitioning of present‐day data sets are, therefore, commonly used to test the predictive performance of models. However, these approaches suffer from the problems of spatial and temporal autocorrelation in the calibration and validation sets. Using observed distribution shifts among 116 British breeding‐bird species over the past ～20 years, we are able to provide a first independent validation of four envelope modelling techniques under climate change. Results showed good to fair predictive performance on independent validation, although rules used to assess model performance are difficult to interpret in a decision‐planning context. We also showed that measures of performance on nonindependent data provided optimistic estimates of models' predictive ability on independent data. Artificial neural networks and generalized additive models provided generally more accurate predictions of species range shifts than generalized linear models or classification tree analysis. Data for independent model validation and replication of this study are rare and we argue that perfect validation may not in fact be conceptually possible. We also note that usefulness of models is contingent on both the questions being asked and the techniques used. Implementations of species–climate envelope models for testing hypotheses and predicting future events may prove wrong, while being potentially useful if put into appropriate context.     

Using unclassified continuous remote sensing data to improve distribution models of red-listed plant species

Remote sensing (RS) data may play an important role in the development of cost-effective means for modelling, mapping, planning and conserving biodiversity. Specifically, at the landscape scale, spatial models for the occurrences of species of conservation concern may be improved by the inclusion of RS-based predictors, to help managers to better meet different conservation challenges. In this study, we examine whether predicted distributions of 28 red-listed plant species in north-eastern Finland at the resolution of 25 ha are improved when advanced RS-variables are included as unclassified continuous predictor variables, in addition to more commonly used climate and topography variables. Using generalized additive models (GAMs), we studied whether the spatial predictions of the distribution of red-listed plant species in boreal landscapes are improved by incorporating advanced RS (normalized difference vegetation index, normalized difference soil index and Tasseled Cap transformations) information into species-environment models. Models were fitted using three different sets of explanatory variables: (1) climate-topography only; (2) remote sensing only; and (3) combined climate-topography and remote sensing variables, and evaluated by four-fold cross-validation with the area under the curve (AUC) statistics. The inclusion of RS variables improved both the explanatory power (on average 8.1 % improvement) and cross-validation performance (2.5 %) of the models. Hybrid models produced ecologically more reliable distribution maps than models using only climate-topography variables, especially for mire and shore species. In conclusion, Landsat ETM+ data integrated with climate and topographical information has the potential to improve biodiversity and rarity assessments in northern landscapes, especially in predictive studies covering extensive and remote areas.     

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.     

Uncertainty of bioclimate envelope models based on the geographical distribution of species

Aim We explored the effects of prevalence, latitudinal range and spatial autocorrelation of species distribution patterns on the accuracy of bioclimate envelope models of butterflies. Location Finland, northern Europe. Methods The data of a national butterfly atlas survey (NAFI) carried out in 1991–2003 with a resolution of 10 × 10 km were used in the analyses. Generalized additive models (GAM) were constructed, for each of 98 species, to estimate the probability of occurrence as a function of climate variables. Model performance was measured using the area under the curve (AUC) of a receiver operating characteristic (ROC) plot. Observed differences in modelling accuracy among species were related to the species’ geographical attributes using multivariate GAM. Results Accuracies of the climate–butterfly models varied from low to very high (AUC values 0.59–0.99), with a mean of 0.79. The modelling performance was related negatively to the latitudinal range and prevalence, and positively to the spatial autocorrelation of the species distribution. These three factors accounted for 75.2% of the variation in the modelling accuracy. Species at the margin of their range or with low prevalence were better predicted than widespread species, and species with clumped distributions better than scattered dispersed species. Main conclusions The results from this study indicate that species’ geographical attributes highly influence the behaviour and uncertainty of species–climate models, which should be taken into account in biogeographical modelling studies and assessments of climate change impacts.     

Bet hedging in a warming ocean: predictability of maternal environment shapes offspring size variation in marine sticklebacks

Bet hedging at reproduction is expected to evolve when mothers are exposed to unpredictable cues for future environmental conditions, whereas transgenerational plasticity (TGP) should be favoured when cues reliably predict the environment offspring will experience. Since climate predictions forecast an increase in both temperature and climate variability, both TGP and bet hedging are likely to become important strategies to mediate climate change effects. Here, the potential to produce variably sized offspring in both warming and unpredictable environments was tested by investigating whether stickleback (Gasterosteus aculeatus) mothers adjusted mean offspring size and within‐clutch variation in offspring size in response to experimental manipulation of maternal thermal environment and predictability (alternating between ambient and elevated water temperatures). Reproductive output traits of F1 females were influenced by both temperature and environmental predictability. Mothers that developed at ambient temperature (17 °C) produced larger, but fewer eggs than mothers that developed at elevated temperature (21 °C), implying selection for different‐sized offspring in different environments. Mothers in unpredictable environments had smaller mean egg sizes and tended to have greater within‐female egg size variability, especially at 21 °C, suggesting that mothers may have dynamically modified the variance in offspring size to spread the risk of incorrectly predicting future environmental conditions. Both TGP and diversification influenced F2 offspring body size. F2 offspring reared at 21 °C had larger mean body sizes if their mother developed at 21 °C, but this TGP benefit was not present for offspring of 17 °C mothers reared at 17 °C, indicating that maternal TGP will be highly relevant for ocean warming scenarios in this system. Offspring of variable environment mothers were smaller but more variable in size than offspring from constant environment mothers, particularly at 21 °C. In summary, stickleback mothers may have used both TGP and diversified bet‐hedging strategies to cope with the dual stress of ocean warming and environmental uncertainty.     

Does dispersal capacity matter for freshwater biodiversity under climate change?

相似文献   

Large-scale spatial and temporal dynamics of the vulnerable and highly mobile superb parrot

Aim To investigate the spatial and temporal dynamics of the vulnerable and highly mobile superb parrot (Polytelis swainsonii) across its range in south‐eastern mainland Australia. Location South‐eastern Australia (27°–37° S latitude and 141°–151° E longitude). Methods We used generalized additive models (GAMs) to model time‐specific bird atlas occurrence data against time‐specific plant productivity data, plus a range of environmental predictor variables. We then examined the effects of environmental variables on the temporal and spatial patterns of predicted abundance and distribution of the superb parrot using a correlative mapping approach. Results Key findings from GAM analysis were: (1) there was a strong positive relationship between abundance and plant productivity in all regions, but (2) the response of abundance to other predictor variables often differed between regions. Correlative mapping predictions of the abundance and distribution of the superb parrot also indicated that: (1) predicted abundance varied through time and space, (2) predicted abundance sometimes decreased in all regions, but at other times some regions had high abundance when others had low, and (3) changes in plant productivity (and therefore climate) were associated with this variation. Main conclusion The superb parrot favours productive landscapes that are also favoured for agriculture. Movements appear to be associated with seasonal and year to year climate variability. Thus, variation in the recorded abundance of the superb parrot may mask population trends, suggesting that existing population estimates are unreliable. Also, high abundances in some areas, and at some times, may reflect deteriorating habitat conditions elsewhere rather than species recovery. Temporal variability in the distribution of the superb parrot makes it difficult to identify specific drought refugia. Consequently, through time, as key habitat continues to deteriorate, the species will become increasingly vulnerable and threatened. Whole‐landscape habitat conservation and restoration strategies are therefore needed to sustain superb parrot populations in the long‐term.     

Predicting insect phenology across space and time

Many species appear to be undergoing shifts in phenology, arising from climate change. To predict the direction and magnitude of future changes requires an understanding of how phenology depends on climatic variation. Species show large‐scale spatial variation in phenology (affected by differentiation among populations) as well as variation in phenology from year‐to‐year at the same site (affected predominantly by local plasticity). Teasing apart spatial and temporal variation in phenology should allow improved predictions of phenology under climate change. This study is the first to quantify large‐scale spatial and temporal variation in the entire emergence pattern of species, and to test the relationships found by predicting future data. We use data from up to 33 years of permanent transect records of butterflies in the United Kingdom to fit and test models for 15 butterfly species. We use generalized additive models to model spatial and temporal variation in the distribution of adult butterflies over the season, allowing us to capture changes in the timing of emergence peaks, relative sizes of peaks and/or number of peaks in a single analysis. We develop these models using data for 1973–2000, and then use them to predict phenologies from 2001 to 2006. For six of our study species, a model with only spatial variation in phenology is the best predictor of the future, implying that these species have limited plasticity. For the remaining nine species, the best predictions come from a model with both spatial and temporal variation in phenology; for four of these, growing degree‐days have similar effects over space and time, implying high levels of plasticity. The results show that statistical phenology models can be used to predict phenology shifts in a second time period, suggesting that it should be feasible to project phenologies under climate change scenarios, at least over modest time scales.     

The transparency,reliability and utility of tropical rainforest land‐use and land‐cover change models

Land‐use and land‐cover (LULC) change is one of the largest drivers of biodiversity loss and carbon emissions globally. We use the tropical rainforests of the Amazon, the Congo basin and South‐East Asia as a case study to investigate spatial predictive models of LULC change. Current predictions differ in their modelling approaches, are highly variable and often poorly validated. We carried out a quantitative review of 48 modelling methodologies, considering model spatio‐temporal scales, inputs, calibration and validation methods. In addition, we requested model outputs from each of the models reviewed and carried out a quantitative assessment of model performance for tropical LULC predictions in the Brazilian Amazon. We highlight existing shortfalls in the discipline and uncover three key points that need addressing to improve the transparency, reliability and utility of tropical LULC change models: (1) a lack of openness with regard to describing and making available the model inputs and model code; (2) the difficulties of conducting appropriate model validations; and (3) the difficulty that users of tropical LULC models face in obtaining the model predictions to help inform their own analyses and policy decisions. We further draw comparisons between tropical LULC change models in the tropics and the modelling approaches and paradigms in other disciplines, and suggest that recent changes in the climate change and species distribution modelling communities may provide a pathway that tropical LULC change modellers may emulate to further improve the discipline. Climate change models have exerted considerable influence over public perceptions of climate change and now impact policy decisions at all political levels. We suggest that tropical LULC change models have an equally high potential to influence public opinion and impact the development of land‐use policies based on plausible future scenarios, but, to do that reliably may require further improvements in the discipline.