The Stationary Distribution and Extinction of Generalized Multispecies Stochastic Lotka-Volterra Predator-Prey System

This paper is concerned with the existence of stationary distribution and extinction for multispecies stochastic Lotka-Volterra predator-prey system. The contributions of this paper are as follows. (a) By using Lyapunov methods, the sufficient conditions on existence of stationary distribution and extinction are established. (b) By using the space decomposition technique and the continuity of probability, weaker conditions on extinction of the system are obtained. Finally, a numerical experiment is conducted to validate the theoretical findings.


Introduction
The dynamic relationship between the predators and the preys has long been and will continue to be one of the dominant themes in both ecology and mathematical ecology due to its universal existence and importance [1].The classic predator-prey model is the Lotka-Volterra model, governed by the following differential equation: where () and () denote the prey and predator population size, respectively, at time .For the prey component, the parameters  and  are the fixed growth and mortality rates, respectively.For the predator component, the parameters  and  are the fixed growth and mortality rates, respectively.Since then, variants of the two-species Lotka-Volterra system have been frequently investigated to describe population dynamics with predator-prey relations; see, for example, [2][3][4].
Recently, the multispecies predator-prey systems have received a great deal of research attention since they took the differences among individual growth and mortality into account (see [5][6][7][8]).In order to understand the nature of the competitive interactions and relationships between predator and prey, Yang and Xu [8] considered the following periodic -prey and −-predator Lotka-Volterra differential system with periodic coefficients: where   (),  = 1, . . ., , denotes the density of prey species at time  and   (),  =  + 1, . . ., , denotes the density of predator species at time .Under the assumption that   () > 0 ( = 1, . . ., ),   () ≥ 0 ( =  + 1, . . ., ),   () > 0 ( = 1, . . ., ),   ≥ 0 ( ̸ = ) are continuous periodic functions with a common periodic  > 0, a set of sufficient conditions on the existence and global attractiveness of the periodic solution to system (2) are obtained.Recently, Chen and Shi [5] further considered the almost periodic case of more complicated systems than system (2) under the almost periodic case.By constructing a suitable Lyapunov function, they obtained a set of sufficient conditions which guarantees the existence of a unique globally attractive positive almost periodic solution to the corresponding system.
On the other hand, from the biological point of view, population systems in the real world are inevitably affected by environmental noise.In the past decades, the dynamics of stochastic populations and related topic have received a great deal of research attention (see [9][10][11][12][13][14][15][16][17][18][19][20][21]), since they have been successfully used in a variety of application fields, including biology (see [22][23][24][25][26][27][28]), epidemiology (see [29,30]), and neural networks (see [31][32][33]).More recently, the asymptotic properties of stochastic predator-prey systems have received a lot of attention; the readers can refer to [10,11,34] and the references therein.For example, the dynamics of the density dependent stochastic predator-prey system with different functional response have been studied by Ji and Jiang in [10,11].Vasilova [34] has investigated a stochastic Gilpin-Ayala predator-prey model with timedependent delay, and certain asymptotic results regarding the long-time behavior of trajectories of the solution and sufficient criteria for extinction of species for a special case of the considered system are given.
In this paper, considering the effect of environmental noise, we introduce stochastic perturbation into the growth rate of the prey and the predator in system (2) and assume that parameters   and   are constant.Then we obtain the following -prey and  − -predator stochastic Lotka-Volterra system with constant coefficients: where () = ( 1 (),  2 (), . . .,   ()) is an -dimensional Brownian motion and  2  will be called the noise intensity.Throughout this paper, we always assume that the following hypothesis holds: In the study of stochastic population systems, extinction and existence of stationary distribution are two important and interesting properties, respectively, meaning that the population system will die out or the distribution of the solution converges weakly to the probability measure in the future, which have received a lot of attention (see [12,[35][36][37]).Then one question arises naturally: under what condition can system (3) have a stationary distribution and become extinct, respectively?This issue constitutes the first motivation of this paper.
In addition, the existing literatures (see [10,35]) show clearly that if the noise intensity of every prey species is more than twice the corresponding intrinsic growth rate, the population will become extinct exponentially.Then one interesting question is as follows: What will happen if the noise intensity equals twice the intrinsic growth rate?Thus, the second purpose of this paper is to solve this interesting problem.
The organization of the paper is as follows.Section 2 describes some preliminaries.The main results are stated in Sections 3 and 4. In Section 3, sufficient conditions are obtained under which there is a stationary distribution to system (3).By utilizing some novel stochastic analysis techniques, sufficient criteria for ensuring the extinction of system (3) are obtained in Section 4. Section 5 provides some numerical examples to check the effectiveness of the derived results.Conclusion is made in Section 6.

Notation
Throughout this paper, unless otherwise specified, let (Ω, F, {F  } ≥0 , P) be a complete probability space with a filtration {F  } ≥0 satisfying the usual conditions (i.e., it is increasing and right continuous while F 0 contains all P-null sets).Let () = ( 1 (),  2 (), . . .,   ()) be an n-dimensional Brownian motion defined on the probability space.If  ∈  × is symmetric, its largest and smallest eigenvalues are denoted by  max () and  min ().Let  * = ( * 1 ,  * 2 , . . .,  *  ) be the positive equilibrium of the corresponding deterministic predator-prey system to system (3), that is, the solution to the following equation: *  = 0,  =  + 1, . . ., . ( In the same way as Zhu and Yin [38] and Liu et al. [39] did, we can also show the following result on the existence of global positive solution.

Lemma 1. Suppose that condition (4) holds; then one has the following assertions:
(i) For any given initial value  0 ∈   + , there is a unique solution (,  0 ) to system (3) and the solution will remain in   + with probability 1; namely, for any  0 ∈   + .(ii) For any given initial value  0 ∈   + and any  > 0, almost every sample path of   (,  0 ) is locally but uniformly Holder continuous.

Stationary Distribution
In this section, we mainly show that system (3) has a stationary distribution.Let us give a lemma that will be used in the following proof.Let () be a homogeneous Markov process in   ⊂   described by the following stochastic equation: The diffusion matrix is Lemma 3 (see [41]).One assumes that there is a bounded open subset  ⊂   with a regular (i.e., smooth) boundary such that its closure  ⊂   , and (i) in the domain  and some neighborhood therefore, the smallest eigenvalue of the diffusion matrix () is bounded away from zero; (ii) if  ∈   \ , the mean time  at which a path issuing from  reaches the set  is finite, and sup ∈    < +∞ for every compact subset  ∈   and throughout this paper one sets inf 0 = ∞.
One then has the following assertions: (1) The Markov process () has a stationary distribution (⋅) with density in   .
(2) Let () be a function integrable with respect to the measure (⋅).Then Remark 4. The proof is given by [41] in detail.Precisely, the existence of a stationary distribution with density is obtained in Theorem 4.3 on pp.117.The ergodic property is referred to Theorem 4.2 on pp.110.To validate (i), we can directly show that  min {()} > 0. To validate (ii), it suffices to prove that there is some neighborhood  and a nonnegative  2function  such that, for any  ∈   \ , L() is negative (for details refer to [42], pp.1163).

Extinction
Extinction is one of the most basic questions that can be studied in the population dynamics, which means the population system will die out.Most of the time we need to know the extinction rate of the species for which we have to make a suitable policy in advance and to make useful measures to protect them from becoming extinct.
Proof.Let   () =   (,  0 ) for simplicity.To make the proof clear, we are going to divide it into four steps.The first step and the third step are to show the least upper bound of exponential extinction rate for the top  preys and the predators of system (3), respectively.The second step is to show the extinction for the bottom  −  preys of system (3) in the case of  2  = 2  ,  =  + 1, . . ., .The fourth step is to accomplish the proof of assertions (i)-(iii) based on the proof of Steps 1-3.
For  = 1, 2, . . ., , letting  ∞ in (21) yields that lim sup Step 2. The main aim of this step is to show lim  → ∞   () = 0 a.s.for  =  + 1, . . ., .When  2  = 2  , (21) turns to be the following form: According to the convergence of ∫ ∞ 0   (), we can decompose the sample space into two exclusive events spaces as follows: On the other hand we can divide the sample space into three mutually exclusive events as follows: From the above, the proof of lim  → ∞   () = 0 a.s. is equivalent to show  ,1 ⊂  ,3 a.s. and  ,2 ⊂  ,3 a.s.Now we give the process in two parts.
Remark 8.Note that, for  = , system (3) becomes the following classic stochastic Lotka-Volterra competitive systems, which have received much attention (see [12,36,38]): And condition (4) becomes the following form: Thus, by Theorem 7, we have the sufficient conditions on extinction for system (45) as Corollary 9.
By using some novel stochastic analysis techniques, we point out that species  is still extinct when  2 = 2  .In comparison with Theorem 4.1 in [12], the conditions imposed on the extinction are weaker.

Example and Simulations
In this section, we present a numerical example to illustrate the usefulness and flexibility of the theorem developed in the previous section.
But we cannot see the value of the extinction rate of the three populations clearly.So we give Figure 3 to show (log  ())/ for  = 1, 2, 3.According to Theorem 7, we can compute that, for  = 1, the growth rate is −0.05220 (it is said that the extinction rate is 0.05220).By the same method, we can know that the extinction rate is 0.00000 for  = 2 and the extinction rate is 0.10125 for  = 3.

Conclusion
This paper is devoted to the existence of stationary distribution and extinction for multispecies stochastic Lotka-Volterra predator-prey system.Firstly, by applying Lyapunov methods, sufficient conditions for ensuring the existence of stationary distribution of the system are obtained.Secondly, some novel techniques have been developed to derive weaker sufficient conditions under which the system is extinctive.Finally, numerical experiment is provided to illustrate the effectiveness of our results.