Complex dynamics of a tri-trophic predator-prey system with double fear effects

Ting Yang
Department of Mathematics, Jiangxi University of Finance and Economics, Nanchang 330013, China
Yue Cai
Department of Mathematics, Jiangxi University of Finance and Economics, Nanchang 330013, China

Article ID: 4577

Keywords: predator-prey system; double fear effects; equilibrium; stability; Hopf bifurcation; coexistence; chaos

Abstract

In this paper, we investigate a tri-trophic predator-prey system with double fear effects, where the prey is affected by fear induced by both the intermediate predator and the top predator. The top predator is assumed to feed on both the basal prey and the intermediate predator. We establish the positivity and boundedness of solutions and derive sufficient conditions for the global asymptotic stability of axial equilibria. The possible interior equilibria are characterized algebraically, and their local stability is analyzed. We further derive the Hopf bifurcation conditions and compute the first Lyapunov coefficient to determine the direction of the bifurcation and the stability of the bifurcating periodic orbit. Numerical simulations reveal several forms of complex dynamics, including a period-three orbit in the associated return dynamics, the coexistence of a stable equilibrium and a stable limit cycle, the coexistence of two stable limit cycles, and chaotic behavior. The corresponding multistability is further illustrated by basin-of-attraction diagrams. Numerical evidence for chaotic attractors in parameter regions admitting different numbers of interior equilibria is provided by positive largest Lyapunov exponents, parameter bifurcation diagrams, and bounded long-time trajectories. These results show that, within the parameter regimes examined, the interaction between the two fear mechanisms and the tri-trophic feeding structure can generate pronounced dynamical transitions, multistability, and strong sensitivity to initial population densities.

Published

2026-08-21

How to Cite

Yang, T., & Cai, Y. (2026). Complex dynamics of a tri-trophic predator-prey system with double fear effects. Advances in Differential Equations and Control Processes, 33(3). https://doi.org/10.59400/adecp4577

Download Citation

Issue

Vol. 33 No. 3 (2026): In Progress

Section

Article

This work is licensed under a Creative Commons Attribution 4.0 International License.

References

[1]Lotka AJ. Elements of Physical Biology. Williams & Wilkins; 1925.

[2]Volterra V. Variations and Fluctuations in the Number of Individuals in Coexisting Animal Species. Leonardo da

[3]Vinci Printing Company; 1926. (in Italian)

[4]Holling CS. The components of predation as revealed by a study of small-mammal predation of the European pine

[5]sawfly. Canadian Entomologist. 1959; 91(5): 293–320.

[6]May RM. Stability and Complexity in Model Ecosystems. Princeton University Press; 2001.

[7]Zhang L, Wang ZC. Spatial dynamics of a diffusive predator–prey model with stage structure. Discrete and

[8]Continuous Dynamical Systems–Series B. 2015; 20(6): 1831–1853.

[9]Guin LN, Pal PJ, Alzahrani J, et al. Influence of Allee effect on the spatiotemporal behavior of a diffusive

[10]predator–prey model with Crowley–Martin type response function. Scientific Reports. 2023; 13: 4710.

[11]Pandey S, Ghosh U, Das D, et al. Rich dynamics of a delay-induced stage-structure prey–predator model

[12]with cooperative behaviour in both species and the impact of prey refuge. Mathematics and Computers in Simulation.

[13]; 216: 49–76.

[14]Shao Y, Yu H, Jin C, et al. Dynamics analysis of a predator–prey model with Allee effect and harvesting effort.

[15]Electronic Research Archive. 2024; 32(9): 5682–5716.

[16]Ge J, Ji W, Liu M. Threshold behavior of a stochastic predator–prey model with fear effect and regime-switching.

[17]Applied Mathematics Letters. 2025; 164: 109476.

[18]Paul B, Sikdar GC, Ghosh U. Effect of fear and non-linear predator harvesting on a predator–prey system in presence

[19]of environmental variability. Mathematics and Computers in Simulation. 2025; 227: 442–460.

[20]Qi H, Jin J, Liu B, et al. Stationary distribution of a stochastic reaction–diffusion predator–prey model with additional

[21]food, fear effect and anti-predator behavior. Applied Mathematics Letters. 2025; 169: 109612.

[22]Hastings A, Powell T. Chaos in a three-species food chain. Ecology. 1991; 72(3): 896–903.

[23]Klebanoff A, Hastings A. Chaos in three-species food chains. Journal of Mathematical Biology. 1994; 32(5):

[24]–451.

[25]Ramasamy S, Banjerdpongchai D, Park P. Chaos control of a delayed tri-trophic food chain model with fear and its

[26]carry over effects. Symmetry. 2023; 15(2): 484.

[27]Qin W, Dong Z. The impact of resource limitation on the pest-natural enemy ecosystem with anti-predator behavior

[28]and fear effect. Advances in Continuous and Discrete Models. 2024; 2024: 10.

[29]Brown JS, Laundré JW, Gurung M. The ecology of fear: Optimal foraging, game theory, and trophic interactions.

[30]Journal of Mammalogy. 1999; 80(2): 385–399.

[31]Preisser EL, Bolnick DI, Benard MF. Scared to death? The effects of intimidation and consumption in predator–prey

[32]interactions. Ecology. 2005; 86(2): 501–509.

[33]Zanette LY, White AF, Allen MC, et al. Perceived predation risk reduces the number of offspring songbirds produce

[34]per year. Science. 2011; 334(6061): 1398–1401.

[35]Wang X, Zanette L, Zou X. Modelling the fear effect in predator–prey interactions. Journal of Mathematical Biology.

[36]; 73(5): 1179–1204.

[37]Kumar A, Dubey B. Modeling the effect of fear in a prey–predator system with prey refuge and gestation delay.

[38]International Journal of Bifurcation and Chaos. 2019; 29(14): 1950195.

[39]Wang X, Zou X. Modeling the fear effect in predator–prey interactions with adaptive avoidance of predators. Bulletin

[40]of Mathematical Biology. 2017; 79(6): 1325–1359.

[41]Pal S, Pal N, Samanta S, et al. Fear effect in prey and hunting cooperation among predators in a Leslie–Gower model.

[42]Mathematical Biosciences and Engineering. 2019; 16(5): 5146–5179.

[43]Panday P, Pal N, Samanta S, et al. Stability and bifurcation analysis of a three-species food chain model with fear.

[44]International Journal of Bifurcation and Chaos. 2018; 28(1): 1850009.

[45]Cong P, Fan M, Zou X. Dynamics of a three-species food chain model with fear effect. Communications in Nonlinear

[46]Science and Numerical Simulation. 2021; 99: 105809.

[47]Kumar V, Kumari N. Controlling chaos in three-species food chain model using fear effect. AIMS Mathematics.

[48]; 5(2): 828–842.

[49]Rana S, Bhattacharya S, Samanta S. Complex dynamics of a three-species food chain model with fear and Allee

[50]effect. International Journal of Bifurcation and Chaos. 2022; 32(10): 2250084.

[51]Hsu SB, Ruan S, Yang TH. Analysis of three-species Lotka–Volterra food web models with omnivory. Journal of

[52]Mathematical Analysis and Applications. 2015; 426(1): 659–687.

[53]Kratina P, LeCraw RM, Ingram T, et al. Stability and persistence of food webs with omnivory: Is there a general

[54]pattern? Ecosphere. 2012; 3(6): 50.

[55]Arancibia-Ibarra C, Flores J. Dynamics of a Leslie–Gower predator–prey model with Holling type II functional

[56]response, Allee effect and a generalist predator. Mathematics and Computers in Simulation. 2021; 188: 1–22.

[57]Sk N, Mondal B, Sarkar A, et al. Chaos emergence and dissipation in a three-species food web model with intraguild

[58]predation and cooperative hunting. AIMS Mathematics. 2024; 9(1): 1023–1045.

[59]Wang T, Wang H. Stoichiometry-dependent fear effect in a food chain model. SIAM Journal on Applied Mathematics.

[60]; 84(3): 783–807.

[61]Han R, Guin LN, Dai B. Consequences of refuge and diffusion in a spatiotemporal predator–prey model. Nonlinear

[62]Analysis: Real World Applications. 2021; 60: 103311.

[63]Pal S, Karmakar S, Pal S, et al. Impact of fear and group defense on the dynamics of a predator–prey system.

[64]International Journal of Bifurcation and Chaos. 2024; 34(2): 2450019.

[65]Cui W, Zhao Y. Bifurcation analysis of a predator–prey model with alternative prey and prey refuges. International

[66]Journal of Bifurcation and Chaos. 2024; 34(2): 2450021.

[67]Din Q, Naseem RA, Shabbir MS. Predator–prey interaction with fear effects: Stability, bifurcation and two-parameter

[68]analysis incorporating complex and fractal behavior. Fractal and Fractional. 2024; 8(4): 221.

[69]Saha S, Pal S, Melnik R. Analysis of the impact of fear in the presence of additional food and prey refuge with

[70]nonlocal predator–prey models. Ecological Modelling. 2025; 505: 111103.

[71]Li TY, Yorke JA. Period three implies chaos. The American Mathematical Monthly. 1975; 82(10): 985–992.

Early fault signals of the rolling bearing in the rotor are weak and present the characteristics of non-periodic and non-stationary; it is more difficult to carry out fault diagnosis on it. In this regard, this paper proposes a weak rolling bearing fault diagnosis algorithm based on whale optimization algorithm, simplistic geometry mode decomposition, and maximum correlated kurtosis deconvolution (WOA-SGMD-MCKD). Firstly, the vibration signal of the rotor platform is obtained, and the Symmetric Geometric Mode Decomposition (SGMD) is used to reconstruct the vibration signal. To obtain the best decomposition effect of the SGMD and overcome modal aliasing, the Whale Optimization Algorithm (WOA) is used to optimize the embedding dimension. Secondly, for the reconstructed vibration signal, the Maximum Correlated Kurtosis Deconvolution (MCKD) is used to extract its impulse component, and the WOA is used to optimize the filter length and deconvolution period of the MCKD so that the frequency envelope spectrum of the vibration signal can be obtained, which can provide the basis for the fault diagnosis of rolling bearings. Finally, the effectiveness and feasibility of the algorithm proposed are verified by a non-periodic and non-stationary simulation platform and rotor maneuvering platform in this paper.

Numerical solution of a 3D mathematical model for the progression of tumor angiogenic factor in a tissue

In this work, the movement of tumor angiogenic factor in a three-dimensional tissue is obtained by the Method of Lines. This method transforms a partial differential equation into a system of ordinary differential equations together with the initial and boundary conditions. The more the number of lines is increased, the more the accuracy of the method increases. This method results in very accurate numerical solutions for linear and non-linear problems in contrast with other existing methods. We present Matlab-generated figures, which are the movement of tumor angiogenic factor in porous medium and explain the biological importance of this progression. The computer codes are also provided.

Mathematical modelling and controllability analysis of fractional order coal mill pulverizer model

This paper investigates the controllability of nonlinear dynamical systems and their applications, with a focus on fractional-order systems and coal mill models. A novel theorem is proposed, providing sufficient conditions for controllability, including constraints on the steering operator and nonlinear perturbation bounds. The theorem establishes the existence of a contraction mapping for the nonlinear operator, enabling effective control strategies for fractional systems. The methodology is demonstrated through rigorous proof and supported by an iterative algorithm for controller design. Additionally, the controllability of a coal mill system represented as a nonlinear differential system, is analyzed. The findings present new insights into the interplay of fractional dynamics and nonlinear systems, offering practical solutions for real-world control problems.