Exploring consensus in 21st century projections of climatically suitable areas for African vertebrates

Africa is predicted to be highly vulnerable to 21st century climatic changes. Assessing the impacts of these changes on Africa's biodiversity is, however, plagued by uncertainties, and markedly different results can be obtained from alternative bioclimatic envelope models or future climate projections. Using an ensemble forecasting framework, we examine projections of future shifts in climatic suitability, and their methodological uncertainties, for over 2500 species of mammals, birds, amphibians and snakes in sub‐Saharan Africa. To summarize a priori the variability in the ensemble of 17 general circulation models, we introduce a consensus methodology that combines co‐varying models. Thus, we quantify and map the relative contribution to uncertainty of seven bioclimatic envelope models, three multi‐model climate projections and three emissions scenarios, and explore the resulting variability in species turnover estimates. We show that bioclimatic envelope models contribute most to variability, particularly in projected novel climatic conditions over Sahelian and southern Saharan Africa. To summarize agreements among projections from the bioclimatic envelope models we compare five consensus methodologies, which generally increase or retain projection accuracy and provide consistent estimates of species turnover. Variability from emissions scenarios increases towards late‐century and affects southern regions of high species turnover centred in arid Namibia. Twofold differences in median species turnover across the study area emerge among alternative climate projections and emissions scenarios. Our ensemble of projections underscores the potential bias when using a single algorithm or climate projection for Africa, and provides a cautious first approximation of the potential exposure of sub‐Saharan African vertebrates to climatic changes. The future use and further development of bioclimatic envelope modelling will hinge on the interpretation of results in the light of methodological as well as biological uncertainties. Here, we provide a framework to address methodological uncertainties and contextualize results.     

Partitioning and mapping uncertainties in ensembles of forecasts of species turnover under climate change

Forecasts of species range shifts under climate change are fraught with uncertainties and ensemble forecasting may provide a framework to deal with such uncertainties. Here, a novel approach to partition the variance among modeled attributes, such as richness or turnover, and map sources of uncertainty in ensembles of forecasts is presented. We model the distributions of 3837 New World birds and project them into 2080. We then quantify and map the relative contribution of different sources of uncertainty from alternative methods for niche modeling, general circulation models (AOGCM), and emission scenarios. The greatest source of uncertainty in forecasts of species range shifts arises from using alternative methods for niche modeling, followed by AOGCM, and their interaction. Our results concur with previous studies that discovered that projections from alternative models can be extremely varied, but we provide a new analytical framework to examine uncertainties in models by quantifying their importance and mapping their patterns.     

Integrating bioclimate with population models to improve forecasts of species extinctions under climate change

Climate change is already affecting species worldwide, yet existing methods of risk assessment have not considered interactions between demography and climate and their simultaneous effect on habitat distribution and population viability. To address this issue, an international workshop was held at the University of Adelaide in Australia, 25–29 May 2009, bringing leading species distribution and population modellers together with plant ecologists. Building on two previous workshops in the UK and Spain, the participants aimed to develop methodological standards and case studies for integrating bioclimatic and metapopulation models, to provide more realistic forecasts of population change, habitat fragmentation and extinction risk under climate change. The discussions and case studies focused on several challenges, including spatial and temporal scale contingencies, choice of predictive climate, land use, soil type and topographic variables, procedures for ensemble forecasting of both global climate and bioclimate models and developing demographic structures that are realistic and species-specific and yet allow generalizations of traits that make species vulnerable to climate change. The goal is to provide general guidelines for assessing the Red-List status of large numbers of species potentially at risk, owing to the interactions of climate change with other threats such as habitat destruction, overexploitation and invasive species.     

Strengthening forecasts of climate change impacts with multi‐model ensemble averaged projections using MAGICC/SCENGEN 5.3

Climate output from general circulation models (GCMs) is being used with increasing frequency to explore potential climate change impacts on species’ distributional range shifts and extinction probability. However, different GCMs do not perform equally well in their ability to hindcast the key climatic factors that potentially influence species distributions. Previous research has demonstrated that multi‐model ensemble forecasts perform better than any single GCM in simulating observed conditions at a global scale. MAGICC/SCENGEN 5.3 is a freeware climate model ‘emulator’ that generates multi‐model ensemble forecasts, conditional on regional and/or global performance, for up to twenty GCMs. In combination with a new application ‘M/SGridder’, this software can be used to produce down‐scaled ensemble forecasts, which minimize climate‐model‐related uncertainty, for a range of ecological problems.     

Conservation management and sustainable harvest quotas are sensitive to choice of climate modelling approach for two marine gastropods

Aim

To establish the robustness of two alternative methods for predicting the future ranges and abundances for two wild‐harvested abalone species (Haliotis rubra Donovan 1808 and H. laevigata Leach 1814): single atmosphere–ocean general circulation model (GCM) or ensemble‐averaged GCM forecasts.

Location

South Australia.

Methods

We assessed the ability of 20 GCMs to simulate observed seasonal sea surface temperature (SST) between 1980–1999, globally, and regionally for the Indian and Pacific Oceans south of the Equator. We used model rankings to characterize a set of representative climate futures, using three different‐sized GCM ensembles and two individual GCMs (the Parallel Climate Model and the Community Climate System Model, version 3.0). Ecological niche models were then coupled to physiological information to compare forecast changes in area of occupancy, population size and harvest area based on forecasts using the various GCM selection methods, as well as different greenhouse gas emission scenarios and climate sensitivities.

Results

We show that: (1) the skill with which climate models reproduce recent SST records varies considerably amongst GCMs, with multimodel ensemble averages showing closer agreement to observations than single models; (2) choice of GCM, and the decision on whether or not to use ensemble‐averaged climate forecasts, can strongly influence spatiotemporal predictions of range, abundance and fishing potential; and (3) comparable hindcasting skill does not necessarily guarantee that GCM forecasts and ecological and evolutionary responses to these forecast changes, will be similar amongst closely ranked models.

Conclusion

By averaging across an ensemble of seven highly ranked skilful GCMs, inherent uncertainties stemming from GCM differences are incorporated into forecasts of change in species range, abundance and sustainable fishing area. Our results highlight the need to make informed and explicit decisions on GCM choice, model sensitivity and emission scenarios when exploring conservation options for marine species and the sustainability of future harvests using ecological niche models.

     

A comparison of absolute performance of different correlative and mechanistic species distribution models in an independent area

To investigate the comparative abilities of six different bioclimatic models in an independent area, utilizing the distribution of eight different species available at a global scale and in Australia. Global scale and Australia. We tested a variety of bioclimatic models for eight different plant species employing five discriminatory correlative species distribution models (SDMs) including Generalized Linear Model (GLM), MaxEnt, Random Forest (RF), Boosted Regression Tree (BRT), Bioclim, together with CLIMEX (CL) as a mechanistic niche model. These models were fitted using a training dataset of available global data, but with the exclusion of Australian locations. The capabilities of these techniques in projecting suitable climate, based on independent records for these species in Australia, were compared. Thus, Australia is not used to calibrate the models and therefore it is as an independent area regarding geographic locations. To assess and compare performance, we utilized the area under the receiver operating characteristic (ROC) curves (AUC), true skill statistic (TSS), and fractional predicted areas for all SDMs. In addition, we assessed satisfactory agreements between the outputs of the six different bioclimatic models, for all eight species in Australia. The modeling method impacted on potential distribution predictions under current climate. However, the utilization of sensitivity and the fractional predicted areas showed that GLM, MaxEnt, Bioclim, and CL had the highest sensitivity for Australian climate conditions. Bioclim calculated the highest fractional predicted area of an independent area, while RF and BRT were poor. For many applications, it is difficult to decide which bioclimatic model to use. This research shows that variable results are obtained using different SDMs in an independent area. This research also shows that the SDMs produce different results for different species; for example, Bioclim may not be good for one species but works better for other species. Also, when projecting a “large” number of species into novel environments or in an independent area, the selection of the “best” model/technique is often less reliable than an ensemble modeling approach. In addition, it is vital to understand the accuracy of SDMs' predictions. Further, while TSS, together with fractional predicted areas, are appropriate tools for the measurement of accuracy between model results, particularly when undertaking projections on an independent area, AUC has been proved not to be. Our study highlights that each one of these models (CL, Bioclim, GLM, MaxEnt, BRT, and RF) provides slightly different results on projections and that it may be safer to use an ensemble of models.     

High-resolution bioclimatic dataset derived from future climate projections for plant species distribution modeling

Species distribution modeling is playing an increasingly prominent role in ecology and global change biology, owing to its potential to predict species range shifts, biodiversity losses, and biological invasion risks for future climates. Such models are now well-established as important tools for biological conservation. However, the lack of high-resolution data for future climate scenarios has seriously limited their application, particularly because of the scale gap between general circulation models (GCMs) and species distribution models (SDMs). A recently introduced change-factor downscaling technique provides a convenient way to build high-resolution datasets from GCM projections. Here, we present a high-resolution (10’ × 10’) global bioclimatic dataset (BioPlant) for plant species distribution. The 15 bioclimatic variables we select are considered those most eco-physiologically relevant. They can be easily calculated from climatic variables common to all GCM projections. In addition to the traditional classes of variables regarding temperature and precipitation, the BioPlant dataset emphasizes the interactions between temperature and precipitation, particularly within plant growing seasons. A preliminary visual analysis shows that variations among GCMs are more significant on a species range scale than on a global scale. Thus, the formerly advocated ensemble modeling method should be applied not only to different SDMs, but also to various GCMs. Statistic analysis suggests that divergent behavior among GCM variations for temperature class variables and classes of precipitation variables requires special attention. Our dataset may provide a common platform for ensemble modeling, and can serve as an example to develop higher-resolution regional datasets.     

An ensemble data assimilation modeling system for operational outdoor microalgae growth forecasting

Microalgae have received increasing attention as a potential feedstock for biofuel or biobased products. Forecasting the microalgae growth is beneficial for managers in planning pond operations and harvesting decisions. This study proposed a biomass forecasting system comprised of the Huesemann Algae Biomass Growth Model (BGM), the Modular Aquatic Simulation System in Two Dimensions (MASS2), ensemble data assimilation (DA), and numerical weather prediction Global Ensemble Forecast System (GEFS) ensemble meteorological forecasts. The novelty of this study is to seek the use of ensemble DA to improve both BGM and MASS2 model initial conditions with the assimilation of biomass and water temperature measurements and consequently improve short-term biomass forecasting skills. This study introduces the theory behind the proposed integrated biomass forecasting system, with an application undertaken in pseudo-real-time in three outdoor ponds cultured with Chlorella sorokiniana in Delhi, California, United States. Results from all three case studies demonstrate that the biomass forecasting system improved the short-term (i.e., 7-day) biomass forecasting skills by about 60% on average, comparing to forecasts without using the ensemble DA method. Given the satisfactory performances achieved in this study, it is probable that the integrated BGM-MASS2-DA forecasting system can be used operationally to inform managers in making pond operation and harvesting planning decisions.     

Conserving the Brazilian semiarid (Caatinga) biome under climate change

The Caatinga is a semiarid biome of the northeast of Brazil with only 1?% of its territory currently conserved. The biome’s biodiversity is highly threatened due to exposure to land conversion for agricultural and cattle ranch. Climate forecasts predict increases in aridity, which could pose additional threats to the biome’s biodiversity. Here, we ask if the remnants of natural vegetation in Caatinga biome, where endemic terrestrial vertebrate species occur, are likely to retain more climatic suitability under climate change scenarios than other less pristine areas of the biome. In order to assess changes in climate suitability across individual species ranges, ensemble forecasting was used based on seven bioclimatic envelope models, three atmosphere–ocean general circulation models, and two greenhouse emission gas scenarios for 2020, 2050, and 2080. We found that most species will gain climatic suitability in the natural vegetation remnants of the Caatinga. Such gains are even greater than the expected to occur within random sets of areas with size similar to the natural vegetation remnants. Our results suggest that natural vegetation remnants will likely play a role of climate refuges for endemic vertebrate species, so efforts should be concentrated in these regions.     

Land-cover data improve bioclimatic models for anurans and turtles at a regional scale

Aim We investigated whether accounting for land cover could improve bioclimatic models for eight species of anurans and three species of turtles at a regional scale. We then tested whether accounting for spatial autocorrelation could significantly improve bioclimatic models after statistically controlling for the effects of land cover. Location Nova Scotia, eastern Canada. Methods Species distribution data were taken from a recent (1999–2003) herpetofaunal atlas. Generalized linear models were used to relate the presence or absence of each species to climate and land‐cover variables at a 10‐km resolution. We then accounted for spatial autocorrelation using an autocovariate or third‐order trend surface of the geographical coordinates of each grid square. Finally, variance partitioning was used to explore the independent and joint contributions of climate, land cover and spatial autocorrelation. Results The inclusion of land cover significantly increased the explanatory power of bioclimatic models for 10 of the 11 species. Furthermore, including land cover significantly increased predictive performance for eight of the 11 species. Accounting for spatial autocorrelation improved model fit for rare species but generally did not improve prediction success. Variance partitioning demonstrated that this lack of improvement was a result of the high correlation between climate and trend‐surface variables. Main conclusions The results of this study suggest that accounting for the effects of land cover can significantly improve the explanatory and predictive power of bioclimatic models for anurans and turtles at a regional scale. We argue that the integration of climate and land‐cover data is likely to produce more accurate spatial predictions of contemporary herpetofaunal diversity. However, the use of land‐cover simulations in climate‐induced range‐shift projections introduces additional uncertainty into the predictions of bioclimatic models. Further research is therefore needed to determine whether accounting for the effects of land cover in range‐shift projections is merited.     

NeuralEnsembles: a neural network based ensemble forecasting program for habitat and bioclimatic suitability analysis

NeuralEnsembles is an integrated modeling and assessment tool for predicting areas of species habitat/bioclimatic suitability based on presence/absence data. This free, Windows based program, which comes with a friendly graphical user interface, generates predictions using ensembles of artificial neural networks. Models can quickly and easily be produced for multiple species and subsequently be extrapolated either to new regions or under different future climate scenarios. An array of options is provided for optimizing the construction and training of ensemble models. Main outputs of the program include text files of suitability predictions, maps and various statistical measures of model performance and accuracy.     

Does seasonal fine‐tuning of climatic variables improve the performance of bioclimatic envelope models for migratory birds?

We examined the influence of 'seasonal fine-tuning' of climatic variables on the performance of bioclimatic envelope models of migrating birds. Using climate data and national bird atlas data from a 10 × 10 km uniform grid system in Finland, we tested whether the replacement of one 'baseline' set of variables including summer (June–August) temperature and precipitation variables with climate variables tailored ('fine-tuned') for each species individually improved the bird-climate models. The fine-tuning was conducted on the basis of time of arrival and early breeding of the species. Two generalized additive models (GAMs) were constructed for each of the 63 bird species studied, employing (1) the baseline climate variables and (2) the fine-tuned climate variables. Model performance was measured as explanatory power (deviance change) and predictive power (area under the curve; AUC) statistics derived from cross-validation. Fine-tuned climate variables provided, in many cases, statistically significantly improved model performance compared to using the same baseline set of variables for all the species. Model improvements mainly concerned bird species arriving and starting their breeding in May–June. We conclude that the use of the fine-tuned climate variables tailored for each species individually on the basis of their arrival and critical breeding periods can provide important benefits for bioclimatic modelling.     

Do we need land‐cover data to model species distributions in Europe?

Aim To assess the influence of land cover and climate on species distributions across Europe. To quantify the importance of land cover to describe and predict species distributions after using climate as the main driver. Location The study area is Europe. Methods (1) A multivariate analysis was applied to describe land‐cover distribution across Europe and assess if the land cover is determined by climate at large spatial scales. (2) To evaluate the importance of land cover to predict species distributions, we implemented a spatially explicit iterative procedure to predict species distributions of plants (2603 species), mammals (186 species), breeding birds (440 species), amphibian and reptiles (143 species). First, we ran bioclimatic models using stepwise generalized additive models using bioclimatic variables. Secondly, we carried out a regression of land cover (LC) variables against residuals from the bioclimatic models to select the most relevant LC variables. Finally, we produced mixed models including climatic variables and those LC variables selected as decreasing the residual of bioclimatic models. Then we compared the explanatory and predictive power of the pure bioclimatic against the mixed model. Results (1) At the European coarse resolution, land cover is mainly driven by climate. Two bioclimatic axes representing a gradient of temperature and a gradient of precipitation explained most variation of land‐cover distribution. (2) The inclusion of land cover improved significantly the explanatory power of bioclimatic models and the most relevant variables across groups were those not explained or poorly explained by climate. However, the predictive power of bioclimatic model was not improved by the inclusion of LC variables in the iterative model selection process. Main conclusion Climate is the major driver of both species and land‐cover distributions over Europe. Yet, LC variables that are not explained or weakly associated with climate (inland water, sea or arable land) are interesting to describe particular mammal, bird and tree distributions. However, the addition of LC variables to pure bioclimatic models does not improve their predictive accuracy.     

Avoidable errors in the modelling of outbreaks of emerging pathogens, with special reference to Ebola

As an emergent infectious disease outbreak unfolds, public health response is reliant on information on key epidemiological quantities, such as transmission potential and serial interval. Increasingly, transmission models fit to incidence data are used to estimate these parameters and guide policy. Some widely used modelling practices lead to potentially large errors in parameter estimates and, consequently, errors in model-based forecasts. Even more worryingly, in such situations, confidence in parameter estimates and forecasts can itself be far overestimated, leading to the potential for large errors that mask their own presence. Fortunately, straightforward and computationally inexpensive alternatives exist that avoid these problems. Here, we first use a simulation study to demonstrate potential pitfalls of the standard practice of fitting deterministic models to cumulative incidence data. Next, we demonstrate an alternative based on stochastic models fit to raw data from an early phase of 2014 West Africa Ebola virus disease outbreak. We show not only that bias is thereby reduced, but that uncertainty in estimates and forecasts is better quantified and that, critically, lack of model fit is more readily diagnosed. We conclude with a short list of principles to guide the modelling response to future infectious disease outbreaks.     

Bioclimatic Modeling of Southern African Bioregions and Biomes Using Bayesian Networks

This study uses Bayesian networks (BNs) to simulate the spatial distribution of southern African biomes and bioregions using bioclimatic variables. Two Tree-Augmented Naïve (TAN) BN models were parameterized from 23 bioclimatic variables using the expectation-maximization (EM) algorithm. Using sensitivity analyses, the relative influence of each variable was determined using the mutual information from which six bioclimatic variables were selected for the final models. Precipitation of the warmest quarter and extra-terrestrial solar radiation was found to be the most influential variables on both bioregion and biome distributions. Isothermality was the least influential bioclimatic variable at both bioregion and biome levels. Overall correspondence was very high at 93.8 and 87.1% for biomes and bioregions, respectively, whereas classification errors were obtained in transition areas indicating the uncertainties associated with vegetation mapping around margins. The findings indicate that southern African bioregions and biomes can be classified and mapped according to key bioclimatic variables. Spatio-temporal, in particular, monthly and quarterly variations in both precipitation and temperature are found to be ecologically significant in determining the spatial distribution of biomes and bioregions. The findings also reflect the hierarchical relationship of biomes and bioregions as a function of local bioclimatic gradients and interactions. The results indicate the ecological significance of bioclimatic conditions in ecosystem science and offer the opportunity to utilize the models for predicting future responses and sensitivities to climatic changes.     

Effects of climate change on an emperor penguin population: analysis of coupled demographic and climate models

Sea ice conditions in the Antarctic affect the life cycle of the emperor penguin (Aptenodytes forsteri). We present a population projection for the emperor penguin population of Terre Adélie, Antarctica, by linking demographic models (stage‐structured, seasonal, nonlinear, two‐sex matrix population models) to sea ice forecasts from an ensemble of IPCC climate models. Based on maximum likelihood capture‐mark‐recapture analysis, we find that seasonal sea ice concentration anomalies (SIC_a) affect adult survival and breeding success. Demographic models show that both deterministic and stochastic population growth rates are maximized at intermediate values of annual SIC_a, because neither the complete absence of sea ice, nor heavy and persistent sea ice, would provide satisfactory conditions for the emperor penguin. We show that under some conditions the stochastic growth rate is positively affected by the variance in SIC_a. We identify an ensemble of five general circulation climate models whose output closely matches the historical record of sea ice concentration in Terre Adélie. The output of this ensemble is used to produce stochastic forecasts of SIC_a, which in turn drive the population model. Uncertainty is included by incorporating multiple climate models and by a parametric bootstrap procedure that includes parameter uncertainty due to both model selection and estimation error. The median of these simulations predicts a decline of the Terre Adélie emperor penguin population of 81% by the year 2100. We find a 43% chance of an even greater decline, of 90% or more. The uncertainty in population projections reflects large differences among climate models in their forecasts of future sea ice conditions. One such model predicts population increases over much of the century, but overall, the ensemble of models predicts that population declines are far more likely than population increases. We conclude that climate change is a significant risk for the emperor penguin. Our analytical approach, in which demographic models are linked to IPCC climate models, is powerful and generally applicable to other species and systems.     

Predicting the future of species diversity: macroecological theory, climate change, and direct tests of alternative forecasting methods

Accurate predictions of future shifts in species diversity in response to global change are critical if useful conservation strategies are to be developed. The most widely used prediction method is to model individual species niches from point observations and project these models forward using future climate scenarios. The resulting changes in individual ranges are then summed to predict diversity changes; multiple models can be combined to produce ensemble forecasts. Predictions based on environment-richness regressions are rarer. However, richness regression models, based on macroecological diversity theory, have a long track record of making reliable spatial predictions of diversity patterns. If these empirical theories capture true functional relationships between environment and diversity, then they should make consistent predictions through time as well as space and could complement individual species-based predictions. Here, we use climate change throughout the 20th century to directly test the ability of these different approaches to predict shifts of Canadian butterfly diversity. We found that all approaches performed reasonably well, but the most accurate predictions were made using the single best richness-environment regression model, after accounting for the effects of spatial autocorrelation. Spatially trained regression models based on macroecological theory accurately predict diversity shifts for large species assemblages. Global changes provide pseudo-experimental tests of those macroecological theories that can then generate robust predictions of future conditions.     

Review of the physiology of human thermal comfort while exercising in urban landscapes and implications for bioclimatic design

This review comprehensively examines scientific literature pertaining to human physiology during exercise, including mechanisms of heat formation and dissipation, heat stress on the body, the importance of skin temperature monitoring, the effects of clothing, and microclimatic measurements. This provides a critical foundation for microclimatologists and biometeorologists in the understanding of experiments involving human physiology. The importance of the psychological aspects of how an individual perceives an outdoor environment are also reviewed, emphasizing many factors that can indirectly affect thermal comfort (TC). Past and current efforts to develop accurate human comfort models are described, as well as how these models can be used to develop resilient and comfortable outdoor spaces for physical activity. Lack of suitable spaces plays a large role in the deterioration of human health due to physical inactivity, leading to higher rates of illness, heart disease, obesity and heat-related casualties. This trend will continue if urban designers do not make use of current knowledge of bioclimatic urban design, which must be synthesized with physiology, psychology and microclimatology. Increased research is required for furthering our knowledge on the outdoor human energy balance concept and bioclimatic design for health and well-being in urban areas.     

Nowcasting and Forecasting the Monthly Food Stamps Data in the US Using Online Search Data

We propose the use of Google online search data for nowcasting and forecasting the number of food stamps recipients. We perform a large out-of-sample forecasting exercise with almost 3000 competing models with forecast horizons up to 2 years ahead, and we show that models including Google search data statistically outperform the competing models at all considered horizons. These results hold also with several robustness checks, considering alternative keywords, a falsification test, different out-of-samples, directional accuracy and forecasts at the state-level.     

The future of invasive African grasses in South America under climate change

Climate change will promote substantial effects on the distribution of invasive species. Here, I used an ensemble of bioclimatic envelope models (Gower Distance, Chebyshev Distance, and Mahalanobis Distance) to forecast climatically suitable areas of South America for 13 invasive African grass species under future climate conditions (year 2050). Under current climatic conditions, the areas with the potential for the highest invasive species richness are located mostly in the tropical climates of South America, except for the Amazon region. In the year 2050, the overall pattern of invasive species richness will not change considerably, and increases in northeastern Amazon and portions of the temperate regions of South America are predicted.