A mathematical model for Zika virus transmission dynamics with a time-dependent mosquito biting rate

Mathematical modeling has become a tool used to address many emerging diseases. One of the most basic and popular modeling frameworks is the compartmental model. Unfortunately, most of the available compartmental models developed for Zika virus (ZIKV) transmission were designed to describe and reconstruct only past, short-time ZIKV outbreaks in which the effects of seasonal change to entomological parameters can be ignored. To make an accurate long-term prediction of ZIKV transmission, the inclusion of seasonal effects into an epidemic model is unavoidable.

We developed a vector-borne compartmental model to analyze the spread of the ZIKV during the 2015-2016 outbreaks in Bahia, Brazil and to investigate the impact of two vector control strategies, namely, reducing mosquito biting rates and reducing mosquito population size. The model considered the influences of seasonal change on the ZIKV transmission dynamics via the time-varying mosquito biting rate. The model was also validated by comparing the model prediction with reported data that were not used to calibrate the model.

We found that the model can give a very good fit between the simulation results and the reported Zika cases in Bahia (R-square = 0.9989). At the end of 2016, the total number of ZIKV infected people was predicted to be 1.2087 million. The model also predicted that there would not be a large outbreak from May 2016 to December 2016 due to the decrease of the susceptible pool. Implementing disease mitigation by reducing the mosquito biting rates was found to be more effective than reducing the mosquito population size. Finally, the correlation between the time series of estimated mosquito biting rates and the average temperature was also suggested.

The proposed ZIKV transmission model together with the estimated weekly biting rates can reconstruct the past long-time multi-peak ZIKV outbreaks in Bahia.