Heteroclinic orbits indicate overexploitation in predator-prey systems with a strong Allee effect.

Species establishment in a model system in a homogeneous environment can be dependent not only on the parameter setting, but also on the initial conditions of the system. For instance, predator invasion into an established prey population can fail and lead to system collapse, an event referred to as overexploitation. This phenomenon occurs in models with bistability properties, such as strong Allee effects. The Allee effect then prevents easy re-establishment of the prey species. In this paper, we deal with the bifurcation analyses of two previously published predator-prey models with strong Allee effects. We expand the analyses to include not only local, but also global bifurcations. We show the existence of a point-to-point heteroclinic cycle in these models, and discuss numerical techniques for continuation in parameter space. The continuation of such a cycle in two-parameter space forms the boundary of a region in parameter space where the system collapses after predator invasion, i.e. where overexploitation occurs. We argue that the detection and continuation of global bifurcations in these models are of vital importance for the understanding of the model dynamics.