#### Abstract

Complex dynamics of a four-species food web with two preys, one middle predator, and one top predator are investigated. Via the method of Jacobian matrix, the stability of coexisting equilibrium for all populations is determined. Based on this equilibrium, three bifurcations, i.e., Hopf bifurcation, Hopf-Hopf bifurcation, and period-doubling bifurcation, are analyzed by center manifold theorem, bifurcation theorem, and numerical simulations. We reveal that, influenced by the three bifurcations, the food web can exhibit very complex dynamical behaviors, including limit cycles, quasiperiodic behaviors, chaotic attractors, route to chaos, period-doubling cascade in orbits of period 2, 4, and 8 and period 3, 6, and 12, periodic windows, intermittent period, and chaos crisis. However, the complex dynamics may disappear with the extinction of one of the four populations, which may also lead to collapse of the food web. It suggests that the dynamical complexity and food web stability are determined by the food web structure and existing populations.

#### 1. Introduction

The concept and research method of the food web was firstly proposed as early as in 1930s by Elton, who initiated exploring community structure in ecosystems from the viewpoint of network [1]. After that, the food webs have been widely studied in various aspects of characteristics, involving stability, complexity, dynamics, etc. Via analyzing the characteristics of food webs, many research works have played an important role in understanding the properties of natural ecosystems [2–4].

Food web systems generally are composed of different combinations of relationships between interacting species, among which predation relationship is a basic one [5, 6]. Understanding the interactions between prey and predator is a significant topic of longstanding interest in biology. The prey-predator models with various functional responses, such as Holling I, II, III, and IV, Leslie-Gower, Michaelis-Menten, ratio-dependent, Beddington-DeAngelis, and so on, have been extensively discussed [7–10]. Many other ecological mechanisms have been also included and investigated in the prey-predator models. For example, Raw et al. proposed and analyzed a tri-trophic prey-predator model with two predators and one prey exhibiting group defense mechanism and found occurrence of complex dynamical behaviors [7].

Over past decades, mathematical models have been investigated to study the dynamics of three or more interacting species in food webs. Early studies mainly focused on three-species food web models, demonstrating the complexity of the dynamics for the food webs [6, 11]. The food web systems may converge not only to stable equilibria and periodic orbits, but also to more complex behaviors such as quasiperiodic oscillations and chaos induced by different bifurcations [11]. As studied by Gakkhar and Gupta [6], generalized Hopf bifurcation can exhibit in a food web system with commensalism and Bogdanov-Takens bifurcation point is detected only for a certain value of commensal coefficient. Gupta and Chandra [11] investigated a prey-predator-scavenger model with quadratic harvesting and found that all three species have unpredictable and irregular behavior which may be controlled with the help of an appropriate number of harvesting so that the system can be dynamically balanced.

Currently, research on food web dynamics slipped into the cases where four species or more species exist and interact [2, 12, 13]. Rich dynamic behaviors of four-species or more species food webs, such as equilibria, bifurcations, periodic orbits, torus, chaotic attractors, and so on, have been found to be contributing to better comprehending on the complexity of species coexistence and community stability. In Gakkhar et al. [2], a four-species food web model was developed and complex dynamic behaviors including strange attractor and chaos were shown. In the food web studied by Wei [13], Wei proved that the food web system possesses a line segment of degenerate steady states for the parameter values on a bifurcation line in the bifurcation diagram and found diverse routes to chaos.

Investigation has demonstrated that the topological structure of food web can play an important role in affecting the food web dynamics. Generally, with the variation of variable diversity and trophic links, the structure of food webs may be different. In Fussmann and Gerd [14], the architecture of a food web with three primary producers, two consumers, one primary predator, one secondary predator, and one tertiary predator was emphasized. Many basic food webs that differ in complexity and in the number of trophic levels and omnivorous links can be generated as substructures of this web. According to the difference of network structure, Kuijper et al. [15] classified the four-species food webs into three types, i.e., competition, omnivory, and food-chain, and found that the architecture of the food web depends on the feeding behaviors of the highest trophic level. Based on the research works of Fussmann and Gerd [14] and Kuijper et al. [15], this research focuses on tri-trophic level food web which includes two-preys, one-middle predator feeding on the two-preys, and one-top predator feeding on all three other populations. The study on the bifurcations and complex dynamics of this food web is still not documented in literature available to the authors.

In this research, we modeled the above food web on the basis of the Hasting-Powell [16] model with application of Holling II functional response, and the exploration is arranged as follows. In Section 2, we introduce the four-species food web model and analyze local stability of the coexisting equilibrium. In Section 3, the bifurcation theorem and center manifold theorem are used to determine the conditions for Hopf bifurcation of the food web. Section 4 investigates the Hopf-Hopf bifurcation and performs simulations to demonstrate the quasiperiodic dynamics. Section 5 numerically studies the period-doubling bifurcation and demonstrates the bifurcation diagram in which a route to chaos emerges. In Section 6, discussion and conclusions are provided.

#### 2. Mathematical Model and Stability Analysis

In 2002, Fussmann and Gerd [14] have studied the food web complexity and chaotic population dynamics, with proposition of architecture of food web models with two, three, and four to five trophic levels. Their results suggested that natural food webs possess architectural properties that may inherently lower the maybe of chaotic community dynamics [14]. Based on their research work, we choose to study the architecture of food web models with three trophic levels including two preys, a middle predator, and a top predator. With the application of Hasting-Powell [16] model and Holling II functional response, a four-species food web model can be described by a dynamical system of four ordinary differential equations: withIn the above equations, , , , and are population densities of the four interacting species and is time; and represent the intrinsic growth rates of species and ; and are the carrying capacities; represents the conversion rate of prey to predator; and are the death rates of species and , respectively. represents a functional response of the Holling-II type; and are parameters used in the functional response: describes the effects of capture rate and handling time; determines how fast per capita feeding rate approaches its saturation rate. From the viewpoint of ecology, parameters , , , , , , and should all be positive constants.

It should also be noticed that population densities are nonnegative. Hence, the state space of the food web system (1) must be strictly restricted in the region . In order to reduce the number of parameters in (1), and to determine which combinations of parameters control the behaviors of the system, system (1) is changed with the application of the following expressions: Using (3) to replace the corresponding variables in (1), then we obtain a simpler form of the food web system: in which

The equilibriums of system (4) can be obtained by soling equations , , , and . Direct calculation yields that system (4) has following six kinds of equilibriums:

(i) Trivial equilibrium always exists, describing the absence of populations.

(ii) Axial equilibriums and indicate the prey population or increases to carrying capacity in the absence of predation. Likewise, equilibrium shows the case where prey populations and increase to the carrying capacity simultaneously in the absence of predation.

(iii) When one predator and one prey are absent simultaneously, system (4) degenerates to be a normal Holling II predator-prey system. In such cases, we can find the equilibriums for coexistence of predator and prey, including , , , and .

(iv) With the consideration of reducing one predator population, system (4) degenerates to be a two-trophic-level food web system with two preys and one predator. The equilibriums of such type of food web have been studied by Misra et al. [17].

(v) The equilibriums for the case with absence of one prey population have been investigated by Hsu et al. [18]. One can refer to that literature for detailed calculations of equilibriums.

(vi) This research focuses on the positive interior equilibrium , which represents the coexistence of the four species in the food web. From system (4), it is easy to know that is a positive root of the following nonlinear algebraic equations:

To determine the local asymptotical stability of the equilibrium , the method of Jacobian matrix is utilized. The Jacobian matrix associated to system (4) at any point can be described as follows: The characteristic equation of the Jacobian matrix at can be written aswhereand is the element of Jacobian matrix at the th row and th volume.

DefineAccording to the Routh-Hurwitz criterion, when the following conditionsare satisfied simultaneously, the equilibrium is locally and asymptotically stable. Otherwise, this equilibrium becomes unstable.

A numerical example is provided here to illustrate the emergence of stable coexistent equilibrium of the food web system. We set parametric conditions of system (4) as follows:

Figure 1 shows the dynamics of the food web when the parameter values are given in (12). In the graph of time series (Figure 1(a)), variables , , , oscillate at the beginning with decreasing amplitude and become steady in long-term behaviors. In the and phase diagrams (Figures 1(b) and 1(c)), it is indicated that the coexistent equilibrium is a stable focus. It can be verified that From determinative condition (11), we also know that the coexistent equilibrium is locally and asymptotically stable.

**(a)**

**(b)**

**(c)**

#### 3. Hopf Bifurcation

In this section, Hopf bifurcation theorem and center manifold theorem [19] are used to determine the occurrence conditions of Hopf bifurcation for the four-species food web model. To study the Hopf bifurcation, bifurcation parameter should be chosen at first. Among all parameters of system (4), parameter is basic and represents the growth rate of prey , controlling one of energy inputs of the food web and therefore we choose it as the bifurcation parameter.

According to the Hopf bifurcation theorem, when system (4) undergoes Hopf bifurcation, the characteristic equation should have a pair of conjugate pure imaginary roots, i.e.,Substituting the eigenvalues (14) into the characteristic equation (8), and separating the real and imaginary parts, then we haveReducing the unknown quantity , the above equations (15) and (16) turn to be Equation (17) determines the relationship of the parameters which exactly leads to the emergence of Hopf bifurcation. With the other parameter values given, we can calculate the threshold value of the bifurcation parameter . In the following, we denote this threshold value of Hopf bifurcation point as .

When the value of parameter changes around the Hopf bifurcation point, the two conjugate imaginary eigenvalues will change to be complex eigenvalues, written as . Substituting into (8), then we obtain the following [12]:

Through (19), we find can be expressed by . Replacing such expression of into (18) results inwhere . Differentiating with respect to and putting , we getIf system (4) undergoes Hopf bifurcation around the interior equilibrium , we should have .

System (4) undergoing Hopf bifurcation at should be also under the satisfaction of the last determinative value. We consider two conjugate imaginary eigenvalues and two conjugate complex eigenvalues for the characteristic equation (8). To calculate the determinative value, we further transform the equations of system (4) by applying the following translation:Substituting (22) into the equations of system (4) and separating the linear part and nonlinear part. Then system (4) can be rewritten aswherein which is a fourth-order polynomial functions in regard to independent variable , and

Via introducing an invertible transformation, i.e., wherewhere the elements of the matrix are provided in Appendix, then the normal form of system (23) can be given bywherein which , , , and are coefficients of the polynomial function .

Next, the center manifold of system (28) at the origin is calculated. On the basis of center manifold theorem, a center manifold exists and can be described as follows: where and are assumed as By computation, we get

Accordingly, the dynamics restricted to the center manifold can now be considered. Combining with (31), system (28) is simplified aswherein which

We introduce the following partial derivative notation: in which subscript demonstrates the partial derivative of the first variable; subscript is the partial derivative of the second variable. On account of (34), we have , , , and . Moreover, system (33) has two eigenvalues as . Consequently, the final determination for the system (4) experiencing Hopf bifurcation is calculated by the following [20]: Direct calculation yields

Summarizing the above analysis, if condition (38) is satisfied, then system (4) undergoes Hopf bifurcation at equilibrium . Furthermore, when establishes, if and , the system experiences supercritical bifurcation, and a stable orbit generates with radius as about ; if and , subcritical Hopf bifurcation occurs and results in the emergence of an unstable orbit with radius as about .

For illustrating the occurrence of Hopf bifurcation in the food web system, an example is provided here. First, parameter is still set as the bifurcation parameter and the other parameter values are given in (12). Under such parametric conditions, the bifurcation critical value of for occurrence of Hopf bifurcation in system (4) is determined as . Simultaneously, the coexistent equilibrium at the bifurcation critical point is and the corresponding eigenvalues are , , and , which show two pairs of imaginary numbers. Moreover, we can derive and the determinative value is calculated as . Consequently, we know that system (4) experiences supercritical Hopf bifurcation at the coexistence equilibrium , i.e., when , the dynamics of the food web converg to the coexistent equilibrium; when , the coexistent equilibrium is weakly attracting; when , the dynamics of the food web are attracted to a limit cycle (Figures 2 and 3).

**(a)**

**(b)**

**(c)**

**(a)**

**(b)**

Figure 2 is plotted to show the dynamics of the food web when . By Routh-Hurwitz criterion, we obtain suggesting that the coexistent equilibrium is unstable. In this case, , , , and will keep periodically oscillating with long-term behavior (Figure 2(a)). In the phase diagram, the periodically oscillating behavior is clearly manifested as a limit cycle (Figures 2(b) and 2(c)).

Figure 3 further demonstrates the Hopf bifurcation diagrams with the change of parameter . As the value of parameter increases, the system dynamics vary from stable equilibrium to limit cycle, forming bowl-shaped surfaces in Figure 3, which vividly and clearly demonstrate the occurrence of the Hopf bifurcation in the food web system.

#### 4. Hopf-Hopf Bifurcation

Beyond the Hopf bifurcation, the food web system can still further undergo Hopf-Hopf bifurcation with parameter variation. Hopf-Hopf bifurcation is a complicated codimension-two bifurcation, which only appears in four-species food web system or higher-species food web systems. Hopf-Hopf bifurcation can lead to the emergence of a torus solution, appearing like tyre. Hopf-Hopf bifurcation in physical engineering field was fully studied [19, 21, 22]. However, the research on Hopf-Hopf bifurcations of food web system is still not documented in the literature available to the authors.

For studying the Hopf-Hopf bifurcation, we choose two parameters and as the bifurcation parameters. Depending on the variation of bifurcation parameters, the food web system may undergo Hopf-Hopf bifurcation. At the bifurcation point, equilibrium has two pairs of different purely imaginary eigenvalues, denoted aswhere .

Then (8) can be rewritten asBalancing the right side of both (8) and (41), we easily obtain

From the conditions and , we can get the threshold values of the two bifurcation parameters as and . Moreover, if the conditions , , and are satisfied, the values of and can be also obtained, i.e., In this critical case, the food web system undergoes Hopf-Hopf bifurcation at the equilibrium . To demonstrate the Hopf-Hopf bifurcation numerically, we set parameter values of the food web system as follows:

Under the parametric conditions of (44), the threshold values of bifurcation parameters at the Hopf-Hopf bifurcation point are determined as , . Meanwhile, the coexisting equilibrium is obtained as and its eigenvalues are and , which verify the occurrence of the Hopf-Hopf bifurcation in the four-species food web system. As shown in Figure 4, the dynamics of the system converge to a stable invariant torus in such case, which represents quasiperiodic dynamics of the food web system.

**(a)**

**(b)**

**(c)**

**(d)**

Figure 5 demonstrates the dynamic transition around the Hopf-Hopf bifurcation point with the change of parameter . As the value of parameter increases, the system dynamics vary from limit cycle to torus (Figure 5(a)). In the Poincaré map, this phenomenon reflects as a dynamic variation from a stable fixed point to a limit cycle, as shown in Figure 5(b).

**(a)**

**(b)**

Figure 6 further exhibits a bifurcation diagram when the food web system experiences Hopf-Hopf bifurcation. The bifurcation diagram is plotted with the successive maximum of as the parameter varies in the range . As shown in Figure 6, with the value increasing around the critical value , the system dynamics experience transition from a stable limit cycle to a torus. After that the food web system comes into the zone of quasiperiodic dynamics, in which a few periodic windows can still take place, such as periodic windows with period 7, 8 orbits (Figures 7(b) and 7(d)). These periodic dynamics are intermittent and appear and disappear abruptly, resulting in alteration of periodic and quasiperiodic dynamics with the increase of bifurcation parameter. At the end of the quasiperiodic zone, the system dynamics turn to be periodic, as demonstrated by period 4 orbit in Figure 7(f).

**(a)**= 1.12

**(b)**= 1.148

**(c)**= 1.15

**(d)**= 1.175

**(e)**= 1.23

**(f)**= 1.27#### 5. Period-Doubling Bifurcation

Based on the Hopf bifurcation, we find that the food web can further undergo period-doubling bifurcation. Under the influence of the period-doubling bifurcation, the food web system may experience period-doubling cascade, eventually leading to the occurrence of complex periodic behaviors and chaotic dynamics.

In order to numerically display the period-doubling bifurcation and corresponding system dynamics, we give the parametric conditions as follows: Here, we also choose parameter as the bifurcation parameter, and the initial condition is set as (0.5, 0.5, 0.5, 0.5). The numerical simulation results are shown in Figures 7–11. With the change of parameter , the dynamics of all species exhibit period-doubling cascade, leading to dynamical transition from period to chaos. When , the system dynamics exhibit simple periodic oscillation, which appears as a limit cycle in phase portrait (Figure 8(a)) and a single point in Poincaré map (Figure 9(a)). When the value of parameter increases to , the period of the periodic oscillation doubles, i.e., period 2 orbit. This reflects as two cycles in Figure 8(b) and two points in the Poincaré map of Figure 9(b). Consequently with the increase of value, we find periodic orbits with period 4, period 8, and period 16 for , , and , respectively. In the phase portrait and Poincaré map, these periodic dynamical behaviors represent as 4, 8, and 16 cycles and points (see Figures 8(c)-8(e) and 9(c)-9(e)). As a result of the period-doubling cascade, chaotic attractors are discovered in long-term behaviors when , as shown in Figures 8(f) and 9(f).

**(a)**= 1.15

**(b)**= 1.25

**(c)**= 1.3425

**(d)**= 1.3452

**(e)**= 1.3495

**(f)**= 1.45

**(a)**= 1.15

**(b)**= 1.25

**(c)**= 1.3425

**(d)**= 1.3452

**(e)**= 1.3495

**(f)**= 1.45**(a)**

**(b)**

**(a)**

**(b)**

**(c)**

To further determine the chaotic behaviors of the food web, sensitivity analysis is performed via applying two close initial conditions and then observing the system dynamics in long-term behaviors. In Figure 10, red and blue waves show the dynamical behaviors with reference to two close initial conditions, namely,* x*_{0} = 0.5,* y*_{0} = 0.5,* z*_{0} = 0.5, = 0.5 and* x*_{0} = 0. 5001,* y*_{0} = 0.5,* z*_{0} = 0.5, = 0.5. We find that a tiny change (0.0001 in magnitude of* x*_{0}) in initial conditions can result in a large difference of* z* and dynamics as the time progresses. It indicates that the solutions under such parametric conditions are highly sensitive to initial conditions, suggesting the emergence of chaotic dynamics for the food web.

To explicitly demonstrate the dynamic transition between nonchaos and chaos with the increase of bifurcation parameter, bifurcation diagrams are plotted, as shown in Figure 11. Firstly, we find that the period-doubling bifurcation starts a route to chaos (Figure 11(a)). Secondly, when the food web goes into the chaotic zone, periodic windows can repeatedly take place. In the periodic windows, period-doubling cascade also emerges and leads to multiperiodic behaviors (Figure 11(b)). Figure 12 displays the system dynamics in a periodic window with period 3 orbit and period-doubling cascade in orbits of period 3, 6, and 12. Thirdly, when the periodic windows end, the food web may experience chaos crisis, which leads to an abrupt transition from chaotic subsets to a chaotic band (Figure 11(c)). Lastly, the chaos disappears when parameter increases to higher values (such as ), and the food web dynamics return to be periodic.

**(a)**= 1.418

**(b)**= 1.424

**(c)**= 1.4242#### 6. Discussion and Conclusion

In this research, we have studied the dynamics of a four-species food web with Holling II functional response. Three bifurcations, Hopf bifurcation, Hopf-Hopf bifurcation, and period-doubling bifurcation, are found. Various dynamical behaviors are observed, including stable equilibrium, periodic oscillations, quasiperiodic oscillations, period-doubling cascade, periodic window, chaotic interior crisis, and chaotic attractors. With the variation of bifurcation parameter, the food web system may even experience the transitions between equilibrium and period, period and quasiperiod, or nonchaos and chaos.

The complex dynamics occurring in the food web result from the interactions among all populations. If one of the populations is extinct, possibly due to catastrophic shifts of environmental conditions, the system dynamics may experience dramatic change. In order to exhibit the contribution of populations on the dynamics of the food web system, we here choose three cases of typical dynamical behaviors, period (Figure 2), quasiperiod (Figure 4), and chaos (Figure 8(f)) and then remove one population and observe the dynamics change. Figure 13 shows the situation where the prey population is extinct. We find that with the absence of one prey, the middle predator or top predator is going to be extinct in the three cases of Figure 13 and that the food web system degenerates a predator-prey system. Similar phenomenon also occurs in the situation when prey or predator is extinct. As shown in Figure 14, it is discovered that, with the absence of predator , the prey population always tends to extinction and that populations and keep periodically oscillating in long-term behavior for all three cases. This result demonstrates that the extinction of one population can lead to instability and collapse of the food web.

**(a)**

**(b)**

**(c)**

**(a)**

**(b)**

**(c)**

When we consider that the middle predator is extinct, the periodic, quasiperiodic, and chaotic behaviors of the food web system all change to be a stable equilibrium (Figure 15). Although the left three populations are still coexisting, the dynamical complexity of the food web is reduced in these cases. The results in Figures 14 and 15 suggest that the omnivorous population can balance the population abundance and enhance the stability of the food web system. Moreover, Figure 15 also implies that the food web has unpredictable and irregular dynamic behaviors which may be controlled with the help of appropriate number of predation so that the food web can be dynamical equilibrium.

**(a)**

**(b)**

**(c)**

#### Appendix

The expression of elements is shown as follows: