Next Article in Journal
Complexity and Chaos Analysis for Two-Dimensional Discrete-Time Predator–Prey Leslie–Gower Model with Fractional Orders
Previous Article in Journal
Predicting Sit-to-Stand Body Adaptation Using a Simple Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dynamic Analysis of a Stochastic Delayed SEIRS Epidemic Model with Lévy Jumps and the Impact of Public Health Education

School of Mathematics and Statistics, Xinyang Normal University, Xinyang 464000, China
*
Author to whom correspondence should be addressed.
Axioms 2023, 12(6), 560; https://doi.org/10.3390/axioms12060560
Submission received: 8 May 2023 / Revised: 31 May 2023 / Accepted: 3 June 2023 / Published: 6 June 2023

Abstract

:
This study presents a class of the stochastic time-delayed susceptible–educated–infective–recovered–susceptible (SEIRS) epidemic model incorporating both public health education and Lévy jumps. We prove that the system has a unique global positive solution. We also provide derived conditions sufficient for both extinction and persistence in the mean. The verification of the findings and conclusions is performed through parameter sensitivity analysis and numerical simulations. This study concludes that public health education, stochastic noises, vaccination, increased disease recovery levels, and reduced patient contact significantly contribute significantly to disease prevention and control.

1. Introduction

Infectious diseases have continually posed a significant threat to both human health and societal development. Thousands of fatalities are attributable to these diseases every year, hence the critical importance of studying their prevention and control mechanisms.
For disease control and prevention, the primary approach is proactive prevention, which underscores the importance of public health education. The latter enhances public awareness of disease prevention and control, thereby mitigating the prevalence of infectious diseases. The recent surge in using new media to disseminate epidemic prevention knowledge has proven efficient. This medium offers several benefits, including high dissemination efficiency, rapid transmission speed, and broad public acceptance [1,2,3,4]. By promoting infectious disease prevention knowledge among residents, the susceptibility of individuals can be mitigated, thereby reducing the risk of infectious diseases and achieving the objective of curbing disease spread. In 2018, Nyang’inja et al. [5] introduced the educated compartment E ( t ) , representing individuals educated about prevention strategies for tungiasis infection, and formulated a mathematical model for the dynamics of jigger infestation incorporating public health education using the systems of ordinary differential equations. Their findings affirmed the effectiveness of public health education as a control measure for eradicating jigger infestation in endemic communities.
Within the context of disease prevention and control, vaccinating susceptible populations has proven to be effective and plays a crucial role in curtailing the spread of epidemics. Vaccinations have been instrumental in controlling diseases such as measles and varicella, serving as a critical preventive measure for infected individuals [6,7,8]. Furthermore, temporary immunity, occurring when an infected individual recovers or a susceptible one gets vaccinated, can also be a significant factor in disease spread. However, this immunity may fade over time, returning the individual to a susceptible state. Many researchers have introduced time delays into their epidemic models to account for this temporary immunity phenomenon [9,10,11,12,13,14,15]. References [9,10,11,12] addressed the temporary immunity gained post-recovery. Kyrychko et al. [11] conducted the analysis of a delayed susceptible–infective–recovered (SIR) model, incorporating temporary immunity, and illustrated the impact of the immunity period—on the long-term dynamics of solutions. In 2010, Xu and Ma [12] researched a delayed susceptible-infective–recovered–susceptible (SIRS) epidemic model, which included saturation incidence and temporary immunity. They observed a correlation between the immunity period—and the amplitude of oscillations near the endemic equilibrium. Li et al. [13] evaluated the temporary vaccination immunity in their epidemic models and provided conditions under which the disease would become extinct. Two cases of temporary immunity, recovery immunity and vaccination immunity, were also considered in the literature [14,15,16]. In [16], Liu and Chen analyzed a deterministic delay SIR epidemic model with vaccination and temporary immunity, suggesting that time-delayed models better mirror actual scenarios. The above study of the literature on temporary immunity shows that considering the effects of time-delay is essential for predicting the extinction and persistence of infectious diseases.
In real-world scenarios, the dynamics of infectious disease models may be impacted by environmental noise during disease transmission, leading to sudden and significant deviations. This has led to the introduction of perturbation into deterministic models. Thus, numerous researchers have explored stochastic infectious disease models [17,18,19]. For instance, the authors of [17] examined a stochastic delayed SIRS epidemic model with seasonal variation, defining the system’s stochastic threshold and observing that the periodic solution’s oscillation intensity is dependent on the noise intensity. In [18], a stochastic cholera model with a saturated recovery rate was discussed, showing that high levels of noise could result in disease extinction. Major environmental disturbances such as floods, earthquakes, volcanoes, tsunamis, and hurricanes suggest the need to introduce Lévy jumps into infectious disease models to account for such abrupt fluctuations [20,21,22,23,24,25]. Zhang Xuegui et al. [20] studied a stochastic two predator–one prey system with Lévy jumps and mixed functional responses, observing that Lévy jumps negatively impacted species’ existence. Liu and Chen [23] proposed a general stochastic non-autonomous logistic model with delays and Lévy jumps, establishing sufficient criteria for extinction, non-persistence in the mean, and the weak persistence of the model. Their findings underscored the significant relationship between time-dependent delay, Lévy noise, and both persistence and extinction.
Building upon the foundation established by previous works, this study seeks to explore the dynamical properties of an innovative stochastic SEIRS model. The objective of this paper is to delve deeper into the factors influencing the spread of the infectious disease. We consider factors such as the effect of public health education, recovery and vaccination rates, contact rate, and time-delay, to analyze their impact on disease persistence and extinction. In particular, we will discuss the effects of discontinuous stochastic interference and its intensity on the average extinction and persistence of diseases.
This paper is structured as follows: Section 2 details the establishment of the stochastic SEIRS model and provides interpretations of relevant parameters. Section 3 demonstrates that the proposed system has a unique global positive solution. In Section 4 and Section 5, we established the sufficient conditions for extinction and persistence in the mean of the system. Section 6 discusses numerical simulations and sensitivity analysis. We conclude the paper with a summary of key findings.

2. Modeling and Preparation

In this paper, we propose a stochastic time-delayed SEIRS epidemic model incorporating Lévy jumps and the influence of public health education. The new system can be written as
d S ( t ) = [ Λ β S ( t ) I ( t ) 1 + p 1 I ( t ) + ρ S ( t τ 1 ) e μ τ 1 + γ I ( t τ 2 ) e μ τ 2 ρ S ( t ) ϵ S ( t ) μ S ( t ) ] d t + σ 1 S ( t ) d B 1 ( t ) + Z γ 1 ( u ) S ( t ) N ˜ ( d t , d u ) , d E ( t ) = ϵ S ( t ) α β E ( t ) I ( t ) 1 + p 2 I ( t ) μ E ( t ) d t + σ 2 E ( t ) d B 2 ( t ) + Z γ 2 ( u ) E ( t ) N ˜ ( d t , d u ) , d I ( t ) = β S ( t ) I ( t ) 1 + p 1 I ( t ) + α β E ( t ) I ( t ) 1 + p 2 I ( t ) ( γ + σ + μ ) I ( t ) d t + σ 3 I ( t ) d B 3 ( t ) + Z γ 3 ( u ) I ( t ) N ˜ ( d t , d u ) , d R ( t ) = [ γ I ( t ) + ρ S ( t ) ρ S ( t τ 1 ) e μ τ 1 γ I ( t τ 2 ) e μ τ 2 μ R ( t ) ] d t + σ 4 R ( t ) d B 4 ( t ) + Z γ 4 ( u ) R ( t ) N ˜ ( d t , d u ) ,
with the initial conditions
S ( θ ) = ϕ 1 ( θ ) 0 , E ( θ ) = ϕ 2 ( θ ) 0 , I ( θ ) = ϕ 3 ( θ ) 0 , R ( θ ) = ϕ 4 ( θ ) 0 , ϕ i ( 0 ) > 0 , i = 1 , 2 , 3 , 4 , θ [ τ , 0 ] ,
where ( ϕ 1 ( θ ) , ϕ 2 ( θ ) , ϕ 3 ( θ ) , ϕ 4 ( θ ) ) L 1 ( [ τ , 0 ] ; R + 4 ) , L 1 ( [ τ , 0 ] ; R + 4 ) denotes the family of Lebesgue integrable functions from [ τ , 0 ] into R + 4 and τ = max { τ 1 , τ 2 } .
Model (1) includes S, E, I, R, representing susceptible, educated, infected, and recovered populations, respectively. The educated population refers to those individuals informed about specific disease prevention strategies. The parameters are defined as follows: Λ signifies the recruitment rate of the total population; β S ( t ) I ( t ) / ( 1 + p 1 I ( t ) ) represents the effective contact rate, with β denoting the maximal effective contact rate between susceptible and infected individuals; ρ is the efficient vaccination rate for susceptible individuals; γ symbolizes the recovery rate; μ is the natural death rate; σ indicates the disease-related death rate within the compartment I; ϵ is the proportion of the educated population; α reflects the impact of public health education on infection reduction, with 0 < α < 1 ; and p 1 , p 2 , τ 1 , and τ 2 are all positive constants, representing the coefficient of saturated incidence rate and the length of temporary immunity period, respectively. B i ( t ) ( i = 1 , 2 , 3 , 4 ) is a standard four-dimensional Brownian motion defined on a complete probability space ( Ω , F , P ) with a filtration { F t } t R + satisfying the usual conditions ( { F t } t R + is right continuous and increasing while F 0 contains all P -null sets), and satisfies B i ( 0 ) = 0 ( i = 1 , 2 , 3 , 4 ) . σ i ( i = 1 , 2 , 3 , 4 ) is the corresponding intensity of white noise. Let N ˜ ( d t , d u ) = N ( d t , d u ) λ ( d u ) d t , N be a Poisson counting measure with a compensator N ˜ and characteristic measure λ on a measurable subset Z of ( 0 , + ) , which satisfies λ ( Z ) < + . S ( t ) , E ( t ) , I ( t ) , R ( t ) are the left limits of ( S ( t ) , E ( t ) , I ( t ) , R ( t ) ) , respectively. γ i ( u ) ( i = 1 , 2 , 3 , 4 ) denotes the intensity of L e ´ vy jumps.
In order to understand the meaning of solutions to stochastic differential system (1) and the Itô formula with L e ´ vy jumps, we present the following definition and lemmas.
Definition 1
([26]). Consider the following stochastic differential equation (SDE):
d x ( t ) = f ( x ( t ) , t ) d t + g ( x ( t ) , t ) d W ( t ) , t 0 t T
with the initial condition x ( t 0 ) = x 0 R n . An R n -valued stochastic process { x ( t ) } t 0 t T is said to be the solution of SDE (3) if:
(i) 
{ x ( t ) } t 0 t T is continuous and { F t } -adapted;
(ii) 
{ f ( x ( t ) , t ) } t 0 t T L 1 ( [ t 0 , T ] ; R n ) , and { g ( x ( t ) , t ) } t 0 t T L 2 ( [ t 0 , T ] ; R n × m ) ;
(iii) 
x ( t ) = x ( t 0 ) + t 0 T f ( x ( s ) , s ) d s + t 0 T g ( x ( s ) , s ) d W ( s ) , for almost all t [ t 0 , T ] with probability one.
Lemma 1
(The L e ´ vy–Itô decomposition theorem [27]). If X is a L e ´ vy process, then there exists b R d , a Brownian motion B and independent Poisson random measure N which is independent of B, such that for each t 0 ,
X ( t ) = b t + B ( t ) + | x | 1 x N ( t , d x ) + | x | < 1 x N ˜ ( t , d x ) ,
where
b = E X ( 1 ) | x | 1 x N ( 1 , d x ) .
Lemma 2
(Itô’s formula with L e ´ vy jumps [28]). Let X ( t ) R n be a solution of the following stochastic differential equation with L e ´ vy jumps:
d X ( t ) = F ( X ( t ) , t ) d t + G ( X ( t ) , t ) d B ( t ) + Z H ( X ( t ) , t , u ) N ˜ ( d t , d u ) ,
where F : R n × R + R n , G : R n × R + R n and H : R n × R + × Z R n are measurable functions.
Given V C 2 , 1 ( R n × R + ; R + ) , we define the operator L V : R n × R + R
L V ( X , t ) = V t ( X , t ) + V X ( X , t ) F ( X , t ) + 1 2 t r a c e [ G ( t ) V X X ( X , t ) G ( X , t ) ] + Z { V ( X + H ( X , t ) , t ) V ( X , t ) V X ( X , t ) H ( X , t ) } λ ( d u ) ,
where
V t ( X , t ) = V X ( X , t ) t ,
V X ( X , t ) = V X ( X , t ) X 1 , V X ( X , t ) X 2 , , V X ( X , t ) X n
V X X ( X , t ) = 2 V X ( X , t ) X i X j n × n .
Then, the generalized Itô formula with L e ´ vy jumps is given by
d V ( X , t ) = L V ( X , t ) d t + V X ( X , t ) G ( X , t ) d B ( t ) + Z { V ( X + H ( X , t ) , t ) V ( X , t ) } N ˜ ( d t , d u ) .
Notice that the first three equations in model (1) do not depend on the last equation, so it can be omitted without loss of generality. Thus, in the following text, we only discuss the simplified model
d S ( t ) = [ Λ β S ( t ) I ( t ) 1 + p 1 I ( t ) + ρ S ( t τ 1 ) e μ τ 1 + γ I ( t τ 2 ) e μ τ 2 ρ S ( t ) ϵ S ( t ) μ S ( t ) ] d t + σ 1 S ( t ) d B 1 ( t ) + Z γ 1 ( u ) S ( t ) N ˜ ( d t , d u ) , d E ( t ) = ϵ S ( t ) α β E ( t ) I ( t ) 1 + p 2 I ( t ) μ E ( t ) d t + σ 2 E ( t ) d B 2 ( t ) + Z γ 2 ( u ) E ( t ) N ˜ ( d t , d u ) , d I ( t ) = β S ( t ) I ( t ) 1 + p 1 I ( t ) + α β E ( t ) I ( t ) 1 + p 2 I ( t ) ( γ + σ + μ ) I ( t ) d t + σ 3 I ( t ) d B 3 ( t ) + Z γ 3 ( u ) I ( t ) N ˜ ( d t , d u ) ,
with the initial value
S ( θ ) = ϕ 1 ( θ ) 0 , E ( θ ) = ϕ 2 ( θ ) 0 , I ( θ ) = ϕ 3 ( θ ) 0 , ϕ 1 ( 0 ) > 0 , ϕ 2 ( 0 ) > 0 , ϕ 3 ( 0 ) > 0 , θ [ τ , 0 ] ,
where ( ϕ 1 ( θ ) , ϕ 2 ( θ ) , ϕ 3 ( θ ) ) L 1 ( [ τ , 0 ] ; R + 3 ) , L 1 ( [ τ , 0 ] ; R + 3 ) denotes the family of Lebesgue integrable functions from [ τ , 0 ] into R + 3 and τ = max { τ 1 , τ 2 } . For the convenience of the subsequent proof, we make the following assumption in this paper.
Assumption 1.
γ ( u ) is a bounded function, γ ( u ) > 1 and | Λ u γ ( u ) | δ , u Z .

3. Existence and Uniqueness of the Positive Solution

In this section, to study the dynamical properties of an epidemic system, we shall prove that the solutions of the system (4) are positive and global.
Theorem 1.
Suppose that Assumption 1 holds, then for any given initial condition (5), the system (4) has a unique solution ( S ( t ) , E ( t ) , I ( t ) ) on t τ , and the solution will remain in R + 3 with a probability of one.
Proof. 
Obviously, since the coefficient of system (4) satisfies the local Lipschitz condition, and then for any given initial condition (5), there is a unique local solution ( S ( t ) , E ( t ) , I ( t ) ) of system (4) on t [ τ , τ e ) , where τ e denotes the explosion time. Now, we prove that the solution is global, that is τ e = a.s. Assume that k 0 1 be sufficiently large such that ( S ( 0 ) , E ( 0 ) , I ( 0 ) ) all lie in 1 k 0 , k 0 . For each integer k k 0 , we define the stopping time
τ k = inf t [ τ , τ e ) : min { S ( t ) , E ( t ) , I ( t ) } 1 k or max { S ( t ) , E ( t ) , I ( t ) } k .
As usual, we define inf = (represents the empty set). Clearly, τ k is strictly increasing when k . Let τ = lim k τ k , thus τ τ e a.s. That is to say, if we show that τ e = a.s., then τ e = and the solution ( S ( t ) , E ( t ) , I ( t ) ) R + 3 for all t τ . Next, we start the proof by contradiction, if this assertion is false, then there exists a pair of constants T > 0 and ϵ ( 0 , 1 ) such that P { τ T } > ϵ . Therefore, there is an integer k 1 k 0 such that for all k k 1 ,
P { τ k T } ϵ .
Define the C 2 -function V: R + 3 R + as follows:
V ( S , E , I ) = S b b ln S b + E 1 ln E + I 1 ln I + ρ e μ τ 1 t τ 1 t S ( s ) d s + γ e μ τ 2 t τ 2 t I ( s ) d s ,
where b is a positive constant to be determined later. Furthermore, the non-negativity of the function can be clearly derived from u 1 ln u 0 for all u > 0 . Let k > k 0 and T > 0 be arbitrary positive constants. Then, by the Itô’s formula, we have
d V = L V d t + σ 1 ( S b ) d B 1 ( t ) + σ 2 ( E 1 ) d B 2 ( t ) + σ 3 ( I 1 ) d B 3 ( t ) + Z [ S γ 1 ( u ) b ln ( 1 + γ 1 ( u ) ) ] N ˜ ( d t , d u ) + Z [ E γ 2 ( u ) ln ( 1 + γ 2 ( u ) ) ] N ˜ ( d t , d u ) + Z [ I γ 3 ( u ) ln ( 1 + γ 3 ( u ) ) ] N ˜ ( d t , d u ) ,
where L V : R + 3 R + is defined by
L V = 1 b S Λ β S I 1 + p 1 I + ρ S ( t τ 1 ) e μ τ 1 + γ I ( t τ 2 ) e μ τ 2 ( ϵ + μ ) S ρ S + 1 2 b σ 1 2 + 1 1 E ϵ S α β E I 1 + p 2 I μ E + 1 2 σ 2 2 + 1 1 I β S I 1 + p 1 I + α β E I 1 + p 2 I ( γ + σ + μ ) I + 1 2 σ 3 2 + Z { [ b γ 1 ( u ) b ln ( 1 + γ 1 ( u ) ) ] + [ γ 2 ( u ) ln ( 1 + γ 2 ( u ) ) ] + [ γ 3 ( u ) ln ( 1 + γ 3 ( u ) ) ] } λ ( d u ) + ρ e μ τ 1 S ( t ) ρ e μ τ 1 S ( t τ 1 ) + γ e μ τ 2 I ( t ) γ e μ τ 2 I ( t τ 2 ) = Λ μ ( S + E + I ) ( σ + γ ) I ρ S b S Λ + b β I 1 + p 1 I ρ b S S ( t τ 1 ) e μ τ 1 γ b S I ( t τ 2 ) e μ τ 2 + b ( ϵ + μ ) + b ρ β S 1 + p 1 I α β E 1 + p 2 I + γ + σ + μ ϵ S E + α β I 1 + p 2 I + μ + ρ S ( t ) e μ τ 1 + γ I ( t ) e μ τ 2 + 1 2 b σ 1 2 + 1 2 σ 2 2 + 1 2 σ 3 2 + Z { [ b γ 1 ( u ) b ln ( 1 + γ 1 ( u ) ) ] + [ γ 2 ( u ) ln ( 1 + γ 2 ( u ) ) ] + [ γ 3 ( u ) ln ( 1 + γ 3 ( u ) ) ] } λ ( d u ) Λ + b ( ϵ + ρ ) + γ + σ + μ ( 2 + b ) + α β p 2 + ( b β γ ( 1 e μ τ 2 ) μ σ ) I ( μ + ρ ( ( 1 e μ τ 1 ) ) S + 1 2 b σ 1 2 + 1 2 σ 2 2 + 1 2 σ 3 2 + 3 m ,
where m = max { Z [ b γ 1 ( u ) b ln ( 1 + γ 1 ( u ) ) ] λ ( d u ) , Z [ γ 2 ( u ) ln ( 1 + γ 2 ( u ) ) ] λ ( d u ) , Z [ γ 3 ( u ) ln ( 1 + γ 3 ( u ) ) ] λ ( d u ) } 0 . Choosing a positive constant b = [ γ ( 1 e μ τ 2 ) + μ + σ ] / β 0 such that b β γ ( 1 e μ τ 2 ) μ σ = 0 , then we have
L V Λ + b ( ϵ + ρ ) + γ + σ + μ ( 2 + b ) + α β p 2 + 1 2 b σ 1 2 + 1 2 σ 2 2 + 1 2 σ 3 2 + 3 m : = K ,
where K > 0 is a constant. Therefore, we obtain
d V K d t + σ 1 ( S b ) d B 1 ( t ) + σ 2 ( E 1 ) d B 2 ( t ) + σ 3 ( I 1 ) d B 3 ( t ) + Z [ S γ 1 ( u ) b ln ( 1 + γ 1 ( u ) ) ] N ˜ ( d t , d u ) + Z [ E γ 2 ( u ) ln ( 1 + γ 2 ( u ) ) ] N ˜ ( d t , d u ) + Z [ I γ 3 ( u ) ln ( 1 + γ 3 ( u ) ) ] N ˜ ( d t , d u ) .
Integrating from 0 to τ k T on both sides of the above inequality yields
0 τ k T d V 0 τ k T K d t + 0 τ k T σ 1 ( S b ) d B 1 ( t ) + 0 τ k T σ 2 ( E 1 ) d B 2 ( t ) + 0 τ k T σ 3 ( I 1 ) d B 3 ( t ) + 0 τ k T Z [ S γ 1 ( u ) b ln ( 1 + γ 1 ( u ) ) ] N ˜ ( d t , d u ) + 0 τ k T Z [ E γ 2 ( u ) ln ( 1 + γ 2 ( u ) ) ] N ˜ ( d t , d u ) + 0 τ k T Z [ I γ 3 ( u ) ln ( 1 + γ 3 ( u ) ) ] N ˜ ( d t , d u ) ,
where τ k T = min { τ k , T } . Consequently, taking the expectation leads to
E [ V ( S ( τ k T ) , E ( τ k T ) , I ( τ k T ) ) ] V ( S ( 0 ) , E ( 0 ) , I ( 0 ) ) + K E [ τ k T ] V ( S ( 0 ) , E ( 0 ) , I ( 0 ) ) + K T .
Let Ω k = { τ k T } , for k k 1 , we can know P ( Ω k ) ϵ . Note that, for each ω Ω k , there exists S ( τ k , ω ) , E ( τ k , ω ) or I ( τ k , ω ) equals either k or 1 / k , then
V ( S ( τ k T ) , E ( τ k T ) , I ( τ k T ) ) ( k 1 ln k ) 1 k 1 + ln k .
Hence,
V ( S ( 0 ) , E ( 0 ) , I ( 0 ) ) + K T E [ I Ω k ( ω ) V ( S ( τ k , ω ) , E ( τ k , ω ) , I ( τ k , ω ) ) ] ϵ ( k 1 ln k ) 1 k 1 + ln k ,
where I Ω k ( ω ) is the indicator function of Ω k . Letting k , it implies that
+ > V ( S ( 0 ) , E ( 0 ) , I ( 0 ) ) + K T = + ,
which is a contradiction. Therefore, we obtain τ = + a.s. The proof is completed. □
Obviously, there is also a unique global positive solution on t τ for system (1) a.s.

4. The Extinction of the Diseases

In this section, we shall discuss the condition of extinction of the diseases. Furthermore, for the sake of convenience, we introduce the following notation: X ( t ) = 1 t 0 t X ( s ) d s , where X ( t ) is an integrable function on [ 0 , + ) . Then, the following lemma and theorem are obtained.
Lemma 3
([29]). Let M = { M t } t 0 be a real-valued continuous local martingale vanishing at t = 0 , then
lim t M , M t = a . s . lim t M t M , M t = 0 a . s . ,
and
lim sup t M , M t t < a . s . lim t M t t = 0 a . s . ,
where M , M t is the quadratic variation of M.
Lemma 4.
Let ( S ( t ) , E ( t ) , I ( t ) ) be the solution of system (4) with any initial condition (5), then
lim t 0 t S ( s ) d B 1 ( s ) t = 0 , lim t 0 t E ( s ) d B 2 ( s ) t = 0 , lim t 0 t I ( s ) d B 3 ( s ) t = 0 a . s .
lim t 0 t Z γ 1 ( u ) S ( s ) N ˜ ( d s , d u ) t = 0 , lim t 0 t Z γ 2 ( u ) E ( s ) N ˜ ( d s , d u ) t = 0 ,
lim t 0 t Z γ 3 ( u ) I ( s ) N ˜ ( d s , d u ) t = 0 a . s .
Lemma 5.
Suppose that ( S ( t ) , E ( t ) , I ( t ) ) is the solution of model (4) with an initial condition (5), then
lim t S ( t ) + I ( t ) + ρ e μ t t τ 1 t e μ s S ( s ) d s + γ e μ t t τ 2 t e μ s I ( s ) d s t = 0 a . s .
Furthermore
lim t S ( t ) t = 0 , lim t I ( t ) t = 0 , lim t e μ t t τ 1 t e μ s S ( s ) d s t = 0 , lim t e μ t t τ 2 t e μ s I ( s ) d s t = 0 a . s .
The proof is similar to that of Lemma 1 in [30], and is omitted in this paper.
Theorem 2.
Suppose that Assumption 1 holds, and ( S ( t ) , E ( t ) , I ( t ) ) can be any solution for system (4) with an initial condition (5), if the following condition
R 1 = Λ μ · β ( 1 + α ) γ + σ + μ + θ 3 = β Λ μ ( γ + σ + μ + θ 3 ) + α β Λ μ ( γ + σ + μ + θ 3 ) < 1
holds, where θ 3 = 1 2 σ 3 2 + Z [ γ 3 ln ( 1 + γ 3 ) ] λ ( d u ) 0 , then the infectious disease of system (4) goes to extinction. That is to say, lim t I ( t ) = 0 a.s.
Moreover
lim t S ( t ) = Λ ϵ + μ + ρ ( 1 e μ τ 1 ) , lim t E ( t ) = ϵ Λ μ ( ϵ + μ + ρ ( 1 e μ τ 1 ) ) a . s .
Proof. 
Applying Itô’s formula to the second equation of system (4), we derive that
d ln I ( t ) = β S 1 + p 1 I + α β E 1 + p 2 I ( γ + σ + μ ) 1 2 σ 3 2 + Z [ ln ( 1 + γ 3 ) γ 3 ] λ ( d u ) d t + σ 3 d B 3 ( t ) + Z [ ln ( 1 + γ 3 ) ] N ˜ ( d t , d u ) = β S 1 + p 1 I + α β E 1 + p 2 I ( γ + σ + μ ) θ 3 d t + σ 3 d B 3 ( t ) + Z [ ln ( 1 + γ 3 ) ] N ˜ ( d t , d u ) ,
where θ 3 = 1 2 σ 3 2 + Z [ γ 3 ln ( 1 + γ 3 ) ] λ ( d u ) . Then, integrating (8) from 0 to t and dividing by t on both sides gives
ln I ( t ) t ln I ( 0 ) t + β S ( t ) + α β E ( t ) ( γ + σ + μ + θ 3 ) + 1 t 0 t σ 3 d B 3 ( s ) + 1 t 0 t Z [ ln ( 1 + γ 3 ) ] N ˜ ( d s , d u ) .
Note that
d S ( t ) + E ( t ) + I ( t ) + ρ e μ t t τ 1 t e μ s S ( s ) d s + γ e μ t t τ 2 t e μ s I ( s ) d s = [ Λ + ρ S ( t τ 1 ) e μ τ 1 + γ I ( t τ 2 ) e μ τ 2 μ ( S + E + I ) ( σ + γ ) I ρ S ] d t + σ 1 S ( t ) d B 1 ( t ) + σ 2 E ( t ) d B 2 ( t ) + σ 3 I ( t ) d B 3 ( t ) + Z [ γ 1 ( u ) S ( t ) + γ 2 ( u ) E ( t ) + γ 3 ( u ) I ( t ) ] N ˜ ( d t , d u ) + μ ρ e μ t t τ 1 t e μ s S ( s ) d s + ρ e μ t · e μ t S ( t ) ρ e μ t · e μ ( t τ 1 ) S ( t τ 1 ) d t + μ γ e μ t t τ 2 t e μ s I ( s ) d s + γ e μ t · e μ t I ( t ) γ e μ t · e μ ( t τ 2 ) I ( t τ 2 ) d t = Λ μ S + E + I + ρ e μ t t τ 1 t e μ s S ( s ) d s + γ e μ t t τ 2 t e μ s I ( s ) d s σ I ρ S d t + σ 1 S ( t ) d B 1 ( t ) + σ 2 E ( t ) d B 2 ( t ) + σ 3 I ( t ) d B 3 ( t ) + Z [ γ 1 ( u ) S ( t ) + γ 2 ( u ) E ( t ) + γ 3 ( u ) I ( t ) ] N ˜ ( d t , d u ) Λ μ S + E + I + ρ e μ t t τ 1 t e μ s S ( s ) d s + γ e μ t t τ 2 t e μ s I ( s ) d s d t + σ 1 S ( t ) d B 1 ( t ) + σ 2 E ( t ) d B 2 ( t ) + σ 3 I ( t ) d B 3 ( t ) + Z [ γ 1 ( u ) S ( t ) + γ 2 ( u ) E ( t ) + γ 3 ( u ) I ( t ) ] N ˜ ( d t , d u ) .
Then, taking integration of the inequality from 0 to t, dividing by t, and taking a superior limit on both sides of (10). Applying the corresponding Lemma 4, we obtain that
lim sup t S ( t ) + E ( t ) + I ( t ) + ρ e μ t t τ 1 t e μ s S ( s ) d s ) + γ e μ t t τ 2 t e μ s I ( s ) d s Λ μ a . s .
Therefore
lim sup t S ( t ) Λ μ , lim sup t E ( t ) Λ μ a . s .
Substituting (11) into (9) gives
lim sup t ln I ( t ) t β · Λ μ + α β · Λ μ ( γ + σ + μ + θ 3 ) = ( γ + σ + μ + θ 3 ) ( R 1 1 ) .
If the condition R 1 < 1 holds, then
lim t I ( t ) = 0 .
Furthermore
d S ( t ) + I ( t ) + ρ e μ τ 1 t τ 1 t S ( s ) d s + γ e μ τ 2 t τ 2 t I ( s ) d s = [ Λ + α β E ( t ) I ( t ) 1 + p 2 I ( t ) ϵ S ( t ) μ S ( t ) ( γ + σ + μ ) I ( t ) ρ S ( t ) + ρ S ( t τ 1 ) e μ τ 1 + γ I ( t τ 2 ) e μ τ 2 + ρ e μ τ 1 S ( t ) ρ e μ τ 1 S ( t τ ) + γ e μ τ 2 I ( t ) γ e μ τ 2 I ( t τ ) ] d t + σ 1 S ( t ) d B 1 ( t ) + σ 3 I ( t ) d B 3 ( t ) + Z γ 1 ( u ) S ( t ) N ˜ ( d t , d u ) + Z γ 3 ( u ) I ( t ) N ˜ ( d t , d u ) = Λ + α β E ( t ) I ( t ) 1 + p 2 I ( t ) ϵ S ( t ) μ S ( t ) ( γ + σ + μ ) I ( t ) ρ S ( t ) + ρ e μ τ 1 S ( t ) + γ e μ τ 2 I ( t ) d t + σ 1 S ( t ) d B 1 ( t ) + σ 3 I ( t ) d B 3 ( t ) + Z γ 1 ( u ) S ( t ) N ˜ ( d t , d u ) + Z γ 3 ( u ) I ( t ) N ˜ ( d t , d u ) ,
then
S ( t ) S ( 0 ) t + I ( t ) I ( 0 ) t + ρ e μ τ 1 t τ 1 t S ( s ) d s ρ e μ τ 1 τ 1 0 S ( s ) d s t + γ e μ τ 2 t τ 2 t I ( s ) d s γ e μ τ 2 τ 2 0 I ( s ) d s t = Λ + α β E ( t ) I ( t ) 1 + p 2 I ( t ) ( ϵ + μ + ρ ρ e μ τ 1 ) S ( t ) ( γ + σ + μ γ e μ τ 2 ) I ( t ) + σ 1 t 0 t S ( s ) d B 1 ( s ) + σ 3 t 0 t I ( s ) d B 3 ( s ) + 1 t 0 t Z γ 1 ( u ) S ( s ) N ˜ ( d s , d u ) + 1 t 0 t Z γ 3 ( u ) I ( s ) N ˜ ( d s , d u ) ,
thus
α β E ( t ) I ( t ) 1 + p 2 I ( t ) = ( ϵ + μ + ρ ρ e μ τ 1 ) S ( t ) Λ + ( γ + σ + μ γ e μ τ 2 ) I ( t ) + φ 1 ( t ) ,
where
φ 1 ( t ) = S ( t ) S ( 0 ) t + I ( t ) I ( 0 ) t + ρ e μ τ 1 t τ 1 t S ( s ) d s ρ e μ τ 1 τ 1 0 S ( s ) d s t + γ e μ τ 2 t τ 2 t I ( s ) d s γ e μ τ 2 τ 2 0 I ( s ) d s t σ 1 t 0 t S ( s ) d B 1 ( s ) σ 3 t 0 t I ( s ) d B 3 ( s ) 1 t 0 t Z γ 1 ( u ) S ( s ) N ˜ ( d s , d u ) 1 t 0 t Z γ 3 ( u ) I ( s ) N ˜ ( d s , d u ) .
In view of Lemmas 4 and 5, one can see that
lim t φ 1 ( t ) = 0 .
Similarly, according to the second equation of (4), we have
d E ( t ) = ϵ S ( t ) α β E ( t ) I ( t ) 1 + p 2 I ( t ) μ E ( t ) d t + σ 2 E ( t ) d B 2 ( t ) + Z γ 2 ( u ) E ( t ) N ˜ ( d t , d u ) ,
we can obtain
E ( t ) E ( 0 ) t = ϵ S ( t ) α β E ( t ) I ( t ) 1 + p 2 I ( t ) μ E ( t ) + σ 2 t 0 t E ( s ) d B 2 ( s ) + 1 t 0 t Z γ 2 ( u ) E ( s ) N ˜ ( d s , d u ) ,
so we have
α β E ( t ) I ( t ) 1 + p 2 I ( t ) = ϵ S ( t ) μ E ( t ) + φ 2 ( t ) ,
here
φ 2 ( t ) = E ( t ) E ( 0 ) t σ 2 t 0 t E ( s ) d B 2 ( s ) 1 t 0 t Z γ 2 ( u ) E ( s ) N ˜ ( d s , d u ) .
Accordingly, Lemma 4, leads to
lim t φ 2 ( t ) = 0 .
Combining (13) and(15) gives
( μ + ρ ρ e μ τ 1 ) S ( t ) + μ E ( t ) = Λ ( γ + σ + μ γ e μ τ 2 ) I ( t ) φ 1 ( t ) + φ 2 ( t ) ,
By (12), (14), and (16), one can easily obtain that
lim t [ ( μ + ρ ρ e μ τ 1 ) S ( t ) + μ E ( t ) ] = Λ .
In addition, owing to
d S ( t ) + ρ e μ τ 1 t τ 1 t S ( s ) d s + γ e μ τ 2 t τ 2 t I ( s ) d s = [ Λ β S ( t ) I ( t ) 1 + p 1 I ( t ) + ρ S ( t τ 1 ) e μ τ 1 + γ I ( t τ 2 ) e μ τ 2 ( ϵ + μ ) S ( t ) ρ S ( t ) + ρ S ( t ) e μ τ 1 ρ S ( t τ 1 ) e μ τ 1 + γ I ( t ) e μ τ 2 γ I ( t τ 2 ) e μ τ 2 ] d t + σ 1 S ( t ) d B 1 ( t ) + Z γ 1 ( u ) S ( t ) N ˜ ( d t , d u ) ,
then
d S ( t ) + ρ e μ τ 1 t τ 1 t S ( s ) d s + γ e μ τ 2 t τ 2 t I ( s ) d s [ Λ + ρ S ( t ) e μ τ 1 + γ I ( t ) e μ τ 2 ( ϵ + μ + ρ ) S ( t ) ] d t + σ 1 S ( t ) d B 1 ( t ) + Z γ 1 ( u ) S ( t ) N ˜ ( d t , d u ) .
Hence,
S ( t ) Λ ϵ + μ + ρ ( 1 e μ τ 1 ) + γ e μ τ 2 ϵ + μ + ρ ( 1 e μ τ 1 ) I ( t ) + φ 3 ( t ) ,
where
φ 3 ( t ) = 1 ϵ + μ + ρ ( 1 e μ τ 1 ) [ S ( t ) S ( 0 ) t + ρ e μ τ 1 t τ 1 t S ( s ) d s ρ e μ τ 1 τ 1 0 S ( s ) d s t + γ e μ τ 2 t τ 2 t I ( s ) d s γ e μ τ 2 τ 2 0 I ( s ) d s t σ 1 t 0 t S ( s ) d B 1 ( s ) 1 t 0 t Z γ 1 ( u ) S ( s ) N ˜ ( d s , d u ) ] .
Furthermore, we have
lim t φ 3 ( t ) = 0 .
Therefore
lim t S ( t ) Λ ϵ + μ + ρ ( 1 e μ τ 1 ) a . s .
Likewise
lim t E ( t ) ϵ Λ μ ( ϵ + μ + ρ ( 1 e μ τ 1 ) ) a . s .
Thus, from (18), we can achieve that
lim t S ( t ) = Λ ϵ + μ + ρ ( 1 e μ τ 1 ) a . s . ,
and
lim t E ( t ) = ϵ Λ μ ( ϵ + μ + ρ ( 1 e μ τ 1 ) ) a . s .
The proof is completed. □
Remark 1.
The expression from R 1 can be divided into two parts, where the second part α β Λ μ ( γ + σ + μ + θ 3 ) indicates the effect of public health education on the basic reproduction number.

5. Persistence in the Mean of the Diseases

Now, we will discuss the sufficient condition of persistence in the mean of the disease after presenting a related definition as follows.
Definition 2.
If lim inf t I ( t ) > 0 a.s., the infective individuals I ( t ) of system (4) have persistence in the mean.
Define
R 2 = 2 β Λ ( ϵ + μ + ρ + θ 1 ) 1 2 ( γ + σ + μ + θ 3 + 2 ) + 2 α β ϵ Λ ( ϵ + μ + ρ + θ 1 ) 1 2 ( μ + θ 2 ) 1 2 ( γ + σ + μ + θ 3 + 2 ) ,
where θ i = 1 2 σ i 2 + Z [ γ i ( u ) ln ( 1 + γ i ( u ) ) ] λ ( d u ) , i = 1 , 2 , 3 .
Theorem 3.
Suppose that Assumption 1 holds, and let ( S ( t ) , E ( t ) , I ( t ) ) be a solution of system (4) with initial condition (5), then, if R 2 > 1 , we can obtain the disease I ( t ) , which is persistence in the mean. Moreover, it satisfies
lim inf t I ( t ) γ + σ + μ + θ 3 + 2 ( a 1 + b 1 ) β + b 2 α β + p 1 + p 2 ( R 2 1 ) > 0 a . s . ,
where
a 1 = β Λ ( ε + μ + ρ + θ 1 ) 3 2 , b 1 = α β ε Λ ( ε + μ + ρ + θ 1 ) 3 2 ( μ + θ 2 ) 1 2 , b 2 = α β ε Λ ( ε + μ + ρ + θ 1 ) 1 2 ( μ + θ 2 ) 3 2 .
Proof. 
Define the C 2 -function H ¯ : R + 3 R + as follows:
H ¯ ( S ( t ) , E ( t ) , I ( t ) ) = ( a 1 + b 1 ) ln S b 2 ln E ln I ,
which a 1 , b 1 , b 2 are positive constants to be determined. Obviously, the function H ¯ ( S , E , I ) has a minimum point ( S ¯ , E ¯ , I ¯ ) inside R + 3 . Thus, define a non-negative C 2 -function H: H ( S , E , I ) = H ¯ ( S , E , I ) H ¯ ( S ¯ , E ¯ , I ¯ ) . Note that using the Itô’s formula for the function yields
d H = L H d t ( a 1 + b 1 ) σ 1 d B 1 ( t ) b 2 σ 2 d B 2 ( t ) σ 3 d B 3 ( t ) Z [ ( a 1 + b 1 ) ln ( 1 + γ 1 ( u ) ) + b 2 ln ( 1 + γ 2 ( u ) ) + ln ( 1 + γ 3 ( u ) ) ] N ˜ ( d t , d u ) ,
which gives that
L H = ( a 1 + b 1 ) Λ S + β I 1 + p 1 I ρ S ( t τ 1 ) e μ τ 1 S γ I ( t τ 2 ) e μ τ 2 S + ( ϵ + μ ) + ρ + σ 1 2 2 + β S 1 + p 1 I α β E 1 + p 2 I + ( γ + σ + μ ) + σ 3 2 2 + b 2 ϵ S E + α β I 1 + p 2 I + μ + σ 2 2 2 + Z ( a 1 + b 1 ) [ γ 1 ln ( 1 + γ 1 ) ] + b 2 [ γ 2 ln ( 1 + γ 2 ) ] + [ γ 3 ln ( 1 + γ 3 ) ] λ ( d u ) β S 1 + p 1 I a 1 Λ S ( 1 + p 1 I ) + α β E 1 + p 2 I b 1 Λ S b 2 ϵ S E ( 1 + p 2 I ) + a 1 ( ϵ + μ + ρ + σ 1 2 2 ) + b 1 ϵ + μ + ρ + σ 1 2 2 + b 2 μ + σ 2 2 2 + γ + σ + μ + σ 3 2 2 + ( a 1 + b 1 ) β I + b 2 α β I + ( 1 + p 1 I ) + ( 1 + p 2 I ) + Z { ( a 1 + b 1 ) [ γ 1 ln ( 1 + γ 1 ) ] + b 2 [ γ 2 ln ( 1 + γ 2 ) ] + [ γ 3 ln ( 1 + γ 3 ) ] } λ ( d u ) 3 a 1 β Λ 3 4 b 1 b 2 α β ϵ Λ 4 + a 1 ( ϵ + μ + ρ + θ 1 ) + b 1 ( ϵ + μ + ρ + θ 1 ) + b 2 ( μ + θ 2 ) + ( γ + σ + μ + θ 3 + 2 ) + [ ( a 1 + b 1 ) β + b 2 α β + p 1 + p 2 ] I ,
where θ i = 1 2 σ i 2 + Z [ γ i ( u ) ln ( 1 + γ i ( u ) ) ] λ ( d u ) , i = 1 , 2 , 3 . Let
a 1 ( ϵ + μ + ρ + θ 1 ) = β Λ ( ϵ + μ + ρ + θ 1 ) 1 2 ,
b 1 ( ϵ + μ + ρ + θ 1 ) = b 2 ( μ + θ 2 ) = α β ϵ Λ ( ϵ + μ + ρ + θ 1 ) 1 2 ( μ + θ 2 ) 1 2 .
It is easy to check that
a 1 = β Λ ( ϵ + μ + ρ + θ 1 ) 3 2 , b 1 = α β ϵ Λ ( ϵ + μ + ρ + θ 1 ) 3 2 ( μ + θ 2 ) 1 2 , b 2 = α β ϵ Λ ( ϵ + μ + ρ + θ 1 ) 1 2 ( μ + θ 2 ) 3 2 .
Substituting (24) into (23), and then deriving that
L H 2 β Λ ( ϵ + μ + ρ + θ 1 ) 1 2 2 α β ϵ Λ ( ϵ + μ + ρ + θ 1 ) 1 2 ( μ + θ 2 ) 1 2 + ( γ + σ + μ + θ 3 + 2 ) + [ ( a 1 + b 1 ) β + b 2 α β + p 1 + p 2 ] I .
We integrate the inequality (22) from 0 to t and divided by t, then substituting (25) into it, thus we can derive that
H ( S ( t ) , E ( t ) , I ( t ) ) H ( S ( 0 ) , E ( 0 ) , I ( 0 ) ) t ( γ + σ + μ + θ 3 + 2 ) ( R 2 1 ) + [ ( a 1 + b 1 ) β + b 2 α β + p 1 + p 2 ] I ( t ) ϕ ( t ) t ,
where
ϕ ( t ) = b 2 0 t σ 2 d B 2 ( s ) + ( a 1 + b 1 ) 0 t σ 1 d B 1 ( s ) + 0 t σ 3 d B 3 ( s ) + 0 t Z [ b 2 ln ( 1 + γ 2 ( u ) ) + ( a 1 + b 1 ) ln ( 1 + γ 1 ( u ) ) + ln ( 1 + γ 3 ( u ) ) ] N ˜ ( d s , d u ) .
The inequality (26) can be rewritten as
I ( t ) 1 ( a 1 + b 1 ) β + b 2 α β + p 1 + p 2 [ ( γ + σ + μ + θ 3 + 2 ) ( R 2 1 ) + ϕ ( t ) t + H ( S ( t ) , E ( t ) , I ( t ) ) H ( S ( 0 ) , E ( 0 ) , I ( 0 ) ) t ] .
Similarly, we have
lim t ϕ ( t ) t = 0 a . s .
Taking the inferior limit of both sides of (27), and by R 2 > 1 yields
lim inf t I ( t ) γ + σ + μ + θ 3 + 2 ( a 1 + b 1 ) β + b 2 α β + p 1 + p 2 ( R 2 1 ) > 0 a . s .
This completes the proof. □

6. Numerical Simulations and Discussion

Model (1) typically applies to diseases found in remote or impoverished areas characterized by limited access to education, insufficient information, poor public health conditions, and extreme weather events. Such diseases include influenza, tungiasis, leptospirosis, and novel coronavirus pneumonia. Numerical simulations in this study use tungiasis as an illustrative example. Tungiasis is a parasitic skin disease caused by the invasion of female sand fleas into the human epidermis, commonly found in tropical, economically disadvantaged regions such as Latin America, the Caribbean, and sub-Saharan Africa. Detailed information about tungiasis can be found in reference [31]. Now, the numerical simulations are performed using the Milstein method [32] and Euler numerical approximation [33]. Assume that the total number of the population is 800, and the initial value is ( S ( 0 ) , E ( 0 ) , I ( 0 ) , R ( 0 ) ) = ( 300 , 200 , 125 , 175 ) . Parameter values are derived from references [5,34,35,36,37] and are presented in Table 1.
For the assumed parameters, we set α = 0.22 , ρ = 0.6 , Z = ( 0 , + ) , λ ( Z ) = 1 , the intensity of stochastic disturbance as follows:
σ 1 2 = σ 2 2 = σ 4 2 = 0.01 , σ 3 2 = 0.02 , γ 1 = 0.01 , γ 2 = 0.02 , γ 3 = 0.08 , γ 4 = 0.04 .
It is easy to calculate that R 1 0.9932 < 1 , and the condition of Theorem 2 is satisfied. The numerical simulation results are given in Figure 1. It can be found that the infected population I ( t ) is extinct in both the deterministic model and the stochastic model. In fact, all values of β which are less than 0.01498 lead to obtaining a value less than one of the basic reproduction numbers R 1 of the system (1), so the disease will always be extinct.
The coefficient β can likewise increase if we increase the contact between the infected population and the susceptible population. Similarly, if we reduce the effective vaccination rate among the infected population, the coefficient ρ will decrease. For this analysis, we retain the parameter values used in Table 1, except for β and ρ , which are chosen to be 0.114 and 0.01, respectively. The values of stochastic noise intensity are as follows:
σ 1 2 = 0.003 , σ 2 2 = σ 3 2 = 0.002 , σ 4 2 = 0.005 , γ 1 = γ 3 = 0.006 , γ 2 = 0.003 , γ 4 = 0.008 .
In this case, R 2 1.2560 > 1 , satisfying the conditions of Theorem 3, and the numerical simulations are shown in Figure 2. We can find that the trend of persistence in the mean of the infected population I ( t ) is consistent with the deterministic model after adding stochastic perturbations. Furthermore, we can calculate lim inf t I ( t ) 0.0973 , which shows the minimum level of infected population in the average sense. By comparing Figure 1 and Figure 2, it becomes evident that the disease contact rate β and vaccination rate ρ significantly influence the extinction and persistence of the disease.
Next, we also explore the effect of stochastic disturbance intensity on disease extinction. Parameters γ and ρ are set at 0.02 and 0.03, respectively, with other parameter values in Table 1 remain unchanged.
The numerical simulations are presented in Figure 3. When we select the intensities of disturbance are small, such as σ 1 2 = σ 2 2 = σ 4 2 = 0.01 , σ 3 2 = 0.02 , γ 1 = 0.01 , γ 2 = 0.02 , γ 3 = 0.08 , γ 4 = 0.04 , it is easy to calculate that R 1 5.0746 > 1 , the infected population I ( t ) is not extinct in both the deterministic model and stochastic models. When we select larger disturbance intensities, σ 1 2 = σ 2 2 = σ 4 2 = 0.1 , σ 3 2 = 0.7 , γ 1 = 0.1 , γ 2 = 0.2 , γ 3 = 0.6 , γ 4 = 0.4 , it is easy to calculate that R 1 0.8879 < 1 , so the condition of Theorem 2 is satisfied. The infected population I ( t ) is extinct in the stochastic models with and without jumps (see Figure 3d,e), which indicates that the large intensities of environmental disturbance have a positive effect on the disease extinction. In addition, by comparing Figure 3b,d and Figure 3c,e, we find that the results of stochastic disturbance with jump are consistent with those of stochastic disturbance without jump. Discontinuous environmental shocks described by L e ´ vy jumps can also suppress the spread of the disease in our model.
To further analyze the specific effects of different parameters on disease extinction and persistence, parameter sensitivity analysis is as follows:
For disease extinction and persistence, the effects of parameters α and β on I ( t ) are considered. At the same time, according to the sensitivity analysis formula in reference [38], we obtain
R 1 α · α R 1 = α 1 + α > 0 , R 2 α · α R 2 = ( α ϵ ) 1 2 2 [ ( μ + θ 2 ) 1 2 + ( α ϵ ) 1 2 ] > 0 ,
R 1 β · β R 1 = Λ ( 1 + α ) · β β Λ ( 1 + α ) = 1 > 0 , R 2 β · β R 2 = ( μ + θ 2 ) 1 2 + ( α ϵ ) 1 2 2 [ ( μ + θ 2 ) 1 2 + ( α ϵ ) 1 2 ] > 0 .
The results of numerical simulations are shown in Figure 4: (a) select α = 0.22 , 0.11 , 0.01 , respectively, so we can see that the smaller α is, the faster population I ( t ) will go extinct; (b) Select α = 0.22 , 0.42 , 0.62 , respectively, so we can see that I ( t ) decreases as α decreases; (c) Select β = 0.01498 , 0.014 , 0.01 , respectively, we can see that the smaller β is, the faster population I ( t ) will go extinct; and (d) select β = 0.114 , 0.124 , 0.134 , respectively, so we can see that I ( t ) decreases as β decreases.
The decrease in parameter α and β will lead to the decrease in R 1 , R 2 , and it can be concluded that parameters α and β have a positive correlation with the basic reproduction number. That is to say, the closer α is to 0, the better the effect of education on disease infection prevention, as well as the smaller β is, and the less contact susceptible people have with infected people, the fewer infected people there are. This indicates that public health education and disease contact rate have a certain role in promoting disease extinction.
For disease extinction and persistence, the effects of parameters γ and ρ on I ( t ) are also considered. According to the sensitivity analysis formulas, we obtain
R 1 γ · γ R 1 = γ γ + σ + μ + θ 3 < 0 , R 2 γ · γ R 2 = γ γ + σ + μ + θ 3 + 2 < 0 ,
R 2 ρ · ρ R 2 = ρ 2 ( ϵ + μ + ρ + θ 1 ) < 0 .
The results of numerical simulations are shown in Figure 5: (a) select γ = 0.427 , 0.527 , 0.627 , respectively, so we can see that the smaller γ is, the more time the population I ( t ) will take to go extinct; (b) select γ = 0.427 , 0.327 , 0.227 , respectively, we can see that I ( t ) increases as γ decreases; (c) select ρ = 0.6 , 0.7 , 0.8 , respectively, we can see that the smaller ρ is, the more time the population I ( t ) will take to go extinct; (d) Select ρ = 0.01 , 0.02 , 0.03 , respectively, we can see that I ( t ) increases as γ decreases.
The decrease in parameters γ and ρ will lead to the increase in R 1 , R 2 , and it can be concluded that parameter γ and ρ have a negative correlation with the basic reproduction number. In other words, the smaller γ is, the fewer the number of infected people will recover, as well as the smaller ρ is, the fewer the number of susceptible people will be vaccinated, such that a greater number of infected people indicates that the disease recovery rate and the vaccination rate have a certain promoting effect on disease extinction.

7. Conclusions

In this paper, we propose a stochastic delayed SEIRS epidemic model considering the impacts of public health education and temporary immunity. We prove the global existence and uniqueness of the positive solution of the system (4), providing sufficient conditions for the extinction and persistence of the disease. We perform numerical simulations to corroborate our results and analyze each parameter’s influence on the disease. In conclusion, factors such as public health education, disease contact rate, disease recovery rate, vaccination rate, and high stochastic noises can all contribute to disease extinction. The higher the levels of education, recovery and vaccination rates, and the level of noise, coupled with a lower rate of contact, the more likely it is for the disease to become extinct. Notably, the basic reproduction numbers R 1 and R 2 are independent of time delay, suggesting that time delay might not significantly affect the disease’s average persistence and extinction.
The study in reference [5] considered the deterministic susceptible–educated–infective–susceptible (SEIS) model impacted by public health education. The disease extinction threshold R 0 = β Λ / ( μ ( γ + σ + μ ) ) obtained in reference [5] shows a strong relationship between disease extinction and the contact rate, recovery rate, and mortality rate. This paper expands on that by also considering the intensity of stochastic disturbance, thus demonstrating that the average persistence and extinction of infectious diseases are influenced by a combination of these factors. Furthermore, large-intensity noise disturbances can lead to disease extinction, as illustrated in Figure 3. Similarly, both continuous interference described by white Gaussian noise and discontinuous environmental shocks represented by L e ´ vy jumps can suppress the spread of the disease in our model.
This work provides a fundamental stochastic delayed dynamical model, which can be instrumental in understanding the impact of public health education and stochastic noises on infectious disease prevention and control. In terms of real-world disease prevention and control efforts, the findings suggest potential measures, such as increasing education awareness, heightening environmental disturbance, minimizing patient contact with others, and encouraging advance vaccination.

Author Contributions

Conceptualization, X.S.; methodology, X.S.; software, D.Z.; validation, D.Z.; formal analysis, D.Z.; writing—original draft preparation, D.Z.; writing—review and editing, X.S. and X.Z.; supervision, X.Z.; funding acquisition, X.S. and X.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This work is sponsored by the Natural Science Foundation of Henan (222300420521) and the Nanhu Scholars Program for Young Scholars of XYNU.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lu, H. Prevention and Control of Infectious Diseases. Highlights Sci. Eng. Technol. 2023, 36, 871–874. [Google Scholar] [CrossRef]
  2. Ghaffar, A.; Rashid, S.F.; Wanyenze, R.K.; Hyder, A.A. Public health education post-COVID-19: A proposal for critical revisions. BMJ Glob. Health 2021, 6, e005669. [Google Scholar] [CrossRef]
  3. Eric, P.; Warner, S.B. A potential new front in health communication to encourage vaccination: Health education teachers. Vaccine 2021, 39, 4671–4677. [Google Scholar]
  4. Cheng, G.; Wang, N.; Qi, L. Modeling the influence of health education on individual behavior in the process of disease transmission. Math. Methods Appl. Sci. 2022, 45, 498–525. [Google Scholar] [CrossRef]
  5. Nyang’inja, R.A.; Angwenyi, D.N.; Musyoka, C.M.; Orwa, T.O. Mathematical modeling of the effects of public health education on tungiasis—A neglected disease with many challenges in endemic communities. Adv. Differ. Equ. 2018, 2018, 426. [Google Scholar] [CrossRef]
  6. Buonomo, B.; Carbone, G.; d’Onofrio, A. Effect of seasonality on the dynamics of an imitation-based vaccination model with public health intervention. Math. Biosci. Eng. 2018, 15, 299–321. [Google Scholar]
  7. Shi, X.; Gao, X.; Zhou, X.; Li, Y. Analysis of an SQEIAR epidemic model with media coverage and asymptomatic infection. AIMS Math. 2021, 6, 12298–12320. [Google Scholar] [CrossRef]
  8. Joshi, H.; Sharma, R.K.; Prajapati, G.L. Stability analysis of a deterministic vaccination model with non-monotonic incidence rate. J. Math. Comput. Sci. 2020, 10, 51–67. [Google Scholar]
  9. Fan, K.; Zhang, Y.; Gao, S.; Wei, X. A class of stochastic delayed SIR epidemic models with generalized nonlinear incidence rate and temporary immunity. Phys. A Stat. Mech. Appl. 2017, 481, 198–208. [Google Scholar] [CrossRef]
  10. Liu, Y.; Zhang, Y.; Wang, Q. A stochastic SIR epidemic model with Lévy jump and media coverage. Adv. Differ. Equ. 2020, 2020, 70. [Google Scholar] [CrossRef] [Green Version]
  11. Kyrychko, Y.N.; Blyuss, K.B. Global properties of a delayed SIR model with temporary immunity and nonlinear incidence rate. Nonlinear Anal. Real World Appl. 2005, 6, 495–507. [Google Scholar] [CrossRef]
  12. Xu, R.; Ma, Z.; Wang, Z. Global stability of a delayed SIRS epidemic model with saturation incidence and temporary immunity. Comput. Math. Appl. 2010, 59, 3211–3221. [Google Scholar] [CrossRef] [Green Version]
  13. Li, R.; Tian, F.; Hu, X. The stochastic SIV epidemic model with varying coefficient and time delay. Math. Appl. 2016, 29, 890–896. [Google Scholar]
  14. Liu, Q.; Jiang, D.; Hayat, T.; Ahmad, B. Analysis of a delayed vaccinated SIR epidemic model with temporary immunity and Lévy jumps. Nonlinear Anal. Hybrid Syst. 2018, 27, 29–43. [Google Scholar] [CrossRef]
  15. Hu, Z.; Cheng, X.; Ma, W. Analysis of an SIS epidemic model with temporary immunity and nonlinear incidence rate. Chin. J. Eng. Math. 2009, 26, 407–415. [Google Scholar]
  16. Liu, Q.; Chen, Q. Stochastic delayed SIR epidemic with vaccination model and temporary immunity. Sci. Sin. Math. 2016, 46, 1745–1756. [Google Scholar]
  17. Shi, X.; Cao, Y.; Zhou, X. Dynamics for a stochastic delayed SIRS epidemic model. Nonlinear Anal. Model. Control. 2020, 25, 705–725. [Google Scholar] [CrossRef]
  18. Zhou, X.; Shi, X.; Wei, M. Dynamical behavior and optimal control of a stochastic mathematical model for cholera. Chaos Solitons Fractals 2022, 156, 111854. [Google Scholar] [CrossRef]
  19. Shi, X.; Cao, Y. Dynamics of a stochastic periodic SIRS model with time delay. Int. J. Biomath. 2020, 13, 2050072. [Google Scholar] [CrossRef]
  20. Zhang, X.; Shao, Y. Analysis of a stochastic predator-prey system with mixed functional responses and Lévy jumps. AIMS Math. 2021, 6, 4404–4427. [Google Scholar] [CrossRef]
  21. Fan, C.; Xiong, Z.; Nie, D. Dynamical analysis of a stochastic ratio-dependent predator-prey system with Lévy jumps. Int. J. Appl. Math. Stat. 2014, 52, 169–177. [Google Scholar]
  22. Kiouach, D.; El-idrissi, S.E.A.; Sabbar, Y. A novel mathematical analysis and threshold reinforcement of a stochastic dengue epidemic model with Lévy jumps. Commun. Nonlinear Sci. Numer. Simul. 2023, 119, 107092. [Google Scholar] [CrossRef]
  23. Liu, Q.; Chen, Q. Analysis of a general stochastic non-autonomous logistic model with delays and Lévy jumps. J. Math. Anal. Appl. 2016, 433, 95–120. [Google Scholar] [CrossRef]
  24. Lu, C.; Ding, X. Dynamical behavior of stochastic delay Lotka–Volterra competitive model with general Lévy jumps. Phys. A Stat. Mech. Its Appl. 2019, 531, 121730. [Google Scholar] [CrossRef]
  25. Liu, R.; Liu, G. Analysis of a stochastic predator-prey population model with Allee effect and jumps. J. Inequalities Appl. 2019, 2019, 60. [Google Scholar] [CrossRef] [Green Version]
  26. Wang, K. Stochastic Biomathematical Models; Science Press: Beijing, China, 2010. (In China) [Google Scholar]
  27. David, A. Lévy Processes and Stochastic Calculus; Cambridge University Press: Cambridge, UK, 2004. [Google Scholar]
  28. Zhao, Y.; Yuan, S. Stability in distribution of a stochastic hybrid competitive Lotka–Volterra model with Lévy jumps. Chaos Solitons Fractals 2016, 85, 98–109. [Google Scholar] [CrossRef]
  29. Liptser, R.S. A strong law of large numbers for local martingales. Stochastics 1980, 3, 217–228. [Google Scholar] [CrossRef]
  30. Liu, Q.; Chen, Q.; Jiang, D. The threshold of a stochastic delayed SIR epidemic model with temporary immunity. Phys. A Stat. Mech. Its Appl. 2016, 450, 115–125. [Google Scholar] [CrossRef]
  31. Ehrenberg, J.P.; Ault, S.K. Neglected diseases of neglected populations: Thinking to reshape the determinants of health in Latin America and the Caribbean. BMC Public Health 2005, 5, 119. [Google Scholar] [CrossRef] [Green Version]
  32. Higham, D.J. An algorithmic introduction to numerical simulation of stochastic differential equations. SIAM Rev. 2001, 43, 525–546. [Google Scholar] [CrossRef]
  33. Protter, P.; Talay, D. The Euler scheme for Lévy driven stochastic differential equations. Ann. Probab. 1997, 25, 393–423. [Google Scholar] [CrossRef]
  34. Mussa, M.; Vivian, M.; Valeria, S.; Hoseenu, P.; Winfrida, J.; Eliud, Y.Y.; Donath, T. Tungiasis infection among primary school children in Northeastern Tanzania: Prevalence, intensity, clinical aspects and associated factors. IJID Reg. 2023, 7, 116–123. [Google Scholar]
  35. CIA. The CIA World Factbook; Skyhorse Publishing: New York, NY, USA, 2011. [Google Scholar]
  36. IndexMundi. Kenya Demographic Profile. 2018. Available online: https://www.indexmundi.com/kenya/demographics_profile.html (accessed on 18 September 2021).
  37. Kong, L.; Li, L.; Kang, S.; Liu, Y.; Feng, W. Dynamic Behavior of a Stochastic Tungiasis Model for Public Health Education. Discret. Dyn. Nat. Soc. 2022, 2022, 4927261. [Google Scholar] [CrossRef]
  38. Samsuzzoha, M.; Singh, M.; Lucy, D. Uncertainty and sensitivity analysis of the basic reproduction number of a vaccinated epidemic model of influenza. Appl. Math. Model. 2013, 37, 903–915. [Google Scholar] [CrossRef]
Figure 1. Numerical simulations of the SEIRS deterministic model and stochastic model of disease extinction, where (a) is the deterministic model and (b) is the stochastic model without jumps and (c) is the stochastic model with jumps.
Figure 1. Numerical simulations of the SEIRS deterministic model and stochastic model of disease extinction, where (a) is the deterministic model and (b) is the stochastic model without jumps and (c) is the stochastic model with jumps.
Axioms 12 00560 g001
Figure 2. Numerical simulations of the SEIRS deterministic model and stochastic model of disease are persistent, where (a) is the deterministic model, and (b) is the stochastic model without jumps and (c) is the stochastic model with jumps.
Figure 2. Numerical simulations of the SEIRS deterministic model and stochastic model of disease are persistent, where (a) is the deterministic model, and (b) is the stochastic model without jumps and (c) is the stochastic model with jumps.
Axioms 12 00560 g002
Figure 3. Numerical simulations of the SEIRS deterministic model and stochastic model of disease extinction, where (a) is the deterministic models, (b,d) are the stochastic models without jumps and (c,e) are the stochastic models with jumps. Here, the disturbance intensity parameters in (b) are σ 1 2 = σ 2 2 = σ 4 2 = 0.01 , σ 3 2 = 0.02 , in (c) are σ 1 2 = σ 2 2 = σ 4 2 = 0.01 , σ 3 2 = 0.02 , γ 1 = 0.01 , γ 2 = 0.02 , γ 3 = 0.08 , γ 4 = 0.04 , in (d) is σ 1 2 = σ 2 2 = σ 4 2 = 0.1 , σ 3 2 = 0.7 , and in (e) are σ 1 2 = σ 2 2 = σ 4 2 = 0.1 , σ 3 2 = 0.7 , γ 1 = 0.1 , γ 2 = 0.2 , γ 3 = 0.6 , γ 4 = 0.4 .
Figure 3. Numerical simulations of the SEIRS deterministic model and stochastic model of disease extinction, where (a) is the deterministic models, (b,d) are the stochastic models without jumps and (c,e) are the stochastic models with jumps. Here, the disturbance intensity parameters in (b) are σ 1 2 = σ 2 2 = σ 4 2 = 0.01 , σ 3 2 = 0.02 , in (c) are σ 1 2 = σ 2 2 = σ 4 2 = 0.01 , σ 3 2 = 0.02 , γ 1 = 0.01 , γ 2 = 0.02 , γ 3 = 0.08 , γ 4 = 0.04 , in (d) is σ 1 2 = σ 2 2 = σ 4 2 = 0.1 , σ 3 2 = 0.7 , and in (e) are σ 1 2 = σ 2 2 = σ 4 2 = 0.1 , σ 3 2 = 0.7 , γ 1 = 0.1 , γ 2 = 0.2 , γ 3 = 0.6 , γ 4 = 0.4 .
Axioms 12 00560 g003
Figure 4. Numerical simulations of the population I ( t ) of disease extinction and persistence under different parameter α and β , where (a) is disease extinction of α , and (b) is disease persistence of α ; (c) is the disease extinction of β ; and (d) is the disease persistence of β .
Figure 4. Numerical simulations of the population I ( t ) of disease extinction and persistence under different parameter α and β , where (a) is disease extinction of α , and (b) is disease persistence of α ; (c) is the disease extinction of β ; and (d) is the disease persistence of β .
Axioms 12 00560 g004
Figure 5. Numerical simulations of the population I ( t ) of disease extinction and persistence under different parameters γ and ρ , where (a) is disease extinction of γ , and (b) is disease persistence of γ , (c) is disease extinction of ρ , and (d) is the disease persistence of ρ .
Figure 5. Numerical simulations of the population I ( t ) of disease extinction and persistence under different parameters γ and ρ , where (a) is disease extinction of γ , and (b) is disease persistence of γ , (c) is disease extinction of ρ , and (d) is the disease persistence of ρ .
Axioms 12 00560 g005
Table 1. Parameter values of the system (1) used in Figure 1 and Figure 2.
Table 1. Parameter values of the system (1) used in Figure 1 and Figure 2.
SymbolParameterValueSource
Λ Recruitment rate of the total population0.44 day 1 [35]
β Effective contact rate0.01498 day 1 [36]
ϵ Dissemination of education rate0.01431 day 1 [5]
α Scaling factor0.22Assumed 1
p 1 , 2 Coefficient of saturated incidence rate0.8 day 1 Assumed 1
ρ Efficient vaccination rate for the susceptible class[0.01–0.6] day 1 Assumed 1
μ Per capita natural death rate0.016 day 1 [35]
γ Recovery rate0.427 day 1 [36]
σ Disease-induced mortality rate0.05 day 1 [36]
τ 1 , 2 Length of temporary immunity period30 day 1 Assumed 1
1 The scale factor α , an indicator of education’s impact, the vaccination rate ρ , saturation incidence coefficients p 1 , 2 are not measured from the actual situation, but estimated from experience. The duration of temporary immunity τ 1 and τ 2 obtained by inoculation and recovery from infection may be different, but for the convenience of numerical simulation, we assume that τ 1 = τ 2 .
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhou, D.; Shi, X.; Zhou, X. Dynamic Analysis of a Stochastic Delayed SEIRS Epidemic Model with Lévy Jumps and the Impact of Public Health Education. Axioms 2023, 12, 560. https://doi.org/10.3390/axioms12060560

AMA Style

Zhou D, Shi X, Zhou X. Dynamic Analysis of a Stochastic Delayed SEIRS Epidemic Model with Lévy Jumps and the Impact of Public Health Education. Axioms. 2023; 12(6):560. https://doi.org/10.3390/axioms12060560

Chicago/Turabian Style

Zhou, Dan, Xiangyun Shi, and Xueyong Zhou. 2023. "Dynamic Analysis of a Stochastic Delayed SEIRS Epidemic Model with Lévy Jumps and the Impact of Public Health Education" Axioms 12, no. 6: 560. https://doi.org/10.3390/axioms12060560

APA Style

Zhou, D., Shi, X., & Zhou, X. (2023). Dynamic Analysis of a Stochastic Delayed SEIRS Epidemic Model with Lévy Jumps and the Impact of Public Health Education. Axioms, 12(6), 560. https://doi.org/10.3390/axioms12060560

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop