Abstract and Applied Analysis

Abstract and Applied Analysis / 2013 / Article
Special Issue

Nonlinear Dynamics in Applied Sciences Systems: Advances and Perspectives

View this Special Issue

Research Article | Open Access

Volume 2013 |Article ID 465183 | https://doi.org/10.1155/2013/465183

Renato Colucci, Daniel Nuñez, "Periodic Orbits for a Three-Dimensional Biological Differential Systems", Abstract and Applied Analysis, vol. 2013, Article ID 465183, 10 pages, 2013. https://doi.org/10.1155/2013/465183

Periodic Orbits for a Three-Dimensional Biological Differential Systems

Academic Editor: Luca Guerrini
Received03 Jun 2013
Accepted07 Jun 2013
Published25 Jun 2013


We study the existence of periodic orbit for a differential system describing the effects of indirect predation over two preys. Besides discussing a generalized version of the model, we present some remarks and numerical experiments for the nonautonomous version of the two models.

1. Introduction

The role of indirect effects in biology has been extensively studied in the last decades (see [17]); for example, in the case of predation (see [8]), the predator can alter the morphology (see [9]) or the behavior of the preys. The preys, in order to reduce the possibility of contacts with the predators, could modify their normal conduct by reducing their activity or by hiding themselves for long time. There are many types of indirect effects (see [3] for a detailed discussion). Another interesting case is the refuge indirect effect (see [10] for an example); anyway it is of great interest trying to describe the indirect interactions in population dynamics. In [11] a model was proposed including indirect effects, modeling the effects of predator Daphnia over two groups of Phytoplankton of different morphology (see [12, 13]), having Phosphorous as resource (see [13] or [7]). The system takes the following form  where , , , , , are positive parameters and where .

In the previous system represents the density of population of a predator (Daphnia or Zooplankton) that predates the preys (Phytoplankton) and that are of different size, in particular being of a smaller size than . The variable represents the amount of resources (Phosphorous) for the preys and . The system admits a constant of motion , and then we can reduce the degree of freedom of the problem by fixing a value of the first integral, . The system can be rewritten in the following way: The dynamics of system (2) have been studied in [14] by using Persistence Theory (see [15, 16]) and in particular an acyclicity approach. It has been shown that in absence (i.e., ) of the terms that describe indirect effects, the system (2) does not admits coexistence of the three populations, while for and for an appropriate choice of the parameters, the system admits coexistence. The paper [14] suggests the importance of indirect effects in describing cases of coexistence with the parameter being a bifurcation parameter for coexistence. In particular under the hypotheses the system is uniformly persistent. Under the previous hypotheses the system admits an attractor that lies at positive distance from and inside the set In particular the boundary fixed point and the segment of fixed point are all instable and their stable and central manifolds are all on . Under the hypotheses (3)-(4) the set always contains an interior fixed point as follows: whose stability character has not been investigated in [14]. The matrix of the linearized system at is where and and . By numerical experiments (see [14]) it was pointed out that the attractor should take the following forms: (1)the interior fixed point ; (2)a periodic orbit; (3)a bidimensional (possibly chaotic) set.

In this work we are concerned with the proof of existence of periodic orbits (see [17, 18] for examples in the classical case), and we will use Hopf bifurcation theory in order to prove existence of periodic orbits. Since the system depends on 7 parameters, the conditions to have Hopf bifurcation are quite complicated and for sake of simplicity we also present a particular case in which all the parameters are fixed but . In general cases we can use the following remark to facilitate our analysis.

Remark 1. The trace of the matrix satisfies the following inequality then at least one eigenvalue has negative real part and as a consequence fold-Hopf bifurcation is not possible.

The rest of the paper is organized as follows. In Section 2 we present general condition to have Hopf bifurcation, while in Section 3 we give an example fixing all the parameters except for . In Section 4 we consider a nonautonomous version of the model and present several numerical experiments in order to study the structure of the attractor. In the last case we discuss a generalized version of the model, both autonomous and nonautonomous, and present some numerical experiments.

2. Condition for Hopf Bifurcation

We consider the characteristic polynomial of the matrix as follows: where A necessary condition for Hopf bifurcation is that each of the previous coefficients has to be positive. It easy to verify that , , if and the last inequality follows from hypothesis (3).

From Hurwitz-Routh criterion we need to study the sign of the term The sign of the previous term depends only on the factor ; in details we have where From the hypothesis (3) we already have ; then if we have that for any , and this implies that the real parts of the eigenvalues are all negative and so the fixed point is asymptotically stable. Moreover, all the persistence solutions converge to the fixed point .

On the contrary if we have that for any . In this case the eigenvalues have all negative real parts and the fixed point is asymptotically stable. If , then becomes negative and the fixed point loses stability. Moreover a periodic orbit appears thanks to Hopf bifurcation theorem (with as the bifurcation value). Recalling that all the boundary fixed points are instable and repulsive (see [14]), we can conclude that the attractor consists of the periodic orbit generated by Hopf bifurcation. Since the system is uniformly persistent under hypotheses (3)-(4), we expect that the periodic orbit is persistent. We summarize as follows.

Proposition 2. Suppose that (3) and (4) hold. Then the matrix has one negative and two purely imaginary eigenvalues if and only if and condition (18) holds.

Thanks to Hopf bifurcation theorem, a periodic orbit appears. In order to illustrate this result, in the next section we will consider an example of existence of limit cycle.

3. An Example of Existence of a Limit Cycle

In this section we present an example of the results of Section 2. We prove the existence of a limit cycle by considering as a bifurcation parameter and fixing the values of the others parameters.

We set We will choose in order to satisfy hypothesis (3) while hypothesis (4) is already satisfied as The fixed point has the following coordinates: and it then exists for which is verified by hypothesis (3). The Jacobian matrix at point is the following: The characteristic polynomial of the matrix takes the form where and where the three inequalities follow from hypothesis (3). We use Hurwitz-Routh criterion to study the sign of the real parts of the eigenvalues. In particular we are interested in the sign of the term which is positive for and negative for . Then Hopf bifurcation occurs at and this gives the existence of a family of periodic orbits.

In particular if satisfies hypothesis (3) and , we have three eigenvalues with negative real part and then the fixed point is stable and (see Figure 1). In the case in which we have one eigenvalue with negative real part and two eigenvalues with positive real part. In this case the attractor is a periodic orbit (see Figure 2). For the eigenvalues are The period of the limit cycle is , where as represented in Figure 3.

We derive the equation of the central manifold in order to investigate the character of the bifurcation. We rewrite the systems with parameters as in (19) as follows: We consider the following change of coordinates: and we rewrite the system in matrix form where the matrix is calculated at . Let , , be the eigenvectors of . Consider the following transformations: and the change of variable where is the matrix with vectors , , as columns. Then the system becomes that is where We consider a 2-dimensional center manifold of the type Using the first equation of the system, we get while using (37) and the equations of the systems, we get Equaling the coefficients, we get with solutions The system restricted on the central manifold takes the form where , . As remarked in the previous section the periodic orbit should persist and should be stable for .

The first Lyapunov coefficient [19, page 99] can be calculated as and this confirm that the Hopf bifurcation is supercritical.

4. The Nonautonomous Case

Since indirect effects can be of seasonal type (see [14]), it is interesting to study the case in which the constant is replaced by the time-dependent term: We consider values of the parameters for which a periodic orbit exists and make the substitution (44) in the indirect effects terms of the system We consider some numerical experiments which show that not only the periodic orbit persists but also sort of period bifurcation phenomena occur.

For simplicity we fix all the parameters but as follows: In Figure 4 we represent the solution for , , while in Figure 5 for , , . We note that as increases, the period of the solution increases, and in Figure 6 we represent the solutions , , as function of time for showing different frequencies of oscillations. It would be interesting to analyze this phenomenon in details together with the possibility to have chaotic motion. In Figure 7 we compare the period and oscillations of the function for .

5. A Generalized Model

In this section we consider a slighter modification of the model, both in the autonomous and nonautonomous cases. If we sum the second and the third equations of the system, we get then the indirect effects terms disappear and so indirect effects do not change the total amount of the preys but only the proportion between them. A simple generalization of the model consists in considering different values of the parameter , that is, and , respectively, in the second and the third equations of the system The dynamical behavior of the system (48) is very similar to the case (see [14] for details). In particular it has the same boundary fixed points with the same stability character. Then we can state the following without giving the details of computation.

Theorem 3. The system (48) admits coexistence for any positive solution if the following two conditions are satisfied:

We note that the previous conditions do not depend on the parameter . The differences between the two models appear in the coordinate of the fixed point as follows: where . In particular, as for the original system, if the hypotheses of the previous theorem are satisfied, the interior fixed point always exists.

The parameter plays a role only in the stability of the point and as a consequence in the structure of the attractor. By the same method of Section 2 it is possible to derive a condition similar to (18) for Hopf bifurcation.

Proposition 4. Suppose that (49) hold. Then the matrix has one negative and two purely imaginary eigenvalues if and only if and where

We note that inequality (51) depends on and .

From the previous discussion we remark that the behavior of the system (48) is qualitative, the same of that of the system (2).

Numerical simulations suggest that the nonautonomous version of system (48) presents peculiar features. We put and , respectively, in the second and third equations of the system. If , the numerical experiments showed a behavior of the system similar to that described in the previous section (i.e., the case ). If , we have that then, since the last term is nonnegative, we expect a more complicated behavior of the solutions.

We present some simulations to described the dynamics of this case. We consider the following values of the parameters: For the attractor is a periodic orbit, and in Figure 8 we represent the case . The most interesting case is when where a new structure appears; see Figures 9 and 10.

For higher value of we obtain an attractor similar to what we have found in the previous section but with a more complicated geometry. In Figures 11 and 12 we represent, respectively, the attractor and the function , , for .

The numerical experiments suggest that indirect effects encourage an oscillatory behavior of the system and, as a consequence, they are useful to describe real cases of coexistence (see [20]). It could be interesting to investigate in details the nonautonomous case, in particular the phenomenon of period doubling and the existence of chaotic motion (see [2022] for examples for the classical case).


This work was partially supported by the research group EMAP (Grupo de investigaciòn en matemàatica y estadistica aplicada) of Pontificia Universidad Javeriana of Cali, Colombia.


  1. B. Bolker, M. Holyoak, V. Křivan, L. Rowe, and O. Schmitz, “Connecting theoretical and empirical studies of trait-mediated interactions,” Ecology, vol. 84, no. 5, pp. 1101–1114, 2003. View at: Google Scholar
  2. D. Cariveau, R. E. Irwin, A. K. Brody, L. S. Garcia-Mayeya, and A. Von Der Ohe, “Direct and indirect effects of pollinators and seed predators to selection on plant and floral traits,” Oikos, vol. 104, no. 1, pp. 15–26, 2004. View at: Publisher Site | Google Scholar
  3. B. A. Menge, “Indirect effects in marine rocky intertidal interaction webs: patterns and importance,” Ecological Monographs, vol. 65, no. 1, pp. 21–74, 1995. View at: Google Scholar
  4. W. E. Snyder and A. R. Ives, “Generalist predators disrupt biological control by a specialist parasitoid,” Ecology, vol. 82, no. 3, pp. 705–716, 2001. View at: Google Scholar
  5. J. T. Wootton, “Indirect effects, prey susceptibility, and habitat selection: impacts of birds on limpets and algae,” Ecology, vol. 73, no. 3, pp. 981–991, 1992. View at: Google Scholar
  6. M. R. Walsh and D.N. Reznick, “Interactions between the direct and indirect effects of predators determine life history evolution in a killifish,” Proceedings of the National Academy of Sciences of the United States of America, vol. 105, no. 2, pp. 594–599. View at: Google Scholar
  7. D. O. Hessen, T. Andersen, P. Brettum, and B. A. Faafeng, “Phytoplankton contribution to sestonic mass and elemental ratios in lakes: implications for zooplankton nutrition,” Limnology and Oceanography, vol. 48, no. 3, 2003. View at: Google Scholar
  8. J. Estes, K. Crooks, and R. Holt, Ecological Role of Predators, vol. 4 of Enciclopedia of Biodiversity, Academic Press, 2001.
  9. V. Lundgren and E. Granéli, “Grazer-induced defense in Phaeocystis globosa (Prymnesiophyceae): influence of different nutrient conditions,” Limnology and Oceanography, vol. 55, no. 5, 2010. View at: Publisher Site | Google Scholar
  10. M. F. Carusela, F. R. Momo, and L. Romanelli, “Competition, predation and coexistence in a three trophic system,” Ecological Modelling, vol. 220, no. 19, pp. 2349–2352, 2009. View at: Publisher Site | Google Scholar
  11. Indirect effects affects ecosystem dynamics, 2011, http://www.ictp-saifr.org/.
  12. R. Margalef, “Life forms of Phytoplanktos as survival alternative in an unstable environment,” Oceanologica Acta, vol. 134, 1978. View at: Google Scholar
  13. O. Sarnelle, “Daphnia as keystone predators: effects on phytoplankton diversity and grazing resistance,” Journal of Plankton Research, vol. 27, no. 12, pp. 1229–1238, 2005. View at: Publisher Site | Google Scholar
  14. R. Colucci, “Coexistence in a one-predator, two-prey system with indirect effects,” Journal of Applied Mathematics, vol. 2013, Article ID 625391, 13 pages, 2013. View at: Publisher Site | Google Scholar
  15. H. L. Smith and H. R. Thieme, Dynamical Systems and Population Persistence, American Mathematical Society, 2011. View at: MathSciNet
  16. G. Butler, H. I. Freedman, and P. Waltman, “Uniformly persistent systems,” Proceedings of the American Mathematical Society, vol. 96, no. 3, pp. 425–430, 1986. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
  17. S. M. Shahruz and D. A. Kalkin, “Limit cycle behavior in three- or higher-dimensional non-linear systems: the Lotka-Volterra example,” Journal of Sound and Vibration, vol. 246, no. 2, pp. 379–386, 2001. View at: Publisher Site | Google Scholar | MathSciNet
  18. P. Van Den Driessche and M. L. Zeeman, “Three-dimensional competitive Lotka-Volterra systems with no periodic orbits,” SIAM Journal on Applied Mathematics, vol. 58, no. 1, pp. 227–234, 1998. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
  19. Y. A. Kuznetsov, Elements of Applied Bifurcation Theory, Springer, 1998. View at: MathSciNet
  20. J. Hulsman and F. J. Weissing, “Biodiversity of plankton by species oscillations and chaos,” Nature, vol. 402, pp. 407–410, 1999. View at: Publisher Site | Google Scholar
  21. J.-M. Ginoux, B. Rossetto, and J.-L. Jamet, “Chaos in a three-dimensional Volterra-Gause model of predator-prey type,” International Journal of Bifurcation and Chaos in Applied Sciences and Engineering, vol. 15, no. 4, pp. 1689–1708, 2005. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
  22. A. Klebanoff and A. Hastings, “Chaos in one-predator, two-prey models: general results from bifurcation theory,” Mathematical Biosciences, vol. 122, no. 2, pp. 221–233, 1994. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet

Copyright © 2013 Renato Colucci and Daniel Nuñez. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles