Dynamics in two nonsmooth predator–prey models with threshold harvesting

Considering a good pest control program should reduce the pest to levels acceptable to the public, we investigate the threshold harvesting policy on pests in two predator–prey models. Both models are nonsmooth and the aim of this paper is to provide how threshold harvesting affects the dynamics of the two systems. When the harvesting threshold is larger than some positive level, the harvesting does not affect the ecosystem; when the harvesting threshold is less than the level, the model has complex dynamics with multiple coexistence equilibria, limit cycle, bistability, homoclinic orbit, saddle-node bifurcation, transcritical bifurcation, subcritical and supercritical Hopf bifurcation, Bogdanov–Takens bifurcation, and discontinuous Hopf bifurcation. Firstly, we provide the complete stability analysis and bifurcation analysis for the two models. Furthermore, some numerical simulations are given to illustrate our results. Finally, it is found that harvesting lowers the level of both species for natural enemy–pest system while raises the densities of both species for the pest–crop system. It is seen that the threshold harvesting policy of the enemy system is more effective than the crop system.

[1]  S Mandal,et al.  Toxin-producing plankton may act as a biological control for planktonic blooms--field study and mathematical modelling. , 2002, Journal of theoretical biology.

[2]  Sanyi Tang,et al.  State-dependent impulsive models of integrated pest management (IPM) strategies and their dynamic consequences , 2005, Journal of mathematical biology.

[3]  Yongzhen Pei,et al.  Evolutionary consequences of harvesting for a two-zooplankton one-phytoplankton system , 2012 .

[4]  Yu. S. Ledyaev,et al.  Nonsmooth analysis and control theory , 1998 .

[5]  Bernt-Erik Sæther,et al.  THRESHOLD HARVESTING FOR SUSTAINABILITY OF FLUCTUATING RESOURCES , 1997 .

[6]  van Dh Dick Campen,et al.  Bifurcation phenomena in non-smooth dynamical systems , 2006 .

[7]  Lili Ji,et al.  Qualitative analysis of a predator–prey model with constant-rate prey harvesting incorporating a constant prey refuge , 2010 .

[8]  Benjamin Leard,et al.  Analysis of predator-prey models with continuous threshold harvesting , 2011, Appl. Math. Comput..

[9]  J. Hale,et al.  Dynamics and Bifurcations , 1991 .

[10]  Dongmei Xiao,et al.  Bifurcations of a Ratio-Dependent Predator-Prey System with Constant Rate Harvesting , 2005, SIAM J. Appl. Math..

[11]  Yongzhen Pei,et al.  Extinction and permanence of one-prey multi-predators of Holling type II function response system with impulsive biological control. , 2005, Journal of theoretical biology.

[12]  Jorge Rebaza,et al.  Dynamics of prey threshold harvesting and refuge , 2012, J. Comput. Appl. Math..

[13]  Remco I. Leine,et al.  Bifurcations of Equilibria in Non-smooth Continuous Systems , 2006 .

[14]  Xinyu Song,et al.  Effect of prey refuge on a harvested predator–prey model with generalized functional response , 2011 .

[15]  Li Zhong,et al.  Stability analysis of a prey-predator model with holling type III response function incorporating a prey refuge , 2006, Appl. Math. Comput..

[16]  Sanyi Tang,et al.  Multiple attractors of host-parasitoid models with integrated pest management strategies: eradication, persistence and outbreak. , 2008, Theoretical population biology.

[17]  Sanyi Tang,et al.  Modelling and analysis of integrated pest management strategy , 2004 .

[18]  Shenquan Liu,et al.  Codimension-two bifurcation analysis in two-dimensional Hindmarsh–Rose model , 2012 .

[19]  Darwin C. Hall,et al.  On the Timing and Application of Pesticides , 1973 .

[20]  Sanyi Tang,et al.  Models for integrated pest control and their biological implications. , 2008, Mathematical biosciences.

[21]  C. Shoemaker.  Optimal timing of multiple applications of pesticides with residual toxicity. , 1979, Biometrics.

[22]  Shaopu Yang,et al.  The codimension-two bifurcation for the recent proposed SD oscillator , 2009 .

[23]  David L. Sunding,et al.  Insect Population Dynamics, Pesticide Use, and Farmworker Health , 2000 .

[24]  V. Marcström,et al.  SUSTAINABLE HARVESTING STRATEGIES OF WILLOW PTARMIGAN IN A FLUCTUATING ENVIRONMENT , 2002 .

[25]  Yang Kuang,et al.  Uniqueness of limit cycles in Gause-type models of predator-prey systems , 1988 .

[26]  Sanyi Tang,et al.  Integrated pest management models and their dynamical behaviour , 2005, Bulletin of mathematical biology.

[27]  Sanyi Tang,et al.  Optimal dosage and economic threshold of multiple pesticide applications for pest control , 2010, Math. Comput. Model..

[28]  P. Sharpe,et al.  Optimal Pesticide Application for Controlling the Boll Weevil on Cotton , 1978 .

[29]  Yongzhen Pei,et al.  Species extinction and permanence in a prey–predator model with two-type functional responses and impulsive biological control , 2008 .