Stability , Chaos and Estimation of the Unknown Parameters of Habitat Destruction Model with Prey-Predator-Top Predator

This paper is devoted to study the problem of stability, chaos behavior and parameters estimation of the habitat destruction model with three-species (prey, predator and top-predator). The mathematical formula of the model and its proposed interactions are presented. Some important special solutions of systems are discussed. The stationary states of the model are derived. Local stability conditions for the stationary states are derived. Furthermore, the chaotic behavior of the model is discussed and presented graphically. Using Liapunov stability technique, the dynamic estimators of the unknown probabilities and their updating rules are derived. It is found that, the control laws are non-linear functions of the species densities. Numerical illustrative examples are carried out and presented graphically.


Introduction
The habitat destruction is very important topic, which has received considerable attention in the last years.It is sure that there is a direct correlation between the preservation of the environment, human civilization and progress, and vice versa.Some believe that the chaos is behind each destruction in the environment.In contrast, the others believe that the chaos provides a chance or space for the great changes and developments.It's for sure that, the chaos is one of the main causes of the environmental destruction, but it helps to achieve the needed balance in the environment.So, there is a controversy about whether the chaos has a positive role in the world or not.The balance of the environment comes also, from the natural and bounded coexistence between the organisms.The future of the habitat should have a great interest.There are some studies that have focused on this, for example see [1].
When the habitat is no longer able to provide appropriate conditions for the life of its organisms, we can say that, the environment reached to the destruction stage, which is considered as an important factor causing extinction.It effects on various species not only directly but also indirectly.The accumulation of the destruction in a local area increases the risk of extinction in a bigger area and the increasing of the habitat destruction caused an imbalance in the densities of the populations of prey and predator.So, we can say that, the causal relation between species extinction and the habitat destruction is very complicated [2] [3] [4].
The variation of the total population size is determined by birth and the mortality rates [5].In a closed population (i.e., no immigration) of size N, the change in population size for a change in time is given by ( ) , where b, d and l = b − d are the birth, death and the growth rates per individual respectively.It is clearly that if l > 0, the population grows; if l < 0, the population heading towards extinction; if l = 0, no change in the population size.Recent models of competition indicate that, the effect of habitat destruction on coexisting between preys and predators is dependent on the ratio of extinction risk due to predation and prey colonization rate [6] [7].
A lattice model represents the motion of a network of particles, where the motion is produced by forces acting between the neighboring particles.Lattice models also are used to simulate the structure of polymers and can exhibit its dynamic behaviors.Mathematically, the motion is governed by a system of ordinary differential equations.The behavior of the particles depends on the precise nature of the interaction.The main problems which have been treated involve the long time behavior of the system [8].
The results of this study are extend to the results that introduced by [2] [3] where it found that the increase of habitat destruction leads to reduction in the chance of coexistence of prey and predator, different patterns of extinction for the species (prey, predator and top predator), and decrease in the oscillations in the densities of both species (prey and predator).Also, this study is considered as a complement for the results of the study [9] that concentrates on two-Species model, (prey-predator) which found that the system in general is unstable.In other hands, the chaotic behavior of the uncontrolled system has seemed clear through the graphs of its limit cycles and some attractors.The dynamic estimators of the unknown parameters and its updating rules over time are derived from the conditions of the asymptotic stability of the system around its steady states.
Al-Mahdi and Khirallah have studied stability and bifurcation analysis of a model of cancer [10] [11].Alwan has studied the stability and behavior for the model of stochastic lattice gas of prey-predator model with pair-approximation.She found that this system has a chaos behavior and she has derived the estimators of the unknown parame-ters and the updating [12] The considered system contains three-species, which including prey, predator and top-predator.The assumed interactions and mathematical form of this model will be presented.Some special solutions of this system will be discussed.Chaos and linear stability will be studied.The dynamic estimators of the unknown parameters of this model and its updating rules will be calculated.
This paper has the following structure.In Section 2, the lattice model and its assumed interactions are discussed and presented graphically.The mathematical formula model will be presented.Section 3 is devoted to study the chaos and linear stability analysis of the system.In Section 4, it is the analysis of some special solutions.Estimation of the unknown parameters and the updating rules are derived in Section 5.In Section 6, numerical solutions are derived and presented graphically.Finally, Conclusions are provided in Section 7.

The Lattice Model
In this section, the assumed interactions of the proposed model will be discussed and the mathematical formula will be presented.
Each site in the lattice can be labeled by X, Y, Z or O, where X, (or Y, Z) is the site occupied by prey (or predators), and O represents the vacant site.The assumed interactions of this model are given by [2]: The above interactions respectively represent two kinds of predation with probability 1, reproduction of prey with probability r 0 and the deaths of the three-species prey, predator and top-predator, with probabilities r 1 , r 2 , r 3 respectively.The transition matrix that presents the previous processes can be written as the following: Also, these interactions can be represented by p c = 0.5 [1].In contrast, when p takes an extremely small value, no barriers may connect with each other.Thus, p can be used as measures of the intensity of habitat destruction.

The Mathematical Formula of the Model
The population dynamics of three-species model is described by the mean-field theory (MFT) [2], as the following form: ( ) ( ) where the dot represents the derivative with respect to time, P 0 , P 1 , P 2 and P 3 are the densities of O, X, Y and Z respectively, r i , i = 0, 1, 2, 3 are four unknown probabilities, which represent the probabilities of reproduction of prey, deaths of the prey, predator and top-predator, respectively.The probability p takes a non-critical value as explained previously.
The state variables of this system are subject to the following relationship Therefore, for simplicity in our study, one can take just three equations instead of four, which are considered sufficient to demonstrate the chaotic behavior of the system and to estimate the unknown probabilities.Accordingly, using the relation in Equation (3), one can reduce the model in Equations ( 2) to the following model:

Chaos and Linear Stability Analysis
This section is devoted to study the local stability and the chaos behavior near the steady states of the three-spice model in Equations ( 4).The stationary states of the system in Equations ( 4) will be derived.The biological conditions for these steady states will be discussed.In addition, the instability conditions will be derived.Furthermore, some numerical solutions for the uncontrolled system and its instability, chaotic behavior and different attractors will be presented graphically.

The Stationary States
In this subsection, the stationary states will be derived the biologically conditions will be discussed.The stability study of the stationary states will be discussed later.
The stationary probabilities of the three-species model in Equations ( 4) can be obtained by setting the following simultaneous equations 0, 1, 2,3 By solving this system analytically one can see the following stationary states.a) The first stationary probability of the system in Equations ( 4) is given by: ( ) This solution corresponds to a vacuum-absorbing state.
b) The second stationary probability of the system is given by: Clearly, the density of the top-predator in this solution is negative.Due to that, this solution is biologically inadmissible.
c) The third stationary probability of the system is where The necessary condition for solution in Equation ( 8) to be biologically admissible is d) The fourth stationary probability of the system (4) is given by: ( ) ( ) The necessary condition for this solution to be biologically admissible is ( ) The fifth stationary probability of the system (4) is given by: The necessary condition for this solution to be biologically admissible is ( ) ( ) It is worth mentioning that, the condition in ( 13) is a sufficient condition to be S 4 and S 3 be biologically admissible because ( ) ( ) ( )

Linear Stability Analysis
In this subsection, we will discuss the stability analysis of the three-species model in Equations ( 4) near its steady states.We will calculate the Jacobian matrix of the system at its steady states and calculate the eigenvalues corresponding to each steady state.Then, we will determine whether these eigenvalues contain at least a positive root as a sufficient condition for instability [13] [14].Also, we will seek to determine the instability conditions for each steady state if needed.The Jacobian matrix L of the system in Equation ( 4) evaluated at the steady state , 1, 2, 3 j P j = is given by: ( ) ( ) ( ) Now let's start by examining the instability of the steady states ,1 1, , 5 i S =  of three-species model.
• The eigenvalues of the Jacobian matrix evaluated at S 1 are given by Clearly, S 1 is unstable if 1 r α > .The steady state S 1 needs further stability analysis if 1 r α < , because in this situation, we will get a critical case.
• The eigenvalues of the Jacobian matrix evaluated at S 2 are given by solving the characteristic equation, which Therefore, the eigenvalues of 2 L are the roots of Equation ( 16) which are given by ( ) ( ) ( ) Due to the existence of a positive root in (17), the steady state S 2 is absolutely unstable.
• The eigenvalues of the Jacobian matrix evaluated at S 3 are given by ( ) For this steady state, if 1 r α < , then 31 0 λ > , which is a sufficient condition for insta- bility of S 3 .but, this condition makes the prey density in S 3 has a negative value, therefore this condition is biologically, not admissible.On other hand, if ( ) 0 λ > .Therefore, this is a sufficient condition in order to be S 3 unstable.
• The eigenvalues of the Jacobian matrix evaluated at S 4 are given by the roots of the following characteristic equation: where ( ) For the linear factor, the necessary condition for its root to be positive value is ( ) ( ) This condition is sufficient for S 4 to be unstable.
It is worth mentioning, for the quadratic polynomial in Equation ( 19), the second coefficient the third coefficient is positive under the condition in (21).Since all the real coefficients of the quadratic polynomial in Equation ( 19) are positive, using the Routh's stability criterion [15], then all the roots of quadratic polynomial are negative if the condition (21) is satisfied.Accordingly, S 4 is unstable under the condition in (21) and the stability decision needs further stability analysis • The eigenvalues of the Jacobian matrix evaluated at S 5 are given by the roots of the following characteristic equation.
( ) where Without solving the polynomial between the rectangular brackets in (22) and using Descartes rule of signs to determine the number of positive zeros of the polynomial, one can find that, the first coefficient is positive, and by simple calculation, it is easy to see that if ) ( ) then, 1 0 B < ; and 3 0 B < .Under this condition, the second coefficient of the polynomial between the rectangular brackets in Equation ( 23) will be positive and the third will be negative.In this case the polynomial has one variation in sign.Therefore, the polynomial has one positive real zeros and the other is negative.So, we can conclude that, S 5 is unstable if the condition in (24) is satisfied.
In the following we will present graphically, the system behavior over time, which illustrate the oscillatory behavior of the system.
In this figure, all state variables of the uncontrolled system in Equations ( 4) that appear in (a, b, c) have an oscillatory behavior.For limit cycles that appear in (d, e, f), all the neighboring trajectories tend to a limit-cycle at time tends to infinity, causing the so-called a stable limit-cycle, which indicates that the system has an oscillatory behavior.
In the following we will present a special behavior of this system which is a periodic orbits which called attractors as in Figure 3. Attractors is a special behavior, which generated by the dynamical system and steadily developing as long as there is no barrier All graphs of uncontrolled densities, limit cycles and attractors in Figure 2 and Figure 3 well agree with the previous linear stability analysis that indicates that the habitat destruction model with three-species has a high chaotic behavior.

Analysis of Some Special Solutions
In this section, we will discuss some special solutions for this model, which give us a clearer perception for this mode.In addition, the absence of one or more of the components of the system has impact on the behavior of the system.• The first special solution occurs when both prey and predator are absent, that is P 1 = 0 and P 2 = 0.In this case, the top-predator decreases exponentially with rate r 3 as indicants in the following growth equation: In Figure 4, the top predator density decreases over time and converges to zero and the empty site density increases over time and converges to one as t converges to infinity.This situation represents a vacuum-absorbing state.• The second special solution occurs in the absence of both prey and predator, that is P 1 = 0 and P 3 = 0.This situation is similar to the previous case, therefore, the predator will be decreased exponentially with rate r 2 , and the empty sites density increases over time to tend to one.This situation represents also, a vacuum-absorbing state.
• Another special solution occurs when both of the predator and top-predator are absent, that is P 2 = 0 and P 3 = 0.In this case the model in Equation ( 4) will be reduced to the contact percolation process [2], that is a one-species model we have where the birth and death processes of species X are assumed, which can be presented as and X O → .In this case the system will take the following reducing form ( ) ( ) By solving this system, we will take the following prey population growth: ( where ( ) Equation ( 27) indicates that, the density of the prey tends to ( ) to infinity.If the death probability of the prey be zero, then the prey population has a logistic growth and all sites of the lattice will be occupied by prey as t tends to infinity [16] [17] to make a prey-absorbing state as shown in Figure 5.
It is useful to notice that, in order the result in (27) be biologically admissible, the following condition must be satisfied: ( ) • Other special solution occurs when the top-predator only, is absent, that is P 3 = 0.
Clearly, in this case the system in Equations ( 4) will be reduced to the subsystem of two-species (Prey-Predator) model, ( ) ( ) which has been studied in [9].It is found that, one absorbing stationary state and two active stationary states.It is found also that, these stationary states are Furthermore, this system showed a chaotic behavior.
• In the absence of prey, P 1 = 0, we will obtain another two-species model (predatortop predator) as given below This subsystem has two steady states, which are ( ) The Jacobian matrix Jsub1 of the subsystem (31) evaluated at the steady sate , 2, 3 k P k = is given by: The Jacobian matrices Jsub11 and Jsub12 at e 11 and e 12 respectively, are given below The eigenvalues of Jsub11 are The eigenvalues of Jsub12 are Clearly, the steady state e 11 is stable, while the steady state e 12 is absolutely unstable.
Therefore, the subsystem (31) is unstable at least in one dimension.
• The following considered special solution occurs when the predator is absent, that is P 2 = 0 In this case, the three-species system will be reduced to a subsystem of twospecies model with prey-top predator as will present below ( ) (37.1) This system has the following two stationary solutions: ( ) The Jacobian matrix Jsub2 of the subsystem in Equation (37) evaluated at its steady states , 1, 3 k P k = is given by: ( ) Therefore, the Jacobian matrices at e 21 and e 22 respectively, are given by: and the eigenvalues of Jsub21 are: and the eigenvalues of Jsub22 are: Clearly, the steady state • The last special solution occurs when no vacant sites in the lattice, that is P 0 = 0.In this case the processes: reproduction of prey, death of prey and death of predator will be disabled, and only the predation process will be remained.As a result of this, the density of predator has the following form Equation (43) shows that, the predator population has a logistic growth and its density converges to one as t converges to infinity, in the same time the prey density converges to zero as t converges to infinity as shown in Figure 6.Therefore all the lattice sites will be occupied by predators to make a predator-absorbing state and the preys will be extinct.

Estimations of the Unknown Probabilities
This section is devoted to derive the estimators of the unknown probabilities of the three-species model of habitat destruction (4).The dynamic estimators of the unknown probability r i , i = 0, 1, 2, 3 and its updating rules will be calculated from the conditions of the asymptotic of the system (4) near its stationary solutions S k , k = 1, 2, •••, 5 with assistance of some feedback variables e 1 , e 2 and e 3 .
At the beginning, we will develop the mathematical formula of the habitat destruction model of three-species with unknown probabilities in (4) to become as follows: ( )( )  The problem now is equivalent to stabilizing of these steady states and determining the updating rules 1 r  of the estimators i r of the system unknown probabilities i r with help of the controllers ( ) For the same values of probabilities, parameters and initial densities above, let us represent the estimators of the unknown probabilities.
In Figure 8, it is shown that, the estimators ˆ, 0,1, 2, r i = converge to the assumed real values r i over time, where the assumed real values are represented by the dotted  lines.It is necessary to know that the very small differences between the assumed real values and the estimated values are be due to the necessary approximations in the numerical method of solution.

Figure 1 .Figure 1 .
Figure 1.The assumed transitions of the habitat destruction model of three-species.

Figure 4 .
Figure 4. Simultaneous evanescence and growth of top-predator and empty sites respectively in the absence case of prey and predator at 3 0 0.33; 3 0.14.P r = =

Figure 5 .
Figure 5. Simultaneous growth of both of prey and empty densities in the absence case of predator at 1 0 0.23, 1 0.05, 1.2 P r α = = = .

Figure 6 .
Figure 6.Simultaneous growth of the predator site density and the prey site density in the absence of empty sites at the initial value 2 0 0.82 P = .

r
are the estimators of the unknown probabilities r i , i = 0, 1, 2, 3 of the system (4).The model in Equations (44) has unstable stationary solution if 1

Figure 7
Figure 7 indicates that the controlled densities of the three species converge to the steady state over time.The steady states are presented by the dotted lines.For the same values of probabilities, parameters and initial densities above, let us represent the estimators of the unknown probabilities.In Figure8, it is shown that, the estimators ˆ, 0,1, 2, r i = converge to the assumed

Figure 7 .
Figure 7.The controlled densities of (a) prey, (b) predator, and (c) top-predator sites of the habitat destruction model of three-species for the following values of probabilities, parameters and initial densities.

Figure 8 .
Figure 8.The estimator ˆi r of the unknown probabilities ri, i = 0, 1, 2, 3 of the three-species model of the habitat destruction system with three-species prey, predator and top-predator.