Abstract

Predator interference refers to the competition effect on individuals of the same species through preventing or limiting access to a resource. In this paper, a novel predator–prey model that includes predator interference is considered. In the novel model, a high prey density is taken into consideration by investigating Crowley–Martin as a numerical response, while the growth rate of prey is logistically described and Holling type II is used as a functional response. In contrast to literature, this model has distinct and notable features: predator and prey species have an intraspecific competition, and this model has a numerical response that generalizes several of the commonly utilized depictions of predator interference. The boundedness of the model is proved. Predator interference effects are theoretically and numerically studied in terms of stability, coexistence, and extinction. Theoretical results show that the system has no effect on coexistence or extinction in these types of models. Numerical simulations present two different dynamics with identical results on the theoretical side of model stabilization for both dynamics. Moreover, the numerical simulations show that there is a change from periodic to steady state dynamics as the value of predator interference increases. This means that the coexistence probability of the model increases. The attained results are explained from both a mathematical and an ecological points of view.

1. Introduction

Recent studies have used differential equations to describe a variety of applications, such as [18]. Population models, a central topic in many sciences due to their significance, study population change over time. One of the first scientists who noticed geometric population growth was Malthus [9], whose model is known by his name and is called the exponential model. After that, Verhulst [10] formulated the logistic model for population growth, which is called the Verhulst model and takes the sigmoidal form. The predator–prey model is one of the most common types of population models used to describe multispecies interactions. Lotka [11] and Volterra [12] independently founded the Lotka–Volterra model, which is considered the fundamental building block of predator–prey interactions and consists of two nonlinear differential equations. In order to make the Lotka–Volterra model more realistic, many ecological processes have been taken into account.

Two of the main components are functional and numerical responses, which connect the equations of predator and prey [13]. There are different kinds of functional and numerical responses. They are classified depending on the presence of predators and prey densities in the formulation of functional and numerical responses [14, 15]. Holling type II is one of the most commonly used types of functional and numerical responses and has been widely used in many studies with different ideas, for example [1621]. Generally, the main findings of these studies addressed the dynamics of the Holling type II models only; they derived their conditions, obtained their theoretical results, and performed simulations that followed certain dynamics such as stability, coexistence, and extinction under some external effects such as allee and refuge effects.

Holling type II is classified as prey-dependent functional and numerical responses because it is subjected to prey density only [22]. However, the Holling type II type was developed by including terms in the denominator to represent predator interference, and thus these types become dependent on the density of prey and predator as in the Beddington–DeAngelis (BD) [23, 24] and Crowley–Martin (CM) [25] types, but the difference between them is that in the case of the Crowley–Martin type, predator interference is taken into account even in a high prey density, whereas in the type of Beddington–DeAngelis, it excludes predator interference in a high prey density. Recent studies, such as those [26, 27], have used the Crowley–Martin type to discuss the dynamic behaviors of prey and predator species. They obtained the stability, coexistence, and extinction conditions of these models.

Predator–prey systems with mixed functional and numerical responses have been used in a few studies to describe some kinds of environment, which study the dynamic behaviors of these models such as [2830]. The authors of [28] used the Holling type II and Beddington–DeAngelis type as mixed functional responses to investigate the stability and find out the persistence condition; while in [29], they used the Holling type IV and Holling type I as mixed functional responses to examine the local and global stability and obtain the persistence condition. Finally, in [30], the authors investigated the stability of the equilibrium points of a discrete Lotka–Volterra predator–prey model with mixed functional responses of Holling type I and II, and the model was studied in terms of its long-term behavior when each control parameter was changed.

However, Holling type II is widely used and well documented in many theoretical and experimental studies; Crowley–Martin was developed based on Holling type II. To the best of our knowledge, Holling type II and Crowley–Martin have not been used as mixed functional responses to investigate the dynamics. This motivates us to investigate the dynamics of the predator–prey model that incorporates mixed functional responses from Crowley–Martin and Holling type II, respectively. We also examine how predator inference affects the dynamic behaviors of this model.

Interference predator is a concept that describes the effect of competition on individuals of the same species through preventing or limiting access to a resource. It is considered one of the intraspecific competition mechanisms [31]. The dynamics of some mathematical models that contain intraspecific competition have been discussed in many research papers [2646], which plays a significant role in systems’ dynamics. Although in both laboratory and natural systems, predator interference has been observed quite frequently, its effects on predator–prey dynamics have only been studied in a limited number of studies that have concentrated on the use of Beddington–DeAngelis as functional and numerical responses [4754], but Beddington–DeAngelis excludes predator interference at high prey densities.

The main contribution of this research is the study of the interference of predators in a comprehensive approach, including a high prey density, which is taken into consideration through an investigation of Crowley–Martin as a numerical response, while the growth rate of prey is logistically described and Holling type II is used as a functional response. This model has remarkable features: the prey and predator species have an intraspecific competition aspect, and this model’s numerical response is considered the most generalized of the commonly utilized depictions of predator interference. This study examines interference predator effects on stability, coexistence, and extinction of this model.

This manuscript is organized as follows: Section 2 introduces the model’s formulation and the proof of boundedness. Then, Kolmogorov conditions are discussed in Section 3. In Section 4, stability analysis and interference predator effects are explained theoretically. Then, numerical simulations are executed to support the theoretical results in Section 5. Finally, conclusions and discussion are given in Section 6.

2. Formulation and Boundedness of the Model

In this section, we introduce a continuous-time model that describes the interactions of predator and prey species and describes the growth rate of prey logistically in the absence of predator using mixed functional and numerical responses. The model is constructed as follows:

With and that represent positive initial conditions. indicates the growth rate of prey and symbolizes the carrying capacity of the prey in the absence of predation. is the accessibility parameter of prey by a predator, is the death rate of the predator, is conversion rate into biomass of the predator, is the handling time and represents the interference predator.

Theorem 1. The solutions on of model (1) in is positive and bounded with any given positive initial values.

Proof. We assume with positive initial conditions is any solution of model (1), so there is a unique solution . We take and observe that where . We consider .
The derivative of the function is as follows:The denominator of equation (1b) can be written as follows:Equation (3) can be written as follows:After simplification of equation (5), then we get the following equation:The inequality (11) obtains the following equation:Consequently, and are ultimately bounded.
We obtain that model (1) has three non-negative equilibrium points that are , , and . and exist without conditions. The existence of coexistence point is gotten by obtaining the positive real root of the following system algebraic equations:To discuss the existence of the coexistence point , we introduce the following analysis:
Through equation (13a), we assume as follows:One gets the following:(i)When , then . As a result, has no conditions attached to it.(ii)When , then . As a result, has no conditions attached to it. is always less than zero because , and the parameters are positive, so .
According to the above analysis, isocline (14) passes through points and , and is a decreasing function of .
Through equation (13b), supposeOne gets the following:(i)When , then . It is clear that under the following condition (17):(ii)As a result, the slope of isocline (16) always remains positive and it passes through the point under condition (17).
It is concluded from the above analysis that the two isoclines (14) and (16) intersect at a unique point when

Proposition 1. System (1) has a unique positive equilibrium point if conditions (17) and (18) are satisfied.

3. Kolmogorov Conditions

In this section, Kolmogorov conditions are applied to ensure that there is a unique nontrivial equilibrium point for model (1). In addition, Kolmogorov conditions are used to obtain the persistence and extinction conditions and validate the model.

Thus, model (1) is written as follows:

The first set of Kolmogorov conditions states that the growth rate of prey and predators decreases when the number of prey is fixed and the number of predators increases. The conditions are expressed as follows:

and are negative since all the variables and parameters are positive, so the growth rate of prey and predators decreases.

The second condition is that the model has a carrying capacity. The condition is indicated as follows:

The third condition is that if the model has a small predator population, it will include a minimum prey value.

The fourth condition is that the prey coexists with predators if the following condition holds:

This implies that

But, if condition (25) is violated, then the predator becomes extinct,

From Kolmogorov analysis, we attain the following theorems:

Proposition 2. Under condition (25), prey coexists with predators.

Corollary 1. The predators will become extinct if condition (26) is satisfied.

Through conditions (25) and (26), we conclude the following corollary:

Corollary 2. Predator interference does not affect coexistence and extinction.

4. Stability Analysis

For population dynamics, the qualitative analysis of differential equations as stability plays an influential role in following the dynamic behaviors of populations [55]. The stability analysis of non-negative equilibrium points is discussed through the following theorems.

We get the Jacobian matrix of model (1) as follows:

The following theorems are obtained to study the stability of positive equilibrium points.

Theorem 2. in model (1) is a saddle point (i.e., unstable).

Proof. By substituting , the Jacobian matrix of the model (1) is as follows:It is noticed that . Therefore, is a saddle point.

Theorem 3. Suppose that holds, then the equilibrium point is stable.

Proof. By substituting , the Jacobian matrix of model (1) is as follows:It has been observed that , so and if and only if . Consequently, is stable ifWe select the values of the parameters and initial conditions, which are as follows:
, , , , and , to plot two of equilibrium points and , while the third equilibrium point, which represents the coexistence equilibrium point , will be considered in detailed in the next section. is a saddle point, as theoretically proved in Theorem 2 in all cases, it is shown in Figure 1 and Figure 2 that the stream of arrows is moving away from point . With the values of parameters, the equilibrium point is represented by , which has two cases as proved theoretically in Theorem 3, in Figure 1, it is observed that is stable, as it heads towards the point , because condition (30) is satisfied in Theorem 3. However, if the value of the death rate of the predator (s) is changed, then condition (30) is not satisfied and point becomes a saddle point (i.e., an unstable case) as shown in Figure 2, where the stream of arrows is moving away from the point .

Theorem 4. The coexistence equilibrium point is locally asymptotically stable if it holds the following conditions:

Proof. The Jacobian matrix of is as follows:The determinant of is as follows:By simplifying the previous determinant, we get the following equation:Thus, if and .
Then, we obtain a trace of as follows:Therefore, if and .
This implies that the coexistence equilibrium point is stable.

Theorem 5. Model (1) has a periodic solution in under the following condition:

Proof. Let,By multiplying the Dulac function by and , we get the following equation:By taking the divertive of according to we have the following equation:Also for we have the following equation:By taking the divertive of according to we have the following equation:As observed, there is a change in the signs ofBy Bendixson–Dulac criterion, model (1) has a periodic solution under condition (37).

5. Numerical Simulations

By selecting parameter values that are biologically feasible, numerical simulations are conducted to present the dynamics of model (1), which show two different dynamics, steady state and periodic dynamics, as obtained in Section 4. In both dynamics, the values of parameters are selected so that the predator and prey coexist, which satisfies Proposition 2. In addition, interference predator effects are explained in terms of the dynamics of predator–prey interactions. Time series and phase portraits with nullclines figures are used to explain these effects and show the dynamics.

The first case is steady state dynamics:

, , , , , , , and . However, the value of the interference predator changes for each figure.

Second case is periodic dynamics:

, , , , , , , , and . However, the value of the interference predator (ν) changes for each figure.

Figures 36 present steady state dynamics, and with increasing the predator interference value for each figure, it is noticed that the dynamic behavior is still in a steady state, but the densities of prey and predator are affected, where the density of prey increases and the density of predator decreases at high values of predator interference, as explained by the intersection point of nullclines. Coexistence and extinction are not affected by predator interference. This corresponds with theoretical results (i.e., Corollary 3.3). In addition, it is clear that nullclines (8) intersect uniquely; this shows that the coexistence equilibrium point is unique.

In the periodic dynamic (Figure 7), it is observed that the size of the oscillation is reduced, as shown in Figure 8, by increasing the value of predator interference. However, as the predator interference value increases, the dynamics change to a steady state, as shown in Figures 9 and 10. This makes the model more stable, in addition to changing the densities of predator and prey as in a steady state dynamic. The change from periodic to steady state dynamic involves the coexistence probability increasing and the extinction probability decreasing as the value of predator interference increases.

6. Discussion and Conclusions

In this paper, Holling type II and Crowley–Martin functional and numerical responses are combined to form a reasonable predator–prey model. This model has remarkable features; the prey and predator species have an intraspecific competition aspect, and this model’s numerical response is considered the most generalized of the commonly utilized depictions of predator interference. Interference predators are included as an influential feature in this model that may affect the dynamics of predator–prey interactions, so this study aims to present the dynamics of this model. In addition, Crowley–Martin was developed based on Holling type II, and it has been widely used in many experimental and theoretical studies in the literature. The boundedness of this model is proved by Theorem 1.

Through Kolmogorov conditions, we conclude some significant results as follows:(i)Based on Kolmogorov’s four conditions (2026), respectively, predator–prey model (1) is endorsed.(ii)Ensure that there is a unique nontrivial equilibrium point for model (1), as illustrated in Proposition 1.(iii)Attain the coexistence and extinction conditions, which are conditions (25) and (26).(iv)Corollary 3.4 is an important result, which shows that the interference predator has no effect on the coexistence and extinction in this kind of models.

Stability analysis of model (1) is introduced to obtain the stability conditions as shown through Theorems 24. According to Theorem 5, model (1) has a periodic solution, so model (1) has two different dynamics: steady state and periodic dynamics.

The numerical simulations present two different dynamics: steady state and periodic dynamics, and the effect of predator interference is shown on the dynamics of predator–prey interactions. The numerical outcomes of predator interference effects are summarized as follows:(i)In the steady state dynamics, there is an increase in the density of prey and a decrease in the density of predators, which can be interpreted as the presence of aggressive behavior by one or more individuals, which affects the rest of the predators [28]. Thus, this means that predator interference has no effect on coexistence or extinction, not even the stability of the system.(ii)In the periodic dynamic, the effect is clear through the changing of the dynamics by making the predator–prey system stabilize with increasing predator interference in this kind of models. This result is consistent with the literature such as [4751].(iii)In addition, a significant result is obtained through the change from periodic to steady state dynamics. This can be interpreted as a meaning that the probability of coexistence increases and the extinction probability decreases in the model when the value of predator interference is increased.

Appendix

Proof. As model (1) is Kolmogorov type that satisfies the Kolmogorov conditions, so Freedman [56] proved that these kinds of systems has a unique nontrivial equilibrium point as referred in proof of a Kolmogorov type theorem (Theorem 5.1), he used the figures and the Poincare–Bendixson theorem to explain his proof.Proposition 2 Under condition (25), prey coexists with predators.Proof. Through Kolmogorov type theorem (Theorem 5.1), which model (1) is Kolmogorov type that satisfies the Kolmogorov conditions, coexistence condition is one of these conditions as referred in proof of a Kolmogorov type theorem (Theorem 5.1), he used the figures and the Poincare–Bendixson theorem to explain his proof.Corollary 1 The predators will become extinct, if the condition (26) is satisfied.Proof. It can show through the analysis of Preposition 1 that the same condition (18) presents when and . Also, Through Kolmogorov type theorem (Theorem 5.1), which model (1) is Kolmogorov type that satisfies the Kolmogorov conditions, the extinction condition is one of these conditions as referred in proof of a Kolmogorov type theorem (Theorem 5.1).Corollary 2 Predator interference does not affect coexistence and extinction.Proof. Since there is no interference predator parameter () in the coexistence and extinction conditions (25) and (26), respectively, then complete the proof of this corollary.

Data Availability

All data generated or analyzed during this study are included in this article.

Conflicts of Interest

The author declares no conflicts of interest.

Acknowledgments

The author extends his appreciation to the deputyship for Research & Innovation, Ministry of Education in Saudi Arabia for funding this research work through the project number (IFP-2020-98).