Next Article in Journal
Domains of Quasi Attraction: Why Stable Processes Are Observed in Reality?
Previous Article in Journal
Approximation of Caputo Fractional Derivative and Numerical Solutions of Fractional Differential Equations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Stochastic Modeling of Three-Species Prey–Predator Model Driven by Lévy Jump with Mixed Holling-II and Beddington–DeAngelis Functional Responses

1
Department of Engineering Mathematics and Computer Science, National School of Applied Sciences, University Hassan First of Sttat, Berrechid 26103, Morocco
2
Department of Mathematics and Computer Sciences, Faculty of Science, Necmettin Erbakan University, Konya 42090, Türkiye
3
Centre for Environmental Mathematics, Faculty of Environment, Science and Economy, University of Exeter, Cornwall TR10 9FE, UK
4
Department of Mathematics, Bartın University, Bartın 74100, Türkiye
*
Author to whom correspondence should be addressed.
Fractal Fract. 2023, 7(10), 751; https://doi.org/10.3390/fractalfract7100751
Submission received: 16 July 2023 / Revised: 9 October 2023 / Accepted: 9 October 2023 / Published: 12 October 2023

Abstract

:
This study examines the dynamics of a stochastic prey–predator model using a functional response function driven by Lévy noise and a mixed Holling-II and Beddington–DeAngelis functional response. The proposed model presents a computational analysis between two prey and one predator population dynamics. First, we show that the suggested model admits a unique positive solution. Second, we prove the extinction of all the studied populations, the extinction of only the predator, and the persistence of all the considered populations under several sufficient conditions. Finally, a special Runge–Kutta method for the stochastic model is illustrated and implemented in order to show the behavior of the two prey and one predator subpopulations.

1. Introduction

Prey–predator dynamics remain a highly fruitful and attractive ecological subject [1,2,3,4]. The basic prey–predator model was proposed by Lotka and Volterra in 1956 [5]; this basic ecological problem was one of the earliest to model competition between prey and predator. One study devoted to the dynamics of fishing illustrated the mathematical and ecological properties of marine fish [6]. Another application of these mathematical tools sought to understand the interaction between wildebeest, zebra, and lion [7] in Kruger National Park in South-Africa. The natural or more realistic growth rate is a logistic suggested by Verhulst [8] for new recruitment of either prey, predators, or both. This logistic growth rate charmed several scientists in various disciplines [9,10,11,12,13,14]. Other studies have been devoted to the type of competition between prey and predators; for example, the basic prey–predator model uses a mass action computational rate [5]. In another study [15], Yavuz and Sene performed a stability analysis for the fractional predator–prey model with harvesting rate. In [16], the authors provided an analysis of the prey–predator model using the Crowley–Martin type functional response. Their main object was to study the stability behavior of the considered population. The Holling-type functional impact on the dynamics of the prey–predator relationship was studied in [17,18]. The starting predator–prey model is based on the assumption that the reaction is between one predator and one prey species. Recently, the authors of [19,20] considered two prey and one predator species under a mixed Holling-II and Beddington–DeAngelis functional response, as follows:
d X 1 d t ( t ) = X 1 ( t ) υ 1 ϱ 1 X 1 ( t ) y ( t ) ι 1 + X 1 ( t ) , d X 2 d t ( t ) = X 2 ( t ) υ 2 ϱ 2 X 2 ( t ) y ( t ) ι 2 + X 2 ( t ) + β y ( t ) , d y d t ( t ) = y ( t ) d ϱ 3 y ( t ) + ϑ 1 X 1 ( t ) ι 1 + X 1 ( t ) + ϑ 2 X 2 ( t ) ι 2 + X 2 ( t ) + β y ( t ) ,
where the two prey and one predator populations are presented by  X 1 X 2 , and y, respectively. The parameters of (1) and their meanings are provided in Table 1.
Deterministic modeling lacks one of the most important properties in prey and predator, that of randomness; for this reason, the stochastic model is preferred in modeling such problems. This is because the stochastic model takes into account the variance aspect, not just the mean, as in the deterministic case. In addition, the stochastic model may offer different predicted results for the same initial data. Certain stochastic prey–predator systems can be established by taking into consideration only the white noise for the studied population [21,22]. Recently, a Brownian random walk perturbation on model (1) was proposed in [20] by studying the persistence and extinction of the prey and predator population. A more basic modification to the modeling technique is to immediately take into account more realistic distributions in order to maintain the model’s straightforward analytic form. This naturally leads to Lévy jumps, a broader class of driving processes [23,24,25,26,27,28,29,30,31,32,33,34,35]. In fact, the dynamical model may encounter abrupt and significant disturbances in predator–prey development due to the intrinsic stochastic features of the related processes [36]. Here, we add the jumps to the model from [20] in response to the results of earlier research. Thus, the model becomes
d X 1 ( t ) = X 1 ( t ) υ 1 ϱ 1 X 1 ( t ) y ( t ) ι 1 + X 1 ( t ) d t + σ 1 X 1 ( t ) d W 1 ( t ) + U Ξ 1 ( u ) X 1 ( t ) N ˜ ( d t , d u ) , d X 2 ( t ) = X 2 ( t ) υ 2 ϱ 2 X 2 ( t ) y ( t ) ι 2 + X 2 ( t ) + β y ( t ) d t + σ 2 X 2 ( t ) d W 2 ( t ) + U Ξ 2 ( u ) X 2 ( t ) N ˜ ( d t , d u ) , d y ( t ) = y ( t ) d ϱ 3 y ( t ) + ϑ 1 X 1 ( t ) ι 1 + X 1 ( t ) + ϑ 2 X 2 ( t ) ι 2 + X 2 ( t ) + β y ( t ) d t + σ 3 y ( t ) d W 3 ( t ) + U Ξ 3 ( u ) y ( t ) N ˜ ( d t , d u ) .
where the left limits of  X 1 ( t ) X 2 ( t ) , and  y ( t )  are denoted by  x ( t ) y ( t ) , and  z ( t ) , respectively, the standard Brownian motion  W i ( t )  is defined on the space  Ω , F , ( F t ) t 0 , P  as a complete probability, and  ( F t ) t 0  is the filtration verifying the usual conditions. We define the finite stationary compensator  ν  on U as the measurable subset of the non-negative half-line; then,  N ( d t , d u )  is the Poisson counting  ν ( d u ) d t σ i  is the intensity of  W i ( t ) , and the jumping intensities are represented by  q i ( u )  for  i = 1 , 2 , 3 .
The remaining parts of this study are organized as follows. In Section 2, we provide several properties of the solution of the model (2). The stochastic extinction of both the prey and the predator is established in Section 3. In Section 4, we prove the persistence of the prey and the extinction of the predator. The stochastic persistence of both the predator and prey is illustrated in Section 5. Finally, Section 6 of this paper is devoted to numerical simulations that validate our theoretical findings.

2. The Well-Posedness of the Solution

The following theorem ensures the well-posedness of the solution of the model (2).
Theorem 1. 
For all initial conditions that are positive, model (2) admits a unique positive global solution a.s.
Proof. 
Using the local Lipschitz of the diffusion and drift, for all initial positive data the problem (2) admits a local unique solution  X 1 ( t ) , X 2 ( t ) , y ( t )  for  t [ 0 , T e ) , with  T e  representing the time of explosion.
Proving the existence globality of the solution amounts to showing that  T e =  a.s. For a given a large number  n 0 > 0 , for any integer  n n 0  the stopping time is defined as
T n = inf { t [ 0 , T e ) / X 1 ( t ) ( 1 n , n ) or X 2 ( t ) ( 1 n , n ) or y ( t ) ( 1 n , n ) } ,
and when  n , the number  T n  increases. Using  T = lim n T n , with  T T e  a.s., it is sufficient to prove that  T = , which is signified  T e = , and  X 1 ( t ) , X 2 ( t ) , y ( t ) R + 3  a.s. Suppose the contrary, i.e.,  T <  a.s. This means that there exists some  T > 0  and  0 < ϵ < 1  such that  P ( T T ) ϵ .
Considering the functional
V X 1 ( t ) , X 2 ( t ) , y ( t ) = X 1 1 log X 1 + X 2 1 log X 2 + y 1 log y ,
from Itô’s formula [37], we obtain
d V X 1 , X 2 , y = L V d t + σ 1 ( X 1 1 ) d W 1 + σ 2 ( X 2 1 ) d W 2 + σ 3 ( y 1 ) d W 3 + U Ξ 1 ( u ) X 1 log 1 + Ξ 1 ( u ) X 1 + U Ξ 2 ( u ) X 2 log 1 + Ξ 2 ( u ) X 2 + U Ξ 3 ( u ) y log 1 + Ξ 3 ( u ) y ,
with
L V = 1 1 X 1 X 1 ( t ) υ 1 ϱ 1 X 1 y ι 1 + X 1 + σ 1 2 2 + 1 1 X 2 X 2 ( t ) υ 2 ϱ 2 X 2 y ι 2 + X 2 + β y + σ 2 2 2 + 1 1 y y d ϱ 3 y + ϑ 1 X 1 ι 1 + X 1 + ϑ 2 X 2 ι 2 + X 2 + β y + ρ σ 3 2 2 + U Ξ 1 ( u ) log 1 + Ξ 1 ( u ) ν ( d u ) + U Ξ 2 ( u ) log 1 + Ξ 2 ( u ) ν ( d u ) + U Ξ 3 ( u ) log 1 + Ξ 3 ( u ) ν ( d u ) ,
then,
L V υ 1 X 1 ϱ 1 X 1 2 + ϱ 1 X 1 + σ 1 2 2 + υ 2 X 2 ϱ 1 X 2 2 + ϱ 2 X 2 + σ 2 2 2 + d + ( ϑ 1 + ϑ 2 ) y ϱ 3 y 2 + σ 2 3 2 + M , L V ( υ 1 + ϱ 1 ) X 1 ϱ 1 X 1 2 + σ 1 2 2 + ( υ 2 + ϱ 2 ) X 2 ϱ 2 X 2 2 + σ 2 2 2 + d + ( ϑ 1 + ϑ 2 ) y ϱ 3 y 2 + σ 3 2 2 + 3 C , L V ( υ 1 + ϱ 1 ) 2 4 ϱ 1 + ( υ 2 + ϱ 2 ) 2 4 ϱ 2 + d + ( ϑ 1 + ϑ 2 ) 2 4 ϱ 3 + σ 1 2 2 + σ 2 2 2 + σ 3 2 2 + M , L V M ,
with
M = 3 max { U Ξ 1 ( u ) log 1 + Ξ 1 ( u ) ν ( d u ) , U Ξ 2 ( u ) log 1 + Ξ 2 ( u ) ν ( d u ) , U Ξ 3 ( u ) log 1 + Ξ 3 ( u ) ν ( d u ) } .
M = ( υ 1 + ϱ 1 ) 2 4 ϱ 1 + ( υ 2 + ϱ 2 ) 2 4 ϱ 2 + d + ( ϑ 1 + ϑ 2 ) 2 4 ϱ 3 + σ 1 2 2 + σ 2 2 2 + σ 3 2 2 + M .
Integrating Equation (3) between 0 and  T n T ,  we obtain
0 E V X 1 ( T n T ) , X 2 ( T n T ) , y ( T n T ) V X 1 ( 0 ) , X 2 ( 0 ) , y ( 0 ) + M T .
Thus, if  n , then
> V X 1 ( 0 ) , X 2 ( 0 ) , y ( 0 ) + M C T = ,
which is in contradiction with the previous supposition. Thus,  T =  and the solution  ( X 1 ( t ) , X 2 ( t ) , y ( t ) ) is unique and global a.s. □

3. Stochastic Extinction

Now, we prove the extinction of the three subpopulations a.s. First, we define
m 1 = υ 1 + U log 1 + Ξ 1 ( u ) Ξ 1 ( u ) ν ( d u ) σ 1 2 2 , m 2 = υ 2 + U log 1 + Ξ 2 ( u ) Ξ 2 ( u ) ν ( d u ) σ 2 2 2 , m 3 = ϑ 1 ϱ 1 υ 1 + U log 1 + Ξ 1 ( u ) Ξ 1 ( u ) ν ( d u ) σ 1 2 2 + ϑ 2 ϱ 2 υ 2 + U log 1 + Ξ 2 ( u ) Ξ 2 ( u ) ν ( d u ) σ 2 2 2 d + U log 1 + Ξ 3 ( u ) Ξ 3 ( u ) ν ( d u ) σ 3 2 2 .
Theorem 2. 
If  max { m 1 , m 2 , m 3 } < 0 ,  then for any initial positive condition we have
lim t + < X 1 > t   = 0 , lim t + < X 2 > t   = 0 , lim t + < y > t   = 0 .
Proof. 
Considering
F X 1 = log X 1 ( t ) ,
from Itô’s formula we have
d F = L F d t + σ 1 d W 1 ( t ) + U log 1 + Ξ 1 ( u ) N ˜ ( d t , d u ) ,
with
L F = υ 1 ϱ 1 X 1 y ι 1 + X 1 σ 1 2 2 + U log 1 + Ξ 1 ( u ) Ξ 1 ( u ) ν ( d u ) ,
then,
L F υ 1 ϱ 1 X 1 σ 1 2 2 + U log 1 + Ξ 1 ( u ) Ξ 1 ( u ) ν ( d u ) ,
and,
log X 1 ( t ) log X 1 ( 0 ) t υ 1 σ 1 2 2 + U log 1 + Ξ 1 ( u ) Ξ 1 ( u ) ν ( d u ) ϱ 1 < X 1 > t + 0 t σ 1 d W 1 ( s ) t d s + 1 t 0 t U log 1 + Ξ 1 ( u ) N ˜ ( d s , d u ) d s ,
therefore,
< X 1 > t   1 ϱ 1 υ 1 + U log 1 + Ξ 1 ( u ) Ξ 1 ( u ) ν ( d u ) σ 1 2 2 + 1 ϱ 1 0 t σ 1 d W 1 ( s ) t d s + 1 t 0 t U log 1 + Ξ 1 ( u ) N ˜ ( d s , d u ) d s 1 ϱ 1 log X 1 ( t ) log X 1 ( 0 ) t .
Let us now take
M t = 0 t σ 1 d W 1 ( s ) ,
then,
lim sup t + < M t , M t > t = lim sup t + σ 1 2 t < .
Using Martingale’s theorem (the strong law of large numbers), we obtain
lim sup t + M t t = 0 .
Thus,
lim t + < X 1 > t   1 ϱ 1 υ 1 + U log 1 + Ξ 1 ( u ) Ξ 1 ( u ) ν ( d u ) σ 1 2 2 .
Following the same technique, we obtain
lim inf t + < X 2 > t   1 ϱ 2 ϱ 2 + U log 1 + Ξ 2 ( u ) Ξ 2 ( u ) ν ( d u ) σ 2 2 2 .
We next consider
V ( t ) = log ( y ) .
Itô’s formula implies that
d V = L V d t + σ 3 d W 3 ( t ) + U log 1 + Ξ 3 ( u ) N ˜ ( d t , d u ) ,
with
L V = d ϱ 3 z + ϑ 1 X 1 ι 1 + X 1 + ϑ 2 X 2 ι 2 + X 2 + β y σ 3 2 2 + U log 1 + Ξ 3 ( u ) Ξ 3 ( u ) ν ( d u ) d ϱ 3 y + ϑ 1 X 1 + ϑ 2 X 2 + U log 1 + Ξ 3 ( u ) Ξ 3 ( u ) ν ( d u ) σ 3 2 2 ,
then,
V ( t ) V ( 0 ) t ϑ 1 < X 1 > + ϑ 2 < X 2 > d ϱ 3 < y > t + U log 1 + Ξ 3 ( u ) Ξ 3 ( u ) ν ( d u ) σ 3 2 2 + 1 t 0 t σ 3 d W 3 ( t ) + U log 1 + Ξ 3 ( u ) N ˜ ( d t , d u ) ,
then,
< y > t   1 ϱ 3 ϑ 1 < X 1 > + ϑ 2 < X 2 > d σ 3 2 2 + U log 1 + Ξ 3 ( u ) Ξ 3 ( u ) ν ( d u ) 1 ϱ 3 V ( t ) V ( 0 ) t + 1 α 3 1 t 0 t σ 3 d W 3 ( t ) + U log 1 + Ξ 3 ( u ) N ˜ ( d t , d u ) .
Similar to  < X 1 > t , we have
lim t + < y > t   1 ϱ 3 ( ϑ 1 ϱ 1 υ 1 + U log 1 + Ξ 1 ( u ) Ξ 1 ( u ) ν ( d u ) σ 1 2 2 + ϑ 2 ϱ 2 υ 2 + U log 1 + Ξ 2 ( u ) Ξ 2 ( u ) ν ( d u ) σ 2 2 2 + U log 1 + Ξ 3 ( u ) Ξ 3 ( u ) ν ( d u ) d σ 3 2 2 ) .

4. Stochastic Extinction of Predator

In order to prove the extinction of the predator, we need to define
m 1 = m 1 2 ϱ 1 υ 1 , m 2 = m 2 2 ϱ 2 υ 2 , m 4 = ϑ 1 + ϑ 2 d + U log 1 + Ξ 3 ( u ) Ξ 3 ( u ) ν ( d u ) σ 3 2 2 .
Theorem 3. 
Suppose that  m 4 0  and  min { m 1 , m 2 } > 0 ; then,
lim t < y > t   = 0 ,
and
lim inf t < X 1 > t   m 1 a . s . lim inf t < X 2 > t   m 2 a . s .
Proof. 
The first Equation of (2) implies that
X 1 ( t ) X 1 ( 0 ) t = 1 t 0 t X 1 ( s ) υ 1 ϱ 1 X 1 ( s ) y ( s ) ι 1 + X 1 ( s ) d s + σ 1 x ( s ) d W 1 ( s ) + 1 t 0 t U Ξ 1 ( u ) X 1 ( s ) N ˜ ( d s , d u ) d s .
Let us consider the following:
F X 1 = log X 1 ( t ) ,
from It’ô’s formula, we have
d F = υ 1 ϱ 1 X 1 σ 1 2 2 + U log 1 + Ξ 1 ( u ) Ξ 1 ( u ) ν ( d u ) y ι 1 + X 1 d t + σ 1 d W 1 ( t ) + U log 1 + Ξ 1 ( u ) N ˜ ( d t , d u ) ,
then,
X 1 ( t ) X 1 ( 0 ) t + F ( t ) F ( 0 ) t υ 1 < X 1 > t 2 ϱ 1 < X 1 > t 1 t 0 t ( X 1 + 1 ) y ι 1 + X 1 d s + υ 1 σ 1 2 2 + 0 t U log 1 + Ξ 1 ( u ) Ξ 1 ( u ) ν ( d u ) d s + 1 t 0 t σ 1 ( 1 + X 1 ( s ) ) d W 1 ( s ) + 1 t 0 t U log 1 + Ξ 1 ( u ) + Ξ 1 ( u ) X 1 ( s ) N ˜ ( d s , d u ) d s
and we know that
lim t + 1 t 0 t ( X 1 + 1 ) y ι 1 + X 1 d s lim t + 1 t 0 t ( X 1 + 1 ) ι 1 + X 1 2 d s 1 / 2 × 0 t y 2 d s 1 / 2 = 0 .
Thus,
lim inf t + < X 1 > t   υ 1 + U log 1 + Ξ 1 ( u ) Ξ 1 ( u ) ν ( d u ) σ 1 2 2 2 ϱ 1 υ 1 .
Using the same technique, we can show that
lim inf t + < X 2 > t   υ 2 + U log 1 + Ξ 2 ( u ) Ξ 2 ( u ) ν ( d u ) σ 2 2 2 2 ϱ 2 υ 2 .
Let us now take
V ( t ) = log ( y ) ,
from Itô’s formula, we have
d V = L V d t + σ 3 d W 3 ( t ) + U log 1 + Ξ 3 ( u ) N ˜ ( d t , d u ) ,
with
L V = d ϱ 3 y + ϑ 1 X 1 ι 1 + X 1 + ϑ 2 X 2 ι 2 + X 2 + β y σ 3 2 2 + U log 1 + Ξ 3 ( u ) Ξ 3 ( u ) ν ( d u ) d ϱ 3 y + υ 1 + υ 2 σ 3 2 2 ,
therefore,
V ( t ) V ( 0 ) t ϑ 1 + ϑ 2 d ϱ 3 < y > t + U log 1 + Ξ 3 ( u ) Ξ 3 ( u ) ν ( d u ) σ 3 2 2 + 1 t 0 t σ 3 d W 3 ( t ) + U log 1 + Ξ 3 ( u ) N ˜ ( d t , d u ) ,
then,
< y > t   1 ϱ 3 ϑ 1 + ϑ 2 d + U log 1 + Ξ 3 ( u ) Ξ 3 ( u ) ν ( d u ) σ 3 2 2 1 ϱ 3 V ( t ) V ( 0 ) t + 1 ϱ 3 1 t 0 t σ 3 d W 3 ( t ) + U log 1 + Ξ 3 ( u ) N ˜ ( d t , d u ) ,
and
lim t + < y > t   1 ϱ 3 ϑ 1 + ϑ 2 d + U log 1 + Ξ 3 ( u ) Ξ 3 ( u ) ν ( d u ) σ 3 2 2 .

5. Stochastic Persistence

In this section, we show that all three subpopulations persist.
Theorem 4. 
For  m 1 > 0 m 2 > 0 , and
1 ϱ 3 ( υ 1 ϑ 1 + ϑ 1 U log 1 + Ξ 1 ( u ) Ξ 1 ( u ) ν ( d u ) ϑ 1 σ 1 2 2 ( ι 1 + υ 1 2 ϱ 1 ) ( υ 1 2 ϱ 1 ) + υ 2 ϑ 2 + ϑ 2 U log 1 + Ξ 2 ( u ) Ξ 2 ( u ) ν ( d u ) ϑ 2 σ 2 2 2 ( ι 2 + υ 2 2 ϱ 2 + β υ 2 2 ϱ 2 ) ( υ 2 2 ϱ 2 ) d + U log 1 + Ξ 3 ( u ) Ξ 3 ( u ) ν ( d u ) σ 3 2 2 ) > 0 ,
which means that all of the subpopulations persist in the mean. In addition,
lim inf t + < X 1 > t   υ 1 + U log 1 + Ξ 1 ( u ) Ξ 1 ( u ) ν ( d u ) σ 1 2 2 2 ϱ 1 υ 1 a . s . lim inf t + < X 2 > t   υ 2 + U log 1 + Ξ 2 ( u ) Ξ 2 ( u ) ν ( d u ) σ 2 2 2 2 ϱ 2 υ 2 a . s . lim inf t + < y > t   1 ϱ 3 ( υ 1 ϑ 1 + ϑ 1 U log 1 + Ξ 1 ( u ) Ξ 1 ( u ) ν ( d u ) ϑ 1 σ 1 2 2 ( ι 1 + υ 1 2 ϱ 1 ) ( υ 1 2 ϱ 1 ) + υ 2 ϑ 2 + ϑ 2 U log 1 + Ξ 2 ( u ) Ξ 2 ( u ) ν ( d u ) ϑ 2 σ 2 2 2 ( ι 2 + υ 2 2 ϱ 2 + β υ 2 2 ϱ 2 ) ( υ 2 2 ϱ 2 ) d + U log 1 + Ξ 3 ( u ) Ξ 3 ( u ) ν ( d u ) σ 3 2 2 ) a . s .
Proof. 
The first Equation of (2) implies that
X 1 ( t ) X 1 ( 0 ) t = 1 t 0 t X 1 ( s ) υ 1 ϱ 1 X 1 ( s ) y ( s ) ι 1 + X 1 ( s ) d s + σ 1 x ( s ) d W 1 ( s ) + 1 t 0 t U Ξ 1 ( u ) X 1 ( s ) N ˜ ( d s , d u ) d s .
Let us consider the following:
F X 1 = log X 1 ( t ) ,
from Itô’s formula, we have
d F = υ 1 ϱ 1 X 1 σ 1 2 2 + U log 1 + Ξ 1 ( u ) Ξ 1 ( u ) ν ( d u ) y ι 1 + X 1 d t + σ 1 d W 1 ( t ) + U log 1 + Ξ 1 ( u ) N ˜ ( d t , d u ) ,
then,
X 1 ( t ) X 1 ( 0 ) t + F ( t ) F ( 0 ) t υ 1 < X 1 > t 2 ϱ 1 < X 1 > t 1 t 0 t ( X 1 + 1 ) y ι 1 + X 1 d s + υ 1 σ 1 2 2 + 0 t U log 1 + Ξ 1 ( u ) Ξ 1 ( u ) ν ( d u ) d s + 1 t 0 t σ 1 ( 1 + X 1 ( s ) ) d W 1 ( s ) + 1 t 0 t U log 1 + Ξ 1 ( u ) + Ξ 1 ( u ) X 1 ( s ) N ˜ ( d s , d u ) d s .
We know that
lim t + 1 t 0 t ( X 1 + 1 ) y ι 1 + X 1 d s lim t + 1 t 0 t ( X 1 + 1 ) ι 1 + X 1 2 d s 1 / 2 × 0 t y 2 d s 1 / 2 = 0 .
Thus,
lim inf t + < X 1 > t   υ 1 + U log 1 + Ξ 1 ( u ) Ξ 1 ( u ) ν ( d u ) σ 1 2 2 2 ϱ 1 υ 1 .
Using the same technique, we can show that
lim inf t + < X 2 > t   υ 2 + U log 1 + Ξ 2 ( u ) Ξ 2 ( u ) ν ( d u ) σ 2 2 2 2 ϱ 2 υ 2 .
Now, we need to take into account that
V ( t ) = log ( y ) ,
where Itô’s formula implies that
d V = L V d t + σ 3 d W 3 ( t ) + U log 1 + Ξ 3 ( u ) N ˜ ( d t , d u ) ,
with
L V = d ϱ 3 y + ϑ 1 X 1 ι 1 + X 1 + ϑ 2 X 2 ι 2 + X 2 + β y σ 3 2 2 + U log 1 + Ξ 3 ( u ) Ξ 3 ( u ) ν ( d u ) ,
and,
V ( t ) V ( 0 ) t ϑ 1 ι 1 + υ 1 2 ϱ 1 υ 1 + U log 1 + Ξ 1 ( u ) Ξ 1 ( u ) ν ( d u ) σ 1 2 2 υ 1 2 ϱ 1 + ϑ 1 ι 2 + υ 2 2 ϱ 2 + β υ 2 2 ϱ 2 υ 2 + U log 1 + Ξ 2 ( u ) Ξ 2 ( u ) ν ( d u ) σ 2 2 2 υ 2 2 ϱ 2 d ϱ 3 < y > t + U log 1 + Ξ 3 ( u ) Ξ 3 ( u ) ν ( d u ) σ 3 2 2 + 1 t 0 t σ 3 d W 3 ( t ) + U log 1 + Ξ 3 ( u ) N ˜ ( d t , d u ) ,
therefore,
< y > t   1 ϱ 3 ϑ 1 υ 1 + ϑ 1 U log 1 + Ξ 1 ( u ) Ξ 1 ( u ) ν ( d u ) ϑ 1 σ 1 2 2 ( ι 1 + N ) ( υ 1 2 ϱ 1 ) + 1 ϱ 3 ϑ 2 υ 2 + ϑ 2 U log 1 + Ξ 2 ( u ) Ξ 2 ( u ) ν ( d u ) ϑ 2 σ 2 2 2 ( ι 2 + υ 2 2 ϱ 2 + β υ 2 2 ϱ 2 ) ( υ 2 2 ϱ 2 ) 1 ϱ 3 d + U log 1 + Ξ 3 ( u ) Ξ 3 ( u ) ν ( d u ) σ 3 2 2 V ( t ) V ( 0 ) t + 1 ϱ 3 1 t 0 t σ 3 d W 3 ( t ) + U log 1 + Ξ 3 ( u ) N ˜ ( d t , d u ) ,
and
lim inf t + < y > t 1 ϱ 3 ( υ 1 ϑ 1 + ϑ 1 U log 1 + Ξ 1 ( u ) Ξ 1 ( u ) ν ( d u ) ϑ 1 σ 1 2 2 ( ι 1 + υ 1 2 ϱ 1 ) ( υ 1 2 ϱ 1 ) + υ 2 ϑ 2 + ϑ 2 U log 1 + Ξ 2 ( u ) Ξ 2 ( u ) ν ( d u ) ϑ 2 σ 2 2 2 ( ι 2 + υ 2 2 ϱ 2 + β υ 2 2 ϱ 2 ) ( υ 2 2 ϱ 2 ) d + U log 1 + Ξ 3 ( u ) Ξ 3 ( u ) ν ( d u ) σ 3 2 2 ) .

6. Numerical Analysis

In this section, we illustrate the algorithm used in the numerical analysis of the considered model (2). First, we use the following equation:
d Z ( t ) = G t , Z ( t ) d t + σ t , Z ( t ) d W t + U Π ( t , u ) H t , Z ( t ) N ˇ ( d t , d u ) .
Thus,
Z ( t ) = Z ( 0 ) + 0 t G ξ , Z ( ξ ) d ξ + 0 t σ ξ , Z ( ξ ) d W ξ Part I + H t , Y ( t ) 0 t U Π ( ξ , u ) N ˇ ( d ξ , d u ) Part II .
Using the Runge–Kutta method, we can provide an approximation of part I in (8):
Z j + 1 = Z j + 1 6 M 1 + 2 M 2 + 2 M 3 + M 4 ,
with
M 1 = G ( t j , Z j ) h t + σ ( t j , Z j ) ( W j + 1 W j ) , M 2 = G ( t j + h t 2 t , Z j + M 1 2 ) h t + σ ( t j + h t 2 , Z j + M 1 2 ) ( W j + 1 W j ) , M 3 = G ( t j + h t 2 , Z j + M 2 2 ) h t + σ ( t j + h t 2 , Z j + M 2 2 ) ( W j + 1 W j ) , M 4 = G ( t j + h t , Z j + M 3 ) h t + σ ( t j + h t , Z j + M 3 ) ( W j + 1 W j ) .
We consider the infinitesimal interval  [ T j , T j + 1 ) ( t j , t j + 1 ) ; then, for the second part, we have the following:
With no jumping on  [ T j , T j + 1 ) , we obtain
H t , Z ( t ) T j T j + 1 U Π ( ξ , u ) N ˇ ( d ξ , d u ) = 0 .
With one jump at the point  t i [ T i , T i + 1 ) ,  we obtain
H t , Z ( t ) T j T j + 1 U Π ( u ) N ˇ ( d ξ , d u ) = H t j , Z ( t j ) Π t j , ς ( t j ) ,
and,
H t , Z ( t ) 0 t U Π ( ξ , u ) N ˇ ( d ξ , d u ) = j = 0 j = n H t i , Z ( t j ) Π t j , ς ( t j ) .
Therefore,
Z j + 1 = Z j + 1 6 M 1 + 2 M 2 + 2 M 3 + M 4
is the Runge–Kutta version of Equation (8), with
M 1 = G ( t j , Z j ) h t + σ ( t j , Z j ) ( W j + 1 W j ) + H ( t j , Z ( t j ) ) Π ( t j , ς ( t j ) ) , M 2 = G ( t j + h t 2 , Z j + M 1 2 ) h t + σ ( t j + h t 2 , Z j + M 1 2 ) ( W j + 1 W j ) + H ( t j + h t 2 , Z ( t j ) + M 1 2 ) Π ( t j + h t 2 , ς ( t j + h t 2 ) ) , M 3 = G ( t j + h t 2 , Z j + M 2 2 ) h t + σ ( t j + h t 2 , Y j + M 2 2 ) ( W j + 1 W j ) + H ( t J + h t 2 , Z ( t j ) + M 2 2 ) Π ( t j + h t 2 , ς ( t j + h t 2 ) ) , M 4 = G ( t j + h t , Z j + M 3 ) h t + σ ( t j + h t , Z j + M 3 ) ( W j + 1 W j ) + f ( t j + h t , Z ( t j ) + M 3 ) Π ( t j + h t , ς ( t j + h t ) ) .
The different values used for the numerical results are provided in Table 2.
The dynamics of the two prey subpopulations  X 1  and  X 2  and the predator subpopulation y can be seen in Figure 1, where the extinction of both prey species and the predator species can be observed.
Figure 2 illustrates the behavior of  X 1 X 2  and y; it can be observed that the prey subpopulations remain strictly positive and the predator subpopulation vanishes, indicates the extinction of the predator and is opposite to the scenario in which the predator and prey populations persist.
The behavior of  X 1 X 2  and y is illustrated in Figure 3, observing that all the considered populations persist, which coincides with our theoretical result in the case of persistence.
In Figure 4, we show a comparison between the stochastic model (2) and deterministic model (1) in the persistence case. It can be observed that the curve of the stochastic case has perturbations around the deterministic curve, which proves that the stochastic model can more efficiently describe the dynamics of the prey–predator model.
Figure 5 shows a comparison of the stochastic model (2) in the persistence case with and without Lévy jumps. It is notable that the curve of the stochastic model with jumping has a number of jumped-up perturbations, while the curve without Lévy jumps has normal perturbations. This observation demonstrates that certain cases can lead to Lévy jumps in the behavior of the prey–predator model.

7. Conclusions and Discussion

This paper has investigated a stochastic prey–predator model driven by Lévy noise in the presence of two prey subpopulations and one predator subpopulation. The model under consideration uses two forms of functional responses, namely, Holling-II  y X 2 ι 2 + X 2 + β y  and Beddington–DeAngelis  y X 1 ι 1 + X 1 , due to the saturation behavior. These two mixed functionals can better describe the reality of the computation between the prey and predator subpopulations. First, we have established the well-posedness of the solution of the model (2). An analysis of persistence and extinction led us to three possible cases: the first is devoted to the extinction of all the considered population when  max { m 1 , m 2 , m 3 } 0 ; the second case is dedicated to the persistence of the two prey species and the extinction of the predator species when  min { m 1 , m 2 } 0  and  m 4 0 ; and in the last case, we show the persistence of all the studied populations under several conditions. In our numerical analysis, we illustrate a special Runge–Kutta method for the stochastic model under Brownian and Lévy jump perturbations as implemented in Matlab R2019b software.
In this way, we considered the three cases with respect to extinction (persistence of the two prey species, extinction of the predator, and persistence of all three species) in simulations in order to validate our theoretical findings. Finally, we present several numerical comparisons. The first is between the deterministic and stochastic cases, which shows the usefulness of considering the stochastic model. The second comparison is between the cases with and without Lévy jumps, which proves that certain cases naturally lead to Lévy jumping in the behavior of the prey–predator model. Future efforts should address remaining issues; one of these is how the presented predator-prey model spreads via Levy jumps. In addition, we intend to study the chaotic behavior of the stochastic predator–prey system. Moreover, it is possible to model the prey–predator relationship using the space–time spectral order Sinc-collocation method (see [36] for example) and to illustrate a new modified stochastic predictor–corrector compact difference method, such as the one suggested in [38].

Author Contributions

Conceptualization, M.Y. (Mehmet Yavuz); Methodology, J.D. and M.Y. (Mustafa Yıldız); Software, J.D.; Validation, M.Y. (Mustafa Yıldız); Formal analysis, M.Y. (Mehmet Yavuz); Investigation, M.Y. (Mustafa Yıldız); Resources, M.Y. (Mehmet Yavuz); Data curation, J.D. and M.Y. (Mehmet Yavuz); Writing—original draft, J.D.; Writing—review & editing, M.Y. (Mehmet Yavuz) and M.Y. (Mustafa Yıldız). All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Our manuscript has no associated data.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Naik, P.A.; Eskandari, Z.; Shahkari, H.E.; Owolabi, K.M. Bifurcation analysis of a discrete-time prey–predator model. Bull. Biomath. 2023, 1, 111–123. [Google Scholar] [CrossRef]
  2. Chatterjee, A.; Pal, S. A predator–prey model for the optimal control of fish harvesting through the imposition of a tax. Int. J. Optim. Control Theor. Appl. 2023, 13, 68–80. [Google Scholar] [CrossRef]
  3. Tripathi, J.P.; Tyagi, S.; Abbas, S. Global analysis of a delayed density dependent predator–prey model with Crowley-Martin functional response. Commun. Nonlinear Sci. Numer. Simul. 2016, 30, 45–69. [Google Scholar] [CrossRef]
  4. Naik, P.A.; Eskandari, Z.; Yavuz, M.; Zu, J. Complex dynamics of a discrete-time Bazykin–Berezovskaya prey-predator model with a strong Allee effect. J. Comput. Appl. Math. 2022, 413, 114401. [Google Scholar] [CrossRef]
  5. Lotka, A. Elements of Mathematical Biology; Dover Publications: New York, NY, USA, 1956. [Google Scholar]
  6. Spencer, P.D.; Collie, J.S. A simple predator–prey model of exploited marine fish populations incorporating alternative prey. ICES J. Mar. Sci. 1996, 53, 615–628. [Google Scholar] [CrossRef]
  7. Fay, T.H.; Greeff, J.C. Lion, wildebeest and zebra: A predator–prey model. Ecol. Model. 2006, 196, 237–244. [Google Scholar] [CrossRef]
  8. Verhulst, P.F. Notice sur la loi que la population suit dans son accroissement. Corresp. MathéMatique Phys. 1838, 10, 113–121. [Google Scholar]
  9. Pearl, R.; Reed, L.J. The logistic curve and the census count of i930. Science 1930, 72, 399–401. [Google Scholar] [CrossRef] [PubMed]
  10. MacLean, M.; Willard, A. The logistic curve applied to Canada’s population. Can. J. Econ. Polit. Sci. 1937, 3, 241–248. [Google Scholar] [CrossRef]
  11. Kot, M. Elements of Mathematical Ecology; Cambridge University Press: Cambridge, UK, 2001. [Google Scholar]
  12. Kingsland, S. The refractory model: The logistic curve and the history of population ecology. Q. Rev. Biol. 1982, 57, 29–52. [Google Scholar] [CrossRef]
  13. Keshet, E.L. Mathematical Models in Biology; McGrawHill: New York, NY, USA, 1988. [Google Scholar]
  14. Ghosh, D.; Santra, P.K.; Mahapara, G.S. A three-component prey-predator system with interval number. Math. Model. Numer. Simul. Appl. 2023, 3, 1–16. [Google Scholar] [CrossRef]
  15. Yavuz, M.; Sene, N. Stability analysis and numerical computation of the fractional predator–prey model with the harvesting rate. Fractal Fract. 2020, 4, 35. [Google Scholar] [CrossRef]
  16. Meng, X.Y.; Huo, H.F.; Xiang, H.; Yin, Q.Y. Stability in a predator–prey model with Crowley-Martin function and stage structure for prey. Appl. Math. Comput. 2014, 232, 810–819. [Google Scholar] [CrossRef]
  17. Ko, W.; Ryu, K. Qualitative analysis of a predator—Prey model with Holling type II functional response incorporating a prey refuge. J. Differ. Equ. 2006, 231, 534–550. [Google Scholar] [CrossRef]
  18. Sugie, J.; Kohno, R.; Miyazaki, R. On a predator-prey system of Holling type. Proc. Am. Math. Soc. 1997, 125, 2041–2050. [Google Scholar] [CrossRef]
  19. Baek, H.; Dongseok, K. Dynamics of a predator-prey system with mixed functional responses. J. Appl. Math. 2014, 2014, 536019. [Google Scholar] [CrossRef]
  20. Zou, X.; Li, Q.; Cao, W.; Lv, J. Thresholds and critical states for a stochastic predator–prey model with mixed functional responses. Math. Comput. Simul. 2023, 206, 780–795. [Google Scholar] [CrossRef]
  21. Liu, Q.; Jiang, D.; Hayat, T.; Alsaedi, A. Dynamics of a stochastic predator–prey model with stage structure for predator and Holling type II functional response. J. Nonlinear Sci. 2018, 28, 1151–1187. [Google Scholar] [CrossRef]
  22. Rihan, F.A.; Alsakaji, H.J. Stochastic delay differential equations of three-species prey-predator system with cooperation among prey species. Discret. Contin. Dyn.-Syst. 2020, 15, 245–263. [Google Scholar] [CrossRef]
  23. Danane, J.; Allali, K.; Hammouch, Z.; Nisar, K.S. Mathematical analysis and simulation of a stochastic COVID-19 Lévy jump model with isolation strategy. Results Phys. 2021, 23, 103994. [Google Scholar] [CrossRef] [PubMed]
  24. Akdim, K.; Ez-zetouni, A.; Danane, J.; Allali, K. Stochastic viral infection model with lytic and nonlytic immune responses driven by Lévy noise. Phys. Stat. Mech. Its Appl. 2020, 549, 124367. [Google Scholar] [CrossRef]
  25. Danane, J. Stochastic predator–prey Lévy jump model with Crowley–Martin functional response and stage structure. J. Appl. Math. Comput. 2021, 67, 41–67. [Google Scholar] [CrossRef]
  26. Zhao, S.; Song, M. A stochastic predator-prey system with stage structure for predator. Abstr. Appl. Anal. 2014, 2014, 518695. [Google Scholar] [CrossRef]
  27. Danane, J. Stochastic Capital–Labor Lévy Jump Model with the Precariat Labor Force. Math. Comput. Appl. 2022, 27, 93. [Google Scholar] [CrossRef]
  28. Choo, S.; Kim, Y.H. Global stability in stochastic difference equations for predator-prey models. J. Comput. Anal. Appl. 2017, 23, 3. [Google Scholar]
  29. Yang, L.; Zhong, S. Global stability of a stage-structured predator-prey model with stochastic perturbation. Discret. Dyn. Nat. Soc. 2014, 2014, 512817. [Google Scholar] [CrossRef]
  30. Zhao, Y.; Yuan, S. Stability in distribution of a stochastic hybrid competitive Lotka–Volterra model with Lévy jumps. Chaos Solitons Fractals 2016, 85, 98–109. [Google Scholar] [CrossRef]
  31. Liu, Q.; Chen, Q.M. Analysis of a stochastic delay predator-prey system with jumps in a polluted environment. Appl. Math. Comput. 2014, 242, 90–100. [Google Scholar] [CrossRef]
  32. Zhang, X.; Li, W.; Liu, M.; Wang, K. Dynamics of a stochastic Holling II one-predator two-prey system with jumps. Phys. Stat. Mech. Its Appl. 2015, 421, 571–582. [Google Scholar] [CrossRef]
  33. Liu, M.; Bai, C.; Deng, M.; Du, B. Analysis of stochastic two-prey one-predator model with Lévy jumps. Phys. A Stat. Mech. Its Appl. 2016, 445, 176–188. [Google Scholar] [CrossRef]
  34. Sabbar, Y. Asymptotic extinction and persistence of a perturbed epidemic model with different intervention measures and standard lévy jumps. Bull. Biomath. 2023, 1, 58–77. [Google Scholar] [CrossRef]
  35. Mouhcine, N.; Sabbar, Y.; Anwar, Z. Stability characterization of a fractional-order viral system with the non-cytolytic immune assumption. Math. Model. Numer. Simul. Appl. 2022, 2, 164–176. [Google Scholar]
  36. Wu, J. Stability of a three-species stochastic delay predator–prey system with Lévy noise. Phys. Stat. Mech. Appl. 2018, 502, 492–505. [Google Scholar] [CrossRef]
  37. Hudson, R.L.; Parthasarathy, K.R. Quantum Ito’s formula and stochastic evolutions. Commun. Math. Phys. 1984, 93, 301–323. [Google Scholar] [CrossRef]
  38. Jiang, X.; Wang, J.; Wang, W.; Zhang, H. A Predictor–Corrector Compact Difference Scheme for a Nonlinear Fractional Differential Equation. Fractal Fract. 2023, 7, 521. [Google Scholar] [CrossRef]
Figure 1. The behavior of the prey and predator subpopulations in the extinction situation.
Figure 1. The behavior of the prey and predator subpopulations in the extinction situation.
Fractalfract 07 00751 g001aFractalfract 07 00751 g001b
Figure 2. The behavior of the prey and predator populations in the extinction–persistence situation.
Figure 2. The behavior of the prey and predator populations in the extinction–persistence situation.
Fractalfract 07 00751 g002aFractalfract 07 00751 g002b
Figure 3. The behavior of the prey and predator populations in the case of persistence.
Figure 3. The behavior of the prey and predator populations in the case of persistence.
Fractalfract 07 00751 g003aFractalfract 07 00751 g003b
Figure 4. Comparison of the behavior of the prey and predator populations the deterministic model and stochastic model in the case of persistence between.
Figure 4. Comparison of the behavior of the prey and predator populations the deterministic model and stochastic model in the case of persistence between.
Fractalfract 07 00751 g004aFractalfract 07 00751 g004b
Figure 5. Comparison of the behavior of the prey and predator populations in the case of persistence with and without Lévy jumps.
Figure 5. Comparison of the behavior of the prey and predator populations in the case of persistence with and without Lévy jumps.
Fractalfract 07 00751 g005aFractalfract 07 00751 g005b
Table 1. The parameters of (1) and their meanings.
Table 1. The parameters of (1) and their meanings.
ParametersDescription
υ 1 The intrinsic growth rate of prey  X 1
υ 2 The intrinsic growth rate for  X 2
ϱ 1 The intra-specific competition for prey  X 1
ϱ 2 The intra-specific competition for prey  X 2
ι 1 The half-saturation due to  X 1
ι 2 The half-saturation due to  X 2
β The impact of the predator interference
δ The death rates of predator
ϱ 3 The rate of intra-species competition for predator
ϑ 1 The average of the transformation of predator from prey  X 1
ϑ 2 The average of the transformation of predator from prey  X 2
dThe natural mortality of the predator y
Table 2. Values of the parameters used in the numerical simulations.
Table 2. Values of the parameters used in the numerical simulations.
ParametersFigure 1Figure 2Figure 3
υ 1 0.7 1.7 2
ϱ 1 233
υ 2 0.65 1.8 2.3
ϱ 2 144
ϑ 1 0.1 0.15 0.2
ϑ 2 0.12 0.17 0.3
ϱ 3 22 1.5
β 0.002 0.0015 0.001
d 0.4 0.35 0.1
σ 1 10 3 10 3 10 2
σ 2 2 × 10 5 2 × 10 4 2 × 10 3
σ 3 2 × 10 5 2 × 10 3 2 × 10 2
Ξ 1 ( u ) 0.006 0.005 0.005
Ξ 2 ( u ) 0.03 0.03 0.03
Ξ 3 ( u ) 0.07 0.06 0.05
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Danane, J.; Yavuz, M.; Yıldız, M. Stochastic Modeling of Three-Species Prey–Predator Model Driven by Lévy Jump with Mixed Holling-II and Beddington–DeAngelis Functional Responses. Fractal Fract. 2023, 7, 751. https://doi.org/10.3390/fractalfract7100751

AMA Style

Danane J, Yavuz M, Yıldız M. Stochastic Modeling of Three-Species Prey–Predator Model Driven by Lévy Jump with Mixed Holling-II and Beddington–DeAngelis Functional Responses. Fractal and Fractional. 2023; 7(10):751. https://doi.org/10.3390/fractalfract7100751

Chicago/Turabian Style

Danane, Jaouad, Mehmet Yavuz, and Mustafa Yıldız. 2023. "Stochastic Modeling of Three-Species Prey–Predator Model Driven by Lévy Jump with Mixed Holling-II and Beddington–DeAngelis Functional Responses" Fractal and Fractional 7, no. 10: 751. https://doi.org/10.3390/fractalfract7100751

APA Style

Danane, J., Yavuz, M., & Yıldız, M. (2023). Stochastic Modeling of Three-Species Prey–Predator Model Driven by Lévy Jump with Mixed Holling-II and Beddington–DeAngelis Functional Responses. Fractal and Fractional, 7(10), 751. https://doi.org/10.3390/fractalfract7100751

Article Metrics

Back to TopTop