Modelling the Influence of Awareness Programs by Media on the Drinking Dynamics

and Applied Analysis 3 R relapsing into drinkers. m represents the growth rate of density of awareness programs which is assumed to be proportional to the number of heavy drinkers, and γ represents the depletion rate of these programs due to ineffectiveness, social problems, and so forth. For instance, as time passes by, some of the campaigns by media lose their influence or fail to attract people. 2.2. Basic Properties 2.2.1. Invariant Region Lemma 1. All feasible regions Ω defined by Ω = { (S (t) , X (t) , A (t) , R (t) ,M (t)) ∈ R 5 + : S (t) + X (t) + A (t) + R (t) ≤ 1; M (t) ≤ m γ } (3) with the initial conditions S(0) ≥ 0,X(0) ≥ 0,A(0) ≥ 0,R(0) ≥ 0, andM(0) ≥ 0 are positively invariant for system (2). Proof. Adding the first four equations of (2), we have dN dt = μ − μN. (4) It follows that N(t) = 1 + N (0) e −μt , (5) whereN(0) is the initial total number of people. Thus, lim t→∞ N(t) = 1. (6) Then S(t) + X(t) + A(t) + R(t) ≤ 1; the subpopulation is represented in proportion to the total number of population. Thus, since the system (2) monitors human population, it is plausible to assume that all its state variables and parameters are nonnegative for all t ≥ 0. Further, from the last equation of system (2), we have


Introduction
Alcohol consumption, particularly heavier drinking, is an important risk factor for many health problems and is a major contributor to the global burden of disease [1].It has been identified as an important risk factor for chronic disease and injury [2].Drinking alcohol causes the liver to become even more damaged in people infected with the hepatitis B virus (HBV) [3].The CAS (Central Authentication Service) findings have shown that alcohol consumption at binge levels and beyond has a significant impact on college students' academic performance, social relationships, risk taking behaviors, and health.This form of drinking is associated with missing class, falling behind in schoolwork, and lower grade point average, a relationship mediated by fewer hours spent in studying [4,5].In addition, students who attended schools with high rates of binge drinking experienced a greater number of secondhand effects, including disruption of sleep or study, property damage, and verbal, physical, or sexual violence, than their peers attending schools with low binge drinking rates [6].Prior studies have indicated that heavy alcohol drinkers are likely to engage in risky sexual behaviours and, thus, more likely to get sexually transmitted infections (STIs) than social drinkers [7].The rapid consumption of large amounts of alcohol, especially by young people, leads to serious antisocial and criminal behavior in urban centres.This phenomenon has grown very rapidly [8].Since mathematical model is a predictive tool which can mimic the process of drinking and provide useful tools to analyze the spread and control of drinking behavior, several different mathematical models for drinking have been formulated and studied recently.Simple models for alcohol treatment are presented by Sanchez et al. [9] which were based on studying binge drinking in a college system and assumed the same "leaving rate." Mulone and Straughan [10] investigated a possible model for binge drinking taking into account admitting and nonadmitting drinkers.But the global analysis of binge drinking model is not discussed in the literature.Huo and Song [11] introduced a more realistic two-stage model for binge drinking problem, where the youths with alcohol problems are divided into those who admit the problem and those who do not admit it.For the other mathematical models for drinking or smoking, we refer to [12][13][14] and the references cited therein.
The media may be the most important source of health information for the general public; it can play a special role in providing a voice for people to express their experiences of illness [15].Jaramillo [16] studied a model presenting an evaluation of the impact on case finding of a mass media health education campaign for TB control in Cali, Colombia.Liu et al. [17] incorporate the psychological impacts of media coverage on disease transmission and emphasize the importance of refining classical mathematical models to reflect the partially self-limiting nature of infectious disease outbreaks, due to the effects of widespread news coverage and fast information flow.Web-based screening and brief interventions that include personalized feedback about their alcohol use have proven to be particularly promising for reducing hazardous drinking among university students [18].Mass media (television, radio, newspapers, billboards, and booklets) have been used as a way of delivering preventive health messages as they have the potential to influence people's behavior.Further, mass media have been deployed in the effort to control and eliminate epidemic diseases [19].Xiang et al. [20] studied a drinking model with public health educational campaigns which can be an efficient option for reducing the spread of disease.Misra et al. [21,22] considered the effects of awareness programs driven by media on the spread of infectious diseases; their results showed that, though awareness programs cannot eradicate infection, they help in controlling the prevalence of disease.The most successful approaches concentrate not on trying to change the drinker's behavior directly but on influencing the community environment around the nondrinker.For instance, education and persuasion may increase knowledge and change attitudes but have no effect on drinking [23,24], so we only consider the interaction from awareness programs between aware individual and nondrinker.
Motivated by the above works, we introduce a mathematical model with the effect of awareness programs on the binge drinking in this paper.Because awareness programs are capable of inducing behavioral changes in susceptible individuals, we introduce a separate class () of those who are aware of risk and refrain from drinking by avoiding contact with the heavy drinkers.This newly formed aware class is assumed to be fully protected from the heavy drinkers due to media-induced isolation and individuals of this class may contact heavy drinkers only if they lose awareness [21,22].Furthermore, we assume that cumulative density of awareness programs () increases at a rate proportional to the number of heavy drinkers.We establish some sufficient conditions for the global stability of equilibria and give some numerical simulations to explain our main result.Our results show that awareness programs are an effective measure in reducing alcohol problems.
The organization of this paper is as follows.In the next section, the effect of awareness programs by media on binge drinking model is formulated.In Section 3, the existence and the stability of equilibria are investigated.Some numerical simulations are given in Section 4. Some discussions are given in Section 5. (), those who drink heavily, denoted by (), and those who are in treatment, denoted by ().Also let () be the cumulative density of awareness programs driven by media.The growth rate of density of awareness programs is assumed to be proportional to the number of heavy drinkers.It is considered that, due to the awareness programs, nondrinkers form a different class and avoid contacting with the heavy drinkers.And we only consider the interaction from awareness programs between aware individual and nondrinker.The total number of population at time  is given by

The Model
The model structure is shown in Figure 1.The transfer diagram leads to the following system of ordinary differential equations: where  is the rate of individuals entering into the system in a given time interval (say, each year).Since we are dealing with youths, we assume that the death rate is negligible and so the leaving rate is also .This is in accordance with an assumption by Sanchez et al. [9], where "leaving rate" presents the departure rate from drinking environment. is the transmission coefficient for the nondrinkers turning to heavy drinkers through peer pressure with the initial conditions (0) ≥ 0, (0) ≥ 0, (0) ≥ 0, (0) ≥ 0, and (0) ≥ 0 are positively invariant for system (2).
Proof.Adding the first four equations of (2), we have It follows that where (0) is the initial total number of people.Thus, lim Then () + () + () + () ≤ 1; the subpopulation is represented in proportion to the total number of population.Thus, since the system (2) monitors human population, it is plausible to assume that all its state variables and parameters are nonnegative for all  ≥ 0. Further, from the last equation of system (2), we have It follows that where (0) represents the initial value of cumulative density of awareness programs.Thus, lim is a positively invariant set for (2).So we consider dynamics of system (2) on the set Ω in this paper.

Positivity of Solutions.
For system (2), to ensure that the solutions of the system with positive initial conditions remain positive for all  > 0, it is necessary to prove that all the state variables are nonnegative, so we have the following lemma.
Proof.Under the given initial conditions, it is easy to prove that the solutions of the system (2) are positive; if not, we assume a contradiction that there exists a first time  1 such that In that case, from the first equation of system (2), we have which is a contradiction meaning that () > 0,  > 0.
Or there exists a  2 such that In that case, from the second equation of system (2), we have which is a contradiction meaning that () > 0,  ≥ 0.
Or there exists a  3 such that In that case, from the third equation of system (2), we have which is a contradiction meaning that () > 0,  ≥ 0.

Analysis of the Model
There are one alcohol free equilibrium  0 and one alcohol present equilibrium  * for system (2).

Alcohol Free Equilibrium and the Basic Reproduction
Number.The model has an alcohol free equilibrium given by  0 = (1, 0, 0, 0, 0) .
In the following, the basic reproduction number of system (2) will be obtained by the next generation matrix method formulated in [25].
Let  = (, , , , )  ; then system (2) can be written as where ) . ( The Jacobian matrices of F() and V() at the alcohol free equilibrium  0 are, respectively, where The model reproduction number denoted by  0 is thus given by

Existence of the Alcohol Present Equilibrium
From the fourth and fifth equations of ( 27), we obtain Substituting  into the third equation of ( 27), we have Then substituting ( 29) and ( 30) into the second equation of ( 27), we get For  ̸ = 0, substituting ( 29), (30), and (31) into the first equation of (27) gives Substituting ( 30) into (32), we get Therefore, there exists a unique positive root in the interval (0, 1) when  0 > 1; there is no positive root in the interval [0, 1] when  0 ≤ 1.
For  0 > 1, we carry out the simulation (see Figure 4) and give the following conjecture.
Remark.Since the global stability of ( 2) is a hard problem, we only carry out simulation (see Figure 4).How to prove the global stability of the alcohol present equilibrium  * of (2) and give the conditions based on model parameters is still an open problem.

Numerical Simulation
To illustrate the analytic results obtained above, we give some simulations using the parameter values in Table 1.Numerical results are displayed in the following figures.It is reported that binge drinking occurs in up to 30 percent of adolescents and alcohol use disorders occur in about 6 percent of this age group; approximately 50 to 60 percent of them remain abstinent at the end of a year's treatment and a majority of those stay dry permanently [29,30].So the initial condition is taken as the lowest by (0) = 0.3 and (0) = 0.15; in addition, we assume that (0) = 0.1, (0) = 0.45, and (0) = 0.002.
Finally, for showing that awareness programs by media can help to control the population of binge drinking, we compare our model with the model which has been studied by Mulone and Straughan [10]; the model in [10] We choose  = 0.032 and the numerical simulation gives  0 = 1.1645 > 1; the alcohol present equilibrium  * of (34) is globally asymptotically stable.The number of the binge drinking population is generally on the rise without effects of awareness programs by media.
For showing the distinction between ( 2) and (34), we take  = 0.032 and compare the variations of the binge drinking population ().We can see that, because of awareness programs by media (Figure 5), the number of people who have alcohol problems decreases due to awareness-induced isolation of susceptible people, so making more people avoid contacting with the heavy drinkers due to the fact that the awareness programs by media can reduce the alcohol problems.
Further, the variations of the binge drinking population () with respect to time  for different values of growth rate of density of awareness programs  are shown in Figure 6.From these figures, it is apparent that as the rate of density of awareness programs  increases the binge drinking population () decreases.Then, the awareness programs by media are an effective measure in alcohol problems.

Discussion
We have formulated the effect of awareness programs on the binge drinking and analyzed their dynamical behaviors.The model exhibits two equilibria, explicitly the alcohol free equilibrium and alcohol present equilibrium.By constructing Lyapunov function, we establish some sufficient conditions for the global stability of the alcohol free and the alcohol present equilibria and give some numerical simulations to explain our main result.In our model, we have assumed that the cumulative density of awareness programs by media is proportional to the population of binge drinking heavily.It may be noted that the number of alcohol consumption cases known to the policy makers is sometimes old and thus the number of awareness depends on data.So it is more plausible to consider a time delay in the growth rate of awareness programs [21].We can modify (2) We leave these works for the future.
is governed by the following system of nonlinear ordinary differential equations: ( + ) .

Figure 6 :
Figure 6: Variations of binge drinking population with time for different values of .
represents the growth rate of density of awareness programs which is assumed to be proportional to the number of heavy drinkers, and  represents the depletion rate of these programs due to ineffectiveness, social problems, and so forth.For instance, as time passes by, some of the campaigns by media lose their influence or fail to attract people.
.  represents the dissemination rate of awareness among nondrinkersdue to the awareness programs, nondrinkers form a different class () and avoid contacting with the heavy drinkers.The constant  denotes the rate of transfer of aware individuals to nondrinker class. is a proportionality coefficient for  going into treatment. is a proportionality coefficient for  relapsing into drinkers.

Table 1 :
Description and estimation of parameters.