Fusing an agent-based model of mosquito population dynamics with a statistical reconstruction of spatio-temporal abundance patterns

The mosquito Aedes aegypti is the vector of a number of medically-important viruses, including dengue virus, yellow fever virus, chikungunya virus, and Zika virus, and as such vector control is a key approach to managing the diseases they cause. Understanding the impact of vector control on these diseases is aided by first understanding its impact on Ae. aegypti population dynamics. A number of detail-rich models have been developed to couple the dynamics of the immature and adult stages of Ae. aegypti. The numerous assumptions of these models enable them to realistically characterize impacts of mosquito control, but they also constrain the ability of such models to reproduce empirical patterns that do not conform to the models' behavior. In contrast, statistical models afford sufficient flexibility to extract nuanced signals from noisy data, yet they have limited ability to make predictions about impacts of mosquito control on disease caused by pathogens that the mosquitoes transmit without extensive data on mosquitoes and disease.Here, we demonstrate how the differing strengths of mechanistic realism and statistical flexibility can be fused into a single model. Our analysis utilizes data from 176,352 household-level Ae. aegypti aspirator collections conducted during 1999-2011 in Iquitos, Peru. The key step in our approach is to calibrate a single parameter of the model to spatio-temporal abundance patterns predicted by a generalized additive model (GAM). In effect, this calibrated parameter absorbs residual variation in the abundance time-series not captured by other features of the mechanistic model. We then used this calibrated parameter and the literature-derived parameters in the agent-based model to explore Ae. aegypti population dynamics and the impact of insecticide spraying to kill adult mosquitoes.The baseline abundance predicted by the agent-based model closely matched that predicted by the GAM. Following spraying, the agent-based model predicted that mosquito abundance rebounds within about two months, commensurate with recent experimental data from Iquitos. Our approach was able to accurately reproduce abundance patterns in Iquitos and produce a realistic response to adulticide spraying, while retaining sufficient flexibility to be applied across a range of settings. Author summaryThe mosquito Aedes aegypti is the vector for a number of the most medically important viruses, including dengue, Zika, chikungunya, and yellow fever. Understanding the population dynamics of this mosquito, and how those dynamics might respond to vector control interventions, is critical to inform the deployment of such interventions. One of the best ways to gain this understanding is through modeling of population dynamics. Such models are often categorized as either statistical or dynamical, and each of these approaches has advantages and disadvantages-for instance, statistical models may more closely match patterns observed in empirical data, while dynamical models are better able to predict the impact of counterfactual situations such as vector control strategies. In this paper, we present an approach which fuses these two approaches in order to gain the advantages of both: it fits empirical data on Aedes aegypti population dynamics well, while producing realistic responses to vector control interventions. Our approach has the potential to inform and improve the deployment of vector control interventions, and, when used in concert with an epidemiological model, to help reduce the burden of the diseases spread by such vectors.