#### Abstract

We develop a nonlinear mathematical model with the effect of awareness programs on the binge drinking. Due to the fact that awareness programs are capable of inducing behavioral changes in nondrinkers, we introduce a separate class by avoiding contacts with the heavy drinkers. 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 stability of the alcohol free and the alcohol present equilibria and give some numerical simulations to explain our main result. Our results show that awareness programs is an effective measure in reducing alcohol problems.

#### 1. 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–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.

#### 2. The Model

##### 2.1. System Description

The total population is divided into five compartments, namely, the susceptible compartment who do not drink or drink only moderately, denoted by , those who are aware of risk and avoid drinking, denoted by , 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 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. 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. 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
**
with the initial conditions , , , , and are positively invariant for system (2).*

*Proof. *Adding the first four equations of (2), we have
It follows that
where is the initial total number of people. Thus,
Then ; 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 . Further, from the last equation of system (2), we have
It follows that
where represents the initial value of cumulative density of awareness programs. Thus,
It implies that the region
is a positively invariant set for (2). So we consider dynamics of system (2) on the set in this paper.

###### 2.2.2. Positivity of Solutions

For system (2), to ensure that the solutions of the system with positive initial conditions remain positive for all , it is necessary to prove that all the state variables are nonnegative, so we have the following lemma.

Lemma 2. *If , , , , and , then the solutions , , , , and of system (2) are positive for all .*

*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 such that
In that case, from the first equation of system (2), we have
which is a contradiction meaning that , .

Or there exists a such that
In that case, from the second equation of system (2), we have
which is a contradiction meaning that , .

Or there exists a such that
In that case, from the third equation of system (2), we have
which is a contradiction meaning that , . Similarly, it can be shown that and for all . Thus, the solutions , , , , and of system (2) remain positive for all .

#### 3. Analysis of the Model

There are one alcohol free equilibrium and one alcohol present equilibrium for system (2).

##### 3.1. Alcohol Free Equilibrium and the Basic Reproduction Number

The model has an alcohol free equilibrium given by 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 and at the alcohol free equilibrium are, respectively, where The model reproduction number denoted by is thus given by

##### 3.2. Global Stability of

Theorem 3. *For system (2), the alcohol free equilibrium is globally asymptotically stable if .*

*Proof. *We introduce the following Lyapunov function [26]:
The derivative of is given by

If , then . As we know, and , so we obtain . Furthermore, only if or . The maximum invariant set in is the singleton . By Lasalle’s invariance principle [27], is globally asymptotically stable in .

##### 3.3. Alcohol Present Equilibrium

###### 3.3.1. Existence of the Alcohol Present Equilibrium

Theorem 4. *If , system (2) has a unique alcohol present equilibrium , where
*

*Proof. *It follows from system (2) that
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 , 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 when ; there is no positive root in the interval when .

For , we carry out the simulation (see Figure 4) and give the following conjecture.*Conjecture*. If , the alcohol present equilibrium of (2) is globally asymptotically stable.*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.

#### 4. 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 and ; in addition, we assume that , , and .

First, we choose and the numerical simulation gives ; then the alcohol free equilibrium is globally asymptotically stable (Figure 2). Second, we choose and the numerical simulation gives ; the alcohol free equilibrium is globally asymptotically stable (Figure 3). Third, we choose and the numerical simulation gives ; the alcohol present equilibrium is globally asymptotically stable (Figure 4).

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] is governed by the following system of nonlinear ordinary differential equations:

We choose and the numerical simulation gives ; 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 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.

#### 5. 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) to the following model with delay: We leave these works for the future.

#### Conflict of Interests

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

#### Acknowledgments

This work was partially supported by the NNSF of China (10961018), the NSF of Gansu Province of China (1107RJZA088), the NSF for Distinguished Young Scholars of Gansu Province of China (1111RJDA003), the Special Fund for the Basic Requirements in the Research of University of Gansu Province of China, and the Development Program for Hong Liu Distinguished Young Scholars in Lanzhou University of Technology.