- Open Access
Recent and projected future climatic suitability of North America for the Asian tiger mosquito Aedes albopictus
Parasites & Vectorsvolume 7, Article number: 532 (2014)
Since the 1980s, populations of the Asian tiger mosquito Aedes albopictus have become established in south-eastern, eastern and central United States, extending to approximately 40°N. Ae. albopictus is a vector of a wide range of human pathogens including dengue and chikungunya viruses, which are currently emerging in the Caribbean and Central America and posing a threat to North America.
The risk of Ae. albopictus expanding its geographic range in North America under current and future climate was assessed using three climatic indicators of Ae. albopictus survival: overwintering conditions (OW), OW combined with annual air temperature (OWAT), and a linear index of precipitation and air temperature suitability expressed through a sigmoidal function (SIG). The capacity of these indicators to predict Ae. albopictus occurrence was evaluated using surveillance data from the United States. Projected future climatic suitability for Ae. albopictus was obtained using output of nine Regional Climate Model experiments (RCMs).
OW and OWAT showed >90% specificity and sensitivity in predicting observed Ae. albopictus occurrence and also predicted moderate to high risk of Ae. albopictus invasion in Pacific coastal areas of the Unites States and Canada under current climate. SIG also well predicted observed Ae. albopictus occurrence (ROC area under the curve was 0.92) but predicted wider current climatic suitability in the north-central and north-eastern United States and south-eastern Canada. RCM output projected modest (circa 500 km) future northward range expansion of Ae. albopictus by the 2050s when using OW and OWAT indicators, but greater (600–1000 km) range expansion, particularly in eastern and central Canada, when using the SIG indicator. Variation in future possible distributions of Ae. albopictus was greater amongst the climatic indicators used than amongst the RCM experiments.
Current Ae. albopictus distributions were well predicted by simple climatic indicators and northward range expansion was predicted for the future with climate change. However, current and future predicted geographic distributions of Ae. albopictus varied amongst the climatic indicators used. Further field studies are needed to assess which climatic indicator is the most accurate in predicting regions suitable for Ae. albopictus survival in North America.
The Asian Tiger mosquito Aedes albopictus Skuse (1894), is an aggressive diurnal-biting insect that is associated with the transmission of over 20 human pathogens including arboviruses and Dirofilaria spp. nematodes ,. Ae. albopictus is listed as one of the top 100 invasive species by the Invasive Species Specialist Group and is considered to be the most invasive mosquito species in the world ,.
Native to South-eastern Asia, Ae. albopictus naturally occurs in a wide range of habitats including coastland, forests, grasslands, urban areas, water courses and wetlands, and has high ecological flexibility being found in densely vegetated rural areas, agricultural areas as well as urban and sub-urban settings. Over the past 30 years this species has been introduced to parts of Europe as well as parts of Africa, Brazil, Central America, the Caribbean, and southern and eastern United States . Its preference for container habitats for breeding, which include used tyres and containers within peri-domestic settings, has promoted its international spread and establishment close to human habitations. International spread is also favoured by its cold-tolerant eggs and capacity to adapt (in terms of diapause of eggs) to temperate environments ,.
In nature and/or the laboratory, Ae. albopictus is a competent vector for a wide range of viral diseases of significance for human health, including those mostly transmitted human-to-human such as dengue and chikungunya viruses, as well as vector-borne zoonoses such as West Nile virus (WNV), Eastern Equine Encephalitis virus, Rift Valley Fever virus, Cache Valley virus and LaCrosse virus . The capacity of Ae. albopictus to feed on a wide range of host species, and to transmit some of these viruses transovarially, enhances its vector potential . The wide range of viruses transmitted by Ae. albopictus means that where it invades it can act as an additional vector of endemic viruses, and permit autochthonous transmission or outbreaks of diseases exotic to the location it has invaded. This capacity is illustrated by the outbreak of chikungunya in Italy in 2007, and autochthonous cases of chikungunya and dengue in a number of locations in Europe ,. Chikungunya and dengue have recently emerged/re-emerged in the Caribbean, to and from where there is considerable trade and travel with North America ,. Even where environmental conditions may be suitable for mosquito vectors and transmission, the introduction of infected mosquitoes or infected people is unlikely to result in sustained transmission of these viruses in most of North America for a range of socio-economic reasons including availability of home air conditioning, urban and building design and human behaviour . However, limited outbreaks or autochthonous cases of these diseases are a possibility where Ae. albopictus populations have become established. Autochthonously-transmitted exotic vector-borne diseases can have a significant public health impact  and there is considerable current need to assess if and where Ae. albopictus populations, and by inference limited outbreaks or autochthonous cases of dengue or chikungunya, could occur.
A number of studies have aimed to predict where Ae. albopictus may be, or invade, under current climate conditions in Europe, Asia and North America, and most identify climate or weather variables (temperature and precipitation) as key determinants of Ae. albopictus distribution ,-. Laboratory-based entomological studies suggest that higher temperatures (accounting for temperature fluctuations ), if not associated with increased desiccation, improve conditions for Ae. albopictus multiplication, survival and activity ,. It is not unreasonable, therefore, to suggest that a warming climate and changes in precipitation in the context of climate change may drive changes in the geographic distribution of Ae. albopictus,. Only a few studies have to date aimed to predict future possible occurrence of Ae. albopictus with climate change ,,, and for North America specifically, to our knowledge there is only one such study, at State-level . Aedes albopictus invaded the southern United States, in Texas in 1985, and genetic analyses suggest that the colonising Ae. albopictus originated in temperate Japan . Ae. albopictus has now thought to have established breeding populations in States that are close to, or border Canada (Figure 1). In Canada, reproducing populations of Ae. albopictus are not known to exist at present (although individual mosquitoes of this species have been found during surveillance ). In this study we aimed to evaluate the potential of this mosquito species to become established more widely in the United States and Canada under current and future climatic conditions, and provide risk of transmission of chikungunya and dengue that are currently exotic to Canada and most of the United States. We assume that introduction of this mosquito into new geographic regions of the United States and Canada is possible either as a consequence of natural expansion of the mosquito’s range from endemic locations in the United States, or by importation by trade, within-North America or with other parts of the world, in products such as house plants and tyres that have been previously implicated as a means by which this species can be introduced ,,.
In this study, indicators of climatic constraints on Ae. albopictus population survival, that have been elucidated previously in studies in Asia and Europe, were used to assess the possible current and future distributions of this mosquito in North America.
Climatological indicators of Ae. albopictus survival
In this study we used three climatological indicators of Ae. albopictus survival to assess current and future climate suitability over the contiguous United States and Canada. These indicators have been used previously to assess current and future climate suitability in Europe for this mosquito .
A criterion used by Caminade et al. that included seasonal activity of Ae. albopictus in Europe was not considered here because series of day length are not uniformly available over the geographic region studied here. The climatological indicators are described as follows:
An indicator (hereafter termed OW) of the suitability of mean air temperature in January (Tjan) for overwinter survival of Ae. albopictus (after Medlock et al. ) as well as the occurrence of sufficient annual precipitation (Pann) for reproduction. OW provides a four point ordinal scale of suitability for Ae. albopictus: (0) very unsuitable if TJan is lower than 0°C and Pann is below 500 mm, followed by the increasing levels of suitability: (1) moderate when 0°C ≤ TJan < 1°C and 500 mm ≤ Pann <600 mm, (2) high when 1°C ≤ TJan < 2°C and 600 mm ≤ Pann <700 mm, and (3) very high when TJan ≥ 2°C and Pann ≥700 mm.
An indicator (after Kobayashi et al.), hereafter termed OWAT, that combines climatic suitability as defined by OW ≥1 with different thresholds of mean annual temperature (Tann). OWAT provided five point ordinal scale: (0) very unsuitable conditions when Tann is below 9°C, (1) low risk when 9°C ≤ Tann < 10°C, (2) moderate risk if 10°C ≤ Tann < 11°C, (3) high risk if 11°C ≤ Tann < 12°C, and (4) totally suitable conditions if Tann ≥ 12°C.
An indicator based on the overwintering and summer temperatures expressed through a sigmoidal function . For this indicator January and summer (June–July–August) temperatures were transformed into an interval ranging between 0 and 255 using sigmoidal functions. For precipitation, suitability was zero when annual precipitation was lower than 450 mm and maximum when precipitation was higher than 800 mm. For summer temperatures, suitability was zero when temperatures were lower than 15°C and higher than 30°C and maximum between 20°C and 25°C. For January temperatures, suitability was zero when temperatures were lower than 2°C, and maximum when temperatures were higher than 3°C. These three parameters that were used to define this indicator were then linearly combined (using the arithmetic average) to define a level of suitability of a set of climatic conditions for Ae. albopictus. The level of suitability according to this indicator (termed SIG hereafter) was finally rescaled to range between 0 and 100.
Observed climate data
The observation-based climatic suitability over North America for the period 1981–2010 was derived using station-based gridded daily observations of temperature and precipitation. To cover both the United States and Canada, two independent high-resolution databases were merged. These databases are the CONUS L2013  for the conterminous United States and ANUSPLIN for Canada south of 60°N ,, and the merged database is termed ANUSPLIN-CONUSL13 in the following.
Geographic and temporal representation of current and future climate data
The geographic domain and horizontal mesh resolution varied in size for the gridded observations, and amongst the climate model outputs. Therefore, for consistency in analyses and comparisons, a common grid projection was designated as a regular latitude-longitude grid with a spatial resolution of 0.44 degrees square, extending over North America between 20°N and 60°N.
For observed and projected climate, the climatic indicators of Ae. albopictus survival were derived from long-term climatic averages of temperature and precipitation over different time windows. Mean temperatures (January, summer [June-July-August] and annual) were calculated for each year and precipitation was aggregated to annual accumulation and 5-year moving averages for each year were calculated for each climate measure to reduce effects of inter-annual variability.
Values for observed or projected temperature and precipitation for each grid cell were used to classify the grid cells into one of the categories of OW and OWAT, and assign the grid cells a value for SIG. Where climate data and climate model output had a different grid projection to the common grid projection (detailed in Table 1), Inverse Distance Weighting (IDW) interpolation (in Climate Data Operators Version 1.6, Max-Planck-Institut für Meteorologie, Hamburg, Germany) was applied to the raw time series of temperature and precipitation. The observed climate data (ANUSPLIN-CONUSL13) interpolated to the common grid are shown in Figure 2.
Validation of indicators and selection of cut-off values
Data on occurrence of Ae. albopictus in surveillance in the United States from 1985–2010 were provided by the United States Centers for Disease Control and Prevention (CDC) at a spatial resolution of county. There is uncertainty in the surveillance data because mosquito surveillance has not been spatio-temporally systematic. While the majority of counties in the United States have undertaken mosquito surveillance during the 1985–2010 period (particularly those of the west coast, the Rocky Mountain region, the upper Midwest, the northeast and the Atlantic coast ) false negative counties likely occur due to absence of mosquito surveillance. The main objective of validation was, therefore, to compare amongst the indicators in their power to predict the observed data and establish if any one of them performed particularly well (suggesting more emphasis should be placed on it) or badly (suggesting it should not be used). The surveillance data were disaggregated into the common grid projection used in all analyses in this study (as described above) and a two- dimensional spatial array of zeros and ones was created: a value of 1 was assigned to each grid cell within a county in which the vector had been detected, and zero if not (Figure 1). Receiver Operator Characteristic (ROC) analysis  was performed in StataSE11 for Windows (College Station, TX, USA) for each of the climatological indicators using presence-absence of Ae. albopictus for each grid cell during 1981–2010 as the outcome. For the ordinal scales of the OW and OWAT indicators, ROC analysis allowed calculation and comparison of the sensitivity and specificity of different categories as cut-off levels for climatic suitability for Ae. albopictus survival. For the continuous SIG indicator, ROC analysis generated an index of the predictive power of SIG (the area under the ROC curve: AUC). ROC analysis also generated specificity and sensitivity values for a range of cut-off values for SIG that could be used to select a cut-off value for classifying suitability of different locations for Ae. albopictus under future projected climate . The Youden index (defined as J = Sensitivity + Specificity – 1 ) was calculated for each of 3092 value points on the ROC curve to provide a single scale of values on which to base selection of cut-off values for climatic suitability using SIG. Two cut-off values for climatic suitability when using the SIG indicator were selected. One value chosen was that giving the maximum value for J, which gave equal weight to sensitivity and specificity of classification given that surveillance data for the United States was not collected systematically. Therefore false negatives could occur, and identification of Ae. albopictus in surveillance could in some circumstances indicate transient individuals rather than permanent reproducing populations of the mosquito. However a second value for J was chosen that gave the highest SIG value for specificity when sensitivity was >90% to provide a more risk-averse assessment of future possible Ae. albopictus distributions.
The disaggregation of county-level surveillance data into grid cells artificially increased the sample size of the validation data available for assessing performance from 3112 counties in the conterminous United States to 5789 grid cells. To assess the extent to which this affected estimates of the performance of the different indicators in predicting occurrence of Ae. albopictus via the area under the ROC curve (AUC), a data set of 3112 grid cells (i.e. a number equal to the number of counties) was selected at random from the total 5789 grid cells. AUC values for OW, OWAT and SIG were then obtained using this reduced 3112 point data set.
Apart from the non-systematic nature of the mosquito surveillance, an additional possible reason for false negative counties in surveillance data is that the mosquito has not yet been introduced even though climatic and other environmental determinants are suitable. To assess the extent to which this possibility could affect estimates of the performance of the different indicators in predicting occurrence of Ae. albopictus via the AUC values, the ROC analyses were repeated for a subset of the data south of 40°N and east of 105°W, which comprises the main region of the United States where Ae. albopictus has been found (Figure 1).
Future projected distributions of Ae. albopictus
To assess effects of climate change on possible future distributions of Ae. albopictus, and evaluate the degree of uncertainty in these projections at a regional scale, projected precipitation and temperature data were obtained from a range of Regional Climate Models (RCMs). The simulations are based on the lateral boundary conditions arising from reanalysis (i.e. a climate or weather model simulation of the past that includes data assimilation of historical observations, see http://reanalyses.org/) for the historical periods or from global climate model (GCM) output for both current and future climate conditions. A time window of 30 years was used to construct climate change scenarios, which, according to availability of RCM output (Table 1), allowed mapping of the climatological indicators of Ae. albopictus survival for the 2020s (2011–2040) and/or 2050s (2041–2070). This approach is consistent with current best practice for studies on impacts of projected climate change . The climate simulations from the nine RCMs used in this study (Table 1) were performed within two main project frameworks. Simulations of four RCMs including CanRCM4 and CRCM5 developed respectively by the Canadian Centre for Climate modelling and analysis division of Environment Canada (CCCma/EC) and by the Centre pour l’Étude et la Simulation du Climat à l’Échelle Régionale (ESCER) at the Université du Québec à Montréal, HIRHAM5 of the Danish Meteorological Institute, and RCA4 of the Swedish Meteorological and Hydrological Institute represent the first available output of an ensemble of RCMs within the Coordinated Regional Climate Downscaling Experiment (CORDEX) project covering the North American Domain . Output from four other RCMs (RegCM3, ECPC, MM5I and WRF) that participated in the North American Regional Climate Change Assessment Program (NARCCAP: ,) and simulations from the CRCM version 4.2.3 , (CRCM4.2.3 runs provided by Ouranos through the CCCma/EC web site) were also considered in our study. In simulations of current climate, RCMs were forced by lateral boundary conditions of the global reanalysis datasets of the European Centre for Medium-Range Weather Forecasts (ECMWF) ERA-40  or ERA-Interim , or the National Centers for Environmental Prediction (NCEP) Department of Energy (DOE) reanalysis II . Climate change projections were driven by six GCM projections under two greenhouse gas (GHG) emission scenarios (Table 1). The new emissions scenarios developed for the Inter-Governmental Panel on Climate Change (IPCC) AR5 were used in the CORDEX project. The emissions, concentrations, and land-cover change projections are described in the Representative Concentration Pathways RCP4.5 and RCP8.5 . The other five RCMs use the previous, but widely-applied emission scenario IPCC SRES A2 . The A2 scenario describes a very heterogeneous world with high population but slower economic growth than in other scenarios. The RCM scenarios were created for plausible change based on near future (2011 to 2040 in the case of the RCP emission scenarios), and mid-term future (2041 to 2070) climatic conditions. The scenarios A2, RCP4.5 and RCP8.5 are quite similar in terms of GHG equivalent concentrations during the first part of the 21st century but GHG concentrations are greater in A2 and RCP8.5 than in RCP4.5 after 2050 ,,.
Values for projected temperature and precipitation for each grid cell were used to classify the grid cells into one of the categories of OW and OWAT, and assign the grid cells a value for SIG. Where RCM output had a different grid projection to the common grid projection IDW interpolation was applied to the raw time series of temperature and precipitation.
Results and discussion
Validation of indicators and selection of cut-off values
Each of the three indicators for climatic suitability for Ae. albopictus performed well in predicting observed Ae. albopictus distributions in the United States, whether using the entire (5789 values) dataset of grid cell values or the reduced (3112 values) dataset (Table 2). For OW, when using a cut-off for prediction of absence and presence between points 0 and 1 on the ordinal scale, sensitivity and specificity of prediction of Ae. albopictus occurrence were both >90%. For OWAT, performance was best when using a cut-off for prediction of absence and presence between points 2 and 3 on the ordinal scale, where sensitivity and specificity of prediction of Ae. albopictus occurrence were 90.2% and 92.2% respectively. For SIG, the area under the ROC curve, when using the reduced dataset, was 0.925 (95% confidence interval 0.913 to 0.936) indicating that SIG was “highly accurate” (using the terminology of Greiner et al. ) in discriminating suitable and non-suitable climate for Ae. albopictus according to the surveillance data. The different indicators of climatic suitability for Ae. albopictus all performed well in comparison with other studies on prediction of Ae. albopictus distributions in other parts of the world and/or using different modelling methods ,, although AUC values for SIG were greater than those for OWAT, and lowest of OW. However, the uncertainties inherent in the surveillance data mean that inter-study comparisons of the ROC AUC values obtained here should be not be over-interpreted.
When using the data from south of 40°N and east of 105°W, sensitivity of all of the indicators was high (>90%) when using the cut-off values described above. However, AUC values were reduced to approximately 0.7 for all indicators due to low specificity (<50%, Table 2). This was anticipated because while the selected area is where most Ae. albopictus-positive locations have been found, it is also an area where there is likely a particularly high number of counties where mosquito surveillance has not occurred and false negative counties occur .While the occurrence of false negative locations may affect AUC values, they should not, however, affect comparisons amongst the indicators. Using the data from south of 40°N and east of 105°W, the order of performance of the indicators changed compared to those obtained using the full data, with OW having the highest AUC values, SIG having the second highest and OWAT the lowest values (Table 2). This suggests that it would be prudent not to reject the possibility of risk of Ae. albopictus becoming established in a particular location on the basis of the findings of only one indicator, and to determine risk using output from all indicators.
Geographic extent of possible current distributions of Ae. albopictus
Climatic suitability maps for Ae. albopictus using OW, OWAT and SIG from observed climate data and one representative RCM (CanRCM4) are shown in Figure 3. The OW and OWAT indicators suggested similar geographic regions of suitable climate including one block of similar, very high climatic suitability extending from Florida to approximately 40° N (with more northerly regions on the Atlantic coast, including Long Island, being suitable) and 105° W in the United States (Figure 3). The western parts of Pacific coast States and Provinces from mid-California to southern coastal British Columbia were also mostly of very high climatic suitability (Figure 3). There were some slight differences with a smaller geographic extent of climatic suitability in States and Provinces bordering the Pacific coast using the OWAT indicator compared to the OW indicator, and the OWAT indicator identified some patches of climatic suitability in States between the mid-west and the Pacific coastal States that were not identified by the OW indicator. The SIG indicator identified a block of climatic suitability extending from Florida to approximately 40° N (with more northerly regions on the Atlantic coast, including Long Island, being suitable) and westwards to 105° W in the United States that was similar to that predicted by OW and OWAT. The SIG indicator also identified western parts of Pacific coast States and Provinces from mid-California to southern coastal British Columbia as climatically suitable although these areas were of geographic limits similar to those predicted by OWAT. In contrast to OW and OWAT, the SIG indicator predicted climatic suitability for regions of more northern States east of approximately 97° W to the Atlantic coast, and northward into southern Ontario, Quebec and the Maritimes in Canada (Figure 3). This was because the values of the SIG indicator are not constrained to zero by sub-zero January temperatures, and because of suitable rainfall and summer temperatures in these regions (Figure 2). With the exception of a few regional differences, the re-analysis outputs of the RCM output for 1989–2010 produced similar patterns of climatological suitability, for each of the three indicators, to values obtained using observed data (Figure 3).
Therefore, in summary, predicted climate suitability using OW and OWAT was similar and conservative, suggesting that the only part of the United States at risk of Ae. albopictus populations is that where this mosquito is known to have become established (in the south-eastern corner of the continent) or on the Pacific coast where the mosquito is known to have become established in the past and actively eradicated (such as Washington State ). In contrast, the SIG indicator suggested additional climatic suitability in northern states in the eastern United States and in southern Ontario, Quebec and the Maritimes in Canada. If the current northern limit of Ae. albopictus populations is well described by the surveillance data in the United States, then SIG likely overestimates risk further north and OW or OWAT are more useful criteria for describing risk of Ae. albopictus populations becoming established. However, if Ae. albopictus populations do occur north of the limit observed to date in surveillance, then SIG may be the more useful criterion on which to assess risk in northern United States and Canada. Some observations in the field suggest that OW and OWAT do not under-estimate risk. In Italy, the equivalent of OWAT cut-off point 2 best described Ae. albopictus population distributions  and in experimental studies in Connecticut just north of the northern limit of climatic suitability according to OW and OWAT criteria, Ae. albopictus populations failed to become established at least in part due to failure to overwinter . At the northern edge of the geographic range of Ae. albopictus according to the surveillance data, at least at one location ongoing field studies support the idea that the presence of Ae. albopictus in surveillance data equates with the presence of reproducing populations of the mosquito . Furthermore, niche modelling studies using global Ae. albopictus distribution data produce a similar pattern of distribution in North America similar to that predicted using the OW and OWAT indicators . However, because surveillance for Ae. albopictus is not, and has not been, systematic in space and time over the United States and in Canada, the true extent of false negative locations in the surveillance data is unknown.
Other factors may limit the predictive power of the methods used here. First, and in common with all “pattern matching” predictive modelling techniques, Ae. albopictus populations are spreading in the United States and the observed surveillance data likely represent the current “realised niche” rather than the full theoretical climatological niche width for this species . Evolution of photoperiodic responses of mosquitoes, including of Ae. albopictus populations in the United States, has been observed and these (particularly development of egg diapause over winter) may be critical to Ae. albopictus population survival in the northern parts of its range in the United States ,. The possibility of such evolutionary adaptations to changing climatic conditions could call into question the validity of assessing future distributions on those observed in the past by changing the climatological niche width. However, northern populations of Ae. albopictus have likely already evolved the over-winter survival-enhancing trait of egg diapause  and that trait should be accounted for in the surveillance data used here. Interactions amongst mosquito species (particularly competition) likely will also impact the realised niche width and the validity of projections here (e.g. ), and are illustrative that the range of ecological determinants of environmental suitability for vectors extends beyond climatic suitability.
Geographic extent of projected future distributions of Ae. albopictus
Projections of potential future climatic suitability for Ae. albopictus using output from CanRCM4 are shown as an illustration in Figure 4. Projections of future climatic suitability from all models are presented in Additional file 1. Projected changes in mean annual temperatures and annual precipitation are shown for reference in Figure 5. Projected northward expansion of the geographic range of Ae. albopictus (by approximately 200–500 km) was modest using OW and OWAT (Figure 4, Additional file 1), at least in comparison to the range spread projected for other arthropod vectors with climate change . The degree of northward range expansion using OW and OWAT was of a similar magnitude for all RCMs although the current and projected final northern limits by 2041–2070 varied amongst the RCMs. The most northerly projected climatically suitable locations were in Montana, North Dakota, Wisconsin, Michigan, Philadelphia, New York State and New England States in the United States and coastal British Columbia, southern Ontario, Quebec and the Maritimes in Canada (Figure 4, Additional file 1). These northern limits were only projected by climate models that projected the most northern possible distribution of climatic suitability under current climate (e.g. CanRCM4, RCA4, Additional file 1). The models that projected more southerly northern limits under current climate (e.g. CRCM4.2.3, ECPC and MM5I, Additional file 1) projected more southerly northern limits (except along the Pacific coast) with only small projected localised areas of climatic suitability in central and eastern regions of northern United States and southern Canada in the future (Figure 4, Additional file 1). A finding of moderate projected range expansion is consistent with region-scale projections for Ae. albopictus in the north-eastern United States using climatological niche predicted using maximum entropy methods .
Using the SIG indicator, the projections for future climate suitability in the United States and Canada eastward from 100° W were more complex compared to projections using the OW and OWAT indicators. Using this indicator northward expansion was more extensive (up to 1000 km) into Canada (Figure 4), a band of unsuitable climate dividing the climatically suitable region of the United States was apparent, and climatic suitability in Texas, Louisiana, Oklahoma and Arkansas was reduced due to a combination of reduced rainfall and increased summer temperature (Figure 5). The SIG indicator also predicted more extensive future climate suitability in States along the United States Pacific coast and southern British Columbia compared to OW and OWAT indicators and predicted a future area of climatic suitability in the eastern foothills of the southern Rocky Mountains (Figure 4, Additional file 1). As for the OW and OWAT indicators, there was variation amongst the RCMs in the extent of the future projected northern limit of climate suitability (Figure 6). It could be argued that the SIG overestimates risk by not accounting for absolute limits on population establishment associated with very cold temperatures in winter and very dry conditions. However, perhaps SIG provide a better indication of the climatic suitability in some urban and sub-urban areas where refugia from extremes of climate, including very low mid-winter temperatures, may exist ,.
Throughout, variations in projections were greater amongst models than amongst emissions scenarios (Figures 4 and 6, Additional file 1). This would be expected as the RCPs are quite similar to one another during first part of the 21st century as described above, while RCM simulations differ due to differences in their sub-grid scale processes or parameterizations, and differences in their GCM driving conditions (e.g. ). Further work is needed to evaluate and quantify the uncertainties arising from variation amongst different RCM outputs, and the power of RCMs to identify suitable temperature and precipitation conditions that are sensitive to more local scale forcing or regional-scale influences not taken into account in the geographic scale of the present study.
In this study, the OW and OWAT indicators currently seem to offer the best fit to existing data on Ae. albopictus distribution in the United States, but predictions of climatic suitability using SIG must be considered as surveillance to date has not been geographically systematic and consistent. If OW and/or OWAT are the most accurate indicators of climatic suitability for Ae. albopictus then the possibility of geographic range expansion of this species under current and mid- and long-term future climate is relatively limited with the possible exception of locations along the Pacific coast. However if SIG is the more accurate indicator of climatic suitability for Ae. albopictus then more geographically widespread expansion of the northern range of this species may occur where other environmental determinants allow the mosquito to become established. Additional systematic field studies and surveillance will be needed, therefore, to identify which climatic indicator is the most accurate at predicting climate suitable for Ae. albopictus and more accurately define the climatic and other environmental determinants of this mosquito  to better model and predict its current and future geographic distributions.
Fifth assessment report of IPCC
Area under the ROC curve
Canadian centre for climate modelling and analysis/environment Canada
Coordinated regional climate downscaling experiment
Department of energy
European centre for medium-range weather forecasts
Centre pour l’Étude et la Simulation du Climat à l’Échelle Régionale
Global climate model
Inverse distance weighting
Inter-governmental panel on climate change
North American Regional climate change assessment program
National centers for environmental prediction
Regional climate model
Representative concentration pathway
Receiver operator characteristic
Special report on emissions scenarios
Gratz NG: Critical review of the vector status of Aedes albopictus. Med Vet Entomol. 2004, 18: 215-227. 10.1111/j.0269-283X.2004.00513.x.
European Environment Agency: The impacts of invasive alien species in Europe, technical report No 16/2012. ᅟ. ᅟ, Source:http://www.eea.europa.eu/publications/impacts-of-invasive-alien-species, [http://www.eea.europa.eu/publications/impacts-of-invasive-alien-species]
Lowe S, Browne M, Boudjelas S, De Poorter M: 100 of the World’s worst invasive alien species. A selection from the global invasive species database. ᅟ. ᅟ, Published by The Invasive Species Specialist Group (ISSG) a specialist group of the Species Survival Commission (SSC) of the World Conservation Union (IUCN), ᅟ, [http://www.issg.org/database/species/reference_files/100English.pdf], [http://www.issg.org/database/species/reference_files/100English.pdf]
European Centre for Disease Prevention and Control:Aedes albopictus. ᅟ. ᅟ, [http://ecdc.europa.eu/en/healthtopics/aedes_albopictus/Pages/index.aspx], [http://ecdc.europa.eu/en/healthtopics/aedes_albopictus/Pages/index.aspx]
Leisnham PT, Towler L, Juliano SA: Geographic variation of photoperiodic diapause but not adult survival or reproduction of the invasive mosquito Aedes albopictus (Diptera: Culicidae) in North America. Ann Entomol Soc Am. 2011, 104: 1309-1318. 10.1603/AN11032.
Benedict MQ, Levine SR, Hawley AW, Lounibos AW: Spread of the tiger: global risk of invasion by the mosquito Aedes albopictus. Vector Borne Zoonotic Dis. 2007, 7: 76-85. 10.1089/vbz.2006.0562.
Charrel RN, De Lamballerie X, Raoult D: Chikungunya outbreaks - the globalization of vectorborne diseases. N Engl J Med. 2007, 356: 769-771. 10.1056/NEJMp078013.
La Ruche G, Souarès Y, Armengaud A, Peloux-Petiot F, Delaunay P, Desprès P, Lenglet A, Jourdain F, Leparc-Goffart I, Charlet F, Ollier L, Mantey K, Mollet T, Fournier JP, Torrents R, Leitmeyer K, Hilairet P, Zeller H, Van Bortel W, Dejour-Salamanca D, Grandadam M, Gastellu-Etchegorry M: First two autochthonous dengue virus infections in metropolitan France. Euro Surveill. 2010, 15: 19676-
Boggild AK, Geduld J, Libman M, Ward BJ, McCarthy AE, Doyle PW, Ghesquiere W, Vincelette J, Kuhn S, Freedman DO, Kain KC: Travel-acquired infections and illnesses in Canadians: surveillance report from CanTravNet surveillance data, 2009–2011. Open Med. 2014, 8: e20-e32.
Khan K, Bogoch I, Brownstein JS, Miniota J, Nicolucci A, Hu W, Nsoesie EO, Cetron M, Creatore MI, German M, Wilder-Smith A: Assessing the origin of and potential for international spread of chikungunya virus from the Caribbean. PLoS Curr. 2014, 6: 6-
Reiter P: Climate change and mosquito-borne disease. Environ Health Perspect. 2001, 109 (Suppl 1): 141-161. 10.1289/ehp.01109s1141.
Berrang-Ford L, Maclean JD, Gyorkos TW, Ford JD, Ogden NH: Climate change and malaria in Canada: a systems approach. Interdiscip Perspect Infect Dis. 2009, 2009: 385487-
Kobayashi M, Nihei N, Kurihara T: Analysis of northern distribution of Aedes albopictus (Diptera: Culicidae) in Japan by Geographical Information System. J Med Entomol. 2002, 39: 4-11. 10.1603/0022-2585-39.1.4.
Medlock JM, Avenell D, Barrass I, Leach S: Analysis of the potential for survival and seasonal activity of Aedes albopictus (Diptera: Culicidae) in the United Kingdom. J Vector Ecol. 2006, 31: 292-304. 10.3376/1081-1710(2006)31[292:AOTPFS]2.0.CO;2.
Roiz D, Neteler M, Castellani C, Arnoldi D, Rizzoli A: Climatic factors driving invasion of the tiger mosquito (Aedes albopictus) into new areas of Trentino, Northern Italy. PLoS One. 2011, 6: e14800-10.1371/journal.pone.0014800.
Wu F, Liu Q, Lu L, Wang J, Song X, Ren D: Distribution of Aedes albopictus (Diptera: Culicidae) in northwestern China. Vector Borne Zoonotic Dis. 2011, 11: 1181-1186. 10.1089/vbz.2010.0032.
Brady OJ, Johansson MA, Guerra CA, Bhatt S, Golding N, Pigott DM, Delatte H, Grech MG, Leisnham PT, Maciel-de-Freitas R, Styer LM, Smith DL, Scott TW, Gething PW, Hay SI: Modelling adult Aedes aegypti and Aedes albopictus survival at different temperatures in laboratory and field settings. Parasit Vectors. 2013, 6: 351-10.1186/1756-3305-6-351.
Alto BW, Juliano SA: Precipitation and temperature effects on populations of Aedes albopictus (Diptera: Culicidae): implications for range expansion. J Med Entomol. 2001, 38: 646-656. 10.1603/0022-2585-38.5.646.
Alto BW, Juliano SA: Temperature effects on the dynamics of Aedes albopictus (Diptera: Culicidae) populations in the laboratory. J Med Entomol. 2001, 38: 548-556. 10.1603/0022-2585-38.4.548.
Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. 2013, Cambridge University Press, Cambridge
Higa Y: Dengue vectors and their spatial distribution. Trop Med Health. 2011, 39 (Suppl 4): 17-27. 10.2149/tmh.2011-S04.
Caminade C, Medlock JM, Ducheyne E, McIntyre KM, Leach S, Baylis M, Morse AP: Suitability of European climate for Asian tiger mosquito Aedes albopictus: recent trends and future scenarios. J R Soc Interface. 2012, 9: 2708-2717. 10.1098/rsif.2012.0138.
Fisher D, Thomas SM, Niemitz F, Reineking B, Beierkuhnlein C: Projection of climate suitability for Aedes albopictus Skuse (Culicidae) in Europe under climate change conditions. Glob Planet Change. 2011, 78: 54-65. 10.1016/j.gloplacha.2011.05.008.
Rochlin I, Ninivaggi DV, Hutchinson ML, Farajollahi A: Climate change and range expansion of the Asian tiger mosquito (Aedes albopictus) in Northeastern USA: implications for public health practitioners. PLoS One. 2013, 8: e60874-10.1371/journal.pone.0060874.
ᅟ: Public Health Ontario Vector-Borne Disease Surveillance Report. ᅟ
Linthicum KJ, Kramer VL, Madon MB, Fujioka K: Introduction and potential establishment of Aedes albopictus in California in 2001. J Am Mosq Control Assoc. 2003, 19: 301-308.
European Center for Disease Prevention and Control: Meeting report, expert meeting on chikungunya modelling. ᅟ. ᅟ, [http://www.ecdc.europa.eu/en/publications/Publications/0804_MER_Chikungunya_Modelling.pdf], [http://www.ecdc.europa.eu/en/publications/Publications/0804_MER_Chikungunya_Modelling.pdf]
European Center for Disease Prevention and Control: Development of Aedes albopictus Risk Maps. ᅟ
Livneh B, Rosenberg EA, Lin C, Nijssen B, Mishra V, Andreadis KM, Maurer EP, Lettenmaier DP: A long-term hydrologically based dataset of land surface fluxes and states for the conterminous United States: update and extensions. J Climate. 2013, 26: 9384-9392. 10.1175/JCLI-D-12-00508.1.
Hutchinson MF, McKenney DW, Lawrence K, Pedlar J, Hopkinson R, Milewska E, Papadopol P: Development and testing of Canada-wide interpolated spatial models of daily minimum/maximum temperature and precipitation for 1961–2003. J Appl Meteorol Climatol. 2009, 48: 725-741. 10.1175/2008JAMC1979.1.
McKenney DW, Hutchinson MF, Papadopol P, Lawrence K, Pedlar J, Campbell K, Milewska E, Hopkinson R, Price D, Owen T: Customized spatial climate models for North America. Bull Am Meteorol Soc. 2011, 92: 1612-1622. 10.1175/2011BAMS3132.1.
United States Geological Survey 2014: Arboviral disease maps. ᅟ. ᅟ, [http://diseasemaps.usgs.gov/], [http://diseasemaps.usgs.gov/]
Hanley JA, McNeil BJ: The meaning and use of the area under a receiver operating characteristic (ROC) curve. Radiology. 1982, 143: 29-36. 10.1148/radiology.143.1.7063747.
Greiner M, Pfeiffer D, Smith RD: Principles and practical application of the receiver-operating characteristic analysis for diagnostic tests. Prev Vet Med. 2000, 45: 23-41. 10.1016/S0167-5877(00)00115-X.
Déandreis C, Pagé C, Braconnot P, Bärring L, Bucchignani E, de Cerff WS, Hutjes R, Joussaume S, Mares C, Planton S, Plieger M: Towards a dedicated impact portal to bridge the gap between the impact and climate communities: lessons from use cases. Clim Change. 2014, 125: 333-347. 10.1007/s10584-014-1139-7.
ᅟ: Coordinated Regional Climate Downscaling Experiment (CORDEX). ᅟ
Mearns LO, Gutowski WJ, Jones R, Leung L-Y, McGinnis S, Nunes AMB, Qian Y: A regional climate change assessment program for North America. EOS. 2009, 90: 311-312. 10.1029/2009EO360002.
Mearns LO: The North American Regional Climate Change Assessment Program Dataset. ᅟ
Laprise R, Caya D, Bergeron G, Giguère M: The formulation of André Robert MC2 (Mesoscale Compressible Community) model. Atmos Ocean. 1997, 35: 195-220. 10.1080/07055900.1997.9687348.
Music B, Caya D: Evaluation of the hydrological cycle over the Mississippi river basin as simulated by the Canadian Regional Climate Model (CRCM). J Hydromet. 2007, 8: 969-988. 10.1175/JHM627.1.
Uppala SM, Kållberg PW, Simmons AJ, Andrae U, Bechtold VDC, Fiorino M, Gibson JK, Haseler J, Hernandez A, Kelly G, Li X, Onogi K, Saarinen S, Sokka N, Allan RP, Andersson E, Arpe K, Balmaseda MA, Beljaars ACM, Berg LVD, Bidlot J, Bormann N, Caires S, Chevallier F, Dethof A, Dragosavac M, Fisher M, Fuentes M, Hagemann S, Hólm EV: The ERA-40 re-analysis. Quart J R Meteorol Soc. 2005, 131: 2961-3012. 10.1256/qj.04.176.
Uppala S, Dee D, Kobayashi S, Berrisford P, Simmons A: Towards a climate data assimilation system: status update of ERA-interim. ECMWF Newsl. 2008, 115: 12-18.
Dee DP, Uppala SM, Simmons AJ, Berrisford P, Poli P, Kobayashi S, Andrae U, Balmaseda MA, Balsamo G, Bauer P, Bechtold P, Beljaars ACM, van de Berg L, Bidlot J, Bormann N, Delsol C, Dragani R, Fuentes M, Geer AJ, Haimberger L, Healy SB, Hersbach H, Hólm EV, Isaksen L, Kållberg P, Köhler M, Matricardi M, McNally AP, Monge-Sanz BM, Morcrette J-J: The ERA-Interim reanalysis: configuration and performance of the data assimilation system. Quart J R Meteorol Soc. 2011, 137: 553-597. 10.1002/qj.828.
Kanamitsu M, Ebisuzaki W, Woollen J, Yang S-K, Hnilo JJ, Fiorino M, Potter GL: NCEP-DOE AMIP-II reanalysis (R-2). Bull Am Meteorol Soc. 2002, 83: 1631-1643. 10.1175/BAMS-83-11-1631.
Moss RH, Edmonds JA, Hibbard KA, Manning MR, Rose SK, van Vuuren DP, Carter TR, Emori S, Kainuma M, Kram T, Meehl GA, Mitchell JF, Nakicenovic N, Riahi K, Smith SJ, Stouffer RJ, Thomson AM, Weyant JP, Wilbanks TJ: The next generations of scenarios for climate change research and assessment. Nature. 2010, 463: 747-756. 10.1038/nature08823.
Nakicenovic N, Swart R: Special Report on Emissions Scenarios: A Special Report of Working Group III of the Intergovernmental Panel on Climate Change. 2000, Cambridge University Press, Cambridge
Solomon S, Qin D, Manning M, Chen Z, Marquis M, Averyt KB, Tignor M, Miller HL: Climate Change 2007: The Physical Science Basis. Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change. 2007, Cambridge University Press, Cambridge
Moore CG:Aedes albopictus in the United States: current status and prospects for further spread. J Am Mosq Control Assoc. 1999, 15: 221-227.
Andreadis TG: Failure of Aedes albopictus to overwinter following introduction and seasonal establishment at a tire recycling plant in the northeastern USA. J Am Mosq Control Assoc. 2009, 25: 25-31. 10.2987/08-5813.1.
Swanson J, Lancaster M, Anderson J, Crandell M, Haramis L, Grimstad P, Kitron U: Overwintering and establishment of Aedes albopictus (Diptera: Culicidae) in an urban La Crosse virus enzootic site in Illinois. J Med Entomol. 2000, 37: 454-460. 10.1603/0022-2585(2000)037[0454:OAEOAA]2.0.CO;2.
Farajollahi A, Nelder MP: Changes in Aedes albopictus (Diptera: Culicidae) populations in New Jersey and implications for arbovirus transmission. J Med Entomol. 2009, 46: 1220-1224. 10.1603/033.046.0533.
Urbanski J, Mogi M, O’Donnell D, DeCotiis M, Toma T, Armbruster P: Rapid adaptive evolution of photoperiodic response during invasion and range expansion across a climatic gradient. Am Nat. 2012, 179: 490-500. 10.1086/664709.
Johnson BJ, Sukhdeo MV: Successional mosquito dynamics in surrogate treehole and ground-container habitats in the northeastern United States: where does Aedes albopictus fit in?. J Vector Ecol. 2013, 38: 168-174. 10.1111/j.1948-7134.2013.12023.x.
Ogden NH, Maarouf A, Barker IK, Bigras-Poulin M, Lindsay LR, Morshed MG, O’Callaghan CJ, Ramay F, Waltner-Toews D, Charron DF: Climate change and the potential for range expansion of the Lyme disease vector Ixodes scapularis in Canada. Int J Parasitol. 2006, 36: 63-70. 10.1016/j.ijpara.2005.08.016.
Zhong D, Lo E, Hu R, Metzger ME, Cummings R, Bonizzoni M, Fujioka KK, Sorvillo TE, Kluh S, Healy SP, Fredregill C, Kramer VL, Chen X, Yan G: Genetic analysis of invasive Aedes albopictus populations in Los Angeles County, California and its potential public health impact. PLoS One. 2013, 8: e68586-10.1371/journal.pone.0068586.
Eum H-I, Gachon P, Laprise R: Developing a likely climate scenario from multiple regional climate model simulations with an optimal weighting factor. Clim Dyn. 2014, 43: 11-35. 10.1007/s00382-013-2021-4.
Petrić D, Bellini R, Scholte EJ, Rakotoarivony LM, Schaffner F: Monitoring population and environmental parameters of invasive mosquito species in Europe. Parasit Vectors. 2014, 7: 187-10.1186/1756-3305-7-187.
We thank the United States CDC Vector-Borne Disease Division, particular Dr Jennifer Lehman, and Dr Chester Moore (now at the Department of Microbiology, Immunology & Pathology, Veterinary Medicine & Biomedical Sciences, Colorado State University) for providing the surveillance data used in this article.
The ANUSPLIN data set has been developed and kindly provided by Dr. Dan McKenney of the Canadian Forest Service, Natural Resources Canada.
We thank the North American Regional Climate Change Assessment Program (NARCCAP) for providing the data used in this paper. NARCCAP is funded by the National Science Foundation (NSF), the U.S. Department of Energy (DoE), the National Oceanic and Atmospheric Administration (NOAA), and the U.S. Environmental Protection Agency Office of Research and Development (EPA).
We acknowledge the Data Access Integration Team for providing the data and technical support. The DAI data download gateway is made possible through collaboration among the Global Environmental and Climate Change Centre (GEC3), the Adaptation and Impacts Research Division (AIRD) of Environment Canada, and the Drought Research Initiative (DRI). The Ouranos Consortium (in Québec) also provides IT support to the DAI team.
The CRCM time series data were generated and supplied by Ouranos’ Climate Simulations Team. We would also like to acknowledge our colleagues from Environment Canada and UQAM for providing the CanRCM4 and CRCM5 data, with a special thanks to Dr. John Scinocca and Dr Yanjun Jiao from CCCma/EC, andPr. Laxmi Sushama and Katja Winger from ESCER/UQAM.
The authors declare that they have no competing interests.
NHO designed the study, undertook statistical analyses and led writing of the article. MR, supported by PG obtained and developed the current climate data, undertook all geomatic manipulations of the data, obtained climate model data and obtained current and future projected values for the climatological indicators from RCM output. CC advised on the original design and contributed to interpreting the findings. All authors read and approved the final version of the manuscript.