Local impact of temperature and precipitation on West Nile virus infection in Culex species mosquitoes in northeast Illinois, USA
© Ruiz et al; licensee BioMed Central Ltd. 2010
Received: 29 October 2009
Accepted: 19 March 2010
Published: 19 March 2010
Models of the effects of environmental factors on West Nile virus disease risk have yielded conflicting outcomes. The role of precipitation has been especially difficult to discern from existing studies, due in part to habitat and behavior characteristics of specific vector species and because of differences in the temporal and spatial scales of the published studies. We used spatial and statistical modeling techniques to analyze and forecast fine scale spatial (2000 m grid) and temporal (weekly) patterns of West Nile virus mosquito infection relative to changing weather conditions in the urban landscape of the greater Chicago, Illinois, region for the years from 2004 to 2008.
Increased air temperature was the strongest temporal predictor of increased infection in Culex pipiens and Culex restuans mosquitoes, with cumulative high temperature differences being a key factor distinguishing years with higher mosquito infection and higher human illness rates from those with lower rates. Drier conditions in the spring followed by wetter conditions just prior to an increase in infection were factors in some but not all years. Overall, 80% of the weekly variation in mosquito infection was explained by prior weather conditions. Spatially, lower precipitation was the most important variable predicting stronger mosquito infection; precipitation and temperature alone could explain the pattern of spatial variability better than could other environmental variables (79% explained in the best model). Variables related to impervious surfaces and elevation differences were of modest importance in the spatial model.
Finely grained temporal and spatial patterns of precipitation and air temperature have a consistent and significant impact on the timing and location of increased mosquito infection in the northeastern Illinois study area. The use of local weather data at multiple monitoring locations and the integration of mosquito infection data from numerous sources across several years are important to the strength of the models presented. The other spatial environmental factors that tended to be important, including impervious surfaces and elevation measures, would mediate the effect of rainfall on soils and in urban catch basins. Changes in weather patterns with global climate change make it especially important to improve our ability to predict how inter-related local weather and environmental factors affect vectors and vector-borne disease risk.
Local impact of temperature and precipitation on West Nile virus infection in Culex species mosquitoes in northeast Illinois, USA.
West Nile virus (WNv), first reported in the United States in 1999, infects many species of birds as well as humans, equids, and other mammals. Though new to North America, the virus had circulated in Africa, Europe and the Middle East for some time prior to 1999; and other outbreaks, including those in southern Russia, Romania and Israel were indications of a change in the range of the virus [1–4]. Mosquitoes transmit the virus between hosts, with Culex species most often implicated as the primary amplification vector  and bridge vector [6, 7]. Many competent avian hosts have been identified both in the lab and field , and recent work in parts of North America have focused on the possible important role of American Robins (Turdus migratorius) in contributing to virus amplification and maintenance in the sylvatic cycle [9–11]. The continued risk for this sometimes severe and even fatal disease prompted establishment of annual surveillance programs of virus infection in mosquito populations. Now that several years have passed since the introduction of the WNv in North America, longitudinal data from testing of mosquitoes and host species (including records of human and equine illness) reported through systematic surveillance are available for development of models of the risk of infection. These records can be used to examine differences in infection between and within years and among locations to better understand the risk of transmission of the virus and to predict the possibility of place and time-specific outbreaks.
Though many published reports characterize associations between climatic and landscape factors and WNv occurrence, broad patterns have remained elusive as inconsistent results make generalization difficult. A review of 15 publications identified common landscape variables used to predict risk of WNv transmission, including distance to riparian corridor, vegetation measures, slope, elevation, human population numbers, housing and road density, type of urban land use, race, income, housing age, and host community structure [12–26]. Each predictor variable had a significant relationship with the WNv response variable in at least one study, but the directions of the effects were inconsistent among studies. From our review, we consider these inconsistencies to result largely from differences in climatic factors and mosquito vectors among geographic locations, heterogeneous temporal and spatial resolution of the analyses, and the different response variables used to measure risk.
In studies focused on spatial or temporal risk of illness from WNv, the response variable typically includes data from human disease cases [13, 14, 16, 20, 27–29], vector abundance or infection [15, 17, 30–34], evidence of infection in non-human hosts [21–24, 26], or a combination of vector and host data [1, 12, 25, 35–37]. When human case data are the response variable, inconsistencies between WNv risk and environmental factors have been attributed to: exposure occurring away from the home residence; time outdoors; use of insect repellent; socioeconomic differences, and even variable case definition [38–41]. An analysis in Indianapolis, IN, with similar climatic features and spatial scale as the current study, identified a relationship with both low precipitation and higher temperatures and increased numbers of human WNV cases .
For analysis of the transmission of mosquito-borne arboviruses, vector abundance and infection data are key elements, and the two are not always proportional in space or time [42, 43]. For this reason, the fine scale patterns relevant to production and dispersal of adult female mosquitoes and the association between WNv infection in mosquitoes and weather in the context of other landscape features are in need of further investigation. When analyses include different vectors, results will often be inconsistent [16, 29]. For example, the population size of Culex pipiens., the primary enzootic and epidemic vector in the eastern U.S. north of 36 degrees latitude, is often impacted negatively by large rain events due to the flushing of catch basins, a primary urban larval habitat, and the reduction of organic content in all ovipositing sites [44, 45]. By contrast, the vector Culex tarsalis generally responds positively to heavy precipitation, which provides the typical larval habitat in rural areas in the western part of the U.S [28, 33]. In semi-permanent wetlands, drought conditions can increase abundance of some vector populations as they result in more larval breeding sites with fewer competitors and mosquito predators . Early season drought with subsequent wetting and low water table depth preceded amplification episodes for both WNv and St. Louis encephalitis virus in peninsular Florida, where Culex nigripalpus functions as the main vector [47–49].
Increased temperature is known to increase growth rates of vector populations , decrease the length of the gonotrophic cycle (interval between blood meals), shorten the extrinsic incubation period of the virus in the vector and increase the rate of virus evolution [50–54]. Kunkel et al.  showed a correlation between the number of days when daily maximum temperature exceeded a threshold (degree days), timing of a seasonal shift to a higher proportion of Culex pipiens among all Culex species, and the onset of the amplification phase of WNv transmission seasonally in Illinois. Other studies considered less proximate weather conditions, such as increased rainfall in the preceding year . Warmer winter temperatures and warmer March and April may lead to larger summer mosquito populations . Temperature has also been linked to the rate of evolution of the virus and warmer temperatures facilitated the displacement of the WNv NY99 genotype by the WN02 genotype . Bertolotti et al.  discovered high genetic variation of WNv at fine temporal and spatial scales, with variation in local temperature offered as one explanation for it. At the same time, in one study, very high temperatures (above 30°C) reduced larval Culex tarsalis survival .
Finally, the choice of geographic scale and units can impact profoundly the outcome of an analysis. Coarse geographic scales such as the county level can obscure fine spatial patterns in a heterogeneous landscape . The census tract provides a somewhat finer scale geographic unit, but while logistically useful, is not a biological meaningful unit. In fact, any use of a spatial unit introduces the modifiable area unit problem [59, 60], where spatial units of analysis become arbitrary or possibly even introduce systematic bias. In addition to the importance of fine spatial scales, coarse temporal scales could obscure some temporal patterns [31, 32].
Our analysis focuses on greater Chicago, Illinois, where WNv infection in mosquitoes, horses and birds was first noted in 2001, and where human illness has been reported every year from 2002 to 2009 [, http://www.idph.state.il.us/envhealth/wnv.htm]. We consider temporal and spatial patterns of infections in mosquitoes for the years 2004 to 2008, years for which comprehensive mosquito testing data were available in Illinois. We have focused, in particular, on the meteorological conditions that precede or are concurrent with amplification, when a sharp increase is seen in the infection rate in mosquito populations. Our three research questions are: 1) Inter-annually: what are the conditions associated with higher mosquito infection in some years compared to others? 2) Intra-annually: what temporal characteristics of rainfall and temperature precede changes in mosquito infection and with what temporal lag? 3) Spatially: can the patterns of rainfall and temperature help explain the differences in mosquito infection across space? We used surveillance data from the Illinois Department of Public Health (IDPH) and publicly available meteorological readings, and consider the heterogeneity of urban land cover through an analysis of digital spatial data to identify and forecast favorable conditions for WNv amplification in the greater Chicago area.
Weather data from archived sources from the two counties and five surrounding counties were retrieved from three sources (Figure 1). The NOAA NESDIS National Climatic Data Center http://cdo.ncdc.noaa.gov/CDO/cdo included from 21-32 stations, depending on the year, and provided both temperature and precipitation. From the U.S. Geological Survey, we retrieved data from 18 to 41 stations per year, all reporting precipitation http://pubs.usgs.gov/wdr/2005/wdr-il-05/start.htm. Additional precipitation data were obtained from a network of 25 stations maintained by the Illinois State Water Survey http://www.isws.illinois.edu/atmos/ccprecipnet.
In addition to the single weeks' precipitation measures, the weekly precipitation was smoothed using a 3-week and a 5-week moving average. These three different weekly measures of precipitation were included to provide several levels of smoothing for this variable.
Other environmental factors and geographic units
For spatial analysis, we created a 2000 m hexagon grid encompassing the two county study area and summarized all data for those units (N = 1,263). The hexagon size was chosen to equal about the same number of units as census tracts in the same area, but in contrast to tracts are of uniform size, providing a more neutral landscape unit. For mosquito infection, precipitation and temperature, we used a local moving average method of geographic interpolation [63, 64] with inverse distance weighting (IDW) from the six closest points to create weekly GIS raster grid (100 m grid size) maps for each week and for each variable. The "Z" values (also known as the support) for the IDW interpolation were based on the calculation of the MIR at the trap locations or the weather variables at the weather stations for which data were available for that week. For this step, we included data from five counties bordering Cook and DuPage counties to better estimate values at the edge of the study area. In this way, we estimated values for each raster grid in the study region. The raster grid-based data were then summarized for the 2000 m hexagons with ArcGIS using the average of all raster grid cells within each single hexagonal geographic unit.
We developed a spatial database of environmental factors included previously in an analysis of human illness risk in this area [13, 65] to assess the spatial differences in MIR relative to both weather-related and other environmental factors. While the MIR and the weather data were summarized by week, other environmental factors did not have a temporal component and their values were static. To account for the social environment, we included the number of people in residence in each hexagon, the percentage of the population defined racially as "Black" in the U.S. census, and the percentage of housing from different decades . We suggest that both mosquito control and mosquito production vary by neighborhood due both to engineering and housing development differences related to age of housing and socio-cultural differences that affect vegetation choices and landscape patterns. Terrain characteristics were measured by the minimum, maximum, range, mean and standard deviation of the elevation from a digital elevation model http://seamless.usgs.gov. These terrain variables provided several ways to measure the degree to which water would flow or remain in an area. For land cover, we measured the mean and standard deviation of imperviousness (landcover.usgs.gov), and the percentage and standard deviation of percentage vegetation (per hexagon) from a digital land cover map http://www.agr.state.il.us/gis/landcover99-00.html. These final variables accounted for the amount of vegetation and the amount of surfaces where soil was covered by pavement. We expected those to be related to local mosquito and bird habitat.
We performed correlations and linear regression as well as classification and regression tree analysis using the R statistical package . For the first and second research questions, we considered the difference of the weekly MIR from the average MIR at that week across years, as well as the difference in temperature and rainfall from their respective long term means. In this way, we accounted for the patterns that were outside the dominant annual cycle of solar radiation and the seasonality component, which remains unchanged in every year, and focused on the components deviating from this fixed annual cycle. To calculate these differences, the expected, or normal, value of mosquito infection for a given week in the Chicago area was needed. It is not easily known since WNv has been in the area for a relatively short time and there can be instability of MIR results . For MIR, we defined the "normal weekly value" as the average for each week in the combined years 2004 to 2007 as measured from the hexagon level dataset. For temperature and precipitation, the 30-year normal weather values are an accepted standard, but we also wanted to consider the specific differences and possible warmer trends of the more recent years in question. Thus, we defined normal weekly weather values both as the average weather values between the years 2004 and 2007, and the 30-year normal values. We then compared the two results. The 30-year normal values were detrended by accounting for the linear change over the years.
Time-lagged correlations were used to explore the association between mosquito infection and weekly temperature and rainfall at a range of temporal lags. Linear regression models were then used to estimate the empirical models that best fit the MIR with weather variables at key temporal lags selected based on the exploratory correlation analysis. Weekly weather variables were considered starting in week 1, while for MIR, we focused on the weeks from 18 to 35 (from about the end of April to the end of August), during the period of increasingly stronger activity and up to the peak values of MIR. We developed linear models with and without an autoregressive term for MIR, so that if MIR was already known, then the prior week's MIR could be included in the estimate for the following weeks. We developed several linear models and used them to reconstruct the weekly pattern of MIR for the years 2004 to 2007. These models were used then to predict the weekly pattern for 2008, thus allowing a test of the model against a year not included in the model development.
For the third research question, we used regression tree analysis (RT) and Random Forests (RF) with weekly temperature and precipitation and the other environmental factors as measured in the hexagon grid. With RT analysis, a set of rules is developed from the independent (or predictor) variables that can best recreate the observed pattern in the response variable . The response variable in this analysis was the average MIR during the period of peak infection, weeks 32 to 34. In this technique the variability of the response variable is partitioned along binary nodes of the predictive covariates that will lead to an average value of the response variable. Nodes of covariates can be nested, with the most basal explaining most of the variability in the response variable. This technique has the advantage over traditional generalized linear models of capturing non-linear relationships between the covariates and response and has no assumptions about spatial or temporal autocorrelations. It has been used to analyze data in other mosquito-borne disease systems . An additional advantage of this technique is that it doesn't rely on the assumptions that are required for parametric statistics and the analysis is not restricted by linearity in predictor or response variables or by multicollinearity in predictor variables. Trees were selected using the "cost-complexity" algorithm, where auxiliary nodes are cut if no significant loss in the mean square error of the predictions is detected. This technique helps to avoid over-parameterization.
RF is a tool for prediction based on several regression trees. Basically, a new response is constructed using the fit of a decision tree using sampling with replacement. A new decision tree is built, and predictions are based on the aggregate outcome of the trees forming the random forest . Here we used fully cross-validated regression trees, which involves leaving out one of the observations when fitting a tree or a forest and measuring the error (difference) between the observed and the predicted value for the observation not used for the fit. Random forests with 1000 trees were used to predict average MIR during the period of high infection, in weeks 32 to 34, for the years 2004, 2005, 2006, and 2007. We analyzed all years together and separately, and used as covariates: 1) all weekly weather and other environmental factors, 2) only weather variables, and 3) only other environmental factors.
Results and Discussion
Research Question 1) Inter-annually: what are the conditions associated with higher mosquito infection in some years compared to others?
Descriptive overview of mosquito infection and weather
West Nile virus infection and weather for Cook and DuPage county from 2002 to 2008.
Average mosquito infection
Percent Culex mosquito pools positive
Number human cases illness
Annual precipitation in cm
(difference from average)
Degree weeks at week 35 base = 22 C
Research Question 2) Intra-annually: what temporal characteristics of rainfall and temperature precede changes in mosquito infection and with what temporal lag?
Temporal model - exploratory phase
Temporal model -regression model development phase
Using the variables indicated from this exploratory analysis, we developed linear regression models to determine the best-fit model for all of the years. We used these to reconstruct the weekly changes in MIR for the years 2004 to 2007 and to simulate weekly changes for the year 2008, which was not included in the model development. We also simulated MIR for the years 2002 and 2003, years for which comparable mosquito data were not available, to learn what patterns would be estimated for a year with high human illness (2002) and a year with few cases (2003). Note that both the independent and dependent variables in the regression models are the differences from the corresponding averages between 2004 and 2007 but for simulated values, the averages were added back in.
The candidate variables for this set of models, based on the exploratory phase, were DWs accumulated from the beginning of the year up to one week prior to the MIR measure (DW), the 5-week-average precipitation 11 weeks prior (prcp5wk_lag11) and the 3-week average precipitation 3 weeks prior (prcp3wk_lag3). Further, considering that in other studies the human WNv incidence was associated with annual precipitation from the preceding year, we included that variable as one of candidate variables (prcp_annual). Finally, we found that the MIR itself is a first order auto-regressive (AR) process, so to simulate the MIR for any specific week, the MIR measured in the previous week was also included as one of the possible explanatory factors. After developing the best-fit models, we also created models that were somewhat less strong, but would allow for prediction earlier than the best-fit models. For these, we considered models both with and without the AR term.
Coefficients and variables from temporal regression models.
prcp3wk (3 wk lag)
prcp3wk (4 wk lag)
prcp5wk (11 wk lag)
prcp_annual (prior year)
DW (1 wk lag)
DW (4 wk lag)
DWC (1 wk lag)
For Models 1 and 2, we found that the precipitation variable prcp5wk_lag11 was not selected as an explanatory factor, though it had been found to be strongly negatively correlated with the MIR based on the correlations. This can be explained by its co-linearity with the DW variable (r = -0.7). Specifically, years with a dry spring season also tended to be hot during that period (Additional File 2: Temporal, Figure A). Thus, adding prcp5wk_lag11 did not significantly improve the model.
As noted above, the MIR in our study area peaked around week 34, after which it decreased rapidly until reaching zero at around week 43. Although the year-to-year contrasts during this declining period were not as striking as those in the growth period, we still found some inter-annual variation. For example, MIR decreased at a lower rate in the fall of 2006 and 2007 than in other years. We found these differences can be explained by the higher than normal temperature during that period in those years. Thus we constructed a linear regression model (Model C in Table 2) for the MIR in its declining phase using the cooling DW (DWC) and the AR term.
Forecasting the MIR temporally
An early warning system for mosquito infection from WNv will be most useful when it is possible to estimate mosquito infection as far in advance as possible. Models 3 and 4 were developed for this purpose. Here, earlier meteorological information was used to forecast the MIR to facilitate decisions on corresponding preventive measures (Table 2). While the highest correlation between the accumulated DW and MIR was at a lag of 1 week, the DW value cumulated to 4 weeks prior was also strongly correlated, at around 0.7. For Model 3, we included the accumulated DW value leading by 4 weeks, the AR term with a lag of 2 weeks and the precipitation from the prior year as forecasting factors for MIR. In the event that only meteorological variables were available and no prior MIR was included, then the 3-week average precipitation at a 4-week lag, the accumulated DW at a 4- week lag and the prior year precipitation provides a less strong but adequate forecasting model as seen in Model 4 (Table 2 & Figure 5B).
Sensitivity analysis of the four models based on both leave-one-out cross validation and an approach using 5-fold cross validation both indicated that the predicted R2 was nearly identical to the model R2. Model residuals were found to be nearly normal in Model 1 and without temporal autocorrelation. Model 1 included a 1st order auto-regressive term on MIR (Additional File 2: Temporal, Figure B), but the pattern of residuals for other models without that auto-regressive term did have evidence of the need for this term. The MIR for a prior time period may be necessary to build robust weekly models. The differences seen in the model diagnostics between years highlights the difficulties in finding general trends across only a few years, in a region where weekly weather patterns can diverge considerably from year to year.
Research Question 3) Spatially: can the patterns of rainfall and temperature help explain the differences in mosquito infection across space
Spatial patterns of MIR
Random forest spatial model results.
Variables Tested and Model Type
Most important variables in model
All variables Random Forest
Precip W33, W19, W28
Precip W23, W29, W33
Precip W28, W21
Precip W29, W18, W19, W33
Precip W15, W28
Weather variables only Random Forest
Precip W33, W19
Precip W33, W23, W29
Precip W32, W28, W19
Precip W33, W29, W19
Precip W15, W21
Non-weather variables only. Random Forest
Pct Black, Human population, Elevation range
Pct Black, Maximum elevation, Minimum elevation
% pre-40's housing, Pct Black, % 50's housing
Pct Black, maximum elevation, % 90's housing
Pct Black, maximum elevation, mean elevation
Regression tree spatial model results.
Variables Tested and Model Type
Most important variables in model
All variables Regression Tree
Precip W17, W19
Precip W29, W28
Precip W24, W26
Temp W35, W27
Temp W17, W22
Weather variables only Regression Tree
Precip W17, W19
Precip W28, W29
Precip W24, W26
Temp W35, W27
Temp W17, W22
Non-weather variables only Regression Tree
Impervious surface, % pre-40's housing, % 50's housing
Minimum elevation, Mean elevation, impervious surface
Maximum elevation, impervious surface, Human population
Maximum elevation, impervious surface, Human population
Maximum elevation, impervious surface, Human population
Finely grained temporal and spatial patterns of precipitation and air temperature have a significant impact on the timing and location of increased mosquito infection across the Chicago study area. Temperature, in particular, mediates the magnitude and timing of increased MIR within a season and it is strongly indicated as a key factor for explaining much of the observable differences between years. The effect of increased temperature on MIR is especially strong within a week. This is consistent with expectations and would suggest increased mosquito productivity through a shortened time for development of the Culex vector and increased rates of oviposition activity , which implies increased biting rates for anautogenous populations. Further, more quickly replicating virus spurred by higher temperatures precedes an increase in mosquito infection. The patterns of association between MIR and precipitation for 2004, 2005, and 2006 are similar to the patterns observed in Florida, where periods of dry-downs and wetting events lead to elevated SLE and WNV infections [48, 49]. In Illinois, however, 2007 demonstrated a strong departure from this pattern and further emphasizes the more important role of temperature. Additional years of data on mosquito infection and weather need to be assessed to better determine the dominant patterns.
The lag seen in the effect from rainfall and differences in this effect among years suggests that the relationship with rainfall and MIR is more complex and more variable than the effects of temperature. The observed lower rainfall about eleven weeks prior to the highest MIR weeks (32 to 34) indicates that the processes in play may be related to subtle changes in soil moisture, near ground humidity, and an increase in desirable oviposition sites that would result from an increase in organic matter in standing water. While dry years were more likely to result in higher MIR, both higher temperatures and less rainfall in tandem were more highly associated with the highest mosquito infection rates and/or with more human cases. This was seen in particular in the years 2002 (only human illness data available) and 2005 (evidence from human illness and mosquito infection).
WNV infection in mosquitoes was negatively correlated with the previous year's precipitation. This is contrary to the findings of Landesman et al.  for the eastern part of the United States, and consistent with those in the west. The mechanisms for the patterns between past weather conditions to influence arbovirus transmission are still unclear. Our results support the importance of using off-season or previous year weather data in prediction models, but details on this must be left to a future analysis with a longer time series . We also suggest that it is worthwhile to consider winter temperature conditions, which might lead to a larger number of overwintering adult Culex spp. mosquitoes and increase the potential for amplification the following year due to larger initial vector populations.
Compared to the temporal patterns, the spatial analysis indicated the relative importance of rainfall, with temperature playing a more muted role. The results for the 2004 random forests for example (the strongest of the models) show that the three most important variables all measured precipitation. The more consistently seen measures of importance indicated that drier and warmer conditions tended to result in higher MIR, but many exceptions also occurred, emphasizing the complex nature of the spatial patterns. The spatial models for a single year were consistently stronger than those where data from all years were combined. From a public health perspective, this would indicate that the infection risk changes from year to year as weather patterns vary. This is consistent with the outcome of the temporal models, which also exhibited different patterns by year. At the same time, the spatial models for all years were modestly robust. From this, we conclude that in balance, in future efforts to determine factors that influence MIR, it is more important to focus on the dynamic weather variables, possibly supplemented by more dynamic measures of vegetation and surface moisture, rather than the more static landscape metrics.
In terms of the environmental features that were found to influence MIR, it is worth noting that the more important non-weather variables were the proportion of impervious surfaces and elevation variables. Both of these factors would mediate the effect of rainfall on soils and vegetation and in catch basins. For example, low-lying landscape patches with some impervious surfaces could be expected to be more productive for mosquitoes, since they can accumulate water and have resources necessary for mosquitoes to thrive. The importance of the spatial variable that measures the percentage of the population identified as "Black" by the U.S. Census Bureau was also notable, as this was a variable of note in prior analyses in this region, which focused on human case data [13, 65]. Though human case data would be influenced by possible biases due to differences in reporting rates in different neighborhoods, mosquito infection rates are a more neutral measure of presence of the virus. While the present analysis did not allow us to fully explore this relationship, the persistence of the statistical association in the Chicago area warrants further investigation.
We further suggest that a finer grained spatial measure of temperature than that used here may also reveal that temperature plays a role that was not detected with the resolution of temperature data available. Especially, given the strong importance of temperature in the temporal models, differences in surface temperatures due to the amount of vegetation relative to the built environment may be more important than what was revealed in the current analysis. This can be addressed by using remote sensing techniques to obtain more complete continuous coverage of temperature than what is possible with data from point-based weather stations.
The ability of the weekly models to simulate MIR based on prior precipitation and temperature and prior MIR is moderately strong and indicates that at least on a weekly basis, the amplification of the WNv in mosquitoes can be forecast for this area when the level of mosquito infection from the prior week is available for the model. However, given the practical and statistical difficulties in using fine-scale time series data, a coarser temporal scale may be more realistic from a public health standpoint.
Local forecasts available to public health personnel would need to have parameters specific to the place of interest. By way of comparison to our own results, increased temperatures in a Culex pipiens- driven system in Israel resulted in a significant two-week lag for mosquito abundance and about a five week lag for human cases of WNV . Maximum urban and grass temperature were strong predictors of human cases of WNV in the Chicago, Illinois region and low precipitation and warm temperature were associated with WNV cases in Indianapolis, Indiana . Each of these examples further emphasizes the need for models that are place-specific. The framework for developing early warning systems needs to be standardized but geographic variability makes it necessary to customize the system for each location using region specific patterns .
Just as with other analyses in other locations, the exact nature of the environmental factors most associated with higher mosquito infection were not consistently and clearly delineated, and further studies on mosquito ecology across finely grained heterogeneous landscapes are urged to fully understand the interaction of weather and landscape in shaping mosquito population and infection dynamics. This is especially important in light of changing weather patterns associated with global climate change, as prior stability cannot be assumed. The use of local weather data at multiple locations and the integration of mosquito infection data from several sources across multiple years are important to the strength of the models presented. They allowed us to look more specifically at spatial and temporal scales that are more in keeping with the scale of virus amplification than prior studies. The ability to provide more precise spatial risk estimates is still limited, but the modeling approach presented here will help to define more clearly hypotheses that can be tested relative to the biological processes that are most closely related to the sylvatic transmission cycle.
This project was funded through NSF Ecology of Infectious Disease Grants 0429124 (West Nile Virus: Eco-Epidemiology of Disease Emergence in Urban Areas) and 0840403 (Eco-epidemiology of Disease Emergence in Urban Areas II) and through funding from the University of Illinois Adaptive Infrastructure Information Systems Initiative. We thank the Village of Oak Lawn for ongoing support of this project. Kenneth Kunkel contributed to early discussions on this topic. Jane Messina, Julie Clennon and Joanna Ganning helped with geocoding and organization of spatial mosquito data. . LFC and UK also thank the Research and Policy for Infectious Disease Dynamics (RAPIDD) program of the Science and Technology Directorate and the Fogarty International Center, National Institutes of Health program on mosquito borne diseases, for helpful discussions in the context of this manuscript. We thank Tavis Anderson and two anonymous reviewers for helpful suggestions for revisions, which have greatly improved the paper.
- Platonov AE, Fedorova MV, Karan LS, Shopenskaya TA, Platonova OV, Zhuravlev VI: Epidemiology of West Nile infection in Volgograd, Russia, in relation to climate change and mosquito (Diptera: Culicidae) bionomics. Parasitology Research. 2008, 103: S45-S53. 10.1007/s00436-008-1050-0.View ArticlePubMedGoogle Scholar
- Han LL, Popovici F, Alexander JP, Laurentia V, Tengelsen LA, Cernescu C, Gary HE, Ion-Nedelcu N, Campbell GL, Tsai TF: Risk factors for West Nile virus infection and meningoencephalitis, Romania, 1996. Journal of Infectious Diseases. 1999, 179: 230-233. 10.1086/314566.View ArticlePubMedGoogle Scholar
- Hayes EB, Komar N, Nasci RS, Montgomery SP, O'Leary DR, Campbell GL: Epidemiology and transmission dynamics of West Nile Virus disease. Emerging Infectious Diseases. 2005, 11: 1167-1173.PubMed CentralView ArticlePubMedGoogle Scholar
- Petersen LR, Roehrig JT: West Nile virus: A reemerging global pathogen. Emerg Infect Dis. 2001, 7: 611-614. 10.3201/eid0704.010401.PubMed CentralView ArticlePubMedGoogle Scholar
- Turell MJ, Dohm DJ, Sardelis MR, O Guinn 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. 10.1603/0022-2585(2005)042[0057:AUOTPO]2.0.CO;2.View ArticlePubMedGoogle Scholar
- Kilpatrick AM, Kramer LD, Campbell SR, Alleyne EO, Dobson AP, Daszak P: West Nile virus risk assessment and the bridge vector paradigm. Emerg Infect Dis. 2005, 11: 425-429.PubMed CentralView ArticlePubMedGoogle Scholar
- Hamer GL, Kitron UD, Brawn JD, Loss SR, Ruiz MO, Goldberg TL, Walker ED: Culex pipiens (Diptera: Culicidae): a bridge vector of West Nile virus to humans. J Med Entomol. 2008, 45: 125-128. 10.1603/0022-2585(2008)45[125:CPDCAB]2.0.CO;2.View ArticlePubMedGoogle Scholar
- Kilpatrick AM, LaDeau SL, Marra PP: Ecology of west nile virus transmission and its impact on birds in the western hemisphere. Auk. 2007, 124: 1121-1136. 10.1642/0004-8038(2007)124[1121:EOWNVT]2.0.CO;2.View ArticleGoogle Scholar
- Kent R, Juliusson L, Weissmann M, Evans S, Komar N: Seasonal Blood-Feeding Behavior of Culex tarsalis (Diptera: Culicidae) in Weld County, Colorado, 2007. J Med Entomol. 2009, 46: 380-390. 10.1603/033.046.0226.View ArticlePubMedGoogle Scholar
- Kilpatrick AM, Daszak P, Jones MJ, Marra PP, Kramer LD: Host heterogeneity dominates West Nile virus transmission. Proc R Soc Lond Ser B-Biol Sci. 2006, 273: 2327-2333. 10.1098/rspb.2006.3575.View ArticleGoogle Scholar
- Hamer GL, Kitron UD, Goldberg TL, Brawn JD, Loss SR, Ruiz MO, Hayes DB, Walker ED: Host Selection by Culex pipiens Mosquitoes and West Nile Virus Amplification. Am J Trop Med Hyg. 2009, 80: 268-278.PubMedGoogle Scholar
- Loss SR, Hamer GL, Walker ED, Ruiz MO, Goldberg TL, Kitron UD, Brawn JD: Avian host community structure and prevalence of West Nile virus in Chicago, Illinois. Oecologia. 2009, 159: 415-424. 10.1007/s00442-008-1224-6.View ArticlePubMedGoogle Scholar
- Ruiz MO, Tedesco C, McTighe TJ, Austin C, U K: Environmental and social determinants of human risk during a West Nile virus outbreak in the greater Chicago area, 2002. International Journal of Health Geographics. 2004, 3: 11-10.1186/1476-072X-3-8.View ArticleGoogle 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: 1539-1545. 10.3201/eid1410.071396.PubMed CentralView 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: 197-206. 10.1089/vbz.2007.0154.View ArticlePubMedGoogle 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. International Journal of Health Geographics. 2008, 7: 19-PubMed CentralView ArticlePubMedGoogle Scholar
- Diuk-Wasser MA, Brown HE, Andreadis TG, Fish D: Modeling the spatial distribution of mosquito vectors for West Nile virus in Connecticut, USA. Vector-Borne Zoonotic Dis. 2006, 6: 283-295. 10.1089/vbz.2006.6.283.View ArticlePubMedGoogle Scholar
- Gu WD, Lampman R, Krasavin N, Berry R, Novak R: Spatio-temporal analyses of West Nile virus transmission in Culex mosquitoes in Northern Illinois, USA, 2004. Vector-Borne Zoonotic Dis. 2006, 6: 91-98. 10.1089/vbz.2006.6.91.View ArticlePubMedGoogle Scholar
- Ozdenerol E, Bialkowska-Jelinska E, Taff GN: Locating suitable habitats for West Nile Virus-infected mosquitoes through association of environmental characteristics with infected mosquito locations: a case study in Shelby County, Tennessee. International Journal of Health Geographics. 2008, 7: 12-10.1186/1476-072X-7-12.PubMed CentralView ArticlePubMedGoogle Scholar
- Brownstein JS, Rosen H, Purdy D, Miller JR, Merlino M, Mostashari F, Fish D: Spatial analysis of West Nile Virus: Rapid risk assessment of an introduced Vector-Borne Zoonosis (vol 2, pg 157, 2002). Vector-Borne Zoonotic Dis. 2003, 3: 155-155. 10.1089/153036603768395861.View ArticleGoogle Scholar
- Cooke WH, Grala K, Wallis RC: Avian GIS models signial human risk for West Nile virus in Mississippi. International Journal of Health Geographics. 2006, 5: 36-10.1186/1476-072X-5-36.PubMed CentralView ArticlePubMedGoogle Scholar
- Pradier S, Leblond A, Durand B: Land cover, landscape structure, and West Nile virus circulation in southern France. Vector-Borne Zoonotic Dis. 2008, 8: 253-263. 10.1089/vbz.2007.0178.View ArticlePubMedGoogle Scholar
- Bradley CA, Gibbs SEJ, Altizer S: Urban land use predicts West Nile Virus exposure in songbirds. Ecol Appl. 2008, 18: 1083-1092. 10.1890/07-0822.1.View ArticlePubMedGoogle Scholar
- Gibbs SEJ, Allison AB, Yabsley MJ, Mead DG, Wilcox BR, Stallknecht DE: West Nile virus antibodies in avian species of Georgia, USA: 2000-2004. Vector-Borne Zoonotic Dis. 2006, 6: 57-72. 10.1089/vbz.2006.6.57.View ArticlePubMedGoogle Scholar
- Allan BF, Langerhans RB, Ryberg WA, Landesman WJ, Griffin NW, Katz RS, Oberle BJ, Schutzenhofer MR, Smyth KN, St. Maurice Ad, Clark L, Crooks KR, Hernandez DE, McLean RG, Ostfeld RS, Chase JM: Ecological correlates of risk and incidence of West Nile virus in the United States. Oecologia. 2009, 158: 699-708. 10.1007/s00442-008-1169-9.View ArticlePubMedGoogle Scholar
- Gomez A, Kilpatrick AM, Kramer LD, Dupuis AP, Maffei JG, Goetz SJ, Marra PP, Daszak P, Aguirre AA: Land use and West Nile Virus seroprevalence in wild mammals. Emerg Infect Dis. 2008, 14: 962-965. 10.3201/eid1406.070352.PubMed CentralView ArticlePubMedGoogle Scholar
- Ruiz MO, Walker ED, Foster ES, Haramis LD, Kitron UD: Association of West Nile virus illness and urban landscapes in Chicago and Detroit. International Journal of Health Geographics. 2007, 6: 10-10.1186/1476-072X-6-10.PubMed CentralView ArticlePubMedGoogle Scholar
- Wimberly MC, Hildreth MB, Boyte SP, Lindquist E, Kightlinger L: Ecological Niche of the 2003 West Nile Virus Epidemic in the Northern Great Plains of the United States. PLoS One. 2008, 3: Article No.: e3744-10.1371/journal.pone.0003744.View ArticlePubMedGoogle Scholar
- Soverow JE, Wellenius GA, Fisman DN, Mittleman MA: Infectious Disease in a Warming World: How Weather Influenced West Nile Virus in the United States (2001-2005). Environmental Health Perspectives. 2009, 117: 1049-1052.PubMed CentralView ArticlePubMedGoogle Scholar
- DeGroote J, Mercer DR, Fisher J, Sugumaran R: Spatiotemporal investigation of adult mosquito (Diptera: culicidae) Populations in an eastern Iowa county, USA. J Med Entomol. 2007, 44: 1139-1150. 10.1603/0022-2585(2007)44[1139:SIOAMD]2.0.CO;2.View ArticlePubMedGoogle Scholar
- Trawinski PR, Mackay DS: Meteorologically conditioned time-series predictions of West Nile virus vector mosquitoes. Vector-Borne Zoonotic Dis. 2008, 8: 505-521. 10.1089/vbz.2007.0202.View ArticlePubMedGoogle Scholar
- Walsh AS, Glass GE, Lesser CR, Curriero FC: Predicting seasonal abundance of mosquitoes based on off-season meteorological conditions. Environmental and Ecological Statistics. 2008, 15: 279-291. 10.1007/s10651-007-0056-6.View ArticleGoogle Scholar
- Reisen WK, Cayan D, Tyree M, Barker CA, Eldridge B, Dettinger M: Impact of climate variation on mosquito abundance in California. J Vector Ecol. 2008, 33: 89-98. 10.3376/1081-1710(2008)33[89:IOCVOM]2.0.CO;2.View ArticlePubMedGoogle Scholar
- Rochlin I, Harding K, Ginsberg HS, Campbell SR: Comparative analysis of distribution and abundance of West Nile and eastern equine encephalomyelitis virus vectors in Suffolk County, New York, using human population density and land use/cover data. J Med Entomol. 2008, 45: 563-571. 10.1603/0022-2585(2008)45[563:CAODAA]2.0.CO;2.View ArticlePubMedGoogle Scholar
- Paz S, Albersheim I: Influence of warming tendency on Culex pipiens population abundance and on the probability of West Nile Fever outbreaks (Israeli case study: 2001-2005). EcoHealth. 2008, 5: 40-48. 10.1007/s10393-007-0150-0.View ArticlePubMedGoogle Scholar
- Liu H, Weng QH, Gaines D: Spatio-temporal analysis of the relationship between WNV dissemination and environmental variables in Indianapolis, USA. International Journal of Health Geographics. 2008, 7: 66-10.1186/1476-072X-7-66.PubMed CentralView ArticlePubMedGoogle Scholar
- Winters AM, Bolling BG, Beaty BJ, Blair CD, Eisen RJ, Meyer AM, Pape WJ, Moore CG, Eisen L: Combining mosquito vector and human disease data for improved assessment of spatial west nile virus disease risk. Am J Trop Med Hyg. 2008, 78: 654-665.PubMedGoogle Scholar
- Eisen L, Eisen RJ: Need for improved methods to collect and present spatial epidemiologic data for vectorborne diseases. Emerg Infect Dis. 2007, 13: 1816-1820.PubMed CentralView ArticlePubMedGoogle Scholar
- Ostfeld RS, Glass GE, Keesing F: Spatial epidemiology: an emerging (or re-emerging) discipline. Trends Ecol Evol. 2005, 20: 328-336. 10.1016/j.tree.2005.03.009.View ArticlePubMedGoogle Scholar
- de la C Sierra B, Kouri G, Guzman MG: Race: a risk factor for dengue hemorrhagic fever. Archives of Virology. 2007, 152: 533-542. 10.1007/s00705-006-0869-x.View ArticlePubMedGoogle Scholar
- Gujral IB, Zielinski-Gutierrez EC, LeBailly A, Nasci R: Behavioral risks for West Nile Virus disease, northern Colorado, 2003. Emerg Infect Dis. 2007, 13: 419-425. 10.3201/eid1303.060941.PubMed CentralView ArticlePubMedGoogle Scholar
- Andreadis TG, Anderson JF, Vossbrinck CR, Main AJ: Epidemiology of West Nile virus in Connecticut: A five-year analysis of mosquito data 1999-2003. Vector-Borne Zoonotic Dis. 2004, 4: 360-378. 10.1089/vbz.2004.4.360.View ArticlePubMedGoogle Scholar
- Hamer GL, Walker ED, Brawn JD, Loss SR, Ruiz MO, Goldberg TL, Schotthoefer AM, Brown WM, Wheeler E, Kitron UD: Rapid amplification of West Nile virus: the role of hatch-year birds. Vector-Borne Zoonotic Dis. 2008, 8: 57-67. 10.1089/vbz.2007.0123.View ArticlePubMedGoogle Scholar
- Koenraadt CJM, Harrington LC: Flushing effect of rain on container-inhabiting mosquitoes Aedes aegypti and Culex pipiens (Diptera: Culicidae). J Med Entomol. 2008, 45: 28-35. 10.1603/0022-2585(2008)45[28:FEOROC]2.0.CO;2.View ArticlePubMedGoogle Scholar
- Geery PR, Holub RE: Seasonal abundance and control of Culex spp. in catch basins in Illinois. J Am Mosq Control Assoc. 1989, 5: 537-540.PubMedGoogle Scholar
- Chase JM, Knight TM: Drought-induced mosquito outbreaks in wetlands. Ecol Lett. 2003, 6: 1017-1024. 10.1046/j.1461-0248.2003.00533.x.View ArticleGoogle Scholar
- Shaman J, Day JF, Stieglitz M: Drought-induced amplification of Saint Louis encephalitis virus, Florida. Emerg Infect Dis. 2002, 8: 575-580. 10.3201/eid0801.010049.PubMed CentralView ArticlePubMedGoogle Scholar
- Shaman J, Day JF, Stieglitz M: Drought-induced amplification and epidemic transmission of West Nile Virus in southern Florida. J Med Entomol. 2005, 42: 134-141. 10.1603/0022-2585(2005)042[0134:DAAETO]2.0.CO;2.View ArticlePubMedGoogle Scholar
- Day JF, Shaman J: Using hydrologic conditions to forecast the risk of focal and epidemic arboviral transmission in peninsular Florida. J Med Entomol. 2008, 45: 458-465. 10.1603/0022-2585(2008)45[458:UHCTFT]2.0.CO;2.View ArticlePubMedGoogle Scholar
- Dohm DJ, O'Guinn ML, Turell MJ: Effect of environmental temperature on the ability of Culex pipiens (Diptera: Culicidae) to transmit West Nile virus. J Med Entomol. 2002, 39: 221-225. 10.1603/0022-2585-39.1.221.View ArticlePubMedGoogle Scholar
- Meyer RP, Hardy JL, Reisen WK: Diel Changes in Adult Mosquito Microhabitat Temperatures and Their Relationship to the Extrinsic Incubation of Arboviruses in Mosquitos in Kern County, California. J Med Entomol. 1990, 27: 607-614.View ArticlePubMedGoogle Scholar
- Reisen WK, Fang Y, Martinez VM: Effects of temperature onthe transmission of West Nile virus by Culex tarsalis (Diptera: Culicidae). J Med Entomol. 2006, 43: 309-317. 10.1603/0022-2585(2006)043[0309:EOTOTT]2.0.CO;2.View ArticlePubMedGoogle Scholar
- Kilpatrick AM, Meola MA, Moudy RM, Kramer LD: Temperature, viral genetics, and the transmission of West Nile virus by Culex pipiens mosquitoes. Plos Pathogens. 2008, 4 (6): e1000092-10.1371/journal.ppat.1000092.PubMed CentralView ArticlePubMedGoogle Scholar
- Monath TP, Tsai TF: St Louis encephalitis: lessons from the last decade. Am J Trop Med Hyg. 1987, 37: 405-595.Google Scholar
- Kunkel KE, Novak RJ, Lampman RL, Gu WD: Modeling the impact of variable climatic factors on the crossover of Culex restauns and Culex pipiens (Diptera: Culicidae), vectors of West Nile virus in Illinois. Am J Trop Med Hyg. 2006, 74: 168-173.PubMedGoogle Scholar
- Landesman WJ, Allan BF, Langerhans RB, Knight TM, Chase JM: Inter-annual associations between precipitation and human incidence of West Nile virus in the United States. Vector-Borne Zoonotic Dis. 2007, 7: 337-343. 10.1089/vbz.2006.0590.View ArticlePubMedGoogle Scholar
- Bertolotti L, Kitron UD, Walker ED, Ruiz MO, Brawn JD, Loss SR, Hamer GL, Goldberg TL: Fine-scale genetic variation and evolution of West Nile Virus in a transmission "hot spot" in suburban Chicago, USA. Virology. 2008, 374: 381-389. 10.1016/j.virol.2007.12.040.View ArticlePubMedGoogle Scholar
- Reisen WK: Effect of temperature on Culex tarsalis (Diptera, Culicidae) from Coachella and San-Joaquin valleys of California. J Med Entomol. 1995, 32: 636-645.View ArticlePubMedGoogle Scholar
- Kitron U: Landscape ecology and epidemiology of vector-borne diseases: Tools for spatial analysis. J Med Entomol. 1998, 35: 435-445.View ArticlePubMedGoogle Scholar
- Openshaw S: 'The modifiable areal unit problem'. Concepts and Techniques in Modern Geography. 1984, 38: 41-Google Scholar
- Biggerstaff BJ: PooledInfRate, Version 3.0: a Microsoft Excel Add-In to compute prevalence estimates from pooled samples. 2006, Ft. Collins, CO: Centers for Disease Control and PreventionGoogle Scholar
- Baker DG, Sharratt BS, Chiang HC, Zandlo JA, Ruschy DL: Base temperature selection for the prediction of European corn = borer instars by the growing degree-day method. Agricultural and Forest Meteorology. 1984, 32: 55-60. 10.1016/0168-1923(84)90028-5.View ArticleGoogle Scholar
- Lam NSN: Spatial interpolation methods - a reivew. American Cartographer. 1983, 10: 129-149. 10.1559/152304083783914958.View ArticleGoogle Scholar
- Wong DW, Yuan L, Perlin SA: Comparison of spatial interpolation methods for the estimation of air quality data. Journal of Exposure Analysis and Environmental Epidemiology. 2004, 14: 404-415. 10.1038/sj.jea.7500338.View ArticlePubMedGoogle Scholar
- Messina J, Brown WM, Ruiz MO: West Nile virus in the greater Chicago area, 2002-2006: changing patterns of human illness and social and environmental determinants of risk. Proceedings of the URISA GIS in Public Health Conference; New Orleans, LA. 2007, [http://www.urisa.org]Google Scholar
- Housing USCoPa: Summary Population and Housing Characteristics: Indiana. Book Summary Population and Housing Characteristics: Indiana (Editor ed.^eds.). 2001, City: Washington: Government Printing OfficeGoogle Scholar
- R Development Core Team: R: A language and environment for statistical computing. 2008, R Foundation for Statistical Computing, Vienna, Austria, [http://www.R-project.org]Google Scholar
- Gu WD, Lampman R, Novak RJ: Problems in estimating mosquito infection rates using minimum infection rate. Journal of Medical Entomology. 2003, 40: 595-596. 10.1603/0022-2585-40.5.595.View ArticlePubMedGoogle Scholar
- Olden JD, Lawler JJ, Poff NL: Machine learning methods without tears: A primer for ecologists. Quarterly Review of Biology. 2008, 83: 171-193. 10.1086/587826.View ArticlePubMedGoogle Scholar
- Hu WB, Tong SL, Mengersen K, Oldenburg B, Dale P: Mosquito species (Diptera: Culicidae) and the transmission of ross river virus in Brisbane, Australia. J Med Entomol. 2006, 43: 375-381. 10.1603/0022-2585(2006)043[0375:MSDCAT]2.0.CO;2.View ArticlePubMedGoogle Scholar
- Breiman L: Random forests. Machine Learning. 2001, 45: 5-32. 10.1023/A:1010933404324.View ArticleGoogle Scholar
- Strickman D: Rate of oviposition by Culex quinguefasciatus in San Antonio, Texas, during three years. J Am Mosq Control Assoc. 1988, 4: 339-344.PubMedGoogle Scholar
- Chaves LF, Pascual M: Comparing Models for Early Warning Systems of Neglected Tropical Diseases. Plos Neglected Tropical Diseases. 2007, 1: 10.1371/journal.pntd.0000033.Google Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.