key: cord-0812674-j8fw9qi4
authors: Ma, Zhihui; Han, Shuyan; Li, Shenghua
title: A stochastic eco-epidemiological system with patchy structure and transport-related infection
date: 2021-11-13
journal: J Math Biol
DOI: 10.1007/s00285-021-01688-x
sha: 311798b641e2165f72ddceae83391c40da7cd4ce
doc_id: 812674
cord_uid: j8fw9qi4

In this paper, a stochastic eco-epidemiological system with patchy structure and transport-related infection is proposed and the stochastic dynamical behaviors are investigated. Firstly, by constructing suitable Lyapunov functions, it is revealed that there is a unique globally positive solution starting from the positive initial value. Secondly, it is proved that the presented system is stochastically ultimately bounded and the average in time of the second moment of solution is bounded. Thirdly, we prove that the large enough stochastic perturbations may lead the predator population and the diseases in the predator to be extinct while it is persistent in the deterministic system. Finally, some numerical simulations are given to test our theoretical results.

Eco-epidemiology is an important branch in mathematical biology which both considers the ecological and epidemiological factors. The influence of diseases in the ecological system is an important issue from the mathematical and ecological perspectives. Many researchers presented eco-epidemiological systems and tended to explore ecological systems being subject to epidemiological aspects (Anderson and May 1986; Venturino 1993; Bhattacharyya and Mukhopadhyay 2010; Venturino 2002; Greenhalgh et al. 2020; Haque and Venturino 2007; Hethcote 2000; Mainul 2010; Gulland 1995; Hsieh and Hsiao 2008; Fenton and Rands 2006) . Anderson and May (Anderson and May 1986) firstly proposed an eco-epidemiological system which incorporates diseases into the interacting populations. After that, many researchers proposed and investigated eco-epidemiological systems for the spread of diseases among interacting species and revealed the spread mechanisms of diseases based on their studies (Bhattacharyya and Mukhopadhyay 2010; Venturino 2002; Greenhalgh et al. 2020; Haque and Venturino 2007) . Some works studied the spread of diseases in predator/prey population in Lotka-Volterra type predation systems (Bhattacharyya and Mukhopadhyay 2010; Venturino 2002; Mainul 2010; Saha and Samanta 2020; Saha et al. 2018; Mondal and Samanta 2021) . Hsieh and Hsiao (2008) and Fenton and Rands (2006) proposed a predator-prey system with disease in both prey and predator populations and considered the system dynamical behaviors. The results showed that the dynamical behaviors of predation systems with diseases are more complex than those without diseases. In fact, a certain disease in predator populations is a very common phenomenon and the biologically relevant examples in ecosystems are found by Gulland (1995) . For examples, rabies and Sarcoptes spp. in foxes (Vulpis vulpis) and coyotes (Canis latrans), where rabbits are prey populations. The Phocine Distemper Virus in both the common seal (Phoca vitulina) and the striped dolphin (Stenella coeruleoalba), Oto-strongylus circumlitis, and Uncinaria lucasi, which affect ringed seals (Phoca hispida) and northern fur seals (Callorhinus ursinus), respectively. In this case, fishes are prey populations. The other example is that common seals (Phoca vitulina) and striped dolphins (Stenella coeruleoaba) are infected with canine distemper virus (PDV). In 1988, PDV infectious disease was found in seals. In 1988, PDV was first identified as the cause of death of 18,000 harbor seals (P. vitulina) and 300 gray seals (Halichoerus grypus) on the Nordic coast. In this case, fish is considered as prey and the disease will not cross the species barrier. The Avian Pox, Newcastle Disease, Influenza, Pasteurella multocida, Apergillus fumigatus and Leukocytozoon simondi are infectious agents and all affect a wide range of avian species, in which prey individuals are mainly insects (Venturino 2002) .

On the one hand, the environment habitats are heterogeneous and mainly formed by a set of discrete patches. This induces populations survive in patchy habitats and are associated with their migrations. The migration rate explicates the evolutionarily stable patterns of the populations, and moreover plays an important role in persistence of the interacting populations. One of the popular way is the classical idea of evolutionary ecology where the less fit phenotypes, which are distinguished based on their migration and/or dispersal properties, are replaced by fitter ones as a result of evolution. Based on this issue, the ideal free distribution (IFD) as one of the popular concept are presented . Hence, the population systems in patchy habitats with migration (dispersion) are more realistic. However, it is inevitable that some infectious diseases may spread form one region to another as populations migrate or transport between regions. The migration or transportation among regions was one of the main factors which affected the outbreak of infectious diseases. The migration or transportation will change the disease's dynamics and induce infection diseases to be endemic even if infectious diseases will go to extinction in each region without transport-related infection. This will induce more serious consequences in population survival and evolution. Such as the influenza and COVID-19 will transmit form one country/city to other countries/cities. Therefore, the patchy model systems with the infectious diseases during transport are more realistic and essential, and have been studied by many researchers (Hamilton and May 1977; Abrams et al. 2007; Liu and Takeuchi 2006; Allen 2011; Cui et al. 2006; Allen and Bokil 2012; Hu et al. 2019; Rudnicki 2003; Ji and Jiang 2011; Bao and Shao 2015) . For example, Liu and Takeuchi (2006) proposed an S I QS system to study the effect of transport-related infection and entry screening. The results revealed that the entry screening may be helpful for disease eradication since it can always have the possibility to eradicate the disease led by transport-related infection, and furthermore have the possibility to eradicate disease even when the disease is endemic in both isolated cities. Cui et al. (2006) presented an S I S epidemic model system to understand the effect of transport-related infection on the disease's spread and the corresponding mechanism. The results showed that the transport-related infection intensifies diseases spread if infectious diseases break out to cause an endemic situation in each region. However, almost all published researches for eco-epidemiological systems do not consider the possibility for individuals which may be infected during their travels. This factor plays a very important role in determining the spread and development of infectious diseases, and it is very essential to strengthen restrictions of passengers once infectious diseases appear. Motivated by this, this work firstly proposed an eco-epidemiological system incorporating the transport-related infection.

On the other hand, the dynamical behaviors of ecological populations are always affected by environmental noises or disturbances (Cui et al. 2006; Allen and Bokil 2012; Hu et al. 2019; Rudnicki 2003; Ji and Jiang 2011; Bao and Shao 2015; Mao 2007; Rebelo and Soresina 2020; Xin and Liu 2019; Arifah and Mao 2004; Li and Mao 2009; Das and Samanta 2018) . Deterministic population models in ecology do not usually incorporate environmental fluctuation. They are often justified by the implicit assumption that stochastic deviations are small enough to be ignored in large populations. Deterministic models will be proved ecologically useful only if the dynamical patterns are still in evidence when stochastic effects are introduced. Uncertain changes of populations are usually considered as an effect of environmental stochasticity. In fact, the physical and biological environment of populations is not totally predictable, such as the population reproduction depends on temperature, humidity, parasites, pathogens, environmental pollution and so all. Hence, the population's development and growth should be considered as a stochastic process rather than a deterministic one. Rudnicki (2003) provided a detailed analysis for stability of a stochastic Lotka-Volterra system. Ji and Jiang (2011) and studied a predator-prey system with modified Leslie-Gower and Holling type II schemes with stochastic perturbations. Bao and Shao (2015) further extended this system into a stochastic predation system with all coefficients exhibiting random fluctuations. The published researches have revealed that the environmental fluctuations will induce the birth rate, death rate, transmission coefficient and other parameters of dynamical systems to exhibit randomness with a greater or lesser extent. Moreover, the random fluctuations could apply to control some infectious diseases. Motivated by this, this paper proposed an eco-epidemiological system incorporating the patchy structure, transport-related inflection and stochastic perturbation.

The paper is organized as follows. In Sect. 2, an stochastic eco-epidemiological system with patchy structure and transport-related infection is proposed. The unique global positivity of solutions of the presented system is given in Sect. 3. Section 4 shows that the presented model system is stochastically ultimately bounded and the average in time of the second moment of solution is bounded. In Sect. 5, we will deduce the condition of extinction of the total predator populations and disease in the predators. Finally, the theoretical results are discussed and some numerical simulations are given to test our theoretical results.

Based on the above realistic issues, suppose that diseases only spread in predator population, and predator population is divided into two classes: the susceptible predator and the infected predator. Let X i , S i and I i represent the number of prey, susceptible predator, infected predator in patch i (i = 1, 2), respectively.

The basic assumptions for system formation are as follows:

(i) In the absence of predators, the prey population grows logistically with r i (1− X i K i ) in patch i (i = 1, 2), respectively. (ii) Holling I functional response is adopted for susceptible predators, and they will lose the predation ability after being infected by diseases. (iii) The standard incidence rate is adopted in the spread of diseases among predator populations, and β i S i I i S i +I i is the total number of new cases per unit time in patch i (i = 1, 2), respectively. However, this paper supposes that the mass-action incidence rate for the infection during transports due to the sufficient and frequent contacts between susceptible and infected predators, that is, γ i α 2 i j S i I i is the total number of new cases per unit time during transport from patch i to patch j, (i, j = 1, 2, i = j). (iv) The birth rate and death rate of predator individuals who are traveling are omitted since the birth and death processes are more fast than the predators' one-way transport between two patches. (v) It is assumed that stochastic perturbations are white noises which are directly proportional to X i (t), S i (t) and I i (t) (i = 1, 2), and influenced on theẊ i (t), S i (t) andİ i (t) in the proposed system (i = 1, 2) respectively. Hence, the random perturbations incorporate into incidence functions in patch 1 and patch 2 are as follows:

Based on the above assumptions and diagram, the following compartment diagram is presented:

Based on the above assumptions and diagram, the following stochastic ecoepidemiological system with patchy structure and transport-related infection is proposed:

(2.1)

The ecological and epidemiological meanings of all parameters of system (2) are given in Table 1 

are all positive and denote the intensities of the white noises (i = 1, 2, 3, 4, 5, 6, 7, 8, j = 1, 2).

From the biological point of view, the term α i j S i means that the susceptible predators leave the patch i and γ i α 2 i j S i I i denotes the susceptible predators become infected The effective contact rate at patch i σ i Natural death rate of the predator at patch i δ i

The individuals' rate of recovery due to natural causes or treatment ω i Disease-induced death rate of the predator at patch i α i j The travel rate from patch i to patch j γ i

The effective contact rate in the transportation from patch i to patch j during their travels from patch i to patch j. Therefore, it is assumed that the term α i j S i − γ i α 2 i j S i I i should be nonnegative and the term 1 − γ i α i j I i is positive throughout this paper.

Defining

, and the initial conditions of system (2.1) are as follow:

, P) be a complete probability space with a filtration {F t } t≥0 satisfying the usual condition (i.e., it is increasing and right continuous while F 0 contains all P-null sets), B i (t)(i = 1, 2, 3, 4, 5, 6, 7, 8) are defined on this complete probability space. Meanwhile, defining R d

Generally, we consider the following d-dimensional stochastic differential equation: B(t) denotes an m-dimensional standard Brownian motion defined on the complete probability space ( , F, {F t } t≥0 , P). The differential operator L of system (2.1) is defined by

R + ) denote the family of all nonnegative functions defined on R d × [t 0 , +∞) and continuously twice differentiable in x and once in t. 

In this section, we will prove that system (2.1) has a unique global positive solution by using Lyapunov analysis method.

, and the solution φ(t) will remain in R 6 + with probability one.

Proof Since the coefficients of system (2.1) are locally Lipschitz continuous, for any given initial value (2.2), there is a unique local solution in [0, τ e ), where τ e denotes the explosion time.

Now, we will prove that the solution is global, i.e. τ e = +∞ a.s. To this end, let k 0 ≥ 1 be sufficiently large such that

For each integer k ≥ k 0 define the stopping time as follow

we denote as the empty set and set inf = ∞ throughout this paper. Clearly, τ k is increasing as k → ∞. Set τ ∞ = lim k→∞ τ k , hence τ ∞ ≤ τ e a.s. If we can show that τ ∞ = ∞ a.s., then τ e = ∞ and φ(t) ∈ R 6 + a.s. for all t ≥ 0. If this assertion is not true, then there exist T > 0 and ∈ (0, 1) such that P{τ ∞ ≤ T } > . For any k ≥ k 0 , let k = {ω ∈ : τ k (ω) ≤ T }. That is, P( k ) > for any k ≥ k 0 . Define the C 6 -function V:R 6

For any 0 ≤ t < (τ k ∧ T ), using I tô s f ormula, we have

Using the inequality u ≤ 2u − 1 − lnu for any u ∈ R + , we have

Furthermore, for any k > k 0 , integrating both sides of (3.1) from 0 to τ k ∧ T , and then take the expectation, we get

By the Gronwall inequality , it is obtained that

(3.2)

Noting that for every ω ∈ k , there is at least one of X 1 (τ k , ω), S 1 (τ k , ω), I 1 (τ k , ω),

It follows from (3.2) and (3.3), we have

where I k (ω) denotes the indicator function of k . Let k → +∞, then +∞ > (V (φ(0)) + GT )e max{λ 1 ,λ 2 }T = +∞, which leads to a contradiction. Hence, it must have τ ∞ = +∞ a.s. as required.

Next, we claim that φ(t) is positive almost surely on [0, ∞). On t ∈ [0, ∞), for the system (2.1), we have

where i, j = 1, 2, i = j, k = 2, 6, m = 4, 8.

where i, j = 1, 2, i = j, k = 3, 7, m = 4, 8. Apparently, for the initial condition (2.2), we have

Therefore, (2.1) with (2.2) has the unique global solution φ(t) on t ∈ [0, ∞) and the solution will remain in R 6 + with probability one. This completes the proof.

In this section, we examine the ultimate boundedness in mean, stochastically ultimate boundedness, and the average in time of the second moment of solution.

where ξ = min{σ 1 , σ 2 }, Further, system (2.1) is ultimately bounded in mean.

Proof By Theorem 3.1, the solution φ(t) will remain in R 6 + for all t > 0 with probability one. Define the C 6 -function V:R 6 + → R by V 1 (t, φ(t)) = e ξ t (X 1 + X 2 + S 1 + S 2 + I 1 + I 2 ).

Applying I tô s f ormula, we have

Integrating both sides of (4.1) from 0 to t and using the above equality, we obtain

(4.2)

Now, multiplying both sides of (4.2) by e −ξ t , and taking the expectation, we have

Noting from the equality |φ(t)| ≤ X 1 (t) + X 2 (t) + S 1 (t) + S 2 (t) + I 1 (t) + I 2 (t), for any φ(t) ∈ R 6 , hence lim sup t→∞ E|φ(t)| ≤ K 1 (r 1 + ξ) 2 4r 1 ξ + K 2 (r 2 + ξ) 2 4r 2 ξ .

This completes the proof.

By Chebyshev's inequality and Theorem 4.1, we have the following conclusion.

Theorem 4.3 Assume that there are positive numbers m i , n i , l i , θ i , such that R i :=

i 2 − θ i max{γ 1 α 2 12 , γ 2 α 2 21 }B > 0, i = 1, 2, where m 1 m 2 = 1, n 1 n 2 = 1, l 1 l 2 = 1, θ 1 θ 2 = 1. Then, for any given initial values (2.2), the solution φ(t) of system (2.1) satisfies lim sup

Proof Firstly, Integrating both sides of d(X 1 (t)) ≤ r 1 X 1 (1 − X 1 K 1 ) + a 1 X 1 d B 1 (t) from 0 to t and taking the expectation, we have

This implies lim sup

By the same calculation, we can obtain that lim sup

By theorem 4.1 and theorem 4.2, we have S 1 (t)+S 2 (t) < B = K 1 (r 1 +ξ) 2 4r 1 ξ + K 2 (r 2 +ξ) 2 4r 2 ξ for any sufficiently small > 0. Define the C 6 -function V:R 6 + → R by

Applying I tô s f ormula and the inequality 2ab ≤ a 2 + b 2 , we have

Integrating both sides of the above inequality from 0 to t and taking the expectation, we get

It is straightforward to see that lim sup

This completes the proof.

In this section, we consider asymptotic pathwise estimation the solution of system (2.1) and the extinction of the total predator and the extinction of ecological epidemics with probability one. Because the predator can travel freely between two patches, the extinction considered in the section must occur simultaneously in two patches. Firstly, we give the proof of the extinction of the total predator with probability one.

where H 1 is defined in (5.1) . Moreover, if H 1 < 0 holds, then the total predator of system (2.1) will tend to zero exponentially with probability one.

V 2 (φ(t)) = ln(S 1 + I 1 + S 2 + I 2 ).

Applying I tô s f ormula, we have

(5.1)

Hence,

Integrating both sides of (5.2) from 0 to t and then dividing t, we get ln(S 1 (t) + I 1 (t) + S 2 (t) + I 2 (t)) t ≤ ln(S 1 (0) + I 1 (0) + S 2 (0) + I 2 (0)) t

By the Strong law of large numbers for martingales, we have lim t→∞ B i (t) t = 0, i = 2, 3, 6, 7.

Thus,

It is straightforward to see that, when H 1 < 0 is satisfied, S 1 + I 1 + S 2 + I 2 will go to zero exponentially with probability one, i.e the total predator of two patches will go to extinction simultaneously exponentially with probability one.

This completes the proof.

Remark 1 On the one hand, it is easy to see that H 1 = − min{σ 1 − p 1 K 1 , σ 1 + ω 1 , σ 2 − p 2 K 2 , σ 2 + ω 2 } − 1 8 min{b 2 1 , c 2 1 , b 2 2 , c 2 2 } is decreasing with b 2 1 , c 2 1 , b 2 2 , c 2 2 increasing. Hence, the total number of predators in two patches will die out exponentially as long as b 2 1 , c 2 1 , b 2 2 and c 2 2 are large enough such that H 1 < 0. This means that the white noise may lead predator population to be extinction while they are persistent in the deterministic system. On the other hand, we can obtain that lim sup

i 2 a.s. (i = 1, 2). That is to say, if the white noise is large enough, the prey population will be extinct in statistic system while they are persistent in the deterministic system. Therefore, the larger white noise will have negative effects on the population coexistence.

Next, we will give the proof of the extinction of the disease with probability one by using the same method as Theorem 5.1.

where H 2 is defined in (5.6). Moreover, if H 2 < 0 holds, then the disease of system (2.1) will go to extinction exponentially with probability one.

V 3 (φ(t)) = ln(I 1 + I 2 ).

Applying I tô s f ormula, we have

Next, by theorem 4.1 and theorem 4.2, we have S 1 (t) + S 2 (t) < B = K 1 (r 1 +ξ) 2 4r 1 ξ + K 2 (r 2 +ξ) 2 4r 2 ξ for any sufficiently small > 0. Hence,

where D = −min{σ 1 +ω 1 +δ 1 , σ 2 +ω 2 +δ 2 }− 1 4 min{c 2 1 , c 2 2 }+max{γ 2 α 2 21 , γ 1 α 2 12 }B.

Integrating both sides of (5.3) from 0 to t, we get ln(I 1 (t) + I 2 (t)) ≤ ln(I 1 (0) + I 2 (0)) + Dt + t 0 1 I 1 + I 2 ( β 1 S 1 I 1

is a continuous local martingale whose quadratic variation is as follows

By virtue of the exponential martingale inequality , it is obtained that

where 0 < c < 1,k is a random integer. Applying Borel-Cantelli lemma leads to that for almost all ω ∈ , there exists a random integer k 0 > 0, such that for k > k 0 , we have

that is to say, for any 0 ≤ t ≤ k, we get . 2 The influence of the white noises when only the prey is present in the system (2.1), where a: a 1 = 0.2, a 2 = 0.3, b: a 1 = 2.2, a 2 = 2

In this section, some numerical simulations for system (2.1) are conducted to illustrate our analytical results. The system parameter values are shown in Table 2 and we take initial values φ(0) = (2.5, 1.5, 0.5, 2.5, 1.5, 0.5). Some of these parameter values have been taken from literatures . In order to investigate the influence of stochastic perturbations on system (2.1), the numerical simulations on the stochastic system and the corresponding deterministic system under the same parameters are compared.

From Fig. 2a shows r 1 − a 2 1 2 = 1.98 > 0 and r 2 − 1 a 2 2 = 1.455 > 0, and reveals that the relatively small white noise make the solutions of system (2.1) fluctuating in a large scale. Figure 2b shows r 1 − a 2 1 2 = −0.42 < 0 and r 2 − 1 a 2 2 = −0.5 < 0, and means that large enough white noise will lead to the extinction of prey population while it is persistent in the deterministic system .

Choosing σ 1 = 0.56 and σ 2 = 0.51 for Fig. 3 , and the other parameters are the same as those in Table 2 . Figure 3a shows H 1 = 0.9871875 > 0, and reveals that the relatively small white noise make the solutions of system (2.1) fluctuating in a large scale, while Fig. 3b illustrates that the large enough white noise will lead to extinction of the predator population while it is persistent in the deterministic system because of H 1 = −0.5853125 < 0, which is consistent with the theoretical results of Theorem 5.1. For the effect of the white noise to disease spread, selecting σ 1 = 0.09, σ 2 = 0.1, ω 1 = 0.02 and ω 2 = 0.04, the other parameters are the same as those in Table 2 . Figure 4a shows that the relatively small white noise make the solutions of system (2.1) fluctuating in a large scale. Figure 4b illustrates that large enough white noise leads diseases in predator population to be extinct while it is persistent in the deterministic system, which is consistent with the theoretical results of Theorem 5.2.

In this paper, a predator-prey system with diseases in predator population only incorporating stochastic perturbations, patchy structure and transport-related inflection has been proposed and investigated. By constructing different Lyapunov functions, we proved that the existence and uniqueness of the global positive solution, stochastically ultimate boundedness of the solution and extinction of diseases of system (2.1). Compared with some published researches, ecological epidemics and white noise are incorporated in our model since infectious diseases and environmental noises are widely exist in natural world, and paly an important effect on population's evolution. Furthermore, the existing evidences show that infectious diseases will be endemic while environmental noises are ignored. However, our study shows that environmental noises are very important factors to control infectious diseases, and infectious diseases will be extinct under a certain strength of environmental noises, which is disagreement with some published works.

Moreover, the outbreak of COVID-19 shows that the transport-related inflection plays an importantly positive role in the spread of infectious diseases, and the transport-related inflection make diseases to be transmitted in a more large scale. Our investigation will also provide a strong evidence for the positive effect of the transportrelated inflection on the transmission of infectious diseases. Based on this issue, our results are agreement with previous works in this filed. Some previous works sug- Fig. 4 The influence of the white noises when the disease exist in the system (2.1), where a: c 1 = 0.03, c 2 = 0.04, η 1 = 0.02, η 2 = 0.014, b: c 1 = 3.15, c 2 = 3.3, η 1 = 4, η 2 = 3.75 gested that the population's transport among different regions should be completely prohibited while infectious diseases outbroke. However, our research suggests that the transport should be prohibited under a certain scale while infectious diseases spread, and the corresponding threshold value is obtained. This controlling measure for infectious diseases could reduce economic and other losses of countries and cities as far as possible. Based on this issue, our conclusions are disagreement with some previous researches.

Furthermore, compared with stochastic models which has presented by previous works (Hu et al. 2019; Xin and Liu 2019) , few of them incorporated the patchy structure into corresponding models. In fact, habitat on which populations survive are fragmented as human activities. The environmental fragmentation is negative to the population's permanence and coexistence, and will cause the loss of species biodiversity. In this paper, the patchy structure is incorporated into a stochastic ecological epidemic model, and two patches are connected by population's migrations. Hence, our model is more realistic compared with some previous researches. Meanwhile, the coexistence threshold values of interacting populations and the extinct threshold values of infectious diseases are found respectively. Our results show that the population's permanence and the disease extinction will be achieved for populations or species which survive in more complex environment. Based on this issue, our study is more useful for the biodiversity conservation and the diseases control.

Although our investigation considered many realistic issues and obtained some reasonable results, There are existing some interesting issues to be investigated furtherly. On the one hand, it is necessary to incorporate the time-delay which measures the length of the time to complete one-way transport between two patches. On the other hand, one may propose some more realistic model systems, such as considering the nonlinear incidence rate and the generalized functional response. These will be considered in future work.

The role of behavioral dynamics in determining the patch distributions of interacting species

An introduction to stochastic processes with applications to biology

The invasion, persistence, and spread of iufectious diseases within animal and plant communites

Stochastic delay Lotka-Volterra model

Permanence and extinction of regime-switching predator-prey models

On an eco-epidemiological model with prey harvesting and predator switching: Local and global perspectives

Spreading disease with transport-related infection

Stochastic prey-predator model with additional food for predator

The impact of parasite manipulation and predator foraging behaviour on predator-prey communities

Eco-epidemiological model with fatal disease in the prey

The impact of infectious diseases on wild animal populations-a review, Ecology of Infectious Diseases in Natural Populations

Dispersal in stable habitats

An ecoepidemiological model with disease in predator: the ratio-dependent case

The mathematics of infectious diseases

Predator-prey model with disease infection in both populations

A comparison of deterministic and stochastic predator-prey models with disease in the predator

Dynamics of a stochastic density dependent predator-prey system with Beddington-DeAngelis functional response

Population dynamical behavior of non-autonomous Lotka-Volterra competitive system with random perturbation

Spread of disease with transport-related infection and entry screening

Dynamics of a stochastic regime-switching predator-prey model with harvesting and distributed delays

Persistence and extinction of a modified Leslie-Gower Holling-type II stochastic predator-prey model with impulsive toxicant input in polluted environments

A predator-prey model with disease in the predator species only

Chichester Mondal S, Samanta G (2021) Time-delayed predator-prey interaction with the benefit of antipredation response in presence of refuge

Coexistence in seasonally varying predator-prey systems with Allee effect

Long-time behaviour of a stochastic prey-predator model

Influence of dispersal and strong Allee effect on a two-patch predator-prey model

Analysis of a predator-prey model with herd behavior and disease in prey incorporating prey refuge

A Prey-predator system with disease in prey and cooperative hunting strategy in predator

A Michaelis-Menten Predator-prey model with strong Allee effect and disease in prey incorporating prey refuge

The influence of diseases on Lotka-Volterra systems

Epidemics in predator-prey models: disease in the predators

Analysis of stochastic Nicholson-type delay system with patch structure

Publisher's Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations

Since β 1 S 1 I 1 (S 1 + I 1 )(I 1 + I 2 ) − 1 − c 2 η 2 1 S 2 1 I 2 1 (S 1 + I 1 ) 2 (I 1 + I 2 ) 2 ≤ β 2 1 2(1 − c)η 2 1 β 2 S 2 I 2 (S 2 + I 2 )(It is obtained thatTherefore, for k − 1 < t < k, one can see thatLetting k → +∞, i.e., t → +∞ and in view of the Strong law of large numbers to the Brownian motion, we derive lim t→∞Letting c → 0, we get lim supIt is straightforward to see that, when H 2 < 0 holds, I 1 + I 2 will go to zero exponentially with probability one, i.e. the disease of two patches will go to extinction simultaneously exponentially with probability one.This completes the proof.Remark 2 It is easy to see that H 2 = −min{σ 1 +ω 1 +δ 1 , σ 2 +ω 2 +δ 2 }− 1 4 min{c 2 1 , c 2 2 }+ min{γ 2 α 2 21 , γ 1 α 2 12 }B + is decreasing in c 2 1 , c 2 2 , η 2 1 , η 2 2 . Hence, the disease in two patches will exponentially die out simultaneously as long as c 2 1 , c 2 2 , η 2 1 and η 2 2 are large enough such that H 2 < 0. This means that the white noise may lead diseases to be extinct while these diseases will be endemic in the deterministic system. Therefore, the white noise is beneficial to the infectious disease control.Funding This work was supported by Natural Science Foundation of Gansu Province (No. 21JR7RA165, 20JR5RA238, 21JR7RA535) and the Fundamental Research Funds for the Central Universities (lzujbky-2021-54).