{'fraction_non_alphanumeric': 0.03746103005362265, 'fraction_numerical': 0.032373113854595334, 'mean_word_length': 5.047662141779789, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 3, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 2, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Geographic shifts in the bioclimatic suitability for Aedes aegypti under climate change scenarios in Colombia

The Dengue, Chikungunya and Zika viruses are arboviruses predominantly transmitted to humans through the bite of the female mosquito Aedes aegypti. Currently, the vector represents a potential epidemiological risk in several Latin American and Pacific countries. However, little is known about the geographical distribution and bioclimatic suitability of this mosquito in the projected climate change scenarios in Colombia. Using a species distribution model of maximum entropy (MaxEnt) based on presence-only records obtained from Global Biodiversity Information Facility (GBIF), land elevation obtained from Shuttle Radar Topography Mission (SRTM) and bioclimatic variables (WorldClim), we produced environmental suitability maps of this mosquito vector for present and future geographic distribution. The future distribution were constructed based on the Community Climate System Model (CCSM4) for the years 2050 and 2070, projected according to the Representative Concentration Pathways (RCP) 2.6, 4.5 and 8.5 described by the Intergovernmental Panel on Climate Change (IPCC). For the current conditions, Colombia has ~140,612.8 square km of areas with the possible presence of the vector; however, for the future, this will be reduced by more than 30%. For the future conditions, the suitable areas for A. aegypti decreased compared to the present, mainly for the year 2070 under RCP scenarios 4.5 and 8.5, however, the probability of mosquito occurrence increases in some departments of Colombia. Areas susceptible to the presence of A. aegypti are affected by climate change. The Caribbean and Andean regions have a high probability of mosquito distribution; therefore, control and epidemiological surveillance are required in these areas. The results can serve as an input to define preventive and control measures, especially in areas with a higher risk of contracting the virus.


Introduction
The climate change (CC) can affect the areas of geographic distribution, seasonal activities, migratory patterns, abundances and interactions of species due to sudden changes in precipitation and temperature (Haines et al., 2006;IPCC, 2014). Recently, changes related to extreme climatic and meteorological phenomena have been evidenced, e.g., extremely cold temperatures, increase in extremely warm temperatures, the rise in maximum sea levels and a significantly increased cases of intense rainfall in various regions (IPCC, 2014), all of which are determinants of the geographical distribution of vectors like A. aegypti.
The current distribution of A. aegypti is broadly defined by regional temperature and precipitation trends, but the global CC has the potential to significantly alter the future geographic range of mosquito vectors (Danis-Lozano et al., 2015;Andriamifidy et al., 2019). The A. aegypti mosquito is highly dependent on climatic conditions; for example, if the temperature of the water rises, the larvae take a much shorter time to mature and consequently, there is a greater capacity to produce more offspring during the transmission period (Rueda et al., 1990;Andriamifidy et al., 2019). Besides, in warmer climates, the adult female mosquito digests blood faster and feeds more frequently, therefore, there is an increased transmission of diseases (Gillies, 1953;Andriamifidy et al., 2019). On the other hand, changes in humidity and precipitation can affect the behavio of adult mosquitoes and the availability of aquatic habitats necessary for larval development . Also, if precipitation increases, the number and quality of sites for mosquito breeding also increase (Githeko et al., 2000;Andriamifidy et al., 2019).
The A. aegypti mosquito is a species distributed most predominantly in the tropical and subtropical areas of the Americas, between lattudes 35 north and 35 south (Nelson, 1986), although it has been reported at latitude 45 too, possibly due to invasions during the warm season (Nelson, 1986). In this sense, A. aegypti represents a potential epidemiological risk in countries of Latin America and the Pacific, since it is one of the most important vectors in the transmission of established and emerging viral diseases such as Dengue virus (DENV, genus Flavivirus), Chikungunya (CHIKV, genus Alphavirus), Zika (ZIKV, genus Flavivirus) and yellow fever virus (genus Flavivirus) (Kantor, 2016;Rückert et al., 2017). For example, in Colombia, from 2007 to the epidemiological week 20 of 2018, there were 745,257 cases of Dengue and severe Dengue, 69, 572 of Chikungunya and 13,590 of Zika (Instituto Nacional de la salud, 2018). According to medical studies, A. aegypti can transmit a combination of the viruses simultaneously (Rückert et al., 2017). There are no antiviral agents to treat the infection by the virus, and there are no authorized vaccines, despite the existence of some promising possibilities for eradication of the diseases in their entirety (Githeko et al., 2000;Weaver et al., 2012;Messina et al., 2015).
The transmission of diseases through vectors such as A. aegypti is a permanent challenge to public health. Besides, the situation is more difficult in countries where mosquito proliferation is abundant, and resources for medical care and comprehensive vector control services are limited (Kouri et al., 2007), as in the case of Colombia. It is therefore vitally important to advance in understanding the geographical distribution of the species at the national level since public health resources can be better allocated in anticipation of the emergence of diseases in the most vulnerable and vector-suitable communities.
Bearing in mind that the presence of A. aegypti is closely related to environmental conditions (Grillet et al., 2010;Rogers et al., 2014), including climate, it is feasible to predict the geographical distribution of A. aegypti under climate change scenarios. Some models or tools allow estimating the possible distribution of species according to environmental requirements (Elith and Leathwick, 2009). Species Distribution Models (SDMs) are tools that correlate records of known occurrences and environmental conditions in the localities of occurrence to obtain ecological and evolutionary knowledge and to predict the distribution of species in the landscape, which sometimes require extrapolation in space and time (Elith and Leathwick, 2009).
Several modeling techniques are developed to predict the influence of climate change on the distribution of species. Most of the models correlate the current occurrence of the species with the climatic variables or through an understanding of the physiological responses of the species to climate change (Pearson and Dawson, 2003). Once the climatic variables that influence the distribution of a species are determined, it is possible to estimate the potential redistribution of the species for future climatic scenarios. The SDMs have been used to predict the current potential geographic distribution in insect populations, including mosquitoes, as well as changes in their distribution influenced by environmental and climatic changes (Peterson et al., 2005;Escobar et al., 2016). Some studies worldwide have used niche modeling techniques to determine the geographic distribution of A. aegypti in the face of climate change. Khan et al. (2016) could corroborate that niche modeling techniques are adequate to predict the distribution of risk by the Zika virus; this was supported by Kraemer et al. (2015), who found that the tropical areas of the world are suitable for the global distribution of vectors A. aegypti and A. albopictus. Recently, Alaniz et al. (2017) modeled the distribution of A. aegypti worldwide and found that the risk of contracting Zika for the Americas is concentrated in the Atlantic coast of South America and on the coastal zones of the Caribbean Sea in Central America and North America.
Despite the existence of studies on the distribution of A. aegypti at a global level, it is necessary to develop research at a national level, regional level, which will allow the identification and prioritization of the administrative divisions of Colombia (departments) and/or areas with the highest probability of the presence of the vector. Besides, in most studies the spatial resolution is low (for example, 5 or 10 km), which makes it impossible to identify with reasonable precision the local distribution of the vector and, therefore, the risk of contracting viruses and the application of mitigation strategies to them. In this sense, the main objective of our research is to predict the potential geographic distribution of the species A. aegypti under climate change scenarios in Colombia, taking into account a relatively precise spatial resolution (1 km 2 ) with respect to other studies. The present modeling and future geographic distribution changes of A. aegypti will provide a useful tool for disease management and health service delivery planning in the country.

Study area
The study area comprised the departments of Colombia except the San Andr es and Providencia.
We excluded The San Andr es and Providencia islands because of their size and the spatial resolution (1km) of the environmental variables used for the present and future scenarios. The exclusion of these islands facilitated the analysis and design of the geographical distribution maps of the mosquito and promoted not to use different spatial scales to visualize the results.
Colombia is located in the northwestern tip of South America, with the Antilles sea in the north, the Pacific Ocean in the west, by Venezuela and Brazil in the east, by Panama in the northwest and by Peru and Ecuador in the south (Ministerio de Asuntos Exteriores y de Cooperaci on, 2017).
Colombia is divided into five geographical regions with different climatic and terrain conditions ( Figure 1B): 1. The Caribbean region located in the north of the country receives an average annual rainfall between 500 and 2000 mm, with the maximum annual average temperature ranging between 28 and 32 C (high and medium Guajira) and between 24 and 28 C for the Sierra Nevada de Santa Martha and its surroundings (IDEAM, 2005). 2. The Andean Region located in the center of the country (northern ranges of the Andean mountain) records an average annual rainfall up to 2000 mm. Along the Eastern range and in the valleys of Alto Magdalena and Alto Cauca, the annual temperature averages between 24 and 28 C (valleys of the main rivers such as Magdalena, Cauca, Patía and Sogamoso). The highlands of Cundinamarca, Boyac a and Nariño, the mountainous area of the center of Antioquia, Cauca and record between 12 and 16 C, 8 C in the high mountains (IDEAM, 2005). 3. The Orinoquía Region located to the east of the country receives on average an annual rainfall between 2000 to 3000 mm (central and eastern zone) and the average annual temperature ranges between 24 and 28 C for the flat area and between 8 and 20 C for the Piedmont area (IDEAM, 2005). 4. The region of the Amazon located in the south of the country receives on average an annual rainfall between 3000 and 4500 mm, with the average annual temperature ranging between 24 and 28 C and in Piedmont between 12 and 20 C (IDEAM, 2005). 5. The Pacific Region located in the western part of the country receives an average annual rainfall between 3000 and 12,000 mm, and the average annual temperature ranges between 24 and 28 C (IDEAM, 2005).

Species occurrences and environmental data
We downloaded 128 occurrence data (with replicas) for A. aegypti at country level from the Global Compendium of A. aegypti occurrence available at the Global Biodiversity Information Facility (GBIF, www .gbif.org, consulted on November 2017). The localities that are not spatially independent would bias the modeling with MaxEnt (Veloz, 2009;Anderson and Gonzalez, 2011;Boria et al., 2014;Phillips et al., 2016), Therefore, we filtered the localities and thus improved the performance of the ecological niche models (Boria et al., 2014). To reduce the spatial bias of the raw localities, the Spatially Rarefied Occurrence tool of the SDMtoolbox 2.0 (Brown et al., 2017), package was used, reducing occurrence locations to a single point within 10 km as the Euclidean distance (Pearson et al., 2007;Anderson and Raza, 2010), which resulted from the analysis of heterogeneity among the environmental variables filtered in this study ( Figure 1 Supplementary material (SM)). To calculate climate heterogeneity the SDMtoolbox package was used following the steps below: 1) main analysis of components of bioclimatic variables, 2) calculation of heterogeneity and 3) selection of heterogeneity class and Euclidean distance (for example, if heterogeneity is classified as high, medium or low the localities of occurrence could be spatially filtered at 5 km 2 , 10 km 2 and 30 km 2 respectively) (Brown et al., 2017).
We used bioclimatic variables of Worldclim for present and future with a spatial resolution of approximately (~) 30 s of arc (~0.86 km 2 ) ( Table 1). The variables for the present corresponding to the average of the years 1950-2000 (Hijmans et al., 2005), (http://worldclim.org/ver sion2, consulted on 2017), for the future medium-term 2050 (average years 2041-2060) and long term 2070 (average years 2061-2080) (Hijmans et al., 2005), (http://www.worldclim.org/cmip5_30s, consulted on 2017). These variables were derived from monthly temperature and precipitation values to generate biologically significant variables (Fick and Hijmans, 2017). For future bioclimatic variables, we used the model of the Community Climate System (CCSM4) for the years 2050 and 2070 projected according to the RCP Representative Concentration Pathways 2.6, 4.5 and 8.5, which represent a radiative forcing of~3 W/m 2 , 4.5 W/m 2 and 8.5 W/m 2 respectively. The Intergovernmental Panel on Climate Change (IPCC) described these scenarios, which is part of the simulations of the fifth phase of the Coupled Model Intercomparison Project Phase 5 (CMIP5). Global Circulation Model CCSM4 was selected because it has been used in several studies to predict the climate in Colombia, yielding excellent results (Rodríguez, 2012;Palomino et al., 2015).
Generally, bioclimatic variables are highly correlated (Feilhauer et al., 2012). As a consequence, it can be a problem for parameter estimation, because it inflates the variance of the regression parameters, leading to the incorrect identification of relevant predictor variables in a statistical model (Dormann et al., 2013). To reduce the highly correlated bioclimatic variables, we used the tool Remove Highly Correlated Variables from the SDMtoolbox 2.0 package (Brown, 2014), with which we performed the multicollinearity test by using the Pearson correlation coefficient (r) to determine a set of independent variables (Sillero, 2011). Variables with a cross-correlation coefficient value equal to or greater than 0.8 were excluded (Yang et al., 2013;Beckham and Atkinson, 2017). It should be noted that to retain variables with a correlation less than or equal to 80%, the tool "Remove highly correlated variables" from the SDMtoolbox package was used (Brown et al., 2017). This tool evaluates correlations between all data in the input environment and then removes the layers that are correlated at the level specified by the user (Brown et al., 2017). This took into account layers that are frequently used in ecology and evolution studies and variables that best represent the original input climate data (as they directly reflect actual measurements) and are not derived from several layers or a subset of the data (Brown et al., 2017).
To explain the altitudinal variation of the A. aegypti distribution, we used altimetry data as another physical variable (~0.86 km 2 resolution) extracted from the Shuttle Radar Topography Mission (SRTM) available at www.srtm.usgs.gov, consulted on 2017). It should be noted that the same altimetry of the land was used for both present and future modeling. Besides, elevation does not directly affect the distribution of the species but instead changes in temperature and air pressure that vary with altitude (Pearson, 2010). The environmental layers were processed in the ArcGIS10.6 software to obtain the information of the area of interest (Colombia).

Model selection, configurations and evaluation parameters
For A. aegypti distribution, we used the Maximum Entropy species distribution model (MaxEnt, version 3.4.1, available at http://biodiversi tyinformatics.amnh.org/open_source/maxent/, consulted on 2018). MaxEnt estimates a probability distribution of species of interest by finding the probability distribution of maximum entropy . The model takes into account georeferenced presence data of a species at a geographical level and the different characteristics (environmental variables) that influence its distribution .
MaxEnt was chosen because it is one of the most frequently and accurate methods of predicting the distribution of species (Elith et al., 2006;Phillips et al., 2006;Wisz et al., 2008;Ortega-Huerta and Peterson, 2008;Tognelli et al., 2009;Gomes et al., 2018) and is one of the most-used to predict the geographic distribution of insects, such as A. aegypti (Cardoso-Leite et al., 2014;Kraemer et al., 2015;Fatima et al., 2016;Khan et al., 2016;Alaniz et al., 2017). Besides, it is one of the models that presents the best prediction capabilities when different sample sizes are used, and even with sizes less than 25 occurrence data (Pearson et al., 2007;Phillips and Dudík, 2008;Wisz et al., 2008).
MaxEnt requires the use of background data or pseudo-presence data, that is, regions where the species does not have points of presence (Barbet-Massin et al., 2012). We used background data for MaxEnt calibration. It is essential that the background data used for training MaxEnt model does not relate to areas where the species is absent due to dispersion limitations or biotic interactions (Anderson and Raza, 2010;Barve et al., 2011). To avoid bias by the randomly constructed background , a sample of background points can be selected within a maximum radial distance of known occurrences (Thuiller et al., 2009). A maximum radial search distance of 200 km was used which covered the area of a Convex Minimum Polygon (MCP) of 800,000 km 2 which was constructed taking into account the occurrence data of the species and similar procedures proposed by other authors (Figure 3 SM) (VanDerWal et al., 2009;Brown and Yoder, 2015). MCP construction and background data used SDMtoolbox package 2.0 (Brown et al., 2017). The implementation of this method reduces the distribution of the species to the appropriate habitat within an area of known occurrence (Brown et al., 2017). The modalities for present and future climates are adjusted to the predetermined settings in MaxEnt, with the exception of the following parameters: random test percentage and training equal to 30% and 70% respectively; the output of the models was logistic, because it improves the calibration of the model and can be interpreted as a response to the probability of finding suitable areas for the species, based on environmental variables; the entity classes were: hinge, linear and quadratic, taking into account the number of filtered occurrences ; the maximum number of background points equal 10 000 (Kramer-Schadt et al., 2013); the maximum number of interactions equal 5000 (Beckham and Atkinson, 2017); Regularization multiplier of 2, because it reduces the level of low levels (Radosavljevic and Anderson, 2014). Since Colombia presents significant spatial variations (for example, elevation of terrain), it is probable that spatial groups of localities are too adjusted to environmental biases. Therefore, the performance values of the model are inflated (Veloz, 2009;Hijmans and Hall, 2012;Boria et al., 2014). This usually the case when using randomly partitioned occurrence data sets (Radosavljevic and Anderson, 2014) and a geographic division is indispensable (Osborne and Suarez-Seoane, 2002). Besides, such evaluations become more critical when modeling the distribution of species through space or time under CC (Anderson, 2013). To reduce the environmental biases related to spatial groups of localities, 100 replicas (Beckham and Atkinson, 2017) of Bootstrap type (replicated sets of samples chosen by sampling with replacement) taking into account the percentage of random and training tests of 30% and 70% respectively. This configuration causes the models to be calibrated and evaluated with different permutations of the points of occurrence and background data in different regions of the study area.
Jackknife test was performed to identify the relative importance of the environmental variables. During the Jackknife test, each variable was excluded in turn, and a model was created with the remaining variables. Then a model was created using each variable in isolation. Also, a model was created using all variables as before (Phillips, 2010). The environmental variable with the highest gain, when used in isolation, made the most significant influence in the mosquito modeled (Phillips, 2010). The environmental variable that decreases the gain the most when it is omitted appears to have the most information that is not present in other variables, therefore, it is likely to be highly influential (Phillips, 2010). In addition, the model showed response curves that relate the probability of finding suitable areas for mosquito distribution, taking into account the environmental variables used in the modeling.
The precision of the species distribution model was evaluated using the Area under Receiver Operating Characteristics (ROC) curve, better known as AUC (Fielding and Bell, 1997). This technique is one of the most used to evaluate the prediction of species distribution models (Elith et al., 2006). It is interpreted as the probability that a randomly chosen presence location has a higher classification than a randomly chosen background point, i.e. the AUC is used to determine how the model distinguishes between presences and absences (Merow et al., 2013). The ROC curve shows the false positive rate (the model predicts that a species will be present where it has not been observed) on the "X" axis and the true positive rate (the predictions may coincide with the observed presences of the species) on the "Y" axis represented in the range of probability threshold values (Phillips, 2010). The closer the ROC curve is to the axis, the higher the area below the curve (AUC) and therefore, the more accurate the model (Phillips, 2010). The AUC was evaluated from 0.0 to 1, and the precision of the model is better when approaching 1, that is, there is more significant discrimination of suitable areas than those not suitable for the species . Models with an AUC value greater than 0.75 are classified as potentially useful for predicting the distribution of a species (Elith et al., 2006). Another way to validate the results was to contrast the maps of the potential distribution of the species with the real cases of Dengue, Chikungunya, and Zika historically reported in the country.
To determine the effect of climate change on the geographic distribution of A. aegypti, changes in distribution between current and future models were quantified. For that, the Distribution Changes Between Binary SDMs tool of the SDMtoolbox 2.0 package was used (Brown, 2014), which required binary maps (result of the 100 replications) of potential distribution with a logistic presence training threshold of 10 percentiles (Cardoso-Leite et al., 2014;Radosavljevic and Anderson, 2014). The distribution changes between the binary distribution models were determined to find the areas of potential distribution, range reduction, range expansion and zones without changes (Brown et al., 2017). We determined the distribution shift by estimating centroid changes in the distribution using the Centroid Changes tool of the SDMtoolbox 2.0 package (Brown et al., 2017). We determined the centroid changes to summarize the distribution changes of the nucleus in the ranges of A. aegypti, that is, by reducing the distribution of the mosquito to a single central point by creating a vector line that represents the magnitude and direction of the change in time (Brown et al., 2017).

Results and discussion
The bioclimatic variables used in the models are the following: average annual temperature (Bio1), Mean diurnal range (Mean of monthly (max temp -min temp)) (Bio2), isothermality ((Bio2/Bio7) Â 100) (Bio3), Annual precipitation (Bio12) and precipitation of driest month (Bio14). These variables were selected because they had a correlation less than or equal to 0.8 (Table 1 SM) (Yang et al., 2013;Beckham and Atkinson, 2017). The variable temperature annual range (Bio 7) showed a correlation less than 80%; however, it was not included in the final modeling because it did not contribute significantly in previous modeling tests, this being a filter used in other studies (Zeng et al., 2016). 54 spatially independent occurrence data were used to power the MaxEnt model (Figure 2: SM). Out of 128 occurrence data downloaded from GBIF were reduced to 115, due to replicas (12) and points outside the study area (1). After applying Spatially Rarefy Occurrence to the points of presence (115), we obtained 54 spatially independent data separated with a minimum radial distance of 10 km 2 (Brown et al., 2017).
The results revealed a good performance by MaxEnt in the model. The model had an excellent relationship with a random model ( Figure 4A, SM) obtaining an AUC (average of 100 replications) of training~0.812 ( Figure 4B, SM) (Hosmer and Lemeshow, 2000). The AUC obtained allows classifying the model as potentially useful to predict the distribution of the species, because it exceeded a value of 0.75 (Elith et al., 2006).
Among the determining variables in the distribution of A. aegypti, elevation stands out, because it presented greater gain when used in isolation and, therefore, seems to have the most useful information in itself. The bioclimatic variable that reduces the gain when it is omitted is Bio3; therefore, it seems to have the most significant amount of information that is not present in the other variables, as evidenced in a distribution study of mosquito-carried out in Brazil (Cardoso-Leite et al., 2014).
The response curves (Figure 2) reflect the environmental suitability of A. aegypti, considering a MaxEnt model created using only the corresponding variable. The ranges of environmental suitability with respect to the bioclimatic variables and the corresponding probability percentages are shown below: (Bio1) Annual Mean Temperature of~28-32.5 C with~64%, the result is comparable with the optimum range for the species ranging from~25 to 32 C (Khormi and Kumar, 2014). It should be noted that if the temperature exceeds 30 C, the adult reproductive capacity is affected (Carrington et al., 2013) and the extrinsic incubation period of the viruses decreases (Watts et al., 1987); however, (Attaway et al., 2016), pointed out that in laboratory tests, the mosquito has shown competition between vectors at 35 C (Lozano-Fuentes et al., 2012). (Bio2) Mean diurnal range (mean of monthly (max temp -min temp) of 14-15 C with~85%, this being consistent with other studies that suggest that the lower temperature range optimal for A. aegypti is 14-18 C) (Khormi and Kumar, 2014). (Bio3) Isothermality ((Bio2/-Bio7) Â 100) of~96-97.5% with~70%; (Bio12) Annual precipitation of 0-500mm with 68%, highlighting the propensity of the mosquito to develop with a precipitation greater than or equal to 55 mm. Our mosquito develops well (Wiwanitkit, 2006) (Bio14) with Precipitation of the driest Month~smaller than 10mm with~45% and~80% when approaching 0mm, this is consistent with what (Wiwanitkit, 2006) stated, which highlighted that the ideal precipitation values are greater than 2mm.
On the other hand, the Figure 2 response curve related to the elevation of terrain shows that A. aegypti can distribute more predominantly at elevations lower than 1000 meters above sea level (~50% probability) and is more suitable close to 0 MASL (~85% probability), this being similar to what other studies showed (Lozano-Fuentes et al., 2012;Santos and Meneses, 2017). Besides, as the elevation increases, the probability of finding suitable sites decreases (Fatima et al., 2016), highlighting that A. aegypti mostly limits itself to 1500 meters above sea level (Su arez and Nelson, 1981); however, in Colombia, it has been reported at elevations close to 2200 meters above sea level (Su arez and Nelson, 1981); according to our model, because at that elevation, there is~30% probability of finding suitable areas for the mosquito. This can be attributed to climate change, specifically to rising temperatures in high altitude areas with current temperature conditions below the innate mosquito thresholds (Lozano-Fuentes et al., 2012).
The geographical distribution of the A. egypti species varied spatially and temporally in the country. In the present conditions (Figure 3), the environmental suitability of the mosquito distribution is more likely in the Caribbean and Andean regions, especially in level areas, as shown by other studies (Alaniz et al., 2017;Tjaden et al., 2017). The distribution probability decreases in the Amazon region, due to the decrease in population density (Alaniz et al., 2017). In the Caribbean Region, the departments of Atl antico, Magdalena, C ordoba, Bolívar, Sucre and Cesar stand out in the distribution probability of the mosquito; as well as me plains of Valle del Cauca, Nariño, Cauca, Santander, Norte de Santander, Tolima, Antioquia, Huila, Cundinamarca, Caldas, and Risaralda in the Andean Region. The departments with a high probability of having suitable zones for the mosquito in the present, are largely consistent with the number of cases of Dengue, Chicungunya and Zika reported by the National Health Institute in Colombia (Figure 4). For example, from 2007 up to the epidemiological week 20 of 2018 in Colombia, 828,419 cases were presented due to the viruses, of which those in the Valle del Cauca, Santander, Tolima, Antioquia, and Norte de Santander accounted for 116,709,98,586,67,978,64,810 and 54,504 cases respectively (Instituto Nacional de la salud, 2018). Also, a study that searched for dengue cases reported in Colombia between 2004 and 2013 showed that dengue outbreaks are reported mainly in departments of Norte de Santander, Santander, Huila, Tolima, Valle del Cauca and Antioquia (Castrill on et al., 2015). The probabilistic map of potential geographic distribution present for A. aegypti in Colombia is shown below.
Below is the total number of cases of Dengue, Chikungunya and Zika reported in Colombia, sorted by department.
For the years 2050 and 2070, the geographical distribution of A. aegypti varied, compared to the present one. When quantifying the distribution areas under the climatic scenarios (Table 3), significant decreases in the ideal areas for the future were noticed. For the present climatic conditions, Colombia has~140,612.8 km 2 areas where the mosquito is possibly present; however, by 2050, they are reduced tõ 92,233,~94,231.9 and~94,670.7 km 2 under climatic scenarios RCP 2.6, 4.5 and 8.5 respectively, while for 2070 under the same scenarios, Figure 3. Probabilistic map of environmental suitability present for A. aegypti in Colombia. The red, yellow and green zones indicate high medium and low probability of mosquito distribution, respectively. the areas of possible presence are~93,809.2,~84,116.2 and~87,985.8 km 2 . It should be inferred that the situation is not straightforward but rather complicated as some areas will see an upsurge in mosquito distribution, while others are likely to see a decline in the future compared to the present. Also, the changes in species distribution area are not always linearly related to the intensity of emission scenarios (RCP) (Wilby and Dessai, 2010;Goberville et al., 2015). Thus, when modeling future distribution in our case, the same was observed for the 2070s RCP 4.5, as identified in other studies (Samy et al., 2016). This may be due to uncertainties arising from GCMs and emission scenarios.
The results indicate that suitable distribution areas for mosquito would be reduced by more than 30% in the future (2050 and 2070). When comparing future scenarios, it is evident that by 2050, the suitable areas for mosquitos increase with warming (RCP 8.5), while by 2070, there were fluctuations. However, there were significant decreases mainly in RCP 4.5 and RCP 8.5 with respect to RCP 2.6. This change is possibly associated with extreme temperature increases in 2070 over 2050. When the warming becomes high, the areas where the mosquitos were present are no longer suitable in the future. Temperature over the optimum levels affects the life cycle of mosquitos and its reproductive capacity (Rueda et al., 1990;Khormi and Kumar, 2014;Andriamifidy et al., 2019). Although, by 2070, the ideal distribution area for the mosquito decreases, the distribution centroid shift tends to move towards the north of the country as warming increased (RCP 8.5) and towards the south when the warming is low (RCP 2.6) (Figures 5 and 6).
Despite the fact that for the future, the ideal areas for the mosquito presented a decreasing tendency with respect to the present, in departments such as Nariño and Cauca, in some areas of the departments such as Norte de Santander, Choc o, Meta, Caquet a, and Amazonas, there was an expansion of areas suitable for the mosquito. In addition, the probability of finding the species increased from~89% in the present distribution to~92% in some areas of the departments that present climatic conditions suitable for A. aegypti. The results are comparable with those of previous studies, where it was concluded that the optimum areas for the development of A. aegypti will decrease by the year 2070, mainly in some countries of South America, and among these, Colombia (Khormi and Kumar, 2014). On the other hand, to analyze the results presented in this study, it is necessary to take into account that when rainfall is high, possibly in RCP 2.6, the number and quantity of sites for mosquito breeding also increases (Githeko et al., 2000), but when the rainfall increases, the temperature (RCP 8.5) decreases. Further, when the temperature increases, there is a higher capacity to produce more offspring, because the larvae take less time to mature and besides, mosquitoes increase during the transmission period (Rueda et al., 1990) and the female A. aegypti feeds more frequently (Gillies, 1953).
Environmental suitability maps and potential distribution changes for A. aegypti for the years 2050 and 2070 are shown in Figures 5 and 6.
The distribution centroid changes of the environmental suitability of the mosquito for the future are shown in Figure 7. The results indicate that the direction and magnitude of the changes are different in the climatic scenarios studied. The distribution centroid of the current habitat is located at 75.230 W longitude and 6.630 N latitude in the department of Antioquia, Colombia. The distribution Overall, the results indicated the core distributional shift expressed a relatively high magnitude (~150 km) towards the south when temperatures are low (RCP 2.6-2050), and towards north when temperatures are high (RCP 8.5-2070) distribute around 100 km respect to the current distribution centroid location 74.823 W, 6.422 N) (Figure 7).
Literature review shows that this is one of the first studies that predict the distribution and potential changes in the geographic range of A. aegypti in Colombia; taking into account current climatic conditions This research has broad implications in the area of health, especially in the prediction of infectious diseases. For example, for the National Institute of Health, in charge of epidemiological surveillance and control in Colombia, these results are of vital importance for the taking of preventive and control measures, especially in the areas with a higher risk of the presence of vectors.
On the other hand, this study is a valuable input for the construction of predictive models with MaxEnt, since it describes in a general way the procedures to improve the performance of the model and to avoid some biases when configuring the parameters of the software, taking into account the environmental predictors and the presence data of the species.

Conclusions
It can be concluded that the modeling with MaxEnt allowed identifying the potential zones of Colombia where A. aegypti can be distributed. The areas susceptible to the presence of A. aegypti are affected by climate change. The Caribbean and Andean regions presented suitable areas with high probability of mosquito distribution in the present and future climatic conditions. In respect of the Caribbean Region, this included some areas of the departments of Atl antico, Magdalena, C ordoba, Bolívar, Sucre and Cesar stand out in the distribution probability of the mosquito; as well as some plains of Valle del Cauca, Nariño, Cauca, Santander, Norte de Santander, Tolima, Antioquia, Huila, Cundinamarca, Caldas, and Risaralda in the Andean Region.
In respect of the future (2050 and 2070), the geographical distribution of the mosquito varied from the present, ruling out significant decreases in the areas ideal for the insect. In the present conditions, Colombia has~140,612.8 square km of areas with possible presence of the vector; however, for the future, this was reduced by more than 30%. On the other hand, when comparing the geographical distribution of the mosquito between 2050 and 2070, it was found that in 2050, there is a tendency for the ideal areas to increase, as the temperature increases (RCP8.5), while for 2070 under the same scenario, the areas suitable for the mosquito decreased.
The results found can serve as an input to take preventive and control measures, especially in areas with a higher risk of contracting the virus. Besides, some departments without evidence of being affected by A. aegypti in the present conditions could become areas of vulnerability in the future, with an active transmission. For this reason, it is necessary to implement or reinforce prevention measures for these areas.
We suggest modeling the distribution of A. aegypti taking into account other environmental variables, for example, population density, economic activities, water network, soil cover and among other variables which are related to the distribution and biology of the mosquito. Also, it is required to incorporate new occurrence data sets, for example, in the Orinoco region of Colombia.
Also, it is necessary to know the dynamics of the viruses transmitted by the vector to widely understand the distribution patterns of the same in Colombia and thus, take the relevant control measures to reduce the risk of transmission.
Finally, based on these models, we propose to create a GIS virtual platform for early warning system against communicable infectious diseases like Dengue, Zika, and Chikungunya in Colombia. Also, generate new predictions considering a more detailed spatial resolution to study the phenomena in the prioritized areas of the country.

Author contribution statement
Cristiam Victoriano Portilla Cabrera: Conceived and designed the experiments; Performed the experiments; Analyzed and interpreted the data; Contributed reagents, materials, analysis tools or data; Wrote the paper.
John Josephraj Selvaraj: Conceived and designed the experiments; Analyzed and interpreted the data; Contributed reagents, materials, analysis tools or data; Wrote the paper.

Funding statement
This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.