American Institute of Mathematical Sciences

doi: 10.3934/dcdsb.2021160
Online First

Online First articles are published articles within a journal that have not yet been assigned to a formal issue. This means they do not yet have a volume number, issue number, or page numbers assigned to them, however, they can still be found and cited using their DOI (Digital Object Identifier). Online First publication benefits the research community by making new scientific discoveries known as quickly as possible.

Readers can access Online First articles via the “Online First” tab for the selected journal.

The qualitative behavior of a plankton-fish interaction model with food limited growth rate and non-constant fish harvesting

 Department of Mathematics, Institute of Science, Banaras Hindu University, Varanasi-221 005, India

*Corresponding author: R. P. Gupta; e-mail: ravipguptaiitk@gmail.com, ravigupta@bhu.ac.in

Received  October 2020 Revised  March 2021 Early access June 2021

Fund Project: The work of the third author (Shivam Saxena) is supported by Council of Scientific and Industrial Research (No.09/013(0721)/2017)

The manuscript aims to investigate the qualitative analysis of a plankton-fish interaction with food limited growth rate of plankton population and non-constant harvesting of fish population. The ecological feasibility of population densities of both plankton and fish in terms of positivity and boundedness of solutions is shown. The conditions for the existence of various equilibrium points and their stability are derived thoroughly. This study mainly focuses on how the harvesting affects equilibrium points, their stability, periodic solutions and bifurcations in the proposed system. It is shown that the system exhibits saddle-node bifurcation in the form of a collision of two interior equilibrium points. Existence conditions for the occurrence of Hopf-bifurcation around interior equilibrium points are discussed. Lyapunov coefficients are examined to check the stability properties of these periodic solutions. We have also plotted the bifurcation diagrams for saddle-node, transcritical and Hopf bifurcations. A detailed algorithm for the occurrence of Bogdanov-Takens bifurcation is derived and finally some numerical simulations are also carried out to validate the theoretical results. This work suggests that the harvesting of fish population can change the dynamics of the system, which may be useful for the ecological management.

Citation: R. P. Gupta, Shristi Tiwari, Shivam Saxena. The qualitative behavior of a plankton-fish interaction model with food limited growth rate and non-constant fish harvesting. Discrete & Continuous Dynamical Systems - B, doi: 10.3934/dcdsb.2021160
References:

show all references

References:
Zero-growth isoclines along with the interior equilibrium states of the system (6) are depicted in this figure. Here, red and green color curves represent the plankton and fish nullclines respectively. For the values of the parameters $r = 0.1$, $k = 1$, $a = 0.3$, $m = 0.03$, $d = 0.07$, $n = 0.23$, $c = 0.4$, we get two equilibrium (A), then the two unique equilibrium collide to reach unique equilibrium (B) and then no equilibrium (C), for $h = 0.09, 0.1116729526$ and $0.13$ respectively. If we take $r = 0.25$, $k = 10$, $a = 0.2$, $m = 0.5$, $d = 0.9$, $h = 0.2$, $c = 0.2$, $n = 0.2$, then the system has unique interior equilibrium point (D)
In (A) $L_{2*}$ is stable for $r = 0.1$, $k = 1$, $a = 0.3$, $m = 0.03$, $d = 0.07$, $h = 0.11$, $c = 0.4$ and $n = 0.25$, if we increase the bifurcation parameter $n$ to $n = n^{[H]} = 0.3068476027$, the system possesses a periodic solution about $L_{2*}$ which is represented in (B) and finally for $n = 0.307$ and keeping rest of the parameters same, the periodic solution collide with the saddle point $L_{1*}$ to give a homoclinic orbit about $L_{2*}$ which is shown in (C) where as $L_{1*}$ remains saddle in each case. Figure (D) represents the bifurcation diagram for $r = 0.1$, $k = 1$, $a = 0.3$, $m = 0.03$, $d = 0.07$, $h = 0.11$ and $c = 0.4$ with respect to the bifurcation parameter $n$
(A) Bifurcation diagram w.r.t. parameter $a$ in case of two interior equilibrium. (B) Bifurcation diagram w.r.t. parameter $a$ in case of one interior equilibrium
(A) For $r = 0.1$, $k = 1$, $a = 0.3$, $m = 0.03$, $n = 0.23$, $d = 0.07$, $c = 0.4$ and $h = 0.09$ there are two interior equilibrium points of the system (6). (B) The model (6) attains a unique instantaneous equilibrium point for $h = h^{[SN]} = 0.1116729526$ and keeping rest of the values same. (C) For $h = 0.13$ and keeping other parameter fixed the system (6) has no equilibrium point. (D) Represents bifurcation curves in $\lambda_1$$\lambda_2$-plane
(A) The components of interior equilibrium are plotted to show their stability. The red curves stand for stable branch and green curves stand for unstable branch. (B) This figure depicts the P-components of $L_1$ and $L_{2*}$ for $r = 2.0$, $k = 22$, $a = 0.01$, $m = 0.02$, $n = 0.1$, $d = 0.08$, $h = 0.05$, $c = 0.2$ as d varies. When $d<1.51$ the equilibrium $L_1$ is stable while the interior equilibrium $L_{2*}$ is unstable and when $d>1.51$, the interior equilibrium $L_{2*}$ is stable while the equilibrium $L_1$ is unstable
(A) For $r = 0.1$, $k = 1$, $a = 0.3$, $m = 0.03$, $d = 0.07$, $c = 0.4$, $h = 0.1116729526$, $n = 0.2313139014$ and $h = 0.13$ the system attains an unique instantaneous equilibrium point which is a cusp of co-dimension 2. (B) Here, the trace and determinant of the variational matrix at $(N_{2*},P_{2*})$ are plotted in green and red color where their intersection is $(n_0, h_0) = (0.2313139014, 0.1116729526)$
In (A) there is a stable point for $r = 0.25$, $k = 10$, $a = 0.2$, $m = 0.5$, $d = 0.9$, $h = 0.2$, $c = 0.2$ and $n = 0.194$, if we increase the bifurcation parameter $n$ to $n = n^{[H]} = 0.1952313043$, the system possesses a periodic solution which is shown in figure (B) and finally for $n = 0.23$ and keeping rest of the parameters same, the periodic solution collide with the saddle point to give a homoclinic orbit which is given in figure (C) and $L_{1*}$ remains saddle in each case. (D) represents the bifurcation diagram for $r = 0.25$, $k = 10$, $a = 0.2$, $m = 0.5$, $d = 0.9$, $h = 0.2$ and $c = 0.2$ with respect to the bifurcation parameter $n$
 [1] K. Q. Lan, C. R. Zhu. Phase portraits of predator--prey systems with harvesting rates. Discrete & Continuous Dynamical Systems, 2012, 32 (3) : 901-933. doi: 10.3934/dcds.2012.32.901 [2] C. R. Zhu, K. Q. Lan. Phase portraits, Hopf bifurcations and limit cycles of Leslie-Gower predator-prey systems with harvesting rates. Discrete & Continuous Dynamical Systems - B, 2010, 14 (1) : 289-306. doi: 10.3934/dcdsb.2010.14.289 [3] Kie Van Ivanky Saputra, Lennaert van Veen, Gilles Reinout Willem Quispel. The saddle-node-transcritical bifurcation in a population model with constant rate harvesting. Discrete & Continuous Dynamical Systems - B, 2010, 14 (1) : 233-250. doi: 10.3934/dcdsb.2010.14.233 [4] Yuzhou Tian, Yulin Zhao. Global phase portraits and bifurcation diagrams for reversible equivariant Hamiltonian systems of linear plus quartic homogeneous polynomials. Discrete & Continuous Dynamical Systems - B, 2021, 26 (6) : 2941-2956. doi: 10.3934/dcdsb.2020214 [5] Soliman A. A. Hamdallah, Sanyi Tang. Stability and bifurcation analysis of Filippov food chain system with food chain control strategy. Discrete & Continuous Dynamical Systems - B, 2020, 25 (5) : 1631-1647. doi: 10.3934/dcdsb.2019244 [6] Jaume Llibre, Marzieh Mousavi. Phase portraits of the Higgins–Selkov system. Discrete & Continuous Dynamical Systems - B, 2022, 27 (1) : 245-256. doi: 10.3934/dcdsb.2021039 [7] Shubo Zhao, Ping Liu, Mingchao Jiang. Stability and bifurcation analysis in a chemotaxis bistable growth system. Discrete & Continuous Dynamical Systems - S, 2017, 10 (5) : 1165-1174. doi: 10.3934/dcdss.2017063 [8] Kerioui Nadjah, Abdelouahab Mohammed Salah. Stability and Hopf bifurcation of the coexistence equilibrium for a differential-algebraic biological economic system with predator harvesting. Electronic Research Archive, 2021, 29 (1) : 1641-1660. doi: 10.3934/era.2020084 [9] Ming Liu, Dongpo Hu, Fanwei Meng. Stability and bifurcation analysis in a delay-induced predator-prey model with Michaelis-Menten type predator harvesting. Discrete & Continuous Dynamical Systems - S, 2021, 14 (9) : 3197-3222. doi: 10.3934/dcdss.2020259 [10] Shengqin Xu, Chuncheng Wang, Dejun Fan. Stability and bifurcation in an age-structured model with stocking rate and time delays. Discrete & Continuous Dynamical Systems - B, 2019, 24 (6) : 2535-2549. doi: 10.3934/dcdsb.2018264 [11] Patrick Martinez, Judith Vancostenoble. Lipschitz stability for the growth rate coefficients in a nonlinear Fisher-KPP equation. Discrete & Continuous Dynamical Systems - S, 2021, 14 (2) : 695-721. doi: 10.3934/dcdss.2020362 [12] Jackson Itikawa, Jaume Llibre. Global phase portraits of uniform isochronous centers with quartic homogeneous polynomial nonlinearities. Discrete & Continuous Dynamical Systems - B, 2016, 21 (1) : 121-131. doi: 10.3934/dcdsb.2016.21.121 [13] Hebai Chen, Xingwu Chen. Global phase portraits of a degenerate Bogdanov-Takens system with symmetry (Ⅱ). Discrete & Continuous Dynamical Systems - B, 2018, 23 (10) : 4141-4170. doi: 10.3934/dcdsb.2018130 [14] Philip Gerlee, Alexander R. A. Anderson. Diffusion-limited tumour growth: Simulations and analysis. Mathematical Biosciences & Engineering, 2010, 7 (2) : 385-400. doi: 10.3934/mbe.2010.7.385 [15] Yoichi Enatsu, Yukihiko Nakata. Stability and bifurcation analysis of epidemic models with saturated incidence rates: An application to a nonmonotone incidence rate. Mathematical Biosciences & Engineering, 2014, 11 (4) : 785-805. doi: 10.3934/mbe.2014.11.785 [16] Patrice Bertail, Stéphan Clémençon, Jessica Tressou. A storage model with random release rate for modeling exposure to food contaminants. Mathematical Biosciences & Engineering, 2008, 5 (1) : 35-60. doi: 10.3934/mbe.2008.5.35 [17] Peng Zhong, Suzanne Lenhart. Optimal control of integrodifference equations with growth-harvesting-dispersal order. Discrete & Continuous Dynamical Systems - B, 2012, 17 (6) : 2281-2298. doi: 10.3934/dcdsb.2012.17.2281 [18] Benjamin Wincure, Alejandro D. Rey. Growth regimes in phase ordering transformations. Discrete & Continuous Dynamical Systems - B, 2007, 8 (3) : 623-648. doi: 10.3934/dcdsb.2007.8.623 [19] Yoora Kim, Gang Uk Hwang, Hea Sook Park. Feedback limited opportunistic scheduling and admission control for ergodic rate guarantees over Nakagami-$m$ fading channels. Journal of Industrial & Management Optimization, 2009, 5 (3) : 553-567. doi: 10.3934/jimo.2009.5.553 [20] Jicai Huang, Sanhong Liu, Shigui Ruan, Xinan Zhang. Bogdanov-Takens bifurcation of codimension 3 in a predator-prey model with constant-yield predator harvesting. Communications on Pure & Applied Analysis, 2016, 15 (3) : 1041-1055. doi: 10.3934/cpaa.2016.15.1041

2020 Impact Factor: 1.327