Ecological consequences of global bifurcations in some food chain models

Food chain models of ordinary differential equations (ode’s) are often used in ecology to gain insight in the dynamics of populations of species, and the interactions of these species with each other and their environment. One powerful analysis technique is bifurcation analysis, focusing on the changes in long-term (asymptotic) behaviour under parameter variation. For the detection of local bifurcations there exists standardised software, but until quite recently most software did not include any capabilities for the detection and continuation of global bifurcations. We focus here on the occurrence of global bifurcations in four food chain models, and discuss the implications of their occurrence. In two stoichiometric models (one piecewise continuous, one smooth) there exists a homoclinic bifurcation, that results in the disappearance of a limit cycle attractor. Instead, a stable positive equilibrium becomes the global attractor. The models are also capable of bistability. In two three-dimensional models a Shil’nikov homoclinic bifurcation functions as the organising centre of chaos, while tangencies of homoclinic cycle-to-cycle connections ‘cut’ the chaotic attractors, which is associated with boundary crises. In one model this leads to extinction of the top predator, while in the other model hysteresis occurs. The types of ecological events occurring because of a global bifurcation will be categorized. Global bifurcations are always catastrophic, leading to the disappearance or merging of attractors. However, there is no 1-on-1 coupling between global bifurcation type and the possible ecological consequences. This only emphasizes the importance of including global bifurcations in the analysis of food chain models.