Shifting suitability for malaria vectors across Africa with warming climates
BMC Infectious Diseases volume 9, Article number: 59 (2009)
Climates are changing rapidly, producing warm climate conditions globally not previously observed in modern history. Malaria is of great concern as a cause of human mortality and morbidity, particularly across Africa, thanks in large part to the presence there of a particularly competent suite of mosquito vector species.
I derive spatially explicit estimates of human populations living in regions newly suitable climatically for populations of two key Anopheles gambiae vector complex species in Africa over the coming 50 years, based on ecological niche model projections over two global climate models, two scenarios of climate change, and detailed spatial summaries of human population distributions.
For both species, under all scenarios, given the changing spatial distribution of appropriate conditions and the current population distribution, the models predict a reduction of 11.3–30.2% in the percentage of the overall population living in areas climatically suitable for these vector species in coming decades, but reductions and increases are focused in different regions: malaria vector suitability is likely to decrease in West Africa, but increase in eastern and southern Africa.
Climate change effects on African malaria vectors shift their distributional potential from west to east and south, which has implications for overall numbers of people exposed to these vector species. Although the total is reduced, malaria is likely to pose novel public health problems in areas where it has not previously been common.
Malaria is a vector-borne anthroponosis, transmitted in large part by Anopheles mosquitoes, that endangers more than 2.5 × 109 humans annually . Its transmission cycle has been modeled in great detail [2, 3], but application of such process-based models has generally been limited to local and regional scales [4, 5], given challenges in spatially explicit parameter estimation. An alternative approach that offers broadest applicability is that of focusing on vector species' geographic distributions via ecological niche modeling techniques now well tested [6–8] within frameworks for reconstructing the geographic dimensions of disease transmission .
Implications of climate change for malaria transmission across Africa have been the subject of numerous commentaries [10–12] and a few attempted analyses [13–18]. One analysis  developed a global picture of climate change effects on malaria transmission, but relied on country-level classifications of risk, and did not discern the fine details of mosquito-species-specific distributional shifts. To date, then, specific climate change scenarios have yet to be applied to derivation of detailed projections of changing patterns of potential for malaria transmission based on mosquito vector species' distributions across Africa.
Ecological niche modeling uses nonrandom associations between known occurrences of species and digital geospatial data layers summarizing relevant environmental parameters to reconstruct the suite of environmental conditions under which the species is able to maintain populations without immigrational subsidy . Modeled niches are generally conservative over considerable time periods [20, 21], offering considerable predictive power regarding distributional ecology of species, including when environmental conditions are changing [8, 22]. Previous studies have used these techniques to illuminate distributions of mosquito species under diverse circumstances [17, 23–28]; perhaps most importantly, mosquito species' distributions are known to be highly sensitive to climatic variations, including on very short time scales . Here, I integrate detailed climate change scenarios for the two most significant Africa-wide mosquito vector species based on ecological niche models and detailed spatial summaries of present-day human population distribution [29, 30] to calculate shifting patterns of potential human exposure to malaria across Africa to be expected over the coming half-century.
Occurrence data for vector species were derived from a recent detailed compilation  that was the basis of previous analyses [16, 28]; in the present study, I focused on the two species with broad, continentwide distributions in the Anopheles gambiae complex: A. gambiae sensu stricto and A. arabiensis, for which 581 and 501 occurrence records were available, respectively. These data are based on sampling that was intense in some countries, and negligible or absent in others – as such, an extrapolative approach to their analysis becomes key. As such, I related the occurrence data to a suite of environmental data including dimensions of climate (annual mean temperature, mean monthly maximum temperature, mean monthly minimum temperature, annual precipitation; from ) and topography (slope, aspect, compound topographic index; from ). The relationship between the occurrence data and the environmental data sets – effectively the basis for the niche model – thus interpolates into unsampled regions based on the environmental characteristics of those regions, and thereby offers improved predictive ability regarding distributions. All environmental data sets were resampled to 0.1° spatial resolution for analysis.
Ecological niche models were generated using the Genetic Algorithm for Rule-Set Prediction (GARP) , a niche-modeling approach that has been the basis for most previous niche modeling applications to disease geography . GARP is an evolutionary-computing method that builds ENMs based on non-random associations between known occurrence points for species and sets of GIS coverages describing the ecological landscape. Occurrence data are used as follows: 50% of occurrence data points are set aside for an independent test of model quality (extrinsic testing data), 25% are used for developing models (training data), and 25% are used for tests of model quality internal to GARP (intrinsic testing data). Distributional data are converted to raster layers by the GARP program; then, by random sampling from areas of known presence (training and intrinsic test data) and areas of 'pseudoabsence' (areas lacking known presences), two data sets are created, each of 1250 points; these data sets are used for rule generation and model testing, respectively.
The first rule is created by applying a method chosen randomly from a set of inferential tools (i.e., atomic rules that specify particular environmental value combinations as suitable, logistic regression, bioclimatic range rules, negated bioclimatic range rules ). The genetic algorithm consists of specially defined operators (e.g., crossover, mutation) that modify the initial rules, and thus the result are models that have "evolved" – after each modification, the quality of the rule is tested (to maximize both significance and predictive accuracy) and a size-limited set of best rules is retained. Because rules are tested based on independent data (intrinsic test data), performance values reflect the expected performance of the rule, an independent verification that gives a more reliable estimate of true rule performance. The final result is a set of rules that can be projected onto a map to produce a potential geographic distribution for the species under investigation.
Following recent recommendations , for each species, I developed 100 replicate random-walk GARP models, and filtered out 90% based on consideration of error statistics, as follows. The 'best subsets' methodology consists of an initial filter removing models that omit (omission error = predicting absence in areas of known presence) heavily based on the extrinsic testing data, and a second filter based on an index of commission error (= predicting presence in areas of known absence), in which models predicting very large and very small areas are removed from consideration. Specifically, in DesktopGARP, I used a "soft" omission threshold of 20%, and 50% retention based on commission considerations; the result was 10 'best subsets' models (binary raster data layers) that were summed to produce a best estimate of geographic prediction. I took as a prediction of suitable conditions those areas for which ≥8 of the 10 replicate models predicted potential for presence, and a prediction of unsuitable conditions those areas for which ≤4 of the 10 replicate models predicted potential for presence.
Future climate conditions were reconstructed based on two general circulation models (GCMs) that have been used to simulate future climates: those of the Hadley Centre (HadCM3)  and the Canadian Center (CGCM1) . From each GCM, I analyzed two greenhouse gas emissions scenarios: the B2 scenario, which is a relatively conservative estimate of climate change, and the A2 scenario, which is more extreme in the climates reconstructed. As they are based on a 30 yr average around 2055, these models do not take into account potential effects of increased climate variability (El Niño events, in particular) on species' distributions. Because these future climate data are provided at a very coarse native spatial resolution (2.5 × 3.75°), I calculated expected changes in temperature (°C) and precipitation (mm) under each scenario from the relatively coarse raw model results; these expected changes were applied to the original IPCC current climate data layers to provide a final pixel resolution of ~30 × 30 km for future-climate data layers. Models developed based on present-day occurrence patterns and environmental variation were projected to these 4 views of likely future climate conditions.
Finally, I calculated human populations living within the potential distributions of these competent malaria vectors across sub-Saharan Africa based on the LandScan database [29, 30], a 1 km resolution summary of present-day human population distributions globally. It should be noted that I make the explicit assumption of stability of human population distributions – although future population projections are available , they are overly coarse spatially to be greatly informative in this analysis, so this point remains as a future challenge for improvement. For each greenhouse gas emissions scenario, I averaged the results of the two climate models to yield single estimates of future potential distribution under that scenario. Intersecting the model predictions with the Landscan dataset, I calculated present-day numbers of people living in areas coinciding with potential distributional areas of each mosquito species, and present-day numbers of people living in areas converting from unsuitable to suitable or vice versa.
Anopheles gambiae and A. arabiensis, the two most important malaria vectors in Africa, have broad distributions across Africa  (Fig. 1). I do not present detailed validations of model predictions across Africa in the present day because I had published such analyses in an earlier study based on precisely the same occurrence data set  – suffice it to say that model predictions were quite robust, offering confidence in the areas identified as suitable versus unsuitable. Projections of the ecological niche model onto likely future climate conditions indicated that both mosquito species are likely to see less suitable conditions for their populations across portions of West Africa, where temperature increases of 1.5–2.7°C are likely to be manifested (Fig. 2). In contrast, both species are likely to see improving conditions in regions of southern Africa, in areas where annual mean temperatures are increasing sufficiently to permit these species to establish populations (Fig. 2).
Present-day human populations coinciding spatially with the inferred geographic distributions of these vector species presently total 389,155,713 for A. gambiae and 520,289,130 for A. arabiensis over 10.5 × 106 and 14.5 × 106 km2, respectively. Under the two greenhouse gas scenarios, >30 million people are living in areas projected to see increasing exposure to A. gambiae, and >14 million people are living in areas projected to see increased exposure to A. arabiensis. At the same time, however, 78–111 million people are living in areas projected to see reduced exposure to A. gambiae, and 135–171 million people are living in areas projected to see reduced exposure to A. arabiensis (Table 1). These changing patterns are distributed unevenly across the continent: reduced exposure across West Africa and the Sahel, but increased potential for exposure in East Africa and southern Africa (Fig. 1).
This study illustrates the complexities of effects of changing climate conditions on spatial patterns of the potential for disease transmission. Ecological niche models were developed for the two most significant and broadly distributed malaria vectors in Africa and projected onto two greenhouse gas emissions scenarios for each of two global climate models. Clearly, a first priority is to expand these analyses from the two most broadly distributed species in the Anopheles gambiae complex to the broader suite of potential vectors across the continent, as, for instance, A. funestus, A. nili, and A. moucheti are competent vectors in Cameroun . The present analysis was of necessity restricted to the A. gambiae complex for lack of access to high-quality, continentwide occurrence data for other species. A further limitation of this study is its analysis based on present-day human population levels and distributions, for lack of detailed information on future status; similarly, I neglect the effects of urbanization and interventions on mosquito vector distributions and malaria transmission rates.
In the previous global analysis , as in my more detailed analysis, increasing malaria risk was projected in East Africa and decreasing malaria risk was projected for West Africa, so the general pictures presented by the two analyses are parallel, but the detail offered in the present study is considerably greater. Nonetheless, my results require careful examination and exploration, because overall reductions in potential for human exposure nonetheless are achieved by reductions in some regions and increases in others.
Ecological niche models reconstruct patterns of association between species' occurrences and environmental variation across space and time , and offer a predictive understanding of distributional patterns that can extend over changing environmental conditions. They require testing and validation prior to use , however, which in the present case has been the subject of a detailed, published analysis demonstrating excellent predictivity of distributional patterns of the vector species even across broad, unsampled regions . These tests, as well as more general validations of the analytical framework [7, 42–44], increase confidence that projections of niche dimensions across periods of environmental change will also have predictive power. That mosquito species are likely to respond to changing climates is quite probable, particularly in light of recent studies demonstrating that mosquito distributions in time and space respond to the same sets of stimuli .
The climate change picture regarding shifting patterns of malaria vector distributions across Africa is complex. In general, future expectations are that, while fewer people live in areas that will be suitable for malaria vectors across West Africa, more people will be living in exposed areas in East Africa and southern Africa. Because West Africa is more densely populated than parts of eastern and southern Africa, a net decrease in potential for human exposure to these vectors results. For A. gambiae, 78–112 million fewer people are projected to be in suitable areas for the vector species, as compared with 33–34 million people living in areas becoming suitable; for A. arabiensis, reductions in potential for exposure are more dramatic: 136–171 million in areas no longer suitable, but only 14 million in areas newly suitable. The spatial distribution of these shifts, however, is perhaps the most intriguing result: malaria transmission risk is likely to shift eastward and southward, establishing new foci in areas presently lacking intense malaria transmission.
Adding still more complexity are comparisons of the two greenhouse gas emissions scenarios. The A2 scenario emphasizes a heterogeneous world, with regional self-reliance and preservation of local identities, which results in continuously increasing population, regional economic development and slower technological change. The B2 scenario, on the other hand, is one of local solutions to economic, social and environmental sustainability, global population increases at rates lower than A2. Curiously, though, the more extreme (in terms of climate change) A2 scenario is associated with 20–30% less human exposure to malaria across Africa, given its more extreme projected effects on Anopheles distributions, particularly in West Africa. The difference between the two scenarios is thus not so much in terms of how many people live in areas newly suitable for malaria vectors, but rather in terms of how many fewer people live in areas no longer suitable.
Scientific results regarding implications of climate change are usually pessimistic in nature. This analysis departs somewhat from that theme – malaria vector species are expected to undergo distributional shifts in the face of changing climates that will leave fewer people overall in areas suitable for the important vector species. The net reduction nonetheless includes increasing malaria vector presence in areas of eastern and southern Africa that are not presently exposed to these species. As such, the picture is one of complex rearrangement of malaria transmission across Africa with climate change.
Guerra CA, Snow RW, Hay SI: Mapping the global extent of malaria in 2005. Trends in Parasitology. 2006, 22: 353-358. 10.1016/j.pt.2006.06.006.
Beck LR, Rodriguez MH, Dister SW, Rodriguez AD, Washino RK, Roberts DR, Spanner MA: Assessment of a remote sensing-based model for predicting malaria transmission risk in villages of Chiapas, Mexico. Am J Trop Med Hyg. 1997, 56: 99-106.
Omumbo J, Ouma J, Rapuoda B, Craig MH, Lesueur D, Snow RW: Mapping malaria transmission intensity using geographical information systems GIS: an example from Kenya. Ann Trop Med Parasitol. 1998, 92: 7-21.
Gemperli A, Vounatsou P, Sogoba N, Smith T: Malaria mapping using transmission models: Application to survey data from Mali. American Journal of Epidemiology. 2006, 163: 289-297. 10.1093/aje/kwj026.
Gemperli A, Sogoba N, Fondjo E, Mabaso M, Bagayoko M, Briet OJT, Anderegg D, Liebe J, Smith T, Vounatsou P: Mapping malaria transmission in West and Central Africa. Tropical Medicine & International Health. 2006, 11: 1032-1046. 10.1111/j.1365-3156.2006.01640.x.
Elith J, Graham CH, Anderson RP, Dudik M, Ferrier S, Guisan A, Hijmans RJ, Huettman F, Leathwick JR, Lehmann A, et al: Novel methods improve prediction of species' distributions from occurrence data. Ecography. 2006, 29: 129-151. 10.1111/j.2006.0906-7590.04596.x.
Araújo MB, Pearson RG, Thuiller W, Erhard M: Validation of species-climate impact models under climate change. Glob Change Biol. 2005, 11: 1504-1513. 10.1111/j.1365-2486.2005.01000.x.
Pearson RG, Dawson TP: Predicting the impacts of climate change on the distribution of species: Are bioclimate envelope models useful?. Global Ecology and Biogeography. 2003, 12 (5): 361-371. 10.1046/j.1466-822X.2003.00042.x.
Peterson AT: Biogeography of diseases: A framework for analysis. Naturwissenschaften. 2008, 95: 483-491. 10.1007/s00114-008-0352-5.
Martens WJM, Niessen LW, Rotmans J: Potential impact of global climate change on malaria risk. Environmental Health Perspectives. 1995, 103: 458-464. 10.2307/3432584.
Martens P, Kovats RS, Nijhof S, Vries Pd, Livermore MTJ, Bradley DJ, Cox J, McMichael AJ: Climate change and future populations at risk of malaria. Global Environmental Change. 1999, 9: S89-S107. 10.1016/S0959-3780(99)00020-5.
Tanser F, Sharp B, Sueur Dl: Potential effect of climate change on malaria transmission in Africa. Lancet. 2003, 362: 1792-1798. 10.1016/S0140-6736(03)14898-2.
Rogers DJ, Randolph SE: The global spread of malaria in a future, warmer world. Science. 2000, 289: 1763-1766. 10.1126/science.289.5478.391b.
Lindblade KA, Walker ED, Wilson ML: Early warning of malaria epidemics in African highlands using Anopheles (Diptera: Culicidae) indoor resting density. Journal of Medical Entomology. 2000, 37: 664-674.
Hay SI, Cox J, Rogers DJ, Randolph SE, Stern DI, Shanks GD, Myers MF, Snow RW: Climate change and the resurgence of malaria in the East Africa highlands. Nature. 2002, 415: 905-909. 10.1038/415905a.
Rogers DJ, Randolph SE, Snow RW, Hay SI: Satellite imagery in the study and forecast of malaria. Nature. 2002, 415: 710-715. 10.1038/415710a.
Moffett A, Shackelford N, Sarkar S: Malaria in Africa: Vector species' niche models and relative risk maps. PLoS ONE. 2007, 2: e824-10.1371/journal.pone.0000824.
Lieshout Mv, Kovats RS, Livermore MTJ, Martens P: Climate change and malaria: analysis of the SRES climate and socio-economic scenarios. Global Environmental Change. 2004, 14: 87-99. 10.1016/j.gloenvcha.2003.10.009.
Grinnell J: Field tests of theories concerning distributional control. American Naturalist. 1917, 51: 115-128. 10.1086/279591.
Wiens JJ, Graham CH: Niche conservatism: Integrating evolution, ecology, and conservation biology. Annual Review of Ecology, Evolution and Systematics. 2005, 36: 519-539. 10.1146/annurev.ecolsys.36.102803.095431.
Peterson AT, Soberón J, Sánchez-Cordero V: Conservatism of ecological niches in evolutionary time. Science. 1999, 285: 1265-1267. 10.1126/science.285.5431.1265.
Araujo MB, Pearson RG: Equilibrium of species' distributions with climate. Ecography. 2005, 28 (5): 693-695. 10.1111/j.2005.0906-7590.04253.x.
Cooper RD, Waterson DGE, Frances SP, Beebe NW, Sweeney AW: Speciation and distribution of the members of the Anopheles punctulatus (Diptera: Culicidae) group in Papua New Guinea. Journal of Medical Entomology. 2002, 39 (1): 16-27.
Sweeney AW, Beebe NW, Cooper RD, Bauer JT, Peterson AT: Environmental factors associated with the distribution and range limits of the malaria vector Anopheles farauti sensu stricto in Australia. Journal of Medical Entomology. 2006, 43: 1068-1075. 10.1603/0022-2585(2006)43[1068:EFAWDA]2.0.CO;2.
Benedict MQ, Levine RS, Hawley WA, Lounibos LP: Spread of the tiger: Global risk of invasion by the mosquito Aedes albopictus. Vector-Borne and Zoonotic Diseases. 2007, 7: 76-85. 10.1089/vbz.2006.0562.
Peterson AT, Martínez-Campos C, Nakazawa Y, Martínez-Meyer E: Time-specific ecological niche modeling predicts spatial dynamics of vector insects and human dengue cases. Transactions of the Royal Society of Tropical Medicine and Hygiene. 2005, 99: 647-655. 10.1016/j.trstmh.2005.02.004.
Levine RS, Benedict MQ, Peterson AT: Distribution of Anopheles quadrimaculatus Say s.l. and implications for its role in malaria transmission in the US. Journal of Medical Entomology. 2004, 41: 607-613.
Levine RS, Peterson AT, Benedict MQ: Geographic and ecologic distributions of the Anopheles gambiae complex predicted using a genetic algorithm. Am J Trop Med Hyg. 2004, 70: 105-109.
Dobson JE, Bright EA, Coleman PR, Durfee RC, Worley BA: LandScan – A global population database for estimating populations at risk. Photogrammetric Engineering & Remote Sensing. 2000, 66: 849-857.
LandScan: LandScanTM Global Population Database. 2006, Oak Ridge, Tennessee: Oak Ridge National Laboratory, [http://www.ornl.gov/landscan/]
Coetzee M, Craig M, le Sueur D: Distribution of African malaria mosquitoes belonging to the Anopheles gambiae complex. Parasitology Today. 2000, 16: 74-77. 10.1016/S0169-4758(99)01563-X.
New M, Hulme M, Jones P: A 1961–1990 Mean Monthly Climatology of Global Land Areas. 1997, Norwich, U.K.: Climatic Research Unit, University of East Anglia
USGS: HYDRO1k Elevation Derivative Database. 2001, Washington, D.C.: U.S. Geological Survey, [http://edc.usgs.gov/products/elevation/gtopo30/hydro/index.html]
Stockwell DRB, Peters DP: The GARP modelling system: Problems and solutions to automated spatial prediction. International Journal of Geographical Information Science. 1999, 13: 143-158. 10.1080/136588199241391.
Anderson RP, Lew D, Peterson AT: Evaluating predictive models of species' distributions: Criteria for selecting optimal models. Ecological Modelling. 2003, 162: 211-232. 10.1016/S0304-3800(02)00349-6.
Pope VD, Gallani ML, Rowntree VJ, Stratton RA: The impact of new physical parametrizations in the Hadley Centre climate model – HadAM3. 2002, Bracknell, Berks, UK: Hadley Centre for Climate Prediction and Research
Flato GM, Boer GJ, Lee WG, McFarlane NA, Ramsden D, Reader MC, Weaver AJ: The Canadian Center for Climate Modelling and Analysis Global Coupled Model and its climate. Climate Dynamics. 1999, 16: 451-467. 10.1007/s003820050339.
Arnell NW, Livermore MJL, Kovats S, Levy PE, Nicholls R, Parry ML, Gaffin SR: Climate and socio-economic scenarios for global-scale climate change impacts assessments: Characterising the SRES storylines. Global Environmental Change. 2004, 14: 3-20. 10.1016/j.gloenvcha.2003.10.004.
Bigoga J, Manga L, Titanji V, Coetzee M, Leke R: Malaria vectors and transmission dynamics in coastal south-western Cameroon. Malaria Journal. 2007, 6: 5-10.1186/1475-2875-6-5.
Peterson AT: Uses and requirements of ecological niche models and related distributional models. Biodiversity Informatics. 2006, 3: 59-72.
Peterson AT: Kansas Gap Analysis: The importance of validating distributional models before using them. Southwestern Naturalist. 2005, 50: 230-236. 10.1894/0038-4909(2005)050[0230:KGATIO]2.0.CO;2.
Pearson RG, Raxworthy C, Nakamura M, Peterson AT: Predicting species' distributions from small numbers of occurrence records: A test case using cryptic geckos in Madagascar. Journal of Biogeography. 2007, 34: 102-117. 10.1111/j.1365-2699.2006.01594.x.
Martínez-Meyer E, Peterson AT: Conservatism of ecological niche characteristics in North American plant species over the Pleistocene-to-Recent transition. Journal of Biogeography. 2006, 33: 1779-1789. 10.1111/j.1365-2699.2006.01612.x.
Martínez-Meyer E, Peterson AT, Hargrove WW: Ecological niches as stable distributional constraints on mammal species, with implications for Pleistocene extinctions and climate change projections for biodiversity. Global Ecology and Biogeography. 2004, 13: 305-314. 10.1111/j.1466-822X.2004.00107.x.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2334/9/59/prepub
I thank Mark Benedict for his interest and assistance with securing key data resources.
The author declares that they have no competing interests.
ATP developed entire manuscript
About this article
Cite this article
Peterson, A.T. Shifting suitability for malaria vectors across Africa with warming climates. BMC Infect Dis 9, 59 (2009). https://doi.org/10.1186/1471-2334-9-59