Degenerate Bogdanov Takens Bifurcations Research Articles

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.

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.

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].

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.

Global Analysis of a Liénard System with Quadratic Damping

In this paper, the global analysis of a Liénard equation with quadratic damping is studied. There are 22 different global phase portraits in the Poincaré disc. Every global phase portrait is given as well as the complete global bifurcation diagram. Firstly, the equilibria at finite and infinite of the Liénard system are discussed. The properties of the equilibria are studied. Then, the sufficient and necessary conditions of the system with closed orbits are obtained. The degenerate Bogdanov-Takens bifurcation is studied and the bifurcation diagrams of the system are given.

Degenerate Bogdanov–Takens bifurcations in a one-dimensional transport model of a fusion plasma

Experiments in tokamaks (nuclear fusion reactors) have shown two modes of operation: L-mode and H-mode. Transitions between these two modes have been observed in three types: sharp, smooth and oscillatory. The same modes of operation and transitions between them have been observed in simplified transport models of the fusion plasma in one spatial dimension. We study the dynamics in such a one-dimensional transport model by numerical continuation techniques. To this end the MATLAB package cl_matcontL was extended with the continuation of (codimension-2) Bogdanov–Takens bifurcations in three parameters using subspace reduction techniques. During the continuation of (codimension-2) Bogdanov–Takens bifurcations in 3 parameters, generically degenerate Bogdanov–Takens bifurcations of codimension-3 are detected. However, when these techniques are applied to the transport model, we detect a degenerate Bogdanov–Takens bifurcation of codimension 4. The nearby 1- and 2-parameter slices are in agreement with the presence of this codimension-4 degenerate Bogdanov–Takens bifurcation, and all three types of L–H transitions can be recognized in these slices. The same codimension-4 situation is observed under variation of the additional parameters in the model, and under some modifications of the model.

Bifurcation analysis on a generalized recurrent neural network with two interconnected three-neuron components

A class of recurrent neural networks is constructed by generalizing a specific class of n-neuron networks. It is shown that the newly constructed network experiences generic pitchfork and Hopf codimension one bifurcations. It is also proved that the emergence of generic Bogdanov–Takens, pitchfork–Hopf and Hopf–Hopf codimension two, and the degenerate Bogdanov–Takens bifurcation points in the parameter space is possible due to the intersections of codimension one bifurcation curves. The occurrence of bifurcations of higher codimensions significantly increases the capability of the newly constructed recurrent neural network to learn broader families of periodic signals.

The saddle-node-transcritical bifurcation in a population model with constant rate harvesting

We study the interaction of saddle-node and transcritical bifurcations in a Lotka-Volterra model with a constant term representing harvesting or migration. Because some of the equilibria of the model lie on an invariant coordinate axis, both the saddle-node and the transcritical bifurcations are of codimension one. Their interaction can be associated with either a single or a double zero eigenvalue. We show that in the former case, the local bifurcation diagram is given by a nonversal unfolding of the cusp bifurcation whereas in the latter case it is a nonversal unfolding of a degenerate Bogdanov-Takens bifurcation. We present a simple model for each of the two cases to illustrate the possible unfoldings. We analyse the consequences of the generic phase portraits for the Lotka-Volterra system.