# A Kinetic Model for the Spread of Rumor in Emergencies

Authors: Guanghua Chen; Huizhang Shen; Teng Ye; Guangming Chen; Naphtali Kerr
Journal: Discrete Dynamics in Nature and Society
 (2013)

---

## Abstract

To study rumor spreading, a novel rumor spreading model with latent, constant recruitment, and varying total population is investigated. We carry out a global qualitative analysis of the rumor spreading model with latent, constant recruitment, and varying total population and studied the existence and stability of the rumor-elimination and prevailing equilibria and the impact different management strategies have on the rumor spreading process. The Runge-Kutta method is used for numerical simulation, and the simulation results imply that increasing the efforts to crack down on rumor spreaders, the timely disclosure of transparent information, the popularization of scientific knowledge, and the improvement of the ignorants’ discrimination may be effective rumor control measures.


---

## Body

## 1. Introduction

Rumor is a typical social phenomenon that runs through the whole evolutionary history of mankind, and its spreading plays a significant role in a variety of human affairs [1–4]. A rumor is “an unverified account or explanation of events circulating from person to person and pertaining to an object, event, or issue in public concern” [2]. Rumor is usually created to defame someone or to spread false information about public events; it not only infringes upon others’ interests, but also poses a threat to social stability, and it appears when a group tries to make sense of an ambiguous, uncertain, or chaotic situation. For example, in 2011, the Fukushima Daiichi nuclear disaster resulted from a 9.0-magnitude earthquake and ensuing tsunami in Japan caused a salt-buying frenzy in China. This salt-buying frenzy came from the rumors that Japan’s nuclear leak would contaminate the sea off east China’s Shandong province, thereby compromising the future safety of salt taken from the sea, and additionally it is believed that iodized salt helps protect people from nuclear radiation or the salt produced from contaminated seawater will be unsafe to eat. With the rumors spreading, this frenzy swept across China and caused social panic as well as instability in just a few days. From the afternoon of March 16, 2011, people in coastal areas such as Zhejiang, Jiangsu, and Shanghai, rushed out to buy up supplies, pushing up the price tenfold, before the panic spread inland. The next day, other inland areas such as Shanxi, Hubei, Anhui, Jiangxi, Guizhou, and Sichuan provinces were also caught up in the panic buying. The situation eased after the government ordered the securing of salt supplies and as experts confirmed that salt was almost useless to counter radiation. Other rumors have triggered similar responses elsewhere. Drugs stores and health food shops in Russia’s Far East and British Columbia, Canada, have reported shortages of iodine pills, despite health officials insisting that potassium iodide is not antiradiation drug. The “salt” farce not only affected market order and people’s lives, but also became an “international joke” that damaged the national image. The spread of rumors often results in financial loss, social instability, and, in extreme cases, even injuries and loss of lives [5]. Therefore, to study the spread of rumors under emergencies is an urgent issue for crisis management, and the cracking down on rumors is of great importance.Rumors can be viewed as an “infection of the mind”, and their spreading shows an interesting resemblance to that of epidemics. Pioneering contributions to their modeling, based on epidemic models, date back to [1, 5–20]. As far as we know, Daley first proposed the similarity between epidemics and rumors using mathematical analysis [6]. A classical rumor spreading model, introduced by Daley and Kendall [7], is conceptually similar to the susceptible-infective-recovered (SIR) model for epidemiology. Maki and Thompson [8] focused on the analysis of the rumor spreading model based on mathematical theory and developed the MT model. Many researchers used the DK model and the MT model extensively for their quantitative studies [10–16]. Kawachi [15] presented deterministic models for the transmission of a rumor. He investigated the age-independent case and introduced four models, which are classified according to whether the population is closed or not and whether the rumor is constant or variable. Zhao and Wang [17] established an ISRW dynamical model considering the medium as a subclass to explore the mechanism of spreading of individual-to-individual and medium-to-individual. Lebensztayn et al. [18] investigated the case where a new uninterested class of people exists. Z.-L. Zhang and Z.-Q. Zhang [19] studied a model for describing the interplay between rumor spreading and emergency.In particular, recently, a susceptible-infected-removed (SIR) rumor spreading model was introduced in [20]. This model considered a closed homogeneously mixing population of N individuals. The entire population is divided into three groups that are S, I, and R. Here S, I, and R stand for the people who are spreading rumor (spreaders, similar to infective), those who never heard rumor (ignorants, similar to susceptible), and the ones who heard rumor but do not spread it (stiflers, similar to removed), respectively.In the last decades, many rumor models have been given for modeling the spread of rumor. However, the above models play a very important role.The rumor spreading model introduced in [20] makes some assumptions about the rules of rumor spreading. For example, when an ignorant contacts a spreader, the ignorant may believe the rumor to become a spreader directly. This assumption ignores the possibility that the ignorant may experience a latent period before becoming spreader. This possibility lies in that, faced with mixed and conflicting messages, many people’s judgment becomes clouded, and some ignorants have difficulty to discriminate the authenticity of information immediately after they heard the rumor. In addition, there is no inflow to the susceptible class or outflow from any of the classes in their model. Still, their model was extremely innovative and is very useful in the modeling and analysis of rumor spreading.In order to make the rumor spreading process fit realistic cases more closely, we extend the above SIR rumor spreading model by considering the situation that when an ignorant contacts a spreader, a certain percentage of the ignorants may experience a latent period before becoming spreaders. In addition, considering the mobility of people in a certain area, we assume that there is inflow to the ignorant class and outflow from each class. Based on the above rumor spreading mechanism, we build a novel rumor spreading model called theXYZR model to further study the dynamics of rumors spreading.The rest of this paper is organized as follows. We describe theXYZR rumor spreading model and discuss the existence of equilibrium in Section 2. In Section 3, we carry out a qualitative analysis of the model. Stability conditions for the rumor-elimination equilibrium and the prevailing equilibrium are derived, respectively. In Section 4, we report the numerical results and discuss the implications of these results. Finally, conclusions are drawn in Section 5.


## 2.XYZR Rumor Spreading Model

We define ourXYZR rumor spreading model in the following way. Let N(t) denote the total population at time t. We divide the population into four classes: the ignorant class, the latent class, the spreader class, and the stifler class, according to the perception and reaction of an individual to a rumor. Population densities at time t are denoted by X(t), Y(t), Z(t), and R(t), respectively. Those who belong to the ignorant class, whom we call ignorants, do not know about the rumor. Those who belong to the latent class, whom we call latents, know about the rumor and cannot discriminate the authenticity of information immediately after he/she heard the rumor and require active effort to verify the information. The latent class also has infectious force due to the possibility that a latent who has the awareness of seeking confirmation with each other may share the information with the ignorants who are the latent’s friends or relatives immediately after he/she heard the rumor; on the other hand, the latent who does not understand the rumor due to the low level of his/her scientific knowledge or forgets about the rumor is probably not going to be sharing the rumor with his/her friends or relatives. In this case, we assume that an ignorant contacts a latent with probability λ∈(0,1]. Those who belong to the spreader class, whom we call spreaders, know about the rumor and spread it actively. With the development of emergency, a spreader may spontaneously lose interest in the rumor if the rumor itself is not absorbing. In that case, we assume that an ignorant contacts one of these spreaders with probability 1-λ. Those who belong to the stifler class, whom we call stiflers, know about the rumor and do not spread it. Each of which we call rumor class. By definition, we have
(1)N(t)=X(t)+Y(t)+Z(t)+R(t).We assume that no transition of rumor class happens unless a spreader and latent contact someone, and there exists infectious force in the latent class and spreader class. TheXYZR rumor spreading model is depicted in the transfer diagram in Figure 1.Figure 1
The flowchart ofXYZR model for the rumor spreading process.First, we consider the spread of rumor with variable population size; we assume that the rumors spread in the population with constant immigration and emigration. All recruitment is into the ignorant class and occurs at a constantε, and k1, k2, k3, and k4 are the emigration rates of those in the ignorant class, the latent class, the spreader class, and the stifler class, respectively. We assume that ε, k1, k2, k3, and k4 are positive constants, and that emigration is independent of rumor class.As shown in Figure1, the rumor spreading rules of our XYZR model can be summarized as follows.(i)
When an ignorant contacts a latent with probabilityλ∈(0,1], the ignorant becomes a latent at a rate β, namely latent rate.(ii)
When an ignorant contacts a spreader with probability1-λ, the ignorant believes the rumor and becomes a spreader directly at a rate μ, namely spreading rate.(iii)
Because some ignorants who have expertise or scientific knowledge (such as scientists or experts in the related fields) will not be affected by the rumors, hence, these ignorants directly become stiflers at a rateσ, namely stifling rate.(iv)
Any spreader may automatically lose interest in spreading, and consequently becomes a stifler at rateθ, namely removed rate.For the meantime, we assume the rumor is “constant”; that is, the same remark is transmitted at all times. We assume that(1-λ)(βY+μZ)X and λX(βY+μZ) are ignorants that change their rumor class and become spreaders and latents, respectively. When the latents contact the spreaders, some of the latents may firmly believe the rumor and spread the rumor to the ignorants, and so, we assume that αYZ latents change their rumor class and become spreaders, where α is a positive constant number representing the proportion of latents changing their rumor class.TheXYZR model was established based on the above assumptions according to the dynamic theory. The differential equations are described as follows:
(2)dXdt=ε-k1X-(βY+μZ)X-σX,dYdt=λ(βY+μZ)X-k2Y-αYZ,dZdt=(1-λ)(βY+μZ)X-(k3+θ)Z+αYZ,dRdt=θZ+σX-k4R,dNdt=dXdt+dYdt+dZdt+dRdt=ε-k1X-k2Y-k3Z-k4R.
From the first equation of system (2), we have that
(3)dXdt≤ε-(k1+σ)X.
Therefore
(4)limsupX≤t→∞εk1+σ.
From the last equation of system (2), we can see that
(5)dNdt=dXdt+dYdt+dZdt+dRdt=ε-k1X-k2Y-k3Z-k4R≤ε-k¯N,
where k¯=min{k1,k2,k3,k4}, thus, it can be obtained that
(6)limsupN≤t→∞εk¯.
The equation for R can be omitted in our analysis as R does not appear in any of the other equations. This shows that the model can be studied in the feasible region.Consider the following:(7)Θ={(X,Y,Z)∈R+3:X≤εk1+σ,X+Y+Z≤εk¯},
where R+3 denotes the nonnegative cone and its lower dimensional faces. It can be verified that Θ is positively invariant with respect to (2).Once the dynamics of(X,Y,Z) are understood, those of R can then be determined from the equation dR/dt=θZ+σX-k4R. Then the reduced limiting dynamical system is given by
(8)dXdt=ε-k1X-(βY+μZ)X-σX=f1(X,Y,Z),dYdt=λ(βY+μZ)X-k2Y-αYZ=f2(X,Y,Z),dZdt=(1-λ)(βY+μZ)X-(k3+θ)Z+αYZ=f3(X,Y,Z).
We first consider the existence of equilibria of system (8).Denote the boundary and the interior ofΘ in R+3 by ∂Θ and Θο, respectively. The rumor-elimination equilibrium occurs when there are no latents and spreaders in the system; that is, Y=0 and Z=0, which means that the rumors die out.For any values of parameters, system (8) always has a rumor-elimination equilibrium ϕ0=(ε/(k1+σ),0,0).A rumor-prevailing equilibriumϕ*=(X*,Y*,Z*) satisfies X*,Y*,Z*>0. To find the rumor-prevailing equilibrium (or positive equilibrium) ϕ*=(X*,Y*,Z*), set
(9)f1(X,Y,Z)=ε-k1X-(βY+μZ)X-σX=0,f2(X,Y,Z)=λ(βY+μZ)X-k2Y-αYZ=0,f3(X,Y,Z)=(1-λ)(βY+μZ)X-(k3+θ)Z+αYZ=0.
From (9), we can show that a unique ϕ* exists with
(10)X*=(k3+θ)λk3β+(k2+α)μ+λ(θ-k2μ).
From the definition of feasible region (7), we know that ε/((k1+σ)X*)≥1.Define(11)R0=ε(k1+σ)X*=ε(λk3β+(k2+α)μ+λ(θ-k2μ))(k1+σ)(k3+θ).
Then R0 is a threshold parameter that determines the number of equilibria.From the mathematical analysis of (9) we can see that(1)
ifR0≤1, then ϕ0 is the only equilibrium in Θ;(2)
ifR0>1, then there are two equilibria, ϕ0 and a unique rumor-prevailing equilibrium ϕ*=(X*,Y*,Z*).In the next section, we will study the property of these equilibria and perform a global qualitative analysis of model (8).


## 3. Dynamical Analysis ofXYZR Rumor Spreading Model

Theorem 1.
SupposeR0<1, the rumor-elimination equilibrium ϕ0 of system (8) is locally asymptotically stable.Proof.
The Jacobian matrix of system (8) at ϕ0=(ε/(k1+σ),0,0) is(12)J(ϕ0)=[-k1-δ-β(εk1+σ)-μ(εk1+σ)0λβ(εk1+σ)-k2λμ(εk1+σ)0(1-λ)β(εk1+σ)(1-λ)μ(εk1+σ)-(k3+θ)].One eigenvalue of J(ϕ0) is λ1=-(k1+δ)<0. The other two eigenvalues λ2, λ3 are eigenvalues of the 2 × 2 matrix. 							(13)L=[λβ(εk1+σ)-k2λμ(εk1+σ)(1-λ)β(εk1+σ)(1-λ)μ(εk1+σ)-(k3+θ)].
Then we obtain
(14)tr(L)=λβ(εk1+σ)-k2+(k3+θ)×[(1-λ)μ(ε/(k1+σ))(k3+θ)-1].
Using our assumption that R0<1, we have
(15)λβ(εk1+σ)-k2<0,(1-λ)μ(ε/(k1+σ))(k3+θ)-1<0.
Thus, we get
(16)tr(L)=λβ(εk1+σ)-k2+(k3+θ)×[(1-λ)μ(ε/(k1+σ))(k3+θ)-1]<0.
Simple calculations show that
(17)det(L)=[λβ(εk1+σ)-k2]×[(1-λ)μ(εk1+σ)-(k3+θ)]-λμ(εk1+σ)(1-λ)β(εk1+σ)=ε(k2+α)(k1+σ)(1-R0).
Therefore, det(L)>0 if and only if R0<1. Hence, by Routh-Hurwitz criterion, the rumor-elimination equilibrium ϕ0 of system (8) is locally asymptotically stable if R0<1. This completes the proof.Next, we focus on investigating the global stability of the rumor-elimination equilibrium and rumor-prevailing equilibrium. Firstly, we consider the global stability of the rumor-elimination equilibrium point.Theorem 2.
The rumor-elimination equilibriumϕ0 is(i)
globally asymptotically stable in the feasible regionΘ if R0≤1,(ii)
unstable ifR0>1.Proof.
To prove the global asymptotic stability ofϕ0 we use the method of Lyapunov functions. Motivated by those used in [21], we take a Lyapunov function ψ of form
(18)ψ=(βk2+μαk2(k3+θ))Y+μk3+θZ.
Then
(19)dψdt=(βk2+μαk2(k3+θ))dYdt+μk3+θdZdt=[(βk2+μαk2(k3+θ))]×[λ(βY+μZ)X-k2Y-αYZ]+μk3+θ×[(1-λ)(βY+μZ)X-(k3+θ)Z+αYZ]=(k1+δεR0X-1)(βY+μZ).
Since R0≤1 and X≤ε/(k1+σ), we have
(20)dψdt=(k1+δεR0X-1)(βY+μZ)≤(R0-1)(βY+μZ)≤0.
Furthermore, dψ/dt=0⇔Y=Z=0 or R0=1 and X=ε/(k1+σ). Therefore the largest invariant set in the closure Θ¯ of Θ where dψ/dt=0 is the singleton {ϕ0}. By LaSalle’s invariance principle [21], ϕ0 is globally asymptotically stable in Θ, and the proof is completed.Next, we focus on investigating the global stability of the rumor-prevailing equilibriumϕ* when R0>1.Theorem 3.
IfR0>1, then the rumor-prevailing equilibrium ϕ* of the system (8) is globally stable in Θο.Proof.
By Routh-Hurwitz criterion, to prove the global stability ofϕ*, we need to prove that all eigenvalues of corresponding Jacobian matrix J of system (8) have negative real parts. Then we want to show, when R0>1, that the Routh-Hurwitz conditions hold, namely, tr(J)<0, det(J)<0, and tr(J)*M-det(J)<0, respectively.
The Jacobian matrix of system (8) at ϕ*  is(21)J(ϕ*)=[-k1-(βY*+μZ*)-δ-βX*-μX*λ(βY*+μZ*)λβX*-k2-αZ*λμX*-αY*(1-λ)(βY*+μZ*)(1-λ)βX*+αZ*(1-λ)μX*-(k3+θ)+αY*].Simple calculations show that
(22)tr(J(ϕ*))=-k1-(βY*+μZ*)-δ+λβX*-k2-αZ*+(1-λ)μX*-(k3+θ)+αY*<0.
Now we calculate
(23)det(J(ϕ*))=(-k1-δ-βY*-μZ*)×(λβX*-k2-αZ*)×(μX*-λμX*-k3-θ+αY*)+(-βX*)(λμX*-αY*)×((1-λ)(βY*+μZ*))+(-μX*)((1-λ)βX*+αZ*)×λ(βY*+μZ*)<0.
Denote M the sum of determinants of second-order principal minors of matrix J(ϕ*), thus we have
(24)M=(-k1-(βY*+μZ*)-δ)×(λβX*-k2-αZ*)+βX*λ(βY*+μZ*)+(-k1-(βY*+μZ*)-δ)×((1-λ)μX*-(k3+θ)+αY*)+μX*(1-λ)(βY*+μZ*)+(λβX*-k2-αZ*)((1-λ)μX*-(k3+θ)+αY*)-(λμX*-αY*)((1-λ)βX*+αZ*).
Hence, according to (22)–(24), we get
(25)tr(J(ϕ*))*M-det(J(ϕ*))<0ifR0>1.
By (22), (23), (25), and Routh-Hurwitz criterion, we know that all eigenvalues of Jacobian matrix J(ϕ*) have negative real parts, and this completes the proof.


## 4. Discussions and Simulations

In this section, we report the numerical results and discuss the implications of these results. We carry out numerical simulations of our model (2) by using the Runge-Kutta method and the software MATLAB 7.0. First, the obtained results of theoretical analysis (Theorems 2 and 3) are numerically validated (Figures 2 and 3).Figure 2
Whenε=3, k1=0.0005, k2=0.0002, k3=0.04, β=0.0005, k4=0.2, μ=0.00005, σ=0.85, α=0.25, θ=0.001, λ=0.06, and R0=0.2725<1.Figure 3
Whenε=3, k1=0.00005, k2=0.002, k3=0.004, β=0.00015, λ=0.6, k4=0.0002, μ=0.00005, σ=0.0015, α=0.0015, θ=0.0001, and R0=28.5484>1.Numerical results carried out for the system (2) show that the rumor “dies out” when the threshold parameter R0<1 (Figure 2), and the rumor persists at a “prevalent” level when R0>1 (Figure 3).Recall the threshold parameterR0, we have
(26)R0=ε(k1+σ)X*=ε(λk3β+(k2+α)μ+λ(θ-k2μ))(k1+σ)(k3+θ).It is straightforward from (26) thatR0increases as μ or β increases, and R0 decreases as σ increases. In reality, if there exist many latents and spreaders, which means that the parameters μ and β will be sufficiently large such that R0>1, then the rumors will persist at a “prevalent” level. Conversely, if the government departments made an immediate response which includes cracking down the rumor spreaders, promulgating authoritative information through authoritative channels, and improving the ignorants’ discrimination, then the parameters μ and β will become sufficiently small and σ becomes sufficiently larger such that R0<1, and the rumor will “die out.”Furthermore, from the perspective of crisis management strategy, it should be noted that the government departments cracking down the rumor spreaders will lead to the decrease of the spreading rateμ, the timely disclosure of transparent information will result in the decrease of the latent rate β, and the popularization of scientific knowledge and the improvement of the ignorants’ discrimination will lead to the increase of the stifler rate σ.The above analysis indicates thatμ, β, and σ that appear in the threshold parameter R0 are three key parameters in the XYZR rumor spreading model. Therefore, to offer useful suggestions to crisis manager, we will carry out numerical simulations by varying the above three parameters to investigate the impact corresponding management strategies have on the rumor spreading process, respectively, and discuss these simulation results with the “Chinese salt-buying frenzy case.”Figure4 shows how the size of the spreader class changes over time for different spreading rates μ with the other parameters fixed. From Figure 4, we can see that the lower the spreading rate μ, the smaller the number of rumor spreaders class. This is expected because as the spreading rate μ decreases, there is lower spreading capacity for spreaders, which means fewer chances for spreaders to be in contact with ignorants, making the spreading of a rumor much harder. Consider the Chinese salt-buying frenzy case, Chinese government and media showed a rapid and effective response in coping with the panic, such as Chinese public security authorities quickly detaining several people for allegedly fabricating and spreading online rumors about “a salt shortage” after tsunami crippled Japan’s Fukushima Daiichi nuclear power plant triggered panic buying and hoarding of table salt in China, which timely controlled the source and spread of the rumors, and the salt-buying frenzy was under control effectively. One possible interpretation on this is that, in reality, if there is no strict legal punishments on those caught spreading rumors, the negative factors will run wild to destroy the social order and even pose a threat to social stability, which will bring great damage to people and society. This reflects the reality that clearing up the roots and causes of the rumors by increasing the efforts to crack down on rumor spreaders is crucial to curbing rumors.Figure 4
Size of the spreader class over time under different spreading rateμ with ε=3, k1=0.00005, k2=0.002, k3=0.004, k4=0.0002, β=0.00015, σ=0.0015, α=0.015, θ=0.0001, and λ=0.6.Figure5 describes how the size of the spreader class changes over time for different latent rate β with the other parameters fixed. As shown in Figure 5, the lower the latent rate β, the smaller the number of rumor spreaders class. Therefore, the new added factors have affected the process of rumor spreading. With a lower β value, fewer ignorants become the latents. Furthermore, it should be noted that some of these latents can transform to the spreaders when they believe the rumors they heard, and this process postpones the time of the spreaders’ final disappearance. This implies that authorities should release information with detailed facts immediately after rumors begin to circulate, because the preciseness and openness of crisis information was helpful to dispel public panic. Consider the Chinese salt-buying frenzy case, in the initial phase of the event, faced with mixed and conflicting messages, many people’s judgment becomes clouded, and people are eager to acquire authentic information through authoritative channels. Moreover, the Chinese government’s transparent and active approach to releasing information regarding the nuclear leak accident in Japan should be praised. The relevant government departments immediately worked out policies to stabilize the salt market by increasing supply and strengthening supervision. This suggests that getting a timely response can minimize the hazards.Figure 5
Size of the spreader class over time under different latent rateβ with ε=3, k1=0.00005, k2=0.002, k3=0.004, k4=0.0002, μ=0.00005, σ=0.0015, α=0.015, θ=0.0001, and λ=0.6.Figure6 illustrates how the size of the spreader class changes over time for different stifling rate σ with the other parameters fixed. From Figure 6, we can see that the higher the stifling rate σ, the smaller the peak value of the spreader class, and the higher the stifling rate is σ, the faster the rumor terminates. As shown in Figure 6, the number of spreaders shows a much greater decline while the value of σ increases from 0.0015 to 0.15. With a higher σ value, more ignorants (such as experts or doctors who know the professional knowledge of nuclear protection) directly go to the stiflers status and do not participate further in the spreading dynamics, which leads to smaller number of the spreaders. It should be noted that the increase of the σ value mainly depends on the popularization of scientific knowledge and the improvement of the ignorants’ discrimination. Moreover, the average level of people’s scientific knowledge is, to some degree, very low, which leads to easily believing in rumors. Only 3.27 percent of Chinese population has basic scientific knowledge, according to a survey conducted by the China Science and Technology Museum in 2010. This reflects the reality that the popularization of scientific knowledge is needed, and related experts and departments should undertake due obligations. The results coincide with the fact of the evolution of hoarding iodized salt. Consider the Chinese salt-buying frenzy case, due to a few people’s lack of basic scientific knowledge about radiation, they rushed to buy iodized salt in panic after the nuclear leak in Japan. However, it is worth mentioning that Chinese government and media showed a rapid and effective response in coping with the panic. Besides the timely and transparent release of the monitoring results, the country’s media outlets also began to spread scientific knowledge about radiation, such as its definition and impact on people’s health. State-owned China Central Television kept inviting specialists almost every day to explain the facts about nuclear power to the public. After the authorities announced that the rumor was absurd and malicious then the rumor perished rapidly. It is worth mentioning that this time the panic buying salt relay panic buying of Radix Isatidis in China, signifying the popularization of scientific knowledge of shoulder heavy responsibilities.Figure 6
Size of the spreader class over time under different stifling rateσwith ε=3, k1=0.00005, k2=0.002, k3=0.004, k4=0.0002, β=0.00015, μ=0.00005, α=0.015, θ=0.0001, and λ=0.6.


## 5. Conclusion

In this paper, we considered some more realistic factors and extended the SIR rumor spreading model in [20] to make it closer to real life. Our stability analysis approach is quite different from that adopted in [20]. We have carried out a global qualitative analysis of a XYZR rumor spreading model with latent and constant recruitment and varying total population and studied the existence and stability of the rumor-elimination and prevailing equilibria and the impact different management strategies have on the rumor spreading process.We derived that the dynamics of system (2) are completely determined by the threshold parameter R0. In terms of the threshold parameter R0, our main results indicate that when R0<1, the rumor-elimination equilibrium ϕ0 is globally asymptotically stable in the feasible regionΘ, and thus that the spread of rumor will die out spontaneously (Figure 2). When R0>1, the rumor-prevailing equilibrium ϕ* exists and is globally stable in Θο, so that the rumor, if initially present, will persist at the unique prevailing equilibrium level (Figure 3).From the perspective of crisis management, the numerical results imply that increasing the efforts to crack down on rumor spreaders, the timely disclosure of transparent information and the popularization of scientific knowledge and the improvement of the ignorants’ discrimination may be effective rumor control measures.



---

*Source: 605854-2013-12-09.xml*