Noise-induced transitions in a non-smooth SIS epidemic model with media alert

Abstract: We investigate a non-smooth stochastic epidemic model with consideration of the alerts from media and social network. Environmental uncertainty and political bias are the stochastic drivers in our mathematical model. We aim at the interfere measures assuming that a disease has already invaded into a population. Fundamental findings include that the media alert and social network alert are able to mitigate an infection. It is also shown that interfere measures and environmental noise can drive the stochastic trajectories frequently to switch between lower and higher level of infections. By constructing the confidence ellipse for each endemic equilibrium, we can estimate the tipping value of the noise intensity that causes the state switching.


Introduction
These days have witnessed the horrible pandemic of COVID-19. As is known to the public, key characteristic of this killer is its silent spreads. Before we fully understand the virus, it has already established. In general, what do we do once a disease has invaded into a population? To mitigate the spread of a disease and eventually get a full control to the disease, the behavior changes at the population level can bend the infection curve. How to promote a massive behavior change? Two typical governing systems function differently. A central government can force a mandatory behavior change without any court challenge, period. If a mandatory preventive measure cannot be exercised or the administration failed to act following a scientific manner, media and social network are able to propagate the positive message in order to persuade behavior changes. A behavior change at the population level is highly media oriented. In fighting COVID-19 pandemic, WHO Director-General Tedros Adhanom Ghebreyesus ever mentioned "We're not just fighting an epidemic; we're fighting an infodemic." For emerging and re-emerging epidemic diseases, research has shown that positive media coverage is an effective measure to control the diseases [1].
Massive media propagation and social network response to an epidemic take off only after an outbreak has been observed, so that the initial infection growth is following the regular epidemic process. The study here focuses on how to mitigate an epidemic through the awareness by media and social network. Media and social network response to an epidemic have to be based upon the number of known infected. That is, the pressure from media and social network impose to the population only until the prevalence reaches certain tipping point. We, therefore, will set up a cutting point for the number of infections, then let the media or the social network play its role in our model. A quick searching across the board of relevant literature, this question has been tackled mainly from the view of deterministic modeling. Scholars have proposed media functions to model the transmission rates [2][3][4][5][6][7]. For instance, Cui et al. [3] analyzed a susceptible-exposed-infected model by using β(1 − e −mI ) to characterize the effect of media coverage. Their results showed that weak media effects may not be helpful because the infections oscillate periodically. Misra et al. [8] formulated an SIS model to study the impact of awareness programs conducted by a media campaign on the spread of an infectious disease and showed that the spread of an infectious disease can be controlled by using awareness programs but the disease remains endemic due to immigration.
Research work in highlighting the role of media coverage and social media network without tipping point can be read in [2,3,[5][6][7]. In reality, because of time delay and asymptotic infections, we consider that only when the number of infected individuals reaches a tipping value, the alerts and the pressure from media and social network can take place [9,10]. Liu et al. [4] set a base line of infection (the tipping point) I c to launch media alert. They consequently proposed a non-smooth but continuous transmission rate: where nonnegative constant p denotes the intensity of the media coverage. By assuming that the susceptible individuals follow the logistic growth and newborns directly enter into the susceptible class, they studied the following model: where S (t) and I(t) are the sizes of susceptible and infectious population at time t, respectively. All parameters are positive. r is the intrinsic growth rate of the susceptible population; a is the carrying capacity of the community in the absence of infection; d denotes the natural death rate; γ represents the recovery rate of the infected individuals; is the disease-induced death rate. Liu et al. [4] conducted a rigorous and complete mathematical analysis to model (2). Deterministic models always give us a good start to catch certain characteristic of a disease transmission. Obviously, these mean-field models are based upon the homogeneous population, uniform individual behavior, and steady environmental settings. In managing to fight an epidemic or pandemic infections, unexpected random factors also play a significant role. A certain false claims, misinformation, conspiracy theories can change the course of epidemic. The rumors such as black people are immune to COVID-19 and coconut oil kills the virus of COVID-19 have misled some communities in battling the virus. Spreading false or misleading information may prevent the timely and effective adoption of appropriate behaviours and of public health recommendations or measures [11]. Occasional political bias or unusual unscientific statements can even create a temporal chaos in favor of the disease spread. In the real world, the transmissions of infectious diseases are inevitably affected by these unexpected factors. The variations of physical environmental factors [12], such as humidity, temperature, etc. can have a significant impact on the transmission of a disease. A well-known work was the experimental result of Lowen et al. [13].
The introduction of environmental noise into deterministic setup could exhibit unexpected behaviors which have no analogue in the deterministic models [14]. In the last few decades, researchers have done a lot of work regarding to this aspect [14][15][16][17][18][19][20][21][22][23][24][25]. In [24], a stochastically forced predator-prey model with environmental toxins was proposed to analyze the dynamical behavior of noise-induced transitions from coexistence to prey-only extirpation in the bistable zone. Nonlinear dynamical models can exhibit various new phenomena, such as stochastic resonance [26], noise-induced transitions [27], noise-induced ordering [28], noise-induced chaos [29], and noise-induced complexity [30]. Studies of the environmental noise effects on systems with multiple attractors have been steady growing [31]. Multi-stable systems of deterministic models can exhibit noise-induced state switching between deterministic attractors [32][33][34][35] when a noise is considered. Using the stochastic sensitivity functions (SSF), Bashkirtseva et al. [36] developed a new technique to construct the analytical description of stochastically forced equilibria and cycles of discrete-time models. In the present paper, we will use the SSF technique to explore the noise-induced state switching between two epidemic equilibria and estimate a threshold value of the noise intensity.
We organize this paper as following. In Section 2, we review the main results of deterministic model (2) and propose its stochastic counterpart. The studies of noise-induced transitions between two epidemic equilibria and construction of the confidence ellipses for equilibria will be presented in Section 3. Finally, we conclude our study by a simple discussion in Section 4.

Model formulation and the global existence and uniqueness of the solutions
We first briefly review the main results of deterministic model (2), then we study the corresponding stochastic counterpart.

A summary to the deterministic model
Notice that model (2) always exists two equilibria: the origin O(0, 0) and the disease-free equilibrium E 0 (a, 0). The basic reproductive number is Here, we highlight the interested results from [4], which will be directly used in the subsequent stochastic analysis.
Lemma 2.1. If R 0 < 1, E 0 (a, 0) is globally asymptotically stable. If R 0 > 1 together with certain extra conditions, in addition to saddles O(0, 0) and E 0 (a, 0), there exist three endemic equilibria E 1 , E 2 and E 3 , where E 1 and E 3 are asymptotically stable and E 2 is a saddle. The separatrices (the stable manifolds of the saddle E 2 ) separate the domain of attractions of E 1 and E 3 . Hence, bistable dynamics could be observed in model (2) if R 0 > 1.
In this paper, we pay more attentions to the epidemiological implications. Model (2) was set to investigate the media effects after the disease has already invaded the population. So the media alert and the pressure from social network have no impact to the initial growth of the disease. That is why the basic reproductive number R 0 is independent of the media alert parameters I c and p. The appearance of the bistable dynamics is the straightforward outcome of the introduction of media impact. A bifurcation diagram is shown in Figure 1. It can be seen that model (2) undergoes a backward bifurcation at R 0 ≈ 4.8705. There is a window (4.3515, 4.8705) for R 0 where model (2) can have two stable equilibria E 1 and E 3 and a saddle E 2 (see also Figure 2 (a) where β = 0.0139). The appearance of the stable endemic state competing with the other stable endemic state is known as the bistability phenomenon. The media is helpful to interfare the infections by forcing (2) to choose smaller endemic E 1 as the attractor.

The stochastic model
In this paper we intend to explore the stochastic version of model (2). Stochastic factors such as the environmental noise, social noise, random misinformation spread though media and social network and political motivated biased information are not easy to handle with deterministic models. We integrate all these random factors into the deterministic model (2) to look into the following stochastic differential equation: where B 1 and B 2 are standard one-dimensional independent Brownian motions. We use σ S and σ I to denote the noise intensities. For the simplicity of discussion in this research, we assume σ S = σ I = σ.  Table 1.
Because system (4) is piecewise defined, we distinguish the noise intensities by σ = σ 1 when I > I c and σ = σ 2 when I ≤ I c .

Global existence and uniqueness of solutions
Since we are more interested in the asymptotic behavior of the model, we first anticipate that all solutions exist for all forward times. Now we prove the existence and uniqueness of the global positive solution of stochastic model (4).
In the following, we study the phenomenon of noise-induced transitions between two stochastic attractors for stochastic model (4).

Analysis of noise-induced transitions
In this section, the phenomenon of noise-induced transitions between two endemic equilibria is studied by constructing confidence ellipses with the technique of stochastic sensitivity functions (SSF).
We take the parameter values in Table 1 Table 1 is shown, where the red-dotted curve Γ (the stable manifold of E 2 ) is the separatrix of two basins of attraction and the gray one (I = I c ) is the separatrix of two subsystems. As one can see, for the deterministic model, the trajectories with initial value inside of Γ tend to E 1 , and the trajectories with initial value outside of Γ approach E 3 . For the stochastic model, we take endemic equilibria E 3 as an example. It is worth noting that endemic equilibria E 3 is located above the threshold I c . When the noise intensity is small, here, we take σ 1 = 0.1, the stochastic trajectories with initial value near E 3 will hover around it (see Figure 3). However, as the noise intensity increases (σ 1 = 0.25), the steady state switches from high epidemic E 3 to low epidemic E 1 (see Figure 4). It is also observed that the noise-induced transitions from low epidemic E 1 to high epidemic E 3 .

Noise-induced state frequent switching
In the following, we construct confidence ellipses for stochastic model (4) by applying the stochastic sensitivity function technique (see the Appendix in [37]). We then further estimate the critical value of the noise intensity that causes state switching.
For an endemic equilibriumĒ = (S ,Ī) of deterministic model (2), let F = f 11 f 12 f 21 f 22 be the Jacobian matrix of ( f, g) with respect to (S , I) atĒ, where We define g 11 =S , g 22 =Ī, and G = g 11 0 0 g 22 . After solving the following system of linear we obtain the stochastic sensitivity matrix W = w 11 w 12 w 21 w 22 . Given that a noise intensityσ and a fiducial probability P, then applying the formula A. Fixing fiducial probability P = 0.95, we next observe how different intensity noise affects the confidence ellipse. When σ 1 = 0.1, the corresponding stochastic states and confidence ellipses of subsystem (5) are shown in Figure 5. As shown in the figure, the random states of the stochastic model are distributed around the deterministic endemic equilibrium E 3 , and located inside of the confidence ellipses with probability 0.95.
We then take two larger noise intensities σ 1 = 0.24 and σ 1 = 0.28, the corresponding confidence ellipses are shown in Figure 6. As the noise intensity increases, the confidence ellipse is growing. It crosses separatrix Γ (the stable manifold of E 2 ) and ultimately arrives at the domain of attraction of the low endemic equilibrium E 1 . The noise intensity value that corresponds to the confidence ellipse intercepts with Γ can be regarded as the critical noise intensityσ 1 . Hereσ 1 ≈ 0.24. A typical example that the trajectory of stochastic model (4) starting near the endemic equilibrium E 3 converges to the other endemic equilibrium E 1 is shown in Figure 4 in which σ 1 = 0.25 >σ 1 . Putting the confidence ellipses for E 1 and E 3 together can give us a better looking of the whole system. For weaker noise (σ 1 = 0.1, σ 2 = 0.02), confidence ellipses of both E 1 (green) and E 3 (blue) are located inside their own attraction basin. They are distinctly separated by the separatrix Γ (in red) (see Figure 7 (a)). The random trajectories starting from inside of the separatrix will fluctuate around the endemic equilibrium E 1 , and the ones starting from outside of the separatrix will fluctuate around the endemic equilibrium E 3 . The random trajectories leaving the unforced deterministic attractors concentrate in their small neighborhoods (see Figure 7 (c)). One can see the weaker noise does not affect the dynamics of the stochastic model (4). We fix σ 2 = 0.02 and increase σ 1 to 0.28, then the confidence ellipse (blue) of E 3 crosses the separatrix (see Figure 8 (a)). The stochastic trajectories starting within the attraction basin of E 3 will cross the separatrix and approach E 1 with high probability (see Figure 8 (c)). This is also shown by the time series (see Figure 8 (b)): the number of infected cases ultimately is stabilized around I 1 (a smaller value). To better understand this, we numerically obtain the stationary probability density (SPD) of model (4). In Figure 8 (d), we can see that the SPD for stochastic model (4) has a maximum value at I 1 , which implies that the sample trajectory will stay for a longer time in the neighborhood of I 1 . In other words, I 1 is stable in the meaning of probability (with a bigger probability). This illustrates that the phenomenon of noise-induced transition occurs as the noise intensity σ 1 increases. We also obtain a symmetric result by fixing σ 1 = 0.1 and increasing σ 2 to 0.04 (see Figure 9). In this scenario, the confidence ellipse of E 1 crosses the separatrix Γ. That is, the number of infected cases ultimately maybe jump to I 3 (a bigger value), thus making the infections worse.
Finally, we let both intensities change, increasing σ 1 to 0.28 and σ 2 to 0.04. As expected that the phenomenon of frequent random hopping occur between two attractors E 1 and E 3 is observed. We can see from Figure 10 Figure 9. Confidence ellipses (a), time series (b), random trajectory (c), and stationary probability density (d) for stochastic model (4) with initial value (146, 77) and noise intensity σ 1 = 0.1, σ 2 = 0.04. and the region that the two confidence ellipses intersect builds a "transition bridge" between basins of attraction. In this overlapped region, noise-induced transitions can occur with greater probability (see Figure 10 (b) and (c) ). In addition, we can also see from Figure 10 (d) that the SPD for stochastic model (4) has two maximum value at I 1 and I 3 . In this case, stochastic force can worse the picture of infection, it can also mitigate an infection.  Figure 10. Confidence ellipses (a), time series (b), random trajectory (c), and stationary probability density (d) for stochastic model (4) with initial value (146, 77) and noise intensity σ 1 = 0.28, σ 2 = 0.04.

Discussion
What do we do "if the horse has left the barn"? We are still optimistic. Never too late to take an action. Scientific oriented media alert and positive energy spread over the social network are able to guide us walking out of the woods. For epidemics associated with emerging and re-emerging infectious diseases, we can bend the infection curve if we are oriented scientifically. We can smartly use the media coverage and social network to help us to fight an epidemic. Using mathematical model we have shown that alerts from media and/or social network play an important role in the battle of fighting the disease spread.
The final size of the epidemic depends on the media and social network coverage parameters p and I c . A clear monotonic patten of the final size can be seen shown in Figure 11 (a). Under the same noise intensity (we set σ 1 = 0.005), the greater the media and social network coverage, the smaller the final size. The same monotonic patten is observed in Figure 11 (b), the final infection scale decreases when the threshold I c decreases. This suggests that the decision makers should set appropriate threshold I c to minimise the infection scale. Since the alerts from media social network are triggered only after an initial epidemic, it does reduce the magnitude of the infection. However, one cannot overstate its role in containing an epidemic, mainly because these alerts do not modify the basic reproductive number. It is out of the reach of these alerts to bring the basic reproductive number less than one.
The magnitude of a confidence ellipse is proportional to the intensity of the noise. The bigger the noise, the bigger the corresponding confidence ellipse. Then we are sure that intensifying the noise will raise the likelihood of endemic state shift. However, in terms steady state mobile attributed to the noise, current model cannot distinguish whether it bring the epidemic downward or upward.