Stochastic Lotka–Volterra mutualism model with jumps

The existence and uniqueness of the global positive solution are proved for the system of stochastic differential equations describing a two-species Lotka–Volterra mutualism model disturbed by white noise, centered and noncentered Poisson noises. For the considered system, sufﬁcient conditions of stochastic ultimate boundedness, stochastic permanence, non-persistence and strong persistence in the mean are obtained


Introduction
In nature we can find many examples where the interaction of two or more species is to the advantage of all.These population systems are described by mutualism models.The simplest two-spices Lotka-Volterra mutualism model and its properties are presented in J.D. Murray [1].A deterministic nonautonomous two-species Lotka-Volterra mutualism model is described by the system dx i (t) = x i (t) r i (t) − a ii (t)x i (t) + a ij (t)x j (t) dt, i, j = 1, 2, i = j, where x i (t), i = 1, 2, denote population densities of each species at time t, r i (t) > 0, i = 1, 2, denote the intrinsic growth rates of species x i (t), i = 1, 2. The carrying capacities of species x i (t) at time t are r i (t)/a ii (t) > 0, i = 1, 2, and coefficients a ij (t) > 0, i, j = 1, 2, i = j , describe the influence of the j -th population upon the i-th population at time t.
In the real world population systems are often subject to environmental noise.Therefore, it is natural to describe such systems by the systems of stochastic differential equations.In the paper by Peiyan Xia et al. [2], the authors consider the stochastic nonautonomous two-species Lotka-Volterra mutualism model of the form dx i (t) = x i (t) r i (t) − a ii (t)x i (t) + a ij (t)x j (t) dt + σ i (t)dw i (t) , (1) i, j = 1, 2, i = j , where r i (t), a ij (t), σ i (t), i = 1, 2, are all positive, continuous and bounded functions on [0, +∞), and w 1 (t), w 2 (t) are mutually independent Wiener processes.The authors show that the stochastic system (1) has a unique global (no explosion in a finite time) solution for any positive initial value and that the p-th moment of the solution is bounded.The sufficient conditions for stochastic permanence, persistence in the mean, nonpersistence and global attractivity of the system (1) are obtained.
In the paper by L. Shaikhet and A. Korobeinikov [3], the authors studied the asymptotic properties of Lotka-Volterra competition and mutualism models driven by autonomous system of stochastic differential equations Here, a 1 and a 2 are per capita rates of growth of populations x(t) and y(t), respectively, b 11 , b 12 , b 21 and b 22 reflect the intraspecific competition (b 11 > 0 and b 22 > 0) and interspecies interaction (b 12 and b 21 ).For competing or symbiotic species, a 1 > 0, a 2 > 0. For competing populations b 12 > 0, b 21 > 0, and for symbiotic populations b 12 < 0, b 21 < 0. The authors showed that solutions to the considered system with positive initial conditions converge to a certain compact region in the model phase space and oscillate around this region thereafter.So, the solutions of the considered stochastic system are bounded and the system is persistent.The necessary condition of the species extinction is obtained, sufficient conditions of the extinction of competing species, for which b 12 > 0, b 21 > 0 are derived, but sufficient conditions of the extinction of symbiotic species, for which b 12 < 0, b 21 < 0, were not obtained in this paper.
If we want to take into account abrupt environmental perturbations, such as epidemics, fires, earthquakes, etc. in the considered models, we must introduce Poisson noises into the population models for describing such discontinuous systems.In the paper by Y. Gao and X. Zhang [4], the authors considered a two-dimensional autonomous stochastic Lotka-Volterra mutualistic parasite-host system with "white" noise and "small" jumps, corresponding to the centered Poisson measure where x i (t − ), i = 1, 2, are the left limits of x i (t), i = 1, 2, w i (t), i = 1, 2, are mutually independent standard one-dimensional Wiener processes, Ñ(t, A) = N(t, A) − tλ(A), A is a Borel set in R, λ(Z) < +∞, N(t, A) is the Poisson measure, which is independent of w i (t), i = 1, 2, r i > 0, a ij > 0, σ i > 0, i, j = 1, 2. The sufficient conditions for extinction and persistent in the mean of species x i (t), i = 1, 2, are obtained.Then, the authors established the sufficient criteria for stability in distribution of the considered parasite-host system.
In the paper by J. Bao et al. [5], the authors studied the stochastic competitive multi-species Lotka-Volterra model under the action of the one-dimensional "white" noise and jumps, generated by a centered Poisson measure.The authors proved that the model admits a unique global positive solution, which has a uniformly finite p-th moment with p > 0. Stochastic ultimate boundedness, existence of invariant measure and long-term behaviors of solutions are discussed.
The paper by Q. Liu et al. [6] is devoted to the study of two-species mutualism model driven by the system of the autonomous stochastic differential equations where w ij (t), i, j = 1, 2, are mutually independent standard one-dimensional Wiener processes, independent of N(t, A).It is shown that the positive solution of the considered system is stochastically ultimate bounded.Then authors establish sufficient and necessary conditions for the stochastic permanence and extinction of the system.The impact of centered and noncentered Poisson noises to the stochastic nonautonomous mutualism model is studied in the paper by Olg.Borysenko and O. Borysenko [7].The authors considered the system of nonautonomous stochastic differential equations where w i (t), i = 1, 2, are mutually independent standard one-dimensional Wiener processes, , are a finite measures on the Borel sets A in R. The existence and uniqueness of the global positive solution to the considered system is proved.The authors obtain sufficient conditions of stochastic ultimate boundedness, stochastic permanence, nonpersistence in the mean, strong persistence in the mean and extinction of the solution to the considered system.
In this paper, we consider the nonautonomous stochastic mutualism model with jumps generated by centered and noncentered Poisson measures.So, the novelty of considered model is following: we investigate the nonautonomous stochastic Lotka-Volterra model and we take into account not only "small" jumps, corresponding to the centered Poisson measure but also the "large" jumps, corresponding to the noncentered Poisson measure.This model is driven by the system of nonautonomous stochastic differential equations where , are finite measures on the Borel sets A in R.
As far as we know, there are no papers devoted to the dynamical properties of the stochastic mutualism model (2), even in the case of a centered Poisson noise.
In the following we will use the notations For the bounded, continuous function f (t), t ∈ [0, +∞), let us denote We prove that system (2) has a unique, positive, global solution for any positive initial value and this solution is stochastically ultimate bounded.The sufficient conditions for stochastic permanence, nonpersistence and strong persistence in the mean of the system are derived.
The rest of this paper is organized as follows.In Section 2, we prove the existence of the unique global positive solution to the system (2).In Section 3, we prove the stochastic ultimate boundedness of the solution to the system (2).In Section 4, we obtain conditions under which the solution to the system (2) is stochastically permanent and strong persistence in the mean.In Section 5 the sufficient conditions for nonpersistence of the system (2) are obtained.

Existence of a global solution
Let ( , F, P) be a probability space, w i (t), i = 1, 2, t ≥ 0, are mutually independent standard one-dimensional Wiener processes on ( , F, P), and ν i (t, A), i = 1, 2, are mutually independent Poisson measures defined on ( , are finite measures on the Borel sets in R. On the probability space ( , F, P) we consider an increasing, right continuous family of complete sub-σ -algebras {F t } t≥0 , where We need the following assumption.In what follows we will assume that Assumption 1 holds.Theorem 1.There exists a unique global solution X(t) to the system (2) for any initial value X(0) = X 0 > 0, and Proof.The idea of proof is taken from [2].The coefficients of the system (2) are locally Lipschitz continuous.Therefore, for any initial value X 0 there exists a unique local solution X(t) = (x 1 (t), x 2 (t)) on [0, τ e ), where sup t<τ e |X(t)| = +∞ (cf.Theorem 6, p. 246, [8]).To show this solution is global, we need to show that τ e = +∞ a.s.Let n 0 ∈ N be sufficiently large for x i0 ∈ [1/n 0 , n 0 ], i = 1, 2. For any n ≥ n 0 we define the stopping time Clearly, τ n is increasing as n → +∞.Set τ ∞ = lim n→∞ τ n , whence τ ∞ ≤ τ e a.s.If we prove that τ ∞ = ∞ a.s., then τ e = ∞ a.s. and X(t) ∈ R 2 + a.s.for all t ∈ [0, +∞).So, we need to show that τ ∞ = ∞ a.s.If this statement is false, there are constants T > 0 and ε ∈ (0, 1), such that P{τ ∞ < T } > ε.Hence, there is For the nonnegative function by the Itô formula, system (2), and definition of the stochastic integral with respect to the noncentered Poisson measure ν 2 (dt, dz) we derive where Using the inequality x 1 x 2 ≤ εx and under conditions of the theorem, there is a constant K > 0, such that L(x 1 , x 2 , t) ≤ K. Hence, from (5) we have Whence taking expectations, we have Set n = {ω ∈ : Then by (4), P( n ) = P{τ n ≤ t} > ε, ∀n ≥ n 1 .Note that for every ω ∈ n there is some i such that x i (τ n , ω) equals either n or 1/n.So, It then follows from (6) that where 1 n is the indicator function of n .Letting n → ∞ leads to the contradiction ∞ > V (X 0 ) + KT = ∞.This completes the proof of the theorem.

Stochastically ultimate boundedness
Definition 1 ([9]).The solution X(t) to the system ( 2) is said to be stochastically ultimately bounded, if for any ε ∈ (0, 1), there is a positive constant χ = χ(ε) > 0, such that for any initial value X 0 ∈ R 2 + , the solution to the system (2) has the property that lim sup Theorem 2. The solution X(t) to the system (2) is stochastically ultimately bounded for any initial value X 0 ∈ R 2 + .Proof.Let τ n be the stopping time defined in Theorem 1. Applying the Itô formula to the process V (t, x i (t)) = e t x p i (t), i = 1, 2, p > 0, we obtain for i, j = 1, 2, i = j , Under Assumption 1 there are constants K i (p) > 0, i = 1, 2, such that From ( 7) and ( 8) for the process where c i > 0, i = 1, 2, some constants, which we will define later, we have where S stoch i , i = 1, 2, are the sums of corresponding stochastic integrals in (9).From the Young inequality we have . (10) So, applying (10), we derive the estimate Using the condition a Due to (11) we have Therefore, there is a constant C(p) > 0 such that L(x 1 , x 2 ) ≤ C(p).So, taking the expectation in (9), we obtain Letting n → ∞ leads to the estimate for any inial value X 0 ∈ R 2 + .Theorem 3. If min i=1,2 inf t≥0 (r i (t) − β i (t)) > 0, where r i (t) and β i (t) are defined respectively in (2) and (3), then the solution X(t) to the system (2) with the initial condition X 0 ∈ R 2 + is stochastically permanent.Proof.Using the same arguments as in the corresponding part of the proof of Theorem 3 ( [7]) and the condition min i=1,2 inf t≥0 (r i (t) − β i (t)) > 0, we can conclude that for sufficiently small 0 < θ < 1 and sufficiently small λ = λ(θ ) > 0 there exists a constant K > 0 such that where U i (t) = 1/x i (t), i = 1, 2. Let τ n be the stopping time defined in Theorem 1.
Then by integrating (13) and taking the expectation we have If n → ∞, then we obtain the estimate From ( 14) we have lim sup From ( 12) and (15) by the Chebyshev inequality we can derive, that for arbitrary ε ∈ (0, 1), there are positive constants H = H (ε) and h = h(ε) such that lim inf The proof is completed.Definition 3 ([11]).The solution X(t) = (x 1 (t), x 2 (t)), t ≥ 0, to the system (2) is said to be strongly persistent in the mean if for every initial data X 0 > 0, we have lim inf t→∞ 1 t t 0 x i (s)ds > 0 a.s., i = 1, 2. Therefore, the solution X(t) to the system (2) with the initial condition X 0 ∈ R 2 + will be strongly persistent in the mean.
Proof.For the system (2) by the Itô formula, we have for i, j = 1, 2, i = j , where the martingales have quadratic variation The rest of the proof is the same as the proof of Theorem 6 in [7].The proof is completed.
inf 11 a inf 22 , we can choose constants c i , i = 1, 2, such that a t→∞