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.

Keywords

Share and Cite:

Alwan, S. and Hamood Odhah, O. (2016) Stability, Chaos and Estimation of the Unknown Parameters of Habitat Destruction Model with Prey-Predator-Top Predator. Journal of Applied Mathematics and Physics, 4, 2254-2271. doi: 10.4236/jamp.2016.412218.

1. 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 parameters 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.

2. 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] :

(1)

The above interactions respectively represent two kinds of predation with probability 1, reproduction of prey with probability r0 and the deaths of the three-species prey, predator and top-predator, with probabilities r1, r2, r3 respectively. The transition matrix that presents the previous processes can be written as the following:

Also, these interactions can be represented by Figure 1. In the lattice boundary, the barriers are putting as a link between neighboring sites probability p. which called barrier density. The interactions (1.1), (1.2) and (1.3) are assumed occur between neighboring lattice points as presented in Figure 1.

Percolation is the case when the largest set of connected barriers encompasses the entire system. The local habitat destruction is shown when p takes a value near unity, almost all barriers are connected (Percolation case). The probability of the percolation takes a nonzero value when p exceeds a critical point pc. This value is assumed to be

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

pc = 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:

(2.1)

(2.2)

(2.3)

(2.4)

where the dot represents the derivative with respect to time, P0, P1, P2 and P3 are the densities of O, X, Y and Z respectively, ri, 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

(3)

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:

(4.1)

(4.2)

(4.3)

3. 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.

3.1. 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

(5)

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:

(6)

This solution corresponds to a vacuum-absorbing state.

b) The second stationary probability of the system is given by:

(7)

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

(8)

where

(9)

The necessary condition for solution in Equation (8) to be biologically admissible is

(10)

d) The fourth stationary probability of the system (4) is given by:

(11)

The necessary condition for this solution to be biologically admissible is

e) The fifth stationary probability of the system (4) is given by:

(12)

The necessary condition for this solution to be biologically admissible is

(13)

It is worth mentioning that, the condition in (13) is a sufficient condition to be S4 and S3 be biologically admissible because

3.2. 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 of the system in Equation (4) evaluated at the steady state is given by:

(14)

Now let’s start by examining the instability of the steady states of three-species model.

・ The eigenvalues of the Jacobian matrix evaluated at S1 are given by

(15)

Clearly, S1 is unstable if. The steady state S1 needs further stability analysis if, because in this situation, we will get a critical case.

・ The eigenvalues of the Jacobian matrix evaluated at S2 are given by solving the characteristic equation, which is

(16)

Therefore, the eigenvalues of are the roots of Equation (16) which are given by

(17)

Due to the existence of a positive root in (17), the steady state S2 is absolutely unstable.

・ The eigenvalues of the Jacobian matrix evaluated at S3 are given by

(18)

For this steady state, if, then, which is a sufficient condition for instability of S3. but, this condition makes the prey density in S3 has a negative value, therefore this condition is biologically, not admissible. On other hand, if, then. Therefore, this is a sufficient condition in order to be S3 unstable.

・ The eigenvalues of the Jacobian matrix evaluated at S4 are given by the roots of the following characteristic equation:

(19)

where

(20)

For the linear factor, the necessary condition for its root to be positive value is

(21)

This condition is sufficient for S4 to be unstable.

It is worth mentioning, for the quadratic polynomial in Equation (19), the second coefficient is positive because of and. Also, 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, S4 is unstable under the condition in (21) and the stability decision needs further stability analysis

・ The eigenvalues of the Jacobian matrix evaluated at S5 are given by the roots of the following characteristic equation.

(22)

where

(23)

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

(24)

then,; and. 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, S5 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.

4. 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.

Figure 2. The uncontrolled densities of habitat destruction of three-species model and its limit cycles for the probabilities r0 = 0.5, r1 = 0.01, r2 = 0.15, r3 = 0.4 and the barrier density p = 0.3 and the initial densities

Figure 3. Different attractors of the habitat destruction model of the three-species for the following rates, parameters and initial densities:

・ The first special solution occurs when both prey and predator are absent, that is P1 = 0 and P2 = 0. In this case, the top-predator decreases exponentially with rate r3 as indicants in the following growth equation:

(25)

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.

Figure 4. Simultaneous evanescence and growth of top-predator and empty sites respectively in the absence case of prey and predator at

・ The second special solution occurs in the absence of both prey and predator, that is P1 = 0 and P3 = 0. This situation is similar to the previous case, therefore, the predator will be decreased exponentially with rate r2, 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 P2 = 0 and P3 = 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. In this case the system will take the following reducing form

(26.1)

(26.2)

By solving this system, we will take the following prey population growth:

(27)

where

(28)

Equation (27) indicates that, the density of the prey tends to as t tends 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:

(29)

・ Other special solution occurs when the top-predator only, is absent, that is P3 = 0. Clearly, in this case the system in Equations (4) will be reduced to the subsystem of two-species (Prey-Predator) model,

Figure 5. Simultaneous growth of both of prey and empty densities in the absence case of predator at.

(30.1)

(30.2)

(30.3)

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 unstable. Furthermore, this system showed a chaotic behavior.

・ In the absence of prey, P1 = 0, we will obtain another two-species model (predator- top predator) as given below

(31.1)

(31.2)

This subsystem has two steady states, which are

(32)

The Jacobian matrix Jsub1 of the subsystem (31) evaluated at the steady sate is given by:

(33)

The Jacobian matrices Jsub11 and Jsub12 at e11 and e12 respectively, are given below

(34)

The eigenvalues of Jsub11 are

(35)

The eigenvalues of Jsub12 are

(36)

Clearly, the steady state e11 is stable, while the steady state e12 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 P2 = 0 In this case, the three-species system will be reduced to a subsystem of two- species model with prey-top predator as will present below

(37.1)

(37.2)

This system has the following two stationary solutions:

(38)

The Jacobian matrix Jsub2 of the subsystem in Equation (37) evaluated at its steady states is given by:

(39)

Therefore, the Jacobian matrices at e21 and e22 respectively, are given by:

(40)

and the eigenvalues of are:

(41)

and the eigenvalues of are:

(42)

Clearly, the steady state is unstable if and stable if and the stability decision needs more stability analysis if, in which one get a critical case. One can conclude that the subsystem (37) is unstable at least in one direction.

・ The last special solution occurs when no vacant sites in the lattice, that is P0 = 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

(43)

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.

Figure 6. Simultaneous growth of the predator site density and the prey site density in the absence of empty sites at the initial value.

5. 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 ri, i = 0, 1, 2, 3 and its updating rules will be calculated from the conditions of the asymptotic stability of the system (4) near its stationary solutions Sk, k = 1, 2, ∙∙∙, 5 with assistance of some feedback variables e1, e2 and e3.

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:

(44.1)

(44.2)

(44.3)

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

(45)

The problem now is equivalent to stabilizing of these steady states and determining the updating rules of the estimators of the system unknown probabilities with help of the controllers.

The Lyapunov function of the system (4) that will be used can be considered as the following positive definite form

(46)

and its time derivative is given by

(47)

By choosing the controllers as follows:

(48.1)

(48.2)

(48.3)

and by substituting (48) in (47), one can get that the appropriate choice for the updating rules for to be a negative definite function are given by

(49.1)

(49.2)

(49.3)

(49.4)

where and are non-negative control constants.

By substituting the choices (48) and (49) in (47), will be a negative definite function, which proves the asymptotic stability of the solution (45) in Lyapunov sense. Further that, this solution is globally asymptotically stable since V is radially unbounded [18] [19] . And by substituting Equation (48) in the controlled system in Equation (44), in addition Equation (9) we get the following closed system of nonlinear deferential equations following form:

(50.1)

(50.2)

(50.3)

(50.4)

(50.5)

(50.6)

(50.7)

This system can be used to study the dynamics of the habitat destruction model of three-species with unknown probabilities. Since this system is nonlinear, its solution will be numerical as will be seen bellow.

6. Numerical Solution

This section presents numerical solution of the controlled nonlinear system of the habitat destruction model with prey-predator-top-predator, and the estimators of the system unknown probabilities to show how the control for this system is possible. Also, numerical examples for controlled model are carried out for various probabilities values and different initial densities. A numerical solution is presented below

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 Figure 8, it is shown that, the estimators converge to the assumed real values ri over time, where the assumed real values are represented by the dotted

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. The estimator 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.

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.

7. Conclusions

This paper is devoted to study the linear stability, chaos and estimation of the unknown parameters of the habitat destruction model with prey-predator and top-predator. Using the linear stability analysis, it is found that, this system has a high chaotic behavior. 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. The controller laws are found as non-linear functions of the species densities. Numerical examples for the controlled system are carried out and presented graphically.

The chaotic behavior of this system gives us a reason to study the optimal control for this system in future.

Conflicts of Interest

The authors declare no conflicts of interest.

 [1] Ergazaki, M. and Ampatzidis, G. (2012) Students’ Reasoning about the Future of Disturbed or Protected Ecosystems & the Idea of the “Balance of Nature”. Research in Science Education, 42, 511-530. https://doi.org/10.1007/s11165-011-9208-7 [2] Nakagiri, N. and Tainaka, K. (2004) Indirect Effects of Habitat Destruction in Model Cosystems. Ecological Modelling, 174, 103-114. https://doi.org/10.1016/j.ecolmodel.2003.12.047 [3] Szwabinski, J. and Pekalski, A. (2006) Effects of Random Habitat Destruction in a Predator-Prey Model. Physica A, 360, 59-70. https://doi.org/10.1016/j.physa.2005.05.079 [4] Nakagiri, N., Tainaka, K. and Tao, T. (2001) Indirect Relation between Extinction and Habitat Destruction. Ecological Modelling, 137, 109-118. https://doi.org/10.1016/S0304-3800(00)00417-8 [5] El-Gohary, A. (2005) Optimal Control of Stochastic Lattice of Prey-Predator Models. Applied Mathematics and Computation, 160, 15-28. https://doi.org/10.1016/S0096-3003(03)00668-4 [6] Klausmeier, C.A. (1998) Extinction in Multispecies and Spatially Explicit Models of Habitat Destruction. The American Naturalist, 152. https://doi.org/10.1086/286170 [7] Swihart, R., Feng, Z., Slade, N., Mason, D. and Gehring, T. (2001) Effects of Habitat Destruction and Resource Supplementation in a Predator-Prey Metapopulation Model. Journal of Theoretical Biology, 210, 287-303. https://doi.org/10.1006/jtbi.2001.2304 [8] Liggtt, T.M. (1985) Interacting Particle Systems. Springer-Verlag, Inc., New York. [9] Alwan, S. and El-Gohary, A. (2011) Chaos, Estimation and Optimal Control of Habitat Destruction Model with Uncertain Parameters. Computers and Mathematics with Applications, 62, 4089-4099. [10] Al-Mahdi, A. and Khirallah, M. (2016) Bifurcation Analysis of a Model of Cancer. European Scientific Journal, 12, 68-83. [11] Khirallah, M. and Al-Mahdi, A. (2014) Stability Analysis of a Model of Cancer Treatment by Immunotherapy. European Journal of Scientific Research, 121, 161-173. [12] Alwan, S. (2016) Chaos Behavior and Estimation of the Unknown Parameters of Stochastic Lattice Gas for Prey-Predator Model Pair-Approximation. Applied Mathematics, 7, 1765-1779. https://doi.org/10.4236/am.2016.715148 [13] Hacker, T. (1970) Flight Stability and Control. American Elsevier Publishing Company, Inc., New York. [14] Cook, P.M. (1994) Nonlinear Dynamical Systems. Prentice Hall International, Upper Saddle River. [15] Ogata, K. (1997) Modern Control Engineering. Prentice-Hall, Inc., Upper Saddle River. [16] El-Gohary, A. and Bukhari, F. (2003) Optimal Control of Stochastic Prey-Predator Models. Applied Mathematics and Computation, 146, 403-415. https://doi.org/10.1016/S0096-3003(02)00592-1 [17] El-Gohary, A. and Al-Ruzaiza, A. (2007) Chaos and Adaptive Control in Two Prey, One Predator System with Non-linear Feedback. Chaos, Solitons and Fractals, 34, 443-453. https://doi.org/10.1016/j.chaos.2006.03.101 [18] Brogan, W.L. (1982.) Modern Control Theory. Prentice-Hall, Inc., Upper Saddle River. [19] El-Gohary, A. (2002) Global Stabilization of a Rotational Motion of a Rigid Body Using Rotors System. Applied Mathematics and Computation, 133, 297-307. https://doi.org/10.1016/S0096-3003(01)00235-1