A Rosenzweig-MacArthur Model with Continuous Threshold Harvesting in Predator Involving Fractional Derivatives with Power Law and Mittag-Leffler Kernel

The harvesting management is developed to protect the biological resources from over-exploitation such as harvesting and trapping. In this article, we consider a predator-prey interaction that follows the fractional-order Rosenzweig-MacArthur model where the predator is harvested obeying a threshold harvesting policy (THP). The THP is applied to maintain the existence of the population in the prey-predator mechanism. We first consider the Rosenzweig-MacArthur model using the Caputo fractional-order derivative (that is, the operator with the power-law kernel) and perform some dynamical analysis such as the existence and uniqueness, non-negativity, boundedness, local stability, global stability, and the existence of Hopf bifurcation. We then reconsider the same model involving the Atangana-Baleanu fractional derivative with the Mittag-Leffler kernel in the Caputo sense (ABC). The existence and uniqueness of the solution of the model with ABC operator are established. We also explore the dynamics of the model with both fractional derivative operators numerically and confirm the theoretical findings. In particular, it is shown that models with both Caputo operator and ABC operator undergo a Hopf bifurcation that can be controlled by the conversion rate of consumed prey into the predator birth rate or by the order of fractional derivative. However, the bifurcation point of the model with the Caputo operator is different from that of the model with the ABC operator.