Resonance properties and chaotic dynamics of a three-dimensional discrete logistic ecological system within the neighborhoods of bifurcation points
Yujiang Chen, Lin Li, Lingling Liu, Zhiheng Yu
TL;DR
The paper addresses the rich dynamics of a three-dimensional discrete-time ecological system, modeled by a map $F:\,\mathbb{R}_+^3\to\mathbb{R}_+^3$ with fixed points $O=(0,0,0)$, $E_1=(1-1/\mu,0,0)$, $E_2=(1/\beta,1-1/\mu-1/\beta,0)$, and $E_3$, by applying complete polynomial discriminant theory, center-manifold reductions, and normal-form analysis. It derives exhaustive codimension-1 bifurcations (transcritical, flip, Neimark-Sacker) near the fixed points and maps the ensuing dynamics, including strong resonances $1:2$, $1:3$, $1:4$ and weak-resonance Arnold tongues, complemented by Marotto chaos criteria. The contributions include rigorous bifurcation proofs, explicit normal-form unfoldings, and detailed resonance sequences, supported by numerical simulations that illustrate invariant circles, period-doubling, and chaotic regimes. The results enhance understanding of multi-species ecological interactions by revealing how parameter changes can drive transitions between extinction, persistent oscillations, and chaotic population dynamics, with implications for ecosystem management and resilience.
Abstract
In this paper, we delve into the dynamical properties of a class of three-dimensional logistic ecological models. By using the complete discriminant theory of polynomials, we first give a topological classification for each fixed point and investigate the stability of corresponding system near the fixed points. Then employing the bifurcation and normal form theory, we discuss all possible codimension-1 bifurcations near the fixed points, i.e., transcritical, flip, and Neimark-Sacker bifurcations, and further prove that the system can undergo codimension-2 bifurcations, specifically 1:2, 1:3, 1:4 strong resonances and weak resonance Arnold tongues. Additionally, chaotic behaviors in the sense of Marotto are rigorously analyzed. Numerical simulations are conducted to validate the theoretical findings and illustrate the complex dynamical phenomena identified.
