Global Hopf bifurcation for two zooplankton-phytoplankton model with two delays

In this paper, we study the global existence of a bifurcating periodic solution for a two zooplankton-phytoplankton model with two delays. First, we demonstrate that the bifurcating periodic solution exists when one delay increases and the other delay remains unchanged. Second, we give simulation to describe the bifurcating periodic solution when one delay increases. Our work answers the question in Sect. 5 (Shi and Yu in Chaos Solitons Fractals 100:62–73, 2017).


Introduction
The dynamics of a plankton model is important for an aquatic system: phytoplankton produce oxygen by photosynthesis and absorb nearly half carbon dioxide, so phytoplankton exert great influence on our ecosystem. There have been many works about the dynamics of a phytoplankton-zooplankton model in recent years. As we know, delays always exist in nature, which could bring rich dynamics in a large range of biological systems. So the delayed phytoplankton-zooplankton model has attracted much interest [1][2][3][4][5][6][7][8][9][10]; the model contains two delays, that is, gestation delay for zooplankton and the maturity delay for toxic phytoplankton, which have been studied by many researchers [1,4,6,[8][9][10] in recent years. In [1,9,10], the authors took the gestation delay as a parameter, and the Hopf bifurcation was discussed. In [4,6,8], the author took the delay required for the maturation of toxic phytoplankton as a parameter, and the dynamics of system was studied. We know that there always exist more than one delay in the world, so the dynamics of multiple delays should be more real and interesting. The dynamics of more than one delay have turned up in recent years [1,6,[11][12][13]. For example, the author [6] discussed Hopf bifurcation of the following phytoplankton-zooplankton system with two delays: where the delays τ 1 and τ 2 stand for the time required for maturity of two toxic phytoplankton. Besides we studied the Hopf bifurcation with two delays caused by the gestation of two zooplankton [1]: where P(t), Z i (t) (i = 1, 2) denote the density of phytoplankton and zooplankton respectively, r is the intrinsic growth rate for phytoplankton, K is the environmental carrying capacity of phytoplankton, μ i denotes the maximum uptake rate for zooplankton. P α i +P denotes the Holling functional response which describes the grazing for zooplankton. β i , ρ i are the ratio of biomass conversion and the rate of toxic substance produced by per unit biomass of phytoplankton, d i is the death rate, g i denotes the intrinsic competition coefficients, α i is the half saturation constant for Holling functional response, apparently β i < μ i holds. We know that Wu [14] gave the condition about the global Hopf bifurcation for a delay differential system. There exist much work about the global Hopf bifurcation [5,[7][8][9][10][11][15][16][17]. Especially there are many works on a phytoplankton-zooplankton system [5,[7][8][9][10]. For example, we take the gestation delay as a parameter, and the global Hopf bifurcation of the following system is discussed [10]: (1.3) By now, there have been little work about the global Hopf bifurcation with two delays. In paper [11], the global existence of a bifurcating periodic solution of a neural network model with two delays was studied by limiting τ 1 ∈ (0, τ 10 ) and taking τ 2 as a bifurcation parameter. Based on the above works, we shall study the global Hopf bifurcation of system (1.2) with delays, which is the second question given in Sect. 5 [1]. In this paper we should answer this question: if two delays coexists for system (1.2), limiting τ 1 ∈ (0, τ 10 ), when τ 2 increases, does the periodic solution exist for ever? This paper is organized as follows: in Sect. 2, we give some preliminary work about system (1.2). In Sect. 3, based on [14], we demonstrate the global existence of a periodic solution, then we give some simulation of the periodic solution of system (1.2).

The mathematical model
Before our discussion, we give the condition for the existence and uniqueness of positive equilibrium.
We know when delays do not exist, dz i dt < 0 as β iρ id i < 0. Hence we assume that β iρ id i > 0 holds. From [4], we know that, for system (1.2), a positive equilibrium exists when h 3 Furthermore, the positive equilibrium is unique when h 3 (P) is monotone on (P 0 , P 1 ), the definition of P 0 , P 1 could be found in [4].
From the description of dynamics of system (1.2) under case (5), we could obtain
In [1], we have demonstrated this theorem, so we omit here. Proof: For system (1.2), suppose that there is a nontrivial periodic solution with period τ 2 , then the following system has a nontrivial periodic solution: (3.6) We know that the P -Z 1 -plane and the Z 1 -Z 2 -plane are invariant manifolds for system (3.6), there are no orbits crossing these two planes for the orbits do not intersect each other. Let t = inf{t : Z 1 (t) = 0}, combining with initial condition (3.5) and the second equation of (3.6), we get from which the orbits do not cross the P -Z 2 -plane, so all the orbits must lie in the first quadrant and equilibrium E * should be in the interior of a periodic solution. From Lemma 2.1 (τ 1 ∈ (0, τ 10 ), τ 2 = 0), E * is stable and unique, so it is globally stable in the first quadrant, the periodic solution does not exist. Thus there is no periodic orbit in the first quadrant. onto the T space is also bounded. By Lemma 3.2, we obtain that the connected com-

Conclusion
In this paper, we study the global Hopf bifurcation of two zooplankton-phytoplankton system with two delays. Fixing the delay τ 1 in its stable interval and taking the delay τ 2 as a parameter, we demonstrate that the bifurcating periodic solution exists as τ 2 increases, where the delays τ 1 , τ 2 are due to gestation of two zooplankton. We conclude the global oscillatory nature for the solution of a delayed system when it exhibits oscillatory behavior. Our work answers the question in Sect. 5 [1].