Degenerate Bogdanov-Takens Research Articles

Bifurcation of a Leslie–Gower Predator–Prey Model with Nonlinear Harvesting and a Generalist Predator

A Leslie–Gower predator–prey model with nonlinear harvesting and a generalist predator is considered in this paper. It is shown that the degenerate positive equilibrium of the system is a cusp of codimension up to 4, and the system admits the cusp-type degenerate Bogdanov–Takens bifurcation of codimension 4. Moreover, the system has a weak focus of at least order 3 and can undergo degenerate Hopf bifurcation of codimension 3. We verify, through numerical simulations, that the system admits three different stable states, such as a stable fixed point and three limit cycles (the middle one is unstable), or two stable fixed points and two limit cycles. Our results reveal that nonlinear harvesting and a generalist predator can lead to richer dynamics and bifurcations (such as three limit cycles or tristability); specifically, harvesting can cause the extinction of prey, but a generalist predator provides some protection for the predator in the absence of prey.

Dynamics of the generalized Rosenzweig–MacArthur model in a changing and patchy environment

In this paper, we propose a generalized Rosenzweig–MacArthur predator–prey model with two patches and Allee effects in predators. In a constant environment, for the single-patch model we show the existence of degenerate Bogdanov–Takens bifurcation with codimension up to 4 and degenerate Hopf bifurcation with codimension up to 3, which indicate that the nilpotent cusp of codimension 4 is the organizing center of the bifurcation set. Our results show that Allee effects in predators can not only induce much more complex bifurcations and dynamics (such as three limit cycles, multitype tristability), but also serve as a destabilizing factor to make predators extinction (deterministically) when the environment is sufficiently enriched, which provides strong support for the so-called “paradox of enrichment”. For the two-patch model, we numerically reveal the existence of eight positive steady states including two new coexisting attractors, which means that population dispersal can help species coexist or increase the population abundance. In a changing environment, for the single-patch model we find that the population can track unstable states, fast positive or slow negative environmental change can be beneficial to predators by delaying or even avoiding extinction; while for the two-patch model we observe some interesting phenomena: fast negative (or slow positive) environmental change can make predators in different patches trend to extinction synchronously, and rate-induced different regime shifts.

Bifurcation Analysis of a Holling–Tanner Model with Generalist Predator and Constant-Yield Harvesting

In this paper, we introduce constant-yield prey harvesting into the Holling–Tanner model with generalist predator. We prove that the unique positive equilibrium is a cusp of codimension 4. As the parameter values change, the system exhibits degenerate Bogdanov–Takens bifurcation of codimension 4. Using the resultant elimination method, we show that the positive equilibrium is a weak focus of order 2, and the system undergoes degenerate Hopf bifurcation of codimension 2 and has two limit cycles. By numerical simulations, we demonstrate that the system exhibits homoclinic bifurcation and saddle–node bifurcation of limit cycles as the parameters are varied. The main results show that constant-yield prey harvesting and generalist predator can lead to complex dynamic behavior of the model.

Effects of whaling and krill fishing on the whale-krill predation dynamics: bifurcations in a harvested predator-prey model with Holling type I functional response.

In the Antarctic, the whale population had been reduced dramatically due to the unregulated whaling. It was expected that Antarctic krill, the main prey of whales, would grow significantly as a consequence and exploratory krill fishing was practiced in some areas. However, it was found that there has been a substantial decline in abundance of krill since the end of whaling, which is the phenomenon of krill paradox. In this paper, to study the krill-whale interaction we revisit a harvested predator-prey model with Holling I functional response. We find that the model admits at most two positive equilibria. When the two positive equilibria are located in the region , the model exhibits degenerate Bogdanov-Takens bifurcation with codimension up to 3 and Hopf bifurcation with codimension up to 2 by rigorous bifurcation analysis. When the two positive equilibria are located in the region , the model has no complex bifurcation phenomenon. When there is one positive equilibrium on each side of , the model undergoes Hopf bifurcation with codimension up to 2. Moreover, numerical simulation reveals that the model not only can exhibit the krill paradox phenomenon but also has three limit cycles, with the outmost one crosses the line under some specific parameter conditions.

Bifurcations and global dynamics of a predator–prey mite model of Leslie type

AbstractIn this paper, we study a predator–prey mite model of Leslie type with generalized Holling IV functional response. The model is shown to have very rich bifurcation dynamics, including subcritical and supercritical Hopf bifurcations, degenerate Hopf bifurcation, focus‐type and cusp‐type degenerate Bogdanov–Takens bifurcations of codimension 3, originating from a nilpotent focus or cusp of codimension 3 that acts as the organizing center for the bifurcation set. Coexistence of multiple steady states, multiple limit cycles, and homoclinic cycles is also found. Interestingly, the coexistence of two limit cycles is guaranteed by investigating generalized Hopf bifurcation and degenerate homoclinic bifurcation, and we also find that two generalized Hopf bifurcation points are connected by a saddle‐node bifurcation curve of limit cycles, which indicates the existence of global regime for two limit cycles. Our work extends some results in the literature.

Degenerate codimension-2 cusp of limit cycles in a Holling–Tanner model with harvesting and anti-predator behavior

In this paper, a Holling–Tanner predator–prey model with constant-yield prey harvesting and anti-predator behavior is investigated. It is shown that the codimension of a nilpotent cusp is 3 and a degenerate Bogdanov–Takens bifurcation of codimension 3 acts as an organizing center for rich dynamical behaviors including saddle–node bifurcation, Hopf bifurcation, saddle–node bifurcation of limit cycles, and homoclinic bifurcation. Moreover, a Hopf bifurcation of codimension 2 is found, a codimension-2 cusp of limit cycles is firstly observed, and the corresponding normal form coefficients are also obtained, which indicate that there exists an acute angle region of three coexistent limit cycles. In particular, a reversed S-shaped saddle–node bifurcation curve of limit cycles is found. In addition, we observe that the anti-predator behavior may cause nearly extinction of the predator population while the prey population reaches its maximum which is less than the carrying capacity due to the effect of harvesting; however, the constant-yield prey harvesting may drive both species to extinction suddenly. Finally, numerical simulations including bifurcation diagrams and phase portraits are given to illustrate the theoretical results.

Bifurcations driven by generalist and specialist predation: mathematical interpretation of Fennoscandia phenomenon.

In this paper, we revisit a predator-prey model with specialist and generalist predators proposed by Hanski et al. (J Anim Ecol 60:353-367, 1991) , where the density of generalist predators is assumed to be a constant. It is shown that the model admits a nilpotent cusp of codimension 4 or a nilpotent focus of codimension 3 for different parameter values. As the parameters vary, the model can undergo cusp type (or focus type) degenerate Bogdanov-Takens bifurcations of codimension 4 (or 3). Our results indicate that generalist predation can induce more complex dynamical behaviors and bifurcation phenomena, such as three small-amplitude limit cycles enclosing one equilibrium, one or two large-amplitude limit cycles enclosing one or three equilibria, three limit cycles appearing in a Hopf bifurcation of codimension 3 and dying in a homoclinic bifurcation of codimension 3. In addition, we show that generalist predation stabilizes the limit cycle driven by specialist predators to a stable equilibrium, which clearly explains the famous Fennoscandia phenomenon.

Bifurcations in Holling-Tanner model with generalist predator and prey refuge

Refuge provides an important mechanism for preserving many ecosystems. Prey refuges directly benefit prey but also indirectly benefit predators in the long term. In this paper, we consider the complex dynamics and bifurcations in Holling-Tanner model with generalist predator and prey refuge. It is shown that the model admits a nilpotent cusp or focus of codimension 3, a nilpotent elliptic singularity of codimension at least 4, and a weak focus with order at least 3 for different parameter values. As the parameters vary, the model can undergo three types degenerate Bogdanov-Takens bifurcations of codimension 3 (cusp, focus and elliptic cases), and degenerate Hopf bifurcation of codimension 3. The system can exhibit complex dynamics, such as multiple coexistent periodic orbits and homoclinic loops. Moreover, our results indicate that the constant prey refuge prevents prey extinction and causes global coexistence. A preeminent finding is that refuge can induce a stable, large-amplitude limit cycle enclosing one or three positive steady states. Numerical simulations are provided to illustrate and complement our theoretical results.

An SIRS model with nonmonotone incidence and saturated treatment in a changing environment.

Nonmonotone incidence and saturated treatment are incorporated into an SIRS model under constant and changing environments. The nonmonotone incidence rate describes the psychological or inhibitory effect: when the number of the infected individuals exceeds a certain level, the infection function decreases. The saturated treatment function describes the effect of infected individuals being delayed for treatment due to the limitation of medical resources. In a constant environment, the model undergoes a sequence of bifurcations including backward bifurcation, degenerate Bogdanov-Takens bifurcation of codimension 3, degenerate Hopf bifurcation as the parameters vary, and the model exhibits rich dynamics such as bistability, tristability, multiple periodic orbits, and homoclinic orbits. Moreover, we provide some sufficient conditions to guarantee the global asymptotical stability of the disease-free equilibrium or the unique positive equilibrium. Our results indicate that there exist three critical values r_1, r_2 and r_3 for the treatment rate r: (i) when rge max {r_1, r_2}, the disease will disappear; (ii) when r<min {r_1, r_3}, the disease will persist. In a changing environment, the infective population starts along the stable disease-free state (or an endemic state) and surprisingly continues tracking the unstable disease-free state (or a limit cycle) when the system crosses a bifurcation point, and eventually tends to the stable endemic state (or the stable disease-free state). This transient tracking of the unstable disease-free state when {mathscr {R}}_0>1 predicts regime shifts that cause the delayed disease outbreak in a changing environment. Furthermore, the disease can disappear in advance (or belatedly) if the rate of environmental change is negative and large (or small). The transient dynamics of an infectious disease heavily depend on the initial infection number and rate or the speed of environmental change. Supplementary InformationThe online version contains supplementary material available at 10.1007/s00285-022-01787-3.

Bifurcation analysis of a predator–prey model with Beddington–DeAngelis functional response and predator competition

In this paper, we consider a predator–prey model with Beddington–DeAngelis functional response and predator competition, which is a five‐parameter family of planar vector field. It is shown that the model can undergo a sequence of bifurcations including focus type degenerate Bogdanov–Takens bifurcation of codimension 3 and Hopf bifurcation of codimension at least 2 as the parameters vary. Our theoretical results indicate that predator competition can cause richer dynamics such as two limit cycles enclosing one or three hyperbolic positive equilibria and three kinds of homoclinic orbits (homoclinic to hyperbolic saddle, saddle‐node, or neutral saddle). Moreover, there exists a threshold value for predator capturing rate , below or equal to which the predators always tend to extinction, above which the predators and preys will coexist in the form of multiple steady states or periodic oscillations for all positive initial populations. Numerical simulations are presented to illustrate the theoretical results.

Stability Analysis of a Leslie–Gower Model with Strong Allee Effect on Prey and Fear Effect on Predator

In this paper, we propose a Leslie–Gower predator–prey model with strong Allee effect on prey and fear effect on predator. We discuss the existence and local stability of equilibria by making full use of qualitative analytical theory. It is shown that the above system exhibits at most two positive equilibria and it can undergo a series of bifurcation phenomena. We indicate that the dynamical behavior of the model is closely related to the fear effect on predator. In detail, when the fear effect parameter [Formula: see text], the system will undergo degenerate Hopf bifurcation. There exist two limit cycles (the inner is stable and the outer is unstable). However, when [Formula: see text], the system will undergo degenerate Bogdanov–Takens bifurcation. Also, by numerical simulation, we conclude that the stronger the fear effect, the bigger the density of prey species. The above shows that fear effect on predator is beneficial to the persistence of the prey species. Our results can be seen as a complement to previous works [González-Olivares et al., 2011; Pal &amp; Mandal, 2014].

Linking bifurcation analysis of Holling–Tanner model with generalist predator to a changing environment

AbstractBifurcation theory has been highly popular in the analysis of mathematical models. However, stability and bifurcation analyses are only for asymptotic dynamics while applied scientists care more about transient dynamics. In this paper, we first rigorously analyze Holling–Tanner model with generalist predators who have alternative food sources, and then discuss transient dynamics via a changing environment. For a constant environment, we provide a complete bifurcation analysis with high codimension. It is shown that the highest codimension of a nilpotent cusp is 3, and the model can undergo degenerate Bogdanov–Takens bifurcation of codimension 3. Moreover, by using resultant elimination to solve the semialgebraic varieties of Lyapunov coefficients, we show that a center‐type equilibrium is a weak focus with order at most 2, and the model can exhibit Hopf bifurcation of codimension 2. Our results indicate that generalist predators can cause not only richer dynamics and bifurcations, but also the extinction of prey for some positive initial densities. Numerical simulations, including the coexistence of a limit cycle and a homoclinic cycle, tristability, two limit cycles, are presented to illustrate the theoretical results. In a changing environment, the populations start along one stable state but can track unstable states or oscillations when the system crosses a bifurcation point, and then tend to another stable state or oscillations. This tracking on transient dynamics predicts regime shifts under environmental changes. When environmental conditions vary, the populations can track unstable states in the constant environment. The rate of environmental change determines how long the system tracks an unstable state although finally the solution under environmental change is attracted to a stable steady state or limit cycle. Finally, we focus on a periodic environment and find that the populations converge to a periodic solution or an invariant torus depending on both the initial environmental capacity and the amplitude of periodic fluctuation.

Degenerate Bogdanov-Takens bifurcation of codimension 4 in Holling-Tanner model with harvesting

In this paper, we revisit the Holling-Tanner model with constant-yield prey harvesting. It is shown that the highest codimension of a nilpotent cusp is 4, and the model can undergo degenerate Bogdanov-Takens bifurcation of codimension 4. Moreover, when the model has a center-type equilibrium, we show that it is a weak focus with order at least 3 and at most 4, and the model can exhibit Hopf bifurcation of codimension 3. Some algebraic methods including resultant elimination and pseudo-division are used to solve the semi-algebraic varieties of normal form coefficients or focal values. Our results indicate that constant-yield prey harvesting can cause not only richer dynamics and bifurcations, but also the coextinction of both populations with some positive initial densities. Finally, numerical simulations, including the coexistence of limit cycle and homoclinic cycle, and three limit cycles, are presented to illustrate the theoretical results.

Nilpotent singularities and bifurcations structure of the Landau–Ginzburg theory of cortex dynamics

The dynamics research of cerebral cortex activity is an important part of neuroscience. A previous work (Proc. Natl. Acad. Sci. USA 115 (2018) E1356-E1365) put forward a cortex dynamics model based on Landau–Ginzburg theory and it generated an oscillation and bistability. However, the full dynamics of cerebral cortex is still ill-defined. In this paper, we rigorously prove that this model can undergo Bogdanov–Takens bifurcation of codimension 2, Bogdanov–Takens bifurcations of codimension 3 with cusp, focus and elliptic types, and doubly degenerate Bogdanov–Takens bifurcation of codimension 4. Our work complements and enriches the results in the field of cerebral cortex activity.

Global analysis in Bazykin's model with Holling II functional response and predator competition

In this paper, we study the well-known Bazykin's model with Holling II functional response and predator competition. A detailed bifurcation analysis, depending on all four parameters, reveals a rich bifurcation structure, including supercritical and subcritical Bogdanov-Takens bifurcation, degenerate Hopf bifurcation of codimension at most 2, and a focus type degenerate Bogdanov-Takens bifurcation of codimension 3, originating from a nilpotent focus of codimension 3 which acts as the organizing center for the bifurcation set. Moreover, some sufficient conditions to guarantee the global asymptotical stability of the semi-trivial equilibrium or the unique positive equilibrium are also given. Our analysis indicates that we can classify the long-time dynamics of the model with a threshold value c0 for the natural mortality rate c of predators, in detail, the following are true. (i) When c≥c0, the prey will persist and predators will eventually go extinct for all positive initial populations. (ii) When c<c0, the prey and predators will coexist, for all positive initial populations, in the form of multiple positive equilibria or multiple periodic orbits. Our results can be seen as a complement to the work by Bazykin et al. [2–5], Hainzl [22,23], Kuznetsov [30].

Degenerate Bogdanov\u2013Takens bifurcations in a bulk viscous cosmology

Using the dynamical system theory we show that the Friedmann–Robertson–Walker (FRW) cosmological model with bulk viscous fluid in the presence of cosmological constant is equivalent to a degenerate two dimensional Bogdanov–Takens normal form. The equation of state parameter, omega , the bulk viscosity coefficient, xi , and the cosmological constant, Lambda , define the necessary parameters for unfolding the degenerate Bogdanov-Takens system. The fixed points of the system are discussed together with the variation of their stability properties upon changing the relevant parameters omega , Lambda and xi . The variation of the stability properties are visualized by the appropriate bifurcation diagrams. Phase portrait for finite domain and global phase portrait are displayed and the issue of the structural stability is discussed. Typical issues such as late acceleration or inflation that can be induced by viscosity and could have relevance to observational cosmology are also discussed.

Host-virus evolutionary dynamics with specialist and generalist infection strategies: Bifurcations, bistability, and chaos.

In this work, we have investigated the evolutionary dynamics of a generalist pathogen, e.g., a virus population, that evolves toward specialization in an environment with multiple host types. We have particularly explored under which conditions generalist viral strains may rise in frequency and coexist with specialist strains or even dominate the population. By means of a nonlinear mathematical model and bifurcation analysis, we have determined the theoretical conditions for stability of nine identified equilibria and provided biological interpretation in terms of the infection rates for the viral specialist and generalist strains. By means of a stability diagram, we identified stable fixed points and stable periodic orbits, as well as regions of bistability. For arbitrary biologically feasible initial population sizes, the probability of evolving toward stable solutions is obtained for each point of the analyzed parameter space. This probability map shows combinations of infection rates of the generalist and specialist strains that might lead to equal chances for each type becoming the dominant strategy. Furthermore, we have identified infection rates for which the model predicts the onset of chaotic dynamics. Several degenerate Bogdanov-Takens and zero-Hopf bifurcations are detected along with generalized Hopf and zero-Hopf bifurcations. This manuscript provides additional insights into the dynamical complexity of host-pathogen evolution toward different infection strategies.

Bifurcation analysis in a host-generalist parasitoid model with Holling II functional response

In this paper we study a host-generalist parasitoid model with Holling II functional response where the generalist parasitoids are introduced to control the invasion of the hosts. It is shown that the model can undergo a sequence of bifurcations including cusp, focus and elliptic types degenerate Bogdanov-Takens bifurcations of codimension three, and a degenerate Hopf bifurcation of codimension at most two as the parameters vary, and the model exhibits rich dynamics such as the existence of multiple coexistent steady states, multiple coexistent periodic orbits, homoclinic orbits, etc. Moreover, there exists a critical value for the carrying capacity of generalist parasitoids such that: (i) when the carrying capacity of the generalist parasitoids is smaller than the critical value, the invading hosts can always persist despite of the predation by the generalist parasitoids, i.e., the generalist parasitoids cannot control the invasion of hosts; (ii) when the carrying capacity of the generalist parasitoids is larger than the critical value, the invading hosts either tend to extinction or persist in the form of multiple coexistent steady states or multiple coexistent periodic orbits depending on the initial populations, i.e., whether the invasion can be stopped and reversed by the generalist parasitoids depends on the initial populations; (iii) in both cases, the generalist parasitoids always persist. Numerical simulations are presented to illustrate the theoretical results.

Approximations of the Heteroclinic Orbits Near a Double-Zero Bifurcation with Symmetry of Order Two. Application to a Liénard Equation

A dynamical system possessing an equilibrium point with two zero eigenvalues is considered. We assume that a degenerate Bogdanov–Takens bifurcation with symmetry of order two is present and, in the parameter space, a curve of heteroclinic bifurcation values emerges from the codimension two bifurcation point. Using a blow-up transformation and a perturbation method, we obtain second order approximations both for the heteroclinic orbits and for the curve of heteroclinic bifurcation values. Applications of our results for the truncated normal form and for a Liénard equation are presented. Some numerical simulations illustrating the accuracy of our results are performed.

Theoretical conditions for the coexistence of viral strains with differences in phenotypic traits: a bifurcation analysis.

We investigate the dynamics of a wild-type viral strain which generates mutant strains differing in phenotypic properties for infectivity, virulence and mutation rates. We study, by means of a mathematical model and bifurcation analysis, conditions under which the wild-type and mutant viruses, which compete for the same host cells, can coexist. The coexistence conditions are formulated in terms of the basic reproductive numbers of the strains, a maximum value of the mutation rate and the virulence of the pathogens. The analysis reveals that parameter space can be divided into five regions, each with distinct dynamics, that are organized around degenerate Bogdanov–Takens and zero-Hopf bifurcations, the latter of which gives rise to a curve of transcritical bifurcations of periodic orbits. These results provide new insights into the conditions by which viral populations may contain multiple coexisting strains in a stable manner.