Dynamics of coral reef models in the presence of parrotfish

The decline of coral reefs characterized by macroalgae increase has been a global threat. We consider a slightly modified version of an ordinary differential equation (ODE) model proposed in Blackwood, Hastings, and Mumby [Theor. Ecol. 5 (2012), pp. 105-114] that explicitly considers the role of parrotfish grazing on coral reef dynamics. We perform complete stability, bifurcation, and persistence analysis for this model. If the fishing effort (f) is in between two critical values f0 and f1, then the system has a unique interior equilibrium, which is stable if f0<f<f1 and unstable if f1<f<f0. If f is less (more) than these critical values, then the system has up to two (zero) interior equilibria. Also, we develop a more realistic delay differential equation (DDE) model to incorporate the time delay and treating it as the bifurcation parameter, and we prove that Hopf bifurcation about the interior equilibria could occur at critical time delays, which illustrate the potential importance of the inherent time delay in a coral reef ecosystem. Recommendations for Resource Managers One serious threat to coral reefs is overfishing of grazing species, including high level of algal abundance. Fishing alters the entire dynamics of a reef (Hughes, Baird, & Bellwood, 2003), for which the coral cover was predicted to decline rapidly (Mumby, 2006). One major issue is to reverse and develop appropriate management to increase or maintain coral resilience. We have provided a detailed local and global analysis of model (Blackwood, Hastings, & Mumby, 2012) and obtained an ecologically meaningful attracting region, for which there is a chance of stable coexistence of coral-algal-fish state. The healthy reefs switch to unhealthy state, and the macroalgae-parrotfish state becomes stable as the fishing effort increases through some critical values. Also, for some critical time delays, a switch between healthy and unhealthy reef states occurs through a Hopf bifurcation, which can only appear in the delay differential equation (DDE) model. Eventually, for large enough time delay, oscillations appear and an unhealthy state occurs.