MODELING ZIKA TRANSMISSION DYNAMICS: PREVENTION AND CONTROL

Abstract

The Zika virus (ZIKV) epidemic is depicted to have high spatial diversity and slow growth, attributable to the dynamics of the mosquito vector and mobility of the human populations. In an effort to understand the transmission dynamics of Zika virus, we formulate a new compartmental epidemic model with a system of seven differential equations and 11 parameters incorporating the decaying transmission rate and study the impact of protection measure on basic public health. We do not fit the model to the observed pattern of spread, rather we use parameter values estimated in the past and examine the extent to which the designed model prediction agrees with the pattern of spread seen in Brazil, via reaction–diffusion modeling. Our work includes estimation of key epidemiological parameters such as basic reproduction number ([Formula: see text], and gives a rough estimate of how many individuals can be typically infected during an outbreak if it occurs in India. We used partial rank correlation coefficient method for global sensitivity analysis to identify the most influential model parameters. Using optimal control theory and Pontryagin’s maximum principle, a control model has been proposed and conditions for the optimal control are determined for the deterministic model of Zika virus. The control functions for the strategies (i) vector-to-human contact reduction and (ii) vector elimination are introduced into the system. Numerical simulations are also performed. This work aimed at understanding the potential extent and timing of the ZIKV epidemic can be used as a template for the analysis of future mosquito-borne epidemics.