key: cord-0999012-yjruu6wt authors: Fatima, BiBi; Alqudah, Manar A.; Zaman, Gul; Jarad, Fahd; Abdeljawad, Thabet title: Modeling the Transmission Dynamics of Middle Eastern Respiratory Syndrome Coronavirus with the Impact of Media Coverage date: 2021-03-24 journal: Results Phys DOI: 10.1016/j.rinp.2021.104053 sha: 5c109fc76c6369d22b278e9a3b9f47f6c078093b doc_id: 999012 cord_uid: yjruu6wt Middle East respiratory syndrome coronavirus has been persistent in the Middle East region since 2012. In this paper, we propose a deterministic mathematical model to investigate the effect of media coverage on the transmission and control of Middle Eastern respiratory syndrome coronavirus disease. In order to do this we develop model formulation. Basic reproduction number [Formula: see text] will be calculated from the model to assess the transmissibility of the (MERS-CoV). We discuss the existence of backward bifurcation for some range of parameters. We also show stability of the model to figure out the stability condition and impact of media coverage. We show a special case of the model for which the endemic equilibrium is globally asymptotically stable. Finally all the theoretical results will be verified with the help of numerical simulation for easy understanding. Infectious diseases are responsible for a quarter of all death in the world annually, such as SARS, MERS, and now COVID-19, that exhibit some distinct features such as rapid spread and visible symptoms [1] [2] [3] . One of the initiatives is to inform individuals through media and education as quickly as possible the right preventive understanding of the disease. We understand from experience that the more preventive the inhabitants understand, the better they can stop the disease from spreading. In recent years, there have been global fears over viruses. Outbreaks of MERS, SARS, or influenza A are deadly and spread fast. It is one area where the media can play a life-saving role. Timely, accurate information, under the umbrella of risk communication or disaster communication, can help curb outbreaks and drive people to seek treatment early [4, 5] . Mass media reports can induce individual behaviour change during a disease outbreak, which has been found to be useful as it reduces the force of infection. In May 2015, South Korea was hit by a massive, deadly outbreak of the MERS virus [6] [7] [8] . It was a little known virus that started with the patient suffering symptoms of a common cold, but could kill within two weeks. The virus had first spread in Saudi Arabia, where 40 percent of those who contracted it died , [9] [10] [11] . In South Korea 186 people were infected and 38 of those died within two months of falling ill. A lack of information about the nature of the illness, and about where the outbreaks were occurring, led to a widespread panic. The research focused on media coverage of the 2015 Middle East respiratory syndrome (MERS) crisis in South Korea. In particular by the country's three major terrestrial television stations: KBS and MBC, both public broadcasters, and SBS, a commercial channel. Coverage of the MERS outbreak by these three television stations between May and July 2015 was examined. The author also interviewed eleven journalists and editors, via email and the messaging app 'Kakao talk' voice function, between 7th April to 7th June 2017. Other interviews were conducted in the UK, via events at the Refuters Institute for the Study of Journalism in Oxford. Mathematical modeling and analysis are used for the dynamics of infectious diseases, see for instance [12] [13] [14] [15] [16] [17] [18] [19] . There have been mathematical modeling studies to analyze the impact of media coverage on the spread and control of infectious disease in a given population. In mathematical epidemiology the role of media communication in alerting the outcome of infectious disease outbreak, continuously having a place. The paper aimed to analyse the flow of information during an epidemic and to understand the impact of media coverage on the transmission of infectious and hospitalized individuals. In [20] , the authors extend the classical SEI model by considering a new incidence functional which reflects the impact of the media coverage to the spreading and control of the disease. The incidence function has been considered to play a key role in ensuring that the model indeed give reasonable qualitative description of the transmission dynamics of the diseases. We consider the model of [20] , by taking the hospitalize class. The paper is arrange as follow: In Section 2, we discuss formulation of the model, disease free equilibria and reproductive number. In Section 3, we discuss endemic equilibria and the existence of backward bifurcation for some range of parameter. In Section 4, we discuss sensitivity analysis of the model. In Section 5, we find local and global stability analysis. In Section 6, we discuss a special case of the model and obtained global asymptotic stability of the proposed model. In Section 7, we discuss numerical simulation of the proposed model. In Section 8, we give discussion on the obtained results. We examine the transmission of MERS CoV in a specified region. We distribute the population into the following compartments. The susceptible S(t). • All new born will goes to the susceptible class only, • The infected I(t) who are infectious, • The hospitalize H(t), • The recovered population R(t). To incorporate the effect of the behavioral changes of the susceptible individuals, we used a non-linear incidence rate. The mathematical model based on SIHR model with the incident of mass action is given by: with In model (1) • π is the recruitment rate of susceptible population. • µ 0 represent natural death rate and κ is death rate occur due to disease. • β 1 , β 3 are the contact rate before media alert. • β 1 − β 2 I m+I , and β 3 − β 4 H m+H are the contact rate after media alert. The model (1) have a disease free equilibrium denoted by E 0 and given by E 0 = (S 0 , 0, 0, 0), where the components are define as; S 0 = π µ 0 . Basic reproductive number a threshold representing how many secondary infections results from the introduction of one infected individuals in a susceptible population. For the basic reproductive number R 0 , we use the method of Driessche and Watmough [21] . then R 0 is the spectral radius of FV −1 , that is The basic reproduction number R 0 of this model consists of two parts, representing the two different transmission routes; i.e., from the infected individuals, from hospitalized individuals before media alert. Where β 1 , β 3 show contact rate before media alert. Proof: The total population is represented by Suppose the left hand side of each differential equation (1) be zero, the endemic (S, I, H, R) satisfies putting the above expression into the first equation of model (1) and after simplification, we have The significance of backward bifurcation in the epidemiological model is that of the classical requirement of the basic reproduction number R 0 to be less than one [22, 23] , while necessary for the elimination of the MERS CoV virus from population. The presence of backward bifurcation in the proposed model suggests that the feasibility of MERS virus elimination, when the basic reproduction number is less than one, depends on the initial size of the sub population of the model. For R 0 = 1, the following result holds. We bring out sensitivity analysis of parameters using in the proposed model. This analysis will make it easy to know the parameters that have a essentially effect on reproductive number. We apply the technic given in [24, 25] and given by, To examine the local and global stability analysis of the model (1) about E 0 = (S 0 , 0, 0, 0), we use the following results. Parameter Sensitivity indices Parameter Sensitivity indices Theorem 2. We take the model (1) with all positive parameters. For R 0 > 1 the model (1) possess a unique endemic equilibrium E * (S * , I * , H * , R * ) and is locally asymptotically stable. For R 0 < 1 the model (1) get a unique disease free equilibrium E 0 (S 0 , I 0 , H 0 , R 0 ) and is globally asymptotically stable. Proof: The Jacobian matrix of the suggested model (1) about the DFE point E 0 is The first two eigen values have already negative real part λ 1 = −µ 0 < 0, λ 2 = −(µ 0 + γ) < 0 for the rest of eigen value we take 2 × 2 matrix, by Routh-Hurwitz criteria [26] , we have to prove that trace of A is negative and determinant of A is positive, if R 0 < 1, Hence det(A) > 0 if R 0 < 1 and β 3 ≤ 0, which implies that det(A) is positive, if R 0 < 1. Therefore, Trac(A) < 0 and det(A) > 0 if and only if R 0 < 1. Thus the disease free equilibrium is locally asymptotically stable at E 0 Let us consider Lypunov function V = I. Differentiating V with the solution of model (1), we obtain Hence E 0 is globally stable at disease free equilibrium point. For R 0 > 1, the jacobian matrix at the equilibrium E * is where where a 1 > 0,a 2 > 0, a 3 > 0, a 4 > 0, also a 1 a 2 a 3 > a 2 3 + a 2 1 a 4 . It follows from Routh Hurwtiz criteria all the eigen values (11) have negative real part if R 0 > 1, which means that E * is locally asymptotically stable. Suppose γ=0 in model (1) we have the following SIHR model The first three equations are independent of the fourth equation in the model (12) . We consider the reduced model as: The model (13) bear disease free equilibrium at E 0 = ( π µ 0 ). For endemic equilibrium point we put right hand side of (13) zero , Theorem 3. We take the model (13) with all positive parameters. For R 0 > 1 the model (13) possess a unique endemic equilibrium E * (S * , I * , H * , R * ) and is locally asymptotically stable. For R 0 < 1 the model (13) get a unique disease free equilibrium E 0 (S 0 , I 0 , H 0 , R 0 ) and is globally asymptotically stable. Proof: Jacobian matrix of the suggested model (13) about the point E 0 is When R 0 < 1, then (q + κ + µ 0 )(R 0 − 1) < 0, complete the proof. Jacobian matrix at the point E * = (S * , I * , H * ) is give by, where The characteristic equation of the above jacobian matrix is Where a 1 = (A + 2µ 0 + γ + q + κ), a 2 = (π + µ 0 + α + q)B + (A(α + µ 0 ) + DM, a 3 = π(q + κ + µ 0 )A + M(α + µ 0 )(R 0 − 1)C, a 1 a 2 = Aπ + γπ + κµ 0 + A(α + µ 0 ) + qκ + (R 0 − 1) + 2Dµ 0 , a 1 > 0, a 2 > 0, a 3 > 0, also a 1 a 2 > a 3 . From Routh Hurwtiz criteria [26] all the eigen values (16) have negative real part if and only if R 0 > 1, which means that E * is locally asymptotically stable. For global stability at DFE, we use Lyapunov function theory [27] . Theorem 4. For R 0 < 1 the disease free equilibrium of model (13) is stable globally , if S = S 0 other wise unstable if R 0 > 1. Proof. We define the following Lyapunov function is given by Where's w i i = 1, 2, 3, 4 are positive constant taking time derivative of (17) , we have By using dF dt is negative if S > S 0 and R 0 < 1 and dF dt = 0 if and only if S = S 0 . By Lasala inverience principle [28, 29] , the disease free equilibrium is globally asymptotically stable. For global stability at endemic equilibrium we used the geometrical approach [30, 31] . Proof: The linearized matrix and second additive compound matrix is denoted by J and J |2| model (13) Consider the function G(χ) = G(S, I, H) = diag S I , S I , S I , then, G −1 (χ) = diag I S , I S , I S , taking derivative of, G f (χ), we get , which can be written as where Let (n 1 , n 2 , n 3 ) be a vector in R 3 its norm . defined as n 1 , n 2 , n 3 = max{ n 1 , n 2 + n 3 }. Now by Martin et al [31] , ℓ(N) ≤ sup{g 1 , g 2 } = sup{ℓ(N 11 ) + N 12 ) , ℓ(N 22 ) + N 21 }, where g i = ℓ(N ii ) + N ij ) for i = 1, 2 and i = j, which implies that H m+H . Therefore g 1 and g 2 becomes, such that, g 1 ≤Ṡ S − 2µ 0 − (q + κ) and ℓ(N)dt < 0. Thus the system (1) around (S * , I * , H * ). is globally asymptotically stable. In this section, we solved the proposed deterministic model by using Runge-Kutta method of order 4th, see for detail [32] . To further understand the dynamical behavior of the proposed model we used numerical simulation to verify our analytical findings. In order to do this, we assumed some value of parameters, and some are taken from publish data given in Table. 1. The choice of numerical values of the parameter are taken in such a way that would be more biologically feasible. We also assume the time interval is 10 days with initial population for susceptible S(t), infected I(t), hospitalised H(t), recovered from MERS-CoV R(t). Moreover, the biological interpretation of these results states that if the value of basic reproductive number is less than one, then the susceptible population decreases , while then become stable and shows that there will be always stable susceptible population, see Fig.3a . The dynamics of I(t), H(t), R(t) reveals that the number of these populations will be decreases and reaches to zero as shown in Fig.3a-Fig.3d , which ensure the stability of the proposed model. One of the important factor is to find the relative impact of the basic reproductive number to epidemic parameters as shown in Fig.3 . We developed a mathematical model to analyze the impact of media coverage to the spread of infectious diseases in a given region. We get the following results from SIHRS, and SIHR model. We calculate basic reproductive number R 0 by the method of next generation method. When β 2 , β 4 = 0, the media coverage does not effect the reproductive number R 0 . We discuss the stability of the proposed model. Stability analysis show that the disease free equilibrium is locally asymptotically stable if R 0 < 1. If R 0 > 1, it is shown that a unique endemic equilibrium appears and bifurcation can occur which cause oscillatory phenomena. We discuss the role of media coverage on the spreading of MERS-CoV. Though the media coverage itself is not a determined fact, to eradicate the infection of the diseases, the analysis of the model indicates that, to certain extent, the more media coverage in a given population, the less number of individuals will be infected. Our analytical results show that the susceptible S(t), infected I(t), hospitalize H(t), recovered R(t) converge to equilibrium point. Which ensure the stability of the proposed model. In future, we are planning to develop an optimal mechanism on the basis of local dynamics and sensitivity analysis. This control strategy will help that how to eradicate the infection from the community. Work on such issue is in progress and will be reported soon in the form of a new article. Evidence for camel-to-human transmission of MERS coronavirus World Health Organization. Novel Coronavirus-China, Disease outbreak news : Update Modeling and analysis of effects of awareness programs by media on the spread of infectious diseases Understanding, compliance and psychological impact of the SARS quarantine experience Korea Centers for Disease Control and Prevention. Middle East respiratory syndrome coronavirus outbreak in the Republic of Korea Mary's Hospital. Epidemiologic features of the first MERS outbreak in Korea: focus on Pyeongtaek St The Characteristics of Middle Eastern Respiratory Syndrome Coronavirus Transmission Dynamics in South Korea Clinical course and outcomes of critically ill patients with Middle East respiratory syndrome coronavirus infection Evidence for camel-to-human transmission of MERS coronavirus Middle East respiratory syndrome coronavirus: a case-control study of hospitalized patients Fractional order mathematical modeling of COVID-19 transmission Stability analysis of a dynamical model of tuberculosis with incomplete treatment Computational and theoretical modeling of the transmission dynamics of novel COVID-19 under Mittag-Leffler power law Stability analysis and numerical solutions of fractional order HIV/AIDS model Haar wavelet collocation approach for the solution of fractional order COVID-19 model using Caputo derivative A fractional order HIV-TB coinfection model with nonsingular Mittag-Leffler Law On fractional order model of tumor dynamics with drug interventions under nonlocal fractional derivative. Results Phys On a comprehensive model of the novel coronavirus (COVID-19) under Mittag-Leffler derivative The impact of media coverage on the dynamics of infectious disease Reproduction numbers and sub threshold endemic equilibria for compartmental models of disease transmission Backward bifurcation and optimal control of a vector borne disease Backward bifurcation and control in transmission dynamics of arboviral diseases Modeling and control of the hepatitis B virus spreading using an epidemic model The transmission dynamic of different hepatitis B-infected individuals with the effect of hospitalization Routh-hurwitz criterion and bifurcation method for stability analysis of tuberculosis transmission model Stability analysis and optimal vaccination of an SIR epidemic model The stability of dynamical system Stability of nonautonomous system Dynamical models of tuberculosis and their applications Logarithmic norms and projections applied to linear differential systems The transmission dynamic and optimal control of acute and chronic hepatitis B This research was funded by the Deanship of Scientific Research at Princess Nourah bint Abdulrahman University through the Fast-track Research Funding Program.