Processing math: 47%
Research article

Role of media coverage in a SVEIR-I epidemic model with nonlinear incidence and spatial heterogeneous environment


  • Received: 10 June 2023 Revised: 16 July 2023 Accepted: 21 July 2023 Published: 28 July 2023
  • In this paper, we propose a SVEIR-I epidemic model with media coverage in a spatially heterogeneous environment, and study the role of media coverage in the spread of diseases in a spatially heterogeneous environment. In a spatially heterogeneous environment, we first set up the well-posedness of the model. Then, we define the basic reproduction number R0 of the model and establish the global dynamic threshold criteria: when R0<1, disease-free steady state is globally asymptotically stable, while when R0>1, the model is uniformly persistent. In addition, the existence and uniqueness of the equilibrium state of endemic diseases were obtained when R0>1 in homogeneous space and heterogeneous diffusion environment. Further, by constructing appropriate Lyapunov functions, the global asymptotic stability of disease-free and positive steady states was established. Finally, through numerical simulations, it is shown that spatial heterogeneity can increase the risk of disease transmission, and can even change the threshold for disease transmission; media coverage can make people more widely understand disease information, and then reduce the effective contact rate to control the spread of disease.

    Citation: Pengfei Liu, Yantao Luo, Zhidong Teng. Role of media coverage in a SVEIR-I epidemic model with nonlinear incidence and spatial heterogeneous environment[J]. Mathematical Biosciences and Engineering, 2023, 20(9): 15641-15671. doi: 10.3934/mbe.2023698

    Related Papers:

    [1] Cunjuan Dong, Changcheng Xiang, Wenjin Qin, Yi Yang . Global dynamics for a Filippov system with media effects. Mathematical Biosciences and Engineering, 2022, 19(3): 2835-2852. doi: 10.3934/mbe.2022130
    [2] Cheng-Cheng Zhu, Jiang Zhu, Xiao-Lan Liu . Influence of spatial heterogeneous environment on long-term dynamics of a reaction-diffusion SVIR epidemic model with relaps. Mathematical Biosciences and Engineering, 2019, 16(5): 5897-5922. doi: 10.3934/mbe.2019295
    [3] Yanni Xiao, Tingting Zhao, Sanyi Tang . Dynamics of an infectious diseases with media/psychology induced non-smooth incidence. Mathematical Biosciences and Engineering, 2013, 10(2): 445-461. doi: 10.3934/mbe.2013.10.445
    [4] Yongli Cai, Yun Kang, Weiming Wang . Global stability of the steady states of an epidemic model incorporating intervention strategies. Mathematical Biosciences and Engineering, 2017, 14(5&6): 1071-1089. doi: 10.3934/mbe.2017056
    [5] Cheng-Cheng Zhu, Jiang Zhu . Spread trend of COVID-19 epidemic outbreak in China: using exponential attractor method in a spatial heterogeneous SEIQR model. Mathematical Biosciences and Engineering, 2020, 17(4): 3062-3087. doi: 10.3934/mbe.2020174
    [6] Wenjie Qin, Jiamin Zhang, Zhengjun Dong . Media impact research: a discrete SIR epidemic model with threshold switching and nonlinear infection forces. Mathematical Biosciences and Engineering, 2023, 20(10): 17783-17802. doi: 10.3934/mbe.2023790
    [7] Andrei Korobeinikov, Philip K. Maini . A Lyapunov function and global properties for SIR and SEIR epidemiological models with nonlinear incidence. Mathematical Biosciences and Engineering, 2004, 1(1): 57-60. doi: 10.3934/mbe.2004.1.57
    [8] Xinyu Liu, Zimeng Lv, Yuting Ding . Mathematical modeling and stability analysis of the time-delayed SAIM model for COVID-19 vaccination and media coverage. Mathematical Biosciences and Engineering, 2022, 19(6): 6296-6316. doi: 10.3934/mbe.2022294
    [9] Baoxiang Zhang, Yongli Cai, Bingxian Wang, Weiming Wang . Dynamics and asymptotic profiles of steady states of an SIRS epidemic model in spatially heterogenous environment. Mathematical Biosciences and Engineering, 2020, 17(1): 893-909. doi: 10.3934/mbe.2020047
    [10] Maoxing Liu, Yuhang Li . Dynamics analysis of an SVEIR epidemic model in a patchy environment. Mathematical Biosciences and Engineering, 2023, 20(9): 16962-16977. doi: 10.3934/mbe.2023756
  • In this paper, we propose a SVEIR-I epidemic model with media coverage in a spatially heterogeneous environment, and study the role of media coverage in the spread of diseases in a spatially heterogeneous environment. In a spatially heterogeneous environment, we first set up the well-posedness of the model. Then, we define the basic reproduction number R0 of the model and establish the global dynamic threshold criteria: when R0<1, disease-free steady state is globally asymptotically stable, while when R0>1, the model is uniformly persistent. In addition, the existence and uniqueness of the equilibrium state of endemic diseases were obtained when R0>1 in homogeneous space and heterogeneous diffusion environment. Further, by constructing appropriate Lyapunov functions, the global asymptotic stability of disease-free and positive steady states was established. Finally, through numerical simulations, it is shown that spatial heterogeneity can increase the risk of disease transmission, and can even change the threshold for disease transmission; media coverage can make people more widely understand disease information, and then reduce the effective contact rate to control the spread of disease.



    In epidemiological theory, environmental heterogeneity has been regarded as an indispensable factor in the transmission of infectious diseases. Due to altitude, temperature, humidity, latitude, climate, life factors and other factors, the spread of epidemics in different environments varies greatly. In the process of COVID-19 transmission, due to differences in population density, medical resources and climatic conditions in different regions, the outbreak degree of different regions is also different. For example, during the first outbreak, the epidemic spread faster and more widely in southern and large parts of eastern China, because these areas were densely populated, developed rapidly, and also easy to move people, resulting in a concentrated performance at the peak of the outbreak. In addition, Smith et al.[1] observed temporary fluctuations in mosquito populations affected by environmental variables such as rainfall, humidity and temperature. These spatio-temporal heterogeneity have a significant impact on the mode of transmission of mosquito-borne diseases. Wang et al.[2] found that cholera transmission is influenced by spatial variations and seasonal fluctuations, so disease control measures should pay more attention to the spread process. Cai et al.[3] showed that combinations of spatial heterogeneity tended to enhance the persistence of influenza disease in the model, in other words, influenza infection risk would be very important if spatial heterogeneity was taken into account. Luo et al.[4] observed that dengue is a climate-sensitive disease, and that climatic factors, such as temperature and rainfall, promote dengue transmission through potential changes in vector mosquito density and human behavior, leading to increased exposure to dengue, especially in areas where temperatures favor dengue and other vector-borne diseases. Therefore, in order to study the effects of heterogeneous environments on disease transmission, more and more reaction-diffusion epidemic models have been developed [4,5,6,7,8,9,10,11,12].

    As we all know, clinical results show that for certain diseases, including hepatitis B (HBV)[13], hepatitis C (HCV)[14], most human tuberculosis [15], herpes virus and other infectious diseases, since the recovered person still carries the pathogen, the recovered person may appear with the reactivation of the infection when the pathogen lurking in the tissue reproduces to a certain extent. For example, herpes simplex virus type 2 (HSV-2) is usually transmitted through close physical or sexual contact and can cause genital herpes [16]. The main incidence of genital herpes is due to its frequent recurrence rate. In one study [17], 89% of HSV-2 patients had at least one recurrence during follow-up, with a mean monthly recurrence rate of 34%, 38% of patients relapsed at least six times in the first year, and 20% relapsed more than 10 times. They concluded that almost all initially symptomatic patients with HSV-2 infection had a relapse of symptoms, more than 35% of such patients relapse frequently. Recurrence rates were particularly high in patients with prolonged onset of first infection, whether or not they received acyclovir antiviral chemotherapy. In fact, relapse has been extensively studied as an important feature of human or animal diseases, see [18,19,20,21,22,23]. Tudor[24] was the first to incorporate recurrence into mathematical models, building a bilinear compartment model of morbidity and constant population size (later called the SIRI model). The results show that the basic reproduction number is a threshold parameter for the stability of the system.

    On the one hand, people's perception of disease risk is influenced by information widely disseminated by the media. In the early stages of an outbreak, due to lack of medical diagnosis and vaccination, effective media coverage can reduce infection peaks [25]. In order to characterize the impact of media coverage on disease transmission and control, more and more mathematical models have been proposed. For example, Cui et al. [26] proposed an infectious disease model with Logistic growth under media coverage, with β(I)=μemI as the effective contact rate and m as the influence parameter of media coverage on the effective contact rate. In [27,28], the model adopted β1β2Im+I as the effective contact rate to reflect the reduced amount of contact rate due to media coverage. In [29], the model uses βi(I)=aibim(I) as the effective contact rate, the results suggest that changes in human behavior in response to media coverage can reduce outbreaks and reduce the speed of disease transmission. In [30], the author assumes that the media reports of infectious function of f(I):=(ββ1Im+I)I, and established an age-structured epidemic model for a class of incomplete vaccination. In [31], to model the impact of sanitation and awareness on infectious disease transmission, the authors selected β(M)=(ββ1k1Mp+k1M) as the effective contact rate of susceptible individuals. They showed through model analysis that health and advocacy programs have the ability to lower epidemic thresholds and thus control the spread of infection.

    However, these models based on ordinary differential equations ignore the important factor of spatial heterogeneity. Undoubtedly, it would be more practical to consider the impact of media coverage on disease transmission in heterogeneous environments, but there has been little research in this area. In a recent study, Song et al. [32] studied a class of diffuse epidemic systems with spatial heterogeneity and lagging effects of media influence, where the media influence function is β(x)em(x)I(x,tr). Inspired by [27,28], we suggest that, as the number of infected people increases, corresponding interventions reported by the media can provide profound psychological cues to susceptible individuals to help them reduce their contact with infected individuals. Therefore, the contact rate between susceptible and infected individuals should be assumed to be a monotonically decreasing function. Motivated by the above works, in this paper, we use the general smooth function β(x)β1(x)f(I) as the effective contact rate under the influence of media coverage in a heterogeneous environment, where f(I) is the media coverage function with saturated psychological effect. Consider that "asymptomatic" individuals can still spread the infection even if they do not show any symptoms (for example, a particular but critical feature of the recent COVID-19 epidemic is that a large proportion of the population infected with SARS-Cov-2 virus originates from asymptomatic individuals [33,34].) The main reason is that asymptomatic cases are often unrecognized during the incubation period and therefore have more contacts than symptomatic cases. Therefore, this paper incorporates the transmission of exposed persons to susceptible populations into the mathematical epidemic model. Additionally, since the latent has no obvious symptoms during the incubation period, the effective contact rate is not affected by media reports.

    On the other hand, we also know that vaccination is an important means of preventing infection and relapse of these diseases[35,36,37]. Hence, the effect of spatial factors on relapse and vaccination should be taken into account. Inspired by the above literature, in this paper, we study a SVEIR-I epidemic response-diffusion model that incorporates the impact of media coverage on disease transmission and considers the spread of asymptomatic infected individuals.

    {S(t,x)t=(D1(x)S)+Λ(x)[β(x)β1(x)f(I)]SIβ2(x)SEp(x)Sμ(x)S,V(t,x)t=(D2(x)V)+p(x)Sσ(x)[β(x)β1(x)f(I)]VIσ(x)β2(x)VEμ(x)V,E(t,x)t=(D3(x)E)+{[β(x)β1(x)f(I)]I+β2(x)E}(S+σ(x)V)α(x)Eμ(x)E,I(t,x)t=(D4(x)I)+α(x)E+ρ(x)Rμ(x)Iη(x)Iδ(x)I,R(t,x)t=(D5(x)R)+δ(x)Iμ(x)Rρ(x)R,(S(0,x),V(0,x),E(0,x),I(0,x),R(0,x))=(ϕ1(x),ϕ2(x),ϕ3(x),ϕ4(x),ϕ5(x)),xΩ, (1.1)

    with the homogeneous Neumann boundary condition(NBC)

    S(t,x)n=V(t,x)n=E(t,x)n=I(t,x)n=R(t,x)n=0,xΩ,t0,

    where ΩRn is a bounded domain with the smooth boundary Ω. n denotes the outward normal derivative on Ω. S(t,x),V(t,x),E(t,x),I(t,x) and R(t,x) represent the population density of susceptible individuals, vaccinated individuals, exposed individuals, infected individuals and recovered individuals at location x and time t, respectively. Di(x)(i=1,2,3,4,5) respectively represent the diffusion rate of corresponding individuals at position x. We assume that Λ(x), Di(x), μ(x), β(x), βi(x)(i=1,2), p(x), σ(x), δ(x), α(x), ρ(x) and η(x) are positive, continuous and bounded on ˉΩ. The meanings of other parameters in the Table 1.

    Table 1.  The meaning of parameters in model (1.1).
    Symbol Meaning
    μ(x) Natural mortality rate at location x.
    p(x) Vaccination rate at position x.
    ρ(x) Relapse rate at position x.
    δ(x) Per-capita recovery (treatment) rate at position x.
    Λ(x) The recruitment rate of S at position x.
    η(x) Disease-related death rate at position x.
    α(x) The conversion rate from exposed to infected individuals.
    β(x) The infection rate of S affected by I at position x.
    β2(x) The infection rate of S affected by E at position x.
    β1(x) The maximum reduced contact rate of S at position x affected by media coverage.
    σ(x) The fraction of V being infected and entering E at position x.

     | Show Table
    DownLoad: CSV

    This paper is structured as follows. In Section 2, we first introduce some spaces, give some symbols and make some assumptions. In Section 3, we present and prove some basic results on the existence and ultimate boundedness of global solutions. In Section 4, we define the basic reproduction number R0 for model (1.1), and give the relationship between the principal eigenvalue and the basic reproduction number R0. In Section 5, the dynamics with R0 as the threshold are established, that is, when R0<1, disease-free steady state is globally asymptotically stable, and when R0>1, the system is uniformly persistent. In Section 6, we consider the global dynamics of a spatially homogeneous model. In Section 7, some numerical examples are given to illustrate our results and show the impact of media coverage and spatial heterogeneity on disease transmission. A brief conclusion is given in Section 8.

    Throughout this paper, we first introduce the following assumptions.

    (A1) f(I) is continuously differentiable on [0,+), and it satisfies f(I)>0, f(0)=0 and limI f(I)=1.

    (A2) The limitations of media coverage indicate β(x)>β1(x) for all t>0 and xˉΩ.

    (A3) The incidence of disease is a strictly monotonically increasing function of I, so {[ββ1f(I)]SI}I>0.

    Let X:=C(ˉΩ,R5), Y:=C(ˉΩ,R3) and Z:=C(ˉΩ,R) are Banach spaces with the supremum norm of continuous functions. Furthermore, X+:=C(ˉΩ,R5+), Y+:=C(ˉΩ,R3+) and Z+:=C(ˉΩ,R+) are the positive cones corresponding to X, Y and Z, respectively. Then, (X,X+), (Y,Y+) and (Z,Z+) are ordered Banach spaces. For any function (x) defined on some nonempty set U, we denote

    :=maxxU(x),:=minxU(x).

    Let Ai(t)(i=1,2,3,4,5):C(ˉΩ,R)C(ˉΩ,R) are the C0 semigroup associated with (Di(x))πi(x) subjects to NBC, where π1=μ(x)+p(x), π2=μ(x), π3=μ(x)+α(x), π4=μ(x)+η(x)+δ(x) and π5=μ(x)+ρ(x). Then, we have

    (Ai(t)ϕ)(x)=ΩGi(t,x,y)ϕ(y)dy,t>0,ϕC(ˉΩ,R), (2.1)

    where Gi(t,x,y) is the Green function with (Di(x))πi(x)(i=1,2,3,4,5) subjects to NBC. Furthermore, based on ([38], Corollary 7.2.3), we know that Ai(t) is compact and strongly positive for all t>0. Thus, there exist constants Mi>0 such that Ai(t)Mieαit for all t0, where αi<0 are the principal eigenvalue of (Di(x))πi(x) under NBC. For any initial value function ϕ=(ϕ1,ϕ2,ϕ3,ϕ4,ϕ5)X+, we denote u(t,,ϕ)=(S(t,,ϕ),V(t,,ϕ), E(t,,ϕ),I(t,,ϕ),R(t,,ϕ))T be the solution of model (1.1). Let

    {F1(ϕ)(x)=Λ(x)[β(x)β1(x)f(ϕ4(x))]ϕ1(x)ϕ4(x)β2(x)ϕ1(x)ϕ3(x),F2(ϕ)(x)=pϕ1(x)σ(x)[β(x)β1(x)f(ϕ4(x))]ϕ2(x)ϕ4(x)σ(x)β2(x)ϕ2(x)ϕ3(x),F3(ϕ)(x)=[β(x)β1(x)f(ϕ4(x))]ϕ1(x)ϕ4(x)+β2(x)ϕ1(x)ϕ3(x)+σ(x)[β(x)β1(x)f(ϕ4(x))]ϕ2(x)ϕ4(x)+σ(x)β2(x)ϕ2(x)ϕ3(x),F4(ϕ)(x)=α(x)ϕ3(x)+ρ(x)ϕ5(x),F5(ϕ)(x)=δ(x)ϕ4(x),

    then model (1.1) can be written in the following form:

    {u(t,,ϕ)=A(t)ϕ+t0A(ts)F(u(s,,ϕ))ds,u(0)=ϕ, (2.2)

    where A(t)=diag(A1(t),A2(t),A3(t),A4(t),A5(t)) and F(ϕ)(x)=diag(F1(ϕ)(x),F2(ϕ)(x),F3(ϕ)(x), F4(ϕ)(x),F5(ϕ)(x))T.

    In this section, we give the existence of global solutions and the ultimate boundedness of system (1.1).

    Lemma 1. For any ϕ=(ϕ1,ϕ2,ϕ3,ϕ4,ϕ5)X+, model (1.1) has a unique non-negative solution u(t,,ϕ)=(S(t,,ϕ),V(t,,ϕ),E(t,,ϕ),I(t,,ϕ),R(t,,ϕ))X+ on [0,τ) and τ. Moreover, this solution is a classical solution.

    Proof. Similar to the proof of in ([10], Lemma 1), it is sufficient to show the following subtangential conditions holds

    limh0+dist(ϕ+hF(ϕ),X+)=0,ϕX+.

    We can prove that for any ϕX+ and a sufficiently small h0,

    ϕ(x)+hF(ϕ)(x)=(ϕ1(x)+h{Λ(x)[β(x)β1(x)f(ϕ4(x))]ϕ1(x)ϕ4(x)β2(x)ϕ1(x)ϕ3(x)}ϕ2(x)+h{p(x)ϕ1(x)σ(x)[β(x)β1(x)f(ϕ4(x))]ϕ2(x)ϕ4(x)σ(x)β2(x)ϕ2(x)ϕ3(x)}ϕ3(x)+h{[(β(x)β1(x)f(ϕ4(x)))ϕ4(x)+β2(x)ϕ3(x)][ϕ1(x)+σ(x)ϕ2(x)]}ϕ4(x)+h[α(x)ψ3(x)+ρ(x)ϕ5(x)]ϕ5(x)+hδ(x)ϕ4(x))(ϕ1(x)h[β(x)β1(x)f(ϕ4(x))]ϕ1(x)ϕ4(x)ϕ2(x)hσ(x)[β(x)β1(x)f(ϕ4(x))]ϕ2(x)ϕ4(x)ϕ3(x)ϕ4(x)ϕ5(x))>0,

    which indicates that ϕ+hF(ϕ)X+. That completes the proof.

    To continue our research, first, consider the following system,

    {ω(t,x)t=(d(x)ω(t,x))+α(x)β(x)ω(t,x),xΩ,t0,ω(t,x)n=0,xΩ,t0, (3.1)

    where d(x), α(x) and β(x) are positive continuous functions defined on ˉΩ, and we have the following lemma.

    Lemma 2. [39] System (3.1) admits a unique positive steady state ω0(x), which satisfies the equation

    {(d(x)ω0(x))+α(x)β(x)ω0(x)=0,xΩ,ω0(x)n=0,xΩ,

    and it is globally asymptotically stable in C(ˉΩ,R+). In addition, ω0(x)=ab when α()a and β()b are positive constants.

    Theorem 1. For any ϕ=(ϕ1,ϕ2,ϕ3,ϕ4,ϕ5)X+, model (1.1) has a unique classical solution u(t,,ϕ) =(S(t,,ϕ), V(t,,ϕ), E(t,,ϕ), I(t,,ϕ),R(t,,ϕ))X+ defined on [0,), and this solution is also ultimately bounded.

    Proof. Using a standard argument (see, e.g., ([40], Theorem 3.2)), it is only necessary to get the boundness of solution in ˉΩ×[0,τ). From the first equation of (1.1), we have

    St(D1(x)S)+ΛpSμS,t[0,τ),xˉΩ. (3.2)

    From the comparison principle [41] and Lemma 2, we can see that there exists a constant M1>0 such that S(t,x)M1 for all t[0,τ) and xˉΩ. From the second equation of (1.1), we have

    Vt(D2(x)V)+pM1μV,t[0,τ),xˉΩ. (3.3)

    Through similar analysis, there exists a constant M2>0 such that V(t,x)M2 for all t[0,τ) and xˉΩ.

    Then we can get

    {E(t,x)t(D3(x)E)+βM1I+β2M1E+σβM2I+σβ2M2E(α+μ)E,I(t,x)t(D4(x)I)+αE+ρRμIηIδI,R(t,x)t(D5(x)R)+δIμRρR,t[0,τ),xΩ,En=In=Rn=0,xΩ,t>0, (3.4)

    Consider the following comparison system:

    {v1t=(D3(x)v1)+βM1v2+β2M1v1+σβM2v2+σβ2M2v1(α+μ)v1,v2t=(D4(x)v2)+αv1+ρv3μv2ηv2δv2,v3t=(D5(x)v3)+δv2μv3ρv3,t>0,xΩ,v1n=v2n=v3n=0,xΩ,t>0. (3.5)

    According to the Krein-Rutman theorem [42], the eigenvalue problem associated with the system (3.5) has a strongly positive eigenfunction ξ=(ξ1,ξ2,ξ3) corresponding to the principal eigenvalue λ. Thus, system (3.5) admits a solution ceλtξ(x) for t0, where c is a positive constant and satisfies cξ=(v1(x,0),v2(x,0),v3(x,0))(E(x,0),I(x,0),R(x,0)) for all xˉΩ. Then, using the principle of comparison

    E(t,x),I(t,x),R(t,x)ceλtξ(x),t[0,τ),xˉΩ.

    This means that there is a constant M3>0 such that

    E(t,x)M3,I(t,x)M3,R(t,x)M3,t[0,τ),xˉΩ.

    Therefore, τ=. This shows the global existence of u(t,,ϕ).

    We next show that the solution is ultimately bounded. According to the comparison principle, from inequality (3.2), (3.3) and Lemma 2, it follows that there exists constants of N1>0, N2>0 and times t1>0, t2>0 such that S(t,x)N1 and V(t,x)N2 for all tmax{t1,t2}, xˉΩ, which means that S(t,x) and V(t,x) are ultimately bounded.

    Let

    P(t)=Ω(S(t,x)+V(t,x)+E(t,x)+I(t,x)+R(t,x))dx,

    Using ([43], Theorem 3.7), we have

    dPdt=Ω[Λ(x)μ(x)S(t,x)μ(x)V(t,x)μ(x)E(t,x)μ(x)I(t,x)η(x)I(t,x)μ(x)R(t,x)]dxΛ|Ω|μP,t[0,τ), (3.6)

    where |Ω| is the volume of Ω. Hence, there exists a constant N3>0 and t3>0 such that P(t)N3 for all tt3. Consequently, we have

    ΩE(t,x)dxN3,ΩI(t,x)dxN3,ΩR(t,x)dxN3. (3.7)

    Then, using the similar method in ([9], Theorem 1), for any tt4=max{t1,t2,t3}, we have

    E(t,x)=A3(t)E(t4,x)+tt4A3(ts){[β(x)β1(x)f(I(s,x))]I(s,x)+β2(x)E(s,x)}[S(s,x)+σ(x)V(s,x)]dsM3eα3(tt4)E(t4,x)+ω3tt4e(α+μ)(ts)[βN1ΩI(s,y)+β2N1ΩE(s,y)+σβN2ΩI(s,y)+σβ2N2ΩE(s,y)]dydsM3eα3(tt4)E(t4,x)+ω3N1N2N3(β+β2+σβ+σβ2)α+μ.

    So, we have

    lim suptE(t,x)Zω3N1N2N3(β+β2+σβ+σβ2)α+μ.

    This shows that E(t,x) is ultimately bounded. Similarly, for I(t,x) and R(t,x) we get

    lim suptI(t,x)Zω4N3(α+ρ)μ+η+δ,lim suptR(t,x)Zω5N3δμ+ρ,

    where ω4>0 and ω5>0 are constants. This indicates that the solution is ultimately bounded.

    Corollary 1. For any ϕX+, the solution semiflow Q(t)ϕ=u(t,,ϕ):X+X+ generated by model (1.1) has a compact global attractor.

    Proof. By Theorem 1, we know that the system (1.1) is ultimately bounded, which means that the solution semiflow Q(t)ϕ=u(t,,ϕ):X+X+ is point dissipative on X+. Furthermore, from ([44], Theorem 2.6), we can see that Q(t) is compact for any t>0. Thus, from ([45], Theorem 1.1.3), we get that Q(t):X+X+ has a compact global attractor on X+.

    According to Lemma 2, the system (1.1) has a unique disease-free steady state D0(x)=(S0(x),V0(x),0,0,0), where S0(x) and V0(x) satisfies

    {(D1(x)S0(x))=Λ(x)p(x)S0(x)μ(x)S0(x),(D2(x)V0(x))=p(x)S0(x)μ(x)V0(x),xΩ,t>0,S0(x)n=V0(x)n=0,xΩ,t>0.

    By linearizing model (1.1) at D0(x), we obtain the following linear system

    {E(t,x)t=(D3(x)E)+[β(x)I+β2(x)E][S0(x)+σ(x)V0(x)](α(x)+μ(x))E,I(t,x)t=(D4(x)I)+α(x)E+ρ(x)Rμ(x)Iη(x)Iδ(x)I,R(t,x)t=(D5(x)R)+δ(x)Iμ(x)Rρ(x)R,xΩ,t>0,En=In=Rn=0,xΩ,t>0. (4.1)

    Let (E(t,x),I(t,x),R(t,x))=(eλtϖ3(x),eλtϖ4(x),eλtϖ5(x)), then we get the following eigenvalue problem

    {λϖ3(x)=(D3(x)ϖ3(x))+[β(x)ϖ4(x)+β2(x)ϖ3(x)][S0(x)+σ(x)V0(x)](α(x)+μ(x))ϖ3(x),λϖ4(x)=(D4(x)ϖ4(x))+α(x)ϖ3(x)+ρ(x)ϖ5(x)μ(x)ϖ4(x)η(x)ϖ4(x)δ(x)ϖ4(x),λϖ5(x)=(D5(x)ϖ5(x))+δ(x)ϖ4(x)μ(x)ϖ5(x)ρ(x)ϖ5(x),xΩ,t>0,ϖ3n=ϖ4n=ϖ5n=0,xΩ,t>0. (4.2)

    From ([38], Theorem 7.6.1), we have the following result.

    Lemma 3. (4.2) has a unique principal eigenvalue λ0=λ0(S0(x), V0(x)) with positive eigenvector (ϖ3(x),ϖ4(x),ϖ5(x)).

    Define

    D=(D3(x)000D4(x)000D5(x)),V=(α(x)+μ(x)00α(x)μ(x)+η(x)+δ(x)ρ(x)0δ(x)μ(x)+ρ(x)),
    F=(β2(x)S0(x)+σ(x)β2(x)V0(x)β(x)S0(x)+σ(x)β(x)V0(x)0000000).

    Then system (4.1) is

    ðt=(D(x)ð)+(FV)ð,

    where ð=(E,I,R)T. Let the positive semigroup P(t):Y+Y+ generated by the following system

    ðt=Bð,

    where B:=DV. Then we know that B is a resolvable positive operator [46] and P(t)Y+Y+ can be obtained. Moreover, we denote the spectral radius of Q as r(Q)=sup{|λ|,λσ(Q)} and the spectral bound of Q as s(Q)=sup{Reλ,λσ(Q)}.

    Denote ψ=(ψ3(x),ψ4(x),ψ5(x)) is initial infection distribution. Based on a similar discussion by Luo et al.[10], we define the following operators

    L(ψ)(x)=+0F(x)P(t)ψ(x)dt=F(x)+0P(t)ψ(x)dt, (4.3)

    where F(x)+0P(t)ψ(x)dt denotes the total distribution of new infections. Then, based on [46,47,48,49], the basic reproduction number for model (1.1) is defined by R0:=r(L). Further, through the argument in [46,47], we get the following lemma.

    Lemma 4. (i) sign(R01)=sign(λ0).

    (ii) If R0<1, then D0(x) of system (1.1) is locally asymptotically stable.

    (iii) If R0>1, then D0(x) is unstable.

    In this section, we study the global stability of the disease-free steady state D0(x) of the system (1.1), and establish the following results.

    Theorem 2. The disease-free steady state D0(x) of the model (1.1) is globally asymptotically stable in X+ when R0<1.

    Proof. According to Lemma 4, we know that D0(x) is locally asymptotically stable when R0<1. Thus we only need to show that D0(x) is globally attractive. From the first equation of model (1.1), we get

    {St(D1(x)S)+Λ(x)μ(x)S,xΩ,t>0,Sn=0,xΩ,t>0.

    Using the comparison principle and Lemma 2, we obtain lim suptS(t,x)S0(x) uniformly for xˉΩ, which means that there exist positive constants t1 and ϵ1 such that S(t,x)S0(x)+ϵ1 for any tt1. Similarly, from the second equation of model (1.1), we have

    {Vt(D2(x)V)+p(x)(S0(x)+ϵ1)μ(x)V,xΩ,t>0,Vn=0,xΩ,t>0.

    and then we have there exist positive constants ϵ2 and t2 such that V(t,x)V0(x)+ϵ2 for any tt2 and xˉΩ. Hence, for any tt3 we have

    {Et(D3(x)E)+β(x)(S0(x)+ϵ1)I+β2(x)(S0(x)+ϵ1)E+σ(x)β(x)(V0(x)+ϵ2)I+σ(x)β2(x)(V0(x)+ϵ2)E(α(x)+μ(x))E,It(D4(x)I)+α(x)E+ρ(x)R(μ(x)+η(x)+δ(x))I,Rt(D5(x)R)+δ(x)I(μ(x)+ρ(x))R,xΩ,En=In=Rn=0,xΩ.

    Let (E(t,x),I(t,x),R(t,x)) be the solution of the following system

    {Et=(D3(x)E)+β(x)(S0(x)+ϵ1)I+β2(x)(S0(x)+ϵ1)E+σ(x)β(x)(V0(x)+ϵ2)I+σ(x)β2(x)(V0(x)+ϵ2)E(α(x)+μ(x))E,It=(D4(x)I)+α(x)E+ρ(x)R(μ(x)+η(x)+δ(x))I,Rt=(D5(x)R)+δ(x)I(μ(x)+ρ(x))R,xΩ,En=In=Rn=0,xΩ. (5.1)

    Using comparison principle, (E(t,x),I(t,x),R(t,x))(E(t,x),I(t,x),R(t,x)). It follows from the Lemma 4 that λ0(S0(x)+ϵ1,V0(x)+ϵ2)<0 when R0<1, where λ0(S0(x)+ϵ1,V0(x)+ϵ2)<0 is the principal eigenvalue of system (5.1). Denote (a3(x),a4(x),a5(x)) be the eigenfunction corresponding to this principal eigenvalue λ0(S0(x)+ϵ1,V0(x)+ϵ2)<0, we can get the following solution for system (5.1),

    (E(t,x),I(t,x),R(t,x))=(a3(x),a4(x),a5(x))eλ0(S0(x)+ϵ1,V0(x)+ϵ2)(tt3),tt3.

    Therefore, (E(t,x),I(t,x),R(t,x))0 uniformly on xˉΩ as t. Further, we get the following limit system

    {St=(D1(x)S)+Λ(x)μ(x)S,xΩ,t>0,Vt=(D2(x)V)+p(x)Sμ(x)V,xΩ,t>0,Sn=Vn=0,xΩ,t>0.

    Based on the theory of asymptotically autonomous semiflows ([50], Corollary 4.3) and Lemma 3, we have S(t,x)S0(x) and V(t,x)V0(x) uniformly for xˉΩ as t. That means D0(x) is globally asymptotically stable.

    In this subsection, we prove that R0 is the threshold for disease persistence. We first give the following results.

    Theorem 3. If R0>1, then (1.1) is uniformly persistent, that is, there is a constant ϱ>0 such that for any initial value ϕ=(ϕ1,ϕ2,ϕ3,ϕ4,ϕ5)X+ with ϕ30, ϕ40 and ϕ50, u(t,,ϕ)=(S(t,,ϕ),V(t,,ϕ),E(t,,ϕ),I(t,,ϕ),R(t,,ϕ)) satisfies

    lim inftS(t,,ϕ)ϱ,lim inftV(t,,ϕ)ϱlim inftE(t,,ϕ)ϱ,lim inftI(t,,ϕ)ϱ,lim inftR(t,,ϕ)ϱ,

    uniformly for xˉΩ. Moreover, model (1.1) has at least one endemic steady state D(x)=(S(x), V(x),E(x),I(x),R(x)).

    Proof. Let

    X0:={ϕ=(ϕ1,ϕ2,ϕ3,ϕ4,ϕ5)X+:ϕ30,ϕ40andϕ50},X0:=X+X0={ϕ=(ϕ1,ϕ2,ϕ3,ϕ4,ϕ5)X+:ϕ3=0orϕ4=0orϕ5=0}.

    and

    M:={ϕX+:Q(t)ϕX0for allt0},M1={D0(x)}.

    First, we give the following two claims:

    Claim 1. ϕMω(ϕ)=M1 for any given ϕM (where ω(ϕ) be the omega limit set of ϕ for Q(t)).

    Proof. Since Q(t)D0(x)=D0(x) for all t0, then M1ϕMω(ϕ). Now, we just have to prove ϕMω(ϕ)M1. Since Q(t)ϕX0 for any given ϕM and t0, then E(t,,ϕ)0 or I(t,,ϕ)0 or R(t,,ϕ)0 for all t0.

    The first case, E(t,x)0. By the third and fourth equations of the system (1.1), we yield

    I(t,x)t+R(t,x)t=(D4(x)I)+(D5(x)R)μ(x)Iη(x)Iμ(x)R,

    by the parabolic maximum principle[41], we immediately obtain I(t,x)=0 and R(t,x)=0 from model (1.1). Thus, we yield

    {S(t,,ϕ)t=(D1(x)S(t,,ϕ))+Λ(x)μ(x)S(t,,ϕ),xΩ,t>0,V(t,,ϕ)t=(D2(x)V(t,,ϕ))+p(x)S(t,,ϕ)μ(x)V(t,,ϕ),xΩ,t>0,S(t,,ϕ)n=V(t,,ϕ)n=0,xΩ,t>0. (5.2)

    According to the comparison principle and Lemma 2, S(t,x)S0(x) and V(t,x)V0(x) uniformly on xˉΩ as t hold. This implies ω(ϕ)=D0(x).

    Second case, I(t,x)0. According to the fourth equation of system (1.1), we know that α(x)E(t,x)+ρ(x)R(t,x)=0. Due to α(x)>0 and ρ(x)>0, we have E(t,x)=0 and R(t,x)=0. Thus, by asymptotically autonomous semiflow theory [50], we get S(t,x)S0(x) and V(t,x)V0(x) uniformly on xˉΩ as t. This also implies ω(ϕ)=D0(x).

    The third case, R(t,x)0. Similar to the first case, we get I(t,x)=0 and E(t,x)=0. Additionally, S(t,x)S0(x) and V(t,x)V0(x) uniformly on xˉΩ as t. This implies ω(ϕ)=D0(x). Thus, we have ϕMω(ϕ)M1, and finally get ϕMω(ϕ)=M1.

    Claim 2. If R0>1, then there exists a constant δ>0 such that for any ϕX0, the solution u(t,,ϕ) of the model (1.1) is satisfied

    lim suptu(t,,ϕ)D0(x)Xδ.

    Proof. Suppose that Claim 2 does not hold, then there exists an enough large T1 such that S0(x)δ<S(t,,ϕ)S0(x)+δ,V0(x)δ<V(t,,ϕ)V0(x)+δ,0E(t,,ϕ)<δ,0I(t,,ϕ)<δ,0R(t,,ϕ)<δ, for all tT1 and xˉΩ. From model (1.1), we have

    {Et(D3(x)E)+[β(x)β1(x)](S0(x)δ)I+β2(x)(S0(x)δ)E+σ(x)[β(x)β1(x)](V0(x)δ)I+σ(x)β2(x)(V0(x)δ)E(α(x)+μ(x))E,It(D4(x)I)+α(x)E+ρ(x)Rμ(x)Iη(x)Iδ(x)I,Rt(D5(x)R)+δ(x)Iμ(x)Rρ(x)R,xΩ,En=In=Rn=0,xΩ. (5.3)

    Study the following corresponding comparison system with (5.3)

    {v3t=(D3(x)v3)+[β(x)β1(x)](S0(x)δ)v4+β2(x)(S0(x)δ)v3+σ(x)[β(x)β1(x)](V0(x)δ)v4+σ(x)β2(x)(V0(x)δ)v3(α(x)+μ(x))v3,v4t=(D4(x)v4)+α(x)v3+ρ(x)v5μ(x)v4η(x)v4δ(x)v4,v5t=(D5(x)v5)+δ(x)v4μ(x)v5ρ(x)v5,xΩ,v3n=v4n=v5n=0,xΩ. (5.4)

    For any initial value ϕX0, it is obtained by the parabolic maximum principle [41] for all tT1 and xˉΩ, there is E(t,x)>0, I(t,x)>0 and R(t,x)>0.

    Next, we consider the following eigenvalue problem

    {λϕϱ3=(D3(x)ϕϱ3)+[β(x)β1(x)](S0(x)δ)ϕϱ4+β2(x)(S0(x)δ)ϕϱ3+σ(x)[β(x)β1(x)](V0(x)δ)ϕϱ4+σ(x)β2(x)(V0(x)δ)ϕϱ3(α(x)+μ(x))ϕϱ3,λϕϱ4=(D4(x)ϕϱ4)+α(x)ϕϱ3+ρ(x)ϕϱ5μ(x)ϕϱ4η(x)ϕϱ4δ(x)ϕϱ4,λϕϱ5=(D5(x)ϕϱ5)+δ(x)ϕϱ4μ(x)ϕϱ5ρ(x)ϕϱ5,xΩ,ϕϱ3n=ϕϱ4n=ϕϱ5n=0,xΩ. (5.5)

    Let λ0(ϵ) be the principal eigenvalue of (5.5). There is a constant that is small enough ϵ0>0, such that λ0(ϵ0)>0 and S0(x)ϵ0>0, V0(x)ϵ0>0 for all xˉΩ. In addition, the eigenvalue λ0(ϵ0) corresponds to a strictly positive eigenfunction (ϕϱ3(x),ϕϱ4(x),ϕϱ5(x)) for xˉΩ. Clearly, system (5.4) has a solution (v3(t,x),v4(t,x),v5(t,x))=eλ0(ϵ0)(tT1)(ϕϱ3(x),ϕϱ4(x),ϕϱ5(x)). Then we can choose a constant c1>0, according to the comparison principle, we get

    (E(t,x),I(t,x),R(t,x))>c1(ϕϱ3(x),ϕϱ4(x),ϕϱ5(x))eλ0(ϵ0)(tT1),tT1.

    Owing to λ0(ϵ0)>0, we have limtE(t,x)=+, limtI(t,x)=+ and limtR(t,x)=+. However, by Theorem 1, we know from the boundedness of (E(t,x),I(t,x),R(t,x)) that this is a contradiction, which shows Claim 2 holds.

    Next, we prove the required persistence. By Theorem 1, the solution u(t,,ϕ) of model (1.1) is ultimately bounded, that is, for a constant M>0 and time T2>0, such that u(t,,ϕ)M for all t>T2 and xˉΩ. Therefore, according to the first and second equations of model (1.1), we get

    {St(D1(x)S)+ΛβSMβ2SMpSμS,xΩ,t>0,Vt(D2(x)V)+pSσβVMσβ2VMμV,xΩ,t>0,Sn=Vn=0,xΩ,t>0,

    for any tT2. It follows from Lemma 2 and the comparison principle [41]

    lim inftS(t,,ϕ)Λ(β+β2)M+μ+p,lim inftV(t,,ϕ)pΛ[(β+β2)M+μ+p][(β+β2)σM+μ],

    which implies that S(t,,ϕ) and V(t,,ϕ) has a positive lower bound. Therefore, the uniform persistence of S(t,,ϕ) and V(t,,ϕ) in model (1.1) is proved.

    Define a continuous function p:X+[0,+) by

    p(ϕ)=min{minxˉΩϕ3(x),minxˉΩϕ4(x),minxˉΩϕ5(x)}.

    Obviously, p1(0,+)X0 and p satisfies the following properties: p(ϕ)=0 and ϕX0 or p(Q(t)ϕ)>0. Based on the definition ([51], Section 3), we know that p is a generalized distance function for semiflow Q(t):X+X+. In addition, as can be seen from the discussion in Claim 1 and Claim 2, all the solution of the model (1.1) tend to D0(x) on boundary X0, which implies that D0(x) is a isolated invariant set in X+, no subset of M1 forms a cycle in X0 and Ws(D0)X0=, where Ws(D0) is the stable set of D0(x). Thus, Ws(D0(x))p1(0,)=. By ([51], Theorem 3), there exists a constant ϱ2>0 such that lim inftp(Q(t)ϕ)ϱ2 for all ϕX0, which implies

    lim inftE(t,,ϕ)ϱ2,lim inftI(t,,ϕ)ϱ2,lim inftR(t,,ϕ)ϱ2.

    Define ϱ1=min{Λ(β+β2)M+μ+p,pΛ[(β+β2)M+μ+p][(β+β2)σM+μ]} and let ϱ=min{ϱ1,ϱ2}, the uniform persistence of the system (1.1) is obtained. From the above proof and Theorem 4.7 in [52], we have the following result.

    Lemma 5. When R0>1, model (1.1) has at least one positive steady state D(x)=(S(x),V(x),E(x), I(x),R(x)).

    In this section, we study the dynamic behavior of model (1.1) in homogeneous space, where all the coefficients of model (1.1) are positive constants except the diffusion coefficient.

    {St=(D1(x)S)+Λ[ββ1f(I)]SIβ2SEpSμS,Vt=(D2(x)V)+pSσ[ββ1f(I)]VIσβ2VEμV,Et=(D3(x)E)+{[ββ1f(I)]I+β2E}(S+σV)αEμE,It=(D4(x)I)+αE+ρRμIηIδI,Rt=(D5(x)R)+δIμRρR,xΩ,t0.Sn=Vn=En=In=Rn=0,xΩ,t0. (6.1)

    Obviously, by Lemma 2, the system (6.1) has a disease-free equilibrium Q^0 = (S^{0}, V^{0}, 0, 0, 0) with S^{0} = \frac{\Lambda}{\mu+p} and V^{0} = \frac{p\Lambda}{\mu(\mu+p)} . By simple computation, one has

    \begin{aligned} R_0& = \frac{\beta_2(\mu \Lambda+\sigma p\Lambda)}{\mu(\mu+p)(\alpha+\mu)}+\frac{\beta\alpha(\mu\Lambda+\sigma p\Lambda)(\mu+\rho)}{\mu(\mu+p)(\alpha+\mu)\big[(\mu+\eta+\delta)(\mu+\rho)-\rho\delta\big]}\triangleq R_{01}+R_{02}. \end{aligned}

    Remark 1. R_{01} and R_{02} can be expressed as the contribution of exposed and infected individuals to the basic reproduction number, respectively. It can be seen that a primary case in the latent population has a rate \frac{\beta_2(\mu\Lambda+\sigma p\Lambda)}{\mu(\mu+p)} of infectious contact with the susceptible population in the expected time \frac{1}{\alpha+\mu} , while a primary case in the infected population has a rate \frac{\beta \alpha(\mu\Lambda+\sigma p\Lambda)(\mu+\rho)}{\mu(\mu+p)} of infectious contact with the susceptible population in the expected time \frac{1}{(\alpha+\mu)[(\mu+\eta+\delta)(\mu+\rho)-\rho\delta]} .

    Theorem 4. When R_0 > 1 , there exists a unique endemic equilibrium for system (6.1).

    Proof. Suppose that there exists an endemic equilibrium Q^*(S^*, V^*, E^*, I^*, R^*) for system (6.1), then Q^* satisfies

    \begin{equation} \left\{ \begin{aligned} &\Lambda-[\beta-\beta_1f(I^*)]S^*I^*-\beta_2S^*E^*-\mu S^*-pS^* = 0, \\ &pS^*-\sigma[\beta-\beta_1f(I^*)]V^*I^*-\beta_2V^*E^*-\mu V^* = 0, \\ &\big\{[\beta-\beta_1f(I^*)]I^*+\beta_2E^*\big\}(S^*+\sigma V^*)-\alpha E^*-\mu E^* = 0, \\ &\alpha E^*+\rho R^*-\mu I^*-\eta I^*-\delta I^* = 0, \\ &\delta I^*-\mu R^*-\rho R^* = 0. \end{aligned} \right. \end{equation} (6.2)

    By a simple computation, we have

    \begin{aligned} S^*& = \frac{\Lambda}{[\beta-\beta_1f(I^*)]I^*+\beta_2E^*+\mu+p}, \; E^* = \frac{[(\mu+\rho)(\mu+\eta+\delta)-\rho\delta ]I^*}{\alpha(\mu+\rho)}, \\ V^*& = \frac{p\Lambda}{\{[\beta-\beta_1f(I^*)]I^*+\beta_2E^*+\mu+p\}\{\sigma[\beta-\beta_1f(I^*)]I^*+\sigma\beta_2E^*+\mu\}}. \end{aligned}

    Since I^*\neq 0 , by the third equation in (6.2), we have

    \begin{aligned} S^*+\sigma V^* & = \frac{(\alpha+\mu)E^*}{[\beta-\beta_1f(I^*)]+\beta_2E^*} = \frac{(\alpha+\mu)[(\mu+\rho)(\mu+\eta+\delta)-\rho\delta]}{\alpha(\mu+\rho)[\beta-\beta_1f(I^*)]+\beta_2[(\mu+\rho)(\mu+\eta+\delta)-\rho\delta]} \end{aligned}

    Let C(I) = \alpha(\mu+\rho)[\beta-\beta_1f(I)]I+\beta_2[(\mu+\rho)(\mu+\eta+\delta)-\rho\delta]I . According to the known formula, we have

    \begin{aligned} S^*+\sigma V^* = & \frac{\Lambda}{[\beta-\beta_1f(I^*)]I^*+\beta_2E^*+\mu+p}\bigg[1+\frac{\sigma p}{\sigma[\beta-\beta_1f(I^*)]I^*+\sigma\beta_2E^*+\mu}\bigg]\\ = &\frac{\alpha\Lambda(\mu+\rho)[\sigma C(I^*)+\alpha(\mu+\rho)(\mu+\sigma p)]}{[C(I^*)+\alpha(\mu+p)(\mu+\rho)][\sigma C(I^*)+\alpha\mu(\mu+\rho)]}. \end{aligned}

    Thus,

    \frac{\alpha\Lambda(\mu+\rho)[\sigma C(I^*)+\alpha(\mu+\rho)(\mu+\sigma p)]}{[C(I^*)+\alpha(\mu+p)(\mu+\rho)][\sigma C(I^*)+\alpha\mu(\mu+\rho)]} = \frac{(\alpha+\mu)[(\mu+\rho)(\mu+\eta+\delta)-\rho\delta]}{\alpha(\mu+\rho)[\beta-\beta_1f(I^*)]+\beta_2[(\mu+\rho)(\mu+\eta+\delta)-\rho\delta]}.

    Define H(I): = \psi (I)-\varphi (I) = 0 , where

    \begin{equation} \begin{aligned} \psi (I)& = \frac{\alpha(\mu+\rho)[\sigma C(I)+\alpha(\mu+\rho)(\mu+\sigma p)]}{[C(I)+\alpha(\mu+p)(\mu+\rho)][\sigma C(I)+\alpha\mu(\mu+\rho)]}, \\ \varphi(I)& = \frac{(\alpha+\mu)\big[(\mu+\rho)(\mu+\eta+\delta)-\rho\delta\big]}{\Lambda\{\alpha(\mu+\rho)[\beta-\beta_1f(I)]+\beta_2[(\mu+\rho)(\mu+\eta+\delta)-\rho\delta]\}}. \end{aligned} \end{equation} (6.3)

    Obviously, \varphi(I) is increasing for I > 0 . Further

    \psi' (I) = \frac{h(I)}{\big\{[C(I^*)+\alpha(\mu+p)(\mu+\rho)][\sigma C(I^*)+\alpha\mu(\mu+\rho)]\big\}^2},

    where

    \begin{aligned} h(I) = &\alpha(\mu+\rho) C'(I)\bigg\{\sigma \big[C(I)+\alpha(\mu+\rho)(\mu+p)\big]\big[\sigma C(I)+\alpha\mu(\mu+\rho)\big]\\ &-\big[\sigma C(I)+\alpha(\mu+\rho)(\mu+\sigma p)\big]\big[\sigma C(I)+\alpha\mu(\mu+\rho)\big]\\ &-\big[\sigma C(I)+\alpha(\mu+\rho)(\mu+\sigma p)\big]\sigma \big[C(I)+\alpha(\mu+\rho)(\mu+p)\big]\bigg\} < 0. \end{aligned}

    Combined with (A_1) , (A_2) and (A_3) , we obtain that function \psi(I) is decreasing for I > 0 , and then H(I) is decreasing for I > 0 .

    In addition, from (6.3) we get

    \begin{equation} H(0) = \frac{(\alpha+\mu)\big[(\mu+\rho)(\mu+\eta+\delta)-\rho\delta\big]}{\Lambda\{\beta\alpha(\mu+\rho)+\beta_2[(\mu+\rho)(\mu+\eta+\delta)-\rho\delta]\}}(R_0-1), \end{equation} (6.4)

    and

    H(I) < \frac{\alpha(\mu+\rho)}{C(I)}-\frac{(\alpha+\mu)\big[(\mu+\rho)(\mu+\eta+\delta)-\rho\delta \big]}{\Lambda\{\alpha(\mu+\rho)[\beta-\beta_1f(I)]+\beta_2[(\mu+\rho)(\mu+\eta+\delta)-\rho\delta]\}}.

    It is easy to see that H(0) > 0 and \lim_{I\to \infty}H(I) < 0 when R_0 > 1 . Based on the monotonicity of H(I) , we can get system (6.1) has unique endemic equilibrium Q^* = (S^*, V^*, E^*, I^*, R^*) when R_0 > 1 . This completes the proof.

    Next we mainly discuss the global stability of disease-free equilibrium Q^0 and endemic equilibrium Q^* . Based on the discussion of Theorem 2, we give the following invariant domain

    \Gamma = \{\phi = (\phi_1, \phi_2, \phi_3, \phi_4, \phi_5)\in X_+: \phi_{1}(x)\leq S^0, \phi_{2}(x)\leq V^0, \;\; x\in\bar{\Omega}\}.

    Theorem 5. (a) If R_0 \le 1 , then disease-free equilibrium Q^0 is globally asymptotically stable in domain \Gamma .

    (b) If R_0 > 1 , then equilibrium Q^0 is unstable.

    Proof. Define

    W(t) = \int_{\Omega}^{}\bigg[S^0G\bigg(\frac{S}{S^0}\bigg)+V^0G\bigg(\frac{V}{V^0}\bigg)+E+\frac{\alpha+\mu}{\alpha}I+\frac{(\alpha+\mu)\rho}{\alpha(\mu+\rho)}R\bigg]dx,

    where G(x) = x-1-\ln x . According to [43, Theorem 3.7], we have

    \begin{equation} \begin{aligned} \frac{\mathrm{d} W(t)}{\mathrm{d} t} = & \int_{\Omega}^{}\bigg\{(1-\frac{S^0}{S})\big[\nabla \cdot(D_{1}(x)\nabla S)+\Lambda-[\beta-\beta_1f(I)]SI-\beta_2SE-(\mu+p)S\big]\\ &+(1-\frac{V^0}{V})\big[\nabla \cdot(D_{2}(x)\nabla V)+pS-\sigma[\beta-\beta_1f(I)]VI-\sigma\beta_2VE-\mu V\big]\\ &+ \nabla \cdot(D_{3}(x)\nabla E)+\big\{[\beta-\beta_1f(I)]I+\beta_2E\big\}(S+\sigma V)-(\alpha+\mu) E\\ &+\frac{\alpha+\mu}{\alpha}\big[\nabla \cdot(D_{4}(x)\nabla I)+\alpha E+\rho R-(\mu +\eta +\delta) I\big]\\ &+\frac{(\alpha+\mu)\rho}{\alpha(\mu+\rho)}\big[\nabla \cdot(D_{5}(x)\nabla R)+\delta I-(\mu +\rho) R\big]\bigg\}dx\\ \le & \int_{\Omega}^{}\bigg\{\mu S^0\bigg(2-\frac{S^0}{S}-\frac{S}{S^0}\bigg)+\mu V^0\bigg(3-\frac{S^0}{S}-\frac{V}{V^0}-\frac{SV^0}{S^0V}\bigg)\\ &+\big[\beta-\beta_1f(I)\big] S^0I+\beta_2S^0E+\sigma \big[\beta-\beta_1f(I)\big] V^0I+\sigma \beta_2V^0E-(\alpha+\mu) E\\ &+\frac{\alpha+\mu}{\alpha}\big[\alpha E+\rho R-(\mu +\eta +\delta) I\big]+\frac{(\alpha+\mu)\rho}{\alpha(\mu+\rho)}\big[\delta I-(\mu +\rho) R\big]\bigg\}dx\\ \le & \int_{\Omega}^{}\bigg[\mu S^0\bigg(2-\frac{S^0}{S}-\frac{S}{S^0}\bigg)+\mu V^0\bigg(3-\frac{S^0}{S}-\frac{V}{V^0}-\frac{SV^0}{S^0V}\bigg)\bigg]dx\\ &+\frac{1}{\alpha\mu(\mu+p)(\mu+\rho)}\int_{\Omega}^{}I\bigg\{\mu(\mu+p)(\alpha+\mu)\big[(\mu+\eta+\delta)(\mu+\rho)-\rho\delta\big](R_0-1)\bigg\}dx. \end{aligned} \end{equation} (6.5)

    Therefore, \frac{dW(t)}{dt}\leq 0 when R_0 < 1 and \frac{dW(t)}{dt} = 0 if and only if S = S^0 , V = V^0 and I(t, x) = 0 , from model (6.1), we can obtain E(t, x) = 0 and R(t, x) = 0 . This means that {Q^0} is the largest invariant set in \{(S, V, E, I, R)\in \Gamma: \frac{d W(t)}{dt} = 0\}. By the invariance principle, we can conclude that Q^0 is globally asymptotically stable for model (6.1).

    Theorem 6. If R_0 > 1 and

    \begin{equation} \begin{aligned} \big[g(I)-g(I^*)\big]\bigg[\frac{g(I)}{I}-\frac{g(I^*)}{I^*}\bigg]\leq 0, \;\;(E-E^*)\bigg(\frac{E}{I}-\frac{E^*}{I^*}\bigg)\leq 0, \end{aligned} \end{equation} (6.6)

    where g(I) = [\beta-\beta_1f(I)]I , then endemic equilibrium Q^* is globally asymptotically stable.

    Proof. Define

    L(t) = \int_{\Omega}^{}\big[W_1(t)+W_2(t)+W_3(t)+W_4(t)+W_5(t)\big]dx

    where

    \begin{aligned} W_1(t) & = S-S^*-S^*\ln\frac{S}{S^*}, \; W_2(t) = V-V^*-V^*\ln\frac{V}{V^*}, \; W_3(t) = E-E^*-E^*\ln\frac{E}{E^*}, \\ W_4(t) & = \frac{\alpha+\mu}{\alpha}(I-I^*-I^*\ln\frac{I}{I^*}), \; W_5(t) = \frac{(\alpha+\mu)\rho}{\alpha(\mu+\rho)}(R-R^*-R^*\ln\frac{R}{R^*}). \end{aligned}

    Define \beta(I) = \beta-\beta_1f(I) , we yield that

    \begin{equation} \begin{aligned} \frac{dL(t)}{dt} = &\int_{\Omega}^{}\bigg\{(1-\frac{S^*}{S})\frac{dS}{dt}+(1-\frac{V^*}{V})\frac{dV}{dt}+(1-\frac{E^*}{E})\frac{dE}{dt}\\ &+\frac{\alpha+\mu}{\alpha}(1-\frac{I^*}{I})\frac{dI}{dt}+\frac{(\alpha+\mu)\rho}{\alpha(\mu+\rho)}(1-\frac{R^*}{R})\frac{dR}{dt}\bigg\}dx\\ = &\int_{\Omega}^{}\bigg\{(1-\frac{S^*}{S})\big[\nabla \cdot(D_{1}(x)\nabla S)+\Lambda-\beta(I)SI-\beta_2SE-pS-\mu S\big]\\ &+(1-\frac{V^*}{V})\big[\nabla \cdot(D_{2}(x)\nabla V)+pS-\sigma\beta(I)VI-\sigma \beta_2VE-\mu V\big]\\ &+(1-\frac{E^*}{E})\big[\nabla \cdot(D_{3}(x)\nabla E)+\beta(I)SI+\beta_2SE+\sigma\beta(I)VI+\sigma\beta_2VE-\alpha E-\mu E\big]\\ &+\frac{\alpha+\mu}{\alpha}(1-\frac{I^*}{I})\big[\nabla \cdot(D_{4}(x)\nabla I)+\alpha E+\rho R-(\mu+\eta+\delta)I\big]\\ &+\frac{(\alpha+\mu)\rho}{\alpha(\mu+\rho)}(1-\frac{R^*}{R})\big[\nabla \cdot(D_{5}(x)\nabla R)+\delta I-(\mu+\rho)R\big]\bigg\}dx.\\ \end{aligned} \end{equation} (6.7)

    From [43, Theorem 3.7], we can deduce that

    \begin{align*} \label{system-4.1} \frac{dL(t)}{dt} = &\int_{\Omega}\bigg\{-S^*D_1(x)\frac{\left \| \nabla S \right \|^2}{S^2}-V^*D_2(x)\frac{\left \| \nabla V \right \|^2}{V^2}-E^*D_3(x)\frac{\left \| \nabla E\right \|^2}{E^2}\\ &-\frac{\alpha+\mu}{\alpha}I^*D_4(x)\frac{\left \| \nabla I \right \|^2}{I^2}-\frac{(\alpha+\mu)\rho}{\alpha(\mu+\rho)}R^*D_5(x)\frac{\left \| \nabla R \right \|^2}{R^2}\\ &+(1-\frac{S^*}{S})\big[\Lambda-\beta(I)SI-\beta_2SE-pS-\mu S\big]+(1-\frac{V^*}{V})\big[pS-\sigma\beta(I)VI-\sigma \beta_2VE-\mu V\big]\\ &+(1-\frac{E^*}{E})\big[\beta(I)SI+\beta_2SE+\sigma\beta(I)VI+\sigma\beta_2VE-\alpha E-\mu E\big]\\ &+(1-\frac{I^*}{I})\big[(\alpha+\mu) E+\frac{(\alpha+\mu)\rho}{\alpha} R-\frac{(\alpha+\mu)(\mu+\eta+\delta)}{\alpha}I\big]\\ &+(1-\frac{R^*}{R})\big[\frac{(\alpha+\mu)\rho\delta}{\alpha(\mu+\rho)} I-\frac{(\alpha+\mu)\rho}{\alpha}R\big]\bigg\}dx\\ \leq & \int_{\Omega}\bigg\{\mu S^*(2-\frac{S^*}{S}-\frac{S}{S^*})+\mu V^*(3-\frac{S^*}{S}-\frac{V}{V^*}-\frac{SV^*}{S^*V})\\ &+\big[\sigma\beta(I^*)V^* I^*+\sigma\beta_2V^*E^*\big](2-\frac{S^*}{S}-\frac{S}{S^*})+(1-\frac{S^*}{S})\big[\beta(I^*)S^*I^*+\beta_2V^*E^*\big]+\beta(I)S^*I\\ &+\beta_2S^*E+(\frac{S^*}{S}-\frac{SV^*}{S^*V})[\sigma\beta(I^*)V^*I^*+\sigma\beta_2V^*E^*]+\sigma\beta(I)V^*I+\sigma\beta_2V^*E-\beta(I)SI\cdot\frac{E^*}{E}\\ &-\beta_2SE^*-\sigma\beta(I)VI\cdot\frac{E^*}{E}-\sigma\beta_2VE^*+(\alpha+\mu)E^*-(\alpha+\mu)E^*\cdot\frac{I}{I^*}-(\alpha+\mu)E^*\cdot\frac{EI^*}{E^*I}\\ &+(\alpha+\mu)E^*+\frac{2(\alpha+\mu)\rho}{\alpha}R^*-\frac{(\alpha+\mu)\rho}{\alpha}R\cdot\frac{I^*}{I}-\frac{(\alpha+\mu)\rho\delta}{\alpha(\mu+\rho)}I\cdot\frac{R^*}{R}\bigg\}dx\\ = & \int_{\Omega}\bigg\{\mu S^*\bigg(2-\frac{S^*}{S}-\frac{S}{S^*}\bigg)+\mu V^*\bigg(3-\frac{S^*}{S}-\frac{V}{V^*}-\frac{SV^*}{S^*V}\bigg)\\ &+\beta(I^*)S^*I^*\bigg(3-\frac{S^*}{S}+\frac{\beta(I)I}{\beta(I^*)I^*}-\frac{S\beta(I)IE^*}{S^*\beta(I^*)I^*E}-\frac{I}{I^*}-\frac{I^*E}{IE^*}\bigg)\\ &+\beta_2S^*E^*\bigg(3-\frac{S^*}{S}-\frac{S}{S^*}+\frac{E}{E^*}-\frac{I}{I^*}-\frac{I^*E}{IE^*}\bigg)\\ &+\sigma\beta(I^*)V^*I^*\bigg(4-\frac{S^*}{S}-\frac{SV^*}{S^*V}+\frac{\beta(I)I}{\beta(I^*)I^*}-\frac{V\beta(I)IE^*}{V^*\beta(I^*)I^*E}-\frac{I}{I^*}-\frac{I^*E}{IE^*}\bigg)\\ &+\sigma\beta_2V^*E^*\bigg(4-\frac{S^*}{S}-\frac{SV^*}{S^*V}-\frac{V}{V^*}+\frac{E}{E^*}-\frac{I}{I^*}-\frac{I^*E}{IE^*}\bigg)\\ &+\frac{(\alpha+\mu)\rho\delta}{\alpha(\mu+\rho)}I^*\bigg(2-\frac{I^*R}{IR^*}-\frac{IR^*}{I^*R}\bigg)\bigg\}dx. \end{align*}

    Let g(I) = \beta(I)I and F(x) = 1-x+\ln x . Obviously, we have F(x)\le 0 for all x > 0 , and F(x) = 0 for x = 1 , then

    \begin{align*} &3-\frac{S^*}{S}+\frac{g(I)}{g(I^*)}-\frac{Sg(I)E^*}{S^*g(I^*)E}-\frac{I}{I^*}-\frac{I^*E}{IE^*}\\ = &3-\frac{S^*}{S}+\frac{g(I)}{g(I^*)}-\frac{Sg(I)E^*}{S^*g(I^*)E}-\frac{I}{I^*}-\frac{I^*E}{IE^*}+\ln \frac{S^*}{S}+\ln \frac{Sg(I)E^*}{S^*g(I^*)E}+\ln \frac{Ig(I^*)}{I^*g(I)}+\ln \frac{I^*E}{IE^*}\\ = &F(\frac{S^*}{S})+F(\frac{Sg(I)E^*}{S^*g(I^*)E})+F(\frac{Ig(I^*)}{I^*g(I)})+F(\frac{I^*E}{IE^*})+\frac{I}{g(I)g(I^*)}[g(I)-g(I^*)]\bigg[\frac{g(I)}{I}-\frac{g(I^*)}{I^*}\bigg], \\ &3-\frac{S^*}{S}-\frac{S}{S^*}+\frac{E}{E^*}-\frac{I}{I^*}-\frac{I^*E}{IE^*}\\ = &F(\frac{S^*}{S})+F(\frac{S}{S^*})+F(\frac{I^*E}{IE^*})+F(\frac{IE^*}{I^*E})+\frac{I}{EE^*}(E-E^*)\bigg(\frac{E}{I}-\frac{E^*}{I^*}\bigg), \\ & 4-\frac{S^*}{S}-\frac{SV^*}{S^*V}-\frac{V}{V^*}+\frac{E}{E^*}-\frac{I}{I^*}-\frac{I^*E}{IE^*}\\ = &F(\frac{S^*}{S})+F(\frac{SV^*}{S^*V})+F(\frac{V}{V^*})+F(\frac{I^*E}{IE^*})+F(\frac{IE^*}{I^*E})+\frac{I}{EE^*}(E-E^*)\bigg(\frac{E}{I}-\frac{E^*}{I^*}\bigg), \\ &4-\frac{S^*}{S}-\frac{SV^*}{S^*V}+\frac{g(I)}{g(I^*)}-\frac{Vg(I)E^*}{V^*g(I^*)E}-\frac{I}{I^*}-\frac{I^*E}{IE^*}\\ = &F(\frac{S^*}{S})+F(\frac{SV^*}{S^*V})+F(\frac{Vg(I)E^*}{V^*g(I^*)E})+F(\frac{Ig(I^*)}{I^*g(I)})+F(\frac{I^*E}{IE^*})+\frac{I}{g(I)g(I^*)}[g(I)-g(I^*)]\bigg[\frac{g(I)}{I}-\frac{g(I^*)}{I^*}\bigg]. \end{align*}

    Therefore, from (6.6) we further have \frac{dL(t)}{dt}\leq 0 for all S > 0 , V > 0 , E > 0 , I > 0 and R > 0 . Moreover, \frac{dL(t)}{dt} = 0 if and only if S = S^* , V = V^* , E = E^* , I = I^* and R = R^* , which means that \{Q^*\} is the largest invariant set for \frac{dL(t)}{dt} = 0 . Using invariance principle, we obtain that Q^* of system (6.1) is globally asymptotically stable.

    For convenience, we pay our attention to \Omega = [0, 10] . The nonlinear term under the influence of media reports is expressed as f(I) = \frac{q_1I}{1+c_1I} (where c_1 reflects the deviation degree of people's and relevant departments' understanding of disease information). At the same time, we take the initial value as follows:

    \begin{equation} \left\{ \begin{aligned} &S(x, 0) = 9700\times 0.92\times e^{-10(x-5)^2}, \; V(x, 0) = 9600\times 0.94\times e^{-10(x-5)^2}, \\ &E(x, 0) = 4000\times 0.04\times e^{-10(x-5)^2}, \; I(x, 0) = 40\times 0.02\times e^{-10(x-5)^2}, \; R(x, 0) = 0.\\ \end{aligned} \right. \end{equation} (7.1)

    Case 1: In order to simulate the result of Theorem 2, we take the parameters in the model (1.1) as shown in Table 2. By the numerical scheme in [47], we get the basic reproduction number R_0 = 0.9484 < 1 . Moreover, we plot the time evolution of exposed individuals E(x, t) (Figure 1(a)(c)) and infected individuals I(x, t) (Figure 1(d)(f)) for model (1.1). Figure 1 shows that the density of individuals E(x, t) and infected individuals I(x, t) converge to zero over time, which indicates the Theorem 2.

    Table 2.  Values of all parameters in model (1.1) for Case 1.
    Parameter Value Parameter Value Parameter Value
    \Lambda(x) 40(1+0.5\sin(2\pi x)) c_1(x) 0.15 q_1(x) 0.08
    \beta(x) 4\times10^{-4}(1+0.5\sin(2\pi x)) \sigma(x) 0.6 \alpha(x) 0.2(1+0.2\sin(2\pi x))
    \beta_1(x) 3.6\times10^{-4}(1+0.5\sin(2\pi x)) \delta(x) 0.3(1+0.5\sin(2\pi x)) \eta(x) 0.2(1+0.2\sin(2\pi x))
    \beta_2(x) 8\times10^{-4}(1+0.5\sin(2\pi x)) p(x) 0.3(1+0.5\sin(2\pi x)) \rho(x) 0.4(1+0.5\sin(2\pi x))
    \mu(x) 0.1(1+0.5\sin(2\pi x)) D_1(x) 0.09+0.005\sin(2\pi x) D_2(x) 0.08+0.005\sin(2\pi x)
    D_3(x) 0.07+0.005\sin(2\pi x) D_4(x) 0.06+0.005\sin(2\pi x) D_5(x) 0.085+0.005\sin(2\pi x)

     | Show Table
    DownLoad: CSV
    Figure 1.  Time evolution of system (1.1) with initial values. (a)–(c): The evolution of the E(x, t) over time; (d)–(f): The evolution of the I(x, t) over time.

    Case 2: To support the conclusion of Theorem 3, we choose \beta(x) = 8\times10^{-4}\times(1+0.5\sin(2\pi x)) , \beta_2(x) = 1\times10^{-3}\times(1+0.5\sin(2\pi x)) , \alpha(x) = 0.1\times(1+0.2\sin(2\pi x)) , \eta(x) = 0.08\times(1+0.2\sin(2\pi x)) and the other parameters remain the same as Table 2. By the numerical scheme in [47], we get the basic reproduction number R_0 = 1.8605 > 1 . Moreover, we plot the time evolution of exposed individuals E(x, t) (Figure 2(a)(c)) and infected individuals I(x, t) (Figure 2(d)(f)) for model (1.1). Indeed, as seen in Figure 2, the densities of exposed individuals E(x, t) and infected individuals I(x, t) converge to a spatially heterogeneous steady state over time, which shows Theorem 3.

    Figure 2.  Time evolution of system (1.1) with initial values. (a)–(c): The evolution of the E(x, t) over time; (d)–(f): The evolution of the I(x, t) over time.

    Case 3: To support the conclusion of Theorems 5 and 6, we first select the parameters as shown in Figure 3, and keep the diffusion coefficient as shown in Table 2. In this case, we calculate R_0 = 0.4808 < 1 . It can be seen from the Figure 3 that the density of exposed individuals E(x, t) and infected individuals I(x, t) tends to zero over time, which indicates Theorem 5. Then we choose \beta = 3.6\times10^{-3} , \beta_1 = 6\times10^{-4} , \beta_2 = 6\times10^{-3} and the other parameters remain the same as in Figure 3. In this case, we calculate R_0 = 4.8081 > 1 . As can be seen from the Figure 4, the density of exposed individuals E(x, t) and infected individuals I(x, t) tends to a steady state over time, which indicates Theorem 6.3.

    Figure 3.  (a)–(d): Evolution of the model (6.1) solutions with time when R_0 = 0.4808 . \Lambda = 10000/(2.79\times 365) , \mu = 0.03 , \beta = 3.6\times10^{-4} , \beta_1 = 6\times10^{-5} , \beta_2 = 6\times10^{-4} , p = 0.4 , \alpha = 0.6 , \eta = 0.1 , \rho = 0.3 , \delta = 0.16 , \sigma = 0.4 , q_1 = 0.3 , c_1 = 0.15 . The initial value is (200, 120, 40, 4, 0) .
    Figure 4.  (a)–(d): Evolution of the model (6.1) solutions with time when R_0 = 4.8081 . \Lambda = 10000/(2.79\times 365) , \mu = 0.03 , \beta = 3.6\times10^{-3} , \beta_1 = 6\times10^{-4} , \beta_2 = 6\times10^{-3} , p = 0.4 , \alpha = 0.6 , \eta = 0.1 , \rho = 0.3 , \delta = 0.16 , \sigma = 0.4 , q_1 = 0.3 , c_1 = 0.15 . The initial value is (200, 120, 40, 4, 0) .

    Basic reproduction number is an important index to assess the extinction or persistence of disease, so it is necessary to study the effect of spatial heterogeneity on R_0 . First, we choose \beta(x) = 4\times10^{-4}\times(1+c\times \cos(2\pi x)) and \beta_2(x) = 8.491\times10^{-4}\times(1+0.5\times \cos(2\pi x)) and keep the other parameters as in Table 2 to study the influence of spatial heterogeneity on R_0 , where c\in[0, 1] reflects the spatial heterogeneity intensity of the environment. Specifically, when c = 0 , the space is homogeneous, while as c increases, the spatial heterogeneity of the environment also increases. Second, we choose \beta(x) = 4\times10^{-4}\times(1+0.5\times \cos(2\pi x)) and \beta_2(x) = 8.49\times10^{-4}\times(1+c\times \cos(2\pi x)) to study the influence of spatial heterogeneity, and the other parameters are the same as in Table 2. It can be seen from Figure 5(a), (b) that with the increase of the intensity of spatial heterogeneity c , the basic reproduction number R_0 increases. Further, we can see that R_0 is minimum in the case of spatial homogeneity of c = 0 , and R_0 is greater than 1 when c exceeds the critical value. On the other hand, c_2 is defined as the spatial heterogeneity of the diffusion coefficient, we choose \beta(x) = 6\times10^{-4}\times(1+0.5\times \cos(2\pi x)) , \beta_2(x) = 9\times10^{-4}\times(1+0.5\times \cos(2\pi x)) and D_4(x) = 0.07\times(1+c_2\times\cos(2\pi x)) , and other parameters are the same as in Example 2. Figure 5(c) shows that the basic reproduction number R_0 decreases as the spatial heterogeneity c_2 of the diffusion coefficient increases, and R_0 is less than 1 when c_2 passes a critical value. In other words, in the case of spatial heterogeneity, diffusion coefficient of the change can change the final state of the disease.

    Figure 5.  Relationship between R_0 and spatial heterogeneity of parameters.

    Further, we study the effect of heterogeneity in treatment rate \delta(x) , recurrence rate \rho(x) and disease-related death rate \eta(x) on R_0 , because the difference of treatment level in different regions will affect the spread of the disease, and the recurrence of people in different regions will also affect the spread of the disease due to environmental heterogeneity. We use c_i(i = 3, 4, 5) to indicate the level of treatment in the area. Take \delta(x) = 0.3\times(1+c_3\times\sin(2\pi x)) , \rho(x) = 0.4\times(1+c_4\times\sin(2\pi x)) and \eta(x) = 0.2\times(1+c_5\times\sin(2\pi x)) , where 0\leq c_i\leq 1 ( i = 3, 4, 5 ), and keep the other parameters as in Table 2. In this case, we obtain Figure 5(d)(f), we get that R_0 is a decreasing function of treatment intensity, and R_0 is an increasing function of relapse intensity. Therefore, a higher treatment rate and a better treatment effect (a lower recurrence rate) can reduce the transmission capacity of the disease. This means that the spread of the disease can be effectively controlled by increasing medical facilities and improving treatment. In addition, we find that the basic reproduction number R_0 is not monotonic for \eta(x) . For some reason, R_0 may increase with the enhancement of spatial heterogeneous \eta(x) .

    In order to study the effect of media coverage on disease transmission, we assume that c_1 = 0.04 means that there are extensive media reports, and people have a full understanding of disease information and a strong awareness of prevention, c_1 = 500 means there is little media coverage and people are completely unaware of disease information. In the case of spatial heterogeneity, other parameters are taken to be the same as those in Table 2, and the changes of the density of infected individuals over time under different media coverage intensity are obtained when R_0 < 1 (see Figure 6), where Figure 6(a), (b) represented strong media coverage and Figure 6(c), (d) represented small media coverage intensity. By comparing Figure 6(b), (d), it can be seen that at the same time t and position x , the stronger the media coverage, the smaller the density of infected individuals I(x, t) , which means that in the case of spatial heterogeneity, increasing the intensity of media coverage can reduce the peak value of infected individuals. Similarly, when R_0 > 1 , the density of infected individuals changes over time under different media coverage intensity (see Figure 7). By comparing Figure 7(b), (d), it can be seen that at the same time t and position x , the stronger the media coverage, the smaller the density of infected individuals I(x, t) . This means that in the case of spatial heterogeneity, increasing the intensity of media coverage can reduce the density of infected individuals and reduce the spread of the disease.

    Figure 6.  The effect of different levels of media coverage on the number of I(x, t) when R_0 < 1 . (a)–(b): c_1 = 0.04 , (c)–(d): c_1 = 500 .
    Figure 7.  The effect of different levels of media coverage on the number of I(x, t) when R_0 > 1 . (a)–(b): c_1 = 0.04 , (c)–(d): c_1 = 500 .

    In the case of spatial homogeneity but heterogeneous diffusion, using the parameters in Figures 3 and 4, we fixed the position x and obtained the curve of the density I(x, t) of infected individuals changing with time t under different media reports. It can be observed from the Figure 8 that when R_0 < 1 , the number of latent and infected people tends to zero under different media coverage degrees. When c_1 is smaller (and more widely reported in the media), the time to peak the number of infected individuals, the peak size and the time to elimination of the disease are also smaller. This means that a reduction in effective exposure to media coverage can reduce the number of infections and accelerate the disappearance of infectious diseases.

    Figure 8.  The effect of different levels of media coverage on the number of E(x, t) and I(x, t) when R_0 < 1 .

    It can be observed from Figure 9 that when R_0 > 1 , the number of latent and infected people tend to be different stable values under different media coverage levels. When c_1 is smaller, which is more widely reported in the media, the peak and final numbers of the number of latent and infected individuals are also smaller. This means that increased media coverage can control the effective exposure rate of the population and thus reduce the spread of the disease.

    Figure 9.  The effect of different levels of media coverage on the number of E(x, t) and I(x, t) when R_0 > 1 .

    In this paper, we propose a reaction-diffusion SVEIR-I model with media coverage in spatially heterogeneous environment. First, we analysis the well-posedness of this model. Second, we define and prove that the basic reproduction number R_0 is a threshold parameter that determines the extinction and persistence of the disease in the heterogeneous case. Then, we analyze the dynamic behavior of the model when the space is homogeneous, obtain the existence and uniqueness of the endemic equilibrium, and the global asymptotic stability of the disease-free equilibrium and the endemic equilibrium is proved by establishing an appropriate Lyapunov function. Finally, some numerical simulation examples are given to illustrate our major results.

    Numerical simulation results show that spatial heterogeneity will increase the risk of disease transmission. First, we use the infection rate of spatial heterogeneity to study the impact of heterogeneity on R_0 . Here, we define c as the heterogeneity strength of the spatial infection rate \beta(x) and \beta_2(x) , and find that R_0 will increase with the increase of heterogeneity intensity c , R_0 is the minimum when c = 0 , and when c exceeds a certain critical value, R_0 is greater than 1 , the disease status would change. From a biological perspective, increasing the contact rate between people would make the disease spread faster, which causes disease break out. Second, we study the effect of diffusion coefficient on R_0 under spatial heterogeneity. Here, we define c_2 as the heterogeneity strength of the spatial diffusion coefficient D_4(x) . We find that R_0 will decrease with the increase of heterogeneity intensity c_2 , and R_0 reach its maximum when c_2 = 0 . When c_2 exceed a certain critical value, R_0 is less than 1 , and the threshold state of the disease would change. According to the above analysis, we conclude that the final state of the disease may be affected by spatial heterogeneity. In addition, we analyze the effect of the other parameters on R_0 , and find that some measures to prevent diseases, such as strengthening the construction of regional medical facilities or improving the medical level, could reduce the spread of diseases to a certain extent. Moreover, we focus on the impact of media coverage on disease transmission, including the case of spatially heterogeneous environments and spatially homogeneous environments. The results show that in both cases, when R_0 < 1 , strengthening media coverage can shorten the time of the peak of the number of infected individuals, the size of the peak and the time to eliminate the disease, and when R_0 > 1 , strengthening media coverage can reduce the size of the final scale of infected individuals. Therefore, as a non-pharmaceutical measure to improve information awareness, media coverage plays a very important role in the public response and implementation of disease control.

    The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.

    This research was supported by the Natural Science Foundation of Xinjiang Uygur Autonomous Region (2022D01C64), Natural Science Foundation of China (12201540), University Scientific Research Projects (XJEDU2021Y001) and the Doctoral Foundation of Xinjiang University (620320024).

    All authors declare no conflicts of interest in this paper.



    [1] D. L. Smith, J. Dushoff, F. E. McKenzie, The risk of a mosquito-borne infectionin a heterogeneous environment, PLoS Biol., 2 (2004), e368. https://doi.org/10.1371/journal.pbio.0020368 doi: 10.1371/journal.pbio.0020368
    [2] X. Wang, X. Q. Zhao, J. Wang, A cholera epidemic model in a spatiotemporally heterogeneous environment, J. Math. Anal. Appl., 468 (2018), 893–912. https://doi.org/10.1016/j.jmaa.2018.08.039 doi: 10.1016/j.jmaa.2018.08.039
    [3] Y. Cai, X. Lian, Z. Peng, W. Wang, Spatiotemporal transmission dynamics for influenza disease in a heterogenous environment, Nonlinear Anal. Real World Appl., 46 (2019), 178–194. https://doi.org/10.1016/j.nonrwa.2018.09.006 doi: 10.1016/j.nonrwa.2018.09.006
    [4] Y. Luo, Z. Teng, X. Q. Zhao, Transmission dynamics of a general temporal-spatial vector-host epidemic model with an application to the dengue fever in Guangdong, China, Discrete Contin. Dyn. Syst. Ser. B, 28 (2023), 134–169. https://doi.org/10.3934/dcdsb.2022069 doi: 10.3934/dcdsb.2022069
    [5] T. Zheng, L. Nie, H. Zhu, Y. Luo, Z. Teng, Role of seasonality and spatial heterogeneous in the transmission dynamics of avian influenza, Nonlinear Anal. Real World Appl., 67 (2022), 103567. https://doi.org/10.1016/j.nonrwa.2022.103567 doi: 10.1016/j.nonrwa.2022.103567
    [6] J. Wang, J. Wang, Analysis of a reaction-diffusion cholera model with distinct dispersal rates in the human population, J. Dyn. Differ. Equations, 33 (2021), 549–575. https://doi.org/10.1007/s10884-019-09820-8 doi: 10.1007/s10884-019-09820-8
    [7] J. Wang, F. Xie, T. Kuniya, Analysis of a reaction-diffusion cholera epidemic model in a spatially heterogeneous environment, Commun. Nonlinear Sci. Numerical Simul., 80 (2020), 104951. https://doi.org/10.1016/j.cnsns.2019.104951 doi: 10.1016/j.cnsns.2019.104951
    [8] C. Zhang, J. Gao, H. Sun, J. Wang, Dynamics of a reaction-diffusion SVIR model in a spatial heterogeneous environment, Phys. A Stat. Mech. Appl., 533 (2019), 122049. https://doi.org/10.1016/j.physa.2019.122049 doi: 10.1016/j.physa.2019.122049
    [9] Y. Luo, S. Tang, Z. Teng, L. Zhang, Global dynamics in a reaction-diffusion multi-group SIR epidemic model with nonlinear incidence, Nonlinear Anal. Real World Appl., 50 (2019), 365–385. https://doi.org/10.1016/j.nonrwa.2019.05.008 doi: 10.1016/j.nonrwa.2019.05.008
    [10] Y. Luo, L. Zhang, Z. Teng, et al., Analysis of a general multi-group reaction-diffusion epidemic model with nonlinear incidence and temporary acquired immunity, Math. Comput. Simul., 182 (2021), 428–455. https://doi.org/10.1016/j.matcom.2020.11.002 doi: 10.1016/j.matcom.2020.11.002
    [11] H. Zhao, Y. Shi, X. Zhang, Dynamic analysis of a malaria reaction-diffusion model with periodic delays and vector bias, Math. Biosci. Eng., 19 (2022), 2538–2574. https://doi.org/10.3934/mbe.2022117 doi: 10.3934/mbe.2022117
    [12] J. Wang, B. Dai, Qualitative analysis on a reaction-diffusion host-pathogen model with incubation period and nonlinear incidence rate, J. Math. Anal. Appl., 514 (2022), 126322. https://doi.org/10.1016/j.jmaa.2022.126322 doi: 10.1016/j.jmaa.2022.126322
    [13] A. Marzano, S. Gaia, V. Ghisetti, et al., Viral load at the time of liver transplantation and risk of hepatitis B virus recurrence, Liver Transplant., 11 (2005), 402–409. https://doi.org/10.1002/lt.20402 doi: 10.1002/lt.20402
    [14] A. Y. Kim, J. Schulze zur Wiesch, T. Kuntzen, J. Timm, D. E. Kaufmann, J. E. Duncan, Impaired hepatitis C virus-specific T cell responses and recurrent hepatitis C virus in HIV coinfection, PLoS Med., 3 (2006), e492. https://doi.org/10.1371/journal.pmed.0030492 doi: 10.1371/journal.pmed.0030492
    [15] M. L. Lambert, E. Hasker, A. Van Deun, D. Roberfroid, M. Boelaert, P. Van der Stuyft, Recurrence in tuberculosis: relapse or reinfection?, Lancet Infect. Dis., 3 (2003), 282–287. https://doi.org/10.1016/S1473-3099(03)00607-8 doi: 10.1016/S1473-3099(03)00607-8
    [16] D. W. Kimberlin, D. J. Rouse, Genital herpes, New Eng. J. Med., 350 (2004), 1970–1977. https://doi.org/10.1056/NEJMcp023065
    [17] J. Benedetti, L. Corey, R. Ashley, Recurrence rates in genital herpes after symptomatic first-episode infection, Annals Int. Med., 121 (1994), 847–854. https://doi.org/10.7326/0003-4819-121-11-199412010-00004 doi: 10.7326/0003-4819-121-11-199412010-00004
    [18] P. Van den Driessche, L. Wang, X. Zou, Modeling diseases with latency and relapse, Math. Biosci. Eng., 4 (2007), 205. https://doi.org/10.3934/mbe.2007.4.205 doi: 10.3934/mbe.2007.4.205
    [19] M. Ghosh, S. Olaniyi, O. S. Obabiyi, Mathematical analysis of reinfection and relapse in malaria dynamics, Appl. Math. Comput., 373 (2020), 125044. https://doi.org/10.1016/j.amc.2020.125044 doi: 10.1016/j.amc.2020.125044
    [20] S. Liu, S. Wang, L. Wang, Global dynamics of delay epidemic models with nonlinear incidence rate and relapse, Nonlinear Anal. Real World Appl., 12 (2011), 119–127. https://doi.org/10.1016/j.nonrwa.2010.06.001 doi: 10.1016/j.nonrwa.2010.06.001
    [21] C. Vargas-De-Leon, On the global stability of infectious diseases models with relapse, Abstraction Appl. Mag., 9 (2014).
    [22] Y. Chen, J. Li, S. Zou, Global dynamics of an epidemic model with relapse and nonlinear incidence, Math. Methods Appl. Sci., 42 (2019), 1283–1291. https://doi.org/10.1002/mma.5439 doi: 10.1002/mma.5439
    [23] A. Lahrouz, H. El Mahjour, A. Settati, A. Bernoussi, Dynamics and optimal control of a non-linear epidemic model with relapse and cure, Phys. A Stat.l Mech. Appl., 496 (2018), 299–317. https://doi.org/10.1016/j.physa.2018.01.007 doi: 10.1016/j.physa.2018.01.007
    [24] D. Tudor, A deterministic model for herpes infections in human and animal populations, Siam Rev., 32 (1990), 136–139. https://doi.org/10.1137/1032003 doi: 10.1137/1032003
    [25] T. K. Kar, S. K. Nandi, S. Jana, M. Mandal, Stability and bifurcation analysis of an epidemic model with the effect of media, Chaos Solitons Fractals, 120 (2019), 188–199. https://doi.org/10.1016/j.chaos.2019.01.025 doi: 10.1016/j.chaos.2019.01.025
    [26] J. Cui, Y. Sun, H. Zhu, The impact of media on the control of infectious diseases. J. Dyn. Differ. Equations, 20 (2008), 31–53. https://doi.org/10.1007/s10884-007-9075-0
    [27] D. K. Das, S. Khajanchi, T. K. Kar, The impact of the media awareness and optimal strategy on the prevalence of tuberculosis, Appl. Math. Comput., 366 (2020), 124732. https://doi.org/10.1016/j.amc.2019.124732 doi: 10.1016/j.amc.2019.124732
    [28] S. M. Salman, Memory and media coverage effect on an HIV/AIDS epidemic model with treatment, J. Comput. Appl. Math., 385 (2021), 113203. https://doi.org/10.1016/j.cam.2020.113203 doi: 10.1016/j.cam.2020.113203
    [29] X. Wang, D. Gao, J. Wang, Influence of human behavior on cholera dynamics, Math. Biosci., 267 (2015), 41–52. https://doi.org/10.1016/j.mbs.2015.06.009 doi: 10.1016/j.mbs.2015.06.009
    [30] L. Wang, Z. Liu, X. Zhang, Global dynamics for an age-structured epidemic model with media impact and incomplete vaccination, Nonlinear Anal. Real World Appl., 32 (2016), 136–158. https://doi.org/10.1016/j.nonrwa.2016.04.009 doi: 10.1016/j.nonrwa.2016.04.009
    [31] R. K. Rai, A. K. Misra, Y. Takeuchi, Modeling the impact of sanitation and awareness on the spread of infectious diseases, Math. Biosci. Eng., 16 (2019), 667–700. https://doi.org/10.3934/mbe.2019032 doi: 10.3934/mbe.2019032
    [32] P. Song, Y. Xiao, Analysis of a diffusive epidemic system with spatial heterogeneity and lag effect of media impact, J. Math. Biol., 85 (2022), 17. https://doi.org/10.1007/s00285-022-01780-w doi: 10.1007/s00285-022-01780-w
    [33] D. P. Oran, E. J. Topol, The proportion of SARS-CoV-2 infections that are asymptomatic: a systematic review, Annals Int. Med., 174 (2021), 655–662. https://doi.org/10.7326/M20-6976 doi: 10.7326/M20-6976
    [34] M. Day, Covid-19: identifying and isolating asymptomatic people helped eliminate virus in Italian village, BMJ British Med. J., 368 (2020).
    [35] L. Wang, Z. Liu, C. Guo, Y. Li, X. Zhang, New global dynamical results and application of several SVEIS epidemic models with temporary immunity, Appl. Math. Comput., 390 (2021), 125648. https://doi.org/10.1016/j.amc.2020.125648 doi: 10.1016/j.amc.2020.125648
    [36] S. Zhao, L. Stone, D. Gao, D. He, Modelling the large-scale yellow fever outbreak in Luanda, Angola, and the impact of vaccination, PLoS Neglected Trop. Dis., 12 (2018), e0006158. https://doi.org/10.1371/journal.pntd.0006158 doi: 10.1371/journal.pntd.0006158
    [37] C. C. Zhu, J. Zhu, X. L. Liu, Influence of spatial heterogeneous environment on long-term dynamics of a reaction-diffusion SVIR epidemic model with relapse, Math. Biosci. Eng., 16 (2019), 5897–5922. https://doi.org/10.3934/mbe.2019295 doi: 10.3934/mbe.2019295
    [38] H. L. Smith, Monotone Dynamical Systems: An Introduction to the Theory of Competitive and Cooperative Systems, in Mathematical Surveys And Monographs, Providence, RI, 1995.
    [39] Y. Lou, X. Q, Zhao, A reaction-diffusion malaria model with incubation period in the vector population, J. Math. Biol., 62 (2011), 543–568. https://doi.org/10.1007/s00285-010-0346-8 doi: 10.1007/s00285-010-0346-8
    [40] T. Zheng, Y. Luo, X. Zhou, L. Zhang, Z. Teng, Spatial dynamic analysis for COVID-19 epidemic model with diffusion and Beddington-DeAngelis type incidence, Commun. Pure Appl. Anal., 22 (2023), 365–396. https://doi.org/10.3934/cpaa.2021154 doi: 10.3934/cpaa.2021154
    [41] M. H. Protter, H. F. Weinberger, Maximum Principles in Differential Equations, Prentice Hall, Englewood Cliffs, 1967.
    [42] H. Amann, Fixed point equations and nonlinear eigenvalue problems in ordered Banach spaces, SIAM Rev., 18 (1976), 620–709. https://doi.org/10.1137/1018114 doi: 10.1137/1018114
    [43] J. Groeger, Divergence theorems and the supersphere, J. Geom. Phys., 77 (2014), 13–29. https://doi.org/10.1016/j.geomphys.2013.11.004 doi: 10.1016/j.geomphys.2013.11.004
    [44] J. Wu, Theory and applications of partial functional differential equations, Springer Science - Business Media, 1996.
    [45] X. Q. Zhao, Dynamical systems in population biology, Springer, 2003.
    [46] H. R. Thieme, Spectral bound and reproduction number for infinite-dimensional population structure and time heterogeneity, SIAM J. Appl. Math., 70 (2009), 29–48. https://doi.org/10.1137/080732870 doi: 10.1137/080732870
    [47] W. Wang, X. Q. Zhao, Basic reproduction numbers for reaction-diffusion epidemic models, SIAM J. Appl. Dyn. Syst., (2012), 1652–1673. https://doi.org/10.1137/120872942
    [48] O. Diekmann, J. A. P. Heesterbeek, J. A. J. Metz, On the definition and the computation of the basic reproduction ratio R_0 in the models for infectious disease in heterogeneous populations, J. Math. Biol., 28 (1990), 365–382. https://doi.org/10.1007/BF00178324 doi: 10.1007/BF00178324
    [49] P. V. Driessche, J. Watmough, Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission, Math. Biosci., 180 (2002), 29–48. https://doi.org/10.1016/S0025-5564(02)00108-6 doi: 10.1016/S0025-5564(02)00108-6
    [50] H. R. Thieme, Convergence results and a Poincare-Bendixson trichoyomy for asymptotically autonomous differential equations, J. Math. Biol., 30 (1992), 755–763. https://doi.org/10.1007/BF00173267 doi: 10.1007/BF00173267
    [51] H. L. Smith, X. Q. Zhao, Robust persistence for semidynamical systems, Nonlinear Anal., 47 (2001), 6169–6179. https://doi.org/10.1016/S0362-546X(01)00678-2 doi: 10.1016/S0362-546X(01)00678-2
    [52] P. Magal, X. Q. Zhao, Global attractors and steady states for uniformly persistent dynamical systems, SIAM J. Math. Anal., 37 (2005), 251–275. https://doi.org/10.1137/S0036141003439173 doi: 10.1137/S0036141003439173
  • This article has been cited by:

    1. Rujie Yang, Hong Qiu, A Valid Transport Related SVEIHR Stochastic Epidemic Model with Coverage and Time Delays, 2024, 23, 2224-2880, 815, 10.37394/23206.2024.23.84
    2. Yicheng Hao, Yantao Luo, Zhidong Teng, Role of limited medical resources in an epidemic model with media report and general birth rate, 2025, 10, 24680427, 522, 10.1016/j.idm.2025.01.001
  • Reader Comments
  • © 2023 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(2090) PDF downloads(202) Cited by(2)

Figures and Tables

Figures(9)  /  Tables(2)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog