UNCERTAINTY IN PROJECTED IMPACTS OF CLIMATE CHANGE ON BIODIVERSITY
A FOCUS ON AFRICAN VERTEBRATES
RAQUEL A. GARCIA University of Copenhagen | 2014
Uncertainty in projected impacts of climate change on biodiversity A focus on African vertebrates
Raquel A. Garcia PhD thesis | January 2014 University of Copenhagen | Faculty of Science
ACADEMIC ADVISORS Prof Miguel B. Araújo CMEC, University of Copenhagen, Denmark Imperial College, Silwood Park, UK National Museum of Natural History–CSIC, Spain InBio/CIBIO, Évora University, Portugal
Dr Mar Cabeza Metapopulation Research Group, University of Helsinki, Finland
SUBMITTED TO The PhD School of The Faculty of Science, University of Copenhagen, Denmark January 2014
ASSESSMENT COMMITTEE Prof Jane Hill University of York, Department of Biology, UK
Dr Richard Pearson University College London, Research Department of Genetics, Evolution and Environment, UK
Dr David Nogués-Bravo CMEC, University of Copenhagen, Denmark
COPYRIGHT © 2014 Raquel A. Garcia (Synopsis, Design) © 2012 Blackwell Publishing Ltd (Chapter I) © 2014 The Authors Journal of Biogeography Published by John Wiley & Sons Ltd (Chapter II) © 2014 The Authors (Chapter III) © 2014 The Authors (Chapter IV) © 2014 The Authors (Chapter V)
Contents | iii
Contents
Preface v Summary / Resumé vii Acknowledgements x
Synopsis 1 The uncertain nature of assessments under climate change 5 Variation in data and model decisions 9 Robustness to ecological assumptions 11 Broadening the scope of assessments 13 Embracing uncertainty 15 References 18
Chapter I 27 Exploring consensus in 21st century projections of climatically suitable areas for African vertebrates
Chapter II 65 Matching species traits to projected threats and opportunities from climate change
Chapter III 87 Conservation implications of omitting rare and threatened species from climate change impact modelling
Chapter IV 115 Multiple dimensions of climate change and implications for biodiversity
Chapter V 147 Do projections from bioclimatic envelope models and climate change metrics match?
Preface | v
Preface
This thesis is the result of a four-year PhD pro- ject based at the Center for Macroecology, Evo- lution and Climate at the University of Copen- hagen, in Denmark, and the National Museum of Natural History (CSIC), in Madrid, Spain. The project was supervised by Prof Dr Miguel Araújo and Dr Mar Cabeza. The thesis work also included various stays at the University of Hel- sinki in Finland, the 'Rui Nabeiro' Biodiversity Chair at Évora University in Portugal, and the South African National Biodiversity Institute in Cape Town, South Africa. The work was funded by a PhD grant from the Portuguese Foundation for Science and Technology. The thesis consists of two parts. The first part is a synopsis describing the background and aims of the thesis, summarising the main findings, and discussing the work presented in a broader perspective. The second part consists of five chapters that form the core of the thesis. They include four analytical chapters and one review chapter. At the time of submission, two were published as scientific articles, and the remaining three were in review.
Summary / Resumé | vii
Summary
Evidence for shifts in the phenologies and dis- need to both integrate uncertainties in assess- tributions of species over recent decades has ments, and reduce or circumvent them where often been attributed to climate change. The possible. Integration of uncertainties is illus- prospect of greater and faster changes in cli- trated in two examples. The first uses ensem- mate during the 21st century has spurred a bles of bioclimatic envelope models for sub- stream of studies anticipating future biodiver- Saharan African vertebrates, built with alterna- sity impacts. Yet, uncertainty is inherent to tive climate data and model algorithms. En- both projected climate changes and their effects semble forecasting provides a means for ex- on biodiversity, and needs to be understood ploring the breadth and spatial variation of before projections can be used. This thesis uncertainties, and for building consensus seeks to elucidate some of the uncertainties among projections. Several consensus method- clouding assessments of biodiversity impacts ologies are compared here, including a newly from climate change, and explores ways to proposed methodology that preserves informa- address them. While the focus is mostly on sub- tion about the variability of projections in the Saharan African vertebrates, the methodologi- ensemble. The second example examines model cal advances and conclusions presented are far- outputs for sub-Saharan African amphibians in reaching and have wider relevance. the light of species' vulnerability to climate Throughout the chapters in this thesis, pro- change. An analytical framework is developed jections under changing climates for sub- for distinguishing between different climatic Saharan African vertebrates, based on biocli- threats and opportunities revealed by the bio- matic envelope models, are shown to be af- climatic envelope models, and analysing how fected by multiple uncertainties. Different they each are altered by the consideration of model algorithms produce different outputs, as specific response-mediating traits. do alternative future climate models and sce- Efforts to reduce uncertainties in biodiver- narios of future emissions of greenhouse gases. sity impact assessments are equally important. Another uncertainty arises due to omission of However, many sources of uncertainty cannot species with small sample sizes, which are easily be reduced, not least the omission of difficult to model. The effect of such bias species that are narrow-ranging, poorly known, against narrow-ranging species is often over- or even unknown to science. This uncertainty looked in assessments of biodiversity impacts, can instead be circumvented through the use of but our results for sub-Saharan African am- alternative approaches to assessing impacts. phibians show that it trickles down to conser- This thesis discusses one candidate approach vation strategies. Finally, assumptions about that is independent of species' data: the use of the climatic tolerance of species, their dispersal climate change metrics. By describing the expo- ability, and other characteristics are also shown sure of regions to multiple changes in the mag- to alter model projections for sub-Saharan nitude, timing, position, or availability of cli- African amphibians. matic conditions, metrics can provide infer- Despite numerous calls to address the un- ences about the potential threats and opportu- certainty challenge, appropriate treatment of nities for the biodiversity in those regions. The uncertainty has yet to be formalised in assess- diversity of existing metrics is reviewed here, ments of biodiversity impacts under climate and the picture that emerges is one of multifac- change. The chapters in this thesis highlight the eted changes in climate, with unequal spatial viii | Summary / Resumé
patterns around the world. To help interpret The uncertainties discussed in this thesis, the diversity of climate change metrics, a con- and many others not covered here, impair the ceptual framework is proposed for using them conservation of biodiversity under changing in biodiversity impact assessments. Early test- climates in Africa and elsewhere. Explicitly ing of this framework, by comparing inferences addressing all uncertainties of projected im- from metrics and from the bioclimatic envelope pacts appears overwhelming. Yet, if model models for sub-Saharan African amphibians, projections are to be useful for conservation suggests that climate change metrics might be a planners, they must be as transparent as possi- useful addition to the biodiversity impact as- ble by including an honest description of their sessment toolbox. level of confidence given the current knowl- edge.
Resumé
De seneste årtiers evidens for ændringer i konklusioner mere vidt-rækkende og har bred arters fænologi og geografiske udbredelser er relevans. ofte blevet tilskrevet klimaforandringer. En gennemgående tendens i denne Udsigten til hurtigere og mere omfattende afhandling er at forudsigelser for afrikanske ændringer i klimaet i det 21. århundrede har hvirveldyr, baseret på bioclimatic envelope affødt en strøm af studier, der søger at models (klima-baserede modeller for en arts klarlægge de forventede konsekvenser for udbredelse) under klimaforandringer, er fremtidens biodiversitet. Dog ligger der en påvirkede af flere usikkerheder. Forskellige iboende usikkerhed i forudsigelser af model-algoritmer giver forskellige resultater, klimaforandringer, og deres indflydelse på og det gør alternative klimamodeller og biodiversiteten, som det er nødvendig at afklare scenarier for fremtidige drivhusgas- før disse forudsigelser kan anvendes. Denne udledninger også. En anden usikkerhed opstår afhandling sigter mod at udrede nogle af de som følge af at arter med små sample-sizes usikkerheder der forplumrer en vurdering af udelukkes, fordi de er svære at modellere. klimaforandringernes indflydelse på Effekten af denne skævvridning væk fra arter biodiversiteten, og udforsker måder at med små udbredelser overses ofte i håndtere dem på. Selv om fokus hovedsagligt vurderinger af konsekvenser for ligger på afrikanske hvirveldyr syd for Sahara, biodiversiteten, men vores resultater for er de metodologiske fremskridt og padder syd for Sahara viser at der er en effekt,
Summary / Resumé | ix
der afspejles helt ned i konkrete strategier for udbredelser, og som man ved meget lidt eller naturbevarelse. Endelig vises det at antagelser intet om. Denne usikkerhed kan istedet omgåes om arters klimatiske tolerance, spredningsevne ved at anvende alternative tilgange. Denne og andre egenskaber ændrer model- afhandling diskuterer én mulig tilgang, som er forudsigelser for afrikanske padder. uafhængig af artsdata: at anvende specifikke Til trods for adskillige opfordringer til at enheder til at måle klimaforandringer. Ved at håndtere den udfordring som usikkerhed beskrive et områdes eksponering for ændringer udgør, så mangler der stadig en passende i klimatiske forholds udbredelse, timing, formaliseret tilgang til at behandle usikkerhed i placering, eller tilgængelighed, kan sådanne vurderingerne af klimaforandringernes måleenheder bruges til at udlede noget om de konsekvens for biodiversiteten. Kapitlerne i potentielle trusler mod, og muligheder for, denne afhandling understreger biodiversiteten i de pågældende regioner. En nødvendigheden af både at integrere gennemgang af den diverse gruppe af usikkerheden i vurderinger, og at reducere eller eksisterende klimaforandringsmål viser et omgå den når det er muligt. Integration af billede af multifaceterede klimaforandringer usikkerhed illustreres med to eksempler. Det med forskellige rumlige mønstre rundt om i første gør brug af 'ensembles' verden. Som hjælp til at fortolke disse (sammenstillinger) af flere bioclimatic envelope forskellige klimaforandringsmål foreslåes her models for afrikanske vertebrater syd for en konceptuel ramme for hvordan de kan Sahara, konstrueret med alternative klimadata bruges til at vurdere påvirkninger på og model-algoritmer. Ensemble-forudsigelse gør biodiversiteten. Indledende tests af denne det muligt at udforske bredden og den rumlige ramme, foretaget ved at sammenligne variation i graden af usikkerhed, og at etablere fortolkninger af klimaforandringsmål med en konsensus mellem forudsigelser. Her bioclimatic envelope models for afrikanske sammenlignes flere tilgange til at etablere padder syd for Sahara, peger på at konsensus mellem forudsigelser, inklusive en klimaforandringsmål kan blive endnu et nyligt foreslået metode der bevarer information brugbart værktøj til at vurdere påvirkningen om variabiliteten mellem de enkelte på biodiversiteten. forudsigelser i ensemblet. Det andet eksempel De usikkerheder som denne afhandling undersøger model-resultater for afrikanske omhandler, og adskillige andre der ikke er padder syd for Sahara, set i lyset af hver arts omtalt her, vanskeliggør beskyttelsen af sårbarhed overfor klimaforandringer. Der biodiversiteten mod klimaforandringer i Afrika udvikles en analytisk ramme for adskillelse af og andre steder. Det kan forekomme de klimatiske trusler og muligheder der påvises overvældende at skulle imødegå alle disse af bioclimatic ensemble models, og for analyse af usikkerheder eksplicit, men hvis forudsigelser hvordan disse ændrer sig når respons- baseret på modeller skal være nyttige for medierende træk inkorporeres. naturbevarelsen skal de være så gennemsigtige Det er også vigtigt at søge at reducere som muligt. Dette kan opnås ved at inkludere usikkerheden i vurderingen af påvirkninger på en beskrivelse af forudsigelsernes grad af biodiversiteten. Der er dog mange kilder til usikkerhed, som er ærlig og baseret på den usikkerhed der ikke let lader sig reducere, ikke nyeste viden. mindst udeladelsen af arter der har små
x | Acknowledgements
Acknowledgements
This was not a lone journey. Throughout the During the four years I was privileged to four years of my PhD there were various people work in five labs and to meet many incredible pointing me to new paths, encouraging me to people and researchers, young and senior alike. look forward, and reminding me to slow down I kept my status of 'nomad', constantly moving when I needed. between labs and trying to catch the summers My deepest gratitude goes to my supervi- of both hemispheres. Yet, all five labs felt like sors Miguel Araújo and Mar Cabeza. Thank you 'home' to me, for there was always a warm both for trusting from the start (and some- welcome from my colleagues in each one of times, along the way, more than myself) that I them. Gracias to my colleagues at the Museum could do this. Miguel, thank you for instilling in in Madrid for their support during my first year me the capacity to work independently. Part of of science (and Spanish language). Tak to the my research developed because your answers people at CMEC in Copenhagen for their conta- to my questions often broadened, instead of gious enthusiasm for science. Obrigada to all at narrowing, the options. If my 'consultancy style' the Biodiversity Chair in Évora for the warm writing has improved to something closer to environment I always found there (outdoors as scientific writing, I owe it to you. Getting back well). Kiitos to the Metapopulation group in the manuscripts all marked in red was often Helsinki, ever dynamic and helpful. And thanks dispiriting, yet always so helpful. Thank you for to all at SANBI in Cape Town, for their inspiring finding the time to discuss science, advise me connection to nature. My gratitude also goes to on career matters, or simply have a cup of tea. the leaders of these groups for keeping the The barbecues in Alpedrete or Évora, and the doors open (Miguel Araújo, Carsten Rahbek, picnics in the parks of Copenhagen will be Ilkka Hanski, Mar Cabeza, and Guy Midgley), greatly missed. Thank you Miguel, Sonia, Ágata and to all the staff for their support. and Quino for your friendship and company Special thanks to my colleagues who have away from science. made this journey easier in so many different Mar, thank you for teaching me so much ways. It would have been a lot harder to learn about science, the ethics of science, and the fun about statistics or GIS, start writing in R, delve of making science. Above all, thank you for into all the new concepts in ecology, and learn being an inspiration to stay in science. I so en- to read and criticise scientific papers without joyed our discussions about species ranges, Diogo Alagador, Hedvig Nenzén, Márcia Bar- traits and climate in your garden, in the green- bosa, Maria Triviño, Carla Pinto Cruz, Regan house, or while walking Ada to sleep. You have Early, Heini Kujala, Laura Meller, Joaquín gently guided me and helped me to grow as a Hortal, Guida Santos, Pep Serra Díaz, François scientist and as a person – thank you for caring Guilhaumon, Dora Neto, Katie Marske, David about my scientific development beyond pub- Nogués-Bravo, Paco Ferri, Raul García-Valdéz, lishing papers, and for pushing me to continue Sara Varela, Marisa Peláez, Isaac Pozo, Jan going after my dreams. Only someone who Dempewolf, and each and every one in the five knows me well could have found the 'yellow labs! A warm thank you to Phoebe Barnard and cottage' for one of my summer stays in Helsinki John Measey for taking me out to the field, to all – the best house I've lived in during these four in the 'Reserve Selection Journal Club' in Hel- years! Thank you Mar, Tero and Ada for being sinki and in CMEC/iBiochange workshops and the best neighbours, cooking palls, and friends. retreats for many fruitful discussions, and to
Acknowledgements | xi
Res Altwegg for the 'Stats Lunches'. Thanks also Andrea, Mar and Tero, Mats and Sofie, Guida to Heini Kujala, Mar Cabeza, Joaquín Hortal, and Joaquín, Agi and Andy, Maya and Rudi, Jon, Regan Early, Hedvig Nenzén, and Katharine Rean and Tania, Gayane, Manuel, and Virginia Marske for helpful comments on the synopsis, for making my work stays in Portugal, Den- and to Anna-Sofie Stensgaard and Michael mark, Finland, Spain, Switzerland, South Africa, Borregaard for translating the summary into US, New Zealand, and Mozambique so much Danish. more enjoyable. For sharing ideas and making my journey Many have been on this journey with me, so much more rewarding, I thank all those I near or far, but always present. Thank you to all have worked with in manuscripts or projects my precious friends and their beautiful fami- during my PhD. Thanks to all the co-authors in lies: Sol, Joana and Murat, Sofia, Ana, Agi, the papers in this thesis: Mar and Miguel, Wen- Naoko, Zé, Cuncas, Sara and João, Carla, Cristi- dy Foden, Neil Burgess, Carsten Rahbek, Chris- na, Pedro, Diogo M., Eilat and Shaul, Jon, Maya tian Hof, Louis Hansen, Phil Platts, Alexander and Rudi, Linda, Emma, Sergio and Elizabeth, Gutsche, and Res Altwegg. It has also been (and João Carlos, Cláudio, Maral and Mark, Tania and will continue being) fun collaborating with Rean, Andrea, Thierry, Kian, Mats, Philip, Burcu, Jessica Forrest, Guy Midgley, Phoebe Barnard, Ania, Firmas and Ana, Anália, Checa, Patrícia V., and all in the PARCC, HarmBio, RLI, and AECID/ Serpa, Manuel, Patrícia and Pedro, Ana Luísa, Mozambique teams. Thanks to Neil Burgess, Guida and Joaquín H., Márcia, Diogo A., Dora, Miguel Araújo and Mar Cabeza for inviting me Joaquín C., Hedvig, Maria, Helen S., Helen C-B., into new collaborations. and many others. Thank you all for being a Besides Miguel and Mar, I had a few other source of encouragement and fun! official trainers. I am grateful to Richard Pear- Stepping out of the journey every now and son, Steven Phillips, Wilfried Thuiller, Nick then into the yoga classes with Esmé de Wet, Zimmermann, Antoine Guisan, Jens-Christian Elena Mironov, and other wonderful instructors Svenning, David Nogués Bravo, Daniel Kissling, was a life saver. Namaste. Being a child again Signe Normand, and Rob Anderson for teaching while playing with all my nephews and nieces me about the illusions and disillusions of mod- was equally relaxing and even more fun! For els; to Luís Carrascal for making statistics seem their smiles and for reminding me of what is less daunting; to José Alexandre Diniz-Filho for important in life, I thank my goddaughters Mia decrypting (a little) phylogenetics; to Joaquín and Ona Adalís. Hortal for teaching me to question the data; and And, always, thanks to my family. Obrigada to Jan and Camilla for opening up the fascinat- aos meus pais e irmã pelo apoio em todas as ing world of pedagogy to me. Then my turn minhas aventuras com ou sem pára-quedas, came. Teaching the practical lectures about dentro ou fora de Portugal. Aos meus pais, models in an international PhD course in Co- Nanda e David, obrigada por me ensinarem penhagen was one of the most gratifying expe- tanto mesmo sem darem conta. Obrigada à Inês riences during the PhD. Thank you, David, for por estar sempre presente, como amiga e irmã. trusting me, and all the students in the 2012 A toda a família, e aos regulares nos almoços de and 2013 classes for making my job so easy and domingo, obrigada por sorrirem e por me for teaching me so much. fazerem rir. As any nomad, I am grateful to many peo- ple for hosting me throughout these four years. This work was funded by the Portuguese Thanks to Inês and Zé, Carla and Luís, Patrícia Foundation for Science and Technology and Pedro, Sol and Jorge, Joana and Murat, (grant SFRH/BD/65615/2009).
'There is something fascinating about science. One gets such wholesale returns of conjecture out of such a trifling investment of fact.'
Mark Twain, Life on the Mississippi, 1883 Synopsis
Uncertainty in projected impacts of climate change on biodiversity - A focus on African vertebrates
RAQUEL A. GARCIA
Synopsis Uncertainty in projected impacts of climate change on biodiversity – A focus on African vertebrates
RAQUEL A. GARCIA1,2,3
1 Center for Macroecology, Evolution and Climate, Department of Biology, University of Copenhagen, Denmark 2 Department of Biogeography and Global Change, National Museum of Natural Sciences, CSIC, Madrid, Spain 3 InBio/CIBIO, University of Évora, Évora, Portugal
Africa's savannahs support a higher diversity of gesting a link between the changes observed in ungulate species than anywhere else in the nature and climate3–6, and attributing it to an- world. It is hard to imagine an Africa without thropogenic climate change7. Faced with the the springbok, the kudu, the giraffe, or the sable prospect of increasing rates of climate change antelope. The Kruger National Park in north- during the 21st century8–11, attention in the eastern South Africa is home to large popula- research community then turned to projecting tions of ungulates, but the numbers of some future potential impacts on biodiversity2. For species began dwindling in about 1988. Re- some of the ungulate populations in the Kruger, searchers Joseph Ogutu and Norman Owen- Ogutu and Owen-Smith projected severe de- Smith examined the data from wildlife census clines by 2016, and even extirpation, should the conducted between 1977 and 1996, and offered same dry season trends persist1,12. As in any a possible explanation for the declines of eleven other field, however, projections are inherently ungulate species1. Extreme reduction in dry uncertain. For projections of the effects of cli- season rainfall, together with warmer tempera- mate change on biodiversity to be improved tures, had led to a decrease in the amount of and judiciously used, it is crucial to appreciate green forage inside the park, while boundary their uncertainty. This thesis is a contribution fencing prevented the animals from finding to understanding some of the uncertainties vegetation elsewhere. clouding familiar assessments of the future of Population declines such as those recorded biodiversity under changing climates, and to for Kruger's ungulates, as well as phenological exploring ways to address them. The focus is shifts and species' distributional changes ob- mostly on sub-Saharan African vertebrates, but served in multiple regions, prompted the first the methodological advances and over-arching studies in climate change ecology2. Together, conclusions are more broadly applicable. these studies created a body of evidence sug-
4 | Synopsis
Box 1 | Glossary Vague or ambiguous language introduces uncertainty in assessments of climate change impacts on biodiversity (Table 1). The terminology of climate change ecology and modelling is constantly evolving, with some terms not yet established and others still contentious. This glossary lists some of these terms, and describes the meaning used for each throughout this thesis, with the aim of enhancing clarity for the reader.
Bioclimatic envelope model: correlative model that uses associations between observations of spe- cies' occurrences and climate to define the sets of climatic conditions under which species are likely to maintain viable populations – the Grinnellian niche of species16,17. There is contention on which term best describes what is modelled given the data used and the underlying assumptions. Alternative terms include 'species distribution model', 'ecological niche model', and 'habitat suitability model'18–22. In Chapter III, the term 'species distribution model' is used as a synonym (but see a proposal to distinguish between the two23). Climate change metric: measure of changes in climatic parameters (e.g., temperature, precipitation) in a given site (cell) or across broad regions (sets of cells). Different dimensions of climate change can be measured, including changes in the magnitude, timing, availability, and position of climatic conditions. Ensemble forecasting: modelling framework in which multiple simulations across more than one set of initial conditions, model classes, parameters, and boundary conditions are generated, each combina- tion representing a possible state of the system being forecast24. Ensembles can be combined to gener- ate consensus projections, with the aim of summarising agreement among the individual projections. Exposure: 'the extent of climate change likely to be experienced by a species or locale'25. Future climate projections applied to bioclimatic envelope models yield estimates of exposure of species to climate change, whereas climate change metrics estimate the exposure of regions to climate change and allow for inferences to be made about the level of exposure of the biodiversity in such regions. Narrow-ranging species: species with small geographical ranges. Here, used to refer to species with small numbers of occurrence records such that statistical modelling becomes unviable. Narrow-ranging species in this sense are likely to include the majority of the species typically regarded as 'rare'26 – this term replaces 'narrow-ranging species' in Chapter III, with rarity in this case defined by a relative cut- off point based on eligibility for modelling. Risk: the magnitude and probability associated with the effect of an adverse event on species or popu- lations, for example, by causing their decline, reduced fitness, genetic loss, or extinction. It requires consideration of both threat and vulnerability27. Threat: extrinsic human and natural adverse events occurring in a given area within a given time27. Opportunities brought by extrinsic events can also be considered28. Climatic threats are the focus of the thesis, but non-climatic threats such as those arising from land-use change would fall under the same definition. Trait: 'any morphological, physiological or phenological feature measurable at the individual level, from the cell to the whole-organism level, without reference to the environment or any other level of organization'29. Often replaced with characteristics of species and their ranges which are expected to summarise relevant (yet unavailable) traits. Uncertainty: The lack of certainty or precision. A state of having limited knowledge where it is impos- sible to exactly describe the existing state, a future outcome, or more than one possible outcome. Vulnerability: the intrinsic susceptibility of species to threats27. Here it includes both i) the degree to which species or populations are affected by changes in climate depending on factors such as ecophysi- ology, life history, and microhabitat preferences, and ii) their capacity to adapt to changes through phenotypic plasticity, genetic diversity, evolutionary rates, life history traits, and dispersal and coloni- zation ability25. Whereas exposure refers to the extrinsic threats facing species, vulnerability captures the intrinsic capacity of species to respond.
Uncertainty in projected impacts of climate change | 5
The uncertain nature of assessments certainty), vague or ambiguous expression of under climate change that knowledge (linguistic uncertainty), and subjective human values and preferences asso- Assessing the potential impacts of climate ciated with the system (decision-making uncer- change on biodiversity relies on different ap- tainty)37,42–44. These uncertainties pervade the proaches spanning a continuum from correla- entire modelling process. They affect the three tive to mechanistic models13. At the correlative components of modelling frameworks defined end of the continuum, bioclimatic envelope by Austin45: decisions regarding the generation models14,15 (see Glossary in Box 1) use the asso- of data used in the modelling (what Austin calls ciation between known species' occurrences the 'data model'), decisions regarding the actual and climate to characterise the sets of suitable modelling process and evaluation (the 'statisti- climatic conditions for species (or the realised cal model' in Austin's terminology as applied to subsets thereof). By projecting this characteri- correlative approaches, or, more generally sation to the future, these models describe here, the 'mathematical model'), and the eco- potential changes in the geographical patterns logical knowledge and theory used (the 'eco- of climatic suitability for species, and provide logical model'). Building upon a number of estimates of species' exposure to climate published classifications of uncertainty change. At the opposite end of the continuum, sources37,41–44, this thesis uses a classification mechanistic models explicitly include processes that is based on Austin's components. Three such as the species' physiological constraints broad sources of uncertainty are thus distin- and plastic acclimation capacity30, and the guished, affecting either the 'data model', the demographic dynamics31 underlying the intrin- 'mathematical model', and the 'ecological sic vulnerability of species to climate change. model' (Table 1). Between the two ends of the continuum, ap- First, the data used in assessments are im- proaches of varying complexity combine bio- perfect, inherently complex, and often carry the climatic envelope models with available data uncertainties of the models that such data, in that are expected to characterise the vulner- turn, are based upon. Bioclimatic envelope ability of species32. For example, data on the models rely heavily on species' distributional dispersal ability or physiological climatic toler- data. The quality of such data is affected by ance of species have been used to constrain errors in the identification or location of spe- bioclimatic envelope model projections33,34. cies, by geographical, temporal or taxonomical Whereas the ecological information needed to biases in the sampling, and by the interpolation parameterise mechanistic models is lacking for or expert-based techniques used to map distri- most known species, climate and species' oc- butions46–49. By contrast, mechanistic models currence data of varying quality have become rely on detailed physiological or life-history widely available. As a result, the last decades data that are lacking for most species, and that have seen a surge of studies using bioclimatic are typically affected by intra-specific variation, envelope models to assess how species might subjective judgement by experts, inconsistent become exposed to ongoing and future changes terminology, and the techniques used to fill in climate2,35,36. data gaps41,50,51. Not only are species' geo- Irrespective of the approach used, assess- graphical distributions and biology poorly un- ing the potential impacts of climate changes on derstood (a problem termed the 'Wallacean biodiversity is fraught with uncertainty37–41. shortfall'52), most species on Earth are yet to be Assessments are uncertain because of incom- described53,54 (the 'Linnean shortfall'52). As- plete knowledge of the system (epistemic un- sessments of climate change impacts on biodi-
6 | Synopsis
Table 1 | Major sources of uncertainty in assessments of the impacts of future climate change on biodiversity. Examples are given for three major sources of uncertainty associated with the data used, the mathematical model, and the underlying ecological knowledge. The list is not meant to be exhaustive, but to capture some of the major sources that are relevant in the context of this thesis. For each example, the arrows point to the extreme of the correlative (left) to mechanistic (right) model continuum that is most affected, and the roman numerals refer to the chapters of this thesis where the examples are dis- cussed. 'Data model' uncertainty: decisions about the biological and climatic data used in the modelling C ↔ M Chapter Error, bias or subjective judgement in species occurrence data ← Error, bias, or subjective judgement in mechanistic data → II Omission of species → III Error in interpolated baseline climate data ↔ Error in future climate projections ↔ I Unknown future greenhouse gas emissions ↔ I Error in downscaled climate data ↔ 'Mathematical model' uncertainty decisions about the calibration, evaluation and projection methods C ↔ M Chapter Choice of algorithm ← I Methods for predictor selection, data partitioning, absence estimation, etc. ↔ Choice of evaluation data and technique ↔ Small numbers of species' occurrence ← III Extrapolation to non-analogue conditions ← I Choice of threshold for deriving binary projections ← 'Ecological model' uncertainty decisions about ecological concepts and assumptions underpinning the models C ↔ M Chapter Missing or inadequate biological components ← II Missing or inadequate climate parameters and dimensions, and response forms ↔ IV Missing non-climatic factors ↔ Species-climate equilibrium assumption ← Spatial or temporal scale mismatch/ inadequacy ← Lack of clarity about the object and output of models ← II
versity are thus also crippled by omission of (Box 2). Uncertainty from projections of future many of Earth's species. The deficiencies and climate pervades all assessments and is associ- gaps in biological data are not geographically or ated with assumptions regarding the future taxonomically55,56 random, with the tropics and concentrations of greenhouse gases, imperfect threatened species particularly afflicted (Box understanding of the relevant physical proc- 2). esses, errors in the modelling process, and Uncertainty in climate data is associated with a natural variability58. A host of alternative At- multitude of sources affecting both observa- mosphere-Ocean General Circulation Models tions and projections57,58. Characterisation of (hereafter climate models) are available61, with the current climate relies on interpolated data multiple realisations and under different sce- from weather stations or satellites, with areas narios of future emissions of greenhouse poorly covered59 and spatially heterogeneous60 gases62,63. Most such models are developed at showing decreased data reliability and lower coarse resolutions, but climate data are increas- performance of bioclimatic envelope models ingly being downscaled to resolutions that are
Uncertainty in projected impacts of climate change | 7
Box 2 | The data challenge in the tropics Tropical regions contain most of the world’s biodiversity, but are also where the knowledge of species is poorest. The tropics hold the highest rate of discoveries of new mammal species over the past two decades65, with large numbers of further descriptions expected for Africa up to 203265. Tropical species are also more often narrow-ranging than species in temperate zones66, and are thus more often difficult to model with familiar correlative methods requiring sufficient sample sizes67. The resulting bias against narrow-ranging species in assessments may even be larger than is readily apparent, because such species are less likely to be described currently68. Data quality for known species also poses a challenge55,69 in the tropics, and across most of Africa. For vertebrates at the sub-Saharan African scale, the database held at the Zoological Museum in Den- mark70 (Chapter I) is the most comprehensive collation from multiple sources, albeit not free from sampling bias towards the most biodiverse areas71. Most regional efforts towards systematic data col- lection and treatment of uncertainty focus on South Africa or Southern Africa. The data generated through the Southern African Bird Atlas Project72 and analysed with occupancy models73 is a case in point. Knowledge of species' traits is also deficient in the tropics, as illustrated by the lower mean num- ber of data entries in the mammalian dataset PanTHERIA that originate in the tropics56, and the dispro- portionate number of tropical amphibians without data in the Global Amphibian Assessment74. The data challenge in the tropics extends to climate observations and projections. Weather station coverage is poorer in tropical regions75. In comparison to mid-latitudes, the future climates of the trop- ics, and Africa in particular, are still not fully understood due to more complex systems with strong ocean-atmosphere coupling11,75,76. Precipitation projections are generally less certain, with those for tropical Africa subject to especially large decadal variability relative to the mean precipitation77. In turn, temperature increases relative to natural internal variability are expected to be larger in the tropics than mid-latitudes11, although with high uncertainty associated with future emissions of greenhouse gases78. Reasons linked to wealth, communication and security have been suggested to explain deficiencies in species' data79, and possibly also explain poor climate data. Lack of basic data is reflected in a biased fingerprint of ecological effects of climate change4–6, where evidence is sparse for the tropics in general, and for most African countries. Likewise, existing efforts to anticipate future impacts in the tropics are crippled by lack of data, with most published studies focusing on temperate regions80.
more relevant for conservation planning64. Yet, also introducing uncertainty in assessments. In downscaled products are still constrained by climate change and biological invasion applica- the uncertainty of the coarse-scale data that tions in particular, bioclimatic envelope model drive them39,64,81, as well as new uncertainties projections are affected when models are introduced by the methodologies and assump- forced to extrapolate to climatic conditions tions of downscaling techniques82. outside the range of conditions used to cali- Second, the different mathematical models that brate the models84,89. Different algorithms vary are used in impact assessments each also carry in how they handle truncated functions de- their own methodological uncertainties. In scribing the association between climate and bioclimatic envelope modelling, decisions re- suitability for species, causing model outputs to garding the choice of algorithm83,84, estimation vary. Whereas mechanistic models are more of absence data85,86, and thresholds for trans- robust to extrapolation90, they are subject to forming modelled probabilities of suitability other uncertainties associated with the meth- into projections of species' presence and ab- ods and parameters required to handle their sence87,88, for example, all lead to considerable increased ecological complexity41. variation in model outputs. Species with small The third source of uncertainty in assess- numbers of occurrence records are dificult to ments is associated with the underlying eco- model with familiar correlative techniques67, logical decisions. Knowledge of the causal
8 | Synopsis
mechanisms of species' occurrences is incom- synopsis, and discussed in the broader context plete, including knowledge of relevant aspects of the uncertainty challenge in the last section. of climate91,92 described at adequate spatial and Specifically, the following questions are ad- temporal scales64,93,94; non-climatic threats such dressed: as land-use change and landscape fragmenta- tion95; population, ecophysiological, and evolu- i) How can bioclimatic envelope model uncer- tionary dynamics96–98; and biotic interactions99. tainty be integrated in assessments of the While many such mechanisms are being inte- potential impacts of climate change on biodi- grated in models, uncertainty-free projections versity? are unrealistic, not least because of the com- Chapter I104 compares the uncertainty arising plexity and stochasticity inherent to the eco- from alternative climate projections and bio- logical37 and climatic58 systems. On the correla- climatic envelope model algorithms in assess- tive end of the model continuum, the focus is on ments for sub-Saharan African vertebrates, and the patterns – not the mechanisms – in the explores ways to summarise uncertainty by association between observed species' distribu- building consensus around ensembles of mod- tions and climate13. However, uncertainty can els. In Chapter II105, the focus turns to the un- arise when omitted mechanisms cause such certainty of model projections for sub-Saharan models to depart from working assumptions of African amphibians due to ecological assump- equilibrium between species' distributions and tions. The chapter presents an analytical climate15,100. In addition, lack of clarity about framework for teasing apart projected climatic the ecological underpinnings of bioclimatic threats and opportunities facing species and envelope models pervades most practical ap- examining their robustness to assumptions plications21,101,102, partly because there is still about species' vulnerability. debate about the terminology for describing the assumptions, methods and outputs18–22,103. Such ii) Can assessments of the potential impacts uncertainty can lead to inadequate interpreta- of climate change become more inclusive of tion and misuse of outputs beyond their in- Earth's biodiversity despite incomplete tended purpose102. knowledge of species? This thesis investigates some of the above The omission of species with small sample sizes sources of uncertainty (Table 1), and discusses in bioclimatic envelope modelling studies has ways to address them. Focusing on bioclimatic hitherto been largely overlooked as a source of envelope models, the work presented here uncertainty when projecting the effects of cli- explores how projections are influenced by mate change on biodiversity. Chapter III dem- imperfect climate data, omission of species, the onstrates the conservation implications of such choice of model algorithm, and ecological as- bias against narrow-ranging species of sub- sumptions. Much of the analytical component of Saharan African amphibians, underscoring the the thesis is centred on sub-Saharan African importance of complementing bioclimatic enve- vertebrate species, illustrating some of the lope models with alternative tools that are crucial challenges for assessments in tropical scalable to all species. One such alternative tool and sub-tropical areas of high biodiversity and is the use of simple metrics of climate change as poor knowledge (Box 2). Nevertheless, the proxies for the threats and opportunities facing uncertainties addressed are far-reaching and wholesale biodiversity. Climate change metrics the work should have wider relevance. The are reviewed and illustrated at the global scale thesis is composed of five chapters, summa- in Chapter IV, and compared to bioclimatic rised in the following three sections of this envelope models for sub-Saharan African ver- tebrates in Chapter V.
Uncertainty in projected impacts of climate change | 9
Variation in data and model decisions data, and thresholds for binary proj- ections88,91,115–121. Although these studies varied Across the wealth of studies assessing the ef- in methodology as well as geographical and fects of future climate change, the uncertainty taxonomic scope, the choice of model algorithm surrounding model projections is often only was consistently identified as the major cause partially addressed or even overlooked. Previ- of uncertainty whenever included in the com- ous continental-scale applications of biocli- parison (with one exception121). None of these matic envelope models to African species106–110 comparative studies applied to Africa or to have considered only limited variation in future multiple taxonomic groups, a gap that was filled climate models, greenhouse gas emissions by Chapter I. Here too, the comparison of mul- scenarios, or bioclimatic envelope model algo- tiple projections of species' temporal rithms. Yet, the effect of some of these and turnover122 resulting from models built under other sources of uncertainty on model projec- different assumptions identified the model tions is likely to be important, as illustrated for algorithms as the largest relative contributor to sub-Saharan African vertebrates throughout overall uncertainty, followed by the climate this thesis (see Box 3 for an overview of uncer- models and emissions scenarios. tainties for amphibians in particular). The bio- The contribution of the input climate data climatic envelope models used in this thesis to the uncertainty around Chapter I projec- (with the exception of Chapter III) are de- tions was larger in southern areas of sub- scribed in Chapter I104. For over 2,500 species Saharan Africa (Box 3). In previous studies of sub-Saharan African birds, mammals, am- projecting the impacts of climate change, phibians, and snakes71, models were built at one treatment of climate-model uncertainty has degree resolution (≈ 111 km at the Equator), varied between building multiple bioclimatic using baseline data111 and projections81,112 for envelope model projections with a suite of three climatic predictors: mean temperatures available climate models109,114 or, less fre- of both warmest and coldest months, and an- quently, applying a multi-model consensus nual precipitation. Seven different algorithms, projection123. The rationale behind consensus is three multi-model climate projections, and that combining an ensemble of projections, three emissions scenarios were used, generat- assumed to be independent and representative ing an ensemble of projections. of the breadth of possible states of the system Ensembles of models built with alternative being forecast, yields lower mean error than 24 assumptions offer a way to quantify the un- any of the individual projections24,124. Building certainty of model projections. The multiple consensus is challenging not only because the projections in an ensemble are seen as plausi- required assumptions are difficult to ble representations of the system under study, meet58,125,126, but also because there is still de- and together express, at least partially, the bate on the best methodology to combine pro- breadth of uncertainty. Ensemble forecasting is jections126. Averaging across entire ensembles common in climatology58, and has more re- of projections is frequently done, but implies cently started to be used in bioclimatic enve- the loss of information about extreme 109,113,114 lope modelling . Previous studies have states126,127. To allay the criticism levelled at examined ensembles of models to compare the simple averages, Chapter I uses a consensus relative importance of uncertainties arising methodology that first groups co-varying cli- from different sources, including the choice of mate models before summarising the central algorithms, climate models, emissions scenar- tendency for each group113. Among the suite of ios, methods for collinearity correction or vari- 17 climate models, those that were closest to able selection, variable sets, species occurrence
10 | Synopsis
Box 3 | Uncertainties for sub-Saharan African amphibians under climate change Throughout this thesis, different sources of uncertainty are shown to affect projections for sub-Saharan African amphibians under climate change (Fig. 1). Model uncertainties are exemplified by the variation in algorithms (Chapter I). The ecological uncertainty is illustrated by the variation between the raw projections and projections that were modified a posteriori by estimates of species' vulnerability accord- ing to their dispersal ability and tolerance to climate changes (Chapter II). Climate data uncertainty is shown by the variation across three averaged clusters of climate models, and across three emissions scenarios (Chapter I). Finally, data uncertainty arising from the omission of narrow-ranging species (Chapter III) is revealed by comparing two hypothetical projections, one assuming loss of climatic suit- ability for all (omitted) species and the other assuming continued suitability. Chapter I presents a spatial assessment of the relative importance of three sources of uncertainty, based on Analysis of Variance (ANOVA) performed for each grid cell116,118 and using the temporal turno- ver of species122 as response variable. The results show that algorithm-related uncertainty was the larg- est contributor to overall uncertainty, followed by the variation in climate models and emissions scenar- ios. When this assessment was updated here to include all five sources of uncertainty examined in this thesis, the ecological uncertainty as estimated here was found to contribute more to overall uncertainty than the variation across algorithms (median [and lower and upper quartile] proportions of the total sum of squares across the study area of 48.0 [30.2–63.4] and 30.4 [17.3–46.7], respectively). Figure 1 | Uncertainties in projections for African am- phibians under climate change. Bioclimatic envelope models for 263 sub-Saharan African amphibians were used to project changes in the probabil- ity of climatic suitability by mid- century. Projected median abso- lute (positive and negative) changes are shown for a multi- model climate projection under emissions scenario A1B, and reflecting the median across seven algorithms (a). The coeffi- cient of variation is shown across projections for alternative algo- rithms (b), assumptions about the climate change vulnerability of species (c), emissions scenari- os (d), climate models (e), and assumptions about the impacts on narrow-ranging species ex- cluded from the modelling (f). the median model in sign, magnitude and spa- input data for the modelling, an advantage that tial pattern128 were grouped together and aver- is not trivial in the face of the increasing num- aged. The result was a set of three consensus ber of climate projections available61. climate projections for use in the modelling In contrast to climate data uncertainty, the exercise, retaining the variability of the initial variation in Chapter I projections obtained ensemble. At the same time, the consensus with different bioclimatic envelope model algo- solution presented in Chapter I (subsequently rithms was largest in the Sahel and the South- used in another application129) simplifies the ern Sahara. Algorithm-related uncertainty has
Uncertainty in projected impacts of climate change | 11
previously been associated with the type of sions (Chapter I), they may also be more or species' occurrence data used (presence only less robust to assumptions about ecological versus presence-absence data) and the assump- processes. Notably, consideration of biotic tions made during extrapolation to non- interactions133,134, the dispersal ability of spe- analogous, or novel, climates84. In the case of cies135,136, and the fundamental climatic toler- the model projections built in Chapter I, algo- ances of species33,34,137 have all been previously rithm-related uncertainty affected areas where shown to influence model projections under future temperatures were expected to fall climate change. Such processes are usually above the range of calibration, and was thus omitted in classical bioclimatic envelope mod- possibly due to model extrapolation. Model els, but they are crucial when assessing biodi- projections for these areas are less reliable, as versity risk under climate change. When direct are likely to be any projections in other tropical incorporation of ecological data in mechanistic and sub-tropical regions where emergence of models is not viable, bioclimatic envelope non-analogous climates is expected130 (see model projections can be contrasted with alter- Chapter IV). native assumptions about the vulnerability of The variation in outputs from different al- species138,139. Exposure and vulnerability esti- gorithms has sometimes been summarised by mates are often also combined into indices that building consensus projections109,113,114, but are meant to describe the climate change risk bioclimatic modellers, like climatologists, are facing species138,140. Chapter II105 presents an still debating what is the best methodology to analytical framework to examine projections of combine projections. Chapter I extends previ- changes in climatic suitability with reference to ous comparisons of consensus meth- estimates of species' vulnerability to such odologies24,131,132 by exploring a comprehensive changes. In contrast to previous studies com- suite of five alternative methodologies that paring exposure to vulnerability, Chapter II either measure the central tendency in the argues that the different extrinsic threats and frequency distribution of all projections (mean, opportunities arising from the exposure of weighted mean, and median), measure the species to climate change should be examined central tendency among groups of co-varying separately, and matched to specific intrinsic projections (as above for climate models), or traits that are likely to mediate species' re- select the projection summarising the highest sponses. amount of variation among projections (using In published studies using bioclimatic en- Principal Components Analysis). These five velope models, the degree of species' exposure methodologies were compared by examining to climate change is commonly inferred from their accuracy in predicting the baseline distri- temporal changes in the total size of areas with bution of species. Consensus projections for suitable climatic conditions107,139,141. Species are sub-Saharan African vertebrates were generally considered 'winners' or 'losers' depending on more accurate – and more similar to each other whether such areas become more or less avail- – than individual projections, with the consen- able across a given region. However, such sus method based on co-varying groups yield- summary measures conceal different climatic ing the highest accuracy. opportunities as well as threats, each imposing specific constraints on species28. At the level of each grid cell, climatic suitability may be lost, Robustness to ecological assumptions gained, or become markedly different from neighbour cells, with consequences not only for Not only are bioclimatic envelope model out- the size, but also the position and level of frag- puts sensitive to data and methodological deci-
12 | Synopsis
Box 4 | Teasing apart climatic threats and opportunities The threats and opportunities arising from exposure of species to climate change are normally inferred from bioclimatic envelope models (Chapter II), but can also be gauged with cli- mate change metrics (Chapters IV and V). With the former, the geograph- ical projection of species' bioclimatic envelope is compared across time, and changes in area occupied, distance shifted, and degree of fragmentation can be measured. In turn, climate change metrics depict the spatio- temporal dynamics of climate across broad regions, independently of spe- cies' data. Metrics can describe chang- es in the area occupied by analogous climatic conditions, or in the distance to such conditions, allowing for indi- rect inferences to be made about the potential implications for the biodi- versity in such regions (Chapter IV). Different combinations of the threats and opportunities, measured with models or metrics, have distinct con- servation implications (Fig. 2).
Figure 2 | Climatic threats and opportunities for species. Geographical representation of changes in areas of suitable climate for species between two time steps, with reference to different scenarios of threats and opportunities from future climate change inferred from bioclimatic envelope models or measured changes in climate (light and dark blue scales respectively, see Chapter V). Combinations of threats and opportunities relating to both the size and position of areas with given climatic conditions result in different risks for species occupying such areas (panels (a-d)). Using size as sole criterion would lead to classifying species (a) and (b) as 'winners', and species (c) and (d) as 'losers'. However, examining also the changes in position reveals that the future suitable areas for species (b) and (d) are displaced, leaving no overlap with the baseline areas. As a result, species (b) and (d) should also be flagged, particularly if they are poor dispersers. These species thus require the protection of new areas and enhanced connectivity to reach them, whereas species (a) and (c) depend on protection in situ. The axis for area could be decomposed into losses and gains, and a third axis could be added for the threat of increased fragmentation.
mentation of species' bioclimatic envelopes areas of climatic suitability. Previous frame- across broader regions (Box 4). These compo- works28,34 separated the threat of loss from nents of exposure are likely to have distinct opportunities for gains. Chapter II adds spatial distributions and implications for con- fragmentation, a threat to the persistence143 servation guidance142, but are seldom teased and expansion144 of species' ranges that has apart. Using the bioclimatic envelope models hitherto been largely overlooked in climate built in Chapter I for sub-Saharan African am- change ecology (but see a recent example145). phibians as illustration, Chapter II teases apart The different threats and opportunities are projected losses, fragmentation and gains of also likely to interact with specific traits. The
Uncertainty in projected impacts of climate change | 13
specificity in the interaction between extrinsic analysis were the narrowest-ranging species in threats and intrinsic vulnerability has been the database, which were not modelled in previously shown for non-climatic threats such Chapter I (see Chapter III). Indeed, the analy- as habitat destruction or hunting146–149, but has sis presented in Chapter II for sub-Saharan not been sufficiently addressed in climate African amphibians illustrates the data uncer- change ecology. Risk indices used to inform tainty that is typical of tropical regions (Box 2). priorities for conservation under future cli- Not only did model projections exclude narrow- mates150, for example, frequently combine ranging species, but data for traits in the strict summary measures of exposure with summary sense were also unavailable and had to be re- estimates of vulnerability according to available placed with proxies, carrying their own uncer- trait data, irrespective of the potential specific tainty to the analysis. interactions between the two. The analytical framework proposed in Chapter II aims to select specific response-mediating traits for Broadening the scope of assessments each climatic threat and opportunity, based on The uncertainties discussed in Chapters I and theoretical and empirical research on climate II applied to projections of the exposure of sub- change vulnerability151 (see also collated refer- Saharan African vertebrates to future climate ences in Chapter II). Thus, for example, pro- changes. As already alluded to above, such jected losses of climatic suitability can be as- projections were not available for all species. sessed for their robustness to traits influencing Excluded were species with fewer than 15 one- the species' physiological climatic tolerance, degree gridded occurrence records, considered whereas projected gains of newly suitable areas too uncertain to model due to their small sam- can be contrasted with traits associated with ple size67. The modelling thus covered only the dispersal ability or reproductive output of 67% of the 4,092 vertebrate species – and a species. mere 38% of the 741 amphibian species – in For the sub-Saharan African amphibians the African vertebrate database held at the analysed in Chapter II, spatially overlaying Zoological Museum in Denmark70. Far from projected losses, fragmentation or gains of unique to the studies presented in this thesis, climatically suitable areas with vulnerability the bias against narrow-ranging species per- classifications derived from selected trait vades most modelling exercises under climate data140,152 illustrates how interpretation of change. This often overlooked source of uncer- outputs from bioclimatic envelope models can tainty is the focus of Chapter III and a motiva- be altered with consideration of species' cli- tion behind Chapters IV and V, where an alter- mate change vulnerability (Box 3). The differ- native approach is proposed to broaden the ences revealed underscore the need to clearly scope of assessments. communicate in assessments what is being Using sub-Saharan African amphibians as a modelled and what the outputs mean. The case study, Chapter III demonstrates the impli- modifying effect of vulnerability on model pro- cations that omitting narrow-ranging species in jections was stronger for the species in the projections of the impacts of climate change can study sample that were most exposed to cli- have for conservation priorities. The analysis is mate change. For example, for the amphibians based on a new modelling exercise, but the with the smallest geographical overlap between results presented hold for the models of Chap- baseline and future climate envelopes, pro- ter I. The bioclimatic envelope modelling per- jected gains may be unrealisable due to poor formed in Chapter III left out the species with dispersal ability (Box 4). Excluded from the fewer than 10 one-degree gridded occurrence
14 | Synopsis
records, amounting to over half of all species in existing schemes and priorities derived for the the dataset and including the vast majority of modelled species decreased into the future, threatened species. Multivariate ordination suggesting that the use of bioclimatic envelope analysis showed that the subset of modelled models to inform conservation priority setting species is not representative of either the pre- may systematically downplay important areas sent climate space of all known sub-Saharan for rare and threatened species. Irrespective of African amphibians or their future exposure to the actual impacts of ongoing climate changes climate changes. In comparison to modelled on such rare species, most of them face other species, those that were omitted occupy topog- non-climatic, more immediate threats that raphically complex areas with cooler, wetter warrant priority in conservation strategies163. and less seasonal climates, which are projected The results of Chapter III underscore the to experience lower climate anomalies. Previ- importance of complementing models with ous studies also found that climatically and alternative, more generally applicable tools. topographically diverse or distinct regions tend One such alternative tool is the use of simple to contain a disproportionate richness of nar- metrics of climate change, which quantify the row-ranging species153,154. Looking into the exposure of geographical areas to changes in future (Box 3), topographical diversity may climate parameters. Metrics have previously confer climatic stability155 and buffer the effects helped address a diversity of ecological ques- of broader scale climate change93. Yet, other tions such as the potential risks from future dimensions of climate change, including the climate changes to biodiversity130,153,164 and risk of disappearing climates, may especially conservation areas165–167, suggesting that they affect species with narrow climatic tolerances can be a useful tool in impact assessments168,169. such as those omitted from the modelling130,153 However, the variety of existing metrics, and (see below for a discussion of the different their ecological implications, have hitherto not dimensions of climate change reviewed in been fully appreciated. Chapter IV reviews Chapter IV). climate change metrics and classifies them into Bioclimatic envelope models have sup- metrics of local change and regional change. ported analyses of the role of existing conserva- Climate anomalies are the most commonly used tion schemes under changing climates156,157. metric, depicting average changes in climate at However, the range size bias introduced in such the locality (grid cell) level130. Other local met- analyses is counter to traditional frameworks rics capture changes in the magnitude of ex- for identifying priority sites for conservation, treme climates164, the timing of climatic which focus precisely on those taxa most likely events170, or the local velocity of climate167. to be omitted from the modelling158,159. Previ- In turn, regional metrics characterise ous studies160 have shown that including in- changes in the distribution of climatic condi- formation about unknown or poorly sampled tions over broader geographical areas (sets of species influences the selection of reserves. In cells). These metrics begin with the characteri- the case of the sub-Saharan African amphibians sation of climatic conditions across a given examined in Chapter III, conservation priori- region, and can then describe temporal changes ties derived using the modelled species subset in the availability of analogous climatic condi- exhibited lower spatial congruence with exist- tions across that region130,171,172, as well as ing conservation priority schemes (Biodiversity changes in the direction to, or distance be- Hotspots161, Endemic Bird Areas158 and The tween, the positions of analogous conditions171. Global 200162) than those inferred from the Chapter IV highlights the contrasting spatial omitted species subset. Congruence between patterns arising at the global scale from differ-
Uncertainty in projected impacts of climate change | 15
ent climate change metrics. Whereas polar or even extinctions, of species' ranges over regions face reductions in the global availability time182. Such changes often match, with time of similar climatic conditions, the tropics and lag, the increases and decreases in the area of hot arid regions are exposed to average different climatic conditions brought about by changes beyond historical inter-annual vari- the Earth's warming and cooling cycles183–185. At ability, emergence of novel climates, and in- the same time, measured shifts in the position creased frequency of extreme climates. of climatic conditions have been used to explain To help interpret the diversity of metrics, the opportunity for species to track climatic Chapter IV proposes a conceptual framework conditions across broad regions in both palaeo- for classifying metrics into common currencies cological185 and recent3,5 times. of climatic threat and opportunity for species. To test the feasibility of the framework The framework is based on principles linking proposed in Chapter IV for guiding the use of the persistence of populations and species to metrics, Chapter V compares climate change local and regional climatic suitability, respec- metrics for sub-Saharan Africa with bioclimatic tively173. It builds upon published empirical envelope models built in Chapter I for verte- studies associating different aspects of climate brates. For all taxa, local climate anomalies, as change with observed ecological changes in projected by the metrics, were diagnostic of recent or palaeocological time, to establish changes in climatic suitability of grid cells for potential links between climate change metrics species, as projected by the models. In turn, and threats and opportunities for biodiversity. changes in the area and position of analogous Some of the empirical associations previously climates, as measured by the metrics, were found are highlighted here, and further de- indicative of changes in the size and position of scribed in Chapter IV. species' bioclimatic envelopes, respectively. At the local level, decreased climatic suit- The agreement found between metrics and ability can have effects on the physiology, mor- models was stronger for the vertebrate species phology or behavior of the organisms in a with narrower climatic breadths, which depend population174, potentially leading to changes in more strongly on tracking climates130,153 than demography. Reduced activity of individuals, species with more generalist climate prefer- population declines, and mortality have been ences. Nevertheless, the results in Chapter V linked to metrics of gradual175,176 or extreme demonstrate that, when carefully implemented changes177 in local climate. Positive cases of and interpreted, metrics can provide first ap- increased fitness or abundance have also been proximations to the same extrinsic threats and reported178,179, although less frequently. In turn, opportunities typically captured by bioclimatic observed shifts in the timing of seasonal activi- envelope models. ties of populations have been associated with measured changes in the seasonality of climate180. Where dispersal is involved, the Embracing uncertainty velocity of temperature changes has been sug- This thesis illustrates a suite of critical sources gested to pace the rate of population shifts of uncertainty in assessments of biodiversity across the local topography181. impacts under changing climates. Unknown At the regional level, the spatio-temporal future climates (Chapter I), omission of species dynamics of climate can affect the availability (Chapter III), the diversity of model algorithms and distribution of climatically suitable areas (Chapter I), and species-specific vulnerability for species173 (see also Box 4). Long-term eco- to climate change (Chapter II), are shown to logical data reveal expansions and contractions, affect bioclimatic envelope model projections
16 | Synopsis
for sub-Saharan African vertebrates. These There are many other ways to integrate uncer- uncertainties, and many others not covered tainty in assessments, which are not discussed here, impair the conservation of biodiversity in this thesis. Comparing the outputs from dif- under changing climates in Africa and else- ferent models along the correlative-mechanistic where. Despite numerous calls to explicitly continuum, each with their own uncertainties, address this challenge39,41,47,83,84,131, appropriate can help highlight areas of convergence and treatment of uncertainty has yet to be formal- offset uncertainties189–191. Precision estimates ised in impact assessments and conservation for the data used in bioclimatic envelope mod- adaptation. The chapters in this thesis highlight els can be visualised in maps accompanying the need to both integrate uncertainties in model projections47, or incorporated directly in assessments, and reduce or circumvent them the modelling process192,193. For the species' where possible. distributional dataset used in Chapter I and Firstly, precise projections are unattain- many other datasets, explicit information on able, and existing uncertainties should be in- the sampling effort associated with species' corporated in assessments and quantified to occurrence records is lacking. But even in such the extent possible41,186. In answering the first cases, alternative approaches based on qualita- question posed at the outset, this thesis shows tive, objective differences in survey intensity examples of how bioclimatic envelope model across broad, nested regions49 can be used to uncertainty can be integrated in assess- detect potential biases. ments. Building ensembles of models with Secondly, understanding the uncertainty in alternative assumptions is illustrated in Chap- projections helps to direct efforts to reduce it ter I. Here, as in other published studies, the where most needed. Reducing uncertainty is an use of ensembles is constrained by the breadth integral part of ecology and any science. Models of possible future states sampled, but can, at are 'educated guesses about the future'39, and least, give a lower bound on the range of uncer- increasing the precision (and accuracy) of the tainty. Probabilistic forecasting, where site- guesses is a key goal in climate change ecology. specific probability density functions for mod- Narrowing the guesses based on bioclimatic elled suitability account for several uncertain- envelope models rests partly on improving the ties, has only rarely been applied in bioclimatic data on species' occurrences. Efforts towards envelope modelling187,188. When ensembles are enhanced, targeted sampling194 are warranted used to build consensus projections, these particularly in tropical regions (Box 2), where should preserve information about uncer- possible combining observational, experimen- tainty41, for example by relying on methodolo- tal and mathematical lines of research25. En- gies that capture the breadth of the ensemble hancing the precision of future climate projec- (Chapter I). tions, in turn, is seen by climatologists195 as Chapter II gives another example of how bringing smaller benefits than understanding uncertainty can be integrated in assessments, the influence of existing uncertainties on pro- by examining a posteriori the effect of estimates jected impacts. In fact, despite continuous de- of species' vulnerability on model projections. velopment of climate models and greater com- The proposed framework carries several sim- putational capacity since the Fourth Assess- plifying assumptions that require further scru- ment Report of the Inter-governmental Panel tiny. Nevertheless, its application can reveal on Climate Change, model spread has not where, in geographical space, the model projec- changed much with the new generation of tions may be less robust to crucial ecological models195–197. factors not contemplated in the modelling.
Uncertainty in projected impacts of climate change | 17
More important are perhaps any efforts to be influential201. Yet, increased realism comes reduce the ecological uncertainty arising from at the expense of applicability to many species. the choice of climatic predictors91. Assessments Mechanistic models have only been empirically ignoring influential climate parameters (e.g., tested on few species with available data, and temperature, precipitation) have previously even simple bioclimatic envelope models are been shown to misjudge potential biotic re- not suited to the narrowest-ranging species sponses92. Equally crucial is to measure the (Chapter III). What is more, both correlative relevant dimensions of change in selected pa- and mechanistic models are blind to the major- rameters (e.g., means, extremes, or timing of ity of species on Earth, still unknown54. climatic events)198, an issue that was empha- The second question posed in this thesis sised in Chapter IV by the diversity of spatial asked how assessments can become more patterns, and their potential links to biotic inclusive of Earth's biodiversity. Describing threats, across dimensions. Changes in climate the distributions and auto-ecology of species, extremes, shifts in the timing of climatic events, and modelling their responses, must remain a and climate change velocity, for example, are priority, but gathering data for all species on less frequently considered in bioclimatic enve- Earth is impractical. What is suggested here is lope modelling exercises, but can be integrated that investment in models be accompanied by in the models themselves198 or complement the development of alternative approaches that assessments. Even in large scale, multi-species circumvent the uncertainty associated with modelling exercises such as that in Chapter I, species' data. Climate change metrics are one more ecologically sound choices can be made candidate approach (Chapter IV). Early testing by first grouping species that are likely to be suggests that metrics can provide inferences influenced by similar predictors. What then about extrinsic climatic threats and opportuni- remains a key challenge – in Chapter I and ties in agreement with bioclimatic envelope most other applications199 – is the mismatch in model projections (Chapter V), and metrics spatial scale between available climate predic- should now be compared with available time- tors and species' occurrences64. Whereas coarse series species' data. The use of metrics does not and uncertain global or regional climate models alleviate the need for considering the intrinsic may be useful to inform broad strategies, con- vulnerability of species, but can give first indi- servation planning at fine scales199 may require cations on the potential extrinsic threats facing greater certainty from models of topographic biodiversity. Clearly, providing conservation microclimate200. planners with assessments that are useful re- In an effort to further reduce the ecological quires the integration of multiple approaches. uncertainty of projections, the field of climate In closing, it must be acknowledged that change ecology is becoming more integrated, explicitly addressing all uncertainties of pro- harnessing new disciplines and types of jected impacts appears overwhelming. Yet, if data25,182 to describe the potential behaviour of model projections are to help manage the fu- the ecological system. Complex models are ture impacts of climate change on Earth's bio- valuable for guiding species-based conserva- diversity, they must be as transparent as possi- tion, while also advancing ecological theory. ble by including an honest description of their Indeed, the new uncertainties carried by com- abilities and limitations. Such description must plex models, and even by simplified frame- guide end-users through the ecological mean- works such as that presented in Chapter II, can ing and purpose of the models, and reveal the stimulate new research to identify the circum- level of confidence associated with the outputs stances under which different mechanisms may (Box 3). Projections thus presented can still be
18 | Synopsis
useful to guide decision-making, while avoiding declines among African savanna ungulates. Ecology Letters 6, 412–419 (2003). the risks of over-interpretation and loss of 2. Nabout, J. C. et al. Trends and Biases in Global credibility should additional data or new mod- Climate Change Literature. Natureza & els contradict them in the future. In a rapidly Conservação 10, 45–51 (2012). evolving field as climate change ecology, such 3. Parmesan, C. & Yohe, G. A globally coherent risks are real and could undermine our efforts fingerprint of climate change impacts across natural systems. Nature 421, 37–42 (2003). to adapt to climate change. Ultimately, factoring 4. Parmesan, C. Ecological and Evolutionary 202 uncertainties into decision-making requires Responses to Recent Climate Change. Annual collaboration with social and political scien- Review of Ecology, Evolution, and Systematics 37, tists. Equipped with ensembles of projections, 637–669 (2006). model users have the option of weighing the 5. Chen, I.-C., Hill, J. K., Ohlemüller, R., Roy, D. B. & Thomas, C. D. Rapid Range Shifts of Species risks associated with alternative decisions39,195. Associated with High Levels of Climate Warming. Risk-averse users may opt for the most pessi- Science 333, 1024–1026 (2011). mistic scenario, whereas the middle-term or 6. Rosenzweig, C. et al. Attributing physical and biological impacts to anthropogenic climate the consensus projection may be adequate in a change. Nature 453, 353–357 (2008). 203 risk-tolerant context . When planning for a 7. Root, T. L., MacMynowski, D. P., Mastrandrea, M. D. range of outcomes is viable, assessments under & Schneider, S. H. Human-modified temperatures multiple, plausible hypotheses about future induce species changes: Joint attribution. Proceedings of the National Academy of Sciences of impacts can lead to better conservation deci- the United States of America 102, 7465–9 (2005). 204 sions and may be preferable. 8. IPCC. Climate Change 2001: The Scientific Basis. The picture that emerges from decades of Contribution of Working Group I to the Third climatic research is one of drastic and multifac- Assessment Report of the Intergovernmental Panel on Climate Change. (Cambridge University eted changes in climate. There is still much to Press, 2001). understand about the climate, the distributions 9. IPCC. Climate change 2007: The physical science and biology of species, and the ability of species basis. Contribution of working Group I to the Fourth Assessment Report of the to cope with climate changes and interacting Intergovernemental Panel on Climate Change. threats163,205. In the face of all uncertainties, (Cambridge University Press, 2007). protecting biodiversity is a daunting task. But 10. Peters, G. P. et al. The challenge to keep global are imprecise projections flawed beyond use? warming below 2°C. Nature Climate Change 3, 4–6 (2013). Uncertainty needs to be understood, yet it 11. IPCC. Summary for Policymakers. in Climate should not deter our efforts to anticipate the Change 2013: The Physical Science Basis. impacts of climate change on biodiversity39,195. Contribution of Working Group I to the Fifth Rather, the projections that are attainable with Assessment Report of the Intergovernmental Panel on Climate Change (Stocker, T. F. et al.) the current knowledge should help direct fur- (Cambridge University Press, 2013). 206,207 ther research to priority regions , and 12. Barnosky, A. D. Heatstronke - Nature in an age of prompt us to adapt conservation strategies to Global Warming. 269 (Island Press, 2009). account for uncertain impacts. Sacrificing some 13. Dormann, C. F. et al. Correlation and process in precision to generality may well be needed if species distribution models: bridging a dichotomy. Journal of Biogeography 39, 2119–2131 (2012). we are to rise to the challenge of protecting 14. Guisan, A. & Zimmermann, N. E. Predictive habitat Earth's biodiversity. distribution models in ecology. Ecological Modelling 135, 147–186 (2000). 15. Peterson, A. T. et al. Ecological Niches and Geographic Distributions. (Monographs in References Population Biology, Princeton University Press, 2011). 1. Ogutu, J. O. & Owen-Smith, N. ENSO, rainfall and temperature influences on extreme population 16. Grinnell, J. The Niche-Relationships of the California Thrasher. The Auk 34, 427–433 (1917).
Uncertainty in projected impacts of climate change | 19
17. James, F. C., Johnston, R. F., Wamer, N. O., Niemi, G. 33. Elith, J., Kearney, M. & Phillips, S. The art of J. & Boecklen, W. J. No TitleThe Grinnellian Niche modelling range-shifting species. Methods in of the Wood Thrush. The American Naturalist 124, Ecology and Evolution 1, 330–342 (2010). 17–47 (1984). 34. Arribas, P. et al. Evaluating drivers of vulnerability 18. Kearney, M. Habitat, environment and niche: what to climate change: a guide for insect conservation are we modelling? Oikos 115, 186–191 (2006). strategies. Global Change Biology 18, 2135–2146 (2012). 19. Elith, J. & Leathwick, J. R. Species Distribution Models: Ecological Explanation and Prediction 35. Siqueira, T., Padial, A. A. & Bini, L. M. Mudanças Across Space and Time. Annual Review of Ecology, climáticas e seus efeitos sobre a biodiversidade: Evolution, and Systematics 40, 677–697 (2009). um panorama sobre as atividades de pesquisa. Megadiversidade 5, 17–26 (2009). 20. Warren, D. L. In defense of “niche modeling.”Trends in Ecology & Evolution 27, 497– 36. Guisan, A. et al. Predicting species distributions for 500 (2012). conservation decisions. Ecology Letters 16, 1424– 1435 (2013). 21. McInerny, G. J. & Etienne, R. S. Ditch the niche – is the niche a useful concept in ecology or species 37. Elith, J., Burgman, M. A. & Regan, H. M. Mapping distribution modelling? Journal of Biogeography epistemic uncertainties and vague concepts in 39, 2096–2102 (2012). predictions of species distribution. Ecological Modelling 157, 313–329 (2002). 22. Hortal, J., Lobo, J. M. & Jiménez-Valverde, A. Basic Questions in Biogeography and the (Lack of) 38. Barry, S. & Elith, J. Error and uncertainty in habitat Simplicity of Species Distributions: Putting Species models. Journal of Applied Ecology 43, 413–423 Distribution Models in the Right Place. Natureza & (2006). Conservação 10.2, 108–118 (2012). 39. Wiens, J. A., Stralberg, D., Jongsomjit, D., Howell, C. 23. Anderson, R. P. Harnessing the world’s A. & Snyder, M. A. Niches, models, and climate biodiversity data: promise and peril in ecological change: Assessing the assumptions and niche modeling of species distributions. Annals of uncertainties. Proceedings of the National Academy the New York Academy of Sciences 1260, 66–80 of Sciences 106, 19729–19736 (2009). (2012). 40. Heikkinen, R. K. et al. Methods and uncertainties in 24. Araújo, M. B. & New, M. Ensemble forecasting of bioclimatic envelope modelling under climate species distributions. Trends in Ecology & change. Progress in Physical Geography 30, 751– Evolution 22, 42–47 (2007). 777 (2006). 25. Dawson, T. P., Jackson, S. T., House, J. I., Prentice, I. 41. Beale, C. M. & Lennon, J. J. Incorporating C. & Mace, G. M. Beyond predictions: biodiversity uncertainty in predictive species distribution conservation in a changing climate. Science 332, modelling. Philosophical transactions of the Royal 53–8 (2011). Society of London. Series B, Biological sciences 367, 247–58 (2012). 26. Gaston, K. J. Rarity. (Chapman and Hall, 1994). 42. Regan, H. M., Colyvan, M. & Burgman, M. A. A 27. Araújo, M. B. & Williams, P. H. Selecting areas for taxonomy and treatment of uncertainty for species persistence using occurrence data. ecology and conservation biology. Ecological Biological Conservation 96, 331–345 (2000). Applications 12, 618–628 (2002). 28. Thomas, C. D. et al. A framework for assessing 43. Ascough II, J. C., Maier, H. R., Ravalico, J. K. & threats and benefits to species responding to Strudley, M. W. Future research challenges for climate change. Methods in Ecology and Evolution incorporation of uncertainty in environmental and 2, 125–142 (2011). ecological decision-making. Ecological Modelling 29. Violle, C. et al. Let the concept of trait be 219, 383–399 (2008). functional! Oikos 116, 882–892 (2007). 44. Kujala, H., Burgman, M. A. & Moilanen, A. 30. Kearney, M. & Porter, W. Mechanistic niche Treatment of uncertainty in conservation under modelling: combining physiological and spatial climate change. Conservation Letters 6, 73–85 data to predict species’ ranges. Ecology letters 12, (2012). 334–50 (2009). 45. Austin, M. P. Spatial prediction of species 31. Keith, D. A. et al. Predicting extinction risks under distribution: an interface between ecological climate change: coupling stochastic population theory and statistical modelling. Ecological models with dynamic bioclimatic habitat models. Modelling 157, 101–118 (2002). Biology Letters 4, 560–563 (2008). 46. Hortal, J., Jiménez-Valverde, A., Gómez, J. F., Lobo, 32. Willis, S. G. et al. A framework methodology for J. M. & Baselga, A. Historical bias in biodiversity integrating climate change vulnerability inventories affects the observed environmental assessments from Species Distribution Models and niche of the species. Oikos 117, 847–858 (2008). Trait-based Approaches for use in systematic conservation planning. In preparation.
20 | Synopsis
47. Rocchini, D. et al. Accounting for uncertainty when Bulletin of the American Meteorological Society 93, mapping species distributions: The need for maps 485–498 (2011). of ignorance. Progress in Physical Geography 35, 62. Nakicenovic, N. & Swart, R. Special Report on 211–226 (2011). Emissions Scenarios: A Special Report of Working 48. Kujala, H., Vepsäläinen, V., Zuckerberg, B. & Group III of the Intergovernmental Panel on Brommer, J. E. Range margin shifts of birds Climate Change. 599 (Cambridge University Press, revisited - the role of spatiotemporally varying 2000). survey effort. Global change biology 19, 420–30 63. Moss, R. H. et al. The next generation of scenarios (2013). for climate change research and assessment. 49. Barbosa, A. M., Pautasso, M. & Figueiredo, D. Nature 463, 747–56 (2010). Species–people correlations and the need to 64. Wiens, J. A. & Bachelet, D. Matching the Multiple account for survey effort in biodiversity analyses. Scales of Conservation with the Multiple Scales of Diversity and Distributions 19, 1188–1197 (2013). Climate Change. Conservation Biology 24, 51–62 50. Nakagawa, S. & Freckleton, R. P. Missing inaction: (2010). the dangers of ignoring missing data. Trends in 65. Ceballos, G. & Ehrlich, P. R. Discoveries of new ecology & evolution 23, 592–596 (2008). mammal species and their implications for 51. Fitzsimmons, J. M. How consistent are trait data conservation and ecosystem services. Proceedings between sources? A quantitative assessment. of the National Academy of Sciences 106 , 3841– Oikos 122, 1350–1356 (2013). 3846 (2009). 52. Lomolino, M. V. Conservation Biogeography. in 66. Wiens, J. J., Graham, C. H., Moen, D. S., Smith, S. a & Frontiers of Biogeography: new directions in the Reeder, T. W. Evolutionary and ecological causes geography of nature (Lomolino, M. V. & Heaney, L. of the latitudinal diversity gradient in hylid frogs: R.) 293–296 (Sinauer Associates, Inc., 2004). treefrog trees unearth the roots of high tropical diversity. The American naturalist 168, 579–96 53. May, R. M. & Beverton, R. J. H. How Many Species? (2006). [and Discussion]. Philosophical Transactions of the Royal Society of London. Series B: Biological 67. Stockwell, D. R. B. & Peterson, A. T. Effects of Sciences 330 , 293–304 (1990). sample size on accuracy of species distribution 54. Scheffers, B. R., Joppa, L. N., Pimm, S. L. & models. Ecological Modelling 148, 1–13 (2002). Laurance, W. F. What we know and don’t know 68. Patterson, B. D. Accumulating Knowledge on the about Earth's missing biodiversity. Trends in Dimensions of Biodiversity: Systematic Ecology & Evolution 27, 501–510 (2012). Perspectives on Neotropical Mammals. Biodiversity Letters 2, 79–86 (1994). 55. Feeley, K. J. & Silman, M. R. The data void in modeling current and future distributions of 69. Kamino, L. H. Y. et al. Challenges and perspectives tropical species. Global Change Biology 17, 626– for species distribution modelling in the 630 (2011). neotropics. Biology letters 8, 324–6 (2012). 56. González-Suárez, M., Lucas, P. M. & Revilla, E. 70. Burgess, N., Fjeldså, J. & Rahbek, C. Mapping the Biases in comparative analyses of extinction risk: distributions of Afrotropical vertebrate groups. mind the gap. Journal of Animal Ecology 81, 1211– Species 30, 16–17 (1998). 1222 (2012). 71. Hansen, L. A. Personal communication. (2013). 57. Beaumont, L. J., Hughes, L. & Pitman, A. J. Why is 72. Harrison, J. A., Underhill, L. G. & Barnard, P. The the choice of future climate scenarios for species seminal legacy of the Southern African Bird Atlas distribution modelling important? Ecology letters project. South African Journal of Science 104, 82– 11, 1135–46 (2008). 84 (2008). 58. Knutti, R. et al. A Review of Uncertainties in Global 73. Bled, F., Nichols, J. D. & Altwegg, R. Dynamic Temperature Projections over the Twenty-First occupancy models for analyzing species’ range Century. Journal of Climate 21, 2651–2663 (2008). dynamics across large geographic scales. Ecology 59. Hijmans, R. J., Cameron, S. E., Parra, J. L., Jones, P. and Evolution 3, 4896–4909 (2013). G. & Jarvis, A. Very high resolution interpolated 74. IUCN Species Survival Commission, Conservation climate surfaces for global land areas. Science, International Center for Applied International Journal of Climatology 25, 1965– Biodiversity, NatureServe & IUCN. Global 1978 (2005). Amphibian Assessment. (2004). at 60. Fernández, M., Hamilton, H. & Kueppers, L. M.
Uncertainty in projected impacts of climate change | 21
76. James, R. & Washington, R. Changes in African 90. Keenan, T., Serra, J. M., Lloret, F., Ninyerola, M. & temperature and precipitation associated with Sabate, S. Predicting the future of forests in the degrees of global warming. Climatic Change 117, Mediterranean under climate change, with niche- 859–872 (2013). and process-based models: CO2 matters! Global Change Biology 17, 565–579 (2011). 77. Hawkins, E. & Sutton, R. The potential to narrow uncertainty in projections of regional 91. Synes, N. W. & Osborne, P. E. Choice of predictor precipitation change. Climate Dynamics 37, 407– variables as a source of uncertainty in continental- 418 (2011). scale species distribution modelling under climate change. Global Ecology and Biogeography 20, 904– 78. Hawkins, E. & Sutton, R. The Potential to Narrow 914 (2011). Uncertainty in Regional Climate Predictions. Bulletin of the American Meteorological Society 90, 92. VanDerWal, J. et al. Focus on poleward shifts in 1095–1107 (2009). species’ distribution underestimates the fingerprint of climate change. Nature Climate 79. Amano, T. & Sutherland, W. J. Four barriers to the Change 3, 239–243 (2013). global understanding of biodiversity conservation: wealth, language, geographical location and 93. Logan, M. L., Huynh, R. K., Precious, R. A. & security. Proceedings of the Royal Society B: Calsbeek, R. G. The impact of climate change Biological Sciences 280, 20122649 (2013). measured at relevant spatial scales: new hope for tropical lizards. Global Change Biology 19, 3093– 80. Felton, A. et al. Climate change, conservation and 3102 (2013). management: an assessment of the peer-reviewed scientific journal literature. Biodiversity and 94. Jackson, S. T., Betancourt, J. L., Booth, R. K. & Gray, Conservation 18, 2243–2253 (2009). S. T. Ecology and the ratchet of events: Climate variability, niche dimensions, and species 81. Tabor, K. & Williams, J. W. Globally downscaled distributions. Proceedings of the National Academy climate projections for assessing the conservation of Sciences 106, 19685–19692 (2009). impacts of climate change. Ecological Applications 20, 554–565 (2010). 95. Brook, B. W., Sodhi, N. S. & Bradshaw, C. J. A. Synergies among extinction drivers under global 82. Pierce, D. W. et al. Probabilistic estimates of future change. Trends in Ecology & Evolution 23, 453– changes in California temperature and 460 (2008). precipitation using statistical and dynamical downscaling. Climate Dynamics 40, 839–856 96. Holt, R. D. The microevolutionary consequences of (2013). climate change. Trends in Ecology & Evolution 5, 311–315 (1990). 83. Thuiller, W. et al. Biodiversity conservation: Uncertainty in predictions of extinction risk. 97. Chevin, L.-M., Lande, R. & Mace, G. M. Adaptation, Nature 430, 145–148 (2004). Plasticity, and Extinction in a Changing Environment: Towards a Predictive Theory. PLoS 84. Pearson, R. G. et al. Model-based uncertainty in Biology 8, e1000357 (2010). species range prediction. Journal of Biogeography 33, 1704–1711 (2006). 98. Hill, J. K., Griffiths, H. M. & Thomas, C. D. Climate Change and Evolutionary Adaptations at Species’ 85. Chefaoui, R. M. & Lobo, J. M. Assessing the effects Range Margins. Annual Review of Entomology 56, of pseudo-absences on predictive distribution 143–159 (2011). model performance. Ecological Modelling 210, 478–486 (2008). 99. Post, E. Ecology of climate change. The importance of biotic interactions. (Princeton University Press, 86. Wisz, M. S. & Guisan, A. Do pseudo-absence 2013). selection strategies influence species distribution models and their predictions? An information- 100. Araújo, M. B. & Pearson, R. G. Equilibrium of theoretic approach based on simulated data. BMC species’ distributions with climate. Ecography 28, Ecology 9, 8 (2009). 693–695 (2005). 87. Jiménez-Valverde, A. & Lobo, J. M. Threshold 101. Nogués-Bravo, D. Predicting the past distribution criteria for conversion of probability of species of species climatic niches. Global Ecology and presence to either–or presence–absence. Acta Biogeography 18, 521–531 (2009). Oecologica 31, 361–369 (2007). 102. Araújo, M. B. & Peterson, A. T. Uses and misuses of 88. Nenzén, H. K. & Araújo, M. B. Choice of threshold bioclimatic envelope modeling. Ecology 93, 1527– alters projections of species range shifts under 1539 (2012). climate change. Ecological Modelling 222, 3346– 103. Peterson, A. T. Uses and requirements of 3354 (2011). ecological niche models and related distribution 89. Williams, J. W. & Jackson, S. T. Novel climates, no- models. Biodiversity informatics 3, 59–72 (2006). analog communities, and ecological surprises. 104. Garcia, R. A., Burgess, N. D., Cabeza, M., Rahbek, C. Frontiers in Ecology and the Environment 5, 475– & Araújo, M. B. Exploring consensus in 21st 482 (2007). century projections of climatically suitable areas
22 | Synopsis
for African vertebrates. Global Change Biology 18, 120. Real, R., Luz Márquez, A., Olivero, J. & Estrada, A. 1253–1269 (2012). Species distribution models in climate change scenarios are still not useful for informing policy 105. Garcia, R. A. et al. Matching species traits to planning: an uncertainty assessment using fuzzy projected threats and opportunities from climate logic. Ecography 33, 304–314 (2010). change. Journal of Biogeography 10.1111/jbi.12257 (2014). 121. Bagchi, R. et al. Evaluating the effectiveness of conservation site networks under climate change: 106. McClean, C. J. et al. African plant diversity and accounting for uncertainty. Global Change Biology climate change. Annals of the Missouri Botanical 19, 1236–1248 (2013). Garden 92, 139–152 (2005). 122. Peterson, A. T. et al. Future projections for 107. Huntley, B. et al. Potential impacts of climatic Mexican faunas under global climate change change upon geographical distributions of birds. scenarios. Nature 416, 626–629 (2002). Ibis 148, 8–28 (2006). 123. Fordham, D. A. et al. Managed relocation as an 108. Thuiller, W. et al. Vulnerability of African adaptation strategy for mitigating climate change mammals to anthropogenic climate change under threats to the persistence of an endangered lizard. conservative land transformation assumptions. Global Change Biology 18, 2743–2755 (2012). Global Change Biology 12, 424–440 (2006). 124. Bates, J. M. & Granger, C. W. J. The Combination of 109. Barbet-Massin, M., Walther, B. A., Thuiller, W., Forecasts. Operational Research Quarterly 20, Rahbek, C. & Jiguet, F. Potential impacts of climate 451–468 (1969). change on the winter distribution of Afro- Palaearctic migrant passerines. Biology Letters 5, 125. Jun, M., Knutti, R. & Nychka, D. W. Spatial Analysis 248–251 (2009). to Quantify Numerical Model Bias and Dependence. Journal of the American Statistical 110. Hole, D. G. et al. Projected impacts of climate Association 103, 934–947 (2008). change on a continent-wide protected area network. Ecology Letters 12, 420–31 (2009). 126. Knutti, R., Furrer, R., Tebaldi, C., Cermak, J. & Meehl, G. A. Challenges in Combining Projections 111. New, M., Lister, D., Hulme, M. & Makin, I. A high- from Multiple Climate Models. Journal of Climate resolution data set of surface climate over global 23, 2739–2758 (2010). land areas. Climate Research 21, 1–25 (2002). 127. Beaumont, L. J., Hughes, L. & Pitman, A. J. Why is 112. Meehl, G. A. et al. in (Solomon S et al.) 747–845 the choice of future climate scenarios for species (Cambridge University Press, 2007). distribution modelling important? Ecology Letters 113. Araújo, M. B., Thuiller, W. & Pearson, R. G. Climate 11, 1135–1146 (2008). warming and the decline of amphibians and 128. Duan, Q. & Phillips, T. J. Bayesian estimation of reptiles in Europe. Journal of Biogeography 33, local signal and noise in multimodel simulations of 1712–1728 (2006). climate change. Journal of Geophysical Research 114. Lawler, J. J. et al. Projected climate-induced faunal 115, D18123 (2010). change in the Western Hemisphere. Ecology 90, 129. Naujokaitis-Lewis, I. R. et al. Uncertainties in 588–597 (2009). coupled species distribution–metapopulation 115. Buisson, L., Thuiller, W., Casajus, N., Lek, S. & dynamics models for risk assessments under Grenouillet, G. Uncertainty in ensemble climate change. Diversity and Distributions 19, forecasting of species distribution. Global Change 541–554 (2013). Biology 16, 1145–1157 (2010). 130. Williams, J. W., Jackson, S. T. & Kutzbach, J. E. 116. Diniz-Filho, J. A. F. et al. Partitioning and mapping Projected distributions of novel and disappearing uncertainties in ensembles of forecasts of species climates by 2100 AD. Proceedings of the National turnover under climate change. Ecography 32, Academy of Sciences of the United States of America 897–906 (2009). 104, 5738–42 (2007). 117. Diniz-Filho, J. A. F. et al. Ensemble forecasting 131. Araújo, M. B., Whittaker, R. J., Ladle, R. J. & Erhard, shifts in climatically suitable areas for Tropidacris M. Reducing uncertainty in projections of cristata (Orthoptera: Acridoidea: Romaleidae). extinction risk from climate change. Global Insect Conservation and Diversity 3, 213–221 Ecology and Biogeography 14, 529–538 (2005). (2010). 132. Marmion, M., Parviainen, M., Luoto, M., Heikkinen, 118. Dormann, C. F., Purschke, O., García Márquez, J. R., R. K. & Thuiller, W. Evaluation of consensus Lautenbach, S. & Schröder, B. Components of methods in predictive species distribution uncertainty in species distribution analysis: a case modelling. Diversity and Distributions 15, 59–69 study of the Great Grey Shrike. Ecology 89, 3371– (2009). 86 (2008). 133. Araújo, M. B. & Luoto, M. The importance of biotic 119. Thuiller, W. Patterns and uncertainties of species’ interactions for modelling species distributions range shifts under climate change. Global Change under climate change. Global Ecology and Biology 10, 2020–2027 (2004). Biogeography 16, 743–753 (2007).
Uncertainty in projected impacts of climate change | 23
134. Kissling, W. D., Field, R., Korntheuer, H., Heyder, U. 148. Murray, K. A., Rosauer, D., McCallum, H. & Skerratt, & Böhning-Gaese, K. Woody plants and the L. F. Integrating species traits with extrinsic prediction of climate-change impacts on bird threats: closing the gap between predicting and diversity. Philosophical Transactions of the Royal preventing species declines. Proceedings of the Society B: Biological Sciences 365 , 2035–2045 Royal Society of London Series B: Biological (2010). Sciences 278, 1515–23 (2011). 135. Midgley, G. F., Hughes, G. O., Thuiller, W. & Rebelo, 149. Fritz, S. A., Bininda-Emonds, O. R. P. & Purvis, A. A. G. Migration rate limitations on climate change- Geographical variation in predictors of induced range shifts in Cape Proteaceae. Diversity mammalian extinction risk: big is bad, but only in and Distributions 12, 555–562 (2006). the tropics. Ecology Letters 12, 538–549 (2009). 136. Engler, R. & Guisan, A. MigClim: Predicting plant 150. Rowland, E., Davison, J. & Graumlich, L. distribution and dispersal in a changing climate. Approaches to Evaluating Climate Change Impacts Diversity and Distributions 15, 590–601 (2009). on Species: A Guide to Initiating the Adaptation 137. Veloz, S. D. et al. No-analog climates and shifting Planning Process. Environmental Management 47, 322–337 (2011). realized niches during the late quaternary: implications for 21st-century predictions by 151. Buckley, L. B. & Kingsolver, J. G. Functional and species distribution models. Global Change Biology Phylogenetic Approaches to Forecasting Species’ 18, 1698–1713 (2012). Responses to Climate Change. Annual Review of 138. Heikkinen, R. K. et al. Assessing the vulnerability Ecology, Evolution, and Systematics 43, 205–226 of European butterflies to climate change using (2012). multiple criteria. Biodiversity and Conservation 19, 152. Foden, W. B. et al. in Wildlife in a Changing World: 695–723 (2009). an analysis of the 2008 IUCN Red List of 139. Triviño, M., Cabeza, M., Thuiller, W., Hickler, T. & Threatened Species. (Vié, J.-C., Hilton-Taylor, C. & Araújo, M. B. Risk assessment for Iberian birds Stuart, S. N.) 77–88 (IUCN, 2008). under global change. Biological Conservation 168, 153. Ohlemüller, R. et al. The coincidence of climatic 192–200 (2013). and species rarity: high risk to small-range species 140. Foden, W. B. et al. Identifying the World’s Most from climate change. Biology Letters 4, 568–572 Climate Change Vulnerable Species: A Systematic (2008). Trait-Based Assessment of all Birds, Amphibians 154. Sandel, B. et al. The influence of Late Quaternary and Corals. PLOS ONE 8, e65427 (2013). climate-change velocity on species endemism. 141. Thuiller, W., Lavorel, S., Araujo, M. B., Sykes, M. T. Science 334, 660–4 (2011). & Prentice, I. C. Climate change threats to plant 155. Fjeldså, J., Ehrlich, D., Lambin, E. & Prins, E. Are diversity in Europe. Proceedings of the National biodiversity `hotspots’ correlated with current Academy of Sciences of the United States of America ecoclimatic stability? A pilot study using the 102, 8245–8250 (2005). NOAA-AVHRR remote sensing data. Biodiversity and Conservation 6, 401–422 (1997). 142. Garcia, R. A. & Araújo, M. B. Planejamento para a Conservação em um Clima em Mudança. 156. Hole, D. G. et al. Projected impacts of climate Natureza&Conservação 8 , 78–80 (English version change on a continent-wide protected area available in Supplementary (2010). network. Ecology Letters 12, 420–431 (2009). 143. Araújo, M. B., Williams, P. H. & Fuller, R. J. 157. Araújo, M. B., Alagador, D., Cabeza, M., Nogués- Dynamics of extinction and the selection of nature Bravo, D. & Thuiller, W. Climate change threatens reserves. Proceedings of the Royal Society of European conservation areas. Ecology letters 14, London. Series B: Biological Sciences 269, 1971– 484–92 (2011). 1980 (2002). 158. Stattersfield, A. J., Crosby, M. J., Long, A. J. & Wege, 144. Hodgson, J. A. et al. Habitat re-creation strategies D. C. Endemic Bird Areas of the World: Priorities for promoting adaptation of species to climate for Biodiversity Conservation. Birdlife change. Conservation Letters 4, 289–297 (2011). Conservation Series no.7. 846 (BirdLife 145. Serra-Diaz, J. M. et al. Bioclimatic velocity: the pace International, 1998). of species exposure to climate change. Diversity 159. Myers, N., Mittermeier, R. A., Mittermeier, C. G., da and Distributions 20, 169–180 (2013). Fonseca, G. A. B. & Kent, J. Biodiversity hotspots 146. González-Suárez, M. & Revilla, E. Variability in life- for conservation priorities. Nature 403, 853–858 history and ecological traits is a buffer against (2000). extinction in mammals. Ecology Letters 16, 242– 160. Bini, L. M., Diniz-Filho, J. A. F., Rangel, T. F. L. V. B., 51 (2013). Bastos, R. P. & Pinto, M. P. Challenging Wallacean 147. Isaac, N. J. B. & Cowlishaw, G. How species respond and Linnean shortfalls: knowledge gradients and to multiple extinction threats. Proceedings of the conservation planning in a biodiversity hotspot. Diversity & Distributions 12, 475–482 (2006). Royal Society of London. Series B: Biological Sciences 271, 1135–1141 (2004).
24 | Synopsis
161. Mittermeier, R. A. et al. Hotspots Revisited: Earth’s 177. Allen, C. D. et al. A global overview of drought and Biologically Richest and Most Endangered heat-induced tree mortality reveals emerging Ecoregions. 390 (CEMEX, 2004). climate change risks for forests. Forest Ecology and Management 259, 660–684 (2010). 162. Olson, D. M. & Dinerstein, E. The Global 200: a representation approach to conserving the Earth’ 178. Chamaillé-Jammes, S., Massot, M., Aragón, P. & s most biologically valuable ecoregions. Clobert, J. Global warming and positive fitness Conservation Biology 12, 502–515 (1998). response in mountain populations of common 163. Tingley, M. W., Estes, L. D. & Wilcove, D. S. lizards Lacerta vivipara. Global Change Biology 12, Ecosystems: Climate change must not blow 392–402 (2006). conservation off course. Nature 500, 271–272 179. Tyler, N. J. C., Forchhammer, M. C. & Øritsland, N. (2013). A. Nonlinear effects of climate and density in the 164. Jiménez, M. A. et al. Extreme climatic events dynamics of a fluctuating population of reindeer. change the dynamics and invasibility of semi-arid Ecology 89, 1675–1686 (2008). annual plant communities. Ecology letters 14, 180. Lane, J. E., Kruuk, L. E. B., Charmantier, A., Murie, J. 1227–35 (2011). O. & Dobson, F. S. Delayed phenology and reduced 165. Iwamura, T., Wilson, K. A., Venter, O. & fitness associated with climate change in a wild Possingham, H. P. A Climatic Stability Approach to hibernator. Nature 489, 554–557 (2012). Prioritizing Global Conservation Investments. 181. Ordonez, A. & Williams, J. W. Climatic and biotic PLOS ONE 5, e15103 (2010). velocities for woody taxa distributions over the 166. Wiens, J. A., Seavy, N. E. & Jongsomjit, D. Protected last 16 000 years in eastern North America. areas in climate space: What will the future bring? Ecology Letters 16, 773–781 (2013). Biological Conservation 144, 2119–2125 (2011). 182. Willis, K. J. & MacDonald, G. M. Long-Term Ecological Records and Their Relevance to Climate 167. Loarie, S. R. et al. The velocity of climate change. Change Predictions for a Warmer World. Annual Nature 462, 1052–1055 (2009). Review of Ecology, Evolution, and Systematics 42, 168. Ohlemüller, R. Running Out of Climate Space. 267–287 (2011). Science 334, 613–614 (2011). 183. Davis, M. B. & Shaw, R. G. Range Shifts and 169. Svenning, J.-C., Fløjgaard, C., Marske, K. A., Nógues- Adaptive Responses to Quaternary Climate Bravo, D. & Normand, S. Applications of species Change. Science 292, 673–679 (2001). distribution modeling to paleobiology. Quaternary Science Reviews 30, 2930–2947 (2011). 184. Blois, J. L., McGuire, J. L. & Hadly, E. A. Small mammal diversity loss in response to late- 170. Burrows, M. T. et al. The Pace of Shifting Climate in Pleistocene climatic change. Nature 465, 771–774 Marine and Terrestrial Ecosystems. Science 334, (2010). 652–655 (2011). 185. Nogués-Bravo, D., Ohlemüller, R., Batra, P. & 171. Ohlemüller, R., Gritti, E. S., Sykes, M. T. & Thomas, Araújo, M. B. Climate predictors of late quaternary C. D. Towards European climate risk surfaces: the extinctions. Evolution 64, 2442–2449 (2010). extent and distribution of analogous and non- analogous climates 1931–2100. Global Ecology 186. Littell, J. S., McKenzie, D., Kerns, B. K., Cushman, S. and Biogeography 15, 395–405 (2006). & Shaw, C. G. Managing uncertainty in climate- driven ecological models to inform adaptation to 172. Ackerly, D. D. et al. The geography of climate climate change. Ecosphere 2, art102 (2011). change: implications for conservation 187. Fronzek, S., Carter, T. R., Räisänen, J., Ruokolainen, biogeography. Diversity and Distributions 16, 476– L. & Luoto, M. Applying probabilistic projections of 487 (2010). climate change with impact models: a case study 173. Jackson, S. T. & Overpeck, J. T. Responses of Plant for sub-arctic palsa mires in Fennoscandia. Populations and Communities to Environmental Climatic Change 99, 515–534 (2010). Changes of the Late Quaternary. Paleobiology 26, 194–220 (2000). 188. Wenger, S. J. et al. Probabilistic accounting of uncertainty in forecasts of species distributions 174. Peñuelas, J. et al. Evidence of current impact of under climate change. Global Change Biology 19, climate change on life: a walk from genes to the 3343–3354 (2013). biosphere. Global Change Biology 19, 2303–2338 189. Hijmans, R. J. & Graham, C. H. The ability of climate (2013). envelope models to predict the effect of climate 175. Foden, W. et al. A changing climate is eroding the change on species distributions. Global Change geographical range of the Namib Desert tree Aloe Biology 12, 2272–2281 (2006). through population declines and dispersal lags. Diversity and Distributions 13, 645–653 (2007). 190. Morin, X. & Thuiller, W. Comparing niche- and process-based models to reduce prediction 176. Sinervo, B. et al. Erosion of Lizard Diversity by uncertainty in species range shifts under climate Climate Change and Altered Thermal Niches. change. Ecology 90, 1301–1313 (2009). Science 328, 894–899 (2010).
Uncertainty in projected impacts of climate change | 25
191. Kearney, M. R., Wintle, B. A. & Porter, W. P. Assessments and Conservation Decisions. Correlative and mechanistic models of species Conservation Biology 27, 1147–1157 (2013). distribution provide congruent forecasts under 204. Meller, L. et al. Ensemble distribution models in climate change. Conservation Letters 3, 203–213 conservation prioritization: from consensus (2010). predictions to consensus reserve networks. 192. Tingley, M. W. & Beissinger, S. R. Detecting range Diversity & Distributions (2014). shifts from historical species occurrences: new doi:10.1111/ddi.12162 perspectives on old data. Trends in Ecology & 205. Hof, C., Levinsky, I., Araújo, M. B. & Rahbek, C. Evolution 24, 625–33 (2009). Rethinking species’ ability to cope with rapid 193. Bierman, S. M., Butler, A., Marion, G. & Kühn, I. climate change. Global Change Biology 17, 2987– Bayesian image restoration models for combining 2990 (2011). expert knowledge on recording activity with 206. Fjeldså, J. How broad-scale studies of patterns and species distribution data. Ecography 33, 451–460 processes can serve to guide conservation (2010). planning in Africa. Conservation biology : the 194. Raxworthy, C. J. et al. Predicting distributions of journal of the Society for Conservation Biology 21, known and unknown reptile species in 659–67 (2007). Madagascar. Nature 426, 837–41 (2003). 207. Kerr, J. T., Kharouba, H. M. & Currie, D. J. The 195. Dessai, S., Hulme, M., Lempert, R. & Pielke, R. Do Macroecological Contribution to Global Change We Need Better Predictions to Adapt to a Solutions. Science 316 , 1581–1584 (2007). Changing Climate? Eos, Transactions American Geophysical Union 90, 111–112 (2009). 196. Knutti, R. & Sedlacek, J. Robustness and uncertainties in the new CMIP5 climate model Other publications or manuscripts in projections. Nature Climate Change 3, 369–373 preparation (2013). 197. Stevens, B. & Bony, S. What Are Climate Models R.A. Garcia and M.B. Araújo (2010). Planejamento para Missing? Science 340 , 1053–1054 (2013). a conservação em um clima em mudança (Conservation planning under changing climates). 198. Zimmermann, N. E. et al. Climatic extremes Natureza & Conservação 8(1): 78-80. improve predictions of spatial patterns of tree species. Proceedings of the National Academy of J. Forrest, N. D. Burgess, R. A. Garcia, A. Carlson, S. Sciences of the United States of America 106, Freeman, C. Loucks, Z. Maritim, S. Olimb, N. 19723–19728 (2009). Olwero, A. Schrag, S. Anstey. An Assessment of the Spatial Vulnerability of Selected Vertebrates to 199. Forrest, J. et al. An Assessment of the Spatial Land Use and Climate Change in the Greater Vulnerability of Selected Vertebrates to Land Use Ruvuma Landscape, Eastern Africa. In and Climate Change in the Greater Ruvuma preparation. Landscape, Eastern Africa. In preparation. S. G. Willis, W. Foden, R.G. Smith, N.D. Burgess, D.J. 200. Bennie, J. et al. Range expansion through Baker, E. Belle, J. Carr, N. Doswald, R. A. Garcia, A. fragmented landscapes under a variable climate. Hartley, C. Hof, T. Newbold, C. Rahbek, P. Visconti, Ecology Letters 16, 921–929 (2013). B. Young , S.H.M. Butchart. A framework 201. Evans, M. R. et al. Predictive systems ecology. methodology for integrating climate change Proceedings of the Royal Society B: Biological vulnerability assessments from Species Sciences 280, 20131452 (2013). Distribution Models and Trait-based Approaches 202. Burgman, M. A., Lindenmayer, D. B. & Elith, J. for use in systematic conservation planning. In Managing Landscapes for Conservation under preparation. Uncertainty. Ecology 86, 2007–2017 (2005). 203. Snover, A. K. et al. Choosing and Using Climate- Change Scenarios for Ecological-Impact
Chapter I
Exploring consensus in 21st century projections of climatically suitable areas for African vertebrates
RAQUEL A. GARCIA, NEIL D. BURGESS, MAR CABEZA, CARSTEN RAHBEK, AND MIGUEL B. ARAÚJO Global Change Biology 18, 1253-1269 (2012)
Exploring consensus in 21st century projections of climatically suitable areas for African vertebrates
RAQUEL A. GARCIA1,2,3, NEIL D. BURGESS2,4, MAR CABEZA1,5, CARSTEN RAHBEK2, and MIGUEL B. ARAÚJO1,2,3
1 Department of Biodiversity and Evolutionary Biology, National Museum of Natural Sciences, CSIC, Madrid, Spain 2 Center for Macroecology, Evolution and Climate, Department of Biology, University of Copenhagen, Denmark 3 Rui Nabeiro Biodiversity Chair, University of Évora, CIBIO, Évora, Portugal 4 WWF US Conservation Science Program, Washington, DC, US 5 Metapopulation Research Group, Department of Biosciences, University of Helsinki, Finland
Global Change Biology 18, 1253-1269 (2012)
Abstract
Africa is predicted to be highly vulnerable to 21st rise agreements among projections from the bio- century climatic changes. Assessing the impacts of climatic envelope models we compare five con- these changes on Africa’s biodiversity is, however, sensus methodologies, which generally increase plagued by uncertainties, and markedly different or retain projection accuracy and provide con- results can be obtained from alternative biocli- sistent estimates of species turnover. Variability matic envelope models or future climate projec- from emissions scenarios increases towards late- tions. Using an ensemble forecasting framework, century and affects southern regions of high spe- we examine projections of future shifts in climatic cies turnover centred in arid Namibia. Two-fold suitability, and their methodological uncertainties, differences in median species turnover across the for over 2,500 species of mammals, birds, amphib- study area emerge among alternative climate ians and snakes in sub-Saharan Africa. To summa- projections and emissions scenarios. Our ensem- rise a priori the variability in the ensemble of 17 ble of projections underscores the potential bias general circulation models, we introduce a con- when using a single algorithm or climate projec- sensus methodology that combines co-varying tion for Africa, and provides a cautious first ap- models. Thus we quantify and map the relative proximation of the potential exposure of sub- contribution to uncertainty of seven bioclimatic Saharan African vertebrates to climatic changes. envelope models, three multi-model climate pro- The future use and further development of bio- jections and three emissions scenarios, and ex- climatic envelope modelling will hinge on the plore the resulting variability in species turnover interpretation of results in the light of methodo- estimates. We show that bioclimatic envelope logical as well as biological uncertainties. Here, we models contribute most to variability, particularly provide a framework to address methodological in projected novel climatic conditions over uncertainties and contextualise results. Sahelian and southern Saharan Africa. To summa-
30 | Chapter I
Introduction considerable time. In face of the increasing number of climate projections and statistical Assessments of the potential effects of 21st functions available, calls have been made to century climatic changes on biodiversity com- explicitly address the methodological uncer- monly rely on bioclimatic envelope models tainties of BEMs so as to quantify the confi- (BEMs). Using correlations between climate dence that can be placed in forecasts (Thuiller and known species occurrences, BEMs estimate et al., 2004a, Araújo et al., 2005b, Araújo et al., future shifts in suitable climate for species. 2006, Pearson et al., 2006, Diniz-Filho et al., Widespread use of BEMs has been accompanied 2009, Wiens et al., 2009, Buisson et al., 2010). by discussions of the biological (e.g. Pearson & The high level of certainty typically required for Dawson, 2003, Araújo & Pearson, 2005, Sinclair policy making can hardly be attained using et al., 2007) and methodological (e.g. Heikkinen correlative models – least of all under unknown et al., 2006, Beaumont et al., 2008) uncertain- future climates – leaving model users with the ties that surround the outputs. BEMs rely on option of exploring the uncertainty in projec- assumptions about the association between tions and weighing the risks associated with climate and species distributions, and their alternative actions (Wiens et al., 2009). One biological realism hinges on additional factors approach to address such uncertainties is to influencing species vulnerability to climatic incorporate several assumptions and explore changes, such as ecophysiological and micro- the resulting range of potential results habitat preferences, phenotypic plasticity, evo- ("ensemble forecasting", Araújo & New, 2007). lutionary rates, dispersal ability (Chevin et al., Here, we provide the most extensive investiga- 2010, Dawson et al., 2011, Hof et al., 2011), and tion to date of methodological uncertainties biotic interactions (Araújo & Luoto, 2007, Suttle associated with ensemble forecasts of climate et al., 2007). In turn, BEM results are sensitive change impacts on sub-Saharan African verte- to the data and statistical functions that are brate species. used to describe the associations between spe- The African continent is predicted to be cies and climate. Alternative algorithms differ one of the most vulnerable to 21st century cli- regarding the data used, variable selection, matic changes (Boko et al., 2007, Collier et al., model parameterisation (Guisan & 2008). Forecasts of warming above the global Zimmermann, 2000, Elith et al., 2006, average (Christensen et al., 2007) are projected Heikkinen et al., 2006), and techniques for to affect African biodiversity and people’s live- extrapolation to novel conditions (Thuiller et lihoods (Velarde et al., 2005, Boko et al., 2007, al., 2004b, Pearson et al., 2006, Elith & Graham, Biggs et al., 2008). Yet, in comparison to well- 2009). BEMs are also sensitive to the green- researched regions such as Europe or North house gas emissions scenarios and climate America, Africa has received limited attention models used to simulate future climates regarding the potential impacts of climate (Beaumont et al., 2008). change on biodiversity (Lovett et al., 2005, Research is ongoing to develop more bio- Felton et al., 2009). Attribution of shifts in spe- logically realistic models (Keith et al., 2008, cies distributions to climate change is difficult Anderson et al., 2009, Brook et al., 2009, in Africa (Chown et al., 2010, Midgley & Huntley et al., 2010, Kearney et al., 2010), but Thuiller, 2011) because changes in water avail- the breath of information and data required to ability – the main determinant of ecological appropriately parameterise them is large. Sim- responses (Hawkins et al., 2003) – are spatially pler approaches such as BEMs are thus likely to complex and difficult to document (MacKellar remain important tools for assessing potential et al., 2007). Increases in temperature, how- impacts of climate change on biodiversity for a
Ensemble forecasting for African vertebrates | 31
ever, have been associated with the observed projected to contract in the Kalahari region, and range extension of the Common Swift (Hockey to shift eastward (Thuiller et al., 2006). For & Midgley, 2009), and are likely to have more birds, forecasts revealed losses in southern and severe impacts for tropical species (Deutsch et eastern Africa for breeding birds (Huntley et al., al., 2008, Wright et al., 2009). 2006) and trans-Saharan migrant passerines Previous studies using BEMs at a continen- (Barbet-Massin et al., 2009), but relatively tal scale have projected substantial geographi- small changes for the former in equatorial and cal shifts in suitable climate for African plants, moist tropical forest habitats and even gains for birds and mammals by late-century (Table 1). the latter in the Sahel and Arabian Peninsula. Specifically, McClean and colleagues (2005) Hole and colleagues (2009) also projected predicted losses of suitable climate for plant higher ensemble turnover for breeding birds in species in the Guineo-Congolian forests of southern Africa, despite high persistence of western and central Africa, and gains in the suitable climate across the network of Impor- surrounding uplands as well as the highlands of tant Bird Areas as a whole. The results pre- Namibia and the South African Drakensberg. sented in these studies are, however, contin- Mammal species ranges around the equatorial gent on the specific BEMs and future climate zone in central Africa were projected to shift projections used. Only one study (Barbet- westward, with contractions in the Congo Ba- Massin et al., 2009) has fitted a range of differ- sin, whereas mammals in southern Africa were ent BEMs. Most of them also spanned a limited
Table 1 | Published continental- and sub-continental-scale studies using correlative models to assess the impacts of climate change on African biodiversity. Species data Extent Res Modelling Future scenarios Ref. olu- approach tion Distributional Sub- 1o . Box model, SGA . 1 GCM: HadCM3 McClean data for 5,197 Saharan and BGA . 1 SRES: B1 et al., plant species (≥2 Africa . 3 time periods: 2025, 2055 2005 records) and 2085
Distributional Sub- 1o . Locally . 3 GCMs: HadCM3, Huntley et data for bird Saharan weighted GFDL_R30, al. 2006 species breeding Africa regression ECHAM4/OPYC3 in Africa . 1 SRES: B2 . 1 time period: 2070-99
Extent of occur- Africa 10’ . GAM . 1 GCM: HadCM3 Thuiller et rence data for 277 . 2 SRES: A2 and B2 al., 2006 mammal species . 2 time periods: 2050 and 2080
Distributional Sub- 1o . CRS and GAM . 3 GCMs: HadCM3, ECHAM4, Hole et al. data for 1,608 Saharan GFDL-R30 2009 bird species (≥5 Africa . 1 SRES: B2a records) . 3 time periods: 2025, 2055 and 2085
Distributional Sub- 0.5o . GLM, GAM, CTA, . 5 GCMs: BCM2, ECHAM5, Barbet- data for 64 bird Saharan ANN, MDA, HadCM3, MIROHIC3_2-HI Massin et species (≥6 re- Africa MARS, GBM, RF, and MK3) al. 2009 cords) MaxEnt . 3 SRES: A1B, B1 and A2 . Consensus . 3 time periods: 2030, 2065 projection and 2100
32 | Chapter I
number of General Circulation Models (GCMs) (Ogawa-Onishi et al., 2010), Australia and emissions scenarios, overlooking the vari- (Crossman & Bass, 2008), and Africa (Barbet- ability among simulations of future climates Massin et al., 2009, Coetzee et al., 2009), and which has been shown to be region-specific and have generally yielded higher accuracy than relatively high for most of Africa south of the single-models. Yet, there is still debate on the Equator (Hawkins & Sutton, 2009, 2011). In- best methodologies for combining BEM projec- vestigation of the level of uncertainty associ- tions and only a few comparisons have been ated with the results was, thus, limited. published (Araújo et al., 2005b, Araújo & New, In this article, we use seven BEM tech- 2007, Marmion et al., 2009). Here, we compare niques to describe the bioclimatic envelopes of five methodologies, including the methodology 284 amphibian, 310 snake, 623 mammal and introduced to combine GCMs. 1,506 bird species in sub-Saharan Africa. To Our projections provide insights into the assess the impacts of climatic changes on the potential exposure of sub-Saharan African ver- modelled species, we project their envelopes to tebrates to 21st century climatic changes. We mid- and late-century climates. We use an en- explicitly address the variability in forecasts of semble of 17 GCMs, under the B1, A1B, and A2 species temporal turnover from alternative emissions scenarios from the Intergovernmen- climate projections and BEMs, or combinations tal Panel on Climate Change (IPCC). As large thereof. More specifically, we investigate: a) the ensembles of projections are difficult to inter- relative contribution of different sources of pret, consensus methodologies can be used to uncertainty in forecasts of species turnover; b) average across ensembles. Multi-model climate the predictive accuracy of forecasts from alter- projections are widely used in climatology, but native BEM consensus methodologies; and c) there is still debate on the best consensus the variation in forecasts with alternative BEM methodologies to combine models (Tebaldi & consensus methodologies and climate projec- Knutti, 2007, Knutti et al., 2010). To retain tions. information about the full variability in our ensemble of GCMs, we introduce a methodology that averages co-varying GCMs based on their Data and Methods similarity in both magnitude and spatial pat- tern. We thus reduce the ensemble of GCMs to Species and climate data three multi-model projections, and obtain, for The study region covered continental sub- each of the 2,723 species modelled, 126 projec- Saharan Africa, south of 20°N. Species occur- tions (seven BEM techniques, three climate rence data for amphibian (Hansen et al., projections, three emissions scenarios, and two 2007a), snake (Rasmussen et al., 2007), mam- time periods), or 343,098 projections overall. mal (Galster et al., 2007) and terrestrial bird To facilitate interpretation of this ensemble, we species (Hansen et al., 2007b) in sub-Saharan also summarise agreements among projections Africa were used from the 1° resolution (≈111 from the seven BEMs. Consensus methodolo- km x 111 km at the Equator) databases held at gies have been used in previous climate change the Zoological Museum within the University of ecology work in Europe (Thuiller, 2004, Araújo Copenhagen in Denmark. This is the most com- et al., 2006, Buisson & Grenouillet, 2009, Araújo prehensive biodiversity dataset for Africa, et al., 2011, Thuiller et al., 2011), the Americas compiled from multiple sources and continu- (Diniz-Filho et al., 2009, Lawler et al., 2009, ously improved over 15 years (Burgess et al., Marini et al., 2009, Roura-Pascual et al., 2009, 1998, Brooks et al., 2001). Data were available Diniz-Filho et al., 2010, Marini et al., 2010), Asia for 741 amphibians, 477 snakes, 1,085 mam-
Ensemble forecasting for African vertebrates | 33
mals and 1,789 birds. Because of statistical bird species (Huntley 2006) and a variety of difficulties with modelling species with limited other species in Africa (Chown, 2010 and numbers of occurrence records (Stockwell & references therein). The temperature-based Peterson, 2002, Wisz et al., 2008), we excluded variables selected in our study further reflect species with fewer than 15 records over the the important effect of seasonal temperatures study area. Our threshold may introduce uncer- on species’ distributions (Huntley et al., 2006). tainty in the analysis, yet the effect of sample Future climate projections were derived size on model accuracy is a question that re- from 17 GCMs downscaled to 10’ resolution quires further study. More conservative (Tabor & Williams, 2010; see Appendix S2). All thresholds for sample size have been suggested downscaled GCMs were from the World Climate (e.g. Harrell et al., 1984), but some algorithms Research Programme’s Coupled Model Inter- have been shown to achieve 90% of their max- comparison Project phase 3 multi-model data- imum accuracy with ten records (Stockwell & set projections (Meehl et al., 2007) and de- Peterson, 2002). In total 2,723 species were biased using the change-factor technique and modelled, accounting for 67% of the available observational data from the CRU. The datasets data (284 amphibians, 310 snakes, 623 mam- comprise monthly mean temperature and pre- mals and 1,506 birds, with median range sizes cipitation projections for the 2041-60 and of 71, 90, 94 and 162 grid cells respectively). 2081-00 time intervals. Simulations from the Baseline climate data averaged for the 17 GCMs were used for three illustrative green- 1961-90 period were obtained from the Cli- house gas emissions scenarios from the IPCC’s matic Research Unit (New et al., 2002) at a Special Report Emissions Scenarios resolution of 10’ (≈18.6 km x 18.6 km at the (Nakicenovic et al., 2000). We used high- to Equator). Monthly precipitation and mean low-end scenarios (A2, A1B and B1) that reflect temperature values were used to compute 21 different assumptions about demographic, variables that are commonly useful in biocli- socio-economic and technological development matic modelling studies (see Appendix S1 in on greenhouse gas emissions. For each GCM Supporting Information). We applied principal projection and scenario we computed the three components analysis (PCA) to identify sets of selected variables over the study area. To uncorrelated variables that represent major match the species data resolution, both base- climatic gradients over the study area. From the line and future climate datasets were re- first axis we selected the variable with the sampled in ArcGIS (ESRI, 2006), using bilinear highest loading, i.e., annual precipitation. An interpolation, to the 1° grid over sub-Saharan additional variable was selected from the first Africa. Data processing and statistical analyses axis with high loading but opposite sign (corre- were performed in R (R Development Core lation -0.45): temperature of the warmest Team, 2010) version 2.11.1. month. The variable with the highest loading on the second axis, i.e., temperature of the coldest Combining ensembles of climate projections month, was the third variable selected. The two We first summarised the general tendencies first axes explained 74.3% of the variation among the 17 selected GCMs. In climatology, (Appendix S1). Together, precipitation and multi-model ensemble averages have often temperature influence water availability, which been shown to improve the outcome of climate controls biological activity in the tropics and simulations (Phillips & Gleckler, 2006, Tebaldi sub-tropics (Hawkins et al., 2003). Both pre- & Knutti, 2007, Gleckler et al., 2008, Reichler & cipitation- and temperature-based variables are Kim, 2008, Pierce et al., 2009, Knutti et al., important determinants of the distributions of 2010, but see Fordham et al., in press). How-
34 | Chapter I
ever, averaging ensembles can result in the loss cate high similarity between a given model and of higher-order variability reflected in extreme the multi-model median. projections (Beaumont et al., 2008). To avoid Second, the ensemble of 17 GCMs was par- this limitation, we used a hybrid consensus titioned into groups of co-varying projections approach (Araújo et al., 2006) – hereafter re- according to the 2081-00 D and R statistics ferred to as ‘central cluster’ – that groups co- obtained. We used k-means, a clustering tech- varying projections before averaging them. nique that assigns data points to the closest When there is great variation in projections – pre-defined centre. These centres were the as it is often the case with precipitation – this median points of single linkage hierarchical approach also avoids averaging projections that clusters based on the Euclidean distance matrix are very different or even contradictory, by (Venables & Ripley, 2003). The significance of placing them in different groups. the differences between clusters was tested To combine the GCMs, we used three steps with Anosim, a non-parametric test of analysis in R (R Development Core Team, 2010, see of similarity (Clarke & Warwick, 1994). Anosim Appendix S3 for the R scripts). First, as a basis was applied to the dissimilarity matrix of D and for identifying co-varying projections under R values to test whether the distances between each emissions scenario, we assessed similari- clusters were greater than the distances within ties among GCM simulations for each variable clusters. The initial number of clusters to ex- projected in the late-century, when inter- tract was selected so as to minimise inter- simulation spread becomes larger (Stott & cluster distances in the hierarchical trees, and Kettleborough, 2002, Meehl et al., 2007, was increased when needed until the Anosim Hawkins & Sutton, 2009). Similarities were test was statistically significant. For each emis- assessed separately for each variable because sions scenario we obtained a set of clusters, the performance of climate models varies for each with a number of co-varying late-century different variables (Lambert & Boer, 2001, climate models. The same clusters were applied Gleckler et al., 2008). We used model perform- to the baseline and 2041-60 time periods. ance metrics to characterise the agreement Third, we generated summaries of the pat- between individual simulations for each vari- terns of central tendency in each cluster ex- able and the multi-model median ensemble for tracted. For each cluster of co-varying GCMs the same variable. These metrics were spatially obtained for different emissions scenarios and aggregated point-wise measures of regional time slices, un-weighted median consensus deviations (Duan & Phillips, 2010). The spatial forecasts were computed on each variable. pattern Pearson correlation (R) reflects spatial There are contrasting views on the use of agreement between individual simulations of a weights to perform climate ensemble averages. given variable and the median simulation of They have been shown not to systematically that variable. The signed standardised anomaly change the results in some cases (e.g. Pierce et (D), in turn, measures signed differences in al., 2009) and to improve them in other cases magnitude between individual simulations and (e.g. Min & Hense, 2006). Which model per- the median simulation of a given variable, stan- formance metrics to use as weights also re- dardised using the standard deviation of all mains an issue of debate (Tebaldi & Knutti, simulations. D thus reflects whether a simula- 2007, Gleckler et al., 2008). Because the opti- tion tends to under- or over-estimate a given mal performance weights for future projections variable in relation to the median of all simula- are unlikely the same as for baseline climate tions of that variable, and by how much. D val- (Tebaldi & Knutti, 2007), we opted for un- ues close to zero and R values close to 1 indi- weighted averages. In summary, for each emis-
Ensemble forecasting for African vertebrates | 35
sions scenario and time period combination we (Chefaoui & Lobo, 2008). Yet, random selection obtained a set of clusters of GCMs, each with has been shown to result in higher predictive the median simulation of each variable com- power than strategies that select pseudo- puted across GCMs. For each set, the cluster absences from low-suitability regions (Wisz & with the average D closest to zero and the high- Guisan, 2009). est average R captured the maximum consen- For each species, the seven models were sus among projections, corresponding to the built using a calibration subset of 75% of the ‘central cluster’, whereas clusters departing sites selected at random and evaluated with the from the multi-model median ensemble cap- remaining 25% of the sites. This data-splitting tured extreme projections. procedure was repeated five times. Projections of the probability of occurrence of species in Bioclimatic envelope modelling each site were converted to binary format (presence/absence) using a threshold maximis- Models were fitted at 1° resolution using seven ing the True Skill Statistic (TSS, Allouche et al., presence-absence BEM techniques in BIOMOD 2006). The models were evaluated based on (Thuiller et al., 2009), a computing platform for median omission and commission errors and ensemble forecasting that operates in R envi- TSS on the cross-validated data. The calibrated ronment (R Development Core Team, 2010). models were used to generate projections of The techniques included three regression species’ bioclimatic envelopes under each GCM methods (generalised linear models (GLM), cluster and emissions scenario for 2041-60 and generalised additive models (GAM), and multi- 2081-00. The projections were based on the variate adaptive regression splines (MARS)), final runs of the models using 100% of the data, three machine-learning methods (artificial as data partitions have been shown to add sig- neural networks (ANN), Breiman and Cutler's nificant uncertainty to forecasts (Araújo et al., random forest for classification and regression 2005b, 2009). Because we were interested in (RF), and generalised boosting models (GBM)), measuring changes in climatic suitability for and one classification method (flexible dis- species rather than interpreting model projec- criminant analysis, FDA). Due to differences in tions as estimates of the changes in observed quality, species occurrence data were treated species distributions, we adopted an ‘unlimited differently across taxa. Estimated range maps dispersal’ scenario, whereby species are as- for mammals and birds, based on numerous sumed to be able to track shifting suitable cli- records of species across multiple countries, mate over the entire study area. To comple- lend themselves to be treated as presence- ment these projections, the areas where higher absence data. For most amphibians and snakes, proportions of species are projected to retain however, the data comprise confirmed speci- climatic suitability over time (corresponding to men locality records from museum collections a ‘no dispersal’ assumption) were also mapped. and thus were considered to be closer to pres- ence-only data. For the latter taxa, pseudo- Combining ensembles of BEM forecasts absences were randomly generated to allow fitting models that assume the data to be in the For each taxon, we explored the agreement form of presences and absences. The process of among projections from the seven BEM tech- generating random pseudo-absences in BIO- niques using five consensus methodologies. The MOD weighs them to achieve a prevalence of first three methodologies provide a synthetic 0.5. There is a debate on how to select pseudo- measurement of the central tendency in the absences, and the choice of selection method is frequency distribution of the projections ob- dependent on the purpose of the study tained from all BEMs. Implemented within the
36 | Chapter I
BIOMOD package, the ‘ensemble mean’ com- consensus forecasts in the baseline and 2041- putes the un-weighted mean (e.g. Buisson & 60 time periods. Grenouillet, 2009, Diniz-Filho et al., 2010), the The BEM consensus projections were built ‘ensemble weighted mean’ (e.g. Marmion et al., using 100% of the data for the same reasons 2009) uses the TSS values as weights, and the cited above for single-BEM projections. To ‘ensemble median’ calculates the second quar- evaluate the consensus projections we applied tile of the frequency distribution of forecasts the same five consensus methodologies to the from all models (e.g. Araújo et al., 2005b, five evaluation datasets (25% of the data) and Lawler et al., 2009). The fourth and fifth meth- computed the median omission error, commis- odologies, in turn, pre-select projections that sion error and TSS of all cross-validations and, best summarise consensus among them. in the case of amphibians and snakes, all With the fourth method – ‘central model’ – pseudo-absence runs. To assess the level of we selected the model summarising the highest consensus among BEMs, we performed PCAs amount of variation among projections. For for each species on the probabilistic projec- each species, PCAs were performed on the tions, both for all seven BEMs and for the ‘cen- projected probabilities within the BIOMOD tral cluster’ BEMs only. The proportion of vari- package, and the ‘central model’ corresponded ance explained by the first principal component to the one with the highest PCA loading in the axis provided a measure of consensus among first (consensus) axis (e.g. Thuiller et al., 2005, projections (Thuiller, 2004, Araújo et al., 2006, Algar et al., 2009). In the fifth method – ‘central Grenouillet et al., 2011). cluster’ – we investigated patterns of central tendency among groups of co-varying projec- Mapping shifts in climatic suitability and associated tions (e.g. Araújo et al., 2005b, Araújo et al., uncertainties 2006). Following the approach used to combine For each emissions scenario and GCM cluster the ensemble of GCMs, we clustered the BEMs combination, and for the five BEM consensus based on the similarities between single-BEM methodologies applied, baseline and future probabilistic projections and the multi-model species richness and the number of contracting median probabilistic projection. We used the and expanding species were computed. Spatial same measures of D and R, computed for each patterns of change were investigated using species. Unlike the other four consensus meth- measures of species temporal turnover per grid odologies, the grouping was performed for the cell (Peterson et al., 2002). The ‘species turn- set of all species in one taxon rather than indi- over rate’ refers to local dissimilarities between vidually for each species. The same procedure baseline and future sets of species for which a of clustering and Anosim testing we used for given area is projected to be climatically suit- GCM projections was followed. The correspond- able, and thus incorporates both losses and ing binary projections for the BEMs in each gains of climate space. In addition, in situ per- cluster were combined using a majority vote sistence of climatic suitability for species was criterion, whereby a species was considered investigated. present in grid points where more than half of To evaluate and map the relative contribu- the BEMs in the cluster predicted presence. The tions of emissions scenarios, future climates ‘central cluster’ was the closest to the multi- and BEMs to the overall uncertainty in fore- model projection. Both the PCA and cluster casts, we performed analyses of variance analysis were applied to the end-of-century (ANOVA) in R (R Development Core Team, scenario, when divergence is expected to be 2010). Following Diniz-Filho and colleagues highest, and the same model(s) used to derive (2009), we performed a three-way ANOVA
Ensemble forecasting for African vertebrates | 37
without replication for each grid cell, using the the general trend across models (Figure 1a; see turnover rate as the response variable and the below). emissions scenarios, future climate projections Disagreement among BEM forecasts in late- and BEM consensus methodologies as factors. century was mainly concentrated in the north- An ANOVA using single-BEMs, climates and ern half of the study area, particularly in Sahe- emissions scenarios as factors, in turn, pro- lian and southern Saharan Africa (Figure 2). vided indications on the relative contribution of These areas were predicted to experience mean individual BEMs to uncertainty in turnover temperatures of the warmest and coldest projections (before combining the ensembles). months above the calibration range (Figure 3a; To explore BEM uncertainties associated with Appendix S5), forcing the models to extrapolate model extrapolation, we mapped the areas beyond known relationships. Late-century non- experiencing future climates beyond the range analogue climates covered up to 50% of the of climate values used to fit the models for each study area for the most severe climate projec- of the three variables. For each GCM cluster and tion. A comparison of grid points with analogue scenario combination, the sum of these areas and non-analogue climates (Figure 3b) re- corresponded to ‘non-analogue’ areas, where vealed significant differences in the proportion projections become statistically less reliable of the total sum of squares attributed to BEMs (Heikkinen et al., 2006, Williams et al., 2007, for all time periods and future climates (Kol- Fitzpatrick & Hargrove, 2009). mogorov-Smirnov tests p-value<0.001). The northern regions of high BEM uncertainty mostly corresponded to areas with non- Results analogue climates.
Relative contribution of different sources of Predictive accuracy of forecasts from alternative BEM uncertainty in forecasts of species turnover consensus methodologies Uncertainty in species turnover forecasts was Using the median TSS across all species as an mainly caused by the variability among BEMs. evaluation criterion, we found that consensus In the point-wise ANOVA using BEMs, GCM projections outperformed or equalled all or six clusters, and emissions scenarios as factors, the of the single-BEMs for amphibians and snakes, median proportions of the total sum of squares and between one and seven single-BEMs for across the study area attributed to BEMs mammals and birds (Figure 4a). The ‘central reached between 76% for mammals and 82% cluster’ methodology provided the most robust for snakes by mid-century (Appendix S4). The projections for all taxa, surpassing all single- relative contribution of BEMs to overall uncer- BEMs. This methodology resulted in four clus- tainty decreased by late-century (to between ters for each taxon, with the ‘central cluster’ 61% for mammals and 69% for snakes) due to combining high-accuracy models (GAM and increased divergence among emissions scenar- GLM for mammals, GAM, GBM, GLM and RF for ios. Variability across BEMs was strongly af- birds, and GBM and GLM for amphibians and fected by RF projections that displayed higher snakes for most GCM cluster and scenario com- losses and gains of suitable climate across cli- binations; Anosim p-value<0.05, see Appendix mates and emissions scenarios, departing from S6 for results on 100% of the data). The ‘central cluster’ projections were also the most consen-
38 | Chapter I
Figure 1 | Proportion of suitable climate projected to be lost and gained by species for the seven biocli- matic envelope models (BEM) (a) and the five BEM consensus projections (b) under alternative climate projections. Values are median percentages of grid cells lost or gained for all species of amphibians (n=284), snakes (n=310), mammals (n=623) and birds (n=1,506) in mid- (open circles) and late-century (solid circles). For a given time period, each circle corresponds to one of the nine combinations of three emissions scenarios and three general circulation model (GCM) clusters. The BEMs are Artificial Neural Networks (ANN), Generalised Additive Models (GAM), Generalised Boosting Model (GBM), Generalised Linear Models (GLM), Multivariate Adaptive Regression Splines (MARS), Flexible Discriminant Analysis (FDA), and Random Forests (RF), and the BEM consensus are ensemble mean (EMean), ensemble weighted mean (EWMean), ensemble median (EMed), central model (CMod), and central cluster (CClus). sual, raising the median levels of consensus (Figure 4c; for each taxon, distribution of com- across all species for the four taxa, although to mission errors across species significantly dif- different degrees depending on the number and ferent from the remaining consensus projec- spread of projections in the cluster (Appendix tions according to Wilcoxon signed rank tests, S7). The ‘central model’ methodology, in turn, see Appendix S9). For birds in particular, the yielded high accuracy projections for amphibi- models combined in the ‘central cluster’ projec- ans and snakes but the lowest accuracy of all tions included RF, which incorporate the notion consensus projections for mammals and birds. of ensemble forecasting (Araújo & New, 2007) Whereas for a large fraction of amphibian and and can yield highly accurate projections snake species the ‘central model’ corresponded (Prasad et al., 2006). The extreme discrepancy to high-accuracy models (GBM and ANN), for in commission error of RF models from the mammal and bird species the selection covered other projections, however, suggests that these a wider range of models (see Appendix S8 for models may have over-fitted the training data results on 100% of the data) with varying levels (see Jimenez-Valverde et al., 2008). For am- of accuracy, yielding lower median TSS across phibians and snakes, the measurements of species. commission error may have been affected by High accuracy of consensus projections random pseudo-absences placed in climatically was linked to both low omission and commis- suitable areas (Anderson et al., 2003, Peterson sion errors (Figure 4b and c). The ‘central clus- et al., 2011). ter’ projections displayed lower numbers of known absence points incorrectly predicted
Ensemble forecasting for African vertebrates | 39
Figure 2 | Spatial variation of the relative contribution of future climate projections and bioclimatic envelope models to the variability in species turnover forecasts for amphibian, snake, mammal and bird species. Values shown correspond to the proportion of the total sum of squares accounted for by the bioclimatic envelope models (BEM), general circulation model clusters (GCMcons), emissions scenarios (SRES), and one interaction factor (GCMcons:SRES) in the three-way analysis of variance (ANOVA) per- formed for each grid cell over the study area (N=1,851) on turnover forecasts for each taxon. Data are shown for mid- and late-century.
Variation in forecasts with alternative BEM terns of species turnover for each taxon, with consensus methodologies and climate projections larger variations for amphibians and snakes between the ‘central cluster’ and the remaining When the BEMs were combined and the rela- consensus projections (Appendix S11). Consis- tive sources of uncertainty re-assessed in a tent across all taxa was a trend towards more point-wise ANOVA, differences emerged across species contracting their bioclimatic envelopes taxa. For amphibians and snakes, mid-century than expanding. Depending on the BEM con- turnover forecasts varied most with the BEM sensus projection, between 54% and 74% of consensus methodologies (median proportions the species in each taxon were consistently of the total sum of squares across the study estimated to lose suitable climate across emis- area of 37% and 43% respectively, decreasing sions scenarios by late-century (Figure 5a). to 24% and 29% by late-century; see Table 2). In contrast, alternative climate projections For these taxa, ‘central cluster’ median esti- had the largest impact on mid-century projec- mates of late-century species turnover were up tions for mammals and birds, with the GCM to 1.6 times higher than those produced by the clusters explaining 34% and 40% respectively most conservative BEM consensus methodolo- of the total sum of squares in the ANOVA (Table gies across all climate projections (Appendix 2). Towards the end of the century, the spread S10). By contrast, only up to 1.2-fold differ- across emissions scenarios increased, becoming ences resulted for mammal and bird species. the major source of uncertainty for all taxa. For For each taxon, however, there were consistent each scenario, the three clusters of GCMs ob- trends among consensus methodologies of tained (Anosim statistics 0.81 (A2), 0.72 (A1B) median losses and gains across climate projec- and 0.81 (B1), P=0.001) reflected a warming tions (Figure 1b). Geographically, the five BEM gradient. The maximum consensus ‘central consensus projections displayed similar pat-
40 | Chapter I
of the coldest month; see Appendix S12). Trends across the GCM clusters were less clear for precipitation forecasts, and sometimes showing contrasting directions of change (Ap- pendix S13), but cluster 3 was consistently the driest for all scenarios. Following the warming gradient, median late-century turnover rates across the study area almost doubled from cluster 1 under B1 to cluster 3 under A2 (Ap- pendix S10). For all taxa, a southern area centred in the arid regions of Namibia emerged with high turnover rates by late-century (Figure 6). How- ever, the geographical extent of this high- turnover effect varied from isolated patches of the Kalahari in Namibia and Botswana and of southern Mozambique for cluster 1 under B1, Figure 3 | Comparison of uncertainties arising to most of inland Namibia, Botswana and from bioclimatic envelope models between ana- southern Mozambique for cluster 3 under A2. A logue and non-analogue climate grid cells, for comparison of the areas projected to remain amphibians, snakes, mammals and birds. The climatically suitable for species over time re- maps (a) show the distribution over the study area (N=1,851) of non-analogue climates (dark vealed more striking differences across taxa grey) for scenarios A2 (N=738 non-analogue grid (Figure 7). For amphibians, it was the West cells), A1B (N=632) and B1 (N=374). The graphs African forests that were projected to remain (b) show the frequency of the distributions of the climatically suitable for the largest proportion proportions of the total sum of squares accounted for by bioclimatic envelope models in the point- of species, irrespective of climate projection. A wise three-way analysis of variance (ANOVA) significant proportion of snake species were performed for each taxon using species turnover also projected to persist in this area, but under projections as response variable. The difference is the B1 scenario also in the Albertine Rift moun- shown between analogue (light grey) and non- analogue (dark grey) climate grid cells over the tain forests and extending around the Congo study area (Kolmogorov-Smirnov tests p- Basin. In comparison, larger climatically stable value<0.001 across taxa). Data refer to the ‘maxi- areas were projected for mammal and bird mum consensus’ general circulation model cluster species, particularly under B1, with the highest (cluster 2) under the three emissions scenarios. proportions of species persisting in the Ethio- pian highland mountain, Albertine Rift, East cluster’ (cluster 2) reflected intermediate levels Africa montane and Eastern Arc forests, as well of warming. The other clusters captured low- as the Angolan scarp and Miombo woodlands (cluster 1) and high-end (cluster 3) tempera- and the Drakensberg and eastern coast of South ture variability across climate models. For the Africa for bird species. Projections for cluster 1 temperature of the warmest month, for exam- under B1 showed the lowest proportions of ple, median values of late-century anomalies species of all taxa consistently predicted to projected over the study area varied between contract across all BEM consensus projections 1.8°C (lower and upper quartiles 1.6-2.0%) for (Figure 5b). the low-end cluster 1 under B1 and 5.1°C (4.6- 5.7°C) for the high-end cluster 3 under A2 (similar patterns emerged for the temperature
Ensemble forecasting for African vertebrates | 41
Figure 4 | Comparison between average True Skill Statistics (TSS) of bioclimatic envelope models (BEM) and BEM consensus projections (a), and density functions of omission error (b) and commission error (c) of BEM consensus projections for all species in each taxon. Omission and commission error and TSS values refer to the median of the validation datasets (and pseudo-absence runs in the case of am- phibians and snakes). The TSS plots (a) show the median (full circles) and the upper and lower quartiles (the extremes of the horizontal lines) of the TSS values of all species for each BEM consensus projection, as well as the median TSS values of all species for each BEM projection (black vertical lines). The consen- sus projections are the ensemble mean (EMean), ensemble weighted mean (EWMean), ensemble median (EMed), central model (CMod) and central cluster (CClus). The CMod and CClus projections shown were built using late-century projections from the GCM cluster 2 under A1B. The individual models are Artifi- cial Neural Networks (ANN), Generalised Additive Models (GAM), Generalised Boosting Model (GBM), Generalised Linear Models (GLM), Multivariate Adaptive Regression Splines (MARS), Flexible Discrimi- nant Analysis (FDA), and Random Forests (RF).
Discussion data, and the thresholds used to convert prob- abilistic to binary projections (Araújo et al., We explored the spread in estimates of climati- 2005b, Nenzén & Araújo, 2011). Yet, the uncer- cally suitable areas for African vertebrates tainty arising from future climate projections using alternative climate projections and emis- and from BEMs, which differ in how they ex- sions scenarios, as well as BEMs and combina- trapolate to non-analogue climates, has particu- tions thereof. Our aim was to address the lar relevance in the climate change context. methodological uncertainty associated with the results. Other methodological factors that are Sources of uncertainty not accounted for in our analysis are likely to Although comparisons of sources of uncertainty add further uncertainty to the projections. They in forecasts depend on the amount of variability include gaps and biases in species occurrence captured by each source (Diniz-Filho et al., data (Hortal et al., 2008), the selection of pre- 2009), our assessment spanned a wide variabil- dictors (Synes & Osborne, 2011), the temporal ity for the three sources considered: seven (Roubicek et al., 2010) and spatial resolution classes of BEM that perform differently (Kriticos & Leriche, 2010) of baseline climate
42 | Chapter I
Table 2 | Relative contribution of alternative general circulation models and bioclimatic envelope model consensus methodologies to overall uncertainty in species turnover projections for amphibian, snake, mammal and bird species. The values are proportions of the total sum of squares from the three-way analysis of variance (ANOVA) performed for each grid cell over the study area (N=1,851) to evaluate the relative contributions of the bioclimatic envelope model consensus methodologies (BEMcons), the gen- eral circulation model clusters (GCMcons), the emissions scenarios (SRES) and interactions among these factors, to the variability in turnover forecasts for each taxon. Values correspond to the median and (in brackets) the lower and upper quartiles of the proportions of the total sum of squares attributed to each factor over the study area, and are shown for both mid- and late-century. Amphibians Snakes Mammals Birds 2041-60 BEMcons 36.5 (20.4-60.5) 42.8 (24.8-60.6) 24.7 (11.8-44.1) 20.6 (11.4-34.9) GCMcons 24.4 (9.2-40.7) 25.7 (11.6-41.3) 34.1 (19.0-50.4) 40.2 (24.5-54.7) SRES 9.8 (4.4-16.3) 10.7 (6.0-16.7) 15.2 (9.5-21.7) 18.2 (11.7-25.7) BEMcons:GCMcons 3.4 (1.7-6.4) 3.1 (1.8-5.2) 3.2 (1.8-5.9) 2.0 (1.1-3.3) BEMcons:SRES 1.9 (0.9-3.8) 1.9 (1.1-3.4) 2.3 (1.3-4.1) 2.0 (0.9-4.7) GCMcons:SRES 4.2 (2-7.6.0) 3.9 (2.1-6.8) 4.3 (2.2-7.5) 4.5 (2.4-8.1) BEMcons:GCMcons:SRES 4.3 (2.2-8.5) 3.7 (2.1-5.9) 4.3 (2.5-7.5) 2.4 (1.3-4.4) 2081-00 BEMcons 23.6 (7.6-48.8) 28.8 (11.3-47) 12.7 (3.6-29.3) 12.4 (3.3-24.7) GCMcons 19.2 (8.9-27.7) 19.8 (11.7-27.2) 24.7 (17.6-31.1) 28.0 (21.3-34.6) SRES 28.5 (13.3-48.7) 32.1 (19.6-48.8) 44.2 (31.4-57.1) 45.8 (32.3-57.3) BEMcons:GCMcons 1.9 (0.9-3.7) 1.5 (0.7-3.2) 1.4 (0.6-2.8) 1.0 (0.4-1.8) BEMcons:SRES 1.8 (0.9-3.7) 1.7 (0.9-3.3) 1.9 (0.8-3.7) 1.4 (0.6-2.7) GCMcons:SRES 4.3 (2.2-8.2) 3.8 (2.0-6.9) 3.6 (2.0-6.1) 4.0 (2.0-6.8) BEMcons:GCMcons:SRES 3.4 (1.5-8.5) 2.6 (1.5-5.9) 2.4 (1.1-4.7) 1.2 (0.7-2.3)
Figure 5 | Percentage of species predicted to lose or gain suitable climate consistently across the three emissions scenarios (a) and across the five BEM consensus methodologies (b). The percentages of spe- cies of amphibians (n=284), snakes (n=310), mammals (n=623) and birds (n=1,506) that are consistent- ly projected to contract (lighter tones, left side of graphs) or expand (darker tones, right side of graphs) their bioclimatic envelopes for the ‘maximum consensus’ general circulation model cluster across all emissions scenarios (A2, A1B and B1) are shown for each consensus methodology (a). The percentages of species in each taxon that are consistently projected to contract or expand across all BEM consensus projections (EMean=ensemble mean, EWMean=ensemble weighted mean, EMed=ensemble median, CMod=central model and CClus=central cluster) are shown for each climate projection (b). Data are shown for mid- and late-century.
Ensemble forecasting for African vertebrates | 43
Figure 6 | Projected late-century species turnover rate (%) for amphibian, snake, mammal and bird species under the alternative emissions scenarios (A2, A1B and B1) and general circulation model clus- ters (Clusters 1, 2 and 3). Data refer to the ensemble median of all bioclimatic envelope model projec- tions.
when projected to the future (Thuiller et al., 2004a, Araújo et al., 2006, Pearson et al., 2006), three emissions scenarios that span a good portion of the range of the six IPCC scenarios (Manning et al., 2010), and 17 of the 23 IPCC fourth Assessment Report GCMs. In line with previous studies (Araújo et al., 2005b, Pearson et al., 2006, Diniz-Filho et al., 2009, 2010, Nenzén & Araújo, 2011), mid-century projec- tions were most affected by the choice of BEM (Figure 2). The relative amount of variation introduced by BEMs, however, decreases over time (Buisson et al., 2010 and this study), as divergence among emissions scenarios in- creases (Stott & Kettleborough, 2002, Hawkins & Sutton, 2009). By late-century, two-fold dif- Figure 7 | Percentage of species predicted to ferences emerged in species turnover projec- retain climatic suitability under each emissions tions with alternative GCM clusters and emis- scenario. The proportion of the total numbers of species of amphibians (n=284), snakes (n=310), sions scenarios. mammals (n=623) and birds (n=1,506) that are Differences among projections from alter- projected to retain climatic suitability in each native climate projections were especially im- location are shown for the median ensemble of all portant in southern Africa (Figure 2), where bioclimatic envelope models and for the ‘maxi- mum consensus’ general circulation model cluster some of the greatest changes in species turn- under the A2, A1B and B1 emissions scenarios. over were projected (Figure 6). In contrast, variability from BEMs was especially high in the northern half of the study area (Figure 2),
44 | Chapter I
coinciding with projected non-analogue cli- likely to be higher (Platts et al., 2008, mates (Figure 3). The problem of model ex- Fitzpatrick & Hargrove, 2009, Elith et al., 2010, trapolation into 21st century novel climates has Araújo et al., 2011, Roberts & Hamann, 2011) been overlooked in most continental-scale and using a range of BEMs that make different studies in Africa and elsewhere, although there assumptions about the responses of species in have been efforts to quantify it at different those areas (Pearson et al., 2006). geographical scales, for example in Europe (Saetersdal et al., 1998, Thuiller et al., 2004b, Role of consensus approaches Araújo et al., 2011), North America (Roberts & Although multi-model averages of GCMs have Hamann, 2011), and Australia (Fitzpatrick & rarely been used by ecologists (but see Hargrove, 2009, Elith et al., 2010). Projected Beaumont et al., 2011, Roberts & Hamann, novel climate conditions are unevenly distrib- 2011), their value for global change studies has uted worldwide (Williams et al., 2007). In Af- been recognised (Beaumont et al., 2008, rica, forecasts of severe warming increase the Fordham et al., in press). Our ensemble of GCMs risk of non-analogue climates, particularly reflects availability, and may thus not sample under high-end emissions scenarios (Williams the full range of uncertainty or guarantee inde- et al., 2007; Appendix S5). pendence of models (Meehl et al., 2007, Tebaldi Solutions to the problem of extrapolation & Knutti, 2007, Knutti et al., 2008), but the into non-analogue climates include classifying large number of models included is expected to such areas as species absences (Austin & minimise potential biases by model choice Meyers, 1996, Thuiller et al., 2004b), excluding (Knutti et al., 2010). Our approach to combine them from the analysis (Saetersdal et al., 1998), groups of co-varying GCMs enabled us to retain or using a larger calibration area before pro- information about the full variability of projec- jecting to the region of interest (Pearson et al., tions, including those that are extreme, and 2006). Yet, such solutions may lead to mislead- minimise the effect of combining often diver- ing results when the bioclimatic envelope of the gent precipitation projections. In face of the species has not been fully captured by the cali- wealth of climate simulations available, this bration data or when the species is not in equi- approach might prove increasingly useful in librium with climate (Svenning & Skov, 2004, ecological modelling studies. Araújo & Pearson, 2005). With long time hori- How alternative methodologies to combine zons, novel climates are expected to become BEMs affect forecasts of species bioclimatic more widespread (Williams et al., 2007), and envelopes has rarely been investigated (but see the problem of extrapolation may persist even Araújo et al., 2005b, Araújo et al., 2006, when using larger calibration areas, and all the Marmion et al., 2009). Measurements of model more when the species response curves are accuracy in the baseline context do not neces- high or increasing where truncated (Anderson sarily provide an indication of the models’ abil- & Raza, 2010, Webber et al., 2011). The ideal ity to transfer into future conditions (Araújo et solution of using ecological or physiological al., 2005a). However, in the absence of inde- knowledge of the species to classify non- pendent data for evaluation of the models, they analogue climate areas as either presences or can be used as benchmark for verification of the absences (Elith et al., 2010) is not feasible for consistency of alternative consensus method- most species where these data are lacking. The ologies. Using a variety of methodologies, we approach used in our study may thus reflect the found that consensus projections displayed current best practice: mapping non-analogue greater consistency in accuracy for amphibians climates to identify where uncertainties are and snakes. Whereas the three consensus
Ensemble forecasting for African vertebrates | 45
methodologies averaging the full ensemble of tat structure further limit the response of ver- BEM projections were generally consistent for tebrates. African habitats are controlled by the all taxa, for mammals and birds they diverged interaction among climate, atmospheric CO2 from those methodologies combining only pre- and disturbances like fire (Scheiter & Higgins, selected models (Figure 4a). This distinction 2009). In the case of grass-dominated savannas, across taxa resulted from the greater spread of the low levels of CO2 that triggered their devel- projection accuracy from the seven BEMs for opment during the last glacial are clearly being mammals and birds. For all taxa, however, the surpassed, with expected changes in tree cover ‘central cluster’ methodology stood out as the (Bond et al., 2003, Kgope et al., 2010) and asso- most accurate projection. In the case of am- ciated fauna (Sirami et al., 2009). The use in phibians and snakes, accurate and constrained BEMs of vegetation predictors derived from baseline ‘central cluster’ projections led to mechanistically-based dynamic vegetation higher species turnover rates than the remain- models (Triviño et al., In press) could capture ing consensus methodologies (Appendix S11). these effects and, to some extent, reduce uncer- tainties arising in non-analogue situations. Interpreting BEM outputs Exposure of species to climatic changes was measured in our study by three climatic The climatically suitable areas identified in this variables. The realism of our models thus de- study are areas where the exposure of species pends on the relationship between species to climatic changes can potentially allow the distributions and these variables. This relation- persistence of vertebrates through time. As- ship is unlikely to be constant through time, not sessing the potential response of species to only for its complexity but also for the changing these changes would necessitate mechanistic correlations among climatic variables (Morin & information about their sensitivity and adap- Lechowicz, 2008). Decoupling of patterns of tive capacity. Correlative studies have used covariance between predictor and proximal such information as a source of data to infer variables may undermine the value of BEMs absence (Elith et al., 2010) or presence points when extrapolating in time (Jackson et al., (Hijmans & Graham, 2006), as model predictors 2009, Elith et al., 2010). Our variables also (Rödder et al., 2009), or as complementary discount climatic variability nested across tem- information (Morin & Thuiller, 2009), but simi- poral and spatial scales. Multi-decadal climatic lar applications to large datasets like ours are means fail to capture fluctuations and rapid limited by data availability. Information on transitions between climate states (Jackson et species’ dispersal capacity, for example, would al., 2009). Projected changes in mean precipita- determine whether new suitable climate space tion in Southern Africa, for example, differ be- is accessible to species. Given the known varia- tween seasons and do not always parallel those tion in dispersal capacity among the four taxa, in extreme precipitation (Shongwe et al., 2009). effective range shifts would diverge more However, it is the interplay between temporal across taxa than the projected gains of suitable variability and species survival thresholds that climate (Figure 1). Mapping areas that remain determines the effects on species (Jackson et al., climatically suitable for most species through 2009). By the same token, the coarse spatial time reflects potential persistence of species in resolution used here overlooks microclimates situ (Figure 7), whereas integrating dispersal provided by topographical or vegetation fea- rate estimates in BEMs (Midgley et al., 2006) or tures (e.g. Shoo et al., 2011). combining BEM and migration model projec- BEM outputs need to be interpreted in the tions (Iverson et al., 2004) would indicate po- light of the methodological uncertainty ex- tential range shifts. Biotic factors such as habi-
46 | Chapter I
plored and the biological limitations discussed roecology, Evolution and Climate, and NDB also above. If carefully implemented, BEMs can, thanks WWF US. therefore, provide a first-order, parsimonious assessment of the changes in the distribution of suitable climate for species (Thuiller et al., References 2008, Jackson et al., 2009, Chevin et al., 2010, Algar AC, Kharouba HM, Young ER, Kerr JT (2009) Huntley et al., 2010). A new generation of mod- Predicting the future of species diversity: els that couple correlative with mechanistic macroecological theory, climate change, and direct tests of alternative forecasting methods. Ecography, approaches (Keith et al., 2008, Anderson et al., 32, 22-33. 2009, Brook et al., 2009, Huntley et al., 2010) is Allouche O, Tsoar A, Kadmon R (2006) Assessing the required to allow predictions of population accuracy of species distribution models: prevalence, persistence that have direct relevance to local kappa and the true skill statistic (TSS). Journal of Applied Ecology, 43, 1223-1232. or regional conservation (Jackson et al., 2009, Anderson BJ, Akcakaya HR, Araújo MB, Fordham DA, Chevin et al., 2010). Yet our ensemble forecast- Martinez-Meyer E, Thuiller W, Brook BW (2009) ing implementation of BEMs provides general Dynamics of range margins for metapopulations insights into the potential exposure of sub- under climate change. Proceedings of the Royal Society B-Biological Sciences, 276, 1415-1420. Saharan African vertebrates to climate change Anderson RP, Lew D, Peterson AT (2003) Evaluating at the continental scale, as well as a critical predictive models of species' distributions: criteria background for those seeking to interpret these for selecting optimal models. Ecological Modelling, 162, 211-232. results and use them as the basis for decision- making at this large scale. Anderson RP, Raza A (2010) The effect of the extent of the study region on GIS models of species geographic distributions and estimates of niche evolution: preliminary tests with montane rodents (genus Nephelomys) in Venezuela. Journal of Biogeography, Acknowledgements 37, 1378-1393. We are grateful to Thomas Phillips for advice Araújo MB, Pearson RG (2005) Equilibrium of species’ distributions with climate. Ecography, 28, 693-695. regarding the comparison of General Circula- Araújo MB, Pearson RG, Thuiller W, Erhard M (2005a) tion Models, and to Diogo Alagador for assis- Validation of species–climate impact models under tance in statistical analyses and R coding. climate change. Global Change Biology, 11, 1504- 1513. Thanks to Hedvig Nenzén and Sara Varela for Araújo MB, Whittaker RJ, Ladle RJ, Erhard M (2005b) commenting on the manuscript, and to three Reducing uncertainty in projections of extinction risk anonymous reviewers for constructive criti- from climate change. Global Ecology and cism. RAG is funded through a PhD studentship Biogeography, 14, 529-538. from the Portuguese Foundation for Science Araújo MB, Thuiller W, Pearson RG (2006) Climate warming and the decline of amphibians and reptiles and Technology with co-funding from the in Europe. Journal of Biogeography, 33, 1712-1728. European Social Fund Araújo MB, Luoto M (2007) The importance of biotic (SFRH/BD/65615/2009). Research by MBA is interactions for modelling species distributions funded through the FCT PTDC/AAC- under climate change. Global Ecology and Biogeography, 16, 743-753. AMB/98163/2008 project. RAG, MC and MBA Araújo MB, New M (2007) Ensemble forecasting of acknowledge the Spanish Research Council species distributions. Trends in Ecology & Evolution, (CSIC) for support, and RAG and MBA also 22, 42-47. thank the ‘Rui Nabeiro/Delta’ Biodiversity Araújo MB, Thuiller W, Yoccoz NG (2009) Reopening the Chair. MBA and MC also received support climate envelope reveals macroscale associations with climate in European birds. Proceedings of the through the RESPONSES project. NDB, CR, RAG National Academy of Sciences, 106, E45-E46. and MBA thank the Danish National Research Araújo MB, Alagador D, Cabeza M, Nogués-Bravo D, Foundation for supporting the Center of Mac- Thuiller W (2011) Climate change threatens
Ensemble forecasting for African vertebrates | 47
European conservation areas. Ecology Letters, 14, Chevin L-M, Lande R, Mace GM (2010) Adaptation, 484-492. Plasticity, and Extinction in a Changing Environment: Towards a Predictive Theory. PLoS Biol, 8, e1000357. Austin MP, Meyers JA (1996) Current approaches to modelling the environmental niche of eucalypts: Chown S, Hoffmann A, Kristensen T, Angilletta MJ, implication for management of forest biodiversity. Stenseth N, Pertoldi C (2010) Adapting to climate Forest Ecology and Management, 85, 95-106. change: a perspective from evolutionary physiology. Climate Research, 43, 3-15. Barbet-Massin M, Walther BA, Thuiller W, Rahbek C, Jiguet F (2009) Potential impacts of climate change Chown SL (2010) Temporal biodiversity change in on the winter distribution of Afro-Palaearctic transformed landscapes: a southern African migrant passerines. Biology Letters, 5, 248-251. perspective. Philosophical Transactions of the Royal Society B-Biological Sciences, 365, 3729-3742. Beaumont LJ, Hughes L, Pitman AJ (2008) Why is the choice of future climate scenarios for species Christensen JH, Hewitson B, Busuioc A et al. (2007) distribution modelling important? Ecology Letters, Regional Climate Projections. In: Climate Change 11, 1135-1146. 2007: The Physical Science Basis. Contribution of Working Group I to the Fourth Assessment Report of Beaumont LJ, Pitman A, Perkins S, Zimmermann NE, the Intergovernmental Panel on Climate Change (eds Yoccoz NG, Thuiller W (2011) Impacts of climate Solomon S, Qin D, Manning M, Chen Z, Marquis M, change on the world's most exceptional ecoregions. Averyt K B, Tignor M, Miller H L) pp 847-940. Proceedings of the National Academy of Sciences, 108, Cambridge University Press, Cambridge, United 2306-2311. Kingdom and New York, NY, USA. Biggs R, Simons H, Bakkenes M, Scholes RJ, Eickhout B, Clarke KR, Warwick RM (1994) Similarity-based testing van Vuuren D, Alkemade R (2008) Scenarios of for community pattern: the two-way layout with no biodiversity loss in southern Africa in the 21st replication. Marine Biology, 118, 167-176. century. Global Environmental Change-Human and Policy Dimensions, 18, 296-309. Coetzee BWT, Robertson MP, Erasmus BFN, Rensburg BJv, Thuiller W (2009) Ensemble models predict Boko M, Niang I, Nyong A et al. (2007) Africa. In: Climate Important Bird Areas in southern Africa will become Change 2007: Impacts, Adaptation and Vulnerability. less effective for conserving endemic birds under Contribution of Working Group II to the Fourth climate change. Global Ecology and Biogeography, 18, Assessment Report of the Intergovernmental Panel on 701-710. Climate Change (eds Parry M L, Canziani O F, Palutikof J P, Van Der Linden P J, Hanson C E), pp Collier P, Conway G, Venables T (2008) Climate change 433-467. Cambridge University Press, Cambridge UK. and Africa. Oxford Review of Economic Policy, 24, 337- 353. Bond WJ, Midgley GF, Woodward FI (2003) The importance of low atmospheric CO2 and fire in Crossman ND, Bass DA (2008) Application of common promoting the spread of grasslands and savannas. predictive habitat techniques for post-border weed Global Change Biology, 9, 973-982. risk management. Diversity and Distributions, 14, 213-224. Brook BW, Akçakaya HR, Keith DA, Mace GM, Pearson RG, Araújo MB (2009) Integrating bioclimate with Dawson TP, Jackson ST, House JI, Prentice IC, Mace GM population models to improve forecasts of species (2011) Beyond Predictions: Biodiversity extinctions under climate change. Biology Letters, 5, Conservation in a Changing Climate. Science, 332, 53- 723-725. 58. Brooks T, Balmford A, Burgess N et al. (2001) Toward a Deutsch CA, Tewksbury JJ, Huey RB, Sheldon KS, Blueprint for Conservation in Africa. Bioscience, 51, Ghalambor CK, Haak DC, Martin PR (2008) Impacts of 613-624. climate warming on terrestrial ectotherms across latitude. Proceedings of the National Academy of Buisson L, Grenouillet G (2009) Contrasted impacts of Sciences, 105, 6668-6672. climate change on stream fish assemblages along an environmental gradient. Diversity and Distributions, Diniz-Filho JAF, Mauricio Bini L, Fernando Rangel T, 15, 613-626. Loyola RD, Hof C, Nogués-Bravo D, Araújo MB (2009) Partitioning and mapping uncertainties in ensembles Buisson L, Thuiller W, Casajus N, Lek S, Grenouillet G of forecasts of species turnover under climate (2010) Uncertainty in ensemble forecasting of change. Ecography, 32, 897-906. species distribution. Global Change Biology, 16, 1145- 1157. Diniz-Filho JAF, Nabout JC, Bini LM, Loyola RD, Rangel TF, Nogues-Bravo D, Araújo MB (2010) Ensemble Burgess N, Fjeldså J, Rahbek C (1998) Mapping the forecasting shifts in climatically suitable areas for distributions of Afrotropical vertebrate groups. Tropidacris cristata (Orthoptera: Acridoidea: Species, 30, 16-17. Romaleidae). Insect Conservation and Diversity, 3, Chefaoui RM, Lobo JM (2008) Assessing the effects of 213-221. pseudo-absences on predictive distribution model Duan Q, Phillips TJ (2010) Bayesian estimation of local performance. Ecological Modelling, 210, 478-486. signal and noise in multimodel simulations of climate
48 | Chapter I
change. Journal of Geophysical Research, 115, Hawkins BA, Field R, Cornell HV et al. (2003) Energy, D18123. water, and broad-scale geographic patterns of species richness. Ecology 84(12), 3105–3117 Elith J, Graham CH, Anderson RP et al. (2006) Novel methods improve prediction of species' distributions Hawkins E, Sutton R (2009) The Potential to Narrow from occurrence data. Ecography, 29, 129-151. Uncertainty in Regional Climate Predictions. Bulletin of the American Meteorological Society, 90, 1095- Elith J, Graham CH (2009) Do they? How do they? WHY 1107. do they differ? On finding reasons for differing performances of species distribution models. Hawkins E, Sutton R (2011) The potential to narrow Ecography, 32, 66-77. uncertainty in projections of regional precipitation change. Climate Dynamics, 37, 407-418. Elith J, Kearney M, Phillips S (2010) The art of modelling range-shifting species. Methods in Ecology and Heikkinen RK, Luoto M, Araújo MB, Virkkala R, Thuiller Evolution, 1, 330-342. W, Sykes MT (2006) Methods and uncertainties in bioclimatic envelope modelling under climate ESRI (2006) ArcGIS. Environmental Systems Research change. Progress in Physical Geography, 30, 751-777. Institute (ESRI). Redlands CA. Hijmans RJ, Graham CH (2006) The ability of climate Felton A, Fischer J, Lindenmayer D et al. (2009) Climate envelope models to predict the effect of climate change, conservation and management: an change on species distributions. Global Change assessment of the peer-reviewed scientific journal Biology, 12, 2272-2281. literature. Biodiversity and Conservation, 18, 2243- 2253. Hockey PAR, Midgley GF (2009) Avian range changes and climate change: a cautionary tale from the Cape Fitzpatrick M, Hargrove W (2009) The projection of Peninsula. Ostrich, 80, 29-34. species distribution models and the problem of non- analog climate. Biodiversity and Conservation, 18, Hof C, Levinsky I, Araújo MB, Rahbek C (2011) 2255-2261. Rethinking species' ability to cope with rapid climate change. Global Change Biology, 17, 2987-2990. Fordham D, Wigley T, Brook B (in press) Multi-model climate projections for biodiversity risk assessments. Hole DG, Willis SG, Pain DJ et al. (2009) Projected Ecological Applications, doi:10.1890/1811- impacts of climate change on a continent-wide 0314.1891. protected area network. Ecology Letters, 12, 420-431. Galster S, Burgess ND, Fjeldså J, Hansen LA, Rahbek C Hortal J, Jiménez-Valverde A, Gómez JF, Lobo JM, Baselga (2007) One degree resolution databases of the A (2008) Historical bias in biodiversity inventories distribution of 1085 species of mammals in Sub- affects the observed environmental niche of the Saharan Africa. On-line data source-Version 1.00. species. Oikos, 117, 847-858. Zoological Museum, University of Copenhagen, Huntley B, Collingham YC, Green RE, Hilton GM, Rahbek Denmark. C, Willis SG (2006) Potential impacts of climatic Gleckler PJ, Taylor KE, Doutriaux C (2008) Performance change upon geographical distributions of birds. Ibis, metrics for climate models. Journal of Geophysical 148, 8-28. Research, 113, D06104. Huntley B, Barnard P, Altwegg R et al. (2010) Beyond Grenouillet G, Buisson L, Casajus N, Lek S (2011) bioclimatic envelopes: dynamic species' range and Ensemble modelling of species distribution: the abundance modelling in the context of climatic effects of geographical and environmental ranges. change. Ecography, 33, 621-626. Ecography, 34, 9-17. Iverson LR, Schwartz MW, Prasad AM (2004) Potential Guisan A, Zimmermann NE (2000) Predictive habitat colonization of newly available tree-species habitat distribution models in ecology. Ecological Modelling, under climate change: an analysis for five eastern US 135, 147-186. species. Landscape Ecology, 19, 787-799. Hansen AJ, Burgess ND, Fjeldså J, Rahbek C (2007a) One Jackson ST, Betancourt JL, Booth RK, Gray ST (2009) degree resolution databases of the distribution of 739 Ecology and the ratchet of events: Climate variability, species of amphibians in Sub-Saharan Africa. On-line niche dimensions, and species distributions. data source-version 1.00. Zoological Museum, Proceedings of the National Academy of Sciences, 106, University of Copenhagen, Denmark. 19685-19692. Hansen LA, Fjeldså J, Burgess ND, Rahbek C (2007b) One Jimenez-Valverde A, Lobo JM, Hortal J (2008) Not as degree resolution databases of the distribution of 1789 good as they seem: the importance of concepts in resident birds in Sub-Saharan Africa. On-line data species distribution modelling. Diversity and source-Version 1.00. Zoological Museum, University Distributions, 14, 885-890. of Copenhagen, Denmark. Kearney MR, Wintle BA, Porter WP (2010) Correlative Harrell FE, Lee KL, Califf RM, Pryor DB, Rosati RA (1984) and mechanistic models of species distribution Regression modelling strategies for improved provide congruent forecasts under climate change. prognostic prediction. Statistics in Medicine, 3, 143- Conservation Letters, 3, 203-213. 152.
Ensemble forecasting for African vertebrates | 49
Keith DA, Akçakaya HR, Thuiller W et al. (2008) Intergovernmental Panel on Climate Change (eds Predicting extinction risks under climate change: Solomon S, Qin D, Manning M, Chen Z, Marquis M, coupling stochastic population models with dynamic Averyt K B, Tignor M, Miller H L), pp 747-845. bioclimatic habitat models. Biology Letters, 4, 560- Cambridge University Press, Cambridge, United 563. Kingdom and New York, NY, USA. Kgope BS, Bond WJ, Midgley GF (2010) Growth Midgley GF, Hughes GO, Thuiller W, Rebelo AG (2006) responses of African savanna trees implicate Migration rate limitations on climate change-induced atmospheric [CO2] as a driver of past and current range shifts in Cape Proteaceae. Diversity and changes in savanna tree cover. Austral Ecology, 35, Distributions, 12, 555-562. 451-463. Midgley G, Thuiller W (2011) Potential responses of Knutti R, Allen MR, Friedlingstein P et al. (2008) A terrestrial biodiversity in Southern Africa to Review of Uncertainties in Global Temperature anthropogenic climate change. Regional Projections over the Twenty-First Century. Journal of Environmental Change, 11, 127-135. Climate, 21, 2651-2663. Min S-K, Hense A (2006) A Bayesian approach to climate Knutti R, Furrer R, Tebaldi C, Cermak J, Meehl GA (2010) model evaluation and multi-model averaging with an Challenges in Combining Projections from Multiple application to global mean surface temperatures Climate Models. Journal of Climate, 23, 2739-2758. from IPCC AR4 coupled climate models. Geophysical Research Letters, 33, 5pp. Kriticos DJ, Leriche A (2010) The effects of climate data precision on fitting and projecting species niche Morin X, Lechowicz MJ (2008) Contemporary models. Ecography, 33, 115-127. perspectives on the niche that can improve models of species range shifts under climate change. Biology Lambert SJ, Boer GJ (2001) CMIP1 evaluation and Letters, 4, 573-576. intercomparison of coupled climate models. Climate Dynamics, 17, 83-106. Morin X, Thuiller W (2009) Comparing niche- and process-based models to reduce prediction Lawler JJ, Shafer SL, White D, Kareiva P, Maurer EP, uncertainty in species range shifts under climate Blaustein AR, Bartlein PJ (2009) Projected climate- change. Ecology, 90, 1301-1313. induced faunal change in the Western Hemisphere. Ecology, 90, 588-597. Nakicenovic N, Alcamo J, Davis G et al. (2000) Special Report on Emissions Scenarios: A Special Report of Lovett JC, Midgely GF, Barnard P (2005) Climate change Working Group III of the Intergovernmental Panel on and ecology in Africa. African Journal of Ecology, 43, Climate Change, 599pp. Cambridge University Press, 167–169. Cambridge, U.K. MacKellar NC, Hewitson BC, Tadross MA (2007) Nenzén HK, Araújo MB (2011) Choice of threshold alters Namaqualand’s climate: recent historical changes projections of species range shifts under climate and future scenarios. Journal of Arid Environments, change. Ecological Modelling, 222, 3346– 3354. 70, 604-614. New M, Lister D, Hulme M, Makin I (2002) A high- Manning MR, Edmonds J, Emori S et al. (2010) resolution data set of surface climate over global land Misrepresentation of the IPCC CO2 emission areas. Climate Research, 21, 1-25. scenarios. Nature Geoscience, 3, 376-377. Ogawa-Onishi Y, Berry PM, Tanaka N (2010) Assessing Marini MÂ, Barbet-Massin M, Lopes LE, Jiguet F (2009) the potential impacts of climate change and their Major current and future gaps of Brazilian reserves conservation implications in Japan: A case study of to protect Neotropical savanna birds. Biological conifers. Biological Conservation, 143, 1728-1736. Conservation, 142, 3039-3050. Pearson RG, Dawson TP (2003) Predicting the impacts of Marini MÂ, Barbet-Massin M, Martinez J, Prestes NP, climate change on the distribution of species: are Jiguet F (2010) Applying ecological niche modelling bioclimate envelope models useful? Global Ecology to plan conservation actions for the Red-spectacled and Biogeography, 12, 361-371. Amazon (Amazona pretrei). Biological Conservation, 143, 102-112. Pearson RG, Thuiller W, Araújo MB et al. (2006) Model- based uncertainty in species range prediction. Marmion M, Parviainen M, Luoto M, Heikkinen RK, Journal of Biogeography, 33, 1704-1711. Thuiller W (2009) Evaluation of consensus methods in predictive species distribution modelling. Diversity Peterson AT, Ortega-Huerta MA, Bartley J, Sanchez- and Distributions, 15, 59-69. Cordero V, Soberón J, Buddemeier RH, Stockwell DRB (2002) Future projections for Mexican faunas under McClean CJ, Jon CL, Küper W et al. (2005) African Plant global climate change scenarios. Nature, 416, 626- Diversity and Climate Change. Annals of the Missouri 629. Botanical Garden, 92, 139-152. Peterson AT, Soberón J, Pearson RG, Anderson RP, Meehl GA, Stocker TF, Collins WD et al. (2007) Global Martínez-Meyer E, Nakamura M, Araújo MB (2011) Climate Projections. In: Climate Change 2007: The Ecological Niches and Geographic Distributions, Physical Science Basis. Contribution of Working Group 328pp. Monographs in Population Biology 49, I to the Fourth Assessment Report of the Princeton University Press.
50 | Chapter I
Phillips TJ, Gleckler PJ (2006) Evaluation of continental under Global Warming. Part I: Southern Africa. precipitation in 20th century climate simulations: Journal of Climate, 22, 3819-3837. The utility of multimodel statistics. Water Resources Shoo LP, Storlie C, Vanderwal J, Little J, Williams SE Research, 42, 10pp. (2011) Targeted protection and restoration to Pierce DW, Barnett TP, Santer BD, Gleckler PJ (2009) conserve tropical biodiversity in a warming world. Selecting global climate models for regional climate Global Change Biology, 17, 186-193. change studies. Proceedings of the National Academy Sinclair SJ, White MD, Newell GR (2007) How Useful Are of Sciences, 106, 8441-8446. Species Distribution Models for Managing Platts PJ, McClean CJ, Lovett JC, Marchant R (2008) Biodiversity under Future Climates? Ecology and Predicting tree distributions in an East African Society, 15, 8. biodiversity hotspot: model selection, data bias and Sirami C, Seymour C, Midgley G, Barnard P (2009) The envelope uncertainty. Ecological Modelling, 218, 121- impact of shrub encroachment on savanna bird 134. diversity from local to regional scale. Diversity and Prasad A, Iverson L, Liaw A (2006) Newer Classification Distributions, 15, 948-957. and Regression Tree Techniques: Bagging and Stockwell DRB, Peterson AT (2002) Effects of sample Random Forests for Ecological Prediction. size on accuracy of species distribution models. Ecosystems, 9, 181-199. Ecological Modelling, 148, 1-13. R Development Core Team (2010) R: A language and Stott PA, Kettleborough JA (2002) Origins and estimates environment for statistical computing. R Foundation of uncertainty in predictions of twenty-first century for Statistical Computing, Vienna, Austria. temperature rise. Nature, 416, 723-726. Rasmussen JB, Hansen LA, Burgess ND, Fjeldså J, Rahbek Suttle KB, Thomsen MA, Power ME (2007) Species C (2007) One degree resolution databases of the Interactions Reverse Grassland Responses to distribution of 467 species of snakes in Sub-Saharan Changing Climate. Science, 315, 640-642. Africa. On-line data source-Version 1.00. Zoological Museum, University of Copenhagen, Denmark. Svenning J-C, Skov F (2004) Limited filling of the potential range in European tree species. Ecology Reichler T, Kim J (2008) How Well Do Coupled Models Letters, 7, 565-573. Simulate Today's Climate? Bulletin of the American Meteorological Society, 89, 303-311. Synes NW, Osborne PE (2011) Choice of predictor variables as a source of uncertainty in continental- Roberts DR, Hamann A (2011) Predicting potential scale species distribution modelling under climate climate change impacts with bioclimate envelope change. Global Ecology and Biogeography, 20, 904- models: a palaeoecological perspective. Global 914. Ecology and Biogeography, DOI: 10.1111/j.1466- 8238.2011.00657.x. Tabor K, Williams JW (2010) Globally downscaled climate projections for assessing the conservation Rödder D, Schmidtlein S, Veith M, Lötters S (2009) Alien impacts of climate change. Ecological Applications, Invasive Slider Turtle in Unpredicted Habitat: A 20, 554-565. Matter of Niche Shift or of Predictors Studied? PLoS ONE, 4, e7843. Tebaldi C, Knutti R (2007) The use of the multi-model ensemble in probabilistic climate projections. Roubicek AJ, VanDerWal J, Beaumont LJ, Pitman AJ, Philosophical Transactions of the Royal Society A: Wilson P, Hughes L (2010) Does the choice of climate Mathematical, Physical and Engineering Sciences, 365, baseline matter in ecological niche modelling? 2053-2075. Ecological Modelling, 221, 2280-2286. Thuiller W (2004) Patterns and uncertainties of species' Roura-Pascual Nr, Brotons Ls, Peterson A, Thuiller W range shifts under climate change. Global Change (2009) Consensual predictions of potential Biology, 10, 2020-2027. distributional areas for invasive species: a case study of Argentine ants in the Iberian Peninsula. Biological Thuiller W, Araújo MB, Pearson RG, Whittaker RJ, Invasions, 11, 1017-1031. Brotons L, Lavorel S (2004a) Biodiversity conservation: Uncertainty in predictions of extinction Saetersdal M, Birks HJB, Peglar SM (1998) Predicting risk. Nature, 430, 145-148. changes in Fennoscandian vascular-plant species richness as a result of future climatic change. Journal Thuiller W, Brotons L, Araújo MB, Lavorel S (2004b) of Biogeography, 25, 111-122. Effects of restricting environmental range of data to project current and future species distributions. Scheiter S, Higgins SI (2009) Impacts of climate change Ecography, 27, 165-172. on the vegetation of Africa: an adaptive dynamic vegetation modelling approach. Global Change Thuiller W, Lavorel S, Araujo MB, Sykes MT, Prentice IC Biology, 15, 2224-2246. (2005) Climate change threats to plant diversity in Europe. Proceedings of the National Academy of Shongwe ME, van Oldenborgh GJ, van den Hurk BJJM, de Sciences of the United States of America, 102, 8245- Boer B, Coelho CAS, van Aalst MK (2009) Projected 8250. Changes in Mean and Extreme Precipitation in Africa
Ensemble forecasting for African vertebrates | 51
Thuiller W, Broennimann O, Hughes G, Alkemade JRM, performance of species distribution models. Diversity Midgley GF, Corsi F (2006) Vulnerability of African and Distributions, 14, 763-773. mammals to anthropogenic climate change under Wright SJ, Muller-Landau HC, Schipper JAN (2009) The conservative land transformation assumptions. Future of Tropical Species on a Warmer Planet. Global Change Biology, 12, 424-440. Conservation Biology, 23, 1418-1426. Thuiller W, Albert C, Araújo MB et al. (2008) Predicting global change impacts on plant species' distributions: Future challenges. Perspectives in Plant Ecology Evolution and Systematics, 9, 137-152. Supporting Information Thuiller W, Lafourcade B, Engler R, Araujo MB (2009) BIOMOD - a platform for ensemble forecasting of Appendix S1. Principal Components Analysis to species distributions. Ecography, 32, 369-373. select predictor variables. Thuiller W, Lavergne S, Roquet C, Boulangeat I, Araújo Appendix S2. General Circulation Models used MB (2011) Consequences of climate change to the Tree of Life in Europe Nature, 470, 531-534. in the study. Triviño M, Thuiller W, Cabeza M, Hickler T, Araújo MB Appendix S3. R scripts used to build the ‘cen- (In press) The Contribution of Vegetation and tral cluster’ consensus projections. Landscape Configuration for Predicting Appendix S4. Sources of uncertainty in the Environmental Change Impacts on Iberian Birds. ensemble of forecasts from individual BEMs. Velarde SJ, Malhi Y, Moran D, Wright J, Hussain S (2005) Valuing the impacts of climate change on protected Appendix S5. Non-analogue climate maps. areas in Africa. Ecological Economics, 53, 21-33. Appendix S6. Single-BEMs selected for the Venables WN, Ripley BD (2003) Modern Applied ‘central clusters’ and Anosim results. Statistics with S, xi+495pp. Springer-Verlag, New Appendix S7. Level of consensus among all York. BEMs and ‘central cluster’ BEMs. Webber BL, Yates CJ, Le Maitre DC et al. (2011) Modelling horses for novel climate courses: insights Appendix S8. Single-BEMs selected as ‘central from projecting potential distributions of native and model’ across species. alien Australian acacias with correlative and Appendix S9. mechanistic models. Diversity and Distributions, 17, Statistical differences in the dis- 978-1000. tributions of TSS and omission and commission Wiens JA, Stralberg D, Jongsomjit D, Howell CA, Snyder error among the five BEM consensus MA (2009) Niches, models, and climate change: projections. Assessing the assumptions and uncertainties. Proceedings of the National Academy of Sciences, 106, Appendix S10. Frequency distribution of spe- 19729-19736. cies turnover for alternative GCM clusters, Williams JW, Jackson ST, Kutzbach JE (2007) Projected emissions scenarios and BEM consensus distributions of novel and disappearing climates by methodologies. 2100 AD. Proceedings of the National Academy of Sciences, 104, 5738-5742. Appendix S11. Late-century species turnover Wisz MS, Guisan A (2009) Do pseudo-absence selection for alternative BEM consensus methodologies. strategies influence species distribution models and Appendix S12. Frequency distribution of cli- their predictions? An information-theoretic approach mate anomalies over the study area. based on simulated data. BMC Ecology, 9, 8. Appendix S13. Climate anomaly maps for the Wisz MS, Hijmans RJ, Li J, Peterson AT, Graham CH, Guisan A (2008) Effects of sample size on the three variables.
52 | Chapter I
Appendix S1: Principal Components Analysis to select predictor variables
21 predictor variables were computed using baseline monthly precipitation and temperature data for 1961-60 from the Climatic Research Unit (CRU, New et al., 2002). A principal components analysis (PCA) was performed on the values of these variables over the study area (N=1,851) in order to select a smaller set of variables and minimise co-linearity. The loadings for the first three principal components (PC1 to 3) are shown. The last three rows show the standard deviation, pro- portion of variance explained and cumulative variance explained by each PC.
Variables PC1 PC2 PC 3 mean temperature of the coldest month 0.022 0.365 0.016 mean temperature of the warmest month 0.277 0.194 -0.023 annual mean temperature 0.196 0.305 -0.028 annual sum of precipitation -0.295 0.161 0.095 annual mean precipitation -0.295 0.161 0.095 precipitation of the driest month -0.194 0.080 -0.558 precipitation of the wettest month -0.231 0.172 0.348 precipitation of the coldest month -0.143 0.171 0.036 precipitation of the warmest month -0.250 0.023 -0.208 mean temperature of the driest month 0.050 0.345 -0.027 mean temperature of the wettest month 0.268 0.161 -0.162 mean temperature of the coldest quarter 0.042 0.364 0.022 mean temperature of the warmest quarter 0.272 0.203 -0.038 precipitation of the driest quarter -0.209 0.097 -0.538 precipitation of the wettest quarter -0.246 0.158 0.343 precipitation of the coldest quarter -0.160 0.190 0.050 precipitation of the warmest quarter -0.270 0.020 -0.108 mean temperature of the driest quarter 0.068 0.355 -0.010 mean temperature of the wettest quarter 0.272 0.165 -0.147 temperature seasonality 0.159 -0.275 -0.057 precipitation seasonality 0.282 -0.010 0.167
Standard Deviation 2.936 2.642 1.301 Proportion of Variance Explained 0.410 0.332 0.081 Cumulative Proportion of Variance 0.410 0.743 0.823
Ensemble forecasting for African vertebrates | 53
Appendix S2: General Circulation Models used in the study
17 General Circulation Models (GCM) were used in this study. Their native resolution is shown in the table below, but downscaled projections at 10’ (Tabor & Williams, 2010) were used here. Model name abbreviations used in our study are also shown. The GCM clusters in which each GCM was included are given in the shaded columns for the three emissions scenarios (A2, A1B and B1).
Source Model Abbrev. Resolution Clusters y x A2 A1B B1 Bjerknes Centre for Climate Research BCCR-BCM2.0 bc2 2.80 2.80 2 1 3 National Center for Atmospheric Research CCSM3 cs3 1.41 1.41 2 1 1 Canadian Centre for Climate Modelling and CGCM3.1(T47) t47 3.75 3.75 3 3 3 Analysis Centre National de Recherches CNRM-CM3 cm3 2.80 2.80 3 3 3 Météorologiques, Météo-France Commonwealth Scientific and Industrial CSIRO-MK3.0 m30 1.90 1.90 1 1 1 Research Organisation (CSIRO) Atmospheric Research Commonwealth Scientific and Industrial CSIRO-MK3.5 m35 1.90 1.90 3 2 2 Research Organisation (CSIRO) Atmospheric Research Max Planck Institute for Meteorology ECHAM5/MPI- eh5 1.90 1.90 3 2 2 OM Meteorological Institute of the University of ECHO-G ecg 3.90 3.90 2 3 3 Bonn, Meteorological Research Institute of the Korea Meteorological Administration (KMA), and Model and Data Group Geophysical Fluid Dynamics Laboratory GFDL-CM2.0 cm2 2.00 2.50 3 2 3 (GFDL), National Oceanic and Atmospheric, Administration (NOAA), U.S. Department of Commerce Geophysical Fluid Dynamics Laboratory GFDL-CM2.1 c21 2.00 2.50 3 2 3 (GFDL), National Oceanic and Atmospheric, Administration (NOAA), U.S. Department of Commerce Goddard Institute for Space Studies (GISS), GISS-ER mer 4.00 5.00 2 3 3 National Aeronautics and Space Administra- tion (NASA) Institute for Numerical Mathematics INM-CM3.0 im3 4.00 5.00 2 1 3 Institut Pierre Simon Laplace IPSL-CM4 cm4 2.50 3.75 3 2 2 Center for Climate System Research (Univer- MIROC3.2 32m 2.80 2.80 1 3 3 sity of Tokyo), National Institute for Environ- (medium reso- mental Studies, and Frontier Research Center lution) for Global Change Meteorological Research Institute MRI-CGCM2.3.2 c23 2.80 2.80 1 1 1 National Center for Atmospheric Research PCM pc1 2.80 2.80 1 1 1 Hadley Centre for Climate Prediction UKMO-HadCM3 hd3 2.50 3.75 3 2 2
54 | Chapter I
Appendix S3: R scripts used to build the “central cluster” consensus projections
Available online.
Appendix S4: Sources of uncertainty in the ensemble of forecasts from individual BEMs
A point-wise three-way analysis of variance (ANOVA) without replication was performed for each grid cell over the study area (N=1,851) to evaluate the relative contributions of bioclimatic enve- lope models (BEM), the General Circulation Model clusters (GCMcons), the emissions scenarios (SRES), and interactions among these factors, to the variability in turnover forecasts for each taxon. The values in the table correspond to the median value and (in brackets) the lower and upper quar- tiles of the proportion of the total sum of squares attributed to each factor and are shown for both mid- and late-century, for each taxon.
Amphibians Snakes Mammals Birds 2041-60 BEMcons 76.5 (61.1-87.4) 81.7 (70.9-88.9) 76.2 (61.8-86.4) 79.8 (61.7-87.9) GCMcons 7.3 (2.7-16.2) 6.8 (3.0-13.4) 7.6 (3.7-16.5) 8.2 (4.2-19.1) SRES 3.3 (1.3-6.7) 3.1 (1.6-5.6) 4.0 (2.1-7.6) 4.0 (2.2-7.7) BEMcons:GCMcons 3.2 (1.6-5.8) 2.2 (1.3-4.1) 2.8 (1.5-5.4) 2.2 (1.2-4.3) BEMcons:SRES 1.4 (0.8-2.6) 1.1 (0.6-1.8) 1.2 (0.7-2.2) 1.0 (0.5-1.7) GCMcons:SRES 1.1 (0.4-2.2) 0.9 (0.4-1.7) 1.0 (0.5-2.2) 0.9 (0.4-1.9) BEMcons:GCMcons:SRES 2.2 (1.2-4.0) 1.6 (0.9-2.8) 1.8 (1-3.4.0) 1.2 (0.6-2.1) 2081-00 BEMcons 63.4 (38.7-83.1) 69.0 (48.2-82.5) 60.6 (37.3-76.6) 62.6 (34.1-76.7) GCMcons 7.4 (2.6-16.2) 7.5 (4.0-14.1) 9.6 (5.1-16.4) 10.2 (6.2-18.2) SRES 12.0 (3.8-30.9) 13.5 (6.1-27.6) 18.2 (9.5-33.7) 17.0 (9.1-35.2) BEMcons:GCMcons 2.4 (1.4-3.9) 1.7 (1.1-2.6) 1.9 (1.2-3.2) 1.8 (1.1-2.7) BEMcons:SRES 2.5 (1.5-4.3) 2.0 (1.2-3.1) 2.3 (1.4-3.8) 2.2 (1.4-3.6) GCMcons:SRES 1.4 (0.5-2.8) 0.9 (0.4-2.0) 1.1 (0.5-2.3) 1.0 (0.5-2.4) BEMcons:GCMcons:SRES 2.4 (1.4-4.2) 1.7 (1.1-2.8) 1.8 (1.1-3.0) 1.4 (0.9-2.3)
Ensemble forecasting for African vertebrates | 55
Appendix S5: Non-analogue climate maps
For both mid- and late-century, non-analogue climates correspond to the sum of the areas where mean temperature of the warmest month, mean temperature of the coldest month and annual precipitation values beyond the observed ranges in the baseline period are projected. Areas in red were projected to experience future climate conditions lacking baseline analogues. The percentage of the total area projected to experience non-analogue climates is indicated in red and the same percentage for each variable (tcm=temperature of the coldest month, tw=temperature of the warmest month, ps=annual precipitation) in black underneath each map. Data are shown for the three General Circulation Model clusters under each emissions scenario (A2, A1B and B1).
56 | Chapter I
Appendix S6: Single-BEMs selected for the ‘central clusters’ and Anosim results
For each taxon and climate projection, the seven bioclimatic envelope models were clustered based on similarities in magnitude and spatial pattern between single-model probabilistic projections and the multi-model median probabilistic projection for late-century. The k-means clustering re- sulted in four clusters, shown below for the three General Circulation Model clusters under each emission scenario (A2, A1B and B1). Cluster 1 is the maximum consensus ‘central cluster’. The analysis of similarity (Anosim) test was used to test the statistical significance of the clustering.
BEM clusters Anosim SRES GCM ANN GAM GBM GLM MARS FDA RF Statistics P-value A2 Cluster 1 4 2 1 1 3 3 3 0.912 0.004 Cluster 2 4 2 1 1 3 3 3 0.765 0.013 Cluster 3 2 4 1 1 3 3 1 0.971 0.004 A1B Cluster 1 4 2 1 2 1 3 1 0.853 0.004 Cluster 2 4 2 1 1 3 3 3 0.853 0.004
Amphibians Cluster 3 2 4 1 1 3 3 3 0.735 0.027 B1 Cluster 1 4 2 1 1 3 3 3 0.824 0.008 Cluster 2 4 2 1 2 1 3 1 0.853 0.004 Cluster 3 4 2 1 1 3 3 3 0.912 0.004
A2 Cluster 1 3 2 1 1 4 4 4 0.853 0.004 Cluster 2 2 4 1 1 3 3 3 0.882 0.004 Cluster 3 2 4 1 1 3 3 1 0.941 0.004
Snakes A1B Cluster 1 3 2 1 1 4 4 4 0.853 0.004 Cluster 2 2 4 1 1 3 3 3 0.853 0.004 Cluster 3 2 4 1 1 3 3 3 0.882 0.004 B1 Cluster 1 3 2 1 1 4 4 4 0.824 0.002 Cluster 2 3 2 1 1 4 4 4 0.853 0.004 Cluster 3 3 2 1 1 4 4 4 0.882 0.004
A2 Cluster 1 4 1 2 1 2 3 2 0.853 0.003 Cluster 2 4 1 2 1 3 2 2 0.941 0.003 Cluster 3 4 2 2 1 3 2 2 0.912 0.003 A1B Cluster 1 4 1 2 1 2 3 2 0.794 0.006
Mammals Cluster 2 4 1 2 1 2 3 2 0.618 0.045 Cluster 3 4 1 1 2 3 3 1 1.000 0.003 B1 Cluster 1 4 1 2 1 2 3 2 0.933 0.034 Cluster 2 4 1 2 1 2 3 2 0.765 0.016 Cluster 3 4 1 2 1 2 3 2 1.000 0.003
A2 Cluster 1 2 1 1 1 4 3 1 0.933 0.032 Cluster 2 2 1 1 1 4 3 1 0.911 0.032
Birds Cluster 3 2 1 1 1 4 3 1 1.000 0.003 A1B Cluster 1 2 1 1 1 4 3 1 0.978 0.032 Cluster 2 2 1 1 1 4 3 1 0.933 0.032 Cluster 3 2 1 1 1 4 3 1 0.882 0.011 B1 Cluster 1 4 1 2 1 2 3 2 0.911 0.032 Cluster 2 2 1 4 1 4 3 4 0.889 0.032 Cluster 3 2 1 1 1 4 3 1 0.933 0.032
Ensemble forecasting for African vertebrates | 57
Appendix S7: Level of consensus among all BEMs and ‘central cluster’ BEMs
To evaluate the degree of consensus achieved in bioclimatic envelope model (BEM) consensus projections, Principal Components Analysis (PCA) was performed for each species on the probabil- istic projections. The proportion of variance explained by the first principal component axis re- flects the degree of consensus among projections. PCAs were performed on baseline and late- century projections, both on the seven single-BEMs and on the ‘central cluster’ BEMs. The graphs show the frequency distributions across all amphibian (n=284), snake (n=310), mammal (n=623) and bird (n=1,506) species, with the degree of consensus increasing on a scale from 0 to 1. For future projections calculations, the ‘maximum consensus’ General Circulation Model cluster (clus- ter 2, the ‘central cluster’) under the emissions scenario A1B was used.
Appendix S8: Single-BEMs selected as ‘central model’ across species
With the ‘central model’ methodology to combine the ensembles of bioclimatic envelope models (BEM), we selected the model summarising the highest amount of variation among projections for each species. A Principal Components Analysis (PCA) was performed for each species on the late- century projected probabilities and the ‘central model’ corresponded to the one with the highest PCA loading in the first (consensus) axis. The graphs show the percentage of species for which each of the seven single BEMs was selected as the ‘central model’. Results are shown for the projections built with 100% of the data.
58 | Chapter I
Appendix S9: Statistical differences in the distributions of TSS and omission and commission error among the five BEM consensus projections
Pair-wise differences in True Skill Statistics (TSS), omission and commission error frequencies for the five bioclimatic envelope model (BEM) consensus projections for amphibian (n=284), snake (n=310), mammal (n=623) and bird (n=1,506) species were tested using the Wilcoxon signed rank test. EMean is the ensemble mean, EWMean the ensemble weighted mean, EMed the ensemble median, CMod the central model, and CClus the central cluster. For each taxon and accuracy meas- ure, Bonferroni-corrected P-values (number of tests=10) are shown above the diagonal and the test statistics below.
Method Omission error Commission error EMean EWMean EMed CMod CClus EMean EWMean EMed CMod CClus EMean 1.70E-02 3.36E-11 1.37E-01 8.15E-09 7.89E-12 5.70E-35 2.63E-19 5.47E-47 EWMean 4049 6.62E-15 1.13E-03 7.08E-05 9519 8.73E-15 6.73E-09 1.57E-45 EMed 13681 2802 2.57E-01 1.33E-18 2689 29204 1.07E+00 4.57E-41 CMod 8373 7593 13164 2.36E-13 32400 28041 21704 2.02E-21
Amphibians CClus 14402 12482 20381 16303 40149 39831 38580 33304
EMean 1.19E+00 3.83E-12 1.81E-01 1.03E-01 3.27E-04 2.30E-29 3.04E-15 2.14E-51 EWMean 9137 3.45E-13 2.96E-02 4.01E-02 16091 1.51E-13 6.91E-09 3.99E-51 EMed 20619 6024 2.63E-02 9.26E-03 5680 35393 4.65E+00 3.97E-45 Snakes CMod 13201 11260 20768 3.46E+00 35288 32048 24302 4.80E-33 CClus 12536 12079 19925 15602 47891 47827 46097 41596
EMean 2.05E-21 2.36E-08 1.84E-41 8.97E-12 6.93E+00 4.45E-15 4.88E-06 2.57E-97 EWMean 9805 5.53E-38 9.47E-57 3.85E+00 77856 1.47E-09 1.18E-02 1.71E-97 EMed 37665 10056 6.88E-31 5.21E-28 49203 110835 8.24E-01 3.60E-99
Mammals CMod 18247 12514 24231 7.71E-62 113297 102932 95942 2.87E-66 CClus 64478 41076 77061 125512 189618 188631 189976 167655
EMean 1.62E-52 9.46E-16 4.79E-82 3.13E-33 4.00E-02 4.20E-27 1.10E-04 9.08E-241 EWMean 79127 1.08E-84 1.32E-132 9.92E-01 435115 3.27E-09 2.67E-01 9.77E-238 Birds EMed 239289 80336 7.06E-65 3.40E-70 294998 591882 5.57E+00 1.19E-239 CMod 136689 73731 156865 1.53E-134 593707 559791 503122 1.39E-189 CClus 404898 274275 530221 738702 1108284 1098301 1104396 1037136
Method TSS EMean EWMean EMed CMod CClus EMean 4.31E-06 7.55E-01 2.78E-02 1.55E-44 EWMean 12609 2.54E-02 9.44E+00 3.89E-43 EMed 17772 15934 1.43E+00 1.56E-43 CMod 24215 19997 22112 1.39E-34
Amphibians CClus 39855 39539 39629 37250
EMean 7.54E-02 4.88E-01 1.11E+00 7.77E-43 EWMean 19213 3.59E+00 9.37E+00 3.64E-35 EMed 20990 22359 5.44E+00 1.43E-39 Snakes CMod 25957 23608 24426 1.56E-30 CClus 46034 43936 45166 41583
EMean 2.11E-47 6.11E-01 2.96E-54 2.63E-77 EWMean 27902 9.20E-39 1.71E-85 1.44E-55 EMed 98747 35386 1.43E-58 2.29E-82
Mammals CMod 25971 7915 23224 1.63E-100 CClus 180356 166971 183086 193325
EMean 1.15E-118 1.15E-02 1.54E-140 6.55E-218 EWMean 159544 2.88E-98 2.10E-208 3.64E-173 Birds EMed 583286 198652 8.45E-141 1.41E-223 CMod 134970 44482 135868 2.66E-241 CClus 1090226 1032510 1102239 1125877
Ensemble forecasting for African vertebrates | 59
Appendix S10: Frequency distribution of species turnover for alternative GCM clusters, emissions scenarios and BEM consensus projections
Frequency distribution of species turnover rates over the study area (N=1851) for amphibians, snakes, mammals and birds, for mid- and late-century. Data are shown for the five bioclimatic en- velope model consensus projections (EMean=ensemble mean, EWMean=ensemble weighted mean, EMed=ensemble median, CMod=central model, and CClus=central cluster), and for the three Gen- eral Circulation Model clusters under emissions scenarios A2, A1B and B1.
60 | Chapter I
Ensemble forecasting for African vertebrates | 61
Appendix S11: Late-century species turnover for alternative BEM consensus methodologies
Projected late-century turnover rate (%) for amphibian, snake, mammal and bird species for alter- native BEM consensus projections. The five consensus projections are: ensemble mean (EMean), ensemble weighted mean (EWMean), ensemble median (EMed), central model (CMod), and central cluster (CClus). Projections refer to the ‘maximum consensus’ General Circulation Model cluster (cluster 2) under the A1B emission scenario.
62 | Chapter I
Appendix S12: Frequency distribution of climate anomalies over the study area
Boxplots of mean temperature of the warmest and coldest month and annual precipitation anoma- lies over the study area (N=1,851). Data are shown for the mid- (white bars) and late-century (grey bars), for the three General Circulation Model clusters under emissions scenarios A2, A1B and B1.
Ensemble forecasting for African vertebrates | 63
Appendix S13: Climate anomaly maps for the three variables
Difference between future (mid- and late-century) projections and baseline data for the three vari- ables (annual precipitation and temperature of the coldest and warmest months), for the different General Circulation Model clusters (1 to 3) under alternative emissions scenarios (A2, A1B and B1). Note that for annual precipitation, the scale (in mm) covers both negative values (blue tones, refer- ring to decreased precipitation in the future) and positive values (red tones, referring to increased precipitation in the future), whereas for the temperature-based variables the scale (in degrees Celsius) has positive values only (blue through to red tones indicating increasingly warmer future temperatures).
References
Clarke KR, Warwick RM (1994) Similarity-based testing Tabor K, Williams JW (2010) Globally downscaled for community pattern: the two-way layout with no climate projections for assessing the conservation replication. Marine Biology, 118, 167-176. impacts of climate change. Ecological Applications, 20, 554-565. Duan Q, Phillips TJ (2010) Bayesian estimation of local signal and noise in multimodel simulations of climate Venables WN, Ripley BD (2003) Modern Applied change. J. Geophys. Res., 115, D18123. Statistics with S, New York, Springer-Verlag. New M, Lister D, Hulme M, Makin I (2002) A high- resolution data set of surface climate over global land areas. Climate Research, 21, 1-25.
64 | Chapter I
Addendum
Figure 7 in Garcia et al. (2012) shows the per- alternative figure might be perceived as provid- centage of species projected to retain climatic ing more useful information for examining the suitability under each emissions scenario. For regional impacts of climate change on climatic each location i, the proportion R of species suitability for sub-Saharan African vertebrates. retaining climatic suitability is given by , where ri is the number of species that retain climatic suitability in loca- tion i under climate change and N is the total number of species in sub-Saharan Africa in the baseline period. Differences across taxa and across emissions scenarios are visible. But because the proportions of retention through- out the region are relative to the same total number of species, it is the areas of high base- line richness that are emphasised. As this figure may lead to misinterpreta- tions, we show below an alternative figure where in situ retention is calculated in relation to local (rather than total) species richness. In this case, the proportion R of species retaining climatic suitability in location i is given by , where ri is the number of species that retain climatic suitability in loca- tion i under climate change and ni is the num- Percentage of species predicted to retain cli- ber of species in location i in the baseline sce- matic suitability under each emissions scenario. nario. Here the variability among scenarios is The proportion of the local numbers of species still clear, as is the greater potential impact of of amphibians, snakes, mammals and birds that climate change on in situ climatic suitability for are projected to retain climatic suitability in amphibians and snakes in comparison to each location are shown for the median ensem- ble of all bioclimatic envelope models and for mammals and birds. However, by normalising the ‘maximum consensus’ general circulation the retention of climatic suitability to local model cluster under the A2, A1B and B1 emis- species richness, the alternative maps empha- sions scenarios. sise the areas where higher proportions of local richness are projected to retain climatic suita- bility, irrespective of total richness. As such, the
Chapter II
Matching species traits to projected threats and opportunities from climate change
RAQUEL A. GARCIA, MIGUEL B. ARAÚJO, NEIL D. BURGESS, WENDY B. FODEN, ALEXANDER GUTSCHE, CARSTEN RAHBEK, AND MAR CABEZA, Journal of Biogeography doi:10.1111/jbi.12257 (2014)
Matching species traits to projected threats and opportunities from climate change
RAQUEL A. GARCIA1,2,3, MIGUEL B. ARAÚJO1,2,3,4 , NEIL D. BURGESS2,5,6, WENDY B. FODEN7,8, ALEXANDER GUTSCHE9, CARSTEN RAHBEK2,4, and MAR CABEZA10
1 Center for Macroecology, Evolution and Climate, Department of Biology, University of Copenhagen, Denmark 2 Department of Biogeography and Global Change, National Museum of Natural Sciences, CSIC, Madrid, Spain 3 InBio/CIBIO, University of Évora, Évora, Portugal 4 Imperial College London, Silwood Park, Ascot, Berkshire, United Kingdom 5 WWF US Conservation Science Program, Washington, DC, US 6 United Nations Environment Programme World Conservation Monitoring Centre, Cambridge, UK 7 Global Species Programme, International Union for Conservation of Nature (IUCN), Cambridge, UK 8 Animal, Plant and Environmental Sciences, University of the Witwatersrand, Johannesburg, South Africa 9 Museum für Naturkunde, Leibniz Institute for Research on Evolution and Biodiversity, Berlin, Germany 10 Metapopulation Research Group, Department of Biosciences, University of Helsinki, Finland
Journal of Biogeography doi:10.1111/jbi.12257 (2014)
Abstract
Aim Climate change can lead to decreased climatic all species, and individually for groups of species suitability within species’ distributions, increased with different combinations of threats and oppor- fragmentation of climatically suitable space, tunities. and/or emergence of newly suitable areas outside Results In the Congo Basin and arid Southern present distributions. Each of these extrinsic Africa, projected losses for wide-ranging amphibi- threats and opportunities potentially interacts ans were compounded by sensitivity to climatic with specific intrinsic traits of species, yet this variation, and expected gains were precluded by specificity is seldom considered in risk assess- poor dispersal ability. The spatial overlap be- ments. We present an analytical framework for tween exposure and vulnerability was more pro- examining projections of climate change-induced nounced for species projected to have their cli- threats and opportunities with reference to traits mate space contracting in situ or shifting to distant that are likely to mediate species’ responses, and geographical areas. Our results exclude the poten- illustrate the applicability of the framework. tial exposure of narrow-ranging species to shrink- Location Sub-Saharan Africa. ing climates in the African tropical mountains. Methods We applied the framework to 195 sub- Main conclusions We illustrate the application of Saharan African amphibians with both available a framework combining spatial projections of bioclimatic envelope model projections for the climate change exposure with traits that are likely mid-21st century and trait data. Excluded were to mediate species’ responses. Although the pro- 500 narrow-ranging species mainly from montane posed framework carries several assumptions areas. For each of projected losses, increased that require further scrutiny, its application adds a fragmentation and gains of climate space, we degree of realism to familiar assessments that selected potential response-mediating traits and consider all species to be equally affected by cli- examined the spatial overlap with vulnerability mate change-induced threats and opportunities. due to these traits. We examined the overlap for
68 | Chapter II
Introduction of local populations (e.g. Sinervo et al., 2010). Second, even where suitable climate space Fingerprints of recent climate change impacts persists, it may become more fragmented. Ar- on species’ distributions are already apparent eas of suitable climate may lose contiguity, with (Parmesan & Yohe, 2003), but predicting future fragments becoming more isolated and smaller climate change impacts is still a major scientific in area. Third, gains of climatic suitability out- challenge (Pereira et al., 2010). Predictions of side current distributions of species bring po- species’ exposure to climate change have to tential opportunities for colonization. These date relied mostly on bioclimatic envelope components of exposure are likely to have models, and are increasingly being comple- distinct spatial distributions and conservation mented with available trait data to estimate implications, but are seldom teased apart (but species’ vulnerability (e.g. Heikkinen et al., see Midgley et al., 2003; Heikkinen et al., 2009; 2009). Assessing the risk of species’ extinction Araújo et al., 2011). or decline requires an understanding of threats, Species also vary in their degree of intrinsic or extrinsic adverse events, and vulnerability, sensitivity and capacity to adapt to exposure or the intrinsic susceptibility of species to (Williams et al., 2008; Chevin et al., 2010), and threats (Araújo & Williams, 2000). It is thus traits can capture this variation. Traits are important to understand the interaction be- meant as ‘any morphological, physiological or tween threats and vulnerability, as specific phenological feature measurable at the individ- traits are likely to mediate species’ responses to ual level, from the cell to the whole-organism different threats (Isaac & Cowlishaw, 2004; level, without reference to the environment or Fritz et al., 2009; Murray et al., 2011; González- any other level of organization’ (sensu Violle et Suárez et al., 2013). For example, under habi- al., 2007, p. 884). Examples are limb or wing tat-modifying processes such as agriculture and length, and clutch size. In many cases, ecologi- logging, small-sized, habitat specialist mam- cal characteristics of species and their geo- mals are most affected, whereas under proc- graphical or environmental ranges are used as esses that directly affect survival, such as hunt- proxies for traits, in the expectation that they ing, the most susceptible are mammals with summarize combinations of traits. Examples large body size and small litter size (González- are species’ geographical range size and cli- Suárez et al., 2013; see also Isaac & Cowlishaw matic breadth. 2004). In climate change risk assessments, Previous research has identified general however, this specificity in the interaction be- traits that predispose species to extinction (e.g. tween threat and vulnerability has not been Purvis et al., 2000), and specific traits that me- addressed sufficiently. diate the effect of particular threats on species In studies using bioclimatic envelope mod- (Isaac & Cowlishaw, 2004; Murray et al., 2011; els, the level of species’ exposure to climate González-Suárez et al., 2013). Under changing change is commonly inferred from temporal climates, poor dispersal ability and habitat or changes in the overall size of species’ climati- climatic specialization, for example, have been cally suitable space (e.g. Thuiller et al., 2005b; suggested to increase vulnerability (Peters & Araújo et al., 2006; Huntley et al., 2006; Feeley Darling, 1985), and have been found to corre- et al., 2012; Triviño et al., 2013). Such summary late positively with empirical data on range measures conceal different opportunities as contractions (Beaumont & Hughes, 2002; Botts well as threats, each imposing specific con- et al., 2013). Generally, traits have accounted straints on species (Thomas et al., 2011). First, for a significant but small amount of the varia- loss of climatic suitability within existing dis- tion in climate change-induced range shifts tributions is expected to affect the persistence
Species traits and climatic threats and opportunities | 69
(Buckley & Kingsolver, 2012). Yet, changes in exacerbate projected losses and increased the size, level of fragmentation and position of fragmentation of climatically suitable areas, or species’ climate space each represent a distinct restrain projected gains of newly suitable areas. threat or opportunity under changing climates, We test the practicality of the framework and thus are likely to interact with particular on wide-ranging sub-Saharan African amphibi- sets of ‘response-mediating’ traits (sensu Luck ans, using available bioclimatic envelope model et al., 2012). and trait data. Our application of the frame- Besides evolution on longer time-scales, work reflects situations where only proxies for three main determinants of climate change traits are available, and where bioclimatic enve- vulnerability can be distinguished (Williams et lope model projections are available only for al., 2008; Chevin et al., 2010). First, traits de- wide-ranging species, in our case leading to the scribing plasticity of individual phenology, exclusion of most narrow-ranging endemics behaviour or physiology affect the potential of typical of biodiverse montane regions. World- individuals to persist in situ under changing wide, amphibian populations are declining due climates. For example, species able to physio- to a multitude of threats that include habitat logically tolerate a wide range of climatic varia- destruction, climate change and the pathogeni- tion (Huey et al., 2012), or adapt their behav- cal fungal disease chytridiomycosis (Blaustein iour to lessen exposure to unsuitable climates & Kiesecker, 2002; Hof et al., 2011; Li et al., (Chown, 2012), have, all else being equal, 2013). Climate change, often in tandem with higher chances of persistence in situ under land-use change, is expected to affect large climate change. Second, traits influencing the areas of tropical Africa in the future (Hof et al., potential of individuals to disperse affect their 2011; Foden et al., 2013). While our aim here is capacity to colonize newly suitable environ- to present a spatially explicit framework for ments (Pöyry et al., 2009) or move between linking threats or opportunities and vulnerabil- fragments of suitable climate. Third, life-history ity in climate-change risk assessments, the traits influencing population growth, although practical application we show can also contrib- not affected by environmental change, may ute to a better understanding of climate change constrain the rate of dispersal or in situ adapta- risks facing wide-ranging sub-Saharan African tion. In the case of reproductive traits, for ex- amphibians. ample, frequent or early reproduction and high fecundity, should increase colonization oppor- tunities (Angert et al., 2011). Material and Methods Here, we present an analytical framework for examining projections of climate change- Study species induced threats or opportunities for species Our study sample consists of 195 amphibian with reference to the vulnerability of species. species (see Appendix S1 in Supporting Infor- Each threat or opportunity – exposure to loss, mation) restricted in their distributions to sub- fragmentation and gain of climate space – is Saharan Africa, and with available projections matched to specific response-mediating traits of species bioclimatic envelope models (Garcia that potentially render species vulnerable, et al., 2012) and trait data (Foden et al., 2008, according to expectations from theory and 2013). From all 695 species in the original empirical evidence (Fig. 1). Areas of spatial species distribution (Hansen et al., 2007) and overlap between threats or opportunities and trait (Foden et al., 2008, 2013) databases, only associated vulnerability are then identified. 272 had model projections available (excluded These are areas where traits can potentially were 423 species with fewer than 15 gridded
70 | Chapter II
Exposure to
seline climate space (left) to -mediating traits associated with spe- across taxa. by dashed lines, and the component of exposure under consider- responses to climate change-induced threats and opportunities. tunities, examples are given of potential response areas are each illustrated with diagrams of the shift from ba imate space represented mography, where supported by theoretical and empirical studies | Examples of traits that potentially mediate species’ Figure 1 loss, increased fragmentation and gain of climatically suitable cl baseline the with lines, closed (right, space climate future ation highlighted in black). For each of these threats or oppor cies’ plasticity, dispersal ability and/or de
Species traits and climatic threats and opportunities | 71
occurrence records), and, of these, 195 also had First, for each species we considered local data for all traits (excluded were 77 species losses in pixels projected to be climatically with some trait data missing). Our sample is, suitable in the baseline period but unsuitable in therefore, restricted to the widest-ranging the future. Second, for fragmentation we used a species in the dataset, and may not be repre- distance-based measure of contagion of suit- sentative of the overall taxonomic and geo- able climate space. Contagion was measured as graphical amphibian diversity patterns in Af- the weighted average of the number of suitable rica (Appendix S2a–c). In particular, our data pixels among a set of ki neighbours of a central exclude most species endemic to biodiverse pixel yi, where the weight given to the grid cell montane areas such as the Cameroon/ Nigerian yj is wij = 1/dij, and dij is the great-circle distance highlands and the Eastern Afromontane biodi- between grid cells yi and yj (Araújo et al., 2002; versity hotspot, and most species in the highest equation 1). Owing to poor dispersal ability of threat categories of the IUCN Red List (Appen- most amphibians, we considered only the first- dix S2d). order neighbours (maximum = 8) adjacent to the central pixel. Changes in contagion were Extrinsic threats and opportunities from climate given by the difference between future and change baseline contagion, with negative values indi- cating reduced contagion, i.e. increased frag- We used published baseline (1961–90) and mentation. mid-century (2041–60) projections of climati- cally suitable areas for our 195 amphibians (1) species, at one degree latitudinal–longitudinal resolution (c. 111 km × 111 km at the equator), according to mean temperatures of warmest and coldest month and annual precipitation Third, local gains corresponded to pixels (for detailed methods see Garcia et al., 2012). projected to be climatically unsuitable in the Future projections were for a multi-model baseline period but suitable in the future. climate ensemble under the A1B greenhouse Fourth, for species projected to gain newly gas emissions scenario (Nakicenovic & Swart, suitable areas in mid-century, we also com- 2000). We used consensus projections obtained puted, for each pixel of newly suitable climate, by computing the median among seven biocli- the great-circle distance to the nearest pixel of matic envelope modelling techniques, in pres- baseline suitable climate. ence–absence format, and assuming unlimited dispersal in future projections. Here, to charac- Intrinsic vulnerability to climate change terize climate change-induced threats and op- We sourced the trait data from the IUCN’s trait- portunities, we computed four types of pro- based climate change vulnerability assessment jected changes in climatic suitability: local loss, for amphibians (Foden et al., 2008, 2013). fragmentation and gain of climatic suitability, These data are mainly ecological characteristics and distance to new areas gained. Each metric of species or their ranges rather than traits in a was quantified per pixel, as described below, strict sense (sensu Violle et al., 2007). From the yielding maps of changes for each species. data available, we selected characteristics that Composite maps for our species sample were are likely to summarize response-mediating also obtained by summing the number of spe- traits under climate change (Fig. 1), and that cies with projected local loss, increased frag- are not strongly correlated. We thus selected mentation, or gain of suitable climate in each tolerance to temperature and precipitation pixel. change, dependence on precipitation cues,
72 | Chapter II
dispersal ability and reproductive output (Ta- traits and exposure (Fig. 1), we matched plas- ble 1; see Foden et al., 2008, 2013, for details). ticity characteristics to local losses, dispersal For each of these characteristics, we used the characteristics to increased fragmentation, and climate change vulnerability classification of both dispersal and demography characteristics species by Foden et al. (2008, 2013): species to gains (see Table 1). For climatic tolerance, were assigned ‘high’, ‘lower’ or ‘unknown’ we considered for each species either the toler- scores of vulnerability based on the ranking of ance to temperature or to precipitation de- all sub-Saharan African species in the IUCN pending on the variable of highest importance database or on pre-defined criteria (see Table in the bioclimatic envelope models (assessed 1). with permutations within the BIOMOD computing platform in R; Thuiller et al., 2009). Spatial overlap between threats or opportunities and We thus identified where projected losses were vulnerability for species with high vulnerability according to plasticity traits, where projected fragmentation We spatially assessed where climate-induced was for species with high vulnerability accord- threats or opportunities, defined with the met- ing to dispersal traits, and where projected rics of species’ exposure to climate change, gains were for species with high vulnerability overlapped with high climate change vulner- according to dispersal or demography traits. ability of species according to our selected The assessment of spatial overlap was first traits. Following expected interactions between
Table 1 | Species traits used as estimates of intrinsic vulnerability to climate change exposure for sub-Saharan African amphibians. Three climate change-induced threats and opportunities for species were defined based on bioclimatic envelope model projections (Garcia et al., 2012): loss, increased frag- mentation, and gain of suitable climate space. For each, different sets of species traits or characteristics of species and their ranges (Foden et al., 2008, 2013) were selected that are likely to mediate species’ responses (see Fig. 1).
Intrinsic vulnerability to climate change
Plasticity . Tolerance to temperature change: average absolute deviation for all cells in species’ refined range for the monthly means; high risk if ≤ 1.06 °C, i.e. 25% of all 704 sub-Saharan African spe- cies in the IUCN dataset with the narrowest tolerance ranges. . Tolerance to precipitation change: average absolute deviation for all cells in species’ refined
Loss range for the monthly means; high risk if ≤ 46.89 mm, i.e. 25% of all 704 sub-Saharan African species in the IUCN dataset with the narrowest tolerance ranges. . Dependence on environmental cues: high risk if dependent on rainfall or increased water avail- ability for mass breeding (excludes species buffered by occurring in forests).
Dispersal . Dispersal ability: high risk if not known to have become established outside their natural ranges, not associated with flowing water, and have small ranges (≤ 4000 km2, i.e. 25% of all 704 sub- tation Saharan African species in the IUCN dataset with the smallest ranges). Fragmen-
Dispersal . Dispersal ability: high risk if not known to have become established outside their natural ranges, not associated with flowing water, and have small ranges (≤ 4000 km2, i.e. 25% of all 704 sub- Saharan African species in the IUCN dataset with the smallest ranges). Gain Demography Extrinsic threats and opportunities from climate change . Reproductive output: high risk if ≤ 50 offspring annually (where known) or viviparous.
Species traits and climatic threats and opportunities | 73
conducted for all species in our sample, and distant areas; and ‘expanding’ species were then individually for groups of species pro- projected to retain most of their baseline suit- jected to experience different combinations of ability and have large and distant gains outside threats and opportunities. To identify these their distributions. groups, we ranked species according to each species’ overall projected losses, gains and distance to gains of suitable climate. For each Results species, the overall local losses of climatic suit- Our results concern 195 sub-Saharan African ability L were quantified using the proportion amphibians with wider ranges, significantly of baseline area of suitable climate (ai) pro- larger temperature tolerances (P-value < 0.05, jected to be lost in the future (li; equation 2). Student’s t-test), and a larger proportion of The overall opportunity for gains G was meas- species with higher reproductive output than ured as the proportion of baseline area of suit- the remaining 500 species in the dataset (Ap- able climate (ai) projected to be gained in the pendix S2c,e). For this subset, the spatial pat- future (gi; equation 3). The surface area across terns of climate change exposure (greyscale the A pixels of the study area was measured maps in Fig. 2) showed widespread local losses taking into account the curvature of the Earth. of climatically suitable space, but a concentra- The overall distance to new areas D for each tion of increased fragmentation in montane species was the mean of the minimum great- areas and local gains in the Congo Basin. Spe- circle distances dib between each pixel gained i cies were also not randomly distributed across and the baseline suitable areas b computed sub-Saharan Africa with regard to climate across the N pixels gained (equation 4). change vulnerability due to selected traits or