Estimating the malaria risk of African mosquito movement by air travel
© Tatem et al; licensee BioMed Central Ltd. 2006
Received: 15 May 2006
Accepted: 14 July 2006
Published: 14 July 2006
The expansion of global travel has resulted in the importation of African Anopheles mosquitoes, giving rise to cases of local malaria transmission. Here, cases of 'airport malaria' are used to quantify, using a combination of global climate and air traffic volume, where and when are the greatest risks of a Plasmodium falciparum-carrying mosquito being importated by air. This prioritises areas at risk of further airport malaria and possible importation or reemergence of the disease.
Monthly data on climate at the World's major airports were combined with air traffic information and African malaria seasonality maps to identify, month-by-month, those existing and future air routes at greatest risk of African malaria-carrying mosquito importation and temporary establishment.
The location and timing of recorded airport malaria cases proved predictable using a combination of climate and air traffic data. Extending the analysis beyond the current air network architecture enabled identification of the airports and months with greatest climatic similarity to P. falciparum endemic regions of Africa within their principal transmission seasons, and therefore at risk should new aviation routes become operational.
With the growth of long haul air travel from Africa, the identification of the seasonality and routes of mosquito importation is important in guiding effective aircraft disinsection and vector control. The recent and continued addition of air routes from Africa to more climatically similar regions than Europe will increase movement risks. The approach outlined here is capable of identifying when and where these risks are greatest.
Throughout history the opening of travel and trade routes between countries has been accompanied by the spread of diseases and their vectors [1, 2]. Air travel has been identified as a prime factor in the global spread of both infectious and vector-borne diseases [3, 4] and represents a longstanding concern . Many recent disease vector invasions are suspected to have resulted from air travel [6–8] and, with continual rapid expansion in global air travel, the threat of future invasions should not be ignored. The public health and economic impacts of past disease vector invasions  are illustrated by Aedes aegypti's invasion of the Americas  and the escape of Anopheles gambiae to Brazil from Africa .
The last 30 years has seen air travel to tropical regions of the World rise dramatically. International tourist arrivals to sub-Saharan Africa (SSA) increased from 6.7 million in 1990 to over 17 million in 2000 . Mosquito species, including Aedes albopictus and An. gambiae s.l., have been shown to survive long haul flights [11–13]. For example, in one 3-week period in 1994, it was estimated that 2000–5000 Anopheline mosquitoes were imported into France at a rate of 8–20 Anopheline mosquitoes per flight . An example of the effects of such importations are the many cases of autochthonous (locally-acquired) malaria, which are principally clustered around international airports: so called "airport malaria" [14, 15]. These occur primarily through the transport of infected Anopheles mosquitoes that can survive for long enough after arrival to transmit malaria . Recent studies suggest that, where used, routine disinsection is proving effective in reducing airport malaria risk , although the number of countries implementing such procedures is in decline [17, 18]
Airport malaria cases are rare, with just two cases per year on average recorded (1969–99), all but one of which were P. falciparum malaria. However, these cases provide important evidence of sufficient traffic volumes and climatic similarity between origin and destination for the survival of malaria-carrying Anopheles mosquitoes. Here we describe an exploratory approach which makes use of this information derived from airport malaria cases to quantify, in terms of global climate and air traffic, which airports have the greatest risks of local P. falciparum malaria transmission through importation from sub-Saharan Africa of infected mosquitoes. We estimate a) the risks based on year 2000 air traffic volumes, b) how this varies throughout the year and c) where the greatest potential future risks would lie through the opening of new routes.
Flight data on total passenger numbers in the year 2000 moving between the World's top 100 airports by traffic (100% aircraft capacity was assumed), were obtained from OAG Worldwide Ltd . The database includes cargo flights. For full geographical coverage, the database also included data on the principal airports of 143 other countries not represented in this top 100. Data on a total of 7129 routes between 278 international airports in the year 2000 were thus available. To provide an estimate of the most likely months of malaria movement, information on the onset and end of the principal malaria transmission seasons for each African country were obtained from maps of malaria seasonality across Africa . These maps were derived only for P. falciparum transmission by An. gambiae s.l., and all results presented here refer only to this parasite and vector combination. The results may be extended to other species, but this paper concerns itself with the most efficient vector of malaria (An. gambiae) and the most pathogenic of the malaria parasites (P. falciparum) [21, 22].
A 10 × 10 minute (~18 × 18 km at the equator) spatial resolution gridded climatology was used to extract mean temperature, rainfall and humidity data for a synoptic year (1961–1990) , and these layers were then linearly rescaled to a common range of values by dividing through by the maximum in each layer. Climatic similarity between origin and destination airports was assumed to determine principally whether mosquitoes from originating airports would survive and have the potential to produce local malaria transmission around the destination airports. The locations of the 278 airports were superimposed onto the monthly climatology surfaces and each 10 × 10 minute spatial resolution grid square covering the airport location identified. To ensure a representative climate was included, the eight land grid squares surrounding each airport square were also identified, forming a three by three grid square centered on the airport. This was not possible for coastal airports or those on small islands, where reduced numbers of land grid squares were extracted. Airports located on islands too small to be represented by the climatology surfaces were eliminated from the analysis, reducing the sample size to 259 airports. Ideally, data for these airports would have been included in the analysis, but the need to use a global climatology meant this was not possible. Data from the grid squares identified in each of the three climatology surfaces thus formed climate 'signatures' for each month, for each airport.
Distance measures, clustering and dendrograms
The lack of sufficient variance in the climate data constituting the majority of signatures, and the location of many airports on small islands, dictated that only simple Euclidean distance could be used as a measure of climatic similarities between origin and destination airports . Euclidean distance is defined as the shortest straight line distance between two points, in this case, the distance between the environmental signature centroids in three-dimensional climatic space, as defined by the temperature, rainfall and humidity values. Euclidean distances between each signature centroid and the centroids of every other signature were calculated to derive separate monthly climate "dissimilarity" matrices. The climate dissimilarity matrices were then subject to hierarchical clustering using an agglomerative algorithm. The clustering results were then translated into dendrograms based on centroid linkage using Phylip v3.63 (University of Washington, USA), for each month of the year. The dendrograms are monthly climate-based phenetic trees, and represent the global air transport network remapped in terms of disease vector suitability [2, 24]. Tests using those signatures which did facilitate the use of more sophisticated distance measures (e.g. divergence, Jefferies-Matusita), revealed that the resultant dendrogram architecture was very similar to those developed using Euclidean distances (results not shown).
Climatic similarity thresholds
To define how similar airport climates need to be to permit the temporary survival of imported Anopheles and the possible transmission of P. falciparum malaria, confirmed examples of where this has occurred were used. These cases of airport malaria confirm that, at the time of year of the case, the climates in the vicinity of the origin and destination airports were similar enough for the survival of malaria-carrying Anopheles, and transmission by them at the destination. In no suspected cases of airport malaria has the origin of imported Anopheles or malaria been confirmed unambiguously, so it was assumed that the origin was the SSA airport or region most climatically similar to the destination and within its primary transmission season. This gives the most conservative estimate of the climatic range within which imported malaria transmission is possible.
The African airports in the database were overlaid on the malaria seasonality maps and classified as either suitable or not for transmission during each month of the year. Those non-malarious airports linked via a dendrogram branch lower than the airport malaria-defined climatic limit to a malarious airport during its transmission season were identified as being sufficiently similar climatically, based on previous airport malaria cases, for there to be a risk of infectious mosquito arrival, survival and possible consequent local malaria transmission.
To obtain a monthly measure of malarial mosquito movement risk to those airports identified as being at risk within the dendrogram in 2000, the monthly Euclidean environmental distances and traffic data were used. The Euclidean distances between malarious airports and other airports were rescaled linearly by dividing through by the maximum Euclidean distance in the climate dissimilarity matrix, and inverted to lie between zero and one. Therefore, the most climatically similar airport pairs had a value close to one and those distinctly different, a value close to zero. Similarly, the air traffic data were also rescaled, resulting in the airports with the most direct traffic running between them in a year taking a value close to one, while those direct routes with little or no traffic a value close to or at zero. It was assumed for this analysis that passenger numbers on each route remained consistent year-round. For each month and route between malarious and other airports, the rescaled climatic and traffic matrices were then multiplied together to provide a simple monthly measure of the likelihood of temporary malarious mosquito importation and possible consequent local transmission.
Sensitivity of results
The results obtained are dependent upon both the accuracy of the seasonality map and the choice of dendrogram threshold. It was therefore essential to assess how sensitive these results were to small alterations in each. The approach, results and related discussion are provided in Additional File 1, and demonstrate that the results presented remain relatively insensitive to small changes in the seasonality map and choice of dendrogram threshold.
Year 2000 situation
Year 2000 top 10 air travel risk routes for P. falciparum infected An. gambiae invasion and subsequent autochthonous transmission.
Annual No. Passengers
Paris Charles de Gaulle
Top 20 climatically closest linked destination airports with malarious-SSA airports within principal transmission season per-month.
Dar Es Salaam
St Vincent and the Grenadines
Dar Es Salaam
Dar Es Salaam
Year 2000 situation
The routes identified in table 1 as the most likely for temporary malarious mosquito invasion and possible consequent autochthonous transmission in 2000 reflect both the structure of the global air traffic network and sufficient climatic similarities between SSA and Europe in July, August and September. There is strong correspondence between the risk routes and suspected airport malaria case location (Figure 2) and timing (Figure 3). The predominance of West African airports in table 1 occurs because the malaria transmission season timing coincides with the European summer. Although many other airports around the World were identified as more similar climatically to malarious SSA airports than the airport malaria examples, the majority of flights originating from SSA have European destinations, with over 100,000 per annum travelling on certain routes (e.g. Abidjan to Paris, see Table 1). Figure 4 shows that European destinations may be sufficiently similar climatically to SSA airports for malarious mosquito survival only 2–4 months a year, yet because almost all air traffic is directed toward Europe, the continent receives the majority of imported P. falciparum cases  and almost all recorded airport malaria episodes (Figure 1) . Such a short window of sufficient climatic similarity suggests that increased vigilance and disinsection in terms of malaria-carrying Anopheles importation by incoming flights from SSA may only be required during the summer months in Europe, especially during unusually hot and humid periods. Disinsection, however, is also undertaken for arbovirus-carrying mosquito control. This paper has not addressed this risk, but other work suggests shipping to be most important in their spread [7, 24].
The likelihood of P. falciparum or An. gambiae establishment through air travel in Europe and many other locations remains low. Unsuitable year-round climate, enforced disinsection measures, An. gambiaes' intolerance of urban areas  and competition from local mosquitoes that are inefficient vectors of P. falciparum all provide barriers to establishment. However, that all destinations identified for 2000 have experienced local malaria transmission within the months shown, is evidence of the existing risks and the need to incorporate climatic information in predictions of future risks. Although a conservative climatic similarity threshold was chosen, the fact that only 19 risk routes (of a possible 1,032) were identified for 2000, shows that the current air traffic network is restricted climatically, with almost all air traffic from SSA directed to climatically dissimilar European airports. This is perhaps one of the reasons why An. gambiae and P. falciparum malaria have escaped so few times from Africa.
Evidence of the importance of shipping traffic volume, particularly container shipping, in vector invasion has been outlined elsewhere [2, 24]. While the focus of this analysis has been on passenger flights, this includes thousands of cargo-laden flights, and in terms of cargo-only flights, database interrogation shows that these make up less than 7% of all scheduled flights. The role of passenger air traffic volume is emphasized here through comparison of Figures 1 and 2. The spatial pattern of airport malaria cases is reflected in the annual incoming passenger numbers, with the adjusted r2 = 0.704 (n = 19, p < 0.01) between the two. The fact that France has received more suspected cases of airport malaria than the UK, despite similar incoming passenger volumes from endemic airports, can be attributed to more routes linking it to West Africa, where transmission season and European summer timings match, though other factors, such as differences in surveillance systems, may play a part. Given that air traffic volumes generally peak in the northern hemisphere summer months, it should be additionally noted that the year round constant traffic volumes assumed here could produce an underestimation of risk levels.
The results show that many airports in regions outside of Europe are climatically more favorable than is Europe for Anopheles survival, and for many more months of the year. The effects of opening up new routes from malaria-endemic SSA countries to non-European destinations could therefore have potentially serious consequences. The accidental introduction of An. gambiae into Brazil in 1930, resulted in over 16,000 malaria deaths and a mosquito control program that cost 3 billion USD equivalent today [9, 27]. Recently, new routes have been opening, with long-haul flights from SSA direct to Washington, Beijing, Hong Kong and Bangkok amongst others. With more routes planned, and larger planes capable of traveling greater distances in shorter times scheduled, the risk of another An. gambiae escape grows . While recent evidence points to the effectiveness of routine aircraft cabin disinsection in flights to the UK from SSA , elsewhere in the World disinsection is in decline. The World Health Organization continues to recommend aircraft disinsection , but fears over the health effects of the pyrethroid insecticides used  and resultant law suits  have led to many airlines and governments ceasing the practice [17, 18].
The destinations in Table 2 with climates matched almost perfectly for the month in question to SSA airports home to P. falciparum infected Anopheles are all in malarious regions, or those declared malaria-free only recently [32, 33]. Should regular flights commence on any of those routes identified in Table 2, the potential for temporary An. gambiae invasion and consequent P. falciparum transmission in those months may be considerable. Though many in Table 2 appear to be unlikely routes in the near future, just outside the top 20 are routes of potential future operation or routes which have recently opened (e.g. Addis Ababa to Washington, Dakar to New York). A route for almost every month of the year is identified in the top 20 alone, highlighting that should SSA become better connected by air, a year-round linkage with climatically similar airports would exist, potentially facilitating the escape of An. gambiae from SSA. Moreover, the opening of flight routes from south-east Asia to SSA may also speed up the global movement of drug resistant malaria (in humans or mosquitoes) and insecticide resistant mosquitoes [34, 35].
Figure 4 shows, as expected, that the climatic regimes of the airports in the tropical regions of the Americas and south-east Asia are more similar to those of malaria-endemic SSA airports than are the latter to those of most European airports. This is especially true of the airports in Central America and the Caribbean which have the closest match of all airports globally, and for the longest period. This explains the predominance of Caribbean and Central American destination airports in Table 2.
Unsurprisingly the periods of greatest similarity of higher latitude climates to those of SSA airports are the summer months of May to September in the northern hemisphere and November to February in the southern hemisphere (Figure 5); it is, therefore, during these months that airport malaria in these regions is most likely. The equivalent periods of similarity for tropical latitudes are much more variable, with little apparent spatial coherency in peak month. Evidence of this can be seen in Table 2, which includes two routes to Miami, USA, one from Dar Es Salaam in East Africa in May, the other from Abidjan in West Africa in October. In the past, high risk locations and months for temporary Anopheles invasion leading to local transmission have been predictably along specific routes to certain European airports, and solely in the summer months. Figure 5 shows that should routes between P. falciparum-endemic SSA countries and other tropical locations continue to open up, the risks, origin and timing of Anopheles invasion will be much less predictable.
The continued occurrence of airport malaria cases and rise in imported P. falciparum malaria cases [15, 25] are indicative of an expanding global air transport network and increased travel to malarious countries . It also reflects a decline in aircraft disinsection [17, 18]. The approach presented highlights routes and months within the current global air transport network at risk of importation and temporary establishment of P. falciparum-carrying mosquitoes. The development of new air travel routes from SSA suggests that the relative risks identified will continue to increase and that monitoring schemes based on climate suitability methods could be used to optimize Anopheles-specific disinsection and control efforts. The analysis presented here is a first step, and future challenges to refine the predictions will involve incorporating information on temporal variations in passenger numbers, flight stopover risks, intra-species competition, human populations at risk, breeding site availability, possible climate change , disinsection and land transport, and quantifying the relative importance of sea and air transport for vectors and diseases.
We thank Abdisalan Noor, Sarah Randolph, Robert Snow and Briony Tatem for comments on earlier drafts of the manuscript. AJT and SIH are funded by a Research Career Development Fellowship (to SIH) from the Wellcome Trust (#069045).
- McNeill WH: Plagues and people. 1976, New York , Anchor PressGoogle Scholar
- Tatem AJ, Rogers DJ, Hay SI: Global transport networks and infectious disease spread. Advances in Parasitology. 2006, 62: 293-343.PubMed CentralView ArticlePubMedGoogle Scholar
- Gezairy HA: Travel epidemiology: WHO perspective. International Journal of Antimicrobial Agents. 2003, 21: 86-88. 10.1016/S0924-8579(02)00365-5.View ArticlePubMedGoogle Scholar
- Grais RF, Ellis JH, Glass GE: Assessing the impact of airline travel on the geographic spread of pandemic influenza. European Journal of Epidemiology. 2003, 18: 1065-1072. 10.1023/A:1026140019146.View ArticlePubMedGoogle Scholar
- Massey A: Epidemiology in relation to air travel. 1933, London , H. K. Lewis and Co. LimitedGoogle Scholar
- Gratz NG, Steffen R, Cocksedge W: Why aircraft disinsection?. Bull World Health Organ. 2000, 78 (8): 995-1004.PubMed CentralPubMedGoogle Scholar
- Lounibos LP: Invasions by insect vectors of human disease. Annual Review of Entomology. 2002, 47: 233-266. 10.1146/annurev.ento.47.091201.145206.View ArticlePubMedGoogle Scholar
- Mangili A, Gendreau MA: Transmission of infectious diseases during commercial air travel. Lancet. 2005, 365: 989-996. 10.1016/S0140-6736(05)71089-8.View ArticlePubMedGoogle Scholar
- Soper FL, Wilson DB: Anopheles gambiae in Brazil: 1930 to 1940. 1943, New York , Rockefeller FoundationGoogle Scholar
- Franco-Parades C, Santos-Preciado JI: Problem pathogens: prevention of malaria in travellers. Lancet Infectious Diseases. 2006, 6: 139-149. 10.1016/S1473-3099(06)70410-8.View ArticleGoogle Scholar
- Misao T, Ishihara M: An experiment on the transportation of vector mosquitoes by aircraft (in japanese). Rinsho to Kenkyu. 1945, 22: 44-46.Google Scholar
- Russell RC: Survival of insects in the wheelbays of Boeing 747B aircraft on flights between tropical and temperate airports. Bulletin of the World Health Organisation. 1987, 65: 659-662.Google Scholar
- Russell RC: Transport of insects of public health importance on international aircraft. Travel Medicine International. 1989, 26-31.Google Scholar
- Isaäcson M: Airport malaria: a review. Bulletin of the World Health Organization. 1989, 67 (6): 737-743.PubMed CentralPubMedGoogle Scholar
- Muentener P, Schlagenhauf P, Steffen R: Imported malaria (1985-95): trends and perspectives. Bulletin of the World Health Organization. 1999, 77 (7): 560-566.PubMed CentralPubMedGoogle Scholar
- Hutchinson R, Bayoh M, Lindsay S: Risk of airport malaria in the UK. European Mosquito Bulletin. 2005, 19: 12-13.Google Scholar
- The plane truth about disinsection. Environmental Health Perspectives. 1999, 107 (8): 397-398.
- Russell RC, Paton R: In-flight disinsection as an efficacious procedure for preventing international transport of insects of public health importance. Bulletin of the World Health Organization. 1989, 67 (5): 543-547.PubMed CentralPubMedGoogle Scholar
- OAG Worldwide Ltd. [http://www.oag.com]
- Tanser FC, Sharp B, le Sueur D: Potential effect of climate change on malaria transmission in Africa. The Lancet. 2003, 362 (9398): 1792-1798. 10.1016/S0140-6736(03)14898-2.View ArticleGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- Snow RW, Guerra CA, Noor AM, Myint HL, Hay SI: The global distribution of clinical episodes of Plasmodium falciparum malaria. Nature. 2005, 434: 214-217. 10.1038/nature03342.PubMed CentralView ArticlePubMedGoogle Scholar
- New M, Lister D, Hulme M, Makin I: A high-resolution data set of surface climate over global land areas. Climate Research. 2002, 21: 1-25.View ArticleGoogle Scholar
- Tatem AJ, Hay SI, Rogers DJ: Global traffic and disease vector dispersal. Proceedings of the National Academy of Sciences. 2006, 103: 6242-6247. 10.1073/pnas.0508391103.View ArticleGoogle Scholar
- Toovey S, Jamieson A: Rolling back malaria: how well is Europe doing?. Travel Medicine and Infectious Diseases. 2003, 1: 167-175. 10.1016/j.tmaid.2003.09.004.View ArticleGoogle Scholar
- Hay SI, Guerra CA, Tatem AJ, Atkinson PM, Snow RW: Urbanization, malaria transmission and disease burden in Africa. Nature Reviews Microbiology. 2005, 3 (1): 81-90. 10.1038/nrmicro1069.PubMed CentralView ArticlePubMedGoogle Scholar
- Killeen GF, Fillinger U, Kiche I, Gouagna LC, Knols BGJ: Eradication of Anopheles gambiae from Brazil: lessons for malaria control in Africa?. The Lancet Infectious Diseases. 2002, 2: 618-627. 10.1016/S1473-3099(02)00397-3.View ArticlePubMedGoogle Scholar
- Powell JR, Coluzzi M: Malarial mosquito: Is Anopheles gambiae plotting an escape. International Herald Tribune. 2004Google Scholar
- World Health Organisation: WHO recommendations on the disinsecting of aircraft. Weekly Epidemiological Record. 1985, 60 (7): 45-47.Google Scholar
- Das R, Cone J, Sutton P: Aircraft disinsection. Bulletin of the World Health Organisation. 2001, 79 (9): 900-901.Google Scholar
- Woodyard C: Fliers fume over planes treated with pesticides. USA Today. 2001Google Scholar
- Guerra CA, Snow RW, Hay SI: Defining the global spatial limits of malaria transmission in 2005. Advances in Parasitology. 2006, 62: 157-179.PubMed CentralView ArticlePubMedGoogle Scholar
- Guerra CA, Snow RW, Hay SI: Mapping the global extent of malaria in 2005. Trends in Parasitology. 2006, in press:Google Scholar
- Okeke I: Stopping the spread of drug-resistant malaria. Science. 2004, 306: 2039-10.1126/science.306.5704.2039c.View ArticlePubMedGoogle Scholar
- Roper C, Pearce R, Nair S, Sharp B, Nosten F, Anderson T: Intercontinental spread of Pyrimethamine-resistant malaria. Science. 2004, 305 (5687): 1124-10.1126/science.1098876.View ArticlePubMedGoogle Scholar
- Loutan L: Malaria: still a threat to travellers. International Journal of Antimicrobial Agents. 2003, 21: 158-163. 10.1016/S0924-8579(02)00367-9.View ArticlePubMedGoogle Scholar
- Tatem AJ, Rogers DJ, Hay SI: T8.5: Traffic in disease vectors: 2005, 2015 and 2030. Commissioned as part of the U.K. Government's Foresight project, Infectious Diseases: preparing for the future. 2006, London, U.K. , Office of Science and InnovationGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.