This study investigates the dynamical interaction between guano biomass density (x), invertebrate biomass density (y), and fish biomass density (z) through a three-compartment trophic model representing a nutrient-based cave ecosystem. The analysis identifies three equilibrium points: the consumer-free equilibrium E0, the predator-free equilibrium E1, and the coexistence equilibrium E2. Local stability analysis shows that the coexistence equilibrium is asymptotically stable, characterized by eigenvalues with negative real parts (1 = -0.519706 and 2,3 = -0.085385 0.188169i). Numerical simulations using the fourth–fifth order Runge–Kutta method (RK45) support these analytical results, showing trajectories that exhibit damped oscillations before converging to the steady state. Furthermore, a bifurcation analysis reveals a critical Branching Point (BP) at the predation rate b2 0.043956. This threshold signifies a transcritical bifurcation where the system transitions from a predator-extinction regime to a stable coexistence regime, highlighting the sensitivity of the food web to energy transfer efficiency. These findings suggest that under the assumed parameter set, the interaction between guano nutrients, invertebrates, and fish can maintain a stable ecological balance through top-down control and nutrient-dependent dynamics.
Copyrights © 2026