Research Article | Open Access
Hopf-Pitchfork Bifurcation in a Phytoplankton-Zooplankton Model with Delays
The purpose of this paper is to study the dynamics of a phytoplankton-zooplankton model with toxin delay. By studying the distribution of the eigenvalues of the associated characteristic equation, the pitchfork bifurcation curve of the system is obtained. Furthermore, on the pitchfork bifurcation curve, we find that the system can undergo a Hopf bifurcation at the positive equilibrium, and we derive the critical values where Hopf-Pitchfork bifurcation occurs.
The study of the dynamical interaction of zooplankton and phytoplankton is an important area of research in marine ecology. Phytoplanktons are tiny floating plants that live near the surface of lakes and ocean. They provide food for marine life, oxygen for human being, and also absorb half of the carbon dioxide which may be contributing to the global warming . Zooplanktons are microscopic animals that eat other phytoplankton. Toxins are produced by phytoplankton to avoid predation by zooplankton. The toxin producing phytoplankton not only reduces the grazing pressure on them but also can control the occurrence of bloom; see Chattopadhyay et al.  and Sarkar and Chattopadhyay . Phytoplankton-zooplankton models have been studied by many authors [4–9]. In , models of nutrient-plankton interaction with a toxic substance that inhibit either the growth rate of phytoplankton, zooplankton, or both trophic levels are proposed and studied. In , authors have dealt with a nutrient-plankton model in an aquatic environment in the context of phytoplankton bloom. Roy  has constructed a mathematical model for describing the interaction between a nontoxic and a toxic phytoplankton under a single nutrient. Saha and Bandyopadhyay  considered a toxin producing phytoplankton-zooplankton model in which the toxin liberation by phytoplankton species follows a discrete time variation. Biological delay systems of one type or another have been considered by a number of authors [10, 11]. These systems governed by integrodifferential equations exhibit much more rich dynamics than ordinary differential systems. For example, Das and Ray  investigated the effect of delay on nutrient cycling in phytoplankton-zooplankton interactions in the estuarine system. In this paper we present a phytoplankton-zooplankton model to investigate its dynamic behaviors. The model we considered is based on the following plausible toxic-phytoplankton-zooplankton systems introduced by Chattopadhayay et al.  where and are the densities of phytoplankton and zooplankton, respectively. , , , , , , and are positive constants. is toxin delay, is the delay kernel and a non-negative bounded function defined on as follows:
For a set of different species interacting with each other in ecological community, perhaps the simplest and probably the most important question from a practical point of view is whether all the species in the system survive in the long term. Therefore, the periodic phenomena of biological system are often discussed [12–16].
The primary purpose of this paper is to study the effects of toxin delay on the dynamics of (1). That is to say, we will take the delay passes through a critical value, the positive equilibrium loses its stability and bifurcation occurs. By studying the distribution of the eigenvalues of the associated characteristic equation, the pitchfork bifurcation curve of the system is obtained. Furthermore, we derive the critical values where Hopf-Pitchfork bifurcation occurs.
The paper is structured as follows. In Section 2, we discuss the local stability of the positive solutions and the existence of Pitchfork bifurcation. In Section 3, the conditions for the occurrence of Hopf-Pitchfork bifurcation are determined.
2. Stability and Pitchfork Bifurcation
In this section, we focus on investigating the local stability and the existence of Pitchfork bifurcation of the positive equilibrium of system (1). It is easy to see that system (1) has a unique positive equilibrium , where where .
Let By the linear chain trick technique, then system (1) can be transformed into the following system: It is easy to check that system (5) has an unique positive equilibrium with provided that the condition holds.
Let , , and ; then system (5) can be transformed into the following system: where
Then linearizing system (6) at is It is easy to see that the associated characteristic equation of system (11) at the positive equilibrium has the following form and thus the characteristic equation of system (5) is given by where Obviously, , , , , and .
From (9), the following lemma is obvious.
Lemma 1. If the condition holds, then is always a root of (9) for all .
Let , , and ; we have the following results.
Proof. Clearly, is a root to (9) if and only if , which means . Substituting into and taking the derivative with respect to , we obtain
Then we can get
For any , by solving (12), we can obtain . If , which means that is a single zero root to (9), and hence the conclusion of (i) follows.
From (11), it follows that Then we get For any , by solving (14), we can obtain . If , which means that is a double zero root to (9), and hence the conclusion of (ii) follows. This completes the proof.
From Lemma 2, we have the following result.
Theorem 3. Suppose that holds if , then, the system (5) undergoes a Pitchfork bifurcation at the positive equilibrium.
3. Hopf-Pitchfork Bifurcation
In the following, we consider the case that (9) not only has a zero root, but also has a pair of purely imaginary roots (), when holds.
Substituting () and into (9) and separating the real and imaginary parts, one can get
It is easy to see from (15) that where
Let . Then (16) can be written as In terms of the coefficient in define by . It is easy to know from the characters of cubic algebraic equation that is a strictly monotonically increasing function if . If and or , but , then has always no positive root. Therefore, under these conditions, (9) has no purely imaginary roots for any and this also implies that the positive equilibrium of system (1) is absolutely stable. Thus, we can obtain easily the following result on the stability of positive equilibrium of system (1).
Theorem 4. Assume that holds and or and or , and . Then the positive equilibrium of system (5) is absolutely stable; namely; is asymptotically stable for any delay .
In what follows, we assume that the coefficients in satisfy the condition , , .Then, according to Lemma 2.2 in , we know that (16) has at least a positive root ; that is, the characteristic equation (9) has a pair of purely imaginary roots . Eliminating in (15), we can get that the corresponding such that (9) has a pair of purely imaginary roots , are given by
Let be the roots of (9) such that when satisfying and . We can claim that In fact, differentiating two sides of (9) with respect to , we get Then where . It follows from the hypothesis that and therefore the transversality condition holds.
Summarizing the previous discussions, we have the following result.
Theorem 6. Suppose that the conditions , , and are satisfied. (i)If and , then the system (1) undergoes a Pitchfork bifurcation at positive equilibrium . (ii)If and , then system (1) can undergo a Hopf-Pitchfork bifurcation at the positive equilibrium .
In this section, we present some particular cases of system (1) as follows: From , we know that the system (24) undergoes a Hopf bifurcation at the positive equilibrium. In this paper, we get the condition that (9) has a zero root and also get the conditions that (9) has double zero roots. Furthermore, we obtain the conditions that (9) has a single zero root and a pair of purely imaginary roots. Under this condition, system (1) undergoes a Hopf-Pitchfork bifurcation at the positive equilibrium. Especially, when , system (24) reduces to We can conclude that the positive equilibrium is locally asymptotically stable in the absence of toxin delay.
The authors are grateful to the referees for their valuable comments and suggestions on the paper. The research of the authors was supported by the Fundamental Research Fund of Henan University (2012YBZR032).
- S. Khare, O. P. Misra, C. Singh, and J. Dhar, “Role of delay on planktonic ecosystem in the presence of a toxic producing phytoplankton,” International Journal of Differential Equations, vol. 2011, Article ID 603183, 16 pages, 2011.
- J. Chattopadhyay, R. R. Sarkar, and A. El Abdllaoui, “A delay differential equation model on harmful algal blooms in the presence of toxic substances,” IMA Journal of Mathematics Applied in Medicine and Biology, vol. 19, no. 2, pp. 137–161, 2002.
- R. R. Sarkar and J. Chattopadhayay, “Occurrence of planktonic blooms under environmental fluctuations and its possible control mechanism—mathematical models and experimental observations,” Journal of Theoretical Biology, vol. 224, no. 4, pp. 501–516, 2003.
- S. Chakraborty, S. Roy, and J. Chattopadhyay, “Nutrient-limited toxin production and the dynamics of two phytoplankton in culture media: a mathematical model,” Ecological Modelling, vol. 213, no. 2, pp. 191–201, 2008.
- K. Das and S. Ray, “Effect of delay on nutrient cycling in phytoplankton-zooplankton interactions in estuarine system,” Ecological Modelling, vol. 215, no. 1–3, pp. 69–76, 2008.
- S. R.-J. Jang, J. Baglama, and J. Rick, “Nutrient-phytoplankton-zooplankton models with a toxin,” Mathematical and Computer Modelling, vol. 43, no. 1-2, pp. 105–118, 2006.
- B. Mukhopadhyay and R. Bhattacharyya, “Modelling phytoplankton allelopathy in a nutrient-plankton model with spatial heterogeneity,” Ecological Modelling, vol. 198, no. 1-2, pp. 163–173, 2006.
- S. Roy, “The coevolution of two phytoplankton species on a single resource: allelopathy as a pseudo-mixotrophy,” Theoretical Population Biology, vol. 75, no. 1, pp. 68–75, 2009.
- T. Saha and M. Bandyopadhyay, “Dynamical analysis of toxin producing phytoplankton-zooplankton interactions,” Nonlinear Analysis: Real World Applications, vol. 10, no. 1, pp. 314–332, 2009.
- J. M. Cushing, Integrodifferential Equations and Delay Models in Population Dynamics, Springer, Berlin, Germany, 1977.
- Y. Kuang, Delay Differential Equations with Applications in Population Dynamics, vol. 191, Academic Press, New York, NY, USA, 1993.
- Y.-H. Fan and L.-L. Wang, “Periodic solutions in a delayed predator-prey model with nonmonotonic functional response,” Nonlinear Analysis: Real World Applications, vol. 10, no. 5, pp. 3275–3284, 2009.
- G.-P. Hu and W.-T. Li, “Hopf bifurcation analysis for a delayed predator-prey system with diffusion effects,” Nonlinear Analysis: Real World Applications, vol. 11, no. 2, pp. 819–826, 2010.
- Z.-P. Ma, H.-F. Huo, and C.-Y. Liu, “Stability and Hopf bifurcation analysis on a predator-prey model with discrete and distributed delays,” Nonlinear Analysis: Real World Applications, vol. 10, no. 2, pp. 1160–1172, 2009.
- X.-P. Yan and W.-T. Li, “Bifurcation and global periodic solutions in a delayed facultative mutualism system,” Physica D, vol. 227, no. 1, pp. 51–69, 2007.
- Z. Zhang and Z. Wang, “Periodic solution for a two-species nonautonomous competition Lotka-Volterra patch system with time delay,” Journal of Mathematical Analysis and Applications, vol. 265, no. 1, pp. 38–48, 2002.
- Y. Song, M. Han, and J. Wei, “Stability and Hopf bifurcation analysis on a simplified BAM neural network with delays,” Physica D, vol. 200, no. 3-4, pp. 185–204, 2005.
Copyright © 2013 Jia-Fang Zhang and Dan Zhang. 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.