Threshold Dynamics of Stochastic Cholera Epidemic Model with Direct Transmission

This paper extends the cholera human-to-human direct transmission model from a deterministic form to a stochastic framework. This is expressed as mixed system of stochastic and deterministic diﬀerential equations. A Lyapunov function is created to investigate the global stability of the stochastic cholera epidemic, which shows the existence of global positivity of the solution using the theory of stopping time. We then ﬁnd the threshold quantity of the extended stochastic cholera epidemic model. We derive a parametric condition (cid:2) R 0 , and for additive white noise, we establish suﬃcient conditions for the extinction and the persistence of the cholera infection. Finally, for a suitable choice of the parameter of the system for (cid:2) R 0 , we perform numerical simulations for both scenarios of extinction and persistence of the dynamic of the cholera infection.


Introduction
Infectious diseases significantly impact human health since they can have terrible effects.The dynamical behavior of infectious diseases must be understood in order to safeguard human health and control an infection [2].Cholera, a diarrhea illness caused by Vibrio Cholera and characterized by diarrhea, is one of the most severe infectious diseases [3].WHO has identified cholera as a public health issue [6].Numerous theoretical and clinical research have been investigated in the literature.Cholera has been documented [7], but it still poses a serious hazard to public health in developing nations.
Mathematical modeling is a powerful tool for understanding how disease spreads [8].Additionally, it provides future outbreak predictions and epidemic control strategies [10].A disease's epidemic model may be deterministic or stochastic.An easy technique to analyze a system of ordinary differential equations is to use a deterministic model [11].Due to their straightforward formulation, these models are commonly utilized.Deterministic models, however, have a number of limitations, including the fact that when the population is small enough, they are less informative and more challenging to estimate and analyze.Additionally, random effects disrupt deterministic models, causing trajectories to deviate from expected noise behavior [12].Deterministic models have their strengths and can be useful for systems with well-defined rules and clear cause-and-effect relationships, but they often fall short in capturing the complexity, randomness, and uncertainty present in many real-world systems.Stochastic systems offer a more flexible and realistic way to model such complexity, making them a valuable tool in various scientific, engineering, and social applications [13,14].
Various techniques can be used to study stochastic epidemic systems [15].The time and state variables show that these processes are different from one another [16].The literature has a variety of stochastic models that investigate how environmental noise affects the spread of infectious diseases.The stochastic models are efficient at calculating asymptotic expressions for the likelihood of an outbreak occurring [17].One of the traditional models is a stochastic SIS epidemiological model.[18] presents such a model whose basic characteristics with vaccination are studied in [19].By employing the threshold quantity of the deterministic system, the authors of [20] investigated a stochastic SIS epidemic model's dynamics and produced conditions for the noise of disease persistence and extinction.Some authors studied the long-term behavior of an SIR stochastic epidemic model [21].They established a threshold condition for the extinction and persistence of the model and supported their theoretical predictions through numerical simulations.[22] used a stochastic SIRS model for studying the noise impact of an infectious disease with a saturated incidence rate and [23] reported a SIQS model for investigating the threshold dynamics for white noise.
Numerous scholars reported different mathematical models to study the cholera epidemic [24][25][26].Depending on the environmental noises, researchers studied a cholera stochastic system with vaccination and derived significant conditions regarding the basic threshold number [27].The role of the aquatic reservoir on cholera disease was reported in [29], where they studied the endemic and epidemic dynamics of the infection.The impact of factors like contaminated water, temperature, rainfall, etc., on the cholera outbreak was investigated through a stochastic model system [30].However, this work is solely concerned with theoretical considerations, and no analysis has been performed.More recently, authors reported a fundamental deterministic model that considers the direct contact transmission mechanism and the dynamics of the cholera sickness [31].The model they used for the said investigation has the following structure: A description of the parameters in system (1), for t > 0 is given below.The term S is the compartment of the susceptible individuals, I is the compartment of the infected populace, T denotes the class of treated people, and R is the class of the recovered class.
The birth rate of the population is represented by Π while the natural rate of mortality is symbolized by ζ .δ denotes the contact rate between the susceptible class and the infected compartment; μ is the rate of treatment while death due to infection is denoted by α while deaths during therapy are described by γ and η is the rate of recovery, the loss rate of immunity is Λ.
Epidemic deterministic models assume even mixing in large populations, but this fails at outbreak start due to few infected people and stochastic transmission.Homogeneous mixing is not accurate for small initial infections, making deterministic models unsuitable.To address this, we introduce random parameter variation using a parametric perturbation approach, reflecting realworld heterogeneity and contact patterns.This enhances our understanding of disease dynamics during the early stages of epidemics.In this work, we extend the notion of the cholera deterministic model presented in [31] into a stochastic problem by taking direct human-to-human transmission into account.For this purpose, we apply the idea of parametric perturbation by choosing parameters from model (1) and transforming it to a random variable [13,22,23].For more results, we refer to the following collection of articles in the book Disease Prevention and Health Promotion in Developing Countries [35].
We outline the environmental changes that affect the parameter δ such that δ → δ + εdW (t).Here, the standard Brownian motion with intensity σ 2 > 0 is denoted by the symbol W (t). Suppose (Ω, F, P) represents an entire probability space with filtration so that {F} t≥0 meets the prerequisites (i.e., F 0 , including all null sets, is assumed to be continuous from the right and growing).Our stochastic model extends the deterministic one presented in [5] and up to the random excitations, it keeps similar extensions as in [4,27].We study the following stochastic model with identic stochastic perturbation for the classes S and I: (2) We will utilize the Lyapunov analysis to determine whether the solution of system ( 2) is stable globally.We first find the corresponding reproduction number in order to analyze the stochastic cholera epidemic dynamics.The calculated threshold number will provide sufficient conditions for the infection to persist and eventually die.In order to validate the analytical results, we will also carry out a few numerical simulations.

Preliminaries
Introduce the notation R 4 + = {(y 1 , y 2 , y 3 , y 4 )|y i > 0, i = 1, 2, 3, 4}.Suppose that (F, Ω, P) illustrates a probability measure space and the filtration {F} t≥0 observing the necessary relations (that is, all null sets include in F 0 , and is continuously increasing).Take into account the underlying m-dimensional stochastic system f and g are locally Lipschitz functions in χ.The differential operator L define by [16,32] associated with (3) as follow Once the operator L is applied to the function ∂y i ∂y j .So we can define Itô formula for y(t) ∈ R m as dN(y, t) = LN(y, t)dt + N y (y, t)g(y, t)dW t .
3 Qualitative analysis of the system (2) The first question to ask when examining dynamic behavior is whether the solution exists globally.Additionally, a populace dynamics system takes into account whether the solution is non-negative.So, in this part, we begin by demonstrating that the solution of model ( 2) is positive as well as global.Locally Lipschitz coefficients ensure a unique local solution for an ordinary differential equation on a given interval by guaranteeing a bounded rate of change for the solution within a neighborhood around each point.This controlled growth prevents solutions from diverging and overlapping, leading to a distinct trajectory for each initial condition, thereby ensuring uniqueness on the interval [28].For this reason, the coefficients of the equation typically need to meet the local Lipschitz condition and the linear growth condition in order to produce a stochastic differential equation that has a distinct global solution for a given initial value [14] .Even if they are locally Lipschitz continuous, the coefficients of the model (2) do not meet the linear growth condition, hence the solution of the system (2) may blow up at a certain point in time.This part of the paper demonstrates that the solution of the system (2) is positive and global utilizing the Lyapunov analysis approach.
Theorem 1.The distinct solution of the stochastic model (2) on t greater than or equal to zero is unique and for a starting point (S 0 , T 0 , I 0 , R 0 ) ∈ R 4 , with probability one will stay in R 4 , specifically (S, T, I, R) ∈ Ω for positive t almost surely (a.s).
Proof.Since the parameters of the system (2) meet the local Lipschitz conditions, then there exists a distinct local solution (I, T, S, R) for an assigned starting value given as (S 0 , T ( 0 , I 0 , R 0 ) ∈ R 4 on [0, τ e ) a.s., with τ e indicates the time of outbursts.Furthermore, to ascertain that solution is global, we require to verify τ e = ∞ as sure.Assume that j 0 > 0 is sufficient large satisfies (I 0 , T 0 , S 0 , R 0 ) remain in the interval [ 1 j 0 , j 0 ].For each integer j ≥ j 0 we describe the time of stopping τ j = inf{ min(T, I, R, S; t ∈ [0, τ e )) ≤ 1 j or max((T, I, R, S) ≥ j}.
In this paper, ∅ represents the empty set such that inf ∅ = ∞.By definition τ j is rising as j → ∞.Take τ ∞ = lim j→∞ τ j , from which τ e ≥ τ ∞ a.s.Now we require to prove that τ ∞ = ∞ a.s, therefore τ e = ∞ and (I, T, S, R) ∈ R 4 a.s.Furthermore, to accomplish this task, we must prove that τ ∞ = ∞ a.s.Suppose that this assumption is false, then one may find some constants B > 0 and any 0 < < 1 for which P{τ ∞ ≤ B} > .As a significance of this, there exist an integer j 1 ≥ j 0 in such a manner Clearly when t ≤ τ j , we can write which can be further written as Now, let us define a Lyapunov function V from R 4 to R as follows: This function is obviously positive as − ln u − 1 + u ≥ 0, ∀ u > 0. Suppose that j ≥ j 0 and for arbitrary B > 0, applying the Itô Integral formula, one may arrive at According to the definition of the operator L, we have After further simplification, one may write Inserting this into (4) leads to Integrating ( 6) from 0 to τ j ∧ B one may arrive at where τ j ∧ B = min(τ n , t).After taking expectations on both sides, it becomes Put Ω j = {τ j ≤ B} for j ≥ j 1 and by (4), P {Ω k } ≥ .It can be noticed that for each ω ∈ Ω j , there corresponds at least one of S(τ j , ω), I(τ j , ω), T (τ j , ω), R(τ j , ω) and that is equal to j or 1 j , and thus V (S(τ j , ω), I(τ j , ω), T (τ j , ω), R(τ j , ω)) which must not minimum than −1 + j − ln j or 1 j − ln 1 j + 1 = 1 j − 1 + ln j.Subsequently one may arrive at Next, it is easy to write from Eqs.( 4) and ( 7) that Note that 1 Ω k is the indicator function of Ω k .Suppose that j → ∞ induces to the contradiction Thus accordingly, it is necessary that

Extinction of disease
For the dynamical behavior in epidemiology, our major concern is how to regulate the illness spread so that the infectious illness vanishes in a long term.For this purpose, the investigation of the sufficient conditions of disease extinction is performed in this section.
To establish sufficient conditions for the cholera disease extinction, we introduce the new notation in this paper: let x(t) = 1 t x(r)dr.
For our planned stochastic model ( 2), the basic reproduction value R0 is described as Proof.By using the Itô integral formula to our proposed model ( 2), then Integrating ( 9) from zero to t and then take ratio of t, we obtain where the non-discontinuous local martingale with property M (0) is equal to zero, represented by Consider Lemma 1, which can be cast into lim t→0 M (t) t = 0, a.s.
By applying the limit superior on both sides of (10)

Case (b):
Integrating system (2) gives the underlying system Further simplification leads to After some manipulation, the last equation may written in the form Consequently, we arrive at Note that the value of ϕ(t) is given by It is apparent that as t → ∞ so as ϕ(t) → 0 a.s.
Now we use the Itô integral formula [28] to the model ( 2), then it attain the form Integrating from zero to t, we reach the expression Dividing both sides by t to obtain This can be further written as where M (t) = ε t 0 I(r)S(r)dW (r) is a continuous and local martingale which satisfies the condition M (0) = 0, see, e.g., [28].By putting ( 12) in ( 15), one may obtain Note that the term θ(t) in ( 16) is given by Next, by taking limit superior on both sides of ( 16) leads to With the help of the above defined condition (b), inequality (18) gives This gives lim t→0 sup ln I t < 0 a.s., as a result of which we deduce that lim t→0 I = 0 a.s.

Persistence of disease
This section of our manuscript is devoted for the construction of important relations for the disease dispersed in the populace.We follow the techniques discussed in [18,20,21].The planned stochastic system (2) stated that persistence in means if lim t→0 inf t 0 I(r)dr > 0, a.s.
To proceed further, we consider it necessary to state the following lemmas: in case one can find non-negative constants π 0 , π and M satisfies , so for any starting value (S 0 , T 0 , I 0 , R 0 ) ∈ η * , the solution (S, T, I, R) of the planned cholera disease system (2) satisfies the property X 1 ≥ lim inf I ≥ X 2 very nearly surely, where .
Proof.From the last inequality of ( 16), we may write Using the value of R from (8), the last result takes the form: Upon some algebraic manipulation, (19) yields: Taking limit superior on both sides and using Lemma 2 together with (13), we arrive at: On the other hand, from ( 12) and ( 14), one may write Thus, we have With the help of ( 8), the last in-equality gives ln I − ln After some algebra, one may observe that (22) leads to the following equation: Applying the inferior limit to each side and utilize Lemma 3 with ( 16) we can write lim Thus from ( 12) and ( 16), we have

Numerical simulation
In order to validate our theoretical results, we perform three numerical tests for the present SIT R stochastic system (2).This coupled stochastic model is derived from the classic SIT R model (1).Thus, we use a random excitation of the population classes S and I.For computational simplicity, we use white noise as a random perturbation.These noises are independent.System (2) is discretized using the first-order stochastic Runge Kutta scheme.For such a stochastic framework, we refer to [33,33].The derivation of this scheme is given as follows: Where Δt n = t n+1 − t n represents the time step and ΔW tn = W t n+1 − W tn stands for the independent increments of the Gaussian Brownian motion generating the white noise.For simplicity, we implement a fixed time step Δt n = Δt for the evolution of all classes.We subdivide the temporal interval into 1,000 equally spaced subintervals.Subject to various random initial conditions that satisfy our theoretical requirements, we numerically solve the SIT R system (2).We run our code in order to generate six simulations.The parameters of these six tests are given in table 1.For tests 4., 5. and 6, we use the same technique for the parameters σ as implemented in [9,34].Moreover, it should be stressed that the choice of our parameters proved in theorems 2 and 3 is relevant for checking the extinctions and persistence scenarios.In tests 1,2 and 3, we simulate the extinction case and in tests 4,5 and , we simulate the persistence case.Obviously, all these parameters follow the conditions of our proven results.For every simulation, we show two samples of numerical realization and the associated numerical mean solution of a large number of solutions.Namely, we generate 1,000 solutions for each test.
Again suppose that the Jacobian matrix of f and v w.r.tI are F and V respectively,which can be write at the no-infection equilibrium point E 0 = ( Π ζ , 0, 0, 0) as, From the last two expressions, it is easy to write The eigenvalue can be found by the Characteristics equation After solving, we get the reproduction number of stochastic model as, .
s., with the help of Lemma 1 can be written as lim t→0 M (t) t = 0 and lim t→0 θ(t) = 0 a.s.

Lemma 1 .
Suppose that a continuous real-valued martingale denoted by M = {M t } t≥0 is disappearing at positive t.So lim < ∞ almost surely, interprets that lim t→0 Mt t = 0 almost surely.Theorem 2. Let (I, S, T, R) be the solution of the planned system (2) with starting value
The initial values of all classes are represented by the values S(0), I(0), T (0), and R(0).