Skip to main content

First assessment of potential distribution and dispersal capacity of the emerging invasive mosquito Aedes koreicus in Northeast Italy



Invasive alien species represent a growing threat for natural systems, economy and human health. Active surveillance and responses that readily suppress newly established colonies are effective actions to mitigate the noxious consequences of biological invasions. However, when an exotic species establishes a viable population in a new area, predicting its potential spread is the most effective way to implement adequate control actions. Emerging invasive species, despite monitoring efforts, are poorly known in terms of behaviour and capacity to adapt to the new invaded range. Therefore, tools that provide information on their spread by maximising the available data, are critical.


We apply three different approaches to model the potential distribution of an emerging invasive mosquito, Aedes koreicus, in Northeast Italy: 1) an automatic statistical approach based on information theory, 2) a statistical approach integrated with prior knowledge, and 3) a GIS physiology-based approach. Each approach possessed benefits and limitations, and the required ecological information increases on a scale from 1 to 3. We validated the model outputs using the only other known invaded area in Europe. Finally, we applied a road network analysis to the suitability surface with the highest prediction power to highlight those areas with the highest likelihood of invasion.


The GIS physiological-based model had the highest prediction power. It showed that localities currently occupied by Aedes koreicus represent only a small fraction of the potentially suitable area. Furthermore, the modelled niche included areas as high as 1500 m a.s.l., only partially overlapping with Aedes albopictus distribution.


The simulated spread indicated that all of the suitable portion of the study area is at risk of invasion in a relatively short period of time if no control policies are implemented. Stochastic events may further boost the invasion process, whereas competition with Aedes albopictus may limit it. According to our analysis, some of the major cities in the study area may have already been invaded. Further monitoring is needed to confirm this finding. The developed models and maps represent valuable tools to inform policies aimed at eradicating or mitigating Aedes koreicus invasion in Northeast Italy and Central Europe.


An increasing number of species is rapidly spreading outside of their original distributional range and invading new territories, gaining the name of invasive species. The factors underpinning invasion processes are numerous and include socio-economical determinants linked to the intensified speed and density of transcontinental commercial and tourist fluxes [1, 2]. Among abiotic factors, anomalous climatic fluctuations [3] and landscape perturbations [4], mainly due to human exploitation of the environment, modify ecological conditions. These altered conditions trigger or facilitate species mixing at various spatial scales, at times resulting in novel ecosystems which are constituted by persistent assemblages of exotic and indigenous species [5, 6]. Among the risks arising from the increased or shifted geographical distribution of species as well as from novel ecosystems [2, 7, 8], the spread rate of infectious diseases is the most pressing for human health [9]. Indeed, many vertebrate and invertebrate species are competent hosts for one or multiple zoonoses –infectious and parasitic diseases transmissible from animals to humans, whose distribution strictly follow the geographical range of their host species [10].

Bloodsucking arthropods, such as mosquitoes, represent the majority of the organisms able to transmit agents of infectious diseases to humans [11]. Indeed they constitute a system able to overpass the skin barrier and deliver the agent of the disease directly into the blood vessels. They use host blood not only as a food source but also to regulate metabolic processes that cause dramatic and key changes in their physiology [12]. Among them, mosquitoes have been the most successful invasive disease-vector group in the 20th century, and are bridge-vectors of infectious pathogens (e.g., arboviruses) which have caused devastating anthropozoonosis. Some arboviruses, such as dengue fever, Rift Valley fever, yellow fever and chikungunya are transmitted by Aedes species. These are highly invasive container-breeding mosquitoes, with a native geographical distribution barycentre located in tropical and subtropical regions [13].

During the last thirty years, Aedes have spread worldwide, recently becoming pests in several non-tropical countries [14]. In Europe, Italy is the most heavily infested country [14]. Here, the tiger mosquito Aedes albopictus (Skuse, 1894) has been recorded for the first time in 1990 [15] and is now well established [16, 17]. This species has been indicated as the primary vector for the first endemic outbreak of Chikungunya in Europe [18]. Furthermore, in France and Croatia Ae. albopictus has been blamed for the transmission of the first autochthonous dengue cases reported in Europe (in 2010 and 2013) in the last 80 years [19, 20]. However, in temperate countries, the distribution of Aedes species is limited by winter temperature [21, 22], and in Italy, Ae. albopictus is mainly present in areas below 600–800 m a.s.l. [23, 24]. In 2011, Aedes koreicus (Edwards, 1917), was found in Italy [25]. This species is native of South Korea, Japan, parts of China and ex-USSR countries [26] and was recorded in 2008 in Belgium for the first time outside its native range [27]. Aedes koreicus forms a monophyletic taxon with Aedes japonicus (Theobald, 1901), which is another emerging invasive mosquito in USA and Europe. Given its ecological plasticity, Ae. koreicus has been proposed as the next global invasive mosquito species [27, 28], with potential impact on human and animal health as the vector of Dirofilaria immitis, a heartworm, endemic in Northern Italy and the Japanese Encephalitis virus, mostly prevalent in Asia [26, 2931].

According to the few data available about its native distribution, Ae. koreicus may be able to tolerate lower winter temperatures than Ae. albopictus. It is also better adapted to urban environments than the forest dwelling species Ae. japonicus [26, 27]. An exploratory analysis performed using nine data points where Ae. koreicus was sampled in its native range (Korea [32, 33]) revealed that its native habitat has a yearly average temperature of 11.5 (sd 0.8) with the minimum average temperature of the coldest month of −9 (sd 1.7; Matteo Marcantonio, personal communication). With the highly competitive species Ae. albopictus well established in Italy, Ae. koreicus may presumably be outcompeted by Ae. albopictus in many areas with mild climate conditions (e.g., through larvae interspecific competition; [34]), but new populations of Ae. koreicus might establish in areas too cold for Ae. albopictus. Following this scenario, a wider geographical range could be colonized by Aedes mosquitoes, potentially widening the spatial distribution of Aedes-borne diseases. Therefore, describing Ae. koreicus potential distribution is critical for proacting ecological management able to promptly respond to the threat posed by this emerging invasive species [35].

The potential distribution of invasive species in a new geographical area can be assessed through mechanistic or correlative algorithms, generally referred to as invasive Species Distribution Models (iSDMs) [3639]. The main challenge with correlative iSDM is that, while many spatial modelling techniques require species to be at equilibrium with their environment, emerging invasive species are by definition in a dynamic transition state in the invaded range [9]. The equilibrium assumption may mislead predictions over broad areas since species’ capacities to colonize previously unoccupied areas may affect reliability of model prediction [40]. Therefore, the reliability of some ecological modelling techniques is disputed in the scientific literature, with hybrid (mechanistic together with correlative approach) and adaptive frameworks being more and more explored [39, 4143]. However, choosing methods is often dictated by more practical reasons such as availability of field and laboratory data, knowledge of species biology, project deadlines and statistical or mathematical expertise. Purely automatic statistical approaches require minimal knowledge about the species’ life history, ecology and physiology, making predictions easy to achieve. By contrast, other approaches make use of such knowledge to select parameters appropriate to model the physiological requirement of the investigated species or build mathematical representations of ecological processes. Such iSDMs require more effort but are more accurate when detailed physiological information are available or in circumstances that require in-depth understanding of survival and spread framework [41, 42].

As part of the LExEM (Laboratory of Excellence for Epidemiology and Modeling, project, we set a network of traps supported by larval searches in northern Italy. Using the collected data, we estimated the potential distribution of Ae. koreicus in Northeast Italy with three different iSDM approaches. First, we applied an automatic statistical approach, represented by the Maximum Entropy (MaxEnt) modelling [44]. Second, we implemented a logistic regression model with Bayesian framework informed by using prior knowledge retrieved from literature on the ecologically similar and better studied species Ae. albopictus. The information derived from field data is therefore mediated by a-priori ecological information [45, 46]. Third, we applied a Geographic Information System (GIS) physiology-based iSDM, solely relying on known environmental constraints of the species (e.g., the species cannot survive cold winter temperatures, etc.). Beyond describing species distributions, iSDMs have become an important and widely used decision making tool for a variety of applications, such as mapping risk of VBDs as well as their host spread, and determining locations that are potentially susceptible to invasion. Here, making use of multiple iSDMs, we aim to reliably estimate Ae. koreicus potential distribution in Northeast Italy, gathering insights into which iSDM should be preferred on the others. Our final goal was to investigate the future expansion of Ae. koreicus in the study area combining the developed habitat suitability map with available information about transportation networks and the observed species dispersal rate. Integrating the current and potential distribution of emerging invasive species with their preferred spread pathways is pivotal in identifying the most appropriate strategy to mitigate and control their invasion. In this paper, we hope to provide useful and validated spatial information about Ae. koreicus spread to decision makers in order to support control strategies and develop proactive public health policies.

Study area

The study area is located in Northeast Italy (Fig. 1; latitude N46.75, S45.59, longitude W10.38, E12.82; Datum WGS84). We investigated the presence of Ae. koreicus in two administrative units, Trentino and Veneto regions (EU NUTS2 code: ITD2 and ITD3). The area comprises the Eastern section of the Alps and the Northeastern part of the Po Valley. The altitude ranges from 0 to 3,900 m a.s.l. The climate ranges from subartic climate (Köppen climate classification: Dwc) in the Northern-most mountainous part to oceanic climate (Cfc) in the central and southern low altitude and flat part of the study area. The annual average temperature ranges from −6.8 to 15.3, the cumulative annual precipitation ranges from 439 to 1660 mm/year, while the total population count was around 1.7 million with a very variable density, spanning from 0 to 1100 people/km2 [47]. The study area includes one of the most developed agricultural, industrial and commercial areas of Italy [25] with dense connections between central and northern European commercial hubs. The area has already been invaded by Ae. albopictus which was first detected in 1991 to then gradually colonize all the suitable habitats.

Fig. 1
figure 1

Study areas: Right side: European map with the two red rectangles showing Ae. koreicus positive areas for Ae. koreicus in Italy (big rectangle) and in Belgium (small rectangle). For the statistical analysis, the Italian study area was used as the training area while the Belgian one as the test area. Left side: Zoom of the Italian study area showing trap locations and major cities, with the shaded digital elevation model as background


Ethical approval?

Data collection

Sampling design

We sampled a total of 394 locations from April to November in 2013 and 2014 using various collection methods such as larval searches, ovitraps, CDC-light traps and BG-sentinel traps (Biogents AG, Regenburg, Germany). We derived part of the data from entomological surveillance supported by Veneto Region. The sampling was implemented in order to acquire data on the current distribution of Ae. koreicus inside and at the boundary of its known invaded range. The sampled locations were checked every 2 weeks, except for larval search. Eggs collected from ovitraps were maintained in water for hatching, and larvae were reared until the fourth instar for identification [48, 49]. Adults caught in CDC-light traps and BG-sentinel traps were identified and stored at −80 for molecular analysis. Expert entomologists identified each sample at species-level and organized data in a geodatabase. A collection was considered positive (presence) if at least one individual was sampled during the sampling period. The presence/absence dataset was spatially aggregated in a 250 m × 250 m grid aligned to the EuroLST grid [50]. The grid cells containing both positive and negative locations were considered as positive. As a result, a total of 306 grid cells, of which 53 positive and 253 negative, were used in the following statistical analysis.

Environmental data sources and model predictors

To model the habitat suitability of Ae. koreicus, we used a set of environmental predictors based on a 10-years long average (2003–2012) derived from remote sensing (satellite) data at a spatial resolution of 250 m. Temperature variables have been derived from the EuroLST bioclim dataset, freely available at the EuroLST website ( Those bioclim variables which integrate temperature together with precipitation were calculated by merging the EuroLST [50] dataset with precipitation data from the Climate Prediction Center (CPC) Morphing algorithm (CMORPH) Version 1.0 [51] which we calibrated against data from the Global Precipitation Climatology Project (GPCP) Version 2.2 [52]. Moreover, two additional temperature predictors were derived; average temperature of the coldest month (TavgCM) and average temperature of mosquito growing season (TavgGS). The two latter variables were considered in order to match the biology of Aedes spp. with temperature variables. Indeed, in the study area, the time period from April to September represents the most favourable conditions for Aedes population growth. Furthermore, cold temperatures under 0 are a limiting factor for diapausing egg survival. On the one hand, in the literature it is reported that Aedes cold-acclimated and diapausing eggs can survive very low temperature (−10/-12) for a brief to moderate period of time [53, 54]. On the other hand, long periods of average cold weather represent a chronic stress which may more strongly limit the fitness (e.g., hatching) of diapausing eggs. Therefore, the average temperature of the coldest month may be more effective to limit mosquitoes potential distribution than the average minimum temperature of the coldest month.

We obtained data on the vegetation biomass from the MODIS Normalized Difference Vegetation Index (NDVI; MOD13Q1) product [55]. Vegetation indices, as NDVI, have been extensively used to describe disease risk and habitat suitability for different species of mosquitoes [5658].

Given that the ecosystem water content might also limit mosquitoes habitat suitability, as water is a key component of their ecological niche, we calculated the Normalized Difference Water Index (NDWI; [59]), derived from the MODIS surface reflectance product (MOD09A1). NDWI and NDVI embed different wavelengths, and so they should be considered as complementary [59]. We averaged NDVI and NDWI values pixel-wise in four seasonal groups over a three month period each (January-March, April-June, July-September, October-December), in order to have a representation of the vegetation coverage and ecosystem water for each of the four seasons.

The initial set of 29 environmental predictors (19 EuroLST/CMORPH bioclim, 2 further temperature-based variables, 4 seasonal NDVI and 4 seasonal NDWI; Table 1) was used in different combinations as input for the modelling approaches described in detail in the next section. All the considered environmental parameters have already been shown as relevant for mosquito iSDMs [60].

Table 1 Description of the predictor variables. We reported source and spatial resolution of each group of predictor variables

Modelling framework

SDM are a set of algorithms which quantitatively describe areas that support the presence of a given species, based on experimental data, known presence/absence data and the associated environmental conditions. These models seek, despite some limitations, the species ecological niche in the Hutchinsonian sense [61]. We made use of three different modelling approaches to estimate Ae. koreicus potential distribution in Northeast Italy. These three techniques rely on automatic statistical methods or on physiological knowledge of the species life history cycle. The target output of all these three modelling techniques was an environmental suitability indicator, expressed as a continuous value from 0 (no suitability) to 1 (complete suitability). We associated each suitability value to its respective EuroLST grid cell, at a resolution of 250 m. Therefore, we visualized the predicted environmental suitability (habitat suitability) in potential distribution maps.

Maximum Entropy (MaxEnt) approach to species distribution modelling

MaxEnt is a common modelling framework used in species distribution modelling [44, 62]. MaxEnt minimizes the relative entropy between the probability density of the predictors estimated from the presence data and the probability density of the predictors estimated from the region of interest (background information). This means that the geographic extent and the number of background samples influence the results. For each predictor, response curves can be generated describing how predictor values are related to the estimated suitability.

After having tested different buffer sizes, we have placed a buffer of 5 km around the presence data, representative of locations accessible for Ae. koreicus via dispersal and which approximate the overall study area environmental conditions [63, 64]. The resultant region has been used as input for MaxEnt modelling. We used MaxEnt as machine learning algorithm, letting it decide which predictors were important through regularization [65]. Therefore, we ran the MaxEnt model without any previous biologically-based selection of the predictor variables. However, since several predictors were highly correlated among them, we performed a correlation analysis in order to limit multicollinearity issues. We excluded all those predictors showing a correlation higher than 0.50 (Pearson’s r). The predictors exclusion was performed selecting the most correlated couple, followed by a random draw to decide what to exclude of the two predictors. We carried out all the analysis in R [66], using dismo [67] package.

Bayesian logistic regression (logBAY) with Markov Chain Monte Carlo simulation

Logit-link Generalized Linear Models (GLMs) are standard regression methods to model habitat suitability in ecology [68]. The presence or absence of a species is transformed in a probability function, real number in the range [0, 1], through a logistic transformation of the presence/absence odds (log(1/1  −  p)). We wrapped the logistic regression in a Bayesian framework using Just Another Gibbs Sampler (JAGS) [69] in combination with rjags [70] and coda [71] R packages. We used presence or absence data as response variable, while as predictor variables we chose those environmental variables with the strongest credibility in shaping Aedes ecological niche, as follows: i) average temperature of mosquito growing season TavgGS, ii) average annual temperature (TavgY), iii) average of the minimum temperature (TavgM); iv)TavgCM; v) cumulative annual precipitation (PcumY) and vi) spring NDWI (NDWIavgS) [14, 21, 22, 27, 28, 49, 54, 7277]. All predictors were scaled using their mean and standard deviation as follows: (x  −  mean(x))/sd(x), where x is the predictor variable. Even though the role of precipitation as limiting factor for container-breeding mosquitoes is controversial [73], we included it among predictors since water availability affects the aquatic stages of the mosquito life cycle [74]. Moreover a preliminary exploratory analysis showed a high correlation between precipitation and Ae. koreicus presence/absence for our dataset. Nevertheless, the observed correlation may be a spurious pattern linked with a different detectability probability in different parts of the precipitation range [78].

We used Gaussian distributed informed priors derived from [22] for the temperature based variables, whereas non-informative priors (normal distribution with mean = 0; precision = 10E  −  12) for all the other variables (Table 2). To select the combination of variables carrying the most information on mosquito distribution, we run all the models possible combining the six aforementioned predictor variables (including models with interactions between TavgGS or TavgY and PcumY and models with TavgY second order polynomial function). Each model was initialized using maximum likelihood estimates for each coefficient and 10000 burn-in Markov Chain Monte Carlo (MCMC) iterations to find a good starting point to sample a representative Posterior Probability Distributions (PPD). The models were therefore ranked using the Deviance Information Criterion (DIC) with 1,000 MCMC iterations and thinning set of 5. The most informative model (lowest DIC) was used to sample 10,000 times with thinning set of 50 the PPD of model parameters and of Ae. koreicus occurrence in each pixel of the study area. The convergence of MCMC chains was monitored using Gelman and Rubin’s convergence diagnostic between two MCMCs [79]. The PPD Highest Density Interval (HDI) was calculated using the function proposed in [80]. The average value of PPD was assigned to the correspondent pixel, resulting in the habitat suitability map for Ae. koreicus. Furthermore, the uncertainty linked to the average pixel prediction was assessed (and mapped; not reported) deriving the 95 % Bayesian confidence interval of the PPD of each pixel. We reported all the steps to reproduce the logBAY model as an R function in the Additional file 1.

Table 2 Average and precision for informed and non informed priors. The precision of a distribution is the inverse of its standard deviation

GIS physiology-based (PHY) suitability modelling

This iSDM approach considers environmental parameters corresponding to physiological constraints of Ae. koreicus. The exact physiological constraints for this species are currently unknown. We used a conservative and parsimonious approach by assuming that the same environmental parameters which represent limiting factors for Ae. albopictus can be applied, most importantly: the average temperature of the coldest month (TavgCM) and the average temperature of the hottest quarter of the year (TavgHQY). The temperature of the coldest month determines overwintering suitability: if the coldest month is under a certain threshold, diapausing eggs will not survive and a persistent population can not be established. If the temperature of the hottest quarter of the year does not reach a certain value, larvae can not develop and adults can not reproduce. Additionally, precipitation can determine habitat suitability, but needs to be treated with caution because irrigation and small anthropogenic water reservoirs can compensate for low precipitation. Other environmental parameters of potential importance are the average temperature of the mosquito growing season (TavgGS) and annual average temperature (TavgY). TavgY has previously been used to model habitat suitability for Ae. albopictus, but can not be linked to a particular physiological constraint. Suitable summer temperatures might be averaged out by cold winter temperatures, and equally cold winter temperatures might be averaged out by hot summer temperatures. The specific threshold for the three environmental parameters (average temperature of the coldest month, average temperature of the hottest quarter of the year and annual precipitation) was estimated from the values observed at sampling locations with presence of Ae. koreicus. We used the lower bound of the 99 % confidence interval as the low threshold for environmental parameters (Table 3). The estimated thresholds were used for habitat suitability modelling resulting in suitability maps. All temperature thresholds were transformed with a sigmoid function such that zero means not suitable, 0.5 corresponds to the actual threshold and 1 means highly suitable (Appendix C). A margin of 4 was applied to the sigmoid function for temperature. Moreover, for annual precipitation (PcumY), a margin of 200 mm/year was applied. Compared to MaxEnt, we defined a priori response curve for relevant environmental parameters, whereas such response curves are a diagnostic result of MaxEnt. The separate suitability indicators were multiplied in order to obtain a single general suitability index where 0 means that any single parameter was 0 (not suitable) and 1 represents that all single parameters were 1 (highly suitable).

Table 3 Predictor variables used in the GIS physiology-based suitability model. The descriptive statistics refers to the location of all the positive traps in the study area

Model performance accuracy

Model performance accuracy was measured assessing the error rate as a percentage (i.e. error rate (%) equals the number of incorrect cases divided by the total number of cases tested), as well as Cohen’s kappa coefficient (k), which is a measure of agreement that takes into account chance effects [81] and True Skill Statistics (TSS; [82]), an accuracy index not sensitive to prevalence. The optimal thresholds to discriminate the continuous model outputs in the presence or absence category were estimated by maximising sensitivity together with specificity.

We performed a further qualitative validation of the models sensitivity predicting Ae. koreicus overall average suitability and standard deviation in the only other known invaded area in Europe, the Maasmechelen municipality in Eastern Belgium (Fig. 1). In this locality, a viable hibernating Ae. koreicus population persists in a homogeneous 6 km2 industrial area since its first detection in 2008 [27]. When dealing with emerging invasive arthropods, absence points have a high likelihood to represent areas where the trap failed to catch entities of the species despite presence in the area, or areas that are inside the species ecological niche, but which have not yet been invaded (i.e., dispersal limitation). Therefore, we emphasized model sensitivity on specificity since it should be considered more effective to assess iSDM predictive power [8385].

Aedes koreicus spread analysis

To estimate the Ae. koreicus spread rate since introduction in the study area, we used information about its presence since 2011, when it was recorded for the first time in a small village near Belluno (Sospirolo village, latitude N46.14; longitude E12.07, datum: WGS84) [25]. In absence of information about the presence of the species before 2011, we assumed this geographical location as the centre of gravity for the introduction point(s).

In addition, we calculated the centre of gravity for the coordinates of positive traps for 2013 and 2014 respectively. We assumed that the range expansion has been constant through time, therefore we divided the Euclidean distance between 2011, 2013, 2014 centres of gravity by 4, representing the years since the introduction, deriving an approximate spreading rate, defined as spread distance in kilometres per year. Afterwards, we built a road network, weighted by the travelling distance between each road segment and the introduction point (root of the network, assumed to be the village Sospirolo). The road network was acquired from the OpenStreetMap project (, cleaned from tertiary roads, tracks and pathways which were assumed to be of low importance for mosquito dispersal. All the unconnected (not connected to the root of the network) road segments were also removed from the network.

The following step was to intersect the weighted road network with the habitat suitability map, to derive a distance-suitability weighted cost network. For the reference suitability map, we chose the one derived by the iSDM with the best predictive performance. This step was carried out to increase the cost for those locations that, despite being spatially close to the introduction location, were ecologically distant from Ae. koreicus ecological niche. We assumed that, for high suitability values (defined using the suitability threshold at which sensitivity plus specificity were maximized), the cost for the spread of the mosquito was the distance from the introduction point divided by the suitability value, while for suitability values below the threshold, the new weighted distance from the introduction location was the original distance from the introduction location divided by the suitability values raised to the power of 1.5 (penalty derived from empirical observations of the invasion process).

Eventually, we split the distance-suitability weighted road network into invasion cost isolines according to the observed Ae. koreicus dispersal rate. This step was performed in order to estimate areas with the same probability to be invaded in a defined temporal span (in years).

All the spatial analysis were performed using GRASS GIS 7 [86] modules (particularly tool set).


The potential distribution maps for Ae. koreicus derived from the three models are reported in Fig. 2.

Fig. 2
figure 2

Ae. koreicus potential distribution maps: The values range from 0: no suitability; to 1: complete suitability. The green triangles represent the centroids of the main cities in the area

MaxEnt modelling results

The output of the correlation analysis, which was a matrix with 9 predictors (Appendix A), was used as an input in the MaxEnt model.

Variable importance can be estimated by different means: percent contribution, permutation importance, and jacknife gain. These three measurements provide different rankings and are reported in Appendix A. According to the rank sum of the three different criteria, the two most important variables for the MaxEnt model were temperature seasonality (bio4) and maximum temperature of the warmest month (bio05). The response curve of bio4 suggests that Ae. koreicus presence probability is relatively constant until it drops suddenly in the localities where seasonality becomes extreme. Bio5 response curve implies a monotonic increase of presence probability in a temperature range between 22 and 28 °C, after which it assumes an asymptotic trend (Appendix A).

The Ae. koreicus potential distribution map derived from the MaxEnt model showed values ranging from 0 to 0.94, with an average suitability of 0.11 (Fig. 2a). According to the MaxEnt model, suitable areas are concentrated along the main Alpine valleys.

Bayesian logistic regression modelling results

The best logBAY model (lowest DIC; Table 4) comprised the average temperature of the growing season (April to September; TavgGS), the minimum temperature of the coldest month (TavgDEC) and the cumulative annual precipitation (PcumY). All the predictors were positively correlated with the presence of Ae. koreicus. The PPD of the model coefficients with their mean and 95 % HDI is showed for the best model in Appendix B. The 95 % HDI of TavgGS and PcumY did not include 0, meaning that the credible values of these model parameters are different than 0 (Appendix B).

Table 4 Model specifications and DIC for the best 15 logBAY models plus the full model

The logBAY suitability map, built using the average of the PPD, showed values ranging from 0 to 0.83, with an average suitability of 0.14 (Fig. 2b). The suitability predicted by logBAY model cut the study area in two distinct sections: high suitability in the southern part, low suitability in the northern mountainous area. The highest suitability was indicated for Pordenone, Treviso provinces and on the surroundings of lake Garda.

Physiology-based modelling results

The PHY suitability surface resembles the one of MaxEnt, with the difference of a much higher absolute suitability value. By definition, the PHY suitability value for most known presence sites is one. The main difference between the PHY, logBAY models and MaxEnt is due to the upper threshold in the MaxEnt response curves, which are mainly composed of parameters derived from temperature, whereas the other two models have not imposed an upper threshold on temperature (logBAY is a linear regression model, while for PHY there is no physiological restriction for maximum monthly temperature in the study area, which rarely reaches more than 35 °C).

The suitability values predicted by PHY ranged from 0 to 1, with an average suitability of 0.32. PHY had the maximum average suitability value among the developed models, with more than 65 % of the studied areas having values higher than 0.50 (Fig. 2c). Very high suitability was assigned to Po, Adige, Valsugana and Sarca valleys. Furthermore Piave, Isarco as well as other minor valleys were characterized by moderate to high suitability values.

Model validation

The suitability threshold at which sensitivity plus specificity were maximized for each model is reported in Table 5 together with Kappa statistics, TSS and percent error rate. The threshold at which sensitivity plus specificity were maximized varies considerably between the three models. After grouping the suitability values in suitable and not suitable classes using these thresholds, MaxEnt reported the highest error rate. On the contrary, logBAY and PHY had TSS and Kappa values indicating from substantial to almost perfect agreement with observed data [87]. PHY showed the highest sensitivity while logBAY the highest specificity.

Table 5 Model performance accuracy. We reported the suitability thresholds at which sensitivity plus specificity were maximized, Kappa statistics, TSS and error rate for each of the model

Applying the discriminant thresholds listed in Table 5 to the suitability maps, we found that 3 %, 26 % and 30 % of the study area was reported as suitable by MaxEnt, logBay and PHY models respectively.

The result of the cross-tabulation between elevation and suitable area is reported in Fig. 3. All the models agreed on the suitable area being concentrated at low altitude (0–800 m). Half of the total area between 0 and 800 m was indicated as suitable according to logBAY and PHY. It is interesting to note that all the profiles in Fig. 3 show a spike in suitability around 400–500 m. Moreover, PHY model predicted as suitable a remarkable percentage (16 %) in higher altitude area (above 800 m).

Fig. 3
figure 3

Altitude profile of suitable area: This figure depicts the percentage of suitable area over the total area for each altitude class for a) PHY; b) logBAY and c) MaxEnt model. The black line represents the percentage of area in each corresponding altitude class

To further validate the model using an independent set of data, we applied each model to the only other area invaded by Ae. koreicus in Europe: Maasmechelen municipality in Belgium (Fig. 1). We reported the descriptive statistics of the predicted suitability distribution in Table 6. PHY model predicted high suitability, whereas logBAY very low suitability.

Table 6 Descriptive statistics for the distribution of suitability values in Maasmechelen municipality, Belgium, for all three models

Aedes koreicus spread analysis

From 2011 to 2014, the average shift of the invaded area centroid was approximately 8 km/year. The analysis performed to predict the spread of Ae. koreicus showed that the most likely dispersal direction was along Valsugana Valley (Fig. 4). Furthermore, the mosquito may be already present in the northern part of Verona province (south-east of the study area; predicted to be invaded within 5–10 years after introduction). According to our results the mosquito spread might be rather fast in the southern part of the study area due to both the dense road connection with the introduction point and to the high habitat suitability. Even though the dispersal along the highly populated Adige Valley was found to be slower than in the southern part of the study area, probably due to its more rough topography, overall it will potentially be at high risk of invasion during the forthcoming decades.

Fig. 4
figure 4

Potential spread of Aedes koreicus predicted through road network analysis: Areas with the same cost of invasion are displayed using a red-green-blue colour scale. The cost of invasion is expressed in years since the species’ introduction (2011). Cost of invasion is a function of the travelling distance from the introduction point based on the observed rate of shift of the invaded range centroid and the predicted habitat suitability. Major cities (green pushpins) and sampling locations (white circles) are also reported


Aedes koreicus is an emerging invasive species in Europe, and is a nuisance and potential vector of infectious diseases [26, 31]. In this study, we assessed its habitat suitability in Northeast Italy, making use of iSDMs and a limited amount of field data supported by prior knowledge from related species. We also investigated the potential pathways and timing of future spread through a road network analysis. The main outcome was that the known distribution of Ae. koreicus is only a fraction of the potentially suitable area. However, we observed a rather variable characterization of the suitable area, highly dependent on the considered iSDM. MaxEnt predicted the smallest suitable area, while PHY the largest. The observed differences are mainly due to the lack of upper temperature threshold limits for logBAY and PHY (Appendix A, C), which allows high suitability in the warmest section of the study area. Despite potential limitations due to the approximation of biological patterns with linear relations in ecological niche modelling, the lack of upper thresholds for temperature should not to be considered an artefact due to the studied species and the climatic characteristic of the study area. Indeed, Ae. koreicusis a temperate-continental taxon and the studied area is characterized by monthly average temperatures rarely exceeding 30 °C (therefore well inside the temperature niche of other Aedes species; [88]). However, high temperatures might affect local abundance of Ae. koreicus (e.g., [89] observed a decrease of Ae. japonicus larvae survivorship at temperatures over 22 °C under laboratory conditions).

All the models were mainly driven by temperature variables, confirming what has already been found in literature for other Aedes species (e.g., [90]). The most important predictors for MaxEnt were bio4 and bio5, the best logBAY model included two temperature variables TavgGS, TavgDEC in addition to PcumY, while PHY was completely constrained by temperature (the included precipitation threshold was always exceeded in the study area). Mosquitoes are small-bodied poikilotherms, meaning that ambient temperature is the main abiotic factor limiting their ecological niche and, therefore, their geographical distribution [77, 91]. High temperature decreases embryonic (e.g., [92]) and larval (e.g., [93]) development time, and the size of adults (e.g., [94]), while cold winter temperatures have a severe impact on the survival of diapausing eggs [54]. The lesser importance of environmental indicators other than temperature can be explained by the considered spatial scale (i.e., extension and grain), where vegetation variability might be less important than climatic conditions [95], and by the autoecology of Aedes mosquitoes, container-breeding species, able to develop independently of the regional precipitation trend and environmental variability. At a finer spatial scale and in urban habitats, vegetation may be more influential on Aedes life cycle. Indeed, in this setting, even small pockets of vegetation favour habitat heterogeneity, allowing mosquitoes to regulate extreme weather conditions [96].

We performed model validation with a two-step analysis, i.e. a classical accuracy assessment with a dependent set of data and a qualitative sensitivity analysis using an independent set of data. The classical accuracy analysis suggested logBAY as the best model. However, PHY indicated a moderate-to-high agreement with observed data, while MaxEnt performed poorly. All the three models predicted presence locations with high accuracy (high sensitivity). On the contrary, model specificity was relatively low for all three models. As extensively reported in the literature, absence data is the Achille’s heel of SDM due to the high uncertainty linked to absence data veracity (see [9799]). This is especially true when dealing with emerging invasive arthropods, whose true absence is hard to identify (see [84, 85]). Therefore, to further assess the model accuracy, we decided to perform a further sensitivity analysis. MaxEnt and logBAY predicted low suitability for Maasmechelen municipality, despite this, the area has hosted a viable Ae. koreicus population at least since 2008 [27]. By contrast, PHY model predicted a moderate to high average suitability. This model, based on the construction of mechanistic overlay functions for climatic constraints, is partially independent from local datasets and thus tends to be more accurate for prediction on an independent dataset.

The percentage of the study area predicted as suitable varied from 3 to 30 %, encompassing different topographic and environmental conditions. To better characterize the predicted suitable area, we cross-tabulated it with a digital elevation model. PHY and logBAY indicated most of the low-altitude areas as highly suitable, while MaxEnt showed a peak in suitability distributed around moderate altitude (400 m). This outcome may be due to the buffer size (5 km) chosen around the presence points to derive background data, which may over-represent hilly areas, influencing the MaxEnt output. However, all three models showed a peak in suitability around 400–500 m, which may indicate optimal ecological conditions for Ae. koreicus. In support of this hypothesis we noted that the trap with the highest Ae. koreicus abundance was located at an altitude of 451 m. Another interesting outcome is that PHY predicted as suitable areas between 600 and 1500 m. This altitudinal range represents a still empty niche for invasive Aedes, as 600 m is the altitudinal limit for Ae. albopictus distribution [24] in Northeast Italy. As a result, the area between 600 and 1500 m suitable for Ae. koreicus should be particularly monitored as, here, the invasion would not be constrained by biotic interactions with species with similar evolutionary traits.

The accuracy assessment indicates PHY as the model with the highest prediction power, being in moderate-to-high agreement with observed data in the study area and predicting high suitability in a positive location with a different environmental setting. Therefore, PHY model was chosen as reference to investigate how Ae. koreicus may further spread in Northeast Italy. The spread analysis was achieved considering the observed dispersal rate since introduction, the preferred dispersal pathways, study area connectivity and habitat suitability. Aedes species have a short flight range, with a flying dispersal capability of 200–300 m radius per week around the hatching location [100, 101]. However, the short active dispersal range is generally compensated for by long distance used-tyres transportation and the plant nursery trade (Dracaena sp.) in the form of drought resistant eggs [102]. The local dispersal in new invaded areas is also boosted by humans, through the movement of garden waste, moist vegetation and water containers that can hold eggs or larvae as well as dispersal in trucks transporting used tyres or private vehicle [102, 103]. As a consequence, it can be inferred that the local dispersal probability in a newly invaded area is a function of the introduction point, local transportation network as well as habitat suitability. From these premises, we derived a distance-suitability weighted road network to predict which areas in Northeast Italy have the highest probability to be invaded. The results revealed how the centroid of the invaded range has been shifting approximately 8 km/year since 2011 (putative introduction year). Assuming a constant invaded range shift and driving it along the shortest road pathway (lowest cost from the introduction point), weighted according to the suitability of each road segment, we built a potential dispersal network which represents a reliable dynamic description of the invaded area evolution in the next decades. The simulated spread predicted all the known presence locations (except one) as invaded in a time frame of 5–10 years since its introduction. Moreover, it showed how the species may have already invaded the two major cities in the southern part of the study area, Treviso and Pordenone. However, a first investigation in July 2015 did not find positive locations in these cities [104]. Furthermore, the simulated spread predicted the north part of Po Valley and the southern Adige Valley as invaded in the next decade. A favourable topography (continuous flat areas), mild climate and dense and congested road network underlie the predicted rapid spread in these parts of the study area. On the contrary, we noticed no predicted spread in the north side of the study area, apart from limited spots such as the southern Isarco and northern Adige Valleys, where temperature hotspots due to towns (Urban Heat Island) as well as high road connectivity may favour Ae. koreicus spread over the next years.

The simulated spread is a reliable approximation of the future expansion of Ae. koreicus distribution range since it integrates a validated suitability surface as well as the most likely dispersal pathways at local scale. A partial validation of the adopted spread analysis comes from a similar study on Ae. albopictus by [93]. The authors found that Ae. albopictus is currently surfing a dispersal wave in Southern France, with occasional “jumps” that did not result in new colonization fronts. However, it should be remembered that the proposed approach may underestimate the dispersal rate. This is due to the choice to consider the centroids as indication of the invaded area shift as well as the deterministic nature of our approach which does not integrate stochastic events such as occasional introductions in spatially distant but ecologically close locations. Stochasticity in species distribution change underlies unpredictable events that sometimes strongly boost species dispersal and colonization of new areas. Besides, in the southern part of the study area, biotic interactions and out-competition by Ae. albopictus, not considered in this study, may slow down Ae. koreicus spread [34], modifying the outcome of the invasion process. Preliminary larval competition experiments suggested that the larval development of Ae. koreicus might be negatively affected by the presence of Ae. albopictus (Frédéric Baldacchino, personal communication).


Despite the rising concern about biological invasions after recent economic and human health issues due to invasive species (e.g., Drosophila suzukii (Matsumura, 1931) and Xylella fastidiosa (Wells et al. 1986) as crop pests and Ae. albopictus as a vector of tropical pathogens; [18, 105, 106]), at present there is no coordinated plan which aims to manage Ae. koreicus in the study area. Multiple control and mitigation strategies are available to eradicate, mitigate or control invasive species [107]. At the beginning of an invasion, as is the case of Ae. koreicus, the most effective control strategy is through inspections followed by destruction of removable breeding sites (e.g., plastic drums) and treatment with larvicidae of fixed sites (e.g., concrete bins). This strategy is time consuming and might be improved in terms of cost-effectiveness by targeting the most productive breeding sites. However, there is often a limited understanding of the biology of emerging invasive species and, consequently, of the hazard they represent [108]. Delays in early mitigation actions result in escalating costs of control, reduced economic returns from management actions and decreased feasibility of management [35, 109, 110]. iSDMs and spread pathway analysis are powerful tools to shed light on the present and future invader distribution and to inform on-ground control of the invasions [111]. We suggest that modelling and mapping the spatial distribution of invasive mosquitoes, validated by entomological surveys, should routinely support the implementation of control actions to limit their expansion. We hope that the results in this study serve as a foundation for design policies aiming to limit Ae. koreicus invasion in Northeast Italy.


  1. Wilson JRU, Dormontt EE, Prentis PJ, Lowe AJ, Richardson DM. Something in the way you move: dispersal pathways affect invasion success. Trends Ecol Evol (Amst). 2009;24:136–44.

    Article  Google Scholar 

  2. Pysek P, Jarošík V, Hulme PE, Kühn I, Wild J, Arianoutsou M, et al. Disentangling the role of environmental and human pressures on biological invasions across Europe. PNAS. 2010;107:12157–62.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Stenseth NC, Ottersen G, Hurrell JW, Mysterud A, Lima M, Chan K-S, et al. Studying climate effects on ecology through the use of climate indices: The North Atlantic Oscillation, El Niño Southern Oscillation and beyond. Proc R Soc B Biol Sci. 2003;270(1529):2087–96.

    Article  Google Scholar 

  4. Foley JA, DeFries R, Asner GP, Barford C, Bonan G, Carpenter SR, et al. Global consequences of land use. Science. 2005;309:570–4.

    Article  CAS  PubMed  Google Scholar 

  5. Hobbs RJ, Arico S, Aronson J, Baron JS, Bridgewater P, Cramer VA, et al. Novel ecosystems: theoretical and management aspects of the new ecological world order. Glob Ecol Biogeogr. 2006;15:1–7.

    Article  Google Scholar 

  6. Hobbs RJ, Higgs ES, Harris JA. Novel ecosystems: concept or inconvenient reality? A response to Murcia et al. Trends Ecol Evol. 2014;29:645–6.

    Article  PubMed  Google Scholar 

  7. Lockwood JL, Cassey P, Blackburn T. The role of propagule pressure in explaining species invasions. Trends Ecol Evol. 2005;20:223–8 [Special Issue: Invasions, Guest Edited by Michael E. Hochberg and Nicholas J. Gotelli].

    Article  PubMed  Google Scholar 

  8. Lenoir J, Svenning J-C. Climate-related range shifts – a global multidimensional synthesis and new research directions. Ecography. 2014;38:15–28.

    Article  Google Scholar 

  9. Crowl TA, Crist TO, Parmenter RR, Belovsky G, Lugo AE. The spread of invasive species and infectious disease as drivers of ecosystem change. Front Ecol Environ. 2008;6:238–46.

    Article  Google Scholar 

  10. Vora N. Impact of anthropogenic environmental alterations on vector-borne diseases. Medscape J Med. 2008;10:238.

    PubMed  PubMed Central  Google Scholar 

  11. Buck C. The Challenge of Epidemiology: Issues and Selected Readings. Washington, D.C., USA: Pan American Health Org; 1988. p. 989.

  12. Lehane M: The Biology of Blood-Sucking in Insects. Cambridge University Press; 2005.

  13. Reinert JF, Harbach RE, Kitching IJ. Phylogeny and classification of Aedini (Diptera: Culicidae), based on morphological characters of all life stages. Zool J Linnean Soc. 2004;142:289–368.

    Article  Google Scholar 

  14. Medlock JM, Hansford KM, Versteirt V, Cull B, Kampen H, Fontenille D, et al. An entomological review of invasive mosquitoes in Europe. Bull Entomol Res 2015;105(6):1–27.

    Article  Google Scholar 

  15. Sabatini A, Raineri V, Trovato G, Coluzzi M. Aedes albopictus in Italy and possible diffusion of the species into the Mediterranean area. Parassitologia. 1990;32:301–4.

    CAS  PubMed  Google Scholar 

  16. Carrieri M, Bacchi M, Bellini R, Maini S. On the competition occurring between Aedes albopictus and Culex pipiens (diptera: culicidae) in Italy. Environ Entomol. 2003;32:1313–21.

    Article  Google Scholar 

  17. Bellini R, Veronesi R, Venturelli G, Angelini P: Guidelines for Surveillance and Control of the Asian Tiger Mosquito (Aedes Albopictus). Regione Emilia-Romagna; 2005;24:1–85.

    Google Scholar 

  18. Rezza G, Nicoletti L, Angelini R, Romi R, Finarelli AC, Panning M, et al. Infection with chikungunya virus in Italy: an outbreak in a temperate region. Lancet. 2007;370:1840–6.

    Article  CAS  PubMed  Google Scholar 

  19. La Ruche G, Souarès Y, Armengaud A, Peloux-Petiot F, Delaunay P, Desprès P, et al. First two autochthonous dengue virus infections in metropolitan France, September 2010. Euro Surveill. 2010;15:19676.

    PubMed  Google Scholar 

  20. Marchand E, Prat C, Jeannin C, Lafont E, Bergmann T, Flusin O, et al. Autochthonous case of dengue in France, October 2013. Euro Surveill. 2013;18:20661.

    Article  CAS  PubMed  Google Scholar 

  21. Neteler M, Metz M. Is Switzerland suitable for the invasion of Aedes albopictus? PLoS ONE. 2013;8:e82090.

    Article  PubMed  PubMed Central  Google Scholar 

  22. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. ECDC: Development of Aedes albopictus risk maps. 2012. Stockholm: European Centre for Disease Prevention and Control, 2009.

  24. Valerio L, Marini F, Bongiorno G, Facchinelli L, Pombi M, Caputo B, et al. Host-feeding patterns of Aedes albopictus (Diptera: Culicidae) in urban and rural contexts within Rome province, Italy. Vector Borne Zoonotic Dis. 2010;10:291–4.

    Article  PubMed  Google Scholar 

  25. Capelli G, Drago A, Martini S, Montarsi F, Soppelsa M, Delai N, et al. First report in italy of the exotic mosquito species Aedes (Finlaya) koreicus, a potential vector of arboviruses and filariae. Parasit Vectors. 2011;4:188.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Tanaka T, Mizusawa K, Saugstad E: Mosquitoes of Japan and Korea. Volume 16. Contributions of the American Entomological Institute; 1979.

  27. Versteirt V, De Clercq EM, Fonseca DM, Pecor J, Schaffner F, Coosemans M, et al. Bionomics of the established exotic mosquito species Aedes koreicus in Belgium, Europe. J Med Entomol. 2012;49:1226–32.

    Article  CAS  PubMed  Google Scholar 

  28. Kaufman MG, Fonseca DM. Invasion biology of Aedes japonicus japonicus (Diptera: Culicidae). Annu Rev Entomol. 2014;59:31–49.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Feng L. Experiments with dirofilaria immitis and local species of mosquitoes in peiping, north china. Ann Trop Med Parasit. 1930;24:347–66.

    Article  Google Scholar 

  30. Miles JAR. Some ecological aspects of the problem of arthropod-borne animal viruses in the Western Pacific and South-East Asia regions. Bull World Health Organ. 1964;30:197–210.

    CAS  PubMed  PubMed Central  Google Scholar 

  31. Montarsi F, Ciocchetta S, Devine G, Ravagnan S, Mutinelli F, Frangipane Di Regalbono A, et al. Development of Dirofilaria immitis within the mosquito Aedes (Finlaya) koreicus, a new invasive species for Europe. Parasit Vectors. 2015;8:177.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Kim HC, Chong ST, O’brien LL, O’guinn ML, Turell MJ, Lee H-C, et al. Seasonal prevalence of mosquitoes collected from light traps in the Republic of Korea in 2003. Entomol Res. 2006;36:139–48.

    Article  Google Scholar 

  33. Kim HC, Chong ST, Nunn PV, Klein TA. Seasonal prevalence of mosquitoes collected from light traps in the Republic of Korea, 2007. Entomol Res. 2010;40:136–44.

    Article  Google Scholar 

  34. Juliano SA, Lounibos LP, O’Meara GF. A field test for competitive effects of Aedes albopictus on A. aegypti in South Florida: differences between sites of coexistence and exclusion? Oecologia. 2004;139:583–93.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Dunn AM, Hatcher MJ. Parasites and biological invasions: Parallels, interactions, and control. Trends Parasitol 2015;31(5):189–99.

    Article  PubMed  Google Scholar 

  36. Guisan A, Zimmermann NE. Predictive habitat distribution models in ecology. Ecol Model. 2000;135:147–86.

    Article  Google Scholar 

  37. Guisan A, Thuiller W. Predicting species distribution: Offering more than simple habitat models. Ecol Lett 2005;8(9):993–1009.

    Article  Google Scholar 

  38. Elith J, Leathwick JR. Species distribution models: ecological explanation and prediction across space and time. Supplemental Material: Annu Rev Ecol Evol Syst. 2009;40:1–20.

    Google Scholar 

  39. Kearney M, Porter W. Mechanistic niche modelling: combining physiological and spatial data to predict species’ ranges. Ecol Lett. 2009;12:334–50.

    Article  PubMed  Google Scholar 

  40. Leroux SJ, Larrivée M, Boucher-Lalonde V, Hurford A, Zuloaga J, Kerr JT, et al. Mechanistic models for the spatial spread of species under climate change. Ecol Appl. 2013;23:815–28.

    Article  PubMed  Google Scholar 

  41. Dormann CF, Schymanski SJ, Cabral J, Chuine I, Graham C, Hartig F, et al. Correlation and process in species distribution models: bridging a dichotomy. J Biogeogr. 2012;39:2119–31.

    Article  Google Scholar 

  42. Uden DR, Allen CR, Angeler DG, Corral L, Fricke KA. Adaptive invasive species distribution models: A framework for modeling incipient invasions. Biol Invasions 2015;17(10):1–20.

    Article  Google Scholar 

  43. Barker C, Montecino D, Marcantonio M. Modeling the spread and control of the Asian tiger mosquito in Los Angeles. San Francisco: American Geophysical Union Fall meeting; 2015. p. 14–8.

    Google Scholar 

  44. Phillips S, Anderson R, Schapire R. Maximum entropy modeling of species geographic distributions. Ecol Model. 2006;190:231–59.

    Article  Google Scholar 

  45. McCarthy MA, Masters P. Profiting from prior information in Bayesian analyses of ecological data. J Appl Ecol. 2005;42:1012–9.

    Article  Google Scholar 

  46. Hooten MB, Wikle CK. A hierarchical Bayesian non-linear spatio-temporal model for the spread of invasive species with application to the Eurasian Collared-Dove. Environ Ecol Stat. 2007;15:59–70.

    Article  Google Scholar 

  47. CIESIN - Columbia University: Gridded Population of the World, Version 4 (GPWv4): Population Density Grid, Future Estimates. 2015.

  48. Becker N. Mosquitoes and their control. Berlin, London: Springer; 2010.

    Book  Google Scholar 

  49. Montarsi F, Martini S, Dal Pont M, Delai N, Ferro Milone N, Mazzucato M, et al. Distribution and habitat characterization of the recently introduced invasive mosquito Aedes koreicus [Hulecoeteomyia koreica], a new potential vector and pest in north-eastern Italy. Parasit Vectors. 2013;6:292.

    Article  PubMed  PubMed Central  Google Scholar 

  50. Metz M, Rocchini D, Neteler M. Surface temperatures at the continental scale: tracking changes with remote sensing at unprecedented detail. Remote Sens. 2014;6:3822–40.

    Article  Google Scholar 

  51. Joyce RJ, Janowiak JE, Arkin PA, Xie P. CMORPH: a method that produces global precipitation estimates from passive microwave and infrared data at high spatial and temporal resolution. J Hydrometeorol. 2004;5:487–503.

    Article  Google Scholar 

  52. Adler RF, Huffman GJ, Chang A, Ferraro R, Xie P-P, Janowiak J, et al. The Version-2 Global Precipitation Climatology Project (GPCP) Monthly Precipitation Analysis (1979–Present). J Hydrometeorol. 2003;4:1147–67.

    Article  Google Scholar 

  53. Hanson SM, Craig GB. Relationship between cold hardiness and supercooling point in Aedes albopictus eggs. J Am Mosq Control Assoc. 1995;11:35–8.

    PubMed  Google Scholar 

  54. Thomas SM, Obermayr U, Fischer D, Kreyling J, Beierkuhnlein C. Low-temperature threshold for egg survival of a post-diapause and non-diapause European aedine strain, Aedes albopictus (Diptera: Culicidae). Parasit Vectors. 2012;5:100.

    Article  PubMed  PubMed Central  Google Scholar 

  55. Neteler M, Metz M. MODIS and vector-borne diseases. IEEE Earthzine; 2014.

  56. Brown H, Duik-Wasser M, Andreadis T, Fish D. Remotely-sensed vegetation indices identify mosquito clusters of West Nile virus vectors in an urban landscape in the northeastern United States. Vector Borne Zoonotic Dis. 2008;8:197–206.

    Article  PubMed  Google Scholar 

  57. Kleinschmidt I, Bagayoko M, Clarke GP, Craig M, Le Sueur D. A spatial statistical approach to malaria mapping. Int J Epidemiol. 2000;29:355–61.

    Article  CAS  PubMed  Google Scholar 

  58. Lourenço PM, Sousa CA, Seixas J, Lopes P, Novo MT, Almeida APG. Anopheles atroparvus density modeling using MODIS NDVI in a former malarious area in Portugal. J Vector Ecol. 2011;36:279–91.

    Article  PubMed  Google Scholar 

  59. Gao BC. NDWI - A normalized difference water index for remote sensing of vegetation liquid water from space. Remote Sens Environ. 1996;58:257–66.

    Article  Google Scholar 

  60. Eisen L, Eisen RJ. Using geographic information systems and decision support systems for the prediction, prevention, and control of vector-borne diseases. Annu Rev Entomol. 2011;56:41–61.

    Article  CAS  PubMed  Google Scholar 

  61. Hutchinson GE. Concluding Remarks. Cold Spring Harb Symp Quant Biol. 1957;22:415–27.

    Article  Google Scholar 

  62. Phillips SJ, Dudı M. Modeling of species distributions with Maxent: New extensions and a comprehensive evaluation. Ecography. 2008;31:161–75.

    Article  Google Scholar 

  63. Saupe EE, Barve V, Myers CE, Soberón J, Barve N, Hensz CM, et al. Variation in niche and distribution model performance: The need for a priori assessment of key causal factors. Ecol Model. 2012;237–238:11–22.

    Article  Google Scholar 

  64. Merow C, Smith MJ, Edwards TC, Guisan A, McMahon SM, Normand S, et al. What do we gain from simplicity versus complexity in species distribution models? Ecography. 2014;37:1267–81.

    Article  Google Scholar 

  65. Elith J, Phillips SJ, Hastie T, Dudík M, Chee YE, Yates CJ. A statistical explanation of MaxEnt for ecologists. Divers Distrib. 2011;17:43–57.

    Article  Google Scholar 

  66. R Core Team: R: A Language and Environment for Statistical Computing [Internet]. Vienna, Austria: R Foundation for Statistical Computing; 2015. Available from:

  67. Hijmans RJ, Phillips S, Leathwick J, Elith J. dismo: Species Distribution Modeling [Internet]. 2015. Available from:

  68. Elith J, Graham C, Anderson R, Dudík M, Ferrier S, Guisan A, et al. Novel methods improve prediction of species’ distributions from occurrence data. Ecography. 2006;29:129–51.

    Article  Google Scholar 

  69. Plummer M. JAGS: A Program for Analysis of Bayesian Graphical Models Using Gibbs Sampling. 2003. In: Proceedings of the 3rd International Workshop on Distributed Statistical Computing, March 20-22, 2003. Vienna, Austria: Technische Universität Wien.

  70. Plummer M. rjags: Bayesian Graphical Models using MCMC [Internet]. 2015. Available from:

  71. Plummer M, Best N, Cowles K, Vines K. CODA: Convergence Diagnosis and Output Analysis for MCMC. R News. 2006;6(1):7–11.

    Google Scholar 

  72. Fischer D, Thomas S, Neteler M, Tjaden N, Beierkuhnlein C. Climatic suitability of Aedes albopictus in Europe referring to climate change projections: Comparison of mechanistic and correlative niche modelling approaches. Eurosurveillance. 2014;19(6).

  73. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  74. Caminade C, Medlock JM, Ducheyne E, McIntyre KM, Leach S, Baylis M, et al. Suitability of European climate for the Asian tiger mosquito Aedes albopictus: recent trends and future scenarios. J R Soc Interface. 2012;9:2708–17.

    Article  PubMed  PubMed Central  Google Scholar 

  75. Neteler M, Roiz D, Rocchini D, Castellani C, Rizzoli A. Terra and aqua satellites track tiger mosquito invasion: modelling the potential distribution of Aedes albopictus in north-eastern Italy. Int J Health Geogr. 2011;10:49.

    Article  PubMed  PubMed Central  Google Scholar 

  76. 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.

    Article  PubMed  Google Scholar 

  77. Alto BW, Juliano SA. Temperature effects on the dynamics of Aedes albopictus (Diptera: Culicidae) populations in the laboratory. J Med Entomol. 2001;38:548–56.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  78. Krishna YC, Krishnaswamy J, Kumar NS. Habitat factors affecting site occupancy and relative abundance of four-horned antelope. J Zool. 2008;276:63–70.

    Article  Google Scholar 

  79. Gelman A, Rubin D. Inference from iterative simulation using multiple sequences. Stat Sci. 1992;7:457–511.

    Article  Google Scholar 

  80. Kruschke JK. Doing Bayesian Data Analysis. 1st Edition. Academic Press; 2010.

  81. Manel S, Ceri Williams H, Ormerod S. Evaluating presence-absence models in ecology: The need to account for prevalence. J Appl Ecol. 2001;38:921–31.

    Article  Google Scholar 

  82. Allouche O, Tsoar A, Kadmon R. Assessing the accuracy of species distribution models: Prevalence, kappa and the true skill statistic (TSS). J Appl Ecol. 2006;43(6):1223–32.

    Article  Google Scholar 

  83. Ward DF. Modelling the potential geographic distribution of invasive ant species in New Zealand. Biol Invasions. 2006;9:723–35.

    Article  Google Scholar 

  84. Cianci D, Hartemink N, Ibáñez-Justicia A. Modelling the potential spatial distribution of mosquito species using three different techniques. Int J Health Geogr. 2015;14:10.

    Article  PubMed  PubMed Central  Google Scholar 

  85. De Clercq EM, Leta S, Estrada-Peña A, Madder M, Adehan S, Vanwambeke SO. Species distribution modelling for Rhipicephalus microplus (Acari: Ixodidae) in Benin, West Africa: Comparing datasets and modelling algorithms. Prev Vet Med. 2015;118:8–21.

    Article  PubMed  Google Scholar 

  86. Neteler M, Bowman MH, Landa M, Metz M. GRASS GIS: A multi-purpose open source GIS. Environ Model Softw. 2012;31:124–30.

    Article  Google Scholar 

  87. Landis JR, Koch GG. The measurement of observer agreement for categorical data. Biometrics. 1977;33:159–74.

    Article  CAS  PubMed  Google Scholar 

  88. Smith G, Eliason D. Use of elevated temperatures to kill Aedes albopictus and Ae. aegypti. J Am Mosq Control Assoc. 1989;4:557–8.

    Google Scholar 

  89. Scott J. The Ecology of the Exotic Mosquito Ochlerotatus (Finlaya) japonicus japonicus (Theobald 1901) (Diptera: Culicidae) and an Examination of its Role in the West Nile Virus Cycle in New Jersey. Ph.D. thesis. Rutgers University; 2003.

  90. Roche B, Léger L, L'Ambert G, Lacour G, Foussadier R, Besnard G, et al. The spread of Aedes albopictus in metropolitan France: contribution of environmental drivers and human activities and predictions for a near future. PLoS ONE. 2015;10:e0125600.

    Article  PubMed  PubMed Central  Google Scholar 

  91. Brady OJ, Johansson MA, Guerra CA, Bhatt S, Golding N, Pigott DM, et al. Modelling adult Aedes aegypti and Aedes albopictus survival at different temperatures in laboratory and field settings. Parasit Vectors. 2013;6:351.

    Article  PubMed  PubMed Central  Google Scholar 

  92. Trpis M, Haufe WO, Shemanchuk JA. Embryonic development of Aedes (o.) sticticus (diptera: culicidae) in relation to different constant temperatures. Can Entomol. 1973;105:43–50.

    Article  Google Scholar 

  93. Teng H-J, Apperson CS. Development and survival of immature Aedes albopictus and Aedes triseriatus (Diptera: Culicidae) in the laboratory: effects of density, food, and competition on response to temperature. J Med Entomol. 2000;37:40–52.

    Article  CAS  PubMed  Google Scholar 

  94. Rueda LM, Patel KJ, Axtell RC, Stinner RE. Temperature-dependent development and survival rates of Culex quinquefasciatus and Aedes aegypti (Diptera: Culicidae). J Med Entomol. 1990;27:892–8.

    Article  CAS  PubMed  Google Scholar 

  95. Austin MP. Spatial prediction of species distribution: An interface between ecological theory and statistical modelling. Ecol Model. 2002;157(2-3):101–18.

    Article  Google Scholar 

  96. Meyer RP, Hardy JL, Reisen WK. Diel changes in adult mosquito microhabitat temperatures and their relationship to the extrinsic incubation of arboviruses in mosquitoes in Kern County, California. J Med Entomol. 1990;27:607–14.

    Article  CAS  PubMed  Google Scholar 

  97. Hirzel AH, Hausser J, Chessel D, Perrin N. Ecological-niche factor analysis: How to compute habitat-suitability maps without absence data? Ecology. 2002;83:2027–36.

    Article  Google Scholar 

  98. Brotons L, Thuiller W, Araújo MB, Hirzel AH. Presence-absence versus presence-only modelling methods for predicting bird habitat suitability. Ecography. 2004;27:437–48.

    Article  Google Scholar 

  99. Drake JM. Range bagging: a new method for ecological niche modelling from presence-only data. J R Soc Interface. 2015;12:20150086.

    Article  PubMed  PubMed Central  Google Scholar 

  100. Turell MJ, Dohm DJ, Sardelis MR, Oguinn ML, Andreadis TG, Blow JA. An update on the potential of north American mosquitoes (Diptera: Culicidae) to transmit West Nile Virus. J Med Entomol. 2005;42:57–62.

    Article  PubMed  Google Scholar 

  101. Marini F, Caputo B, Pombi M, Tarsitani G, della Torre A. Study of Aedes albopictus dispersal in Rome, Italy, using sticky traps in mark-release-recapture experiments. Med Vet Entomol. 2010;24:361–8.

    Article  CAS  PubMed  Google Scholar 

  102. Eritja R, Escosa R, Lucientes J, Marquès E, Roiz D, Ruiz S. Worldwide invasion of vector mosquitoes: present European distribution and challenges for Spain. Biol Invasions. 2005;7:87–97.

    Article  Google Scholar 

  103. Lucientes-Curdi J, Molina-Moreno R, Amela-Heras C, Simon-Soria F, Santos-Sanz S, Sánchez-Gómez A, et al. Dispersion of Aedes albopictus in the Spanish Mediterranean Area. Eur J Public Health. 2014;24:637–40.

    Article  PubMed  Google Scholar 

  104. Montarsi F, Drago A, Martini S, Calzolari M, De Filippo F, Bianchi A, et al. Current distribution of the invasive mosquito species, Aedes koreicus [hulecoeteomyia koreica] in northern italy. Parasit Vectors 2015, in press.

  105. Cini A, Ioriatti C, Anfora G. A review of the invasion of Drosophila suzukii in Europe and a draft research agenda for integrated pest management. Bull Insectol. 2012;65:149–60.

    Google Scholar 

  106. Saponari M, Boscia D, Nigro F, Martelli GP. Identification of Dna Sequences Related to Xylella fastidiosa in Oleander, Almond and Olive Trees Exhibiting Leaf Scorch Symptoms in Apulia (southern Italy). J Plant Pathol. 2013;95:659–68.

    Google Scholar 

  107. Baldacchino F, Caputo B, Chandre F, Drago A, della Torre A, Montarsi F, et al. Control methods against invasive Aedes mosquitoes in Europe: a review. Pest Manag Sci. 2015;71(11):1471–85.

    Article  CAS  PubMed  Google Scholar 

  108. Hulme PE. Beyond control: Wider implications for the management of biological invasions. J Appl Ecol. 2006;43:835–47.

    Article  Google Scholar 

  109. Goedde N. The poisoning of lake davis: weighting the risks. Environs. 1998;21:3–24.

    Google Scholar 

  110. Simberloff D. How much information on population biology is needed to manage introduced species? Conserv Biol. 2003;17:83–92.

    Article  Google Scholar 

  111. Adams VM, Petty AM, Douglas MM, Buckley YM, Ferdinands KB, Okazaki T, et al. Distribution, demography and dispersal model of spatial spread of invasive plant populations with limited data. Methods Ecol Evol. 2015;6(7);782–94.

    Article  Google Scholar 

Download references


We are grateful to Francesca Bussola, Matteo Mazzucato, Silvia Ciocchetta, Entostudio s.r.l., ULSS 1, ULSS 2 and Veneto Region for data collection. The PhD Scholarship of Matteo Marcantonio is supported by FIRS > T (FEM International Research School Trentino). This work has been funded by the Autonomous Province of Trento (Italy), Research funds for Grandi Progetti, Project LExEM (Laboratory of excellence for epidemiology and modelling, and Veneto Region. We are grateful to the NASA Land Processes Distributed Active Archive Center (LP DAAC) for making the MODIS data available.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Matteo Marcantonio.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

Conceived and designed the study: MM, M Metz, AR. Collected data: FB, DA, FM, GC, SC. Analyzed the data: MM, M Metz. Contributed analysis tools: MM, M Metz, MN. Wrote the paper: MM, M Metz. All authors read and approved the final version of the manuscript.

Additional file

Additional file 1:

Title of data: R function to perform logistic regression with Bayesian inference. Description of data: R function to perform logistic regression with Bayesian inference. (PDF 37 kb)


Appendix A: Correlation matrix and MaxEnt response curves

Fig. 5
figure 5

Correlation Matrix: The correlation matrix for the predictors set, input of MaxEnt model

Table 7 Ranking of the 5 most important variables for MaxEnt model. We assigned a score ranging from 5 to 1 to the first 5 predictors for each of the three measurements of variable importance provided by MaxEnt. Afterwards, we summed the rank to provide an overall metric for variable importance
Fig. 6
figure 6

MaxEnt Response Curves: MaxEnt response curve for the two predictors with the highest overall ranking (as calculated in Table 7). a) MaxEnt response curve for bio04. b) MaxEnt response curve for bio05

Appendix B: Posterior Probability Distributions (PPDs) for the best model

Fig. 7
figure 7

Posterior Probability Distribution for the best model parameters: The red dashed lines represent the distribution of the priors while the black horizontal line is the 95 % High Density Interval of the distribution. The distribution average, lower and upper bound of 95 % HDI were also reported in the figure

Appendix C: Sigmoid curves

Fig. 8
figure 8

Sigmoid curves for PHY model: We present the three sigmoid functions used to transform the environmental parameters to separate suitability indices

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Marcantonio, M., Metz, M., Baldacchino, F. et al. First assessment of potential distribution and dispersal capacity of the emerging invasive mosquito Aedes koreicus in Northeast Italy. Parasites Vectors 9, 63 (2016).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: