Ecological niche modeling of mosquito vectors of West Nile virus in St. John’s County, Florida, USA
© The Author(s). 2016
Received: 15 February 2016
Accepted: 15 June 2016
Published: 29 June 2016
The lack of available vaccines and consistent sporadic transmission of WNV justify the need for mosquito vector control and prediction of their geographic distribution. However, the distribution of WNV transmission is dependent on the mosquito vector and the ecological requirements, which vary from one place to another.
Presence/density data of two WNV mosquito vectors, Culex nigripalpus and Cx. quinquefasciatus, was extracted within 5 km buffer zones around seropositive records of sentinel chickens in order to delineate their predicting variables and model the habitat suitability of probable infective mosquito using MaxEnt software. Different correlations between density data of the extracted mosquito vectors and 27 climate, land use-land cover, and land surface terrain variables were analyzed using linear regression analysis. Accordingly, the correlated predicting variables were used in building up habitat suitability model for the occurrence records of both mosquito vectors using MaxEnt.
The density of both WNV mosquito vectors showed variation in their ecological requirements. Eight predicting variables, out of 27, had significant influence on density of Cx. nigripalpus. Precipitation of driest months was shown to be the best predicting variable for the density of this vector (R 2 = 41.70). Whereas, two variables were proven to predict the distribution of Cx. quinquefasciatus density. Vegetation showed the maximum predicting gain to the density of this mosquito vector (R 2 = 15.74), where nestling birds, in particular exotics, are found. Moreover, Jackknife analysis in MaxEnt demonstrated that urbanization and vegetation data layers significantly contribute in predicting habitat suitability of Cx. nigripalpus and Cx. quinquefasciatus occurrence, respectively, which justifies the contribution of the former in urban and the latter in epizootic transmission cycles of WNV. In addition, habitat suitability risk maps were produced for both vectors in response to their predicting variables.
For the first time in the study area, a quantitative relationship between 27 predicting variables and two WNV mosquito vectors within their foraging habitats was highlighted at the local scale. Accordingly, the predicting variables were used to produce a practical distribution map of probable infective mosquito vectors. This substantially helps in determining where suitable habitats are found. This will potentially help in designing target surveillance and control programmes, saving money, time and man-power. However, the suitability risk maps should be updated when serological and entomological data updates are available.
KeywordsWest Nile virus Culex nigripalpus Culex quinquefasciatus Habitat suitability Florida
West Nile virus (WNV) disease was reported for the first time in the United States in New York City during 1999  with a severe mortality of some wild native birds, in particular American crows, and many species of exotic birds . Since that time, increased public health concerns have surfaced over mosquito-borne diseases due to their consistent and sporadic autochthonous transmission in the state of Florida. Eventually, WNV spread to St. John’s County, FL by 2001 . Details concerning the nature of transmission dynamic and predicting variables of WNV disease in the United States remain controversial, especially with the continuous change in global climate and land use-land cover (LULC) including urban expansion and human population growth. Moreover, mosquito vector biology and their ecological requirements play a major role in the variation of the transmission cycle of WNV.
Both Culex nigripalpus Theobald and Cx. quinquefasciatus Say were recorded as the main vectors of WNV in the state of Florida . The former mosquito species was incriminated in the urban transmission cycle of the virus in Peninsular Florida with a feeding preference on human blood [5–8], whereas the latter species was reported to be responsible for an epizootic cycle and sustaining the virus circulation within reservoir host bird(s) [5, 9]. However, Cx. quinquefasciatus alongside with Cx. pipiens complex, were incriminated in the urban transmission cycle of WNV in northeastern and central parts of the USA [10–12]. The biology and ecology of WNV mosquito vectors may show some elasticity from one place to another based on the available resources in the surrounding environments. Accordingly, the transmission cycle of WNV may show some variations in response to the biology and ecology of these mosquito vectors.
Surveillance and control programmes of mosquito vectors are the most effective tools for arbovirus disease prevention . However, these programmes have low priority and lack adequate funds . Prediction models of suitable habitats for arbovirus transmission in terms of mosquito vector distribution and seropositive sentinel chickens could maximize the potentiality of surveillance and control programmes to break the transmission cycle of WNV. Moreover, the lack of available potential vaccines for WNV and consistent development of insecticide resistance for mosquito vectors justify the necessity to model arbovirus disease so that host spot areas can be targeted during surveillance and control activities.
The previous models of WNV transmission dynamics are based solely on hydrological and meteorological data [5, 9, 15–17]. In addition, other models have been based on socio-environmental predictors in terms of vegetation or urban and sub-urban areas . These previous studies addressed three major mosquito vectors in WNV transmission, Cx. pipiens, Cx. quinquefasciatus, and Cx. nigripalpus. The feeding preference of these mosquitoes may show significant variation from one place to another, consequently the disease transmission cycle may vary. Although these models are useful, their findings did not adequately account for the interrelationships between climate and non-climate variables such as LULC and Digital Elevation Models (DEM) and the overall influence of these interrelationships on arbovirus transmission cycles . In addition, some of them either predicted the distribution risk of WNV on regional scale or used data points of WNV seropositive and mosquito vectors instead of region representing the flight range of the mosquito vectors.
Culex nigripalpus and Cx. quinquefasciatus may share some ecological habitats, but their biological and ecological requirements in terms of blood meal source, water habitat quality, vegetation type, and human activities may show significant variation. Distribution models that include climate, LULC, and DEM variables may give detailed information on variation of suitable habitats between both mosquito vectors and the virus transmission dynamics.
To address the limitations of previous studies, we developed a habitat suitability model for probable infective WNV mosquito vectors within their flight range using data records on mosquito vector occurrence and WNV seropositive sentinel chickens. In addition, we addressed the interrelationships between climate, LULC and DEM variables and their overall influence on the habitat suitability of probable infective mosquito vectors using the Maximum Entropy (MaxEnt) tool. The current study was conducted in St. John’s County, Fl. A total of 43 mosquito species were previously reported inhabiting a great diversity of salt and fresh water habitats throughout the County [20, 21]. Eleven of these mosquito species are known vectors of pathogens, notably arthropod-borne viruses (arboviruses) [3, 4, 22, 23]. The Anastasia Mosquito Control District (AMCD) oversees the entire County’s mosquito surveillance and control programmes to alleviate the risk of mosquito-borne disease transmission. Since the local economy is primarily driven by the tourism industry, the management of mosquito-borne diseases has a significant economic impact. The occurrence of arbovirus cases has a negative impact on tourism, making AMCD surveillance and control programmes an important element to the economic health of the county.
St. John’s County is located in northeastern Florida, USA, with a total area of 1588 km2 intermingled between the St. Johns River to the west and the Atlantic Ocean to the east.
Sixty chickens (21-week old) were distributed over 12 locations in the county from the beginning of April through December every year during 2007–2013 to monitor arbovirus activity. Approximately 2.0 ml of blood was taken from each chicken’s wing vein once a week. Blood samples were kept in labeled vacutainers (Fisher Scientific®) and transported back to the laboratory at the AMCD Base Station in St. Augustine Beach where they were centrifuged at 4375 RPM for 15 minutes. Samples were then placed in a labeled and sealed plastic bag, shipped to the Florida State Department of Health (FDOH) Virus Laboratory in Tampa, Florida and tested for WNV. Samples were sent on the first operation day of the week and results reported by the end of the same week. Once a chicken seroconverted, it was removed, destroyed, and replaced with a new one at its respective location.
Mosquito sampling data
Population density (total number of mosquito vector/season/5 km) of adult host seeking mosquito vectors were monitored by AMCD using CDC light traps (John W. Hock Company, Gainesville, FL) baited with dry ice (BioSensory, Inc.) at 68 permanent locations in the County. Traps were placed outdoors at the beginning of March through the last week of November during 2007–2013. Traps were suspended 1 m above ground surface by a shepherd’s hook and operated for 18–20 h using a 12 V battery. Mosquito collections were transported weekly from the field to the AMCD facility for further identification to species level using the taxonomic keys of Darsie & Ward .
Twenty seven variables used in predicting suitable habitats of WNV mosquito vectors in St. John’s County, Fl
Elevation in meters
Annual Mean Temperature
Mean Temperature of Warmest Quarter
Mean Temperature of Coldest Quarter
Precipitation of Wettest Month
Precipitation of Driest Month
Precipitation Seasonality (Coefficient of Variation)
Precipitation of Wettest Quarter
Precipitation of Driest Quarter
Precipitation of Warmest Quarter
Precipitation of Coldest Quarter
Mean Diurnal Range (Mean of monthly (max temp - min temp))
Isothermality (BIO2/BIO7) (* 100)
Temperature Seasonality (standard deviation *100)
Max Temperature of Warmest Month
Min Temperature of Coldest Month
Temperature Annual Range (BIO5-BIO6)
Mean Temperature of Wettest Quarter
Mean Temperature of Driest Quarter
Leaf Area Index
Human population settlements
Slope measures the rate of changes of elevation at surface location, expressed as an angle from 0 (flat) to 90° (high elevation). Aspect ratio indicates the orientation of the slope, eventually this ratio reflects the places of water accumulation and larval mosquito establishment [29, 30]. Aspect ratio ranges from 0 to 360°, however, in this study, aspect ratio was cosine-transformed to represent ranges from +1 (north-facing slope) to −1 (south-facing slope) . Curvature and hill shade reflect depressions and sunlight intensity in/on land surface. The four indicators were generated from a 30 arc-seconds Digital Elevation Model (DEM).
Land use - land cover data
Because the distribution of mosquitoes is greatly affected by the human population as a source of blood meal , urban areas layer was included as a predictor for distribution of human population (Table 1). Since vegetation represents resting places and sugar meal source for adult mosquitoes, it has been extensively used as a predictor for their occurrence [33–35]. Vegetation land cover used in our model was expressed as Leaf Area Index (LAI). The LAI captures vegetation characteristics such as canopy cover and sugar resources. The LAI was derived from Moderate Resolution Image Spectroradiometer (MODIS) imagery from the Terra Satellite. The LAI is defined as the one-sided green leaf area per unit ground area in broadleaf canopies and as one half the total needle surface area per unit ground area in coniferous canopies. The LAI represented a good indicator for LULC in canopies and forest areas, such as the case in St. John’s County rather than Normalized difference vegetation index (NDVI)  that gives a broader depiction about the green areas. Therefore, we used LAI as a potential remote sensing indicator for vegetation cover to assess their predictive power on distribution of both WNV mosquito vectors.
The surface water has been extensively used in previous studies as an indicator for mosquito distribution range . Thereby, data layers on surface water, including water bodies, rivers and streams, were also included in the linear regression model. Types of surface water were imported, classified according to their types, clipped to the buffer zones around the sentinel chickens. Each water body type was extracted as a separate ASCII file to be used in further analysis.
Habitat suitability modeling of WNV vectors
The MaxEnt software v. 3.3 [39–41] was used to model the suitable habitats of both mosquito vectors. The software uses the occurrence (presence of mosquito vector at least one time/trap/season) records of mosquito vectors in association with the predicting variables to generate suitability risk maps.
The contribution of predicting variables in our model was evaluated using Jackknife analysis in MaxEnt. The predicted habitat probability was categorized into five classes using the natural area breaks in ArcGIS utilizing the training omission rate at 10 percentile training presence value produced from MaxEnt: very low (0–0.1), low (> 0.1–0.2), medium (> 0.2–0.4), high (> 0.4–0.6), and very high (> 0.6) using natural breaks in the symbology tools in ArcGIS.
The presence records of both mosquito vectors were randomly partitioned for model evaluation into two subsamples: 75 % of the records were used for training and building up the model, and 25 % of the records were used for testing the model’s accuracy. The duplicate records of WNV mosquito vectors within ~1 km of the same cell size were excluded .
Five replicate runs were assigned in running the model to generate the average, maximum, minimum and median of the distribution range of mosquito vectors. Two thresholds have been used to examine the performance accuracy of our model based on Phillips et al.  and Phillips & Dudik : (i) the extrinsic omission was evaluated at fixed threshold (10 percentile training presence) and (ii) the area under the curve (AUC) of the receiver operating characteristics (ROC) [40, 41].
To increase the potentiality of the habitat suitability model and maximize the sampling effort of extracted mosquito vectors, a separate ASCII file generated from the 5 km buffer zones was included and weighted to the corresponding mosquito vector density . Accordingly, the habitat suitability of each of the mosquito vectors will be predicted in regards to the flight range of probable infective mosquito and response to the predicting variables. Therefore, the distribution map produced for each mosquito vector reflects the presence of probable infective mosquito. The MaxEnt evaluates different correlations between the presence records of extracted mosquito vectors data and their predicting variables within the sampled areas utilizing logistic regression analysis. In our study, we weighted the presence records to density of mosquito vectors 5 km around seropositive sentinel chickens and attached as ASCII bias file. The sampling points were randomly grouped into 75 % and 25 % for data training and testing, respectively. During data training, a matrix of spatial correlations between sampling points and their associated predicting variables were created. Accordingly, the habitat suitability maps were created for sampled and unsampled areas based on the habitat similarity between sampled and unsampled regions. The suitable habitats for unsampled areas were predicted utilizing habitat similarity between sampled and unsampled ones.
Data layers and variables selection
Summary of stepwise linear regression analysis on density of both WNV mosquito vectors in response to bioclimatic, LULC and DEM variables
Ecological niche modeling of west Nile virus vectors
In the current study, we produced a practical habitat suitability map for probable infective WNV mosquito vectors in St. John’s County. We addressed limitations in previous comparable models in regards to (i) use ~1 km resolution remote sensing data layers in terms of climate, ecological, and DEM variables, (ii) ranking the average contribution of predicting variables correlated to the density and presence of WNV mosquito vectors, and (iii) production of practical distribution risk map for the probable infective mosquito vectors confined to their flight range around seropositive records at local scale to prioritize the surveillance and control activities.
Generally, the predicted habitat suitability maps showed a significant variation in distribution between both mosquito vectors. The very high suitable habitats for Cx. nigripalpus was sporadically distributed in the central and southwestern areas. Unlikely, the habitat suitability produced for Cx. quinquefasciatus was largely confined to the western areas of the county (Figs. 4 and 6). However, the predicted distribution of both mosquito vectors shared some suitable habitats in the western regions.
Since the distribution of mosquito vectors and the virus circulation are confined to areas rather than sampling points, we extracted mosquito sampling points within their 5 km flight range around WNV seropositive records. This reflects the association between a probable infective mosquito population, WNV reservoir bird host(s), and water breeding habitats [25, 26]. Both WNV mosquito vectors are well known as exophilic and exophagic and their adult breeding habitats range from ditches, woodland pools, and freshwater marshes of a semi-permanent or permanent nature. In addition, Cx. quinquefasciatus shows preference to polluted water ranging from agricultural drains and sewage canals. Cx. nigripalpus shows feeding preference for humans rather than Cx. quinquefasciatus, which prefers avian blood [6–9, 15–17]. However, the former vector may feed on avian blood during the dry season, when it comes in close vicinity with reservoir host bird(s) seeking water breeding habitats [9, 15–17]. The LRM demonstrated that precipitation during the wettest (R 2 = 38.47) and driest (R 2 = 40.70) months have a significant role in predicting Cx. nigripalpus density during the dry season (Table 2) [5, 22]. Accordingly, MaxEnt confirmed this finding and demonstrated the significant contribution of urbanization (25.5 %) (Table 1) in predicting Cx. nigripalpus, which reflects the importance of this mosquito vector in the urban transmission cycle rather than Cx. quinquefasciatus [5–7, 16].
Although no WNV human cases were recorded during the study period, the transmission potential still exists due to the availability of competent mosquito vectors, reservoir hosts, probably exotic birds, and suitable ecological habitats. In addition, the lack of data on vector survival rate, feeding preference, gonotrophic cycle period, and competence for the study area prevented us including this information in building up our model. Our spatial distribution model shed light on the information gap in vector capacity parameters, which are needed to deeply understand when, precisely, transmission onset occurs in St. John’s County.
The distribution of these mosquito vectors, like other exophilic mosquitoes, are influenced by surrounding climate and ecological variables. In addition, this distribution gives deep insights on the virus circulation between the bridging and main mosquito vector, reservoir bird and human hosts. The findings of some previous models showed great potential in understanding the biology and ecology of WNV mosquito vectors in state of Florida [9, 16, 22], however, they encountered some limitations in emphasizing the interaction between ecological, climate, DEM variables and their overall influence on mosquito vector distribution at a local scale. Other models concluded the association between WNV transmission was either with forested and urban land [11, 44] or socioeconomic status . However, their findings may not be applicable in other landscapes that encounter WNV transmission due to the heterogeneity in climate and landscape variables . Moreover, some of the previous models predicted the spatial distribution risk of virus circulation based on sampling points rather than the flight range area of these mosquito vectors.
Unlike previous models, we tried to model the spatial distribution of probable infective mosquitoes using the mosquito presence records within their flight range around WNV seropositive records of sentinel chickens. Since arbovirus transmission is vector-density dependent, we evaluated the response of both vectors’ density to their predicting variables. The LRM provides reliable information on the influence of the predicting variables on mosquito vectors density. In this regard, the predicting variables were evaluated and selected using AICc values. Statistical LRM potentially resolves the significant discrepancy in predicting variables for both mosquito vectors. This discrepancy reflects variation in ecological requirements between both mosquito species. The potentiality of the current prediction model was proven by the high AUC values produced by MaxEnt, indicating that occurrence records were likely to be assigned a higher probability of presence than background sites.
Unsurprisingly, temperature and precipitation related variables significantly increased the prediction gain of probable infective Cx. nigripalpus in regards to their density and presence records. Their contribution was demonstrated by LRM and MaxEnt tools in terms of R 2 values and percent contribution, respectively. The contribution of precipitation in driest (R 2 = 41.70, 1.9 %) and wettest (R 2 = 38.47, 13.6 %) months reflects their significance in determining the habitat suitability for this mosquito vector during dry and wet seasons. The shared high predictive power by temperature related variables (32.5 %) confirmed the association between this mosquito vector, reservoir hosts, and human population, especially during the dry season. The significance of urbanization in predicting Cx. nigripalpus demonstrates the suitability of this habitat for mosquito vector density (R 2 = 29.52) and presence probability (25.5 %).
Accordingly, precipitation in the driest months, during the dry season in May – June, may play a major role in shifting breeding habitats of this mosquito and bring them in contact with WNV infected host birds . Alternatively, precipitation of wettest months during July and August increases the water habitat suitability and density/occurrence probability of mosquito vectors. Eventually, this mosquito vector starts to get closer to urbanization areas again and bites human [9, 16, 17]. In addition, this may reflect the importance of this mosquito vector in the urban cycle and infection of the human population, which is confirmed by their feeding preference for humans in previous studies [5, 15, 16, 22].
The contribution of temperature in predicting density and distribution of this mosquito vector was highlighted in a mathematical simulation model by Lord & Day , in which the mortality of WNV mosquito vector populations at 22.5 °C was included as a determinant parameter for disease transmission. However, the number of parameters regarding mosquito biology and ecology prevented investigating different correlations in depth. In the current study, both LRM and Jackknife analysis showed these different interrelationships between Cx. nigripalpus and correlated variables within mosquito vector foraging habitats. The temperature of the coldest quarter (< 21 °C) contributed significantly in predicting the occurrence and spatial distribution of this mosquito vector.
Land topography and geomorphology were potentially used in predicting suitable water habitats of mosquitoes [30, 31, 47]. Four potential indicators were highlighted in previous studies to predict suitable water habitats of mosquitoes: aspect ratio, slope, land surface curvature and hill shade [29, 30, 48, 49]. These indicators give detailed information about locations where water flows and accumulates. Although we did not include data on larval mosquitoes, data on host seeking adult mosquitoes included in the current study reflects the vicinity of water habitats to collected samples. In regards to elevation, it defines soil-water gravitational potential energy  and surface water movement within drainage channels as well as throughout the landscape. Similarly, slope has the potential influence on surface water flow velocity, drainage and accumulation of water [29, 49].
In the present study, hill shade was used as a promising remote sensing data that reflects areas with surface water accumulation . Hill shade showed a significant prediction gain with precipitation and temperature related variables. Among five land surface geomorphic predictors representing elevation, slope, aspect ratio, curvature, and hill shade, only the latter was demonstrated as a potential predicting variable in determining the occurrence/density and distribution of Cx. nigripalpus using LRM and MaxEnt. Hill shade was used in previous studies as a predicting precursor for water suitable habitats of malaria vectors .
Density of Cx. quinquefasciatus was evidently correlated with vegetation and mean temperature diurnal range. Vegetation is important for mosquitoes as a source of sugar meal, but also reflects nestling habitats for reservoir bird hosts, probably human hosts in rural areas. Feeding preference of this mosquito vector is reportedly varied ; however, the blood meal identification in fed mosquitoes revealed their feeding preference to avian blood rather than human [6–8]. Accordingly, this mosquito vector may play a major role in sustaining the epizootic cycle of WNV.
Although some previous models showed potential in understanding the empirical correlations between WNV transmission and some ecological and climate variables in the state of Florida, these models encountered some limitations. These models predicted the spatial distribution risk of virus circulation based on sampling points rather than the flight range area of probable infective mosquito vectors at the local scale. Other models may not be applicable in other landscapes due to the heterogeneity in climate and environmental variables. The variation of landscape and climate variables greatly impact the biology and ecology of mosquito vectors in terms of feeding preference, habitat suitability, and the distribution of reservoir bird host(s), especially exotic birds.
In this model, we demonstrated how GIS, remote sensing data, and habitat suitability modeling tools can be used efficiently to analyze and identify suitable habitats for infective WNV mosquito vectors and assess the distribution risk of this disease at a local scale. Our model is very novel in addressing limitations of previous similar models at the local scale. The risk maps produced will effectively help in determining where suitable habitats are found. This will potentially help in targeted surveillance and control programmes, which will save money, time and man-power, especially if the suitability risk maps are updated with serological and entomological data when available.
The availability of consistent temporal data, in terms of mosquito and remote sensing, limit our current spatial model to delineate accurately, when WNV transmission can occur. Further temporal analysis is needed to couple the current spatial model with the updated data on mosquito distribution to give more detailed information on where and when the suitable habitats of WNV transmission may occur. In addition, information on vector capacity in regards to vector competence, survival rate, and gonotrophic cycle period, need to be considered in building up the spatial-temporal models. Eventually this will give deep insights on the competency of mosquito vector populations in transmitting WNV.
AICc, Corrected Akaike information criterion (AICc); AMCD, Anastasia mosquito control district; ASCII, American standard code for information interchange; AUC, area under curve; Bioclim, bioclimate variables; Cx, culex; DEM, digital elevation model; FDOH, Florida Department of Health; hrs, hours; km, kilometer; LAI, leaf area index; LRM, linear regression model; LULC, land use - land cover; MaxEnt, maximum entropy; MODIS, Moderate Resolution Image Spectroradiometer (MODIS); NDVI, Normalized difference vegetation index (NDVI); ROC, receiver operating characteristics; WNV, West Nile Virus
Authors would like to express their deep gratitude to Anastasia Mosquito Control District for the logistic support, especially Mr. Richard Weaver, Mrs. Marcia K. Gaines, Mrs. Jennifer Gibson for their kind help in data preparation. Also, we would like to thank employees at The Urban Entomology Laboratory, University of Florida for their consistent support.
Institute of Food and Agriculture Sciences/IFAS.
Availability of data and materials
Contact the second author: email@example.com
RDX carried out the experimental design and data sampling. MFS: carried out geo-database and model building, data analysis and writing the first draft of the manuscript. All authors contributed equally in validating the model and approved the final version of the manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
The current study has been approved by the ethical review committee of Institute of Food and Agriculture Science, University of Florida, Gainesville, approval no. 201303836_01.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), 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 (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Centers for Disease Control and Prevention. Guidelines for surveillance, prevention and control of West Nile virus infection—United States. MMWR. 2000;49:25–8.Google Scholar
- Steele KE, Linn MJ, Schoepp RJ, Komar N, Geisbert TW, Manduca RM, et al. Pathology of fatal West Nile virus infections in native and exotic birds during the 1999 outbreak in New York City, New York. Vet Pathol Online. 2000;37(3):208–24.View ArticleGoogle Scholar
- Blackmore CG, Stark LM, Jeter WC, Oliveri RL, Brooks RG, Conti LA, et al. Surveillance results from the first West Nile virus transmission season in Florida, 2001. Am J Trop Med Hyg. 2003;69(2):141–50.PubMedGoogle Scholar
- Sardelis MR, Turell MJ, Dohm DJ, O’Guinn ML. Vector competence of selected North American Culex and Coquillettidia mosquitoes for West Nile virus. Emerg Infect Dis. 2001;7(6):1018–22.View ArticlePubMedPubMed CentralGoogle Scholar
- Day JF. Predicting St. Louis encephalitis virus epidemics: lessons from recent, and not so recent, outbreaks. An Rev Entomol. 2001;46(1):111–38.View ArticleGoogle Scholar
- Kilpatrick AM, Daszak P, Jones MJ, Marra PP, Kramer LD. Host heterogeneity dominates West Nile virus transmission. Proc R Soc Biol. 2006;273:2327–33.View ArticleGoogle Scholar
- Molaei G, Andreadis TG, Armstrong PM, Bueno R, Dennett JA, Real SV, et al. Host feeding pattern of Culex quinquefasciatus (Diptera: Culicidae) and is role in transmission of West Nile Virus in Harris County. Texas Am J Trop Med Hyg. 2007;77(1):73–81.PubMedGoogle Scholar
- Turell MJ, Spring AR, Miller MK, Cannon CE. Effect of holding conditions on the detection of West Nile viral RNA by reverse transcriptase-polymerase chain reaction from mosquito (Diptera: Culicidae) pools. J Med Entomol. 2002;39(1):1–3.View ArticlePubMedGoogle Scholar
- Shaman J, Stieglitz M, Stark C, Le Blancq S, Cane M. Using a dynamic hydrology model to predict mosquito abundances in flood and swamp water. Emerg Infect Dis. 2002;8(1):6–13. Epub 2001/12/26.View ArticlePubMedGoogle Scholar
- Brown H, Diuk-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(2):197–206.View ArticlePubMedGoogle Scholar
- Brown HE, Childs JE, Diuk-Wasser MA, Fish D. Ecological factors associated with West Nile virus transmission, Northeastern United States. Emerg Infect Dis. 2008;14(10):1539–45.View ArticlePubMedPubMed CentralGoogle Scholar
- DeGroote JP, Sugumaran R, Brend SM, Tucker BJ, Bartholomay LC. Landscape, demographic, entomological, and climatic associations with human disease incidence of West Nile virus in the state of Iowa, USA. Int J Health Geogr. 2008;7:19.View ArticlePubMedPubMed CentralGoogle Scholar
- Eldridge BF. Strategies for surveillance, prevention, and control of arbovirus diseases in western North America. Am J Trop Med Hyg. 1987;37:77S–86S.PubMedGoogle Scholar
- Gubler DJ. The global emergence/resurgence of arboviral diseases as public health problems. Arch Med Res. 2002;33:330–42.View ArticlePubMedGoogle Scholar
- Day JF, Shaman J. Using hydrologic conditions to track the risk of focal and epidemic arboviral transmission in peninsular Florida. J Med Entomol. 2008;45:458–69.View ArticlePubMedGoogle Scholar
- Shaman J, Day JF, Stieglitz M. Drought-induced amplification of Saint Louis encephalitis virus, Florida. Emerg Infect Dis. 2002;8(6):575–80. Epub 2002/05/25.View ArticlePubMedPubMed CentralGoogle Scholar
- Shaman J, Day JF, Stieglitz M. St. Louis encephalitis virus in wild birds during the 1990 south Florida epidemic: the importance of drought, wetting conditions, and the emergence of Culex nigripalpus (Diptera: Culicidae) to arboviral amplification and transmission. J Med Entomol. 2003;40(4):547–54. Epub 2003/12/19.View ArticlePubMedGoogle Scholar
- Naish S, Mengersen K, Hu W, Tong S. Forecasting the future risk of Barmah forest virus disease under climate change scenarios in Queensland, Australia. PLoS ONE. 2013;8(5):e62843. doi:10.1371/journal.pone.0062843.View ArticlePubMedPubMed CentralGoogle Scholar
- Rogers D, Randolph S. The global spread of malaria in a future, warmer world. Science. 2000;289:1763–6.View ArticlePubMedGoogle Scholar
- AMCD. About AMCD operations. 2012.Google Scholar
- Naranjo D, Qualls W, Jurado H, Perez J, Xue R-D, Gomez E, et al. Vector control programs in Saint Johns County, Florida and Guayas, Ecuador: successes and barriers to integrated vector management. BMC Public Health. 2014;14(1):674.View ArticlePubMedPubMed CentralGoogle Scholar
- Lord CC, Day JF. Simulation studies of St. Louis Encephalitis and West Nile Viruses: The impact of bird mortality. Vector Borne Zoonotic Dis. 2001;1(4):317–29.Google Scholar
- Cupp EW, Klingler K, Hassan HK, Viguers LM, Unnasch TR. Transmission of Eastern Equine Encephalomyelitis virus in Central Alabama. Am J Trop Med Hyg. 2003;68(4):495–500.PubMedPubMed CentralGoogle Scholar
- ESRI. ArcGIS desktop: release 10. Redlands: Environmental Systems Research Institute; 2011.Google Scholar
- DeMeillon B. Studies on insects of medical importance in South Africa. South African Inst Med Res. 1934;VI:249–308.Google Scholar
- Nayar JK, Sauerman Jr DM. A comparative study of flight performance and fuel utilization as a function of age in females of Florida mosquitoes. J Insect Physiol. 1973;19(10):1977–88.View ArticlePubMedGoogle Scholar
- Darsie RF, Ward RA. Identification and geographical distribution of the mosquitos of North America. North of Mexico: University Press of Florida; 2005.Google Scholar
- Hijmans RJ, Cameron SE, Parra JL, Jones PG, Jarvis A. Very high resolution interpolated climate surfaces for global land areas. Int J Climatol. 2005;25:1965–78.View ArticleGoogle Scholar
- Gorsevski PV, Gessler P, Foltz RB. Spatial prediction of landslide hazard using discriminant analysis and GIS. GIS in the Rockies 2000 Conference and Workshop Applications for the 21st Century Denver, Colorado September 25 - 27, 2000. 2000.Google Scholar
- Mushinzimana E, Munga S, Minakawa N, Li L, Feng CC, Bian L, et al. Landscape determinants and remote sensing of anopheline mosquito larval habitats in the western Kenya highlands. Malar J. 2006;5:13. Epub 2006/02/17.View ArticlePubMedPubMed CentralGoogle Scholar
- Nmor JC, Sunahara T, Goto K, Futami K, Sonye G, Akweywa P, et al. Topographic models for predicting malaria vector breeding habitats: potential tools for vector control managers. Parasit Vectors. 2013;6:14.View ArticlePubMedPubMed CentralGoogle Scholar
- Burkett-Cadena ND, McClure CJW, Estep LK, Eubanks MD. Hosts or habitats: What drives the spatial distribution of mosquitoes? Ecosphere. 2013;4:1–16.View ArticleGoogle Scholar
- Gu W, Müller G, Schlein Y, Novak RJ, Beier JC. Natural plant sugar sources of Anopheles mosquitoes strongly impact malaria transmission potential. PLoS One. 2011;6(1):e15996.Google Scholar
- Andersson IH, Jaenson TGT. Nectar feeding by mosquitoes in Sweden, with special reference to Culex pipiens and Cx. torrentium. Med Vet Entomol. 1987;1:59–64.View ArticlePubMedGoogle Scholar
- Foster WA. Mosquito sugar feeding and reproductive energetics. An Rev Entomol. 1995;40:443–74.View ArticleGoogle Scholar
- Zhu G, Ju W, Chen JM, Liu Y. A Novel Moisture Adjusted Vegetation Index (MAVI) to reduce background reflectance and topographical effects on LAI retrieval. PLoS ONE. 2014;9(7):e102560. doi:10.1371/journal.pone.0102560.View ArticlePubMedPubMed CentralGoogle Scholar
- Ammar SE, Kenawy MA, Abdel-Rahman HA, Gad AM, Hamed AF. Ecology of the mosquito larvae in urban environments of Cairo Governorate, Egypt. J Egypt Soc Parasitol. 2012;42(1):191–202.Google Scholar
- SAS Institute Inc. 2013. Using JMP 11. Cary, NC: SAS Institute Inc.Google Scholar
- Phillips SJ, Anderson RP, Schapire RE. Maximum entropy modeling of species geographic distributions. Ecol Modelling. 2006;190(3-4):231–59.View ArticleGoogle Scholar
- Phillips S, Dudik M, Schapire R. A maximum entropy approach to species distribution modeling. In: Proceedings of the 21st International Conference on Machine Learning. New York: ACM Press; 2004.Google Scholar
- Phillips S, Dudik M. Modeling of species distributions with Maxent: new extensions and a comprehensive evaluation. Ecography. 2008;31:161–75.View ArticleGoogle Scholar
- Zhou G, Munga S, Minakawa N, Githeko A, Yan G. Spatial relationship between adult malaria vector abundance and environmental factors in western Kenya highlands. Am J Trop Med Hyg. 2007;77(1):29–35.PubMedGoogle Scholar
- Phillips SJ, Dudík M, Elith J, Graham CH, Lehmann A, Leathwick J, et al. Sample selection bias and presence-only distribution models: implications for background and pseudo-absence data. Ecol Applications. 2009;19(1):181–97.View ArticleGoogle Scholar
- Ruiz MO, Walker ED, Foster ES, Haramis LD, Kitron UD. Association of West Nile virus illness and urban landscapes in Chicago and Detroit. Int J Health Geogr. 2007;6:10. doi:10.1186/1476-072X-6-10.
- Rochlin I, Turbow D, Gomez F, Ninivaggi DV, Campbell SR. Predictive mapping of human risk for West Nile Virus (WNV) based on environmental and socioeconomic factors. PLoS ONE. 2011;6(8):e23280. doi:10.1371/journal.pone.0023280.View ArticlePubMedPubMed CentralGoogle Scholar
- Gibbs SEJ, Wimberly MC, Madden M, Masour J, Yabsley MJ, Stallnecht DE. Factors affecting the geographic distribution of West Nile virus in Georgia, USA: 2002–2004. Vector Borne Zoonotic Dis. 2006;6(1):73–82.View ArticlePubMedGoogle Scholar
- Saxena R, Das MK, Nagpal BN, Srivastava A, Gupta SK, Kumar A, et al. Identification of risk factors for malaria control by focused interventions in Ranchi district, Jharkhand, India. J Vector Borne Dis. 2014;51:276–81.PubMedGoogle Scholar
- Moore ID, Gessler PE, Nielsen GA, Peterson GA. Soil attribute prediction using terrain analysis. Soil Sc Soc Am J. 1993;57(2):443–52.Google Scholar
- Warren SD, Hohmann MG, Auerswald K, Mitasova H. An evaluation of methods to determine slope using digital elevation data. CATENA. 2004;58(3):215–33.View ArticleGoogle Scholar