Modelling the Potential Role of Media Campaigns in Ebola Transmission Dynamics

A six-compartment mathematical model is formulated to investigate the role of media campaigns in Ebola transmission dynamics. Themodel includes tweets or messages sent by individuals in different compartments.Themedia campaigns reproduction number is computed and used to discuss the stability of the disease states. The presence of a backward bifurcation as well as a forward bifurcation is shown together with the existence and local stability of the endemic equilibrium. Results show that messages sent through media have a more significant beneficial effect on the reduction of Ebola cases if they are more effective and spaced out.


Introduction
The world faced one of the most devastating Ebola virus disease (EVD) outbreaks ever in between 2014 and 2015.EVD is caused by a virus called Ebola, which was discovered in the Democratic Republic of Congo in 1976 near a river called Ebola [1].There are five known species of Ebola: Zaire ebolavirus which has caused the 2014 Ebola disease outbreak [2], Sudan ebolavirus, Cote d'Ivoire ebolavirus, Bundibugyo ebolavirus (Uganda), and Reston ebolavirus which has not yet caused disease in humans [3].This virus lives in animals like bats and primates, mostly found in Western and Central Africa.The virus can be transmitted from animals to humans when an individual comes into contact with an infectious animal through handling of contaminated meat, for example, and contamination is also possible among animals.Contamination can occur among humans when they have nonprotected contact with an infectious individual's fluids like faeces, vomit, saliva, sweat, and blood [4].It can also happen in hospitals, where healthcare practitioners paid a heavy price [1].
Symptoms can appear after 2 to 21 days following contamination and the infectious period can last from 4 to 10 days [5].Some contaminated individuals become symptomatic after 21 days [6], whereas others will never develop symptoms and remain asymptomatic [4,7,8].When the virus gets into a human body, it rapidly replicates and attacks the immune system.So, depending on the state of the infected individual's immune system, death can directly follow or recovery can occur after treatment.According to the World Health Organisation (WHO), a suspected case of EVD is any person, alive or dead, suffering or having suffered from a sudden onset of high fever and having had contact with a suspected or confirmed Ebola case, a dead or sick animal, and at least three of the following symptoms: headaches, anorexia, lethargy, aching muscles or joints, breathing difficulties, vomiting, diarrhoea, stomach pain, inexplicable bleeding, or any sudden inexplicable death [9].Confirmed cases of EVD are individuals who would have tested positive for the virus antigen either by detection of virus RNA by Reverse Transcriptase Polymerase Chain Reaction or by detection of IgM antibodies directed against Ebola [9].Ebola seropositive individuals can be either asymptomatic or symptomatic.Post-Ebola survey results showed that 71% of seropositive individuals monitored were asymptomatic [7].Symptomless EVD patients have low infectivity due to their very low viral load whereas the symptomatic cases transmit the disease through their fluids [8].
Media campaigns have been included in mathematical models in recent years.Exponential functions are mostly used to represent their impact on people's behaviour which affects disease evolution [10,11].
A model where media coverage influences the transmission rate of a given disease is presented in [12].An 2 International Journal of Differential Equations exponentially decreasing function is used to describe the media coverage over time.The results show that media coverage has a short-term beneficial effect on the targeted population.A smoking cessation model with media campaign was given in [13] and results showed that the reproduction number is suppressed when media campaigns that focus on smoking cessation are increased.Thus, spreading information to encourage smokers to quit smoking was an effective intervention.The impact of Twitter on influenza was studied in [14].An exponential term was associated to model the effect of Twitter messages on reducing the transmission rate of influenza.It was noted that Twitter can have a substantial influence on the dynamics of influenza virus infection and can provide a good real-time assessment of the current disease condition.
There is no large-scale treatment for EVD as yet, so stopping the transmission chain remains the only viable form of control.Media campaigns publicise the means of contracting the disease and the behaviour to adopt when a suspected or confirmed Ebola case is detected.The potential effect of media campaigns on Ebola transmission dynamics is thus of great interest.This paper is motivated by the work in [14] and was done as an M.S. research work by the first author [15].We use a mathematical model to describe the transmission dynamics of EVD in the presence of asymptomatic cases and the impact of media campaigns on the disease transmission is represented by a linear decreasing function.The efficacy of media campaigns is a state variable in this model and a differential equation describing its variation is given.We examine the long-term dynamics of EVD and evaluate the potential impact of media campaigns on reducing the number of Ebola cases.The paper is arranged as follows: the model formulation is presented in Section 2, and the model properties and analysis are given in Section 3. The numerical simulations are presented in Section 4 and we give concluding remarks in the last section.

Model Formulation
A deterministic model with six independent compartments comprising individuals that are susceptible (), exposed (), infected asymptomatic (  ), infected symptomatic (  ), recovered (), and deceased () is formulated.The total population size  is given by We only consider the Zaire Ebola virus strain which caused the 2014 Ebola outbreak in West Africa.We assume a constant natural death rate  for the whole model.The study is made over a relatively large period so that those who recover from EVD gain permanent immunity against the strain.Recruitment into the susceptibles class is done through birth or migration at a constant rate Λ and susceptible individuals become exposed after unsafe contact with Ebola virus.After contamination, susceptibles move to compartment  and, considering 1/ as the incubation period, individuals leave the exposed compartment at a rate .After  The general objective of media campaigns against a disease is to increase the population's awareness of the disease and correct misperceptions about how it is spread and how it is and is not acquired [18].The efficacy of messages sent through media is thus their ability to produce the intended results.We consider here that Ebola disease related messages are exchanged by individuals from each of the compartments at any time .After receiving tweets or messages related to Ebola disease, the population decides on the means of preventing or even treating the disease.Messages are assumed to get outdated at a rate .() is defined as the fraction of effective messages sent by individuals of the respective classes at any time .Thus, () is the ratio of effective messages to the total messages sent.The contributions to  from the living compartments are, respectively,  1 ,  2 ,  3 ,  4 , and  5 .The use of the campaigns is to reduce EVD transmission.We assume here that media campaigns primarily target the transmission process and 0 < () ≤ 1, ∀ ≥ 0.
The force of infection will be given by where  is the probability that a contact will result in an infection and  is the number of contacts between susceptible and infectious individuals.The parameter  > 1 describes the high infectivity of dead bodies.The flow diagram is presented in Figure 1.

Model Properties and Analysis
3.1.Existence and Uniqueness of Solutions.The right hand side of system (3)-( 9) is made of Lipschitz continuous functions since they describe the size of a population.According to Picard's Existence Theorem, with given initial conditions, the solutions of our system exist and they are unique.
Proof.We first assume that  >  during the modelling time.
Together with  which is already bounded (see proof in Appendix A), we can say that Ω is bounded and at limiting equilibrium lim →∞ () = Λ/.Besides, any sum or difference of variables in Ω with positive initial values will remain in Ω or in a neighbourhood of Ω.Thus, Ω is positively invariant and attracting with respect to the flow of system (3)-(9).

Positivity of Solutions
Theorem 2. The existing solutions of system ( 3)-( 9) are all positive.
Proof.From (3), we can have Solving for (13) yields which is positive given that (0) is also positive.Similarly, from (4), we have so that which thus shows that () is positive since (0) is also positive.Similarly, from (5), we can write from which we obtain Thus,   is positive since   (0) is positive.
The remaining equations yield So,   (), (), and () are all positive for positive initial conditions.Thus, all the state variables are positive.
International Journal of Differential Equations

Steady States Analysis
. This model has two steady states: the disease-free equilibrium (DFE) which describes the total absence of EVD in the studied population and the endemic equilibrium (EE) which exists at any positive prevalence of EVD in the population.This section is dedicated to the study of local and global stability of these steady states.

3.4.
The Disease-Free Equilibrium and   .The diseasefree equilibrium is given by ( * ,  * ,  *  ,  *  ,  * ,  * ,  * ) = (Λ/, 0, 0, 0, 0, 0, Λ 1 /).To compute the media campaigns reproduction number   , we use the next generation method comprehensively discussed in [19].The renewal matrix  and transfer matrix  at DFE are where The media campaigns reproduction number   is the spectral radius of the matrix  −1 and is given by We can rewrite   =  1 +  2 for elucidation purposes where and ] = ( +  1 )/ 2 .Note here that / 1 is the probability that an individual in  moves either to   or to   ./ 3 is the proportion of symptomatic individuals who die from EVD.Thus, the media campaigns reproduction number is a sum of secondary infections due to infectious individuals in   and the deceased in .Notice here the reduction factor 1− * which represents the attenuating effect of media campaigns on the future number of EVD cases.Theorem 3. The DFE is globally asymptotically stable whenever   <    < 1, where    = min((()), (, ])) and (, ]) will be defined later.When    <   < 1, the DFE is locally stable.Otherwise, the DFE is unstable.
Thus,  is a positive definite function at the DFE.
The derivative of  is given by Also, / ≤ 1 and at equilibrium Plugging ( 25) into (24) yields ( ( ()) − 1)   (26) with Thus, V ≤ 0 when (()) ≤ 1 and, particularly, V = 0 only if  =   =   =  = 0. Since () ≥  * for all  > 0 (see proof in Appendix A), we have   < (()).Because the largest invariant set for which V = 0 in Ω is the DFE and V ≤ 0 if (()) ≤ 1, by using the invariance principle of LaSalle [20], we can conclude that the DFE is globally asymptotically stable for   < (()) < 1. Together with the existence of a backward bifurcation later proven, we finally obtain the global stability of the DFE for   <    < 1.  the reproduction number can be found and then used to analyse the disease evolution.For instance, when   = 1, the critical value of media campaigns efficacy   can be determined.Since the behaviour of the system changes when the reproduction number crosses the value one,   can also be used as a threshold parameter that indicates a behavioural change of the system and thus can help in the disease control for any given set of parameter values.

Existence and Stability of the Endemic Equilibrium.
In this section, we show the existence of the endemic equilibrium (EE).We denote the endemic equilibrium by ( * * ,  * * ,  * *  ,  * *  ,  * * ,  * * ,  * * ).At equilibrium, (3)-( 9) give where we obtain the following equation: where with From (31),  * * = 0 corresponds to the DFE discussed in the previous section.The signs of the solutions of the quadratic equation are given in Table 1.
From Table 1, we notice that, for the existence and uniqueness of the endemic equilibrium, ] 0 must be negative.This is only possible if   > 1.Thus, we have the following theorem on the existence of the endemic equilibrium.3)-( 9) has a unique endemic equilibrium.
(iv) If 0 <   <    , then system (3)-( 9) does not admit any endemic equilibrium and only the DFE exists.Provided ] 1 < 0, the existence of two endemic equilibria for   < 1 suggests the existence of a backward bifurcation since the DFE also exists in that particular domain.The coexistence of DFE and endemic equilibrium when   < 1 is a well known characteristic of a backward bifurcation described in [21].Thus, there exists a critical value of   , denoted by    , for which there is a change in the qualitative behaviour of our model.
At the bifurcation point, there is an intersection between the line   =    and the graph of ( * * ).The discriminant Δ is equal to zero at   = (, ]), which is solution of Equation (35) implies Considering as well the threshold value of the reproduction number from Theorem 3, we can conclude that    = min((()), (, ])).So, 0 <   <    , the DFE is globally stable, the DFE is locally stable and two endemic equilibria exist with one which is stable and the other one unstable. (37) The DFE and EE both describe different qualitative behaviours of our epidemic.Let us set  = (1 −  * ) as our bifurcation parameter, so that In order to describe the stability of the endemic equilibrium, we use the theorem, remark, and corollary in [22] which are based on the Centre Manifold Theory, and formulated in Appendix B.
Theorem 5. A unique endemic equilibrium exists when   > 1 and is locally asymptotically stable.
Proof.For model ( 3)-( 9), the DFE ( 0 ) is not equal to zero.According to Remark 1 in [22], we notice that if the equilibrium of interest in Theorem B.1 is a nonnegative equilibrium  0 , then the requirement that  is nonnegative in Theorem B.1 is not necessary.When some components in  are negative, one can still apply Theorem B.1 on condition that  () > 0, if  0 () = 0, if  0 () > 0,  () does not need to be positive, where () and  0 () denote the th component of  and  0 , respectively.Firstly, let us rewrite system (3)-( 9) introducing The equilibrium of interest here is the DFE denoted by  0 = ( * , 0, 0, 0, 0, 0,  * ) and the bifurcation parameter is  * .The linearisation matrix  of our model at The eigenvalues of  are − (twice), −, 0, and the roots of polynomial (42) below: where Our linearisation matrix  will thus have zero as simple eigenvalue.Statement (A1) is verified.We now show that (A2) is satisfied.

(44)
Setting we have We notice that Besides, since  0 ( 1 ) and  0 ( 7 ) are positive,  1 and  7 do not need to be positive according to Remark 1 in [22].So, statement (A2) is verified.
The formulas of the constants  and  are  After multiple derivations, we have Since  < 0 and  > 0, by using the fourth item of Theorem B.1, we can conclude that when  * changes from negative to positive,  0 changes its stability from stable to unstable.Correspondingly, a negative unstable equilibrium becomes positive and locally asymptotically stable and a forward bifurcation appears [21].

Bifurcation Analysis.
The study of the DFE and EE led to the proof of the existence of a backward and a forward bifurcation for our model.Graphically, they are, respectively, represented in Figures 3(a) and 3(b) where   is chosen as the bifurcation parameter.We have shown that a forward bifurcation exists for values of   greater than one.This means that EVD will persist as long as secondary infections will occur and reducing   to values less than one is enough to eradicate EVD.However, the existence of a backward bifurcation makes it difficult to control the epidemic.In fact, the coexistence of the DFE and the EE for   in [   , 1] shows that reducing the number of secondary infections to less than one is not enough to eradicate EVD.Other control measures like quarantine and contact tracing should be implemented together with media campaigns to reach a globally stable DFE and wipe out EVD.
Figure 4 shows time series plots for the force of infection  for varying initial conditions.The trajectories converge to steady states depending on the initial conditions and the values of   .We can observe that when the DFE is asymptotically stable (  <    ), the force of infection reaches zero.When the DFE is locally stable and the EE is unstable (   <   < 1), the force of infection reflects a persistent infection.When the EE is unstable (  > 1), the force of infection is maximal.This confirms the results obtained at the bifurcation analysis and describes the unstable nature of EVD which can easily become an explosive epidemic after a small increase in its force of infection as   passes through 1.

Numerical Simulations
In this section, we use Matlab to carry out simulations for our model.We first verify our theoretic conclusions related to stability analysis of system (3)-( 9) and then we vary our parameters values to better understand how media campaigns influence the prevalence and transmission of EVD.
It is important to note that the figures chosen are for illustrative purposes only, as we endeavour to verify the analytic results.

Parameters' Estimation.
The parameters used in the simulations are either obtained from the literature or estimated.Since the mean infectious period is set to be from 4 to 10 days, the highest recovery rate  2 is set to 1/4.The recovery rate of asymptomatic individuals is assumed to be greater than the one of the symptomatic individuals since the former have stronger resistance to EVD.Without any reliable source for EVD media related data, we assume that individuals can send EVD related messages through media independently of their disease status.At the beginning of the epidemic, there is neither a recovered nor an asymptomatic infected individual since only symptomatic persons transmit the disease.We also assume that messages are transmitted through media at time  = 0, at least for preventive purpose.The setting of the initial conditions is driven by the fact that the population of Nzérékoré, the region where this 2014 Ebola disease outbreak started in Guinea, is estimated to be 1,663,582 individuals [23].We consider the introduction of infectives in the population and high infectivity of dead bodies ( = 1.5).The initial conditions are then Table 2 gives the description of parameters and their values.

Sensitivity Analysis.
In mathematical modelling, parameters whose values are not precisely known are often used and may vary within some ranges.Numerical methods used to solve equations derived from models may introduce numerical errors in the results.The effects of such errors or uncertainties in the model's parameters are quantified through sensitivity analysis.The aim of sensitivity analysis is to quantify the influence of parameters variation on calculated results [24].
Sensitivity indices allow us to measure the relative change in a state variable when a parameter changes.The normalized forward sensitivity index of a variable to a parameter is the ratio of the relative change in the variable to the relative change in the parameter.When the variable is a differentiable function of the parameter, the sensitivity index may be alternatively defined using partial derivatives (see [25]).Definition 6.The normalized forward sensitivity index of a variable, , that depends differentiably on a parameter, , is defined as Media campaigns in this paper contribute to the limitation of the disease transmission.The reproduction number   is an important concept when it comes to the disease transmission, because it helps to determine EVD incidence.The normalized forward sensitivity indices of   with respect to each parameter  in expression (22) are given by  those with the highest absolute values.Negative and positive correlations of the parameters to the reproduction number are indicated by negative and positive signs.The parameters Λ and  1 have the largest absolute negative numerical values with negative sensitivity index values.Thus, increasing their values will decrease EVD incidence.This result can be explained by the fact that any increase in the two parameters leads to an increase in the efficacy of media campaigns.So, the more the media campaigns are efficacious, the less the value of the reproduction number is.Another important parameter with a negative index is , which is an expected result since burials of EVD dead bodies limit the disease transmission due to infected corpses.The outdating rate of media campaigns  has the most positive influence on EVD reproduction number.This means that the more frequently the messages spread by media on EVD are updated, the lower the number of new infections is.The larger the values of  and  are, the less the efficacy of the messages is and the more the disease spreads.The parameters  and  form the transmission rate and their increase will directly contribute to an increase in the number of EVD cases.

Simulations Results and Interpretation.
Figures 5(a) and 5(b) confirm the results on stability analysis.It follows that when   < 1, the epidemic dies out and, for   > 1, EVD becomes endemic.This is a graphical description of the fact that the DFE is locally stable for   < 1 and the EE is locally stable whenever   > 1.
The media campaigns reproduction number   is made of parameters which differently influence its values in a variety of ways.The relationship between those parameters can be evaluated through contour plots.We chose two parameters,  1 and , whose influence on the reproduction number is clearly significant as shown in the expression of   .Figure 6 shows that  1 largely influences   when compared to .Increasing the values of  1 decreases   .Thus, in order to eradicate EVD, the exchange of EVD related messages is critical in eradicating the epidemic.

Discussion and Conclusion
To model the potential effect of media campaigns on Ebola transmission, we used a deterministic model, with compartments comprising individuals with different EVD infection status, who send EVD related messages through media.The effect of media campaigns on people's behaviour is represented by a reduction factor which decreases the number of new EVD cases.Stability analysis was presented in terms of the model reproduction number   .It was shown that the disease-free and the endemic equilibria are locally stable if   < 1 and   > 1, respectively.The inclusion of the asymptomatic infected class resulted in the model exhibiting a backward bifurcation, emphasizing the necessity of intense efforts against EVD as a result of undetected asymptomatic cases.The existence of a backward bifurcation has important implications in the design of policies and strategies to eradicate or control an epidemic.In the presence of a backward bifurcation, classical policies on disease eradication need to be changed as EVD can persist even when the threshold parameter   is less than one.
To be able to control EVD, governments and international stakeholders should implement feasible campaigns taking into account the social, economic, and mainly the cultural realities of the affected countries.Interventions from these campaigns should target the affected populations and help them to understand the disease, comply with control measures, which sometimes seem severe, and change their behaviour in order to stop the disease transmission chain [1].The best way to contain this outbreak is to jointly implement case isolation, contact tracing with quarantine, and sanitary funeral practices as suggested in [26].
This model is not without shortcomings.People's reaction to media campaigns does not always lead to a reduction in the number of future contamination cases.So, a function representing the influence of media campaigns on individuals' behaviour which takes into account the different cultural settings would be an innovative and informative addition to this model.Aspects of quarantine, contact tracing, and case identification initiatives are possible additions that can make this model more reliable.Despite these shortcomings, this model provides a good description of EVD outbreak.The model investigates a very important aspect in disease control in our times, that is, the use of social media in spreading messages.

6 International
Journal of Differential Equations Theorem 4. (i) If   > 1, (34) has a unique positive solution and system (

Figure 4 :
Figure 4: Time series variation of the force of infection for different values of the reproduction number.
Reproduction Number   .  is considered as a reproduction number whose values depend on the fraction of effective messages on EVD at a given time.Assuming  to be constant, Figure2graphically describes the relationship between two concepts: reproduction number and media campaigns efficacy.It shows the reducing effect of media campaigns on the number of EVD infected individuals and indicates as well how we can test the efficacy of Ebola related messages through the pace of the disease transmission.In fact, for each value of , the corresponding value of

Table 2 :
Parameter values and their description.

Table 3
represents the numerical values of the sensitivity indices of the reproduction number   for the parameters used in the model.The most important parameters are International Journal of Differential Equations

Table 3 :
Sensitivity indices for EVD reproduction number.