key: cord-0861906-u8ujwn8i authors: Guo, Haili; Yin, Qian; Xia, Chengyi; Dehmer, Matthias title: Impact of information diffusion on epidemic spreading in partially mapping two-layered time-varying networks date: 2021-08-20 journal: Nonlinear Dyn DOI: 10.1007/s11071-021-06784-7 sha: 6e7ece337bb0c3a2985a469876f6303b3bd7d08c doc_id: 861906 cord_uid: u8ujwn8i We propose a new epidemic model considering the partial mapping relationship in a two-layered time-varying network, which aims to study the influence of information diffusion on epidemic spreading. In the model, one layer represents the epidemic-related information diffusion in the social networks, while the other layer denotes the epidemic spreading in physical networks. In addition, there just exist mapping relationships between partial pairs of nodes in the two-layered network, which characterizes the interaction between information diffusion and epidemic spreading. Meanwhile, the information and epidemics can only spread in their own layers. Afterwards, starting from the microscopic Markov chain (MMC) method, we can establish the dynamic equation of epidemic spreading and then analytically deduce its epidemic threshold, which demonstrates that the ratio of correspondence between two layers has a significant effect on the epidemic threshold of the proposed model. Finally, it is found that MMC method can well match with Monte Carlo (MC) simulations, and the relevant results can be helpful to understand the epidemic spreading properties in depth. Emerging epidemics can create the heavy harm to human health [1, 2] , such as SARS [3] , Ebola [4] and COVID-19 [5] . How to prevent large-scale outbreaks of infectious diseases has always been a challenging topic in the area of disease intervention and public health [6] . Therefore, it is particularly important to use effective methods to prevent and control the diffusion of infectious epidemics. For that matter, experts in various fields begin to focus on the spread of infectious diseases by establishing infectious disease models [7] [8] [9] , which try to characterize the spreading law and present some effective measures to fight against largescale pandemics. Some classical propagation dynamical models were proposed [10] , such as susceptible-infective (SI), susceptible-infective-susceptible (SIS) and susceptible-infective-recovery (SIR) models. An important indicator of epidemic transmission is the critical threshold for the epidemic spreading, which denotes the min-imum transmission rate to induce an extensive outbreak [11, 12] . Among them, Parshani et al. [13] presented an analytic expression of threshold for (SIS) model over a random network. To be notable, multiple factors may affect the threshold of epidemic spreading process within the population. As an example, Wu et al. [14] investigated the effect of global, local and contact awareness on the propagation threshold of the SIS model in scale-free network. They proved that the global awareness decreased the epidemic prevalence, while the other two types of awareness increased the epidemic threshold. After more than 20 years of development, complex network theory has made great achievements [15] [16] [17] [18] [19] , which provides a new framework to favor the study of epidemic transmission. In the epidemic dynamics, the multilayer network structure offers a novel topology foundation to model the epidemic spreading since the spread of epidemic is often coupled with the diffusion of epidemic-related information, such as rumors and knowledges [20] [21] [22] , risk perception [23, 24] and fear [25, 26] , which is usually spread by means of online social networking platforms, (e.g., WeChat, QQ and Twitter) [27, 28] . Any individual will take preventive steps to prevent themselves from being infected after they know the relevant information of the epidemics, for instance, wearing masks, reducing outdoor activities and paying attention to the personal hygiene. Infected individuals also spread information on social networks. Thus, epidemic-related information plays a significant role in the spread of epidemics. For example, individuals who are aware of the epidemics may reduce the probability of being infected, and infected individuals who spread information facilitate the diffusion of information, and thereby, the probability of susceptible individuals being infected is reduced. Some studies have shown that epidemic information can greatly affect epidemic outbreaks, such as reducing the size of outbreaks and raising the epidemic threshold [29] [30] [31] [32] [33] [34] . In particular, Granell et al. [35] put forward a novel UAU-SIS (unaware/aware/unawaresusceptible/infected/susceptible) model to investigate the coupling characteristics between information diffusion and spread of epidemic dynamics, and they found that epidemic threshold is correlated with the information transmission property on the diffusion layer as well as the network topology of epidemic propagation layer. Also, there is a critical value of information transmission rate; starting from this point, the threshold increases and the outbreak size is reduced. Furthermore, they improved their model [36] to take into account the impact of mass media on epidemic transmission, and the tipping point that had appeared in previous studies disappeared under the influence of mass media. Xia et al. [37] further explored the characteristics of the coupled transmission of information and SIR-type epidemics under the influence of mass media, and they found that the increase of the transmission rate of mass media could significantly reduce the size of epidemic outbreak and increase the epidemic threshold. In addition, Wang et al. [38] discussed the epidemic transmission behaviors under positive and negative prevention information, and they proved that the coexistence of positive and negative prevention information can substantially modify the epidemic outbreak size and the critical threshold giving rise to the largescale diffusion of epidemics. As a further step, they improved the model to explore the influence of mass media on the spread of the epidemic, indicating that mass media has an impact on the size and threshold of the epidemics [39] . However, in the related studies mentioned above, it is usually assumed that each node in one layer precisely corresponds to one node in the other layer; that is, there exists one-to-one mapping relationship between two layers for any node pair. It is also hypothesized that any individual can effectively avoid being infected once he obtains the epidemic-related information; at the same time, infected individuals will carry out the information transmission in the information layer. Nevertheless, during the real outbreaks, infected agents may not transmit epidemic messages in time due to various reasons. As an example, AIDS patients are usually not willing to actively report their sexual partners of their HIV infection [40] . Meanwhile, there may exist some individuals who have the information about epidemics, while they are still reluctant to adopt the effective prevention measures, which renders the infectious diseases to be further spread. Taking an example, a considerable number of students in a university in South Africa have knowledge about HIV, but most of them have not adopted the preventive measures to protect themselves [41] . Henceforth, it is necessary to investigate the impact of the mapping relationship between partial node pairs on epidemic transmission, in which epidemics and epidemic-related information are considered simultaneously. It is worth mentioning that we have firstly discussed the impact of partial mapping of node pairs between two-layered networks based on the static multiplex network topology [42] , and it is indicated that ratio of correspondence between two layers of networks has greatly influenced the critical threshold of proposed epidemic model. However, information diffusion and epidemic transmission usually occur in time-varying networks in reality, and then, the static networks often ignore the time-varying characteristics of individual connection patterns. In order to more reasonably describe the transmission dynamics of epidemics in the real world, time-varying networks need to be considered in the modeling of epidemic spread [43] [44] [45] [46] . For this reason, we try to utilize the activity-driven network to describe the network structure of information diffusion and epidemic propagation. Furthermore, beginning from the framework of two-layered networks, we put forward one coupled epidemic model to delineate the interaction between epidemics and information associated with epidemics, where one layer stands for the epidemic-related information diffusion network, while another one is used as the underlying network for the epidemic transmission. However, different from previous works [35] [36] [37] [38] [39] , we only consider that there exists the mapping relationship for partial node pairs between two layers. On the contrary, they will propagate independently if the corresponding nodes have no such mapping relationships. Particularly, we take use of the microscopic Markov chain (MMC) method to analytically acquire the critical threshold of epidemics for the proposed model. Meanwhile, large quantities of MC simulations are used to demonstrate the efficacy of MMC approach. The remaining sections of this paper are organized as follows. First of all, the proposed epidemic model is introduced in detail in Sect. 2. Secondly, in Sect. 3, MMC method is used to analyze the model and derive the analytical expression of the epidemic threshold, and then, extensive simulations are conducted to explore the influence of correspondence rate on the epidemic spreading dynamics. Finally, we conclude the paper and perform some outlooks in Sect. 4. In this work, starting from an activity-driven network model, time-varying network is used to characterize the coupling transmission and interaction of informa-tion diffusion and epidemics transmission. In order to make our introduction to the proposed epidemic model clearer, we list the definitions of key quantities or parameters of the epidemic model in Table 1 . As shown in Fig. 1 , Layer (V) denotes the diffusion of epidemic-related information on virtual online social networking platforms (e.g., WeChat and Weibo), and Layer (R) describes the epidemic transmission in real social networks (say, friends and family) [47] . There is a mapping relationship between partial node pairs, which indicates that only some individuals will change their behavior once they are infected or receive the epidemic-related information. At each layer, the network is generated by the actively driven model [48] , which includes N nodes. Each node has a fixed activity level, denoting the probability of this node being activated to build links with other ones at each time step. On the upper layer (i.e., Layer V), the activity level for every node is expressed as a i = η V x i ; while for the lower layer (that is, Layer R) the activity level of each node is assumed to be b i = η R y i , where x i and y i are a set of sequences following the power-law dis- and meanwhile, they satisfy the conditions ε ≤ x i ≤ 1 and ε ≤ y i ≤ 1, where ε represents the minimum value of x i or y i , we set ε = 5 × 10 −4 , and γ V and γ R represent the activity index of Layer (V) and Layer (R), respectively. The smaller activity index γ V or γ R means the more heterogeneous node activity distribution on Layer V or Layer R. Figure 1 depicts one basic step of time-varying network construction on Layer (V) or Layer (R), which are divided into the following three steps. 1) At the beginning of time step t, both Layer (V) and Layer (R) are composed of N isolated nodes, and there just exist the mapping relationships between partial node pairs for two-layered networks; 2) the node at layer V (R) becomes the active state with probability a i (b i ) and randomly selects other m V (m R ) nodes to establish the links; and 3) at the next time step t + 1, all links on Layer (V) or Layer (R) are removed, the corresponding layer is restored to N separated nodes, and repeat the previous two steps until the propagation of the two layers reaches a stable state. At each iterating time step, the average degree of Layer (V) and Layer (R) is Next, the epidemics and information spreading model is defined in detail. The process of information x i Node i's activity potential on Layer (V) Rescaling factor on Layer (R) Distribution of activity potential for nodes on Layer (V) Distribution of activity potential for nodes on Layer (R) γ V Exponent of activity level on Layer (V) γ R Exponent of actively level on Layer (R) ε The lower limit of activity potential distribution Correspondence rate between two layers Probability that an individual i on Layer (V) is not informed by his neighbors at time step t Probability that an aware individual i is not infected by his neighbors at time step t q U i (t) Probability that an unaware individual i is not infected by his neighbors at time step t transmission on Layer (V) is described by UAU-type dynamics, in which individuals may be at the unaware (U) or aware (A) state. Among them, an unaware agent will acquire the epidemic-related information with the probability λ from his neighbors, and an aware individual will forget the related information with probability δ. In addition, the classical SIS-type model is used to describe the epidemic transmission at Layer (R), where each individual may be susceptible (S) or infected (I). A susceptible individual is transformed into an infected one with probability β, and the infected individual returns to the susceptible state with probability μ. Dif-ferent from the classic SIS model, information transmission will affect the probability of being infected when the two-layered network is coupled. Accordingly, under the case of two-layered network and partial nodepair mapping, there exist 4 possible states for any node pair in the current model, which include US, UI, AS and AI. Furthermore, any one of 4 states can be transformed into other possible ones with a certain probability, and corresponding state transition can be shown in Fig. 2 . In this model, in order to characterize the partial mapping of node pairs between layers, we use the vector L = [l 1 , l 2 , . . . , l N ] to describe the mapping rela- Fig. 1 The UAU-SIS epidemic model with partial node pairs mapping in time-varying network. Layer (V) denotes the epidemic-related information diffusion network, and there are two kinds of node states including aware (A) and unaware (U); Layer (R) network represents the spread of epidemics, and there are susceptible (S) and infected (I) states in the Layer (R). A virtual connection between two layers means that partial node pairs hold the mapping relationship. Individuals in Layer (V) are represented by a circle, while individuals in Layer (R) are represented by a square. The size of the circle and the square indicates the level of individual activity. The larger the size, the higher the level of activity. Firstly, in panel a, two layers are, respectively, composed of N isolated nodes. As an example, node 6 is in the aware state on Layer (V) , and infected state on Layer (R). Secondly, in panel b, nodes 2, 7 and 10 in Layer (V) are active nodes, and they randomly select nodes at Layer (V) to establish two links (m V = 2). In addition, nodes 1, 6 and 9 in Layer (R) are active, and they randomly select nodes at Layer (R) to establish two links (m R = 2). At this step, node 5 in Layer (R) is infected by neighbor node of 6, and the corresponding node in Layer (V) becomes aware. At the same time, node 4 is also infected by neighbor node 6, while the state of node 4 in Layer (V) remains unaware since there is no mapping relationship between layers. In Layer (V), node 10 becomes aware by acquiring the information from its neighbor of node 6. At the next step, all links are deleted and each layer is still composed of N independent nodes, and then, the processes in panel b are repeated until the diffusion dynamics in the two layers tend to be stable Fig. 2 The probability transition tree for four states (US, AS, AI and UI). Nodes that are unaware of the information will not be transmitted by their neighbors with probability r i , and individuals that are aware of the information will forget the information with probability δ. The probability that the node i is unaware of the information will not be infected by its neighbors which is q U i , while the probability is q A i if the individual is aware of information. In addition, the infected node returns to the susceptible state with probability μ. l i denotes whether the node pairs on two layers of networks correspond to each other; when the nodes in two layers correspond to each other, l i = 1; on the contrary, there is no correspondence between the node pairs, and information diffusion and epidemic transmission are independent of each other and so l i = 0 tionship: If there is no mapping between a pair of nodes, l i = 0; otherwise, l i = 1. The rate of infection (β) will decay with an attenuation factor ω i when a susceptible individual i obtains the information related to epidemics, that is, where ω i can be defined in the following way, In the model, β U = β is used to denote the rate of infection for an unaware susceptible to be one infected neighbor. In addition, when l i = 1, we assume that the epidemic-related information will be immediately spread once a US-type individual is infected. However, the individual state on the Layer (V) will not be infected if there is no mapping between node pairs (i.e., l i = 0). At any time step t, P U S i (t), P U I i (t), P AS i (t) and P AI i (t) can be used to denote the probabilities for an individual i at the states of US, UI, AS and AI, respectively. For each node i, the normalized condition , respectively. Further-more, r i (t) denotes the probability that any unaware node i at Layer (V) will not propagate the information at time step t. Similarly, for the node in Layer (R), if the state of this node in Layer (V) is unaware, then the probability that the individual will not be infected by any infected neighbor at step t is q U i (t); Otherwise, the probability is q A i (t). Thus, we can get the following equations, where P A j = P AS j + P AI j , P I j = P U I j + P AI j . Starting from the state transition trees in Fig. 2 , the dynamical equations of four states for the propose model are obtained as Eq. (3). When the system arrives at the steady state, the following equations can be obtained: Epidemic outbreak threshold is an important indicator to characterize whether the disease outbreaks can happen in the field of infection disease modeling. In this section, the dynamic process of Eq. (3) is analyzed by MMC method, and we deduce the epidemic threshold of the proposed model by using Theorem 1. A = [a 1 , a 2 , . . . , a N ] and Proof When β is near the threshold, the probability that the node is infected is almost 0. Hence, let P I i = P AI i + P U I i = i 1, then Eq. (2) can be approximated as 123 Impact of information diffusion on epidemic spreading When β is near epidemic threshold, Eq. (3) at the steady state can be further simplified as follows, Since After that, the first and second equations in Eq. (6) can be further approximated by removing the higherorder terms O( i ) as Henceforth, we can further simplify the last term in Eq. (6) into the following equation, Since P U i + P A i = 1, we can reconstruct Eq. (8) as follows, Taking the average value of all nodes into account, we calculate the infection density to be ρ I = 1 N N j=1 P I j . Then, Eq. (9) can be written as follows, Multiply both sides of Eq. (9) by b i and take the average value of all nodes to obtain the following equation: where . Thus, we can rewrite Eq. (10) and (11) into the matrix form as follows, where The network system satisfies ρ I θ I b = − → 0 when the epidemic is greatly transmitted. When the epidemic breaks out in a networked system, β c is the minimum β value satisfying the formula Eq. (12) such that μ m R β c is equal to Λ max (H ), which is the maximum eigenvalue of matrix H and can be defined as As a consequence, we can characterize the epidemic threshold as It is worth noting that in Eq. (15) , the epidemic threshold is related to the diffusion of the information, the rate of epidemic recovery, the correspondence rate between the Layer (V) and Layer (R) and attenuation factor ω i . In addition, the level of node activity and connectivity in the Layer (V) and Layer (R) will also vary threshold according to Eq. (15). In this subsection, extensive Monte Carlo simulation experiments are carried out to validate the theoretical predictions in Theorem 1. For all experiments, the number of nodes at each layer is set to be N = 10000, and the vector L is used to describe whether there exists the mapping relationship between corresponding pairs. Once the vector L is determined, it will remain unchanged in subsequent experiments. Therefore, we define the mapping rate between two layers as Θ = ( i l i )/N × 100%, where l i ∈ L . In addition, ρ I represents the proportion of infected individuals at the stable state. In the iterative calculation of MMC, the proportion of I-state is written as Fig. 3 , the relative error is 4.1% when Θ = 40%; similarly, they are 3.7%, 3.5% and 4.5% when mapping rates are 60%, 80% and 100%, respectively. Therefore, MMC method can be accurately used to depict the epidemic spreading behaviors of the model described above. Secondly, in Fig. 4 we consider the effect of Θ on ρ I in the case of different combinations of β and λ . There is no doubt that the value of ρ I increases with the rise of β when λ ∈ [0, 1] is constant. However, this increasing tendency weakens when Θ is varied from 40% to 100%. This is because when the value of β is larger, the disease propagation may play a significant role so as to promote the information diffusion. In the case of the higher corresponding rate, more individuals have opportunities to take the preventive measures to avoid being infected once they receive the diseaserelated information. Meanwhile, larger λ leads to the reduction of ρ I , but the level of reduction is not obvious. For example, in the case of Θ = 40%, the value of ρ I changes from 0.5164 to 0.4124 when λ changes from 0 to 1, while it declines from 0.3886 to 0.015 when Θ = 100%. In addition, the effect on ρ I by increasing λ is inferior to change the value of Θ. For instance, the value of ρ I in Fig. 4a is around 0.4831, 0 .4562, 0.4349 and 0.4234 for different values of λ = 0.2, 0.4, 0.6 and 0.8 when we fix the infection rate β = 0.5. However, ρ I in Fig. 4b-d is 0.4141, 0.3412 and 0 .2688 when β = 0.5 and λ = 0.2. Therefore, we may obtain the conclusion that the increase in the correspondence rate has a greater impact on information dissemination than increasing λ during the epidemic outbreaks. In the proposed model, the following three parameters can affect the structure of the two-layered time- varying networks: activity distribution index γ , contact ability m and the correlation coefficient C of node activity between two layers. Specifically, the activity index γ represents the heterogeneity of individual activity distribution, and the smaller the value of γ , the stronger the heterogeneity of nodes activity distribution. Meanwhile, the contact capacity means the average number of links established by an active node, and the activity correlation between two layers describes how nodes at different activity levels in Layer (V) and Layer (R) are coupled to each other. Next, we will explore the impact of these three factors on epidemic transmission in a partially mapping two-layered time-varying networks. The influence of activity heterogeneity of the node at Layer (V) with different correspondence rate Θ is plotted in Fig. 5 . As can be seen from Fig. 5 , the value of Θ can affect the epidemic size and epidemic threshold. For example, ρ I in Fig.5a is 0.40, 0.31, 0.24 and 0.14 when Θ is varied from 40% to 100%, and the outbreak thresholds in Fig.5b are around 0.16, 0.18, 0.21 and 0.25 when γ V = 3.0. This may be because in the case of higher value of Θ, more people have chances to take preventive measures once they know the epidemic-related information. At the same time, they choose to spread the information after they are infected, thus effectively inhibiting the spread of the epidemics and modestly reducing the threshold of epidemic outbreak. The activity heterogeneity of Layer (V) has little effect on epidemics when Θ is lower. For example, when γ V = 2.0, 2.5, 3.0, 3.5, ρ I at the steady state is 0.3871, 0.3956, 0.4034, 0.4057, respectively, and epidemic threshold β c is 0.1658, 0.1618, 0.1600, 0.1569, respectively, when Θ = 40%, while ρ I is 0.1078, 0.1247, 0.1478, 0.1579, and the threshold β c is 0.2925, 0.1618, 0.1600, 0.1569, respectively, when Θ = 100%. This may be because in the case of low correspondence rate, individuals who know the information and take effective measures become fewer. In this case, increasing node activity can affect and promote the further spread of information, but there are still fewer individuals who take preventive measures, and thus, the inhibition effect on epidemic transmission is not obvious. Moreover, with the activity heterogeneity of nodes in Layer (V) being strengthened (say, decreases of γ V ) , the epidemic size decreases and the threshold increases for different correspondence rates. Because the decrease of γ V leads to the increase of the average activity potential x of nodes in Layer (V), which further induces the mean degree of nodes in Layer (V) to be increased, thereby the epidemic-related information diffusion is promoted and the chance of individuals taking preventive measures is also enlarged to inhibit the spread of epidemics. Figure 6 depicts the effect of activity heterogeneity of the node at Layer (R) at different correspondence rates Θ. We find that, for different corresponding rates Θ, the change of the activity heterogeneity of nodes in Layer (R) has a substantial influence on the density of infected individuals at the steady state and the threshold of epidemic outbreaks. For instance, the steady-state values of ρ I are 0.32, 0.24, 0.15 and 0.06 in Fig. 6a , and the outbreak thresholds are 0.19, 0.22, 0.27 and 0.33 in Fig. 6b when Θ varies from 40% to 100% and γ R = 3.0, which may be because when the mapping rate is high, infected individuals are willing to actively spread this information, and at the same time, more individuals will know the epidemic-related information and take preventive measures, thus slowing down the spread of the epidemic. In addition, in the case of different correspondence rates Θ, changing the activity heterogeneity of nodes in Layer (R) has a certain influence on epidemic transmission. As an example, when γ R = 2.0, 2.5, 3.0, 3.5, the density of infected individuals is 0.3941, 0.3563, 0.3236, 0.2935, and threshold is 0.1615, 0.1784, 0.1917, 0.2052, respectively, when Θ = 40%. Likewise, the ρ I is 0.1189, 0.0851, 0.0588, 0.0365, and the threshold is 0.2802, 0.3090, 0.3335, 0.3556, respectively, when Θ = 100%. We can know that the activity heterogeneity of node in Layer (R) is enhanced, the average activity potential y of nodes in Layer (R) increases, and leading to an increase in the average degree of the network in Layer (R), it is conducive to epidemic transmission in Layer (R). Meanwhile, when the corresponding rate Θ gradually decreases, the influence on the trend of epidemic transmission does not weaken. Figure 7 describes the changes of ρ I and β c in the case of different correspondence rates Θ with the number of links m V established by the active node in Layer (V). The change of m V at different correspondence rates has a great difference as far as the density of infected individuals and the threshold are concerned. For instance, when m V changes from 1 to 8, ρ I decreases from 0.4439 to 0.3524 when Θ = 40%, while it declines from 0.2766 to 0.004 when Θ = 100%. Furthermore, β c increases from 0.1289 to 0.1800 and from 0.1291 to 0.4434, respectively. With m V continuing to increase, the change of ρ I and β c gradually decreases. The increase of m V is conducive the k V to increase, leading to the spread of epidemic information in Layer (V), and thus effectively inhibits the spread of epidemic. When the correspondence rate is higher, increasing the number of links created by an active node in Layer (V) can rapidly reduce the infection density and increase the threshold of epi- demic outbreaks, while, for the lower corresponding rate Θ, the change of m V has little effect on the threshold. As an example, when m V is varied from 1 to 8, the threshold is slightly increased from 0.1289 to 0.1903 at Θ = 40%, while it rises from 0.1291 to 0.6587 for Θ = 100%. Therefore, the increase of correspondence Next, we explore the impact of Θ on the epidemic spreading when m R in Layer (R) is varied, and related results are shown in Fig. 8 . Similarly, m R is varied from 1 to 8, and ρ I is increased from 0 to 0.5354 when Θ = 40%, but it increases from 0 to 0.3058 if Θ = 100%. Meanwhile, β c is reduced from 0.6292 to 0.0786 and from 0.9361 to 0.1173, respectively. As m R increases, the change of ρ I and β c gradually decreases and tends to be stable. The increase of m R leads to k R increase, which is conducive to the spread of epidemic in Layer (R). Moreover, the corresponding rate Θ is smaller, and the influence of infection density is greater when m R is increased. For example, m R changes from 1 to 20, and the infection density changes from 0 to 0.6481 when Θ = 40%, while ρ I changes from 0 to 0.4353 when Θ = 100%. Since the correspondence rate Θ is lower, fewer infected individuals in Layer (R) have chances to spread information, which renders the fewer individuals in Layer (V) to know the information. Thus, fewer individuals will take effective prevention measures, which leads to the further transmission of the epidemic. When m R increases up to a certain value, the change of m R has little impact on the threshold, and the difference between different correspondence rates Θ becomes smaller. As an example, when m R = 20, the correspondence rate Θ changes from 40% to 100%, and the threshold values are 0.0315, 0.0354, 0.0403 and 0.0468, respectively. The higher the corresponding rate Θ, the stronger the inhibitory effect on the epidemic when m R changes. Finally, in order to explore the influence of the proportion of positive or negative correlation node pairs on epidemic size and outbreak threshold in different corresponding rate, we use C to describe the proportion of positive or negative correlation between node pairs. For example, C = 0.5 means that there are 5000 nodes in Layer (V) whose activity ability is positively mapping with nodes of higher activity ability in Layer (R). However, C = −0.5 represents that there are 5000 nodes with higher activity capacity in Layer (V) corresponding to nodes with lower activity capacity in Layer (R). And the remaining 5000 nodes are randomly mapping. As shown in Fig. 9 , the influence of correlation rate C on epidemic size and threshold is increased with the corresponding rate Θ being raised. For instance, the infection density ρ I and threshold β c change from 0.3928 to 0.3846 and from 0.1606 to 0.1695, respectively, when Θ = 40%, while for Θ = 100%, ρ I changes from 0.1379 to 0.0671 and β c changes from 0.2544 to 0.3307. Meanwhile, when the value of C is positive, the infection density ρ I decreases and the threshold β c increase as C increases. However, the changes of ρ I and β c are distinct from those foregoing phenomena when the value of C is negative. It is worth mentioning that the change of infection density ρ I and threshold β c is very obvious at the positive and negative boundary; nevertheless, they do not change much more when C is far from zero. It can be explained that the number of nodes with strong activity is limited, the mapping relationship between nodes with the relatively lower activity ability has little influence on the final results, and the positive correlation promotes the spread of information and inhibits the spread of epidemic. Here, the individuals with high activity ability both in Layer (V) and Layer (R) can enhance the transmission of epidemic-relative information and inhibit spread of the epidemic. As indicated in real life, the behaviors of key individual may determine the evolution of epidemics. In summary, based on the framework of two-layered time-varying networks, we present an improved epidemic model to study the interactions between epidemicrelated information and epidemics, where only partial node pairs have the coupling or mapping relationships on two-layered networks. In the proposed model, the UAU model is used to depict the epidemic-relative information diffusion in the first Layer (V), and the SIS model is used to describe epidemic transmission in the second Layer (R). In this model, only partial individuals are considered to have the mapping relationship between two layers, and so only some individuals would take preventive measures after knowing the epidemic-related information or respond to the state of another layer after being infected. We utilize MMC method to analyze and derive the epidemic threshold of the model, in which the threshold is not only related to the information transmission of Layer (V), but also closely correlated with the ratio of correspondence of node pairs between two layers. A great quantity of Monte Carlo simulations are also performed to validate the theoretical predictions. The current model indicates that the outbreak threshold and infection density of Layer (R) are related to the information diffusion dynamics of Layer (V), and high correspondence rates are beneficial to the suppression of epidemic transmission. Our results are also qualitatively similar to these cases in static network. Different from the static two-layered network propagation dynamics, we investigate how the unique activity properties of the time-varying networks of Layer (V) and Layer (R) and the coupling relationship between the two layers affect the epidemic transmission at dif- ferent corresponding rates. It is found that nodes in Layer (V) with higher heterogeneity, which is conducive to information diffusion, have an inhibitory effect on epidemic, and this phenomenon is more obvious when the correspondence rate is higher. The activity heterogeneity of node in Layer (R) tends to reduce the threshold and enhance the infection density, and the trend is consistent for different corresponding rates. In addition, we also find that the increase of the contact capacity of nodes in Layer (V) can reduce the infection density, and the influence of contact capacity on the spread of epidemic is significant when the correspondence rate is higher. However, the increase of the contact capacity of nodes in Layer (R) is conducive to the transmission of epidemics, but the corresponding rate has little effect on the threshold after the contact capacity reaches a certain value. Finally, the influence of correlation between two layers on epidemic transmission is explained, it is found that there is an abrupt transition between positive and negative correlation coefficient C, and the positive correlation is conducive to the suppression of epidemic transmission. In the future, we will further explore the role of individual heterogeneity in the awareness diffusion under the framework of multiple time-varying networks, which can be helpful to realistically characterize its spreading course. In addition, how to model the impact of mass media on the awareness diffusion will also become very necessary since the media information can confront the whole population. Meanwhile, the disease propagation and corresponding prevention measures may greatly change the population topology, and then, probing into the influence of epidemic spreading on the awareness diffusion will also make a difference from the current work. Emerging infectious diseases: threats to human health and global stability Infectious disease threats in the twenty-first century: strengthening the global response Isolation and characterization of a bat SARS-like coronavirus that uses the ACE2 receptor Ebola haemorrhagic fever Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and coronavirus disease-2019 (COVID-19): the epidemic and the challenges Changing patterns of infectious disease Infectious disease threats: a rebound to resilience The mathematics of infectious diseases Statistical physics of vaccination Infectious Diseases of Humans: Dynamics and Control Thresholds for epidemic spreading in networks Threshold behaviour and final outcome of an epidemic on a random network with household structure Epidemic threshold for the susceptible-infectious-susceptible model on random networks The impact of awareness on epidemic spreading in networks Epidemic processes in complex networks Epidemic outbreaks in complex heterogeneous networks Epidemic dynamics of influenza-like diseases spreading in complex networks An SIR model with infection delay and propagation vector in complex networks A novel epidemic model considering demographics and intercity commuting on complex dynamical networks Modeling the competitive diffusions of rumor and knowledge and the impacts on epidemic spreading The dynamics analysis of a rumor propagation model in online social networks The analysis of an seir rumor propagation model on heterogeneous network Effect of risk perception on epidemic spreading in temporal networks Understanding the community risk perceptions of the covid-19 outbreak in south Korea: infodemiology study Modeling individual fear factor with optimal control in a disease-dynamic system Epidemics and fear Close and ordinary social contacts: how important are they in promoting large-scale contagion Friend recommendation with content spread enhancement in social networks Epidemic dynamics on information-driven adaptive networks Dynamics analysis of epidemic and information spreading in overlay networks Analysis of epidemic outbreaks in two-layer networks with different structures for information spreading and disease diffusion Coupling dynamics of epidemic spreading and information diffusion on complex networks Strategy to suppress epidemic explosion in heterogeneous metapopulation networks Hybrid multiscale coarse-graining for dynamics on complex networks Dynamical interplay between awareness and epidemic spreading in multiplex networks Competing spreading processes on multiplex networks: awareness and epidemics A new coupled disease-awareness spreading model with mass media on multiplex networks Epidemic propagation with positive and negative preventive information in multiplex networks Co-evolution spreading of multiple information and epidemics on two-layered networks under the influence of mass media If the patients decide not to tell what can we do?-TB/HIV counsellors' dilemma on partner notification for HIV HIVprevention measures on a university campus in South Africa-perceptions, practices and needs of undergraduate medical students Interplay between epidemic spread and information diffusion on twolayered networks with partial mapping Epidemic spreading on weighted adaptive networks Epidemic dynamics on an adaptive network Epidemic spreading in modular time-varying networks Epidemic spreading in time-varying community networks Effective degree household network disease model Activity driven modeling of time varying networks Availability of data and material All data and materials can be obtained under the proposed model and algorithms. The authors declare that they have no conflict of interest.Code availability All algorithms are implemented in Java, and the code is available once the request to the corresponding author is reasonable.