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 modelling of within-host Chikungunya virus dynamics with adaptive immune response

Dynamics of a periodic Chikungunya model with temperature and rainfall effects

Mathematical Modeling and Analysis of the Dynamics of Chikungunya in Bangladesh

Mathematical analysis of an HIV model with latent reservoir, delayed CTL immune response and immune impairment

Sensitivity Analysis in an Immuno-Epidemiological Vector-Host Model

Mathematical Modeling and Analysis of the Dynamics of RNA Viruses in Presence of Immunity and Treatment: A Case Study of SARS-CoV-2

Modeling the dynamical behaviour of the interaction of T-cells and human immunodeficiency virus with saturated incidence

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

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

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

Modeling and optimal control analysis of Zika virus with media impact

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

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

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

Stability and Optimality Criteria for an SVIR Epidemic Model with Numerical Simulation

An Immuno-Epidemiological Vector-Host Model with Within-Vector Viral Kinetics

Stability and Hopf Bifurcation of a Cytokine-Enhanced HIV Infection Model with Saturation Incidence and Three Delays

Continuous and discrete dynamics of a deterministic model of HIV infection

Modeling HIV-HBV Co-infection Dynamics: Stochastic Differential Equations and Matlab Simulation with Euler-Maruyama Numerical Method

In Vivo HIV Dynamics, Modeling the Interaction of HIV and Immune System Via Non-Integer Derivatives

Global co-dynamics of viral infections with saturated incidence