Abstract:Zika virus disease is a mosquito-borne disease which is mainly transmitted by the female Aedes mosquitoes. Initially, since the 1950, the Zika virus disease (or Zika fever) was known to occur only in a narrow equatorial belt from Africa to Asia but from 2007 to 2016, the virus spread eastward (i.e. across the Pacific Ocean to the Americas) and that led to the 2015–2016 Zika virus epidemic. Now, it is creating a viral threat to societies all over the World because it spreads very fast in a globally connected world. The Zika virus is a member of the virus family Flaviviridae and its dynamics is governed by the day time activities of Aedes mosquitoes (i.e. A. aegypti and A. albopictus, etc.). However, there also exist new evidences of its transmission dynamics being governed via sexual contact and that has been reported from different nations (i.e. Argentina, Canada, Chile, France, Italy, New Zealand, Peru, Portugal, and USA). This paper introduces a mathematical model to investigate the transmission dynamics of Zika virus disease through the inclusion of the impact of media. Our proposed model involves both the possible routes of transmission (i.e. mosquito borne and sexual transmission). The existence and stability of different equilibria of this model are presented and analyzed in detail. The basic reproduction number R0\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$R_0$$\end{document} of the model is computed and it is found that for R0<1\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$R_0<1$$\end{document}, the disease free equilibrium of the model is locally asymptotically stable. When R0>1\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$R_0>1$$\end{document}, the endemic equilibrium point exists and is locally asymptotically stable under some restriction on parameters. For R0<1\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$R_0<1$$\end{document}, backward bifurcation occurs, which suggests that lowering R0\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$R_0$$\end{document} below one is not enough to eliminate the disease from the population. Finally, we extend our model to optimal control model by introducing two types of control parameters and compare the results of both the models using numerical simulation. Our presented results highlight the efficacy of the role of media and show that the proposed model is efficient in the analysis of transmission dynamics of the Zika virus disease. The sensitivity analysis of the parameters suggests that the key parameters are the mosquito biting rate and the rate of transmission due to the sexual contact between human to human. A small increase in any of these parameters can lead to great increase in the basic reproduction number R0\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$R_0$$\end{document} and this causes a large outbreak of this disease. So, the control policies should aim for controlling both the routes of transmission of this disease.

Mathematical modeling on co-infection: transmission dynamics of Zika virus and Dengue fever

Mathematical Modeling of the Co‐Infection Dynamics of Dengue and Malaria Using Delay Differential Equations

A Preliminary Mathematical Model for the Dynamic Transmission of Dengue, Chikungunya and Zika

Mathematical modelling of the first HIV/ZIKV co-infection cases in Colombia and Brazil

A new mathematical model for Zika virus transmission

Mathematical modeling of Zika virus with vertical transmission in the presence of Wolbachia-infected mosquitoes

A mathematical model of Zika virus transmission with saturated incidence and optimal control: A case study of 2016 zika outbreak in Puerto Rico

Mathematical models of SIR disease spread with combined non-sexual and sexual transmission routes

The Mathematical Modeling of the Host–Virus Interaction in Dengue Virus Infection: A Quantitative Study

An optimal control model for COVID‐19, zika, dengue, and chikungunya co‐dynamics with reinfection

The analysis of a new fractional model to the Zika virus infection with mutant

Analysis of a fractional-order model for dengue transmission dynamics with quarantine and vaccination measures

Zika Virus Model with the Caputo–Fabrizio Fractional Derivative

Fractional dynamics of a two-strain dengue model with co-infections and saturated incidences

Mathematical analysis of a three-strain model for dengue with vaccination

Beyond the biting - limited impact of explicit mosquito dynamics in dengue models

Some models for epidemics of vector-transmitted diseases

Modeling and optimal control analysis of Zika virus with media impact

A mathematical model to study the 2014-2015 large-scale dengue epidemics in Kaohsiung and Tainan cities in Taiwan, China

Mathematical Modelling to Investigate a Wolbachia Intervention to Reduce Dengue Transmission

A mathematical model for varicella-zoster and HIV co-dynamic supported by numerical simulations