Abstract

Allee effect, as an important biological phenomenon, has been considered in many ecosystems, whereas, its influence on interactions of three or more species is little investigated. In this paper we modify a three-species Leslie-Gower type food chain system by incorporating the strong Allee effect into the source. Our results show that the existence of Allee effect contributes to the occurrence of more complex dynamics of the system, including Hopf, saddle-node, transcritical, saddle-node-Hopf, period-doubling, and period-halving bifurcations and chaos.

1. Introduction

Stable coexistence of species populations has always been investigated as a focus issue in ecosystems. Recently, several studies reported that in the low density populations the stable coexistence could be influenced easily by Allee effect [1]. Under the mechanism, species have underlying danger to undergo a catastrophic collapse [2, 3]. Allee effect, a biological phenomenon of density dependence, has been gradually becoming acknowledged and recognized as a ubiquitous and important mechanism in many theoretical researches of population and community ecology. Therefore, there has been a growing interest in investigating Allee effect [2, 47]. Main attention of some theoretical works related to Allee dynamics is focused on the question that whether interventions of Allee effect to a system have a remarkable influence on dynamics and corresponding evolutions [811].

To the best of our knowledge, most of the known literatures devoted to the Allee dynamics are based on the single species or two-component ecosystem [2, 12] due to the complexity of high dimension models. However, interactions among community ecology usually involve multispecies. Thus, dynamics investigation of multispecies suffering Allee effects is a meaningful supplement from the point of view of both ecology and mathematics, also for the purpose of exploring the stable coexistence strategies of ecosystem. Research related to Allee effect in three species populations firstly is investigated by [13]. Their efforts were focused on the qualitative aspect. In other words, they primarily address the relativity of the emergent Allee effect and the size-specific interactions within the prey population. However, the influences of Allee effect on the property of bifurcations and chaos in the multiple population systems have never been considered in detail. Therefore, in the present paper, we will devote our efforts to quantitative aspects of three-species food chain with Allee effects.

Specifically, Allee effect is characterized by the positive correlation between the population density and the total individual fitness usually quantified by the per capita population growth rate and population size or density [1, 7]. Biological mechanisms giving rise to an Allee effect generally involve social or cooperative processes between multiple individuals [14], for example, mate search difficulty, social facilitation of reproduction, and increased predation risk due to schooling behavior [1, 2, 7, 14, 15]. There is also another view on the cause of Allee effect. In [14, 16], investigators considered Allee effect as emergent consequences of stage-specific predation. They pointed out that Allee effect came from the shift which a predator induces in the population stage distribution of its prey [14]. Generally, when a single species is subject to the Allee effects, its dynamics is governed by the following system:where is population density at time , satisfying . presents the per capita growth rate of , and there is a range of in the interval , such thatholds. Generally, (1) has the trivial equilibrium and maximum biological equilibria and also admits, at most, one other biological equilibrium between and . Following the definition given by [17], system (1) is identified to have weak Allee effects provided that (2) and hold. Furthermore, a strong Allee effect can occur in system (1) if (2) and hold.

Some theoretical models obeying the definition of (1) are presented and studied solely [17, 18] or integrated into the population community [8, 1921]. Through above works one can see that the difference between the weak Allee effect and the strong Allee effect is the population growth rate, and the former rises with increasing density; however, the latter does not, which commonly exhibits the alternative stable states or the catastrophic shifts in structure and dynamics corresponding to the bistable or the monostable stages, respectively. Particularly, when system is at the bistable stage, the strong Allee effect introduces a critical population density (or size) above which the particular trophic level of the population will always be present, whereas below which the population will decrease to a low trophic level and even go extinct because the growth rate of the species might become negative [1, 3, 14, 22]. When the species suffers a strong Allee effect, some researches indicated that it is apt to have less ability to overcome these additional negative competition effects [23] and thus to be prone to extinction than other species [2, 24]. Obviously, failing to take the strong Allee effect into account could have detrimental consequences for the management and conversation of some natural and exploited populations. This is the exact reason why most works related to Allee dynamics emphasized the investigation of the strong Allee effect [68, 20, 25]. These motivate us to study the dynamics induced by strong Allee effect in a multiple species system.

The remaining part of this paper is organized as follows. In the next section, we derive a food chain model with demographic Allee effect [7] caused by considering mate search difficulty of the prey. In Section 3, the existence and the local stability of equilibria of system (5) are discussed. In Section 4, we investigate the bifurcation patterns that system (7) can undergo and chaos behaviors by using the Routh-Hurwitz rule and some numerical simulations. Finally, a simple discussion is given to conclude this paper.

2. Model

In recent times, several two-component prey predator systems in a homogeneous environment with the Allee effects and Leslie-Gower term have been investigated [2, 8, 20, 26], which take the form aswhere obeys the definition of (1) and is a logistic growth rule with an Allee-type function. is the predator functional response to prey. is called the Leslie-Gower scheme; it is assumed that the growth of the predator with the per capita maximum fertility rate also follows the logistic type but its environmental carrying capacity depends on the available prey quantity. Such mechanism is also adopted widely in ecosystem modelling [2732].

In the low density populations, the Allee-type function in (3) can be considered the demographic Allee effect; the cause of that is owing to a reduced ability of females to find mates, which leads to the species growth at a lower rate; for instance, the growth of Glanville fritillary butterfly is reduced due to the difficulty of mating limitation [33]. Mathematically, the biological phenomenon can be displayed bywhere is used to quantify the strength of Allee effect; and are the per capita maximum fertility and death rates of the source. They satisfy ; represents the intracompetition intensity of the source population [19]. The growth of the species is assumed to obey the logistic rule but imposed a Michaelis-Menten type function . Biologically, the Allee effect presents a density threshold of population at which the species reaches half of its maximum fertility. By the definitions in (1) and (2), it is easy to infer that system (4) has the strong Allee effect, but not the weak Allee effect. Such Allee effects firstly were presented and used in predator-prey system by [19], who investigate the influence of Allee effect on the ecological and evolutionary conditions that allow for continuously stable strategy and evolutionary branching.

To study Allee dynamics in multiple species, we develop (3), a Leslie-Gower food chain system, by introducing a top predator which consumes the intermediate predator to maintain growth itself. The detailed form of the developed systems is described by the following differential equations:where denote the source, the consumer, and the predator, respectively. Parameters , , , , , , , , , , and are positive, where is a quantification of the food supplied by the source; and are the capture rate and conversion efficiency.

A similar system as (5) but without suffering Allee effect has been investigated by [34]. They showed a path along with which system translates its dynamics from Hopf bifurcation to chaos. Compared with the models of [19, 34], system (5) introduces more biological mechanisms; thus, it is more natural to model the interactions of multiple populations, especially for the case of the source at low density. One also can see that system (5) may have more biological steady states due to the existence of the Allee effect. Thus, it is likely to exhibit more complex dynamics which is unclear and needs to be investigated further. On another side, it is different from the works of [19, 34], and we will devote ourselves to address the dynamics analysis of system (5) caused by the Allee mechanism. Our results show that the Allee effect introduced by the low fitness of the source individual facilitates the occurrences of complex bifurcation and chaotic behaviors.

3. Existence and Stability of Equilibria

In this section, we focus on the local dynamical behaviors of system (5). In view of biological meaning, we consider system (5) only in or . Clearly, with , the solutions of system (5) are positive for all .

For simplicity, we begin by taking the scalingto (5), and getting rid of all the bars then the following nondimensional system is obtained:where

Then the biological equilibria of system (7) are obtained as follows.

Lemma 1. Assume . In , system (7) may have four boundary equilibria denoted by , , , and if and also may admit two positive and if , , , and , where

Remark 2. From Lemma 1, we know that system (7) displays multiple equilibria owing to the existence of the strong Allee effect in the source, which is different from the results of [34] (three nonnegative equilibria). In the following analysis, one can see multiple steady states contributing to the occurrence of saddle-node bifurcation. Thus, the system is prone to undergo the catastrophic collapses.

The Jacobian matrix of system (7) at any nonnegative equilibrium takes the form as follows:where .

By recurring to estimate the eigenvalues of these equilibria, we have the following results.

Lemma 3. Provided , , then , , and are saddles. However, may be a saddle or node or saddle-focus or focus-node for . More precisely, we get the following:(i) is saddle or saddle-focus when , , and is stable node or focus-node when , .(ii) is saddle or unstable focus-node when , , and is saddle or saddle-focus when , .

, , and and detailed proof of Lemma 3 can be seen in the Appendix.

Lemma 4. Assume , , and ; then is unstable, and is a stable node or stable node-focus for .

Proof. The Jacobian matrix of system (7) at equilibrium iswith From Lemma 1, we know that if exists, then . Computation shows that is a monotonously decreasing function of , in which, together with , we have for all . Hence, the determinant of denoted by is negative, which implies that is a saddle or a saddle-focus or an unstable node.
Similarly, the Jacobian matrix of system (7) at equilibrium iswith The characteristic equation of iswhere In order to use the Routh-Hurwitz criterion, we get By computation, we infer that for and forTherefore, , , , and if and (18) hold, which indicates that is stable node or stable node-focus.

4. Results

4.1. Hopf Bifurcation Analysis

From Lemma 4 and theory developed by Liu [35], we derive that system (7) can undergo the Hopf bifurcation at if there exists and satisfies that

Firstly, we notice that if holds. Therefore, if exists, then . can be rewritten as a quadratic function of , andA simple computation shows that for ; thus, . That is, (21) implies that exists and satisfies , , , and . Also notice that , such that .

Secondly, the transverse condition Without losing the generality, we set ; thus, together with the above analysis and Lemma 4, we have the following result.

Theorem 5. Set , , and ; then there is a , so that equilibrium is local asymptotically stable when and unstable when . Furthermore, a Hopf bifurcation occurs when .

Remark 6. In Theorem 5, we choose as bifurcation parameter. When approaches , a stable periodic solution is bifurcated from equilibrium . From Figure 1(a), we can see that the value of the source , except for its stability, does not vary with changing continuously. However, if we take the Allee effect as bifurcation parameter, in Figure 1(b), we discover that the value and stability of will change with varying. These indicate that the Allee effect has the bigger impacts than on the dynamics of system (7).

4.2. Dynamics Evolution

To completely know the dynamics induced by the combined action of the Allee effect and other parameters, in the following, we investigate the dynamical behavior evolution of system (7) by recurring to the numerical methods.

Firstly, the Hopf bifurcation as shown in Figure 1(a) occurs at , which forces three species to shift their states from the stable coexistence to periodical cycle as decreasing continuously. However, the period-1 behavior will be translated to period-2 cycle with a period-doubling bifurcation appearing at . Subsequently, the period-doubling cascade happens with decreasing which may lead to the chaos behavior of system (7); for instance, when , system (7) exhibits a chaos attractor (see Figure 2(a)).

Similar result of Hopf bifurcation also is achieved by the Allee effect . As shown in Figure 1(b), we see that system (7) undergoes the Hopf bifurcation when . In addition, from Lemma 4, we also see that under the condition , , and two equilibria and coexist, one is stable and the other is unstable. Furthermore, and collide when , which signifies that a saddle-node bifurcation occurs (see Figure 1(b)). Similar arguments can be implemented to equilibria and ; there also is a saddle-node bifurcation at .

System parameters are chosen as , , , , , , , , and such that system (7) has a positive equilibrium as the saddle-focus with a stable one-dimensional manifold and an unstable two-dimensional manifold and a boundary equilibrium as the saddle; then we can see that there is a stable limit cycle bifurcating from the homoclinic bifurcation as shown in Figure 3.

Finally, from Lemmas 3 and 4, we can exhibit the existence of transcritical bifurcation associated with equilibria and . Figure 4(a) shows that and get close to each other with increasing and coincide and exchange their stabilities at and subsequently separate into two different steady states.

Remark 7. From the above numerical examples, one can see that system (7) may display Hopf, saddle-node, homoclinic, and transcritical bifurcations with parameters and changing. Moreover, system (7) also can undergo other degenerated bifurcations, besides the saddle-node bifurcation, for example, saddle-node-Hopf bifurcation. It is easy to check that the characteristic polynomial corresponding to the Jacobian matrix (13) has a zero eigenvalue when . If we take then matrix (13) also has a pair of purely imaginary eigenvalues. These imply that system (7) may result in the saddle-node-Hopf bifurcation which is verified in Figure 4(b) (the intersection of the Hopf bifurcation curve and the saddle-node curve).

4.3. Chaotic Behavior

It is known from the published theoretical works that ecological chaos may have prominent influences on the population and community dynamics. Thus, how to identify the factors potentially facilitating or inhibiting chaotic behavior is an important but challenging subject. In what follows, extensive numerical simulations are carried out for various values of parameters and for different sets of initial conditions. We choose , , and as bifurcation parameters, respectively, and keep other parameters fixed.

Firstly, to get a deep insight into the complex dynamics induced by the Allee effect, we implement the numerical simulation in plane (). Figure 5 exhibits a route along with which system (7) transforms its dynamical behaviors from period- attractors () to chaos motion and to period- attractor through some period-doubling and period-halving sequences, respectively. More precisely, with increasing, a period-doubling bifurcation occurs at , which forces system (7) to go to a period-4 cycle. Subsequent period-doubling cascade brings the system to the chaotic motion with increasing continuously.

On the contrary, the chaotic motion can inversely switch to period- state through several period-halving bifurcations. It is known that these chaotic behaviors displayed in system (7) heavily depend on the initial states of three species or the perturbation coming from the external stimulations. That is, the different dynamical regimes are likely to be exhibited at the same parameter value. For instance, choosing , system (7) may undergo a period-2 cycle when or a quasiperiodic behavior when (see Figure 6) which is the most intuitive definition of chaos.

To deeply understand the combined effect of the Allee effect and other parameters, similar effects are produced by an increase in parameters and , respectively (see Figures 7 and 9).

5. Conclusions and Discussion

The complexities of ecosystems, such as multiple bifurcations and chaotic behaviors, have been a focus of intense research in the ecologic and mathematics fields for the purpose of investigating how to keep the persistent coexistence of species populations in community [30, 36]. Such studies of population interaction based on mathematical models have provided valuable insight into ecological dynamics. Recently, several works were devoted to address the influence of Allee effect on dynamics of population to explore the underlying mechanism of regulating stable coexistence of species in community.

In the present paper, we consider the Allee effect for the prey (the resource) in a modified Leslie-Gower food chain system. Some main properties including Hopf, period-doubling, saddle-node, transcritical, and saddle-node-Hopf bifurcations and chaotic behaviors are studied qualitatively and quantitatively. Results based on our analysis and numerical simulations show that the Allee effect not only increases the complexity of the system structure but also promotes the occurrence of chaotic behaviors. As shown in Figures 7 and 8, system displays period-doubling cascade, chaos motion, and period-halving cascade when interfered by the Allee effect (see Figure 7), whereas it only keeps period-2 and period-4 cycles when there is no disturbance of the Allee effect (see Figure 8).

Appendix

The Proof of Lemma 3

Proof. Substituting and into (10), respectively, one can show that both and have eigenvalues and ; thus, and are saddles.
Jacobian matrix at iswhere .
Obviously, one eigenvalue of is ; others are determined by the matrixComputation shows that the determinant of (A.2) satisfies with the conditions and , which implies that matrix (A.2) has two real eigenvalues but has opposite signs. That is, is a saddle.
Similarly we also carry out investigation for the eigenvalues of , one is , and others may be two real roots with the same signs or be a pair of conjugate complex numbers since the determinant of satisfies that . The trace of denoted by can be regarded as a continuous function of variable for . Hence, ifthen there exists at least one so that . Equation (A.3) is equivalent toNotice that for ; hence (A.4) can be rewritten asDefine . Without loss of generality, set and . It follows from (A.5) that is unstable if and stable if and . To sum up the above analysis, we can attain the results of Lemma 3.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors would like to thank the editor and the anonymous reviewers for their constructive suggestions and comments, which improved the presentation of the paper.