Spatial modeling of two mosquito vectors of West Nile virus using integrated nested Laplace approximations

The abundance of Culex restuans and Culex pipiens in relation to ecological predictors is poorly understood in regions of the United States where their ranges overlap. It is suspected that these species play different roles in spreading West Nile virus (WNV) in these regions, but few studies have modeled these species separately or accounted for spatial correlation using Bayesian methods. We used mosquito surveillance data collected by the Pennsylvania Department of Environmental Protection from 2002 to 2016 and integrated nested Laplace approximations with the stochastic partial differential equation approach to predict C. restuans and C. pipiens abundance in relation to several ecological predictors. We then made a predictive risk surface of abundance for each species at locations that were not sampled. Explanatory variables in the models included ecological variables previously described to be important predictors of the abundance of these mosquito species. Developed habitat, temperature, and precipitation were important predictor variables for the abundance of C. restuans, whereas developed habitat, snow water equivalent, and normalized difference water index were important predictor variables for the abundance of C. pipiens. The abundance of C. restuans had a negative relationship with developed habitat in contrast to C. pipiens abundance, which had a positive relationship with developed habitat. Julian date was modeled as a temporal trend for both species and showed C. restuans to be more abundant from late April through late June and C. pipiens to be more abundant from July through September. A seasonal crossover was observed between these two species on Julian day 185, 4 July. We observed different spatial patterns of abundance in the predictive risk maps of each of the species. Our results indicate that modeling the abundance of these species spatially and separately in regions where these two mosquito vectors coexist can help gain further insight into understanding the epidemiology of WNV in human and susceptible animal populations.