Abstract

An epidemic model with time delay has been proposed and analyzed. In this model the effect of awareness programs driven by media on the prevalence of an infectious disease is studied. It is assumed that pathogens are transmitted via direct contact between the susceptible and the infective populations and further assumed that the growth rate of cumulative density of awareness programs increases at a rate proportional to the infective population. The model is analyzed by using stability theory of differential equations and numerical simulations. Both equilibria have been proved to be globally asymptotically stable. The results we obtained and numerical simulations suggest the increasing of the dissemination rate and implementation rate can reduce the proportion of the infective population.

1. Introduction

Infectious diseases that cause mortality, disability, and social and economic disruption are a major threat to mankind, which are responsible for a quarter of all deaths annually in the world [13]. Once an infectious disease appears and spreads in a region, the Center of Disease Control and Prevention will do its best to stop the propagation of the disease [4]. One of the measures is to tell people the appropriate preventive knowledge of diseases as soon as possible through media and education which make people take precautions to reduce their chances of being infected [5]. As the awareness disseminates, people will change their behaviors to alter their susceptibility. There is much evidence that media coverage can play an important role in the spread and control of infectious diseases [611]. In [8] Cui and others established a framework of transmission coefficient with media coverage which is a decreasing function of the number of the infective individuals and they observed a classic threshold-type behavior, with the disease becoming extinct when and going to a globally asymptotically stable equilibrium when , and they concluded that media coverage is critical in disease eradication.

The models studying the spread of infectious diseases are very useful in evaluating strategies to control the diseases in populations. Recently some authors studied the impacts of media coverage and education on the spread of infectious diseases in a given region [1214]. In [12] Misra et al. proposed a nonlinear mathematical model for the effects of awareness programs on the spread of infectious diseases and assumed the growth rate of awareness programs impacting the population is proportional to the number of infective individuals. The model analysis showed that the spread of an infectious disease can be controlled by using awareness programs but the disease remains endemic due to immigration. Yorke and London [15] proposed an SIS type compartmental model for sexually transmitted infections with the assumption that the whole population is aware of risk but only a certain proportion choose to respond by limiting their contacts with the infective population. As a result the spread of infection is controlled, leading to a reduction of the number of individuals becoming infected. A nonlinear mathematical model with delay to capture the dynamics of effect of awareness programs on the prevalence of any epidemic is proposed and analyzed [16], which assumed it increases at a rate proportional to the number of the deaths of infective individuals.

Some scholars focus on the contact rate and most of them assume that media campaigning will aid in modifying the contact rate between susceptible and infective individuals [1727]. To prevent the unboundedness of the contact rate, Capasso [17] used a saturated incidence rate of the form ; Liu and coworkers [23, 24] used a nonlinear incidence rate given by to incorporate the effect of behavioral changes. In [25], the authors study how media coverage influences the dynamics of infectious disease by using SIRS model with the contact rate , where , and are positive constants. Tchuenche and Bauch used an exponentially decreasing function to reveal the force of media [27]. Cui et al. used a similar function as [13] and developed an model using incidence rate to investigate the impact of media coverage on the transmission [28]. Stability analysis of the model has shown that Hopf bifurcation can occur which causes oscillatory phenomena when is sufficiently small. Numerical simulations suggested that the media impact was stronger when the basic reproduction number . Liu et al. [29] described the impact of media coverage using the transmission coefficient , where is the number of hospitalized individuals. And this impact leads to the change of avoidance and contact patterns at both individual and community levels. Liu and Xiao [30] introduce a segmented function to describe the saturated media impact when formulating an epidemic model. A Filippov epidemic model with media coverage is proposed to describe the real characteristics of media impact in the spread of infectious diseases by incorporating a piecewise continuous transmission rate in [31]. Mathematical and bifurcation analysis with regard to the local, global stability of equilibria and local sliding bifurcations are performed. Bhunu et al. [32] and Tchuenche and Bauch[27] focused on the different types of population in their work.

In order to better describe population mixed conditions, some authors study infectious diseases models on the different networks [33, 34]. Funk et al. have overlaid the model of information spread of a contagious disease on two, not necessarily identical networks, with more informed individuals acting to reduce their susceptibility [35]. Liu et al. took into the random perturbation [36]. In [37], the authors extended the classical SIRS epidemic model from a deterministic framework to a stochastic differential equation, and then they gave the conditions of existence of unique positive solution and the stochastic extinction and discussed the exponential p-stability and global stability.

Most of the articles, such as [12], assume that, due to awareness programs, driven by media, some susceptible individuals will avoid their contacts with the infectious individuals resulting in formation of a new class in the population and this newly formed aware class may contract infection only if they lose awareness. But we regard it is unreasonable. In fact sometimes even if persons are conscious of diseases, they will also be infected. Therefore we propose a mathematical model for predicting the future course of any epidemic by considering this newly formed aware class into the model.

The rest of this paper is organized as follows. In the next section, to capture the dynamics of effect of cumulative density of awareness programs on the prevalence of any epidemic a mathematical model is proposed and analyzed. Then we analyze the local and global stability of the disease-free and the unique endemic equilibrium in Section 3. Furthermore, in Section 4 we perform some numerical examples to validate the analytical findings in Section 3 and introduce the importance of the dissemination rate and implementation rate in disease control. In Section 5 we discuss the above content.

2. Mathematical Model and Equilibrium Analysis

In this paper due to the awareness programs about the disease, susceptible individuals are rarely in contact with the infective ones and form a different class, namely, aware susceptible class; thus the total population is divided into three classes, the susceptible population, the aware population, and the infective population, and the proportions of them in the total population are , , and . Assuming that at time the cumulative density of awareness programs driven by media in that region is , which increases at a rate proportional to the infective population and consumes due to causes like inefficiency and psychological barriers at , thus

In fact people cannot take measures to protect themselves in time after the media reports the disease, so we introduce a time delay that represents the interval between the report time and the time of taking measures. We assume that a proportion of infected individuals recover through treatment and, after recovery, a fraction of recovered people will become aware and join the aware susceptible class whereas the remaining fraction will become unaware susceptible. Keeping the above facts in mind, the dynamics of a model are governed by the following systems of nonlinear ordinary differential equations: Here , ,  and   .

Assume diseases spread due to the direct contact between susceptible and infective individuals only. In the above model, the rate of immigration of susceptible population is . is the contact rate of aware susceptible with infective population and is the contact rate of unaware ones. The constant represents the dissemination rate of awareness among unaware susceptible due to which they form a different class; then denotes the rate of transfer of aware susceptible to unaware class. The constants , represent the recovery rate and the natural death rate, respectively. All the above constants are assumed to be positive. Using the fact that , the system (2) becomes as follows: Now it is sufficient to study system (3) in detail rather than system (2).

For the analysis of system (3), we need the region of attraction which is given by the set: , which attracts all solutions initiating in the interior of the positive orthant.

Define the basic reproduction number . There are two equilibria, the disease-free equilibrium and the endemic equilibrium ; the existence of is trivial; then we prove the existence of in detail. In the values of , are obtained by solving the following algebraic equations (for ): Using (5) and (6), we get Further, using (7) and (8) in (4), we obtain a quadratic equation in as where Solving (9) we get We obtain and when and get for .

Remark. From the expression of , it is easy to note that and , which shows that the equilibrium number of infective individuals decreases as the rate of dissemination and the implementation rate of awareness programs increase.

3. Stability Analysis

In this section we present the local and global stability of and .

3.1. The Stability of the Disease-Free Equilibrium

Theorem 1. When , the disease-free equilibrium is locally asymptotically stable if and becomes unstable if .

Proof. The Jacobian matrix corresponding to system (4) is given as follows: where is the eigenvalue. Then the characteristic equation is We get ,  and   . So when ; when .

Theorem 2. When , the disease-free equilibrium is globally asymptotically stable in if .

Proof. To establish the global stability of the disease-free equilibrium , we use Lyapunov’s method and consider the following positive definite function without : Now differentiating with respect to , we get When , . The largest compact invariant set in when is the singleton . Then LaSalle’s invariance principle implies that is globally asymptotically stable in .

3.2. The Stability of the Endemic Equilibrium

Linearizing system (4) about , let , , and and get where We have the disease-free equilibrium and the endemic equilibrium , where the stability of about system (3) is corresponding with . The characteristic equation of the above system at is where is the eigenvalue and The stability of the endemic equilibrium of system (3) is stated in the following theorems.

Theorem 3. When , the endemic equilibrium is locally asymptotically stable if .

Proof. When , the characteristic equation is of the form It is easy to see that From , we can get ; thus . According to Hurwitz criterion, we can know all the ’s have negative real parts; then is locally asymptotically stable.
When , notice that (18) does not have nonnegative real roots. If it has roots with nonnegative real parts they must be complex and should have been obtained from a pair of complex conjugate roots which cross the imaginary axis. Consequently, (18) must have a pair of purely imaginary solutions for some . Assume that ( ) is a root of (18) without loss of generality. That is the case if and only if satisfies the equation Separating the real and imaginary parts, we have the following system, satisfied by : To eliminate the trigonometric functions we square both sides of each equation above and we add the squared above equations to obtain the following forth order equation in : To reduce this fourth order equation in to a quadratic equation let and denote the coefficients by where So all the coefficients of (25) are positive numbers. Then according to Lemma in [38], (25) has no positive real roots; that is, we may not get any positive value of , which satisfy the transcendental equation (18). So all the ’s have negative real parts for all values of the delay ; then is locally asymptotically stable. Thus when , is locally asymptotically stable if .

Theorem 4. When , the endemic equilibrium is globally asymptotically stable in if .

Proof. Using Lyapunov’s method, we consider the following positive function: The derivative of along the system is given by The largest compact invariant set when is the singleton . Then LaSalle’s invariance principle implies that is globally asymptotically stable; that is, is globally asymptotically stable in .

4. Numerical Simulations and Results

To check the feasibility of our analysis of , we present some numerical computations in this section using Matlab by choosing the following set of parameter values: , , , , , , , , , ,  and   when . Let ; it may be checked that the condition of existence of the endemic equilibrium . The equilibrium values for this data are obtained as The basic reproduction number , for the above set of parameter values, is found to be .

For the above parameter values, we select five sets of different initial starts; then the computer generated graphs of aware population, infective population, and cumulative density of awareness programs, respectively, have been drawn in Figure 1, which shows that all the trajectories initiating inside the region of attraction approach towards and , respectively. Both of the equilibria and are locally asymptotically stable for given set of parameter values that numerical simulations support the analysis given in Section 3. In fact, they are globally asymptotically stable in as we have proved.

In the following, we research the relationships of , and the dissemination rate , the implementation rate separately. We make , change from to and get the performances of and , the trajectories of which with respect to time for different and are shown in Figures 2 and 3, respectively. And there are no awareness programs when and are equal to zero. As shown in Figures 2 and 3, both increase and both reduce as the increase of and , which proves the conclusions of the remark. And , are greater influence on than , which state awareness programs have a positive effect on prevention of diseases. In addition the reason why ( ) has a similar trend as the variations of and is that and have a similar influence on ( ). From the figure we obtain that can postpone the time of the balance of equilibrium; thus we can have more time to formulate measures to prevent diseases. There really is an effort here to make it clear that and (awareness programs) play a key role in the prevention and control of diseases.

5. Discussion

The media is widely acknowledged as a key tool for influencing people’s behaviors towards the disease to devise proper policies for controlling the epidemic. Awareness programs through media make people be aware about the disease and take various precautions to reduce their chances of being infected. In this paper, we propose and analyze a mathematical model to study the effect of awareness programs driven by media and the delay on the prevalence of an infectious disease. It is assumed that pathogens are transmitted via direct contact between the susceptible and the infective populations. Assume further that cumulative density of awareness programs increases at a rate proportional to the infective population. The model exhibits two equilibria; the disease-free equilibrium has been shown to be stable for basic reproduction number . For , it becomes unstable, which leads to the existence of an endemic equilibrium. The endemic equilibriums are globally asymptotically stable. The delay has no effect on the stability of the system. The numerical simulations and results that prove the stability of equilibria suggest that if we want to reduce the proportions of the infective population and increase the aware population, we can increase the dissemination rate and implementation rate . They are conducive to controlling the spread of diseases.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This work is supported partially by the National Natural Science Foundation of China (no. 10901145) and the Natural Science Foundation of Shanxi (2012011002-1).