A Stochastic Two Species Competition Model : Nonequilibrium Fluctuation and Stability

The object of this paper is to study the stability behaviours of the deterministic and stochastic versions of a two-species symmetric competitionmodel. The logistic parameters of the competitive species are perturbed by colored noises or Ornstein-Uhlenbeck processes due to random environment. The Fokker-Planck equation has been used to obtain probability density functions. Here, we have also discussed the relationship between stability behaviours of this model in a deterministic environment and the corresponding model in a stochastic environment.


Introduction
In recent years, scientists have become increasingly aware of the fact that most natural phenomena do not follow strictly deterministic laws but rather oscillate randomly about some average behaviour.This is especially true in the ecological models where environmental influences should be taken as stochastic.Many models in theoretical ecology take the Volterra-Lotka model of interacting species as a starting point.In an improved analysis, the influence of the environment has to be taken into account.This is often done by arbitrarily augmenting the deterministic equations with stochastic terms or taking the environmental parameters as time dependent and rapidly varying.In both cases, the corresponding stochastic properties have to be postulated 1 .In 1995, Renshaw mentioned that the most natural phenomena do not follow strictly deterministic laws but rather oscillate randomly about some average so that the deterministic equilibrium is not an absolutely fixed state; instead, it is a "fuzzy" value around which the biological system fluctuates.In fact, randomness or stochasticity plays a vital role in the structure and function of biological systems.The environmental factors are time dependent, randomly varying and should be taken as stochastic.In ecology, we have two types of stochasticity, namely, the demographic stochasticity and the environmental stochasticity 2, 3 .Both types of stochasticity play a significant part in the realistic dynamical modelling of ecosystems.Several systematic procedures have been developed to obtain the stationary probability density function for the response of general single-degree-of-freedom nonlinear oscillators and nonlinear dynamical systems under parametric and external Gaussian white noise excitations 4-7 .
In this paper, we have studied the relationship between stability behaviours of a twospecies symmetric competition model in which environmental parameters are prescribed constant deterministic environment and the corresponding model in which these parameters have an element of random fluctuation stochastic environment , where it is assumed that the fluctuations in the environment manifest themselves as random fluctuations in the logistic parameters about some mean value, as they usually do in nature 8, 9 .

Two-Species Symmetric Competition: Stochastic Differential Equations
Two or more species-populations compete for the same limited food source or in some way the growth of each is inhibited by members both of its own and of the other species.For example, competition may be for territory which is directly related to food resources.Some interesting phenomena have been observed from the study of practical competition models 10, 11 .Here, we will discuss a competition model which demonstrates a fairly general principle which is observed to hold in nature.Assume that the two-species symmetric competition model in a deterministic environment satisfies the following deterministic differential equations: where N 1 t , N 2 t denote the number or density of individuals of two competing species.Thus, it is assumed that the per capita growth rate of each population at any instant is a linear function of the sizes of the two competing populations at that instant.Each population would grow logistically if it were alone with a constant environmental parameter K 0 > 0 and a parameter α > 0 which measures the symmetric competition between the two species i.e., it measures the degree to which the presence of one species affects the growth of the other species , we assume that α < 1.
Let us rewrite the system of differential equations 2.1 in the following form:

2.2
The deterministic equilibrium populations N 1 , N 2 are given by The 2 × 2 interaction matrix in the neighbourhood of this equilibrium is The two eigenvalues of this interaction matrix are λ 1 −K 0 , λ 2 −K 0 1 − α / 1 α , and the deterministic stability criterion is satisfied since α < 1.
To take into account the random environment, we extend model 2.2 to the form of the following stochastic differential equations: where the fluctuations u i t are colored noises or Ornstein-Uhlenbeck processes which are more realistic noises than white noises.These are extremely useful to model rapidly fluctuating phenomena, because it can be seen by studying their spectra that thermal noises in electrical resistance and climate fluctuations, disregarding the periodicities of astronomical origin and so forth, are colored noises to a very good approximation.These examples support the usefulness of the colored noise idealization in applications to natural systems.The mathematical expectations and correlation functions of the processes u i t are given by where , δ −1 0 are, respectively, the intensity and the correlation time of the colored noise and • represents the average over the ensemble of the stochastic process.The u i t are the solutions of the stochastic differential equations 12 where η i dW i /dt i 1, 2 are standard white noises characterized by where δ t denotes the Dirac delta function.
Let y i ln N i ; therefore, N i e y i .Hence, dy i dt K 0 u i t − e y i − αe y j , i 1, 2; j 3 − i .

2.9
Put n i y i − y ; therefore, y i n i y , where y ln N .Hence, Therefore, From 2.7 and 2.11 , we have

Stationary Probability Density
The reduced Fokker-Planck equation is given by 13 where x i are the state variables representing the possible values of X i , p s is the stationary probability density, the drift and diffusion coefficients are where

3.3
After some simplifications using the results of Cai and Lin 13 , we have

3.5
By integrating over all values of x 3 , we obtain Similarly, we have

Concluding Remarks
In this paper, we have studied the relationship between stability behaviours of a two-species symmetric competition model for interacting species, in which environmental parameters are prescribed constants deterministic environment , and the corresponding model in which these parameters have an element of random fluctuation stochastic environment .In the stochastic environment, the environmental parameter K 0 has been perturbed by colored noises or Ornstein-Uhlenbeck processes, which are more realistic noises than white noises.Now, in the deterministic environment, the maximum of the eigenvalues of the 2 × 2 interaction matrix is −Λ −K 0 1 α / 1 − α .Therefore, the stability determining quantity is Λ, and the deterministic stability criterion is satisfied, since Λ > 0 as α < 1 .In the stochastic environment, whose random fluctuation has intensity , the stability provided by the population interaction dynamics is again characterized by Λ.It is no longer enough that Λ > 0, for if Λ, population fluctuations are relatively small; in this case, the probability cloud persists for an appreciable time, and the environment is effectively deterministic.For < Λ, but not much less, populations are likely to undergo significant fluctuations, even though they persist for long times.Finally, if Λ, population exhibit large fluctuations in this case, the interaction dynamics provides an ever-weaker stabilizing influence to offset the randomizing σ 2 1 σ 2 2 , which rapidly lead to extinction.These results are in good agreement with those of Samanta 14 , Samanta and Maiti 15 , and Maiti and Samanta 16 .
The correlation time of the colored noise is τ cor δ −1 0 .The limiting process for τ cor → 0 of the colored noise is the Gaussian white noise, and in this case we obtain the well-known results of May 8 .Now, for τ cor → ∞, we see that the population fluctuations are very much small, and in this case, the environment is effectively deterministic.