Electronic supplementary material

RCM mean bias correction

Rainfall and temperature have been mean-bias corrected for each individual RCM with respect to the 1990-2009climatology. Namely for each year:

For temperature:

For rainfall:

i and j denote a specific grid point per latitude / longitude, and m denotes a specific month. The correction is applied for each month within the season. The reference climatology is based on the 1990-2009 time period.

Overwintering criterion

Figure S1a shows the suitable regions for A. albopictus to survive European winters based on mean observed climatic variables for the 1960-1989 period. The western and southern coasts of England, Wales, western France, part of Belgium and Netherlands, the north-western coasts of Spain (Galicia region), Portugal, northern and the south-western coasts of Italy, the western coasts of the Balkans and Greece appears to be suitable areas for the overwintering of A. albopictus (for Rain_ann>700mm and T_jan > 2°C). The delineation of overwintering areas with T_jan > 0°C and Rain_ann > 500m has been generally accepted in former published studies (Burgess, 1995; Mitchell, 1995; Knudsen et al., 1996; Rodhain, 1996; Fisher et al., 2011). The regions that do not fit this criterion (unsuitable) are depicted by the light grey shading on Figure S1. Considering these standard climatic thresholds, most of western Europe including France, northern and the eastern coast of Spain, Portugal, most of the UK excepting the highland regions of Scotland, the Benelux, most of Italy and the coast of the Adriatic appear to be possible overwintering areas for A. albopictus over the period 1960-1989 while most of eastern Europe and the Scandinavian countries are of low or no suitability as the winters are too extreme or the conditions are too dry for the mosquito to survive. Within the context of climate change, annual rainfall has generally increased over northern Europe (more marked during the winter season) whereas this has decreased over southern Europe and the Mediterranean basin especially (more pronounced in summer) during the last two decades. A significant decrease in the number of frost days in winter has also been observed. Figure S1b shows the regions for the overwintering of A. albopictus over the last 20 years (1990-2009), for which temperature and the rainfall amounts have significantly changed.Aedes albopictus suitable areas (white shading) have spread across the Benelux, the western coasts of Denmark, and western Germany. In Spain, conditions have become more unsuitable for the mosquito. Considering the standard overwintering criterion (light gray shading), it is notable that unsuitable areas have moved eastward over north-western central Europe (Germany), and that part of Croatia, Slovenia and Bosnia-Herzegovina have become more suitable for the overwintering mosquito during the last two decades. As the vector has been introduced recently into these countries through importation, recent climate change conditions might have favoured its survival in the winter conditions.

The future scenario based on regional climate model projections highlight similar overwintering trends for the period 2030-2050 (FigS2). Regions become more suitable over north-western central Europe (Benelux, western Germany and western Denmark), as the unsuitable zones retreat eastward with respect to the recent climatic context. This is related to both an increase in simulated temperatures and increased moister conditions over northern Europe (especially in winter, not shown). Future climate projections highlight that A. albopictus overwintering is also favoured over the Balkans. A larger area of Spain becomes more unsuitable as simulated drier conditions (occurring especially in summer) could limit the southward expansion of the vector. There are no significant future changes over the UK and France with respect to the recent climatic context (1990-2009). However, these results are based on the multi-model ensemble mean. There are significant differences across the different model projections (FigS3). Some RCMs projections are close to the ensemble mean pattern (SMHIRCA, KNMI-RACMO2 models). Some models (C4IRCA4, ETHZ-CLM) simulate a large eastward retreat of the unsuitable areas over north-western and central Europe (including most of Germany and Denmark), while for others (MPI-M-REMO, OURANOSMRRCC) the simulated future overwintering pattern is relatively close to the one which is based on recent climate observations.

Multimedia

An animated movie depicting the climate suitability of A.albopictusbased on Model 1(Kobayashi et al., 2002) has been createdby the authors and is available on youtube at this address:

The present period (1950-2009) is based on gridded observation data (EOBS), while the future projections (2010-2050) are based on the regional climate model ensemble mean carried out within the ENSEMBLES project (SRESA1B) scenario. Trends are depicted as a 5 year running average is applied to the data. The RCMs outputs are mean bias corrected with respect to the observations.

Figure S1: Overwintering capacity of A. albopictus to survive based on mean observed January temperatures and annual rainfall. This is carried out for the period 1960-1989 (a) and 1990-2009 (b) based on the EOBS dataset. White areas depict regions for which climate is suitable for the mosquito. The different grey and black shadings depict regions of unsuitability according to different climatic thresholds. The light gray shading depicts strictly unsuitable areas.

Figure S2: Futureoverwintering capacity of A. albopictus to survive based on mean observed January temperatures and annual rainfall. This is carried out for the period 2030-2050 based on the SimA1B scenario experiment (ensemble mean of all RCM experiments). White areas depicts regions for which climate is suitable for the mosquito. The different grey and black shadings depict regions of unsuitability according to different climatic thresholds. The light gray shading depicts strict unsuitability areas (lower scenario).

Figure S3:Overwintering capacity of A. albopictus to survive based on mean simulated January temperatures and annual rainfall. This is carried out for the period 2030-2050 for each Regional Climate Model projections. The RCMs data are mean bias corrected with the mean reference climate based on the EOBS data for the period 1990-2009 (the delta estimated from the models is added to the mean observed climate).
Figure S4: Climatic parameters employed to drive the different methods. Observed annual rainfall climatology for the period a) 1960-89 and b) 1990-2009. Mean observed January temperatures for the period c) 1960-89 and d) 1990-2009. Observed annual temperature climatology for the period e) 1960-89 and f) 1990-2009. Mean observed summer (June-July-August) temperatures for the period g) 1960-89 and h) 1990-2009.
Figure S5:Aedes albopictussuitability based on the standard overwintering criterion and annual mean temperatures (Model 1). This is carried out for the period 2030-2050 for each Regional Climate Model projections. The RCMs data are mean bias corrected with the mean reference climate based on the EOBS data for the period 1990-2009 (the delta estimated from the models is added to the mean observed climate).
Figure S6:Aedes albopictussuitability based on the MCDA method developed in ECDC, 2009 (Model 2). This is carried out for the period 2030-2050 for each Regional Climate Model projections. The RCMs data are mean bias corrected with the mean reference climate based on the EOBS data for the period 1990-2009 (the delta estimated from the models is added to the mean observed climate).
Figure S7:Weeks of adult mosquito activity based on Medlock et al, 2006 (Model 3). This is carried out for the period 2030-2050 for each regional climate model projections. The RCMs data are mean bias corrected with the mean reference climate based on the EOBS data for the period 1990-2009 (the delta estimated from the models is added to the mean observed climate).
a)Model 1 / b)Model 2 / c)Model 3
/ /
Figure S8:Receiver Operating Characteristic (ROC) curves of the different modelperformances employed in the study (driven by observed climate data for 1990-2009). The analysis is carried out for the grid points where absence and presence of the mosquito has been reported (ECDC/VBORNET dataset).All model outputs have been rescaled between 0 and 1, and all data have been interpolated to the common 0.25° grid.
Country / Year of detection / Reference
USA / 1985 / Hawley et al., 1987
Brazil / 1986 / Forattini, 1986
Mexico / 1988 / Anonymous, 1989
Barbados / 1993 / Reiter, 1998
Dominican republic / 1993 / Peña, 1993
Cuba / 1995 / Broche Borja, 1999
Honduras / 1995 / Anonymous, 1996
Guatemala / 1995 / Ogata Amayoa, 1996
El Salvador / 1996 / Anonymous, 1996
Bolivia / 1997 / Anonymous, 1997
Argentina / 1998 / Rossi et al., 1999
Colombia / 1998 / Velez et al., 1998
Paraguay / 1999 / Anonymous, 1999
Panama / 2002 / Anonymous, 2002
Uruguay / 2003 / Anonymous, 2003
Nicaragua / 2003 / Lugo et al., 2005
South Africa / 1990 / Cornel Hunt, 1991
Nigeria / 1991 / Savage et al., 1992
Cameroon / 2000 / Fontenille Toto, 2001
Equatorial Guinea / 2001 / Toto et al., 2003
Gabon / 2006 / Krueger Hagen, 2007
Lebanon / 2003 / Anonymous, 2003
Israel / 2003 / Pener et al., 2003
Syria / 2005 / Anonymous, 2006
Albania / 1979 / Adhami Reiter, 1998
Italy / 1990 / Dalla Pozza et al., 1994
France / 1999 / Schaffner Karch, 2000
Belgium / 2000 / Schaffner et al., 2004
Montenegro / 2001 / Petric et al., 2001
Switzerland / 2003 / Flacio et al., 2004
Greece / 2003 / Samanidou-Voyadjoglouet al., 2005
Spain / 2004 / Aranda et al., 2006
Croatia / 2004 / Klobucaret al., 2006
Netherlands / 2005 / Sholte et al., 2007
Bosnia and Herzegovina / 2005 / ECDC, 2006
Slovenia / 2005 / Petric et al., 2006
Germany / 2007 / Pluskotaet al., 2008
Malta / 2009 / Gatt et al., 2009
Table 1: Countries where A. albopictus was detected from 1979 to 2011.
Related references:
Adhami, J. and P. Reiter, 1998. Introduction and establishment of Aedes albopictus in Albania. J. Am. Mosq. Contr. Assoc. 14: 340-343.
Anonymous, 1989. Update: Aedes albopictus infestation—United States, Mexico. MMWR38:445–446.
Aranda, C., R. Eritja, and D. Roiz, 2006. First record and establishment of the mosquito Aedes albopictus in Spain. Med. Vet. Entomol. 20: 150-152.
Broche R.G. and Borja E.M., 1999. Aedes albopictus in Cuba. J. Am. Mosq. Control. Assoc. 15:569–570.
Burgess, N.R.H. 1995. Aedes albopictus: a potential problem in the United Kingdom. Parasitologia,37: 121-122.
Cornel A., Hunt R.H., 1991. Aedes albopictus in Africa? First records of live specimens in imported tires in Cape Town. J. Am. Mosq. Control. Assoc. 7:107–108.
Dalla Pozza, G.L., R. Romi, and C. Severini. 1994. Source and spread of Aedes albopictus in the Veneto region of Italy. J. Am. Mosq. Contr. Assoc. 10: 589-592.
ECDC, 2006. Consultation on Chikungunya risk assessment for Europe, Stockholm, 30 March 2006, 9pp.[ Accessed 6 June 2006.
Fisher D., S.M. Thomas, F. Niemitz, B. Reineking, C. Beierkuhnlein, 2011. Projection of Climate suitability for Aedes albopictus Skuse (Culicidae) in Europe under climate change conditions. Global and Planetary Change, 78:54-65.
Flacio, E., P. Lüthy, N. Patocchi, F. Guidotti, M. Tonolla, and R. Peduzzi. 2004. Primo ritrovamento di Aedesalbopictus in Svizzera. Boll. Soc. Ticinese Sc. Nat. 92: 141-142.
Fontenille D. and Toto J.C., 2001. Aedes (Stegomyia) albopictus (Skuse), a potential new dengue vector in southern Cameroon. Emerg. Infect. Dis. 7:1066–1067.
Forattini O.P., 1986. Identificacao de Aedes (Stegomyia) albopictus (Skuse) no Brasil. Rev Saude Pub de Sao Paulo20:244–245.
Gatt, P., Deeming, J.P. & Schaffner F., 2009. First record of Aedes (Stegomyia) albopictus (Skuse) (Diptera: Culicidae) in Malta. Eur. Mos. Bull.,27: 56-64.
Hawley, W.A., P. Reiter, R.S. Copeland, C.B. Pumpuni, and G.B. Craig., 1987. Aedes albopictus in North America: Probable introduction in used tires from Northern Asia. Science236: 1114-1115.
Klobucar, A., E. Merdic, N. Benic, Z. Baklaic, and S. Krcmar. 2006. First record of Aedes albopictus in Croatia. J. Am. Mosq. Cont. Assoc. 22: 147-148.
Knudsen, A.B., R. Romi, and G. Majori. 1996. Occurrence and spread in Italy of Aedes albopictus, with implications for its introduction into other parts of Europe. J. Am. Mosq. Contr. Assoc., 12: 177-183.
Kobayashi, M., N. Nihei, and T. Kurihara, 2002. Analysis of northern distribution of Aedes albopictus in Japan by Geographical Information Systems. J. Med. Entomol., 39: 4-11.
Krueger A. and R.M. Hagen, 2007. Short communication: First record of Aedes albopictus inGabon, Central Africa. Tropical Medicine and International Health12: S. 1105–1107.
Lugo Edel C., Moreno G., Zachariah M.A., López M.M., et al., 2005. Identification of Aedes albopictus in urban Nicaragua. J. Am. Mosq. Control. Assoc. 21:325–327.
Mitchell C., 1995. Geographic spread of Aedes albopictus and potential for involvement in arbovirus cycles in the Mediterranean Basin. J. Vector Ecol.,20:44–58.
Ogata K., Samayoa A.L., 1996. Discovery of Aedes albopictus in Guatemala. J. Am. Mosq. Control Assoc. 12:503–506.
Peña CJ. 1993. First report of Aedes (Stegomyia) albopictus (Skuse) from the Dominican Republic. Vector Ecol. Newsl., 24.
Pener, H., A. Wilamowski, H. Schnur, L. O rshan, U. Shalom, and A. Bear. 2003. Letter to the Editors. Eur. Mos. Bull. 14: 32.
Petric, D., I. Pajovic, A. Cupina, and M. Zgomba. 2001. Aedes albopictus in Entomofauna of Yugoslavia. Symposia ofSerbian Entomologist’s 2001, Entomological Societyof Serbia, Goc. 26-29 September 2001.
Petric D., Zgomba M., Ignjatovic Cupina A., Pajovic I., Merdic E., Boca I., Landeka N., 2006. Invasion of the Stegomyia albopicta to a part of Europe, Presentation at the 15th European Society for Vector Ecology Meeting, 10-14 April 2006, Serres, Greece.
Pluskota, B., Storch, V., Braunbeck, T., Beck, M. & Becker N., 2008. First record of Stegomyia albopicta (Skuse) (Diptera: Culicidae) in Germany. Eur. Mos. Bull.,26: 1-5.
Reiter P., 1998. Aedes albopictus and the world trade in used tires, 1988–1995: The shape of things to come? J. Am. Mosq. Control. Assoc., 14:83–94.
Rodhain, F. 1996. Problèmes posés par l’expansion d’Aedes albopictus. Bull. Soc. Pathol. Exot., 89: 137-141.
Rossi G.C., Pascual N.T., Krsticevic F.J., 1999. First record of Aedes albopictus (Skuse) from Argentina. J. Am. Mosq. Control. Assoc. 15:422.
Samanidou-Voyadjoglou, A., E. Patsoula, G. Spanakos, and N.C. Vakalis. 2005. Confirmation of Aedes albopictus (Skuse) (Diptera: Culicidae) in Greece. Eur. Mos. Bull. 19: 10-12.
Savage H.M., Ezike V.I., Nwankwo A.C.N., Miller B.R., 1992. First record of breeding populations of Aedes albopictus in continental Africa: Implications for arboviral transmission. J. Am. Mosq. Control. Assoc. 8:101–103.
Schaffner F. and S. Karch. 2000. Première observation d’Aedes albopictus en France métropolitaine. Comptes Rendus de l’Académie des Sciences, Paris, Sciences de la vie / Life Sciences323: 373-375.
Schaffner F., W. Van Bortel, and M. Coosemans. 2004. First record of Aedes (Stegomyia) albopictus in Belgium. J. Am. Mosq. Contr. Assoc. 20: 201-203.
Scholte E.J., Jacobs F., Linton Y.M., Dijkstra E., Fransen J. And Takken W., 2007. First record of Aedes (Stegomya) albopictus in the Netherlands. Eur. Mos. Bull., 22:5-9.
Toto J.C., Abaga S., Carnevale P., Simard F., 2003. First report of the oriental mosquito Aedes albopictus on the West African island of Bioko, Equatorial Guinea. Med. Vet. Entomol. 17:343–346.
Velez Iñ, Quiñones ML, Suarez M, Olano V, et al., 1998. Presencia de Aedes albopictus en Leticia, Amazonas,Colombia. Biomedica18:182–198.