Analysis of Stochastic COVID-19 and Hepatitis B Co-infection Model with Brownian and Lévy Noise ================================================================================================ * Michael A. Pobbi * Stephen E. Moore ## Abstract In this article, we formulate and analyze a mathematical model for the coinfection of HBV and COVID-19 that incorporates the effects of Brownian and Lévi noise. We studied the dynamics and effects of these diseases in a given population. First, we establish the basic reproduction number of the disease-free equilibrium point of the stochastic model by means of a suitable Lyapunov function. Additionally, we provided sufficient conditions for the stability of the model around the disease-free equilibrium points. Finally, using a few simulation studies, we demonstrate our theoretical results. In particularly, we derived threshold values for HBV only ![Graphic][1], COVID-19 only, ![Graphic][2], and coinfection ![Graphic][3]for the stochastic model around disease-free equilibrium point. Next, the conditions for stability in the stochastic sense for HBV only, COVID-19 only submodels, and the full model are established. Furthermore, we devote our concentrated attention to sufficient conditions for extinction and persistence using each of these reproductive numbers. Finally, by using the Euler–Murayama scheme, we demonstrate the dynamics of the coinfection by means of numerical simulations. Keywords * COVID-19 * Hepatitis * Stochastic model * Lévi noise * Global Stability * Local Stability ## 1 Introduction Both Coronavirus Disease 2019 (COVID-19) and Hepatitis B virus (HBV) infections remain major public health threats [4]. COVID-19 is caused by Severe Acute Respiratory Syndrome Corona virus 2 (SARS-CoV-2) and affect the functionality of the respiratory system. It is spread from one person to another by coming into direct contact with respiratory droplets from an infected individual. It can also spread by touching surfaces contaminated with the virus and then touching the face. On the other hand, HBV disease is caused by the hepatitis B virus and is associated with serious liver infection, leading to liver malfunction, including cirrhosis, liver fibrosis, and hepatocellular carcinoma. Both COVID-19 and HBV are life-threatening diseases that have been declared global health problems [18]. For example, [32] reported that more than 3 billion people are infected with HBV globally. An estimated 7 million COVID-19 deaths and 700 million infections have been reported since December 2019 [18]. These statistics demonstrate how devastating HBV and COVID-19 diseases are worldwide. Presently, global efforts have been made to develop and distribute COVID-19 vaccines, aiming to reduce severe illness, hospitalizations, and deaths. Over the years, mathematical models have provided frameworks for explaining transmission dynamics and for making necessary interventions for several infectious diseases; see, e.g. [2, 11, 28, 5, 7, 16]. Several models of COVID-19, HBV and other infectious disease transmission dynamics have been developed in the literature. The transmission dynamics of COVID-19, HBV, and other diseases have been extensively studied, mostly in isolation. Both deterministic and stochastic models have been developed to understand the spread of these diseases [13, 11, 9]. In [11], the authors deemed the susceptible-exposed-infectious-recovered (SEIR) model introduced by [34] to describe the transmission dynamics and forecast the national and global spread of COVID-19 disease to be simplistic and largely unrealistic because it excluded the relapse from recovered class to susceptible class. Like [34], authors including [22, 30] also made an assumption of no relapse from the recovered to the susceptible class, an assumption that is critiqued in [28]. [26] examine the extinction and persistence for a HBV only model and also conducted simulations on the model. Amidst inconclusive findings from clinical studies on whether one disease increases the severity of another disease in the event of coinfection (see e.g. in [4, 29]), there has been a spike in interest for both clinical and theoretical scientific research on disease coinfections in recent times. The theoretical study of coinfection through mathematical modelling will be valuable to shaping our understanding of transmission dynamics of these diseases. A number of studies that have been done on co-dynamics of COVID-19 including Dengue [23], cholera [10], malaria, [20, 1, 21, 3, 19] and co-dynamics s TB, see, e.g. [24, 14]. [29] employed a deterministic model and examined the effects of optimal control strategies on the Covid-19, HBV coinfection transmission. On the other hand, [8] examined COVID-19 and HBV coinfection by employing Gaussian noise to explain the perturbation in observed data for the COVID-19 and Malaria coinfection model. Stochastic coinfection models are very scarce in literature, and more so for COVID-19 and HBV Stochastic coinfection models. While mathematical models, see, e.g., [17, 28], have been proposed to describe the dynamics of diseases including COVID-19 and HBV, there is a dearth in studies imprint the effect stochastic perturbations on COVID-19-HBV coinfection. In this paper, we employ both Gaussian noise and Non-Gaussian noise to explain the noise in the Coinfection of COVID -19 and HBV in the presence of COVID-19 Vaccination, and determine conditions for extinction and persistence for the coinfection model. Stochastic coinfection SIRS models unlike classic or deterministic SIRS models, incorporate the effects of random fluctuations in the transmission dynamics of disease and recovery rates, as well as the possibility of coinfection by different diseases. In [35, 15], the authors formulated techniques for examining existence, uniqueness and persistence of solution. Understanding and analysing the dynamics of infectious illness transmission in both humans and animals can benefit from mathematical modelling. The persistence or extinction of the diseases depends on several factors, such as the basic reproduction numbers, the coinfection rates, and the stochastic perturbations; hence, determining the conditions for the persistence or extinction of the diseases is a cumbersome task. One of the methods proposed to analyze the persistence or extinction of stochastic coinfection SIRS models is the use of Lyapunov functions, which are scalar functions that measure the distance from an equilibrium point [33]. For example, [33] employed Lyapunov functions to show that if the basic reproduction numbers of both diseases are less than one, then the disease-free equilibrium is globally asymptotically stable for a stochastic coinfection SIRS model with nonlinear incidence rates and cross-immunity effects. In this paper we develop condition for extinction and persistence for coinfection model. In this article, we will formulate a mathematical model which builds on (2) of [8] by firstly assuming relapse to susceptible state. The new model includes COVID-19 vaccination state, and also assume the presence of both standard random fluctuation and massive disturbances which arise due to sudden environmental shocks in the transmission dynamics. The novel coinfection model will open up some new avenues for further research in stochastic modelling dynamics of coinfection of diseases. Specifically, our study seeks to investigate the following objectives: 1. To develop a compartmental model for explaining the transmission dynamics of HBV-COVID-19 coinfection. 2. To establish the Invariance and Positivity of the stochastic system for coinfection Model. 3. To determine the basic Reproductive number for Stochastic HBV Only, COVID-19 Only and Coinfection models respectively. 4. To determine conditions for stability at Disease Free Equilibrium for Stochastic HBV Only, COVID-19 Only and coinfection models respectively using Brownian Noise and Poisson Lévi Noise. 5. To present numerical simulations and discussions of the model. The rest of the paper is organised as follows; Section 2 discusses the Mathematical Model and the underlying assumptions. In Section 3, we analyze the COVID-19, and HBV sub-models as well as the Coinfection model. Section 4 is devoted to discussing the numerical method and numerical simulation of the model. Finally, we present the conclusion in Section 5. ## 2 Mathematical Model Formulation In this section, we consider a stochastic compartmental mathematical model that uses differential equations, see Figure 1. ![Figure 1:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2024/08/14/2024.08.12.24311861/F1.medium.gif) [Figure 1:](http://medrxiv.org/content/early/2024/08/14/2024.08.12.24311861/F1) Figure 1: Compartmental diagram for the COVID-19 Hepatitis B coinfection cohort. The sub-populations of the human population at time *t*, denoted by *N* (*t*), include susceptible individuals *S*(*t*), COVID-19 vaccinated individuals *V* (*t*), COVID-19 infected individuals *I**C*(*t*), Hepatitis-B infected individuals *I**H* (*t*), individuals infected with both COVID-19 and Hepatitis-B *I**HC*(*t*), and recovered individuals *R*(*t*). We will consider the following assumptions; (*A*1) The total human population (*N* (*t*)) is a sum of all the different compartments, i.e. ![Formula][4] (*A*2) Covid-19 Vaccination wears off after some period of time. (*A*3) There is relapse to susceptible state, (*A*4) There is no vertical transmission. ![Formula][5] with initial conditions *S*(0) *>* 0, *I**C*(0) ≥0, *I**H* (0) ≥ 0, *I**HC*(0) ≥0, *R*(0) ≥0, *V* (0) ≥ 0. *β**C*, *β**H* are the effective contact rate for COVID-19, and HBV respectively. Two popular ways to introduce stochastic factors into epidemic models are: (i) to assume some small and standard random fluctuation, (ii) to assume massive disturbances caused by sudden environmental shocks. The first, is described by Gaussian white noise while the second is described by the Lévy noise. It is assumed that fluctuations in the environment will manifest mainly as fluctuations in the parameter *β*, i.e. *β* → *β* + *ϵ**i**dB**i*(*t*), (*i* = 1, 2, …, 6), where *dB**i*(*t*) is a one dimensional standard Brownian motion with *B**i*(0) = 0, and *ϵ**i* is the intensity of the white noise. Also a multiplicative noise is considered in our model, as random term depends on the state space, e.g. [8]. The Lévy jump process, thus, is divided into linear drift term, Brownian motion and compensated Poisson Process. The stochastic model of the system (2) takes the following form: ![Formula][6] where *dB**i*(*t*) is independent of ![Graphic][7]. 𝒩 is a Poisson counting measure with compensating martingale ![Graphic][8] and characteristic measure *v* on a measurable subset 𝒰 ∈ (0, ∞) where *v*(𝒰) ≤ ∞ and Ψ*i* : 𝒰 → ℜ (*i* = 1, 2, 3, 4, 5, 6) denotes the effects of random jumps and is bounded and continuous. It is assumed that *v* is a Levy measure so that ![Graphic][9]. We note that *S*(*t*−), *I**H* (*t*−), *I**C*(*t*−), *I**HC*(*t*−), *R*(*t*−), *V* (*t*−) denote the left limits of *S*(*t*), *I**H* (*t*), *I**C*(*t*), *I**HC*(*t*), *R*(*t*), *V* (*t*) respectively. The following definitions are necessary for the analysis of the stochastic model. Definition 1 (Itô-Lévi [26]). *First, set* ![Graphic][10]. *Suppose a complete probability space* (Ω, *F, F**t* ≥ 0, *P*) *with filtration F**t* ≥ 0, *which satisfies the usual conditions. B**i*(*t*) *is defined on this probability space*. *We assume an Itô -Lévi process, X*(*t*) ∈ ℜ+, *of the form* ![Formula][11] *where, F* : ℜ*n* × ℜ+ × *S* → ℜ*n*, *G* : ℜ*n* × ℜ+ × *S* → ℜ*n*, *and* Ψ : ℜ*n* × ℜ+ × *S* × *Y* → ℜ*n* *are measurable functions and X*(*t*−) *denotes the left limit of X*(*t*). *F* (*X*(*t*), *t*) *represents the linear drift term, G*(*X*(*t*), *t*) *the Brownian noise*, ![Graphic][12] (*dt, du*) *represents the compensated Poisson random measure, and* Ψ(*X, t, u*) *represents the intensity of jumps. It is assumed that the condition*![Graphic][13] Definition 2 (Itô-Lévi Formula [26]). *We consider the process X expressed by* (4) *and let* 𝒱 ∈ *C*2,1(ℜ*n* × ℜ+ × *S*; ℜ+) *such that* 𝒴 ≡ 𝒱(*t, X*(*t*)). *Then*, 𝒴(*t*) *is again an Itô-Lévi process and* ![Formula][14] *Which is the representation of generalized Itô’s formula with jumps*. ![Formula][15] To gain insights into the underlying mechanisms that drive the spread of infectious diseases to identify the long-term behaviour of disease transmission we conduct qualitative analysis of the stochastic model (3). ## 3 Qualitative Analysis of Model In this section, we present the Positivity of the Solution to (3), determine the Reproductive Number for Single Disease Models and Coinfection model, determine the conditions for Local and Global equilibrium for COVID-19 only, HBV only, and for Coinfection model. Further we determine the conditions for Extinction and Persistence for all three models. ### 3.1 Positivity of Model A first step is to determine the positivity of solution to the stochastic system (3) as a negative state will not be biologically meaningful. We will show the positivity of the first equation in(3). Theorem 1 *Let* ((*S*(*t*), *I**H* (*t*), …, *V* (*t*)) *be a solution of system* (3), *given initial values* ((*S*(0), *I**H* (0), …, *V* (0)) ∈ Ω, *where* Ω ∈ ℜ+, *then* lim*t→∞* *S*(*t*) ≥ 0, lim*t→∞* *I**H* (*t*) ≥ 0, …, lim*t→∞* *V* (*t*) ≥ 0. *Proof*. Let us consider (3). By eliminating terms and rearrangement, we obtain ![Formula][16] where *r* = *β**H* *I**H* (*t*)+*β**C**I**C*(*t*)+*β**HC**I**HC*(*t*) − (*µ*+*θ*). Next, we take a function *g*(*t, S*(*t*)) = ln(*t, S*(*t*)) and apply Itô’s formula to get, ![Formula][17] Combining (6) and (7), solving for *S*(*t*) and taking the limit we obtain, ![Formula][18] Hence *S*(*t*) is positive definite. Similarly, lim*t→∞* *I**H* (*t*) ≥ 0, lim*t→∞* *I**C*(*t*) ≥ 0, lim*t→∞* *I**HC*(*t*) ≥ 0, lim*t→∞* *R*(*t*) ≥ 0, lim*t→∞* *V* (*t*) ≥ 0. ### 3.2 Invariant Region of Stochastic COVID-19 and Hepatitis B Coinfection Model An investigation of the invariant region of the model will help us determine the set of states for which the disease will eventually either persist or die out infinitely. Theorem 2 *Let N(t) be the population of system* (3), *given initial values N* (0) *>* 0, *and* ![Graphic][19], *then* ![Graphic][20], *and* (*S*(*t*), *I**H* (*t*), …, *V* (*t*)) ∈ Ω. *Proof*. By adding and simplifying terms in (3) and choosing *ϵ* = *max*(*ϵ**i*) we obtain; ![Formula][21] Applying (5) we have, ![Formula][22] We let ![Graphic][23]. The complementary solution to (3)) was obtained as ![Formula][24] Next, we determine the particular solution of (8) and obtain the general solution of (8) as ![Formula][25] We then find the lim*t→∞* of the general equation of (8) and apply comparison to obtain ![Formula][26] ### 3.3 Disease Free Equilibrium Point We consider the first equation of (3). For disease free equilibrium, we have *I**C*(*t*) = 0, *R*(*t*) = 0, and let ![Graphic][27], solve for *S*(*t*), and apply comparison to obtain, ![Formula][28] Since 0 *< ϵ*1(*dB*1) + ∫ 𝒰 Ψ1(*u*) *< m <*∞, where m is a positive constant, the final inequality is obtained by comparison. Hence, the disease free equilibrium is given by ![Formula][29] ### 3.4 Reproductive Number of COVID-19 only Model Theorem 3 *Given D*.*F*.*E. point* ![Graphic][30], *the Reproductive Number of COVID-19 only model* ![Graphic][31]. *Proof*. Applying Itô-Lévi formula to the COVID-19 infected class of (3), we obtain ![Formula][32] Considering initial infections *F* and secondary infections *V*, we determine the basic reproduction number by means of the next generation matrix as follows ![Formula][33] Taking the spectral radius of *F V* −1, i.e., *ρ*(*F V* −1), we obtain, ![Formula][34] #### 3.4.1 Local Stability of COVID-19 Only Model Theorem 4 *If* ![Graphic][35], *then for any initial values of* ![Graphic][36], *I**C*(*t*) *satisfies* ![Graphic][37] *Proof*. By integrating both sides of (11) and evaluating at disease free-equilibrium point point, we have, ![Formula][38] where ![Formula][39] are a martingales given ![Graphic][40]. Dividing through by t and taking lim sup*t→∞* of both sides, we get, ![Formula][41] The final inequality holds if ![Graphic][42], which implies,![Graphic][43]. #### 3.4.2 Global Stability of COVID-19 Only Model Theorem 5 *If* ![Graphic][44], *then E**C* *is globally asymptotically stable in* Ω. *Proof*. We consider the Lyapunov function, ![Formula][45] For ![Graphic][46] *and* **Q** =(*γ* + *μ* + *ρ**C*) ![Formula][47] since ![Graphic][48]. Hence *L*1 is globally stable in the domain Ω. ### 3.5 Basic Reproduction Number for HBV Model We consider only HBV transmission in human. The model of the system is given as (3). Theorem 6 *Given D*.*F*.*E. point* ![Graphic][49] *then the Reproductive Number of HBV only Model* ![Graphic][50]. *Proof*. Applying Itô-Lévi formula to the HBV infected class of (3), we obtain, ![Formula][51] Considering initial infections *F* and secondary infections *V*, we determine the basic reproduction number by means of the next generation matrix as follows; ![Formula][52] Taking the spectral radius of *F V* −1, i.e., *ρ*(*F V* −1), we obtain, ![Formula][53] #### 3.5.1 Local Stability of HBV Only Model Theorem 7 *If* ![Graphic][54], *then for any initial values of* ![Graphic][55] *satisfies* ![Graphic][56]. *Proof*. We integrate both sides of (12) and evaluate at D.F.E. point to obtain, ![Formula][57] Next, we divide through by t, find lim sup *→ ∞* and apply the strong law of Martingales as follows, ![Formula][58] which completes the proof. ![Graphic][59] as (*α* + *µ* + *ρ**H*) *>* 0, which implies, ![Graphic][60] where ![Graphic][61], and ![Graphic][62] are a martingales given ![Graphic][63]. #### 3.5.2 Global Stability of HBV Only Model Theorem 8 *If* ![Graphic][64], *then E**H* *is globally asymptotically stable in* Ω *Proof*. We consider the Lyapunov function *L*2, ![Formula][65] For ![Graphic][66],and **P** = (*α*+ *μ*+ *ρ**H*). ![Formula][67] hence *L*2 is globally stable in the domain Ω. ### 3.6 Reproductive Number of Coinfection Model Theorem 9 *Given D*.*F*.*E. point* ![Graphic][68], *the Reproductive Number of Coinfection* ![Graphic][69]. *Proof*. Applying Itô-Lévi formula to the Coinfected equation of (3), we obtain, ![Formula][70] Considering initial infections *F* and secondary infections *V*, and since ![Graphic][71] ![Formula][72] Taking the spectral radius of *FV* −1, i.e., ![Graphic][73] we obtain, ![Formula][74] #### 3.6.1 Local Stability of Disease Free Equilibrium point in the Coinfection Model Theorem 10 *If* ![Graphic][75], *then for any initial values of* (*S*(0), *I**H* (0), *I**C*(0), *I**HC*(0), *R*(0), *V* (0)) ∈ ℜ6, *I**HC*(*t*) *satisfies* ![Graphic][76]. *Proof*. Integrating (14), eliminating terms, choosing constants, and setting **G** = (*η* + *µ* + *ρ**HC*) yields, ![Formula][77] Next, we divide through by t, find lim sup and apply the strong law of Martingales as follows, ![Formula][78] The final inequality holds if ![Graphic][79] as (*η* + *µ* + *ρ**C*) *>* 0, which implies, ![Graphic][80] where ![Graphic][81], and ![Graphic][82] are a martingales given ![Graphic][83]. #### 3.6.2 Global Stability of Coinfection Model Theorem 11 *If* ![Graphic][84], *then E**HC* *is globally asymptotically stable in* Ω *Proof*. We consider the Lyapunov function, Net we differential and employ comparision to determine conditions for Global Stability of Coinfection. ![Formula][85] hence *L*3 is globally stable in the domain Ω. ### 3.7 Extinction of HBV-COVID-19 Coinfection disease Next, we determine conditions under which the disease will eventually die out in the population with a probability of one. In this section we study the conditions of extinction for the coinfection model. Definition 3 *[8] We define* ![Graphic][86], *and* ![Graphic][87]. Lemma 1 *[25] Let* (*S*(*t*), *I**C*(*t*), *I**H* (*t*), *I**HC*(*t*), *R*(*t*), *V* (*t*)) *be the positive solution of system*(4) *with given initial condition* ![Graphic][88], *Let also X(t) be the positive solution of equation* (4) *with given initial condition given condition X*(0) = *N* (0) = *S*(0) + *I**C*(0) + *I**H* (0) + *I**HC*(0) + *R*(0) + *V* (0) ∈ ℜ+. *Then* 1. ![Graphic][89] 2. ![Graphic][90] 3. ![Graphic][91],*and* 4. ![Graphic][92] Theorem 12 *Let* (*S*(*t*), *I**H* (*t*), …, *V* (*t*)) *be the solution of* (3) *with initial values* (*S*(0), *I**H* (0), …, *V* (0) ∈ Ω), *the coinfection disease of model* (3) *goes extinct almost surely* (lim*t→∞* *I**HC*(*t*) = 0) *a*.*s if one of the following assumptions holds:* 1. ![Graphic][93] *and* ![Graphic][94], 2. ![Graphic][95] *and* ![Graphic][96], 3. ![Graphic][97]. *Proof*. By integrating (3), dividing by t, adding terms, followed by some algebraic manipulations we obtain the following equations: ![Formula][98] ![Formula][99] where ![Formula][100] Next, we integrate both sides (14). Since ![Graphic][101] we get, ![Formula][102] We divide through (17) by *t*. Next we substituting (15) and (16) into (17) we obtain, ![Formula][103] ![Formula][104] Next, we find lim sup*t→∞* and apply a Lemma 2 and Lemma 1 to obtain following results, ![Formula][105] ![Formula][106] hence ![Formula][107] Similarly, we can prove ![Graphic][108] when ![Graphic][109], and ![Graphic][110] when ![Graphic][111]. ### 3.8 Persistence in mean Persistence of a disease means that the disease will remain endemic in the population with positive probability. We study the disease persistence for the system reported in (3) and derive that the disease persists under certain some conditions. Definition 4 *[8] Persistent for Model* (3) *will hold if*, ![Formula][112] Lemma 2 *[6] Set g* ∈ ℂ [[0, ∞] × Ω, (0, ∞)], *assume that there exist ξ*, *ξ >* 0 *such that* ![Formula][113] *such that G* ∈ (ℂ[[0, ∞] × Ω, (0, ∞)]) *satisfying* ![Graphic][114]. *Then* ![Graphic][115]. Theorem 13 *Let* (*S*(*t*), *I**H* (*t*), *I**C*(*t*), *I**HC*(*t*), *R*(*t*), *V* (*t*)) *be a solution of system* (3) *with initial values* (*S*(0), *I**H* (0), *I**C*(0), *I**HC*(0), *R*(0), *V* (0)) ∈ Ω. 1. *If* ![Graphic][116] *and* ![Graphic][117], *the disease I**H* *persists in mean. In addition, I**H* *holds* ![Formula][118] 2. *If* ![Graphic][119], max (![Graphic][120]*and* ![Graphic][121]) *<* 1, *then the disease I**C* *is persistent in mean. In addition, I**C* *satisfies* ![Formula][122] *Proof*. We begin with the first statement of the theorem. By rearranging terms in (15), we get: ![Formula][123] Next, we integrate both sides of (12), and eliminate terms to obtain, ![Formula][124] We divide (24) by t, substitute (23) into (24), and find lim inf*t→*+*∞* to get, ![Formula][125] Then, by Lemma 2 and since ![Graphic][126], we conclude, ![Formula][127] which concludes the proof. Next we prove the persistence in mean for COVID-19. Let us integrate both sides of (11), and then eliminate terms to obtain, ![Formula][128] We divide (25) by t, substitute (23) into (25), and find lim inf*t→*+*∞* to get, ![Formula][129] Then, by Lemma 2 and since ![Graphic][130], we conclude, ![Formula][131] Finally in a similar way, we can prove persistence for coinfection of diseases; ![Formula][132] ## 4 Numerical Results In this section, we employ a standard numerical procedure In this section, we employ a standard numerical procedure to obtain numerical simulation results for the HBV-only, COVID-19-only, and HBV-COVID-19 coinfection models. The objective here is to verify the theoretical results obtained earlier in this paper. The parameter values used along with their sources are displayed in Table 2. A few values were assumed. As we sought to verify theoretical results, we selected initial values of *S*(0) = 3000, *I**H* (0) = 10, *I**C*(0) = 10, *I**HC*(0) = 10, *R*(0) = 0, *V* (0) = 40 The Euler-Maruyana scheme [12] was employed in conducting the numerical simulations for all models. The results from numerical simulation for HBV only, COVID-19 only and Confection model are presented in Figures 2 - 8. View this table: [Table 1:](http://medrxiv.org/content/early/2024/08/14/2024.08.12.24311861/T1) Table 1: Parameters and description View this table: [Table 2:](http://medrxiv.org/content/early/2024/08/14/2024.08.12.24311861/T2) Table 2: Parameters Values of Source ![Figure 2:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2024/08/14/2024.08.12.24311861/F2.medium.gif) [Figure 2:](http://medrxiv.org/content/early/2024/08/14/2024.08.12.24311861/F2) Figure 2: Solution for HBV Only Model. ### 4.1 Hepatitis B Virus Only Model By using fixed values of all parameters shown in Table 2 we plot the curves representing the dynamics of the susceptible, infectious, and recovered classes in the HBV-only model (3). The graphs of the compartments are presented in Figure 2 reveal. From the parameter values in Table 2, we obtained our basic reproduction numbers ![Graphic][133]. We observe initially, from Figure 2 that as number of infectious population increase in time, the susceptible population decrease. However, in the long run, the solution curves for both the susceptible and infective compartments behave like decreasing functions and approach a fixed point. A decrease in the number of susceptible individuals occurred within a relatively short period of the epidemic, after which there was no marked change over time. On the other hand, the decline in the number of infective compartments tends to increase but eventually declines to a fixed point. Similarly, the recovered class also approach a fixed point in the long run. In other words, all compartments will reach to the respective endemic fixed point of model (3) within a finite time and in long run. The numerical outcomes verifies the theorem statement that when ![Graphic][134], then the disease free equilibrium is locally and globally asymptotically stable. The impact of various parameter on the transmission dynamics of HBV is explored. Numerical results are presented in Figure 3 and Figure 4. We proceeded by first investigating the effects of effective contact rate on population of infected individuals. This was done using selected values of effective contact rate (*β**H* = 0.001, 0.002, 0.003), while all remaining parameter values of the model were held constant. From figure 3, we observe that as values of contact rate increases, there is an increased possibility for the population to be infected by HBV, and vice versa. ![Figure 3:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2024/08/14/2024.08.12.24311861/F3.medium.gif) [Figure 3:](http://medrxiv.org/content/early/2024/08/14/2024.08.12.24311861/F3) Figure 3: HBV Infection with varied effective contact rate. ![Figure 4:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2024/08/14/2024.08.12.24311861/F4.medium.gif) [Figure 4:](http://medrxiv.org/content/early/2024/08/14/2024.08.12.24311861/F4) Figure 4: Recoveries with varied effective recovery rate. Next, we investigated the effects of recovery rate on population of Recovered individuals (*R*(*t*)). Similarly, we held all other parameters constant, and selected varied recovery rate (*ρ**H* = 0.05, 0.06, 0.07) for this investigation. We observe Figure 4, that more individuals recover with increasing value in *ρ**H*, thus, the reducing the number of infectious individuals. From this we can conclude that by treating the infectious individuals the infected population goes to recovered state and HBV is eliminated from the community. ### 4.2 COVID-19 Only Model In Figure 5 simulation results representing the dynamics of Susceptible, Infectious, Vaccinated and Recovered classes in the COVID-19 only model of(3) are displayed. Figure 5 is plotted using combination of all compartments and by taking fixed values of parameters in Table 2. We observe initially, from Figure 5 that as number of infectious population increase in time, susceptible population sharply decrease. However, in the long run, the solution curves for both Susceptible and infective compartments behave like decreasing functions and approach reach to a fixed point. The decline in the number of infective compartment occurs during the after the twentieth time step of the epidemic, after which there is no remarkable change over time, but for fluctuations around the deterministic solution. ![Figure 5:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2024/08/14/2024.08.12.24311861/F5.medium.gif) [Figure 5:](http://medrxiv.org/content/early/2024/08/14/2024.08.12.24311861/F5) Figure 5: Numerical Solution for COVID-19 Only Model. Similarly, the remaining classes, the vaccinated, and recovered populations also approach fix points in the long run. Thus, all compartments will reach to the respective endemic fixed point of the COVID-19 only model in (3) within a finite time and in long run. The numerical outcomes verifies the theorem’ statement that when ![Graphic][135], then the disease free equilibrium is locally and globally asymptotically stable. The impact of various parameter on the transmission dynamics of COVID-19 was explored. Numerical results are presented in Figure 6 and Figure 7. We proceeded by first investigating the effects of effective contact rate on population of infected individuals. This was done using selected values of effective contact rate (*β**C* = 0.005, 0.006, 0.007), while all remaining parameter values of the model were held constant. Figure 6, shows that as the contact rate increases, the probability that the population will be infected by diseases increases, and vice versa. ![Figure 6:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2024/08/14/2024.08.12.24311861/F6.medium.gif) [Figure 6:](http://medrxiv.org/content/early/2024/08/14/2024.08.12.24311861/F6) Figure 6: COVID-19 Infection for different effective contact rates. ![Figure 7:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2024/08/14/2024.08.12.24311861/F7.medium.gif) [Figure 7:](http://medrxiv.org/content/early/2024/08/14/2024.08.12.24311861/F7) Figure 7: COVID-19 recoveries for different treatment rates. Next, we examine the effects of recovery rate on population of COVID-19 recoveries *R**C*(*t*). Again, while we held all other parameters constant, we used various recovery rate (*ρ**C* = 0.3, 0.4, 0.5). We observe in Figure 7, that more individuals recover with increasing value in *ρ**HC*, thus, the reducing the number of COVID-19 infections in the population. ### 4.3 Coinfection Model Figure 8 reveals that the curves representing the dynamics of Susceptible, Infectious, Vaccinated and Recovered classes in the coinfection model (3). Figure 8 is plotted using fixed values of all parameters from Table 2 we obtained our basic reproduction numbers ![Graphic][136]. ![Figure 8:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2024/08/14/2024.08.12.24311861/F8.medium.gif) [Figure 8:](http://medrxiv.org/content/early/2024/08/14/2024.08.12.24311861/F8) Figure 8: Numerical Solution for coInfection Model. We observe initially, from Figure 8 that as number of infectious population increase in time, susceptible population will decrease, however, in the long run, the solution curves for Susceptible and infective compartments *I**C*, *I**H*, and *I**HC* behave like decreasing functions and approach reach to a fixed point over time, but for oscillation around their deterministic path due to the presence of the noise. Similarly, the all remaining classes either decrease to zero or also approach to the fixed points in the long run. In other words, all compartments will reach to the respective endemic fixed point of model (3) within a finite time and in long run. The numerical outcomes verifies Theory 10 which states that when ![Graphic][137], then the disease free equilibrium is locally and globally asymptotically stable. The impact of various parameter on the transmission dynamics of COVID-19 -HBV coinfection is explored. Numerical results are presented in Figure 9 and Figure 10. We proceed by first investigating the effects of effective contact rate on population of infected individuals. This was done using selected values of effective contact rate (*β**HC* = 0.001, 0.002, 0.003), while all remaining parameter values of the model were held constant. From figure 9, we observe that as values of contact rate increases, there is an increased chance for the population to be infected by both diseases, and vice versa. ![Figure 9:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2024/08/14/2024.08.12.24311861/F9.medium.gif) [Figure 9:](http://medrxiv.org/content/early/2024/08/14/2024.08.12.24311861/F9) Figure 9: Coinfection for different effective contact rates. ![Figure 10:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2024/08/14/2024.08.12.24311861/F10.medium.gif) [Figure 10:](http://medrxiv.org/content/early/2024/08/14/2024.08.12.24311861/F10) Figure 10: Coinfection recoveries for different treatment rates. Next, we examine the effects of recovery rate on population of Coinfected recoveries *R**HC*(*t*). Again, while we held all other parameters constant, we made use of varied recovery rate (*ρ**HC* = 0.1, 0.15, 0.2). We observe in Figure 10, that more individuals recover with increasing value in *ρ**HC*, thus, the reducing the number of coinfectious in the population. ## 5 Conclusion The majority of issues in the real world are not deterministic. Because of its proximity to ambient sounds, the stochastic models allow us to model epidemic diseases in a more realistic manner. In this paper, two distinct diseases were used to study a stochastically perturbed SIRS-type model. Emphasis was given to explaining the effect of small fluctuations and perturbations due to sudden environmental shocks in the coinfection transmission dynamics, hence, the inclusion of the white noise and the Lévy noise. To determine the conditions for the local asymptotic stability of the equilibria, we first computed the reproductive number and the equilibria for the underlying stochastic model (3). Furthermore, the conditions of extinction and persistence of HBV only, COVID-19 only, and coinfection model we investigated, and it was shown that persistence for all models were dependent on the intensity of the white noise and the values of epidemic parameters involved in disease transmission. Additionally, we demonstrated the influence of different parameters on the infectious compartments and demonstrated a broad variety of simulation results for all three stochastic models. ## Data Availability All data produced in the present work are contained in the manuscript ## Footnotes * This version of the manuscript has been revised to update the following: i) Number of Authors ii) definition 1 and 2 iii) R\_{0C}, R\_{0H}, R_{0HC} notations in Theorem 12 iv) inclusion of two new citations. * Received August 12, 2024. * Revision received August 14, 2024. * Accepted August 14, 2024. * © 2024, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NoDerivs 4.0 International), CC BY-ND 4.0, as described at [http://creativecommons.org/licenses/by-nd/4.0/](http://creativecommons.org/licenses/by-nd/4.0/) ## References 1. [1]. A. I. Abioye, O. J. Peter, H. A. Ogunseye, F. A. Oguntolu, T. A. Ayoola, and A. O. Oladapo. A fractional-order mathematical model for malaria and COVID-19 co-infection dynamics. Healthcare Analytics, 4:100210, 2023. 2. [2]. L. J. Allen. An introduction to stochastic processes with applications to biology. CRC press, 2010. 3. [3]. P. Boonyarangka, K. Phontham, S. Sriwichai, K. Poramathikul, K. Harncharoenkul, W. Kuntawunginn, N. Maneesrikhum, S. Srisawath, C. Seenuan, C. Thanyakait, K. Inkabajan, S. Pludpiem, K. Pidtana, S. Demons, B. Vesely, M. Wojnarski, J. S. Griesenbeck, and M. Spring. Co-infection with plasmodium vivax and covid-19 in thailand. Tropical Medicine and Infectious Disease, 7(8), 2022. 4. [4]. X. Chen, Q. Jiang, Z. Ma, J. Ling, W. Hu, Q. Cao, P. Mo, L. Yao, R. Yang, Gao S., X. Gui, W. Hou, Y. Xiong, J. Li, and Y. Zhang. Clinical characteristics of hospitalized patients with sars-cov-2 and hepatitis b virus co-infection. Virol Sin., 35:842–845, 2020. 5. [5]. G. Contreras, M. R. Mora, and P. J. Gómez. Estimation of the epidemiological evolution through a modelling analysis of the covid-19 outbreak. Microbiology, 3:152–158, 2020. 6. [6]. P. Das, R. K. Upadhyay, A. K. Misra, F. A. Rihan, P. Das, and D. Ghosh. Mathematical model of covid-19 with comorbidity and controlling using non-pharmaceutical interventions and vaccination. Nonlinear Dyn., 106:1213–1227, 2021. 7. [7]. M. L. Diagne, H. Rwezaura, S. Y. Tchoumi, and J. M. Tchuenche. A Mathematical Model of COVID-19 with Vaccination and Treatment. Computational and Mathematical Methods in Medicine, 2021:1250129, 2021. 8. [8]. A. Din, S. Amine, and A. Allali. A stochastically perturbed co-infection epidemic model for covid-19 and hepatitis b virus. Nonlinear Dynamics, 111(2), 2023. 9. [9]. M. F. Hama, R. R. Q. Rasul, Z. Hammouch, K. A. H. Rasul, and J. Danane. Analysis of a stochastic seis epidemic model with the standard brownian motion and lévy jump. Results in Physics, 37:105477, 2022. 10. [10]. I. M Hezam, A. Foul, and A. Alrasheedi. A dynamic optimal control model for COVID-19 and cholera co-infection in yemen. Adv. Differ. Equ., 2021(1), December 2021. 11. [11]. E. A. Iboi, C. N. Ngonghala, and A. B. Gumel. Will an imperfect vaccine curtail the covid-19 pandemic in the us? Infectious Disease Modelling, 5:510–524, 2020. 12. [12]. T. Jahnke, T. Udrescu, and M. Mikl. Numerical methods in mathematical finance, 2020. 13. [13]. C. Ji, D. Jiang, and N. Shi. The behavior of an sir epidemic model with stochastic perturbation. Stochastic Analysis and Applications, 30(5):755–773, 2012. 14. [14]. K. G. Mekonen and L. L. Obsu. Mathematical modeling and analysis for the co-infection of COVID-19 and tuberculosis. Heliyon, 8(10):e11195, October 2022. 15. [15]. E. F. Mohamed, E. K. Mohamed, G. Richard, L. Aziz, and T. Regragui. Stationary distribution and threshold dynamics of a stochastic sirs model with a general incidence. Physica A-statistical Mechanics and Its Applications, 534:120696, 2019. 16. [16]. S. E. Moore, H. L. Nyandjo-Bamen, O. Menoukeu-Pamen, J. K. K. Asamoah, and Z. Jin. Global stability dynamics and sensitivity assessment of covid-19 with timely-delayed diagnosis in ghana. Computational and Mathematical Biophysics, 10(1):87–104, 2022. 17. [17]. M. M. Ojo and D. F. E. Goufo. The impact of covid-19 on a malaria dominated region: A mathematical analysis and simulations. Alexandria Engineering Journal, 65:23–39, 2023. 18. [18].World Health Organization. COVID-19 vaccination in humanitarian settings. [https://www.who.int/publications/i/item/9789240079434](https://www.who.int/publications/i/item/9789240079434), 2023. [Online; accessed 25 October 2023]. 19. [19]. M. A. Pobbi, S. M. Naandam, and S. E. Moore. Mathematical modelling and analysis of stochastic malaria and covid-19 co-infection model. doi:10.21203/rs.3.rs-3591875/v1, 2023. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.21203/rs.3.rs-3591875/v1&link_type=DOI) 20. [20]. A. Pusparani, J. Henrina, and A. Cahyadi. Co-infection of covid-19 and recurrent malaria. The Journal of Infection in Developing Countries, 15(05):625–629, 2021. 21. [21]. M. Ray, A. Vazifdar, and S. Shivaprakash. Co-infection with malaria and coronavirus disease-2019. Journal of Global Infectious Diseases, 12(3), 2020. 22. [22]. J. M. Read, J. R. E. Bridgen, D. A. T. Cummings, A. Ho, and C. P Jewell. Novel coronavirus 2019-ncov (covid-19): early estimation of epidemiological parameters and epidemic size estimates. Philosophical Transactions of the Royal Society B, 376(1829):20200265, 2021. 23. [23]. A. Rehman, R. Singh, and P. Agarwal. Modeling, analysis and prediction of new variants of covid-19 and dengue co-infection on complex network. Chaos Solitons Fractals, 150(111008):111008, September 2021. 24. [24]. H. Rwezaura, M. L. Diagne, A. Omame, A. L. de Espindola, and J. M. Tchuenche. Mathematical modeling and optimal control of sars-cov-2 and tuberculosis co-infection: a case study of indonesia. Modeling Earth Systems and Environment, 8:5493–5520, 2022. 25. [25]. Y. Sabbar. Mathematical Analysis of Some Stochastic Infectious Disease Models with White Noises and Lévy Jumps. PhD thesis, Université Sidi Mohamed Ben Abdellah de Fès (Maroc), 2021. 26. [26]. S. M. A. Shah, Y. Nie, A. Din, and A. Alkhazzan. Dynamics of hepatitis b virus transmission with a lévy process and vaccination effects. Electron. J. Differ. Equ., 12:16–45, 2024. 27. [27]. L. Tang, Y. Zhou, L. Wang, S. Purkayastha, L. Zhang, J. … He, and P. X. K. Song. A review of multi-compartment infectious disease models. International Statistical Review, 88, 2020. 28. [28]. S. Y. Tchoumi, M. L. Diagne, H. Rwezaura, and J. M. Tchuenche. Malaria and covid-19 co-dynamics: A mathematical model and optimal control. Applied mathematical modelling, 99:294–327, 2021. 29. [29]. S. W. Teklu. Analysis of hbv and covid-19 coinfection model with intervention strategies. Comput Math Methods Med., 2023. 30. [30]. A. W. Tesfaye and T. S. Satana. Stochastic model of the transmission dynamics of covid-19 pandemic. Advances in Difference Equations, 2021:1–21, 2021. 31. [31].UNICEF. How COVID-19 caused malaria to soar – and what countries and the global community should continue to do in response. [https://data.unicef.org/data-for-action/how-covid-19-caused-malaria-to-soar](https://data.unicef.org/data-for-action/how-covid-19-caused-malaria-to-soar), 2023. [Online; accessed 25 October 2023]. 32. [32]. J. Wang and X. Tian. Global stability of a delay differential equation of hepatitis b virus infection with immune. Electron. J. Differ. Equ., 94:1–11, 2013. 33. [33]. X. Wang, C. Wang, and K. Wang. Extinction and persistence of a stochastic sica epidemic model with standard incidence rate. Advances in Difference Equations, 2021. 34. [34]. J. T. Wu, K. Leung, and G. M. Leung. Nowcasting and forecasting the potential domestic and international spread of the 2019-ncov outbreak originating in wuhan, china: a modelling study. The Lancet, 395(10225):689–697, 2020. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/s01406736(20)30260-9&link_type=DOI) 35. [35]. Z. Zhang, A. Zeb, S. Hussain., and E. Alzahrani. Dynamics of COVID-19 mathematical model with stochastic perturbation. Advances in Difference Equations, 2020(1):451, 2020. [1]: /embed/inline-graphic-1.gif [2]: /embed/inline-graphic-2.gif [3]: /embed/inline-graphic-3.gif [4]: /embed/graphic-2.gif [5]: /embed/graphic-3.gif [6]: /embed/graphic-4.gif [7]: /embed/inline-graphic-4.gif [8]: /embed/inline-graphic-5.gif [9]: /embed/inline-graphic-6.gif [10]: /embed/inline-graphic-7.gif [11]: /embed/graphic-5.gif [12]: /embed/inline-graphic-8.gif [13]: /embed/inline-graphic-9.gif [14]: /embed/graphic-6.gif [15]: /embed/graphic-7.gif [16]: /embed/graphic-8.gif [17]: /embed/graphic-9.gif [18]: /embed/graphic-10.gif [19]: /embed/inline-graphic-10.gif [20]: /embed/inline-graphic-11.gif [21]: /embed/graphic-11.gif [22]: /embed/graphic-12.gif [23]: /embed/inline-graphic-12.gif [24]: /embed/graphic-13.gif [25]: /embed/graphic-14.gif [26]: /embed/graphic-15.gif [27]: /embed/inline-graphic-13.gif [28]: /embed/graphic-16.gif [29]: /embed/graphic-17.gif [30]: /embed/inline-graphic-14.gif [31]: /embed/inline-graphic-15.gif [32]: /embed/graphic-18.gif [33]: /embed/graphic-19.gif [34]: /embed/graphic-20.gif [35]: /embed/inline-graphic-16.gif [36]: /embed/inline-graphic-17.gif [37]: /embed/inline-graphic-18.gif [38]: /embed/graphic-21.gif [39]: /embed/graphic-22.gif [40]: /embed/inline-graphic-19.gif [41]: /embed/graphic-23.gif [42]: /embed/inline-graphic-20.gif [43]: /embed/inline-graphic-21.gif [44]: /embed/inline-graphic-22.gif [45]: /embed/graphic-24.gif [46]: /embed/inline-graphic-23.gif [47]: /embed/graphic-25.gif [48]: /embed/inline-graphic-24.gif [49]: /embed/inline-graphic-25.gif [50]: /embed/inline-graphic-26.gif [51]: /embed/graphic-26.gif [52]: /embed/graphic-27.gif [53]: /embed/graphic-28.gif [54]: /embed/inline-graphic-27.gif [55]: /embed/inline-graphic-28.gif [56]: /embed/inline-graphic-29.gif [57]: /embed/graphic-29.gif [58]: /embed/graphic-30.gif [59]: /embed/inline-graphic-30.gif [60]: /embed/inline-graphic-31.gif [61]: /embed/inline-graphic-32.gif [62]: /embed/inline-graphic-33.gif [63]: /embed/inline-graphic-34.gif [64]: /embed/inline-graphic-35.gif [65]: /embed/graphic-31.gif [66]: /embed/inline-graphic-36.gif [67]: /embed/graphic-32.gif [68]: /embed/inline-graphic-37.gif [69]: /embed/inline-graphic-38.gif [70]: /embed/graphic-33.gif [71]: /embed/inline-graphic-39.gif [72]: /embed/graphic-34.gif [73]: /embed/inline-graphic-40.gif [74]: /embed/graphic-35.gif [75]: /embed/inline-graphic-41.gif [76]: /embed/inline-graphic-42.gif [77]: /embed/graphic-36.gif [78]: /embed/graphic-37.gif [79]: /embed/inline-graphic-43.gif [80]: /embed/inline-graphic-44.gif [81]: /embed/inline-graphic-45.gif [82]: /embed/inline-graphic-46.gif [83]: /embed/inline-graphic-47.gif [84]: /embed/inline-graphic-48.gif [85]: /embed/graphic-38.gif [86]: /embed/inline-graphic-49.gif [87]: /embed/inline-graphic-50.gif [88]: /embed/inline-graphic-51.gif [89]: /embed/inline-graphic-52.gif [90]: /embed/inline-graphic-53.gif [91]: /embed/inline-graphic-54.gif [92]: /embed/inline-graphic-55.gif [93]: /embed/inline-graphic-56.gif [94]: /embed/inline-graphic-57.gif [95]: /embed/inline-graphic-58.gif [96]: /embed/inline-graphic-59.gif [97]: /embed/inline-graphic-60.gif [98]: /embed/graphic-39.gif [99]: /embed/graphic-40.gif [100]: /embed/graphic-41.gif [101]: /embed/inline-graphic-61.gif [102]: /embed/graphic-42.gif [103]: /embed/graphic-43.gif [104]: /embed/graphic-44.gif [105]: /embed/graphic-45.gif [106]: /embed/graphic-46.gif [107]: /embed/graphic-47.gif [108]: /embed/inline-graphic-62.gif [109]: /embed/inline-graphic-63.gif [110]: /embed/inline-graphic-64.gif [111]: /embed/inline-graphic-65.gif [112]: /embed/graphic-48.gif [113]: /embed/graphic-49.gif [114]: /embed/inline-graphic-66.gif [115]: /embed/inline-graphic-67.gif [116]: /embed/inline-graphic-68.gif [117]: /embed/inline-graphic-69.gif [118]: /embed/graphic-50.gif [119]: /embed/inline-graphic-70.gif [120]: /embed/inline-graphic-71.gif [121]: /embed/inline-graphic-72.gif [122]: /embed/graphic-51.gif [123]: /embed/graphic-52.gif [124]: /embed/graphic-53.gif [125]: /embed/graphic-54.gif [126]: /embed/inline-graphic-73.gif [127]: /embed/graphic-55.gif [128]: /embed/graphic-56.gif [129]: /embed/graphic-57.gif [130]: /embed/inline-graphic-74.gif [131]: /embed/graphic-58.gif [132]: /embed/graphic-59.gif [133]: /embed/inline-graphic-75.gif [134]: /embed/inline-graphic-76.gif [135]: /embed/inline-graphic-77.gif [136]: /embed/inline-graphic-78.gif [137]: /embed/inline-graphic-79.gif