Dynamic Analysis of a Predator-Prey Model with Holling-II Functional Response

Jiemin Mo¹, Wanying Li¹, Donghuan He¹, Songbo Wang^2*, Xiaoliang Zhou¹

¹School of Mathematics and Statistics, Lingnan Normal University, Zhanjiang, China.
²School of Business, Lingnan Normal University, Zhanjiang, China.

DOI: 10.4236/jamp.2023.1110188                

Abstract

A predator-prey model with linear capture term Holling-II functional response was studied by using differential equation theory. The existence and the stabilities of non-negative equilibrium points of the model were discussed. The results show that under certain limited conditions, these two groups can maintain a balanced position, which provides a theoretical reference for relevant departments to make decisions on ecological protection.

Keywords

Predator-Prey Model, Holling-II Functional Response, Equilibrium Point, Stability

1. Introduction

Predator-prey system refers to the relationship in which one species is preyed upon by another species in a community of organisms. This relationship is an important part of the biome and can affect the number, distribution and composition of species, thereby altering the stability and integrity of the environment.

The Holling-II model (also known as the Lotka-Volterra model or food chain model) is a classic dynamical model used to study predator-prey interactions. The model is based on the assumption that as the number of prey increases, so does the number of predators; however, as an increase in predators, the likelihood of prey being preyed on increases, leading to a decrease in the number of prey. In turn, the decrease in prey populations inhibits the growth of predator populations. From the assumption Holling-II model, it can be seen that the quantitative relationship between predator and prey is nonlinear. In recent years, many researchers have studied predator-prey dynamics model from different angles (see [1] - [11] ). In the literature [10] , the authors studied a Holling-II predatory functional response in a predator-prey model as follows

$\{\begin{array}{l}\frac{\text{d}x}{\text{d}t}=x\left(a-x-\frac{y}{1+x}\right),\hfill \\ \frac{\text{d}y}{\text{d}t}=y\left(\frac{bx}{1+x}-d_1-d_{2}y\right).\hfill \end{array}$ (1)

They investigated the stabilities of the equilibrium point and the existence of limit cycle. In literature [11] , the authors considered a Leslie-Gower predator-prey model with time delay and linear harvesting term

$\{\begin{array}{l}\frac{\text{d}{x}_1}{\text{d}t}=\left(1-x_1\right)\left(x_1-m\right)-q{x}_1{x}_2+h_1{x}_1,\hfill \\ \frac{\text{d}{x}_2}{\text{d}t}=S{x}_2\left(1-\frac{x_2\left(t-{\tau }_2\right)}{x_1\left(t-{\tau }_1\right)}\right)+h_2{x}_2\hfill \end{array}$ (2)

where $h_1{x}_1$ and $h_2{x}_2$ denote the linear capture terms caused by external interference, and $h_1<0$ and $h_2<0$ . They studied the positive steady-state properties of the model (2) and gave the conditions for Hopf bifurcation near the positive equilibrium point. Mainly inspired by literatures [10] and [11] , combining models (1) and (2), we propose the following predator-prey model with Holling-II type functional response and a linear capture term

$\{\begin{array}{l}\frac{\text{d}u}{\text{d}t}=u\left(a-u-\frac{v}{1+u}\right)+h_{1}u,\hfill \\ \frac{\text{d}v}{\text{d}t}=v\left(\frac{bu}{1+u}-d_1-d_{2}v\right)+h_{2}v\hfill \end{array}$ (3)

where u and v represent the densities of prey and predator populations at time t respectively, a, b, d₁ and d₂ are positive parameters, and h₁ and h₂ are negative parameters. We use differential equation theory to discuss the existence and the stabilities of non-negative equilibrium points of the model (3). The results show that under certain limited conditions, these two populations can maintain balance.

2. The Equilibria of the System Model

In view of the practical significance of the system, we only consider the case where u and v satisfy $E=\left\{\left(u,v\right)|u\ge 0,v\ge 0\right\}$ . Define

$\{\begin{array}{l}Q\left(u,v\right):=u\left(a-u-\frac{v}{1+u}\right)+h_{1}u,\hfill \\ R\left(u,v\right):=v\left(\frac{bu}{1+u}-d_1-d_{2}v\right)+h_{2}v.\hfill \end{array}$

Let

$Q\left(u,v\right)=0,R\left(u,v\right)=0$ (4)

Then, we have three equilibrium points:

$P_1\left(0,0\right)$ , $P_2\left(a+h_1,0\right)$ , $P_3\left(u^{*},v^{*}\right)$

where $u^{*},v^{*}$ satisfy

$\{\begin{array}{l}{l}_1:a-u^{*}-\frac{v^{*}}{1+u^{*}}+h_1=0,\hfill \\ l_2:\frac{b{u}^{*}}{1+u^{*}}-d_1-d_2{v}^{*}+h_2=0.\hfill \end{array}$

We easily know that $u^{*}>0$ and $v^{*}>0$ when

$0<\frac{d_1-h_2}{b-d_1+h_2}<a+h_1$ .

3. Stabilities of the Equilibrium Points

The linear approximation at any point $P\left(u_0,v_0\right)$ of the model (3) is as follows

$\{\begin{array}{l}\frac{\text{d}u}{\text{d}t}={Q^{\prime }}_u\left(u_0,v_0\right)\cdot \left(u-u_0\right)+{Q^{\prime }}_v\left(u_0,v_0\right)\cdot \left(v-v_0\right),\\ \frac{\text{d}v}{\text{d}t}={R^{\prime }}_u\left(u_0,v_0\right)\cdot \left(u-u_0\right)+{R^{\prime }}_v\left(u_0,v_0\right)\cdot \left(v-v_0\right)\end{array}$ (5)

where

${Q^{\prime }}_u\left(u_0,v_0\right)=a-u_0-\frac{v_0}{1+u_0}+u_0\left[-1+\frac{v_0}{{\left(1+u_0\right)}^2}\right]+h_1$ , ${Q^{\prime }}_v=-\frac{u_0}{1+u_0}$ ,

${R^{\prime }}_u=\frac{v_{0}b}{{\left(1+u_0\right)}^2}$ , ${R^{\prime }}_v\left(u_0,v_0\right)=\frac{b{u}_0}{1+u_0}-2d_2{v}_0-d_1+h_2$ .

3.1. Stability of Equilibrium Point $P_1\left(0,0\right)$

In this case, we have

${Q^{\prime }}_u\left(0,0\right)\cdot \left(u-0\right)=\left(a+h_1\right)u$ , ${Q^{\prime }}_v\left(0,0\right)\cdot \left(v-0\right)=0$ ,

${R^{\prime }}_u\left(0,0\right)\cdot \left(u-0\right)=0$ , ${R^{\prime }}_v\left(0,0\right)\cdot \left(v-0\right)=\left(-d_1+h_2\right)v$ .

So the linear approximation Equation (5) is as follows:

$\{\begin{array}{l}\frac{\text{d}u}{\text{d}t}=\left(a+h_1\right)u,\hfill \\ \frac{\text{d}v}{\text{d}t}=\left(h_2-d_1\right)v.\hfill \end{array}$ (6)

The coefficient matrix of Equation (5) is $B:=(\begin{array}{c}a+h_1\\ 0\end{array}\begin{array}{c}0\\ h_2-d_1\end{array})$ . Then the characteristicroots are ${\lambda }_1=a+h_1$ and ${\lambda }_2=h_2-d_1$ . Because ${\lambda }_2<0$ , the following conclusion is obvious.

Theorem 1. If $a+h_1>0$ , then $P_1\left(0,0\right)$ is an unstable saddle point. If $a+h_1<0$ , then $P_1\left(0,0\right)$ is an asymptotically stable node.

3.2. Stability of Equilibrium Point $P_2\left(a+h_1,0\right)$

In this case, we have

${Q^{\prime }}_u\left(a+h_1,0\right)=-\left(a+h_1\right)$ , ${Q^{\prime }}_v\left(a+h_1,0\right)=-\frac{a+h_1}{1+a+h_1}$ ,

${R^{\prime }}_u\left(a+h_1,0\right)=0$ , ${R^{\prime }}_v\left(a+h_1,0\right)=\frac{b\left(a+h_1\right)}{1+a+h_1}-d_1+h_2$ .

So the linear approximation Equation (5) is:

$\{\begin{array}{l}\frac{\text{d}u}{\text{d}t}=\left(-a-h_1\right)u-\frac{a+h_1}{1+a+h_1}\cdot v+{\left(a+h_1\right)}^2,\hfill \\ \frac{\text{d}v}{\text{d}t}=\left[\frac{b\left(a+h_1\right)}{1+a+h_1}-d_1+h_2\right]v.\hfill \end{array}$ (7)

The coefficient matrix of Equation (7) is

$B:=\left(\begin{array}{cc}-a-h_1& -\frac{a+h_1}{1+a+h_1}\\ 0& \frac{b\left(a+h_1\right)}{1+a+h_1}-d_1+h_2\end{array}\right).$

The characteristic roots are

${\lambda }_1=-a-h_1=-\left(a+h_1\right)<0$ , ${\lambda }_2=\frac{\left(a+h_1\right)\left(b-d_1+h_2\right)-\left(d_1-h_2\right)}{1+a+h_1}$ .

Therefore, the following conclusion is obtained:

Theorem 2. If $\left(a+h_1\right)\left(b-d_1+h_2\right)>\left(d_1-h_2\right)$ , then ${\lambda }_2>0$ , so $P_2\left(a+h_1,0\right)$ is an unstable saddle point. If $\left(a+h_1\right)\left(b-d_1+h_2\right)<\left(d_1-h_2\right)$ , then ${\lambda }_2<0$ , so $P_2\left(a+h_1,0\right)$ is an asymptotically stable node.

3.3. Stability of Plane Equilibrium Point $P_3\left(u^{\ast },v^{\ast }\right)$

From

$v^{*}=\left(1+u^{*}\right)\cdot \left(a-u^{*}+h_1\right)=\frac{b{u}^{*}}{d_2\left(1+u^{*}\right)}-\frac{d_1-h_2}{d_2}$ ,

we have that

${Q^{\prime }}_u\left(u^{*},v^{*}\right)=u^{*}\left(\frac{a-u^{*}+h_1}{1+u^{*}}-1\right)$ , ${Q^{\prime }}_v\left(u^{*},v^{*}\right)=-\frac{u^{*}}{1+u^{*}}$ ,

${R^{\prime }}_u\left(u^{*},v^{*}\right)=\frac{b\left(a-u^{*}+h_1\right)}{1+u^{*}}$ , ${R^{\prime }}_v\left(u^{*},v^{*}\right)=\frac{-b{u}^{*}}{1+u^{*}}+d_1-h_2$ .

Therefore, the coefficient matrix of Equation (5) for $P_3\left(u^{*},v^{*}\right)$ is

$B:=\left(\begin{array}{cc}{u}^{*}\left(\frac{a-u^{*}+h_1}{1+u^{*}}-1\right)& -\frac{u^{*}}{1+u^{*}}\\ \frac{b\left(a-u^{*}+h_1\right)}{1+u^{*}}& \frac{-b{u}^{*}}{1+u^{*}}+d_1-h_2\end{array}\right)$

and its characteristic equation is

$\left|B-\lambda E\right|=\left[u^{*}\left(\frac{a-u^{*}+h_1}{1+u^{*}}-1\right)-\lambda \right]\cdot \left(\frac{-b{u}^{*}}{1+u^{*}}+d_1-h_2-\lambda \right)+\frac{b{u}^{*}\left(a-u^{*}+h_1\right)}{{\left(1+u^{*}\right)}^2}=0.$

Let

$a_1=\frac{a-u^{*}+h_1}{1+u^{*}}$ , $b_1=\frac{b{u}^{*}}{1+u^{*}}$ , (8)

then we have

${\lambda }^2+\left(\left(1-a_1\right)u^{*}+b_1-d_1+h_2\right)\lambda +u^{*}\left(a_1-1\right)\left(-b_1+d_1-h_2\right)+a_1{b}_1=0$ . (9)

Let

$p=\left(1-a_1\right)u^{*}+b_1-d_1+h_2$ , $q=u^{*}\left(a_1-1\right)\cdot \left(-b_1+d_1-h_2\right)+a_1{b}_1$ . (10)

Then the characteristic roots of Equation (9) satisfy

${\lambda }_1+{\lambda }_2=-p$ , ${\lambda }_1{\lambda }_2=q$ .

From (8) and (10), we have

$\begin{array}{c}p=\left(1-a_1\right)u^{*}+b_1-d_1+h_2=\left(1-\frac{a-u^{*}+h_1}{1+u^{*}}\right)u^{*}+\frac{b{u}^{*}}{1+u^{*}}-d_1+h_2\\ =\frac{2u^{{*}^2}+\left(1-a-h_1+b-d_1+h_2\right)u^{*}-d_1+h_2}{1+u^{*}}.\end{array}$

Define

$f\left(s\right)=2s^2+\left(1-a-h_1+b-d_1+h_2\right)s-d_1+h_2$ ,

then, its discriminant is

${\Delta }_1={\left(1-a-h_1+b-d_1+h_2\right)}^2+8\left(d_1-h_2\right)>0$ .

Thus, equation $f\left(s\right)=0$ has two roots:

$s_{1,2}=\frac{-\left(1-a-h_1+b-d_1+h_2\right)\mp \sqrt{{\Delta }_1}}{4}$ .

Obviously, we have

$s_1<0$ , $s_2>0$ .

Therefore, we can get

(i) as $0<u^{*}<s_2$ , we have $f\left(u^{*}\right)<0$ , which implies $p<0$ ;

(ii) as $u^{*}>s_2$ , we have $f\left(u^{*}\right)>0$ , which implies $p>0$ ;

(iii) as $u^{*}=s_2$ , we have $f\left(u^{*}\right)=0$ , which implies $p=0$ .

From (8) and (10), we also have

$\begin{array}{l}q=\frac{u^{*}}{{\left(1+u^{*}\right)}^2}(-2\left(-b+d_1-h_2\right)u^{*}{}^2+\left[\left(a+h_1\right)\left(-b+d_1-h_2\right)-3\left(d_1-h_2\right)\right]u^{*}\\ +\left(a+h_1\right)\left(d_1-h_2+b\right)-d_1+h_2).\end{array}$

Define

$\begin{array}{c}g\left(s\right)=-2\left(-b+d_1-h_2\right)s^2+\left[\left(a+h_1\right)\left(-b+d_1-h_2\right)-3\left(d_1-h_2\right)\right]s\\ +\left(a+h_1\right)\left(d_1-h_2+b\right)-d_1+h_2.\end{array}$

then, its discriminant is

$\begin{array}{c}{\Delta }_2={\left[\left(a+h_1\right)\left(-b+d_1-h_2\right)-3\left(d_1-h_2\right)\right]}^2\\ +8\left(-b+d_1-h_2\right)\cdot \left[\left(a+h_1\right)\left(d_1-h_2+b\right)-d_1+h_2\right].\end{array}$

From condition of $0<\frac{d_1-h_2}{b-d_1+h_2}<a+h_1$ , we can know that

$-b+d_1-h_2<0$ , $\left(a+h_1\right)\left(b-d_1+h_2\right)>d_1-h_2$ .

Therefore, the following relationships between q and ${\Delta }_2$ hold,

(a) If ${\Delta }_2<0$ , then $g\left(s\right)>0$ for all $s\in R$ , and thus $q>0$ ;

(b) If ${\Delta }_2=0$ , then $g\left(s\right)=0$ has two equal roots:

$s_3=s_4=\frac{-\left[\left(a+h_1\right)\left(-b+d_1-h_2\right)-3\left(d_1-h_2\right)\right]}{-4\left(-b+d_1-h_2\right)}$ ,

and

$s_3=s_4>0$ .

Therefore, $g\left(s\right)>0$ for either $0<s<s_3$ or $s>s_3$ , and thus $q>0$ ;

(c) If ${\Delta }_2>0$ , Then $g\left(s\right)=0$ has two real roots:

$u_{5,6}=\frac{-\left[\left(a+h_1\right)\left(-b+d_1-h_2\right)-3\left(d_1-h_2\right)\right]\mp \sqrt{{\Delta }_2}}{-4\left(-b+d_1-h_2\right)}$

and

$s_5>0$ , $s_6>0$ .

Therefore, as $s_5<u^{*}<s_6$ , we have $g\left(u^{*}\right)<0$ , which implies $q<0$ ; as $0<u^{*}<s_5$ or $s_6<u^{*}$ , we have $g\left(u^{*}\right)>0$ , which implies $q>0$ .

We give some assumptions as the follows.

(H1) $\begin{array}{l}{\left(\left(a+h_1\right)\left(-b+d_1-h_2\right)-3\left(d_1-h_2\right)\right)}^2\\ >8\left(b-d_1+h_2\right)\left(\left(a+h_1\right)\left(d_1-h_2+b\right)-d_1+h_2\right)\end{array}$ ;

(H2) $\begin{array}{l}{\left(\left(a+h_1\right)\left(-b+d_1-h_2\right)-3\left(d_1-h_2\right)\right)}^2\\ =8\left(b-d_1+h_2\right)\left(\left(a+h_1\right)\left(d_1-h_2+b\right)-d_1+h_2\right)\end{array}$ ;

(H3) $\begin{array}{l}{\left(\left(a+h_1\right)\left(-b+d_1-h_2\right)-3\left(d_1-h_2\right)\right)}^2\\ <8\left(b-d_1+h_2\right)\left(\left(a+h_1\right)\left(d_1-h_2+b\right)-d_1+h_2\right)\end{array}$ .

Theorem 3. Suppose conditions (H1) and $s_5<u^{*}<s_6$ hold, then $P_3\left(u^{*},v^{*}\right)$ is an unstable saddle point.

Proof If condition (H1) holds, we know ${\Delta }_2>0$ . Combined with condition $s_5<u^{*}<s_6$ and above discussion, we get $q<0$ . It implies that Equation (9) has two real roots which have different signs, thus $P_3\left(u^{*},v^{*}\right)$ is an unstable saddle point.

Theorem 4. Suppose conditions (H1) holds, we have

(A) if $u^{*}<\min \left\{s_2,s_5\right\}$ or $s_6<u^{*}<s_2$ , then $P_3\left(u^{*},v^{*}\right)$ is an unstable node or unstable focus;

(B) if $s_2<u^{*}<s_5$ or $\max \left\{s_2,s_6\right\}<u^{*}$ , then $P_3\left(u^{*},v^{*}\right)$ is a stable node or stable focus;

(C) if $u^{*}=s_2$ , then $P_3\left(u^{*},v^{*}\right)$ is a center.

Theorem 5. Suppose conditions (H2) holds, we have

(A) if $u^{*}<s_2$ and $u^{*}\ne s_3$ , then $P_3\left(u^{*},v^{*}\right)$ is an unstable node or unstable focus;

(B) if $u^{*}>s_2$ and $u^{*}\ne s_3$ , then $P_3\left(u^{*},v^{*}\right)$ is a stable node or stable focus;

(C) if $u^{*}=s_2$ , then $P_3\left(u^{*},v^{*}\right)$ is a center.

Theorem 6. Suppose conditions (H3) holds, we have

(A) if $u^{*}<s_2$ , then $P_3\left(u^{*},v^{*}\right)$ is an unstable node or unstable focus;

(B) if $u^{*}>s_2$ , then $P_3\left(u^{*},v^{*}\right)$ is a stable node or stable focus;

(C) if $u^{*}=s_2$ , then $P_3\left(u^{*},v^{*}\right)$ is a center.

4. Conclusion

This paper investigates a Holling-II functional response predator-prey model with linear capture terms. The existence of non-negative equilibrium point of the model was discussed, and the stability of the equilibrium point was analyzed. It is clear from this paper that human capture intensity can maintain an equilibrium position for both populations under certain constraints. However, when the intensity of human capture exceeds a certain limit, the two populations will lose equilibrium, which will have a negative impact on species diversity and environmental stability. In other words, human capture behavior, within certain limits, can not only maintain human needs but also achieve sustainable development of ecological resources. This paper provides a theoretical reference for the decision-making of relevant departments.

Acknowledgements

This work has been supported by the NSF of China (Grant No. 11961021), the NSF of Guangdong Province (Grant No. 2022A1515010964, 2022A1515010193), the Innovation and Developing School Project (Grant No. 2019KZDXM032) and the Key Project of Science and Technology Innovation of Guangdong College Students (Grant No. pdjh2022b0320, pdjh2023b0325).

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

References

[1]	Huang, Y.L., Shi, W.Y., Wei, C.J. and Zhang, S.W. (2020) A Stochastic Predator-Prey Model with Holling-II Increasing Function in the Predator. Journal of Biological Dynamics, 15, 1-18. https://doi.org/10.1080/17513758.2020.1859146
[2]	Mishra, P. and Raw, S.N. (2019) Dynamical Complexities in a Predator-Prey System Involving Teams of Two Prey and One Predator. Journal of Applied Mathematics and Computing, 61, 1-24. https://doi.org/10.1007/s12190-018-01236-9
[3]	Wei, Z. (2021) Impact of the Fear Effect in a Prey-Predator Model Incorporating a Prey Refuge and Partial Closure for the Populations. Applied Mathematics A Journal of Chinese Universities (Ser. A), 36, 277-284. (In Chinese)
[4]	Zou, X., Lv, J. and Wu, Y. (2020) A Note on a Stochastic Holling-II Predator-Prey Model with a Prey Refuge. Journal of the Franklin Institute, 357, 4486-4502. https://doi.org/10.1016/j.jfranklin.2020.03.013
[5]	Liu, Q., Zu, L. and Jiang, D. (2016) Dynamics of Stochastic Predator-Prey Models with Holling-II Functional Response. Communications in Nonlinear Science and Numerical Simulation, 37, 62-76. https://doi.org/10.1016/j.cnsns.2016.01.005
[6]	Feng, G., Zhao, H. and Zhang, X. (2013) Qualitative Analysis and Optimal Harvesting Strategy of Predator-Prey Model with Holling II Functional Response. Mathematics in Practice and Theory, 43, 191-197. (In Chinese)
[7]	Wei, C., Luo, Z. and Du, M. (2009) Qualitative Analysis of Two Species Predator-prey Systems with Holling’s Type II Functional Response. Journal of Wenzhou University (Natural Science Edition), 30, 6-10. (In Chinese)
[8]	Duan, M. and Ma, J. (2023) Dynamical Behavior of Discrete Ratio-Dependent Predator-Prey System. Open Journal of Applied Sciences, 13, 396-413. https://doi.org/10.4236/ojapps.2023.133032
[9]	Wang, B., Shi, S. and Luo, P. (2021) Positive Periodic Solutions of Nonautonomous Delayed Predator-Prey System with Pulse Controls. Journal of Applied Mathematics and Physics, 9, 1452-1467. https://doi.org/10.4236/jamp.2021.97098
[10]	Wu, F. (2016) Discussion on the Stability of Singular Point and the Existence of Limit Cycle for a Class of Predator-Prey System with Holling II Functional Response. Journal of Taishan University, 38, 18-29. (In Chinese)
[11]	Wang, F., Lv, T. and Zhou, L. (2022) Hopf Bifurcation of Leslie-Gower Predator-Prey System with Time Delay and Linear Harvesting Term. Natural Science Journal of Harbin Normal University, 38, 1-9. (In Chinese)