Epidemics on Networks: Reducing Disease Transmission Using Health Emergency Declarations and Peer Communication

TL;DR

This study examines how official health declarations and peer communication influence disease spread on different static social networks, revealing optimal intervention thresholds vary by network type and impact epidemic size.

Contribution

It introduces a stochastic model analyzing the effects of public health information dissemination on epidemic dynamics across various network structures.

Findings

Optimal awareness threshold P* reduces epidemic size in Erdos-Renyi and Small-world networks.

Threshold P* does not minimize epidemic size in Scale-free networks.

The basic reproduction number and threshold P* influence epidemic outcomes differently.

Abstract

Understanding individual decisions in a world where communications and information move instantly via cell phones and the internet, contributes to the development and implementation of policies aimed at stopping or ameliorating the spread of diseases. In this manuscript, the role of official social network perturbations generated by public health officials to slow down or stop a disease outbreak are studied over distinct classes of static social networks. The dynamics are stochastic in nature with individuals (nodes) being assigned fixed levels of education or wealth. Nodes may change their epidemiological status from susceptible, to infected and to recovered. Most importantly, it is assumed that when the prevalence reaches a pre-determined threshold level, P*, information, called awareness in our framework, starts to spread, a process triggered by public health authorities. Information is assumed to spread over the same static network and whether or not one becomes a temporary informer, is a function of his/her level of education or wealth and epidemiological status. Stochastic simulations show that threshold selection P* and the value of the average basic reproduction number impact the final epidemic size differentially. For the Erdos-Renyi and Small-world networks, an optimal choice for P* that minimize the final epidemic size can be identified under some conditions while for Scale-free networks this is not case.

Figures

Tables

Table 1. Table 1 : Parameters and their baseline values assumed for simulations, unless stated otherwise.

	Parameter	Description	Unit	Baseline
Network Parameters	N	Number of nodes in the Network	People	$10000$
	GE	Erdős-Rényi random network	–	$G(N,0.001)$
	GW	Small-world network	–	$G(N,10,0.5)$
	GS	Scale-free random network	–	$P\left(k\right)\sim k^{-2.11}$
	$\overline{C}$	Average number of contact per neighbor per unit time	contact/time	$0.75$
	$X\left(\overline{x}\right)$	Level of education distribution (its mean)	–	$\beta (5,2)\left(0.7\right)$
Infection Parameters	$\beta$	Probability of infection transmission per contact	1/contact	$0.011$
	$1/\gamma$	Average time to recover without treatment	days	$10$
Awareness Parameters	$\tilde{\beta }$	Average probability of awareness transmission per contact	1/contact	$0.3$
	$P^{\ast }$	Prevalence threshold	1	$0.1$
	$1/\tilde{\gamma }$	Average time of behavior change for susceptible individuals	days	$7$
	$\kappa$	Saturation factor in $\sigma$ function	1	$0.85$
	$\theta$	Parameter in which, half maximum $\sigma$ function is obtained	1	$0.5$
	$x^{\ast }$	Level of education threshold	1	$0.4$

Table 2. Table 2 : Table of notation for a conventional network G 𝐺 G in algorithms.

Notation	Description
$S_t$	Set of susceptible nodes at time t
$U_t$	Set of unaware nodes at time t
$A_t$	Set of aware nodes at time t
$C_t$	Set of careless nodes at time t
$I_t$	Set of infected nodes at time t
$Q_t$	Set of quarantine nodes at time t
$F_t$	Set of free nodes at time t
$R_t$	Set of recovered nodes at time t
$G.N\left(\text{k}\right)$	Set of neighbors of node k in Network G
$IP$(k)	Infection period for infected node k
$AP$(k)	Awareness period for aware node k
$ern\left(\alpha \right)$	Exponential random number with average $\alpha$ for $\alpha >0$
$urn$	Uniform random number in $[0,1]$
$c_{kj}$	Probability of having contact between two neighbors k and j
$A\underset{}{\overset{\text{k}}{\to }}B$	Element k moves from set A to set B

Full text

Keywords: Awareness spread, Behavior change, Outbreak and epidemic threats, Erdős-Rényi network, Small-world Network, Scale-free Network,

$\leavevmode\nobreak\ {}^{1}$ School of Human Evolution and Social Change; Simon A. Levin Mathematical Computational Modeling Science Center, Arizona State University, Tempe, AZ 85281

$\leavevmode\nobreak\ {}^{2}$ Division of Applied Mathematics, Brown University, Providence RI, 02906

$\leavevmode\nobreak\ {}^{3}$ Sciences and Mathematics Faculty, College of Integrative Sciences and Arts, Arizona State University, Mesa, AZ 85212, USA

1 Introduction

The decisions that individuals make over an epidemic outbreak depend on multiple factors. Here, they are assumed to depend on available information, misinformation, and the income/education of those making them [9, 15, 5, 36, 30]. There are multiple possible scenarios that consider the decisions that individuals may make over the course of an outbreak. Individuals may modify their behaviors in order to reduce their environmental susceptibility to a disease by washing their hands frequently, avoiding handshakes and avoiding kissing salutes, not taking public transportation during rush hours, using masks and more. The frequency and effectiveness of these decisions may depend on the perceived risk of infection, a function of what each individual “knows”. In short, individual responses to new circumstances are adaptive and may depend on real or perceive risks of infection. Current disease prevalence, may become a marker or a tipping point, that when crossed, triggers individual or policy decisions. Whether or not individuals follow public health officials’ recommendations may be a function of individuals’ economic/educational status. The determination to make a drastic decision may be weakened or reinforced by each individual’s networks of friends. Responses are altered by the opinions of work-related connections. Personal needs play a role, and they include the need to use public transportation or the desire to attend a social event.

The landscape where behavioral decisions take place is not fixed. On the contrary, it may be altered from within (individual decisions) or by the use of preventive or active public health policy decisions or recommendations, some “obvious” like vaccination or quarantine, others drastic like mandated social distancing, as in the 2009-10 flu pandemic in Mexico [14]. Now, whether the dynamical changes experienced by the socio-epidemiological landscape are slow or fast, will depend on many factors. It is within this, often altered, complex adaptive dynamical system, that individual decisions such as individually-driven social distancing, the use of face masks, frequent hand washing, increased condom use, and or the routine use of non-pharmaceutical interventions (NPIs), takes place [39]. It has been theoretically documented [11, 29, 12, 16, 23, 37] that massive behavioral changes can impact the patterns of infection spread, possibly playing a critical role in efforts to prevent or ameliorate disease transmission. Today policies are implemented regardless of our knowledge of who are the “drivers” responsible for inducing behavioral changes. Modeling frameworks exist that allow for the systematic exploration of possible scenarios. Through a systematic exploratory analyses of appropriately selected classes of scenarios, it is possible to identify possibly effective (model-evaluated) public health policies. The use of highly detailed models including individual-based models has some advantages since they can incorporate individuals’ awareness of risk based on available of local information. The evaluations carried out on the effectiveness of changes, at the individual level, can be used to assess, for example, the impact of NPIs in reducing disease prevalence.

Models that couple disease dynamics and awareness to levels of infection risk have been proposed. These models have been used to explore the impact of behavioral changes on the spread of infection. In the review paper of Wang et al. [39] classify models as Rule-Based Models (RBMs), those where individuals make their decision about changing behavior independently of others, and Economic-Epidemiology models (EE models), that is, models where individuals change their behavior in order to maximize their own utility function (what they value) subject to available resources. The EE models account for the responses that individuals take in response to infection risks on disease prevalence at the population level. The modeling and results reported in this manuscript are more closely related to those used in RBMs.

Ruled-Based Models [16, 7, 34, 29, 18, 23, 31, 32, 33] vary from compartmental ODE models [16, 23, 37] to individual-based network models [41, 11, 13, 22]. These models have been used to study the dynamics of highly diverse diseases including, for example, influenza and HIV [31, 10], or in the study of generic infections [29, 18, 23, 33, 32]. Compartmental models (often using a phenomenological approach) categories designed to capture levels of awareness of infection. Such approach that can be used to incorporate ‘awareness’ in network models, is the objective of this manuscript. Some models assume that “awareness” spreads along with the invading disease, that is, through identical contact networks. Here, it is assumed that the disease and information spread over the same social network (a drastic simplification). The possibility that awareness and responses to the presence of a new infection among a subset of the population at risk, may significantly alter regular temporal patterns of disease prevalence (lower highs) have been studied. Studies have also shown that epidemic thresholds can be altered [34, 33] [29] in response to the effectiveness of NPI. Here, the focus is on the role of policy decisions/recommendations in altering disease dynamics, possibly the final epidemic size, within a model where awareness (generated by official actions) spreads among those susceptible to infection and their ‘friends’.

In this paper, we explore the impact of various, prevalence-dependent pre-selected thresholds (decisions made by public health officials) as triggers of possibly temporarily behavioral change. The time to a triggering event is assumed to depend on the prevalence of infection– a decision taken by health authorities (declaring some level of health emergency). We carry out simulations to explore the impact of variations on triggering prevalence-driven levels on the final epidemic size of a non-fatal infection under three distinct fixed artificial social structures modeled as, Erdos-Raynei, Small-world, and Scale-free networks.

2 Modeling Framework

The model captures changing disease transmission dynamics by incorporating, after a health emergency has been declared. The spread of risk information, that play by assumption over the same social network. Risk information is transmitted at different rates depending on the economic/education levels of individuals. We focus on single outbreaks, within a Susceptible–Infected–Recovered (SIR) framework. The dissemination of awareness (risk information) among the susceptible population, a process triggered by the declaration of a health emergency by the public health authorities, generates temporary changes on the knowledge and understanding of risk of infection, thus altering their activity and effectiveness as communicators of risk. Becoming aware of the risk of infection and inciting our ability to communicate risk, takes place at a rate modeled as a function of the economic/educational level of individuals in the network. Economic/educational levels are preassigned qualities to individuals in the network from pre-selected distributions. A susceptible individual may be a member of three sub-classes of susceptible: Unaware, $S^{u}$, Aware, $S^{a}$ and Indifferent, $S^{i}$. Unaware individuals that become aware of the risk of infection may (depending on their economic/educational level), change their state to Aware, that is, be ready to convey information, with different levels of enthusiasm, on the risks of the infection to neighbors in the network. Aware, is assumed to be a temporary state and so a transition to the Indifferent state is included. The triggers that determine or drive awareness state are a function of the disease prevalence level and not the time since start of the outbreak. The case when delays in transmitting real-time prevalence information to public health officials takes place, the most likely scenario is not considered. However, we know that it brings additional consequences (see [38]).

The activation process of public awareness takes place within a time framework determined by the actual (reported in general) disease prevalence and so, it is a function of the epidemic outbreak growth rate. We let $P^{*}$ denote the pre-selected level of prevalence identified by public health officials as the triggering threshold after which, a public awareness process campaign is initiated via mass media, word of mouth, social media and the like. The timing is selected when, according to public health officials. A disease outbreak may be about to reach unacceptable levels, with $P^{*}$, being the pre-selected tolerance level before a health emergency is declared. The $P^{*}$ determines the time $t^{*}$ when the information campaign starts (see Figure 1). The $t^{*}=t^{*}(P^{*})$ identifies the public health system threshold that forces public health action aimed at reducing outbreak consequences. The actual changes in behavior are determined by an stochastic process, and so, the shift in behaviors, is a function of “chance”.

Before $t^{*}$, the dynamics are governed only by the disease transmission process. At the time $t^{*}$, the model activates information process that models to behavioral change that reduce disease transmission. Therefore, $t^{*}$ denotes the time when prevalence reaches $P^{*}$ for the first time, Figure(1).

Sub-section 2.1 introduces the details of our network model, starting with the network structure and node attributes that are used to represent individuals in the network.

2.1 Network Structure

The network is denoted by $\textbf{G}=(\textbf{V},\textbf{E})$ which includes a set $N=|\textbf{V}|$ of nodes representing individuals $\textbf{V}=\{\textbf{i}|\textbf{i}=1,2,3,\dots N\}$ together with a binary adjacency relation defined by the set of edges $\textbf{E}=\{\textbf{ij}|\textbf{i,j}=1,2,3,\dots N\}$, where ij denotes the edge between individual (node) i and individual (node) j. We make use of three network structures, namely, Erdős-Rényi , Small-world, and Scale-free networks -defined in the Appendix- as models of our social landscape, the place where infection and awareness spread and behavioral change take place.

After generating the network [8, 26, 3], we assign some attributes to the nodes (individuals) as follows: each node i is associated with a random variable $x_{i}\in\left[0,1\right]$, generated from a beta-distribution, denoting the level of education of individual i, values closer to [math] corresponding to higher levels of education while those close to $1$ indicate limited education, an individual i is assigned a random number from a beta distribution with shape parameters $\alpha$ and $\beta$.

The underlying network G is weighted, where the weight $0<c_{ij}\leq 1$ for edge ij is the probability, of physical contact between neighbors i and j on any specific day. We assume that the values $c_{ij}$ are randomly generated from a uniform distribution $\mathcal{U}[0.5,1]$. For example if $c_{\textbf{ij}}=\frac{5}{7}$, that means that these two neighbors meet with probability $\frac{5}{7}$ and do not meet with probability $\frac{2}{7}$ on a specific day. We also assume that there are no birth or deaths in the population of nodes or edges, during the epidemic outbreak.

2.2 Dynamic of awareness among susceptible individuals

At any given time $t$, a person k is classified as susceptible, $S_{\textbf{k}}(t)$, infected, $I_{\textbf{k}}(t)$, or recovered, $R_{\textbf{k}}(t)$, to the disease under consideration. Disease risk awareness among susceptible individuals spreads after public health officials make the deliberate decision to stress vigorously the risk of infection, as it may be the case under an influenza epidemic or other health emergency. The model assumes that the decision to promote risk prevention starts once a pre-determined prevalence level $P^{*}$, at time $t^{*}$, has been reached. Health authorities start to spread information about a pathogen’s outbreak, information that gets communicated by those individuals, actively aware, in the network, with different degrees of efficiency. The spread of information is a function of chance (probability model) and economical/educational level. The state of active information-communicator is assumed to be temporary, with individuals moving into a“passive” state, indifferent, after a short period of time. Before time $t^{*}$ every susceptible individual is unaware of the information campaigns but for times $s>t^{*}$, members of $S_{\textbf{{k}}}(s)$ may change status from unaware ($S_{\textbf{k}}^{u}(s)$) to aware ($S_{\textbf{k}}^{a}(s)$), then moving, with some probability, into the indifferent class ($S_{\textbf{k}}^{i}(s)$).

We define $\widetilde{\lambda_{jk}}$ as the probability that $S_{\textbf{k}}^{a}(s)$ informs $S_{\textbf{j}}^{u}(s)$ at time $s$, and causes $S_{\textbf{j}}^{u}(s)\overset{\widetilde{\lambda_{jk}}}{\rightarrow}S_{\textbf{j}}^{a}(s+1)$. Thorough the probability $\widetilde{\gamma_{k}}$ -generated from an exponential distribution with parameter $\widetilde{\gamma}$, the average awareness period is therefore $1/\widetilde{\gamma}$ days. The members of $S_{\textbf{k}}^{a}(s)$ leave the awareness class into indifferent class, that is, $S_{\textbf{k}}^{a}(s)\overset{\widetilde{\gamma_{k}}}{\rightarrow}S_{\textbf{k}}^{i}(s+1)$, with probability $\widetilde{\gamma}_{k}$, after spending, on the average $1/\widetilde{\gamma}$ days in the aware class; the mean of waiting time distribution for the class $S_{\textbf{k}}^{a}(s)$.

2.2.1 Awareness transmission probability

$\widetilde{\lambda}_{jk}(s)$ denotes the probability that an aware person $S_{\textbf{k}}^{a}(s)$ will inform an unaware neighbor $S_{\textbf{j}}^{u}(s)$ at time $s$ via their contact: $\widetilde{\lambda}_{jk}(s)=c_{jk}\widetilde{\beta}$. The value $c_{jk}$ denotes the probability of contact between neighbors j and k, with $\widetilde{\beta}$ denoting the average probability that risk information will pass from an aware individual to an unaware neighbor. $\widetilde{\beta}$ is the average of a first increasing and then decaying and waning function $\widetilde{\beta(s,s^{*})}$ over the awareness period, that is, $\widetilde{\beta}=\widetilde{\gamma}\int_{s^{*}}^{s*+\frac{1}{\widetilde{\gamma}}}\widetilde{\beta(s,s^{*})}ds$, where $s^{*}\geq t^{*}$ is the first day that a person becomes aware. $\widetilde{\beta(s,s^{*})}$ is a function of $s^{*}$, because before this time the person was unaware and therefore, incapable of spreading risk information. Figure (2) helps visualize the probability function $\widetilde{\beta(s,t^{*})}$ for $t^{*}=10$, a function acting on initially aware individual, that is, individuals who became aware on the first day of awareness spread $t^{*}$. After $10$ days, information enters the system, people learn that an epidemic is taking place and the risk and severity of infection. The higher the value of $P^{*}$ the less effective the campaign in reducing the impact of information on the outbreak.

3 Simulations

In this Section, we perform our analysis and simulations on three different networks: Erdős-Rényi random network $\textbf{G}_{E}$, Small-world network $\textbf{G}_{W}$, and Scale-free network $\textbf{G}_{S}$. These networks are defined in the appendix. To generate the networks, we use Networkx111We use the NetworkX124 https://networkx.github.io/ open software platform to generate and analyze the network. and Erdős-Rényi -graph algorithm to generate $\textbf{G}_{E}$ with probability of connection $p=0.001$. Using the same package, the Watts-Strogatz algorithm is used to generate Small-world network $\textbf{G}_{W}$ with probability of rewiring $p=0.5$. We also use Barabasi–Albert’s preferential attachment model in Networkx to generate a Scale-free network $\textbf{G}_{S}$ with degree distribution $P(k)\sim k^{-2.11}$ [2]. The size of all networks is 10000 and the average number of neighbors for all networks is equal to 10. However, the networks differ in terms of degree distributions and other properties due to their different structures. Simulations start by infecting the most connected node (index case) to reach a balanced initial condition [1]. They are came out with the model baseline parameters in Table (1), unless stated otherwise.

3.1 Time series of infection

To determine the effectiveness of awareness spread, we compare the prevalence over time in the absence and presence of awareness. Figure (4) shows the results for network structures $\textbf{G}_{E}$, $\textbf{G}_{W}$, and $\textbf{G}_{S}$. We observe that all networks, with or without awareness, support an outbreak for the chosen parameters (that is, the probability of extinction is low or the results are conditioned on non-extinction). Due to a lack of epidemic threshold for Scale-free network [25, 4, 21], we always observe an outbreak with severity that depends on the initial infected index. For the other network structures, we observe an epidemic threshold depends on the network topology; specifically, on the average number of neighbors and the levels of heterogeneity in the number of neighbors (mean and variance of degree distribution) [19]. For example, for the case of Small-world networks Moore et al. [24] derived an analytic expression for the percolation threshold $p_{c}$, above which there will be an outbreak. For the case when the network is homogeneous (Erdős-Rényi network), the epidemic threshold is proportional to the average number of neighbors (average degree) [27]. Hence, for our simulation for both networks $\textbf{G}_{E}$ and $\textbf{G}_{W}$ we approximated the basic reproduction number using epidemic take-offs, when we had an outbreak, that is, whenever $\mathcal{R}_{0}>1$ ( probability of extinction for the parameters used seemed to be negligible).

The spread of infection on Scale-free network $\textbf{G}_{S}$ is faster than for the other two network topologies. Moreno et al. [25] , modeled infection with immunity on Scale-free and Small-world networks, they observed that the spread of infection on Scale-free complex networks is faster than that of Small-world networks due to a lack of epidemic threshold for network $\textbf{G}_{S}$; large connectivity fluctuations (heterogeneity in degree) on this network causes stronger outbreak incidence [25]. Our result is also related to [4], in which Chowell et al. tested the severity of an outbreak using an SIR model over a family of Small-world networks and an Scale-free network found out that the worst case scenarios (highest infection rate) are observed in the most heterogeneous network, namely, Scale-free networks.

The speed for infection has an impact on efficiency of awareness spread: for Scale-free network $\textbf{G}_{S}$ awareness reduces the peak of infection by roughly $6\%$, but for $\textbf{G}_{E}$ and $\textbf{G}_{W}$ this reduction is around $19\%$. Since in $\textbf{G}_{S}$ network infection spreads faster, hubs (individuals with many neighbors) get infected faster and loosing the chance of transmitting or receiving awareness. Wu et al. [41] showed a similar result: global awareness (behavior change because of higher prevalence in the population) on Scale-free networks cannot be as effective as local awareness in reducing infection.

3.2 Initiation of awareness

The presence and spread of awareness are coupled to the presence of infection, with awareness spreading, by design, after the infection reaches the level $P^{*}$. Therefore, there is a time-lag between starting time of infection and awareness spread. What should the value of $P^{*}$ be, if the goal is to reduce the final epidemic size? When the infection disperses fast- high basic reproduction number- people must be informed immediately, or the policy will have no effect. As noted, there is no threshold condition for $\textbf{G}_{S}$ network. We did not derive expressions for epidemic thresholds. Instead, we used epidemic take-offs to estimate the basic reproduction number as $\mathcal{R}_{0}=1.8$ and $2.01$, numbers corresponding to $\textbf{G}_{E}$ and $\textbf{G}_{S}$, respectively. Numbers that are not too high- such as that of the $\mathcal{R}_{0}$ for measles [17]- we are giving awareness a chance to be effective, since an epidemic with low $\mathcal{R}_{0}$ will reach the value $P^{*}$ much slower.

Proceed to quantify the sensitivity of the final size of the outbreak as a function of the prevalence threshold $P^{*}$ over the network structures GE, GW, and GS, Figure (5), for the simulated GE network that representing homogeneous random mixing and for the simulated GW network, we identify an optimal point for $P^{*}$ (the point that final size is minimized), that means awareness spread is most effective in reducing final size for these values of $P^{*}$. For larger values of $P^{*}$, the awareness policy is useless while for smaller values of $P^{*}$ is less effective, see subfigures (5(a)), and (5(b)).

For the simulated networks $\textbf{G}_{S}$ - in which network structure follows a power-law distribution - we observed that the sooner awareness gets started, the smaller the total number of infected individuals will be, see subfigures (5(c)) (see also [25, 4]).

To investigate the impact of awareness on the previous result and on network $\textbf{G}_{E}$, we vary the value of the awareness basic reproduction number $\widetilde{\mathcal{R}_{0}}$ via changes in the average time of behavior change $1/\tilde{\gamma}$. Figure (6) shows that the value of the optimal prevalence threshold $P^{*}$ increases as the $\widetilde{\mathcal{R}_{0}}$ increases. Accelerating awareness spread- increasing $\widetilde{\mathcal{R}_{0}}$- does not seem to change final awareness size trend (Figure (5)). Nevertheless, it makes the prevalence threshold $P^{*}$ more effective at reducing infection final size.

4 Discussion

We used an agent-based network model to model the spread of infection and risk information within a population of individuals in different epidemiological states, education levels and information states. In our model infection spreads within an SIR framework while awareness disseminates only among susceptible individuals. Specifically, when a fraction of the infected people reaches a pre-determined threshold, public health officials start a massive information campaign on the risk and severity of infection. Information is transmitted by those not infected only. Awareness dynamics are triggered by public health authorities with the objective of conveying information about risk and severity of infection via aware susceptible individuals to unaware susceptible neighbors in the network. The level of awareness achieved is tied in to the the educational level of each individual and the disease prevalence. Aware susceptible individuals can become indifferent, stop propagating information on the the risk of infection, after sometime.

Simulations of the spread of infection or awareness are carried out on different network topologies. The patterns vary even when preserving some network properties such as the mean degree [4, 35, 25]. We simulated the spread of awareness and disease over three different network structures, namely Erdős-Rényi , Small world and Scale-free networks.

The results of the infection-awareness model show that Scale-free networks support a fastest rate of infection with the impact of awareness having the lowest effect on transmission. Small-world and Erdős-Rényi networks became less effective at disease transmission, and the role of information via aware individuals had a stronger impact in reducing transmission, Figure (4). These result are somewhat similar to those obtained by Chowell et al. [4] and Moreno et al. [25] as they studied the rate of infection rate of growth on an SIR model over Small-world and Scale-free networks, in the absence of awareness. These researchers observed that the spread of infection on Scale-free complex networks is faster than in a Small-world one. Pastor et al. [28] and Dezső et al. [6] focused on identifying the best strategy that can be used to control (reduce) an epidemic peak on a scale-free network. They concluded that to in order to control infection public health workers needed to immunize the hubs at an early stage. Our results are somewhat related, that is, we find that in order to control an infection, public health workers need to start awareness spread soon, that is, before hubs play the role of super-spreaders.

Infection final size is a function of network topology and prevalence threshold, $P^{*}$, and the basic reproduction number $\mathcal{R}_{0}$. We observed that for Erdős-Rényi and Small-world networks under small value of $\mathcal{R}_{0}$, an optimal $P^{*}$ can be found for which the infection final size is minimized. This optimal point is the value that provides the best outcome under the prevalence $P^{*}$ at the time, $t^{*}$, the time when temporary protection was promoted by aware susceptible individuals, Figures (5(a)) and (5(b)). For Scale-free network, for which there is no epidemic threshold, a monotonic increasing function final size is observed, a function of the prevalence threshold $P^{*}$, there is no optimal point $P^{*}$ that minimizes infection final size, Figure (5(c)).

Infection network and awareness network: One of the big assumptions in our model is that the infection and awareness spread over the same network, ignoring the fact that individuals may have multiple sources of information outside their physical contact network. In fact, we know that information does flow through virtual neighbors such as facebook friends. How should we incorporate the role of physical network and virtual information? We are exploring possibilities.

Static network: Again, since we are focusing on single epidemic outbreak, we have assumed a short time frame in our simulation. We also assumed that the context network is static. The need to extend the model to a dynamic network where individuals are allowed to change neighbors is important in the study of disease spread over the longer time scale. The impact of awareness in such extended model would be prone to change as well.

For sexually transmitted infections on Scale-free network [20], waiting for enough people to become infected people before starting a campaign about risk of infection is not reasonable. For measles in children with basic reproduction numbers very high [17], waiting to reach pre-selected prevalence thresholds may also not be great idea. Our model results, highlight the importance of campaigns that warns a population about the risk and severity of infection for diseases that do not spread too fast, possibly including some flu infections or the severe acute respiratory syndrome (SARS). For disease with a basic reproduction number that it is is not high [17], and assuming that these infections spread over a random homogeneous networks such as Erdős-Rényi , it may not be unreasonable the existence of an optimal prevalence threshold to start and starting an awareness dynamic campaign.

Appendix

In this appendix, first we define some models for the network structures we used for our simulation, then we present the pseudocode of our model.

Definition 4.1

Erdős-Rényi random network Suppose we have $n$ disconnected nodes, and then we make edge between two arbitrary nodes with probability $p$ independent from every other edge. The graph constructed via this model is called Erdős-Rényi random network $G(n,p)$. The parameter $p$ in the model is a weighting function $p\in\left[0,1\right]$, where for $p$ closer to one the generated graph is more likely to construct graphs with more edges [8].

Definition 4.2

Small-world network Suppose we have $n$ nodes over a ring, where each node in the ring is connected to its $k$ nearest neighbors in the ring. Then with probability $p$ and independent from any other pair of edges we select two edges ${i_{1}j_{1}}$ and ${i_{2}j_{2}}$ such that $j_{k}$ is from the k nearest neighbor of $i_{k}$ for $k=1,2$, and rewire them. Then the constructed network is called Watts-Strogatz Smal-world network [40] $G(n,k,p)$.

Definition 4.3

*Preferential attachment Scale-free network Suppose we want to generate a network of $n$ nodes. We start with one node and follow the following procedure $n$ times: the probability of attaching a new node to the existing ones is proportional to their current degree:

$P(\textnormal{connect a new node i to existing node j})=\frac{deg(j)}{\sum_{k}(deg(k))}$, where the $deg(k)$ is the number of neighbors for node $k$, then with this probability we make new edge $ij$. The constructed network via this model is called Barabasi–Albert Scale-free network [2]. The degree distribution of this network follows a power-law distribution.*

Now, we define assumptions and algorithm of our model. Person $i$ at any given day $t$ has one of the states: $S^{u},S^{a},S^{i},I^{q},I^{f}$, or $R$. The fraction of infected people at dat $t$ is $P(t)$, and $P^{*}$ is a positive fraction $\in[0,1]$. We define $t^{*}=min_{t}\{P(t)=P^{*}\}$ which is the first time that $P(t)$ reaches $P^{*}$. We assume Education level for each person $k$, $x_{k}$ does not change by time. Also prior to time $t^{*}$ every susceptible person is at state $S^{u}$, and finally the time scale for each update is day, in which each person can have at most one update. The following algorithms show the details.