# Effects of awareness diffusion and self-initiated awareness behavior on epidemic spreading - an approach based on multiplex networks

## Abstract

In this paper, we study the interplay between the epidemic spreading
and the diffusion of awareness in multiplex networks. In the model,
an infectious disease can spread in one network representing the
paths of epidemic spreading (contact network), leading to the diffusion of awareness in the other network (information network), and then the diffusion
of awareness will cause individuals to take social distances, which in turn affects the epidemic spreading. As for the diffusion of awareness, we
assume that, on the one hand, individuals can be informed by other aware neighbors in information network, on the other hand, the susceptible individuals can be
*self-awareness* induced by the infected neighbors in the
contact networks (local information) or mass media (global
information). Through Markov chain approach and numerical
computations, we find that the density of infected individuals and the
epidemic threshold can be affected by the structures of the two
networks and the effective transmission rate of the
awareness. However, we prove that though the
introduction of the *self-awareness* can lower the
density of infection, which cannot increase the epidemic threshold
no matter of the local information or global information. Our
finding is remarkably different to many previous results–local information
based behavioral response can alter the epidemic threshold.

###### keywords:

Multiplex netwoeks, Infectious diseases, Awareness diffusion, Epidemic threshold.elsarticle-num

^{1}

## 1 Introduction

The outbreaks of diseases can involve the diffusion of information in regard to the diseases, including the risk of
infection, rumors, fears and so on, which can stimulate
individuals to take spontaneous behavioral responses to protect
themselves, thereby bring profound impacts on the spreading of
disease [1, 2, 3, 4, 5, 6]. For example, recent outbreaks of the H1N1 flu, the bird flu, and the severe acute respiratory syndrome
(SARS) have brought the reduction of going out and the plenty of
people wearing face masks. For
this reason, there has been an increasing focus on the development
of formal models aimed at investigating the interplay of epidemic
spreading and information-based behavioral responses[7, 8, 9, 10]. Such as, based
on the assumption that the probability of susceptible individual
going to the alter state is proportional to the number of infected
neighbors, Sahneh*et al.* extended the SIS
(Susceptible-Infected-Susceptible) model to a
Susceptible-Alter-Infected-Susceptible (SAIS) model [11, 12],
and they found that the way of behavioral response can enhance the
epidemic threshold; Meloni *et al.* studied a meta-population
model that incorporates several scenarios of self-initiated
behavioral changes into the mobility patterns of individuals, and
they found that such behavioral changes do not alter the epidemic
threshold, but may produce a negative impact on disease, i.e., the
density of infection is increased [13]; In Refs. [14, 15, 16], authors
investigated the effects of the information-based behavioral
responses on the epidemic dynamics by designing the transmission
rate as a function of the local infected density or the global
infected density.

Though the effects of information-based behavioral responses on the
epidemic dynamics have been studied by many authors, most of works
assumed the spreadings of information and epidemic are in the same
network. As we know, with the development of technology, information
can fast diffuse through many different channels, such as, the word
of mouth, news media, online social networks, and so on. In view of
this, recent well-studied multiplex network theory has been used
to mimic the interplay of information or related awareness and the
epidemic dynamics[17, 18, 19, 20].
For instance, Sahneh *et al.* have shown
that the information dissemination spread in another network can
help boost the resilience of the agents’ population against the
spreading and found optimal information dissemination for different
topologies [21]; Wang *et al.* investigated the
interplay of the epidemic dynamics and the information dynamics in
multiplex network based on the SIR (R-recovery) model, and focused on
the two fundamental quantities underlying any spreading process:
epidemic threshold and the final epidemic prevalence [22]; Granell
*et al.* established an SIS-UAU model to investigate the
competing effects of the spreading of awareness and the epidemic
dynamics in multiplex with the transmission rate of awareness as
well as the structure of information network [23]. More recently, they
further generalized their model by reducing the probability of
infected individuals becoming awareness and including the effect
of a mass broadcast of awareness (mass media) on the epidemic
dynamics [24].

Existing works either assume that individuals are self-initiated [13, 12],
that is, they become aware because their neighbors are infected, or
individuals can only be informed and become aware by other aware
neighbors [22, 23], the combine effects of the two factors have not been well
studied. In reality, individuals can obtain the disease information through many ways. Inspired these factors, in the current work, we study the
interplay between the diffusion of awareness by incorporating the
*self-awareness* effects and the epidemic dynamics under the
framework of multiplex network. In the model, an infectious disease
first spreads among population represented by the contact network,
and then the outbreak of the disease stimulates some people
(infected or susceptible individuals) become aware of the risk of
infection, and they take some protections to reduce the probability of
infection. Meanwhile, unaware individuals can be informed by other
aware individuals through the information network or become *self-awareness*
induced by the infected neighbors in contact network or mass media. The finding
indicates that the additional self-initiated awareness mechanism can
reduce the density of infection, however, which can not alter the epidemic threshold. The results are
verified by the Monte-Carlo simulations and the microscopic Markov
chain approach (MMCA).

The layout of the paper is as follows: we introduce the model in
Sec. 2. The simulation results and theoretical
analysis are presented in Sec. 3. Finally,
Conclusions and discussions are presented in
Sec. 4. The results for the global information-based *self-initiated * awareness are given in Appendix.

## 2 Model

In this work, we generalize the model of Ref. [23]. In that model, a multiplex includes two layers, one is physical layer representing the spreading of epidemic (contact network), and the other is information layer on where the diffusion of the awareness evolves (information network). All nodes represent the same individuals in both layers, but the connectivity is different in each of them. In the contact layer, a Susceptible - Infected - Susceptible (SIS) model is used to mimic the epidemic dynamics. That is to say, a susceptible node can be infected by one infected neighbor with certain probability, and the infected node can return to susceptible state with probability . On the information layer, the dynamical process of awareness is assumed to be similar to the SIS model, that is, an unaware node (U) can be informed by an aware neighbor (A) with probability , and the aware node can loss awareness and back to unaware state with probability . The interplay of the two processes is modelled as follows: once an individual is infected, s/he will certainly become aware, that is, the probability is . In addition, to distinguish the protective behaviors between the aware individuals and unaware individuals, let and (here . If , the aware individuals are completely immune to the infection.) be the probabilities of unaware and aware susceptible nodes to get infected, respectively.

From the description of the model, one can find that, on the one hand,
the authors assumed that the infected individuals will
*automatically* become aware and are willing to inform the disease information. As we know, in many cases, infected individuals are
unwilling to tell others since they can be discriminated or isolated
by others once others know they are infected by one certain disease. So we make a progress and
assume that infected individuals becoming aware with probability
. On the other hand, in the model, individuals can
*only* be informed by their neighbors through the information
network, that is to say, one individual has no chance to become aware
once their neighbors are unaware. However, as proposed in many
previous works, individuals can become *self-awareness* once their friends are
infected or they are informed by the mass media. Thus, we further
generalize the model as: a susceptible individual can go to aware
state by *self-initiated* response with probability
when contacting one infected friend. Therefore, the probability of the susceptible
individual becoming awareness increases with the number of infected
neighbors in the contact network [12]. Note that, for the
original model in Ref. [23], awareness cannot break out if the probability , so the roles of awareness cannot be played,
however, in our model, the awareness can diffuse among population
even though since susceptible individuals can become awareness by their
*self-initiated* responses. When
and , our model returns to the original model in
Ref. [23].

According to this scheme, an individual can be in four different states: susceptible and unaware(SU), susceptible and aware(SA),infected and unaware(IU), infected and aware(IA). The flow diagram of the model is given in Fig. 1.

## 3 Main results

### 3.1 Theoretical analysis

Denoting and be the adjacency matrices that support the SIS and UAU processes, respectively.
The probability of in one of four states at time is denoted
by , , and
respectively. Assuming the probability of susceptible
(infect) node *not* being *informed* by any neighbors is
(), and the probability of unaware
(aware) susceptible node *not* being *infected* is
(). They are described as [23]:

(1) | |||||

(2) | |||||

(3) | |||||

(4) |

where and (Note: to simplify the model, we do not distinguish the infectivity of IA and IU, meanwhile, the diffusion capabilities of awareness for SA and IA are also the same.).

For each possible state at time , which may give rise to four possible states at time with certain probability, the transition probability trees for node are illustrated in Fig. 2. According to the scheme in Fig. 2, one can easily write the Markov Chain Approach (MMCA) equations [25, 26] for each state, for example, along the top branch of the four trees, one can read the probability of SU at as:

(5) |

Similarly, the other three MMCA equations can be written as:

(6) | |||||

(7) | |||||

(8) | |||||

where . When the system goes to the stationary state, we have for SU state and equivalently for the SU, SA and IU states.

Since the epidemic threshold determines whether the epidemic can outbreak or die out, it is vital to analyze the effects of the different parameters on the epidemic threshold . As stated in Ref. [23], near the threshold, the probability of nodes being infected is very low, i.e., . Consequently, and . Further approximate and by assuming , and then substitute both of them into Eq. (5) and Eq. (6), we obtain

(9) |

(10) |

Combing Eqs. (7–10), then a simple formula is obtained:

(11) |

With , , and , then Eq. (11) can be rewritten as:

(12) |

here if ; otherwise, .

Defining matrix with elements:

(13) |

Eq. (12) can be read as

(14) |

where with be the vector transportation.

The non-trivial solutions of Eq. (14) are eigenvectors of ,
whose eigenvalues are equal to . Therefore, the onset of
the epidemics is given by the largest eigenvalue of ,*i.e.*,
,

(15) |

### 3.2 Numerical simulations

To verify our theoretical results, as in Ref. [23], we build a configuration network with degree distribution and network size as the contact network, and for the information network, which is generated by adding 800 extra random links in the contact network. and represent the adjacency matrices of the contact network and the information network, respectively. All simulation results are obtained by averaging 20 realizations.

We first compare the results from MMCA with Monte-Carlo simulation in Fig. 3 to check the effectiveness of our analysis based on MMCA, from Fig. 3, one can observe that the results based on the two approaches are in good agreement. So in the next figures, our main results are obtained from MMCA.

Then we investigate the effects of the two main parameters of the model— and on the epidemic threshold and the density of infected individuals. Here, we will present the results for , meaning that and . Obviously, once the value of is increased, the epidemic threshold is decreased and the density of infected individuals in enhanced correspondingly.

Fig. 4 plots the density of infection [, see Fig. 4(a)] and aware
individuals [, see Fig. 4(b)] as a function of for different values of
, respectively. Observing
the figure, one can see that though the larger value of can cause more individuals become aware and reduce the density of infection. However, one can find that the increasing of has no influence on the epidemic threshold. The result is remarkably different from many previous results
which claim that the local information-based behavioral response in the *single-layer network* can alter the epidemic threshold. How to understand this nontrivial result? Since the UAU awareness dynamic is the same to the SIS epidemic process. When the epidemic has not broken out, the density of awareness in information network (*i.e.*, ) is only determined by the transmission rate of awareness, , recovery rate and the structure of information network, but is not related to the value of . Namely, near or below the epidemic threshold point, increasing the value of only means that the initial number of the aware individuals is increased, which cannot affect the density of aware individuals at stationary state. In this case, the value of is independent of the value of or , which gives rise to the same value of [see Fig. 4(b)]. Thus, according to Eqs. (13) and (15), the epidemic threshold is invariable owing to the same value of . One should note that,
for the case of single-layer network, the local information-based behavioral response can *directly* reduce the transmission rate of epidemic, leading to the change of the epidemic threshold. For our model, the self-awareness behavior first diffuses through the information network, and then the epidemic process happens in contact network. Thus, the effect of the awareness behavior on the transmission rate is *indirect*.

The density of infected individuals and aware nodes as functions of for different values of parameter are also shown in Fig. 5(a) and Fig. 5(b), respectively. One can see that, similar to Fig. 4, varying the value of has no effect on the epidemic threshold. However, differ to the above case, from Fig. 5(a) we find that the value of also has negligible effect on the density of infection, even in the extreme cases where infected unaware individuals remain unaware of its sickness () or certainly become aware of it (). The result in Fig. 5 is in accord with the Fig. 3 in Ref. [24].

In order to systematically study the effects of and on the , we further explore the full phase diagram () in Fig. 6. Overall, we can see that is not influenced by when is smaller than the epidemic threshold, since epidemic will die out by itself. Once overpasses the epidemic threshold, decreases with for different values of or . More specifically, by comparing Fig. 6(a) with Fig. 6(b) (or comparing Fig. 6(c) with Fig. 6(d)), we can see that is not remarkably influenced by the value of . Likewise, by comparing the Fig. 6(a) with Fig. 6(d) (or comparing Fig. 6(b) with Fig. 6(c)), decreases with , especially for the large value of .

## 4 Conclusions and discussions

Recognizing that, on the one hand, outbreak of an epidemic through a physical-contact network can trigger the spreading of information awareness through other different channels, such as on-line social networks, mass media; on the other hand, an individual can not only be informed by other aware neighbors but also can become self-awareness once some friends in contact network are infected. By introducing the *self-awareness* mechanism for susceptible individuals, we have investigated interplay between the spreading of epidemic and the diffusion of awareness based on the framework of the multiplex networks. We mainly studied the two parameters and characterizing the self-awareness probability of susceptible individuals and infected individuals, respectively. Analysis based on the Markov chain approach as well as the extensive computations reveal that the density of infection can be reduced once the two parameters are increased, however, we found that the impact of self-awareness behavior for susceptible individuals on inhibiting the spreading of epidemic is much better than the self-awareness of the infected individuals, since self-awareness from susceptible individuals can directly reduce their probabilities of being infected. Meanwhile, we found that the self-awareness behavior cannot alter the epidemic threshold no matter of the local or global information, which are in stark contrast with the results obtained from the single-layer networks.

The challenges of studying the intricate interplay between social and biological contagions in human populations are generating interesting science. In this work, we consider the effects of the *self-awareness* behavior based on the multiplex networks on the density of infection and the epidemic threshold, our result implies that the conclusions obtained from single-layer networks may need to be re-examined when they are extended to multiplex networks.

## 5 Acknowledgments

This work is funded the National Natural Science Foundation of China (Grant Nos. 61473001,11331009) and the Doctoral Research Foundation of Anhui University (Grant No. 02303319).

## 6 Apendix: Global information-based awareness

In recent work, Granell *et al.* considered the effect of the mass media on the epidemic process and awareness diffusion [24]. In the model, each individual
becomes aware with probability by assuming that they are informed by a broadcast or mass media. Thus, it can be regarded as a global information-based awareness. One questionable point is that the probability of being awareness is irrelevant to the density of infection. As a result, even the epidemic is almost eliminated, individuals also have the fix probability of being aware. In reality, becoming awareness often means that individuals need to take some protective measures, such as, washing hands, wearing masks or reducing outgoings. These measures indicate certain inconveniences or some cost [27, 28]. Thus, a more realistic situation is that the probability of being awareness should be related to the density of infection. To mimic this case, here we assume the probability of aware from global information is given as: with is the density of infection at time , which indicates that the probability of being awareness adaptively varies with the density of infection.

For this case, we only need to slightly change the local model described in subsection 3.1. We only need to change and as follows:

(16) |

(17) |

### Footnotes

- Corresponding author: haifengzhang1978@gmail.com

### References

- S. Funk, E. Gilad, C. Watkins, V. A. Jansen, The spread of awareness and its impact on epidemic outbreaks, Proc. Natl. Acad. Sci. U.S.A. 106 (16) (2009) 6872–6877.
- S. Kitchovitch, P. Liò, Risk perception and disease spread on social networks, Procedia Comput. Sci. 1 (1) (2010) 2345–2354.
- H.-F. Zhang, Z.-X. Wu, M. Tang, Y.-C. Lai, Effects of behavioral response and vaccination policy on epidemic spreading-an approach based on evolutionary-game dynamics, Sci. Rep. 4 (2014) 5666.
- Y. Bu, S. Gregory, H. L. Mills, Efficient local behavioral-change strategies to reduce the spread of epidemics in networks, Phys. Rev. E 88 (4) (2013) 042801.
- S. Funk, M. Salathé, V. A. Jansen, Modelling the influence of human behaviour on the spread of infectious diseases: a review, J. R. Soc. Interface 7 (50) (2010) 1247–1256.
- H. Yang, M. Tang, H.-F. Zhang, Efficient community-based control strategies in adaptive networks, New. J. Phys. 14 (12) (2012) 123017.
- F. Bagnoli, P. Liò, L. Sguanci, Risk perception in epidemic modeling, Phys. Rev. E 76 (6) (2007) 061904.
- L. Wang, Y. Zhang, Z. Wang, X. Li, The impact of human location-specific contact pattern on the sir epidemic transmission between populations, Int. J. Bifurcat. Chaos 23 (05).
- L. Wang, Z. Wang, Y. Zhang, X. Li, How human location-specific contact patterns impact spatial transmission between populations?, Sci. Rep. 3 (2013) 1468.
- Z. Ruan, M. Tang, Z. Liu, Epidemic spreading with information-driven vaccination, Phys. Rev. E 86 (3) (2012) 036117.
- F. D. Sahneh, C. M. Scoglio, Epidemic spread in human networks, in: Decision and Control and European Control Conference (CDC-ECC), Orlando, Florida, 2011, pp. 3008–3013.
- F. D. Sahneh, F. N. Chowdhury, C. M. Scoglio, On the existence of a threshold for preventive behavioral responses to suppress epidemic spreading, Sci. Rep. 2 (2012) 632.
- S. Meloni, N. Perra, A. Arenas, S. Gómez, Y. Moreno, A. Vespignani, Modeling human mobility responses to the large-scale spreading of infectious diseases, Sci. Rep. 1 (2011) 62.
- Q. Wu, X. Fu, M. Small, X.-J. Xu, The impact of awareness on epidemic spreading in networks, Chaos 22 (1) (2012) 013101.
- Q. Wu, H. Zhang, G. Zeng, Responsive immunization and intervention for infectious diseases in social networks, Chaos 24 (2) (2014) 023108.
- H.-F. Zhang, J.-R. Xie, M. Tang, Y.-C. Lai, Suppression of epidemic spreading in complex networks by local information based behavioral responses, Chaos 24 (4) (2014) 043106.
- S. Gomez, A. Diaz-Guilera, J. Gomez-Gardeñes, C. J. Perez-Vicente, Y. Moreno, A. Arenas, Diffusion dynamics on multiplex networks, Phys. Rev. Lett. 110 (2) (2013) 028701.
- Z. Wang, A. Szolnoki, M. Perc, Interdependent network reciprocity in evolutionary games, Sci. Rep. 3 (2013) 1183.
- M. De Domenico, A. Solé-Ribalta, E. Cozzo, M. Kivelä, Y. Moreno, M. A. Porter, S. Gómez, A. Arenas, Mathematical formulation of multilayer networks, Phys. Rev. X 3 (4) (2013) 041022.
- S. Boccaletti, G. Bianconi, R. Criado, C. Del Genio, J. Gómez-Gardeñes, M. Romance, I. Sendiña-Nadal, Z. Wang, M. Zanin, The structure and dynamics of multilayer networks, Phys. Rep. 544 (1) (2014) 1–122.
- F. D. Sahneh, C. M. Scoglio, Optimal information dissemination in epidemic networks., in: CDC, 2012, pp. 1657–1662.
- W. Wang, M. Tang, H. Yang, Y. Do, Y.-C. Lai, G. Lee, Asymmetrically interacting spreading dynamics on complex layered networks, Sci. Rep. 4 (2014) 5097.
- C. Granell, S. Gomez, A. Arenas, Dynamical interplay between awareness and epidemic spreading in multiplex networks, Phys. Rev. Lett. 111 (12) (2013) 128701.
- C. Granell, S. Gómez, A. Arenas, Competing spreading processes on multiplex networks: Awareness and epidemics, Phys. Rev. E 90 (2014) 012808.
- S. Gómez, A. Arenas, J. Borge-Holthoefer, S. Meloni, Y. Moreno, Discrete-time markov chain approach to contact-based disease spreading in complex networks, EPL (Europhysics Letters) 89 (3) (2010) 38009.
- Y. Wang, D. Chakrabarti, C. Wang, C. Faloutsos, Epidemic spreading in real networks: An eigenvalue viewpoint, in: Reliable Distributed Systems, 2003. Proceedings. 22nd International Symposium on, IEEE, 2003, pp. 25–34.
- H.-F. Zhang, Z. Yang, Z.-X. Wu, B.-H. Wang, T. Zhou, Braess’s paradox in epidemic game: Better condition results in less payoff, Sci. Rep. 3 (2013) 3292.
- F. Fu, D. I. Rosenbloom, L. Wang, M. A. Nowak, Imitation dynamics of vaccination behaviour on social networks, P. Roy. Soc. B-Biol. Sci. 278 (1702) (2011) 42–49.