Dynamical analysis and optimal harvesting of conformable fractional prey-predator system with predator immigration.

Aim of this work is to study the four species various fractional-order prey-predator or Lotka-Volterra (LV) system with both immigration and harvesting effects. The existence and uniqueness, uniform boundedness, persistence, permanence, and extinction of this system solution are analyzed. The stability behavior of the system is obtained with the help of the Routh-Hurwitz (RH) stability criterion. The small changes in fractional-order values can produce a significant impact on the stability of the system is confirmed. This work verifies that the small amount of immigration effect can change the dynamic nature of the LV system. Numerical results are given to illustrate the obtained theoretical results of the stability analysis. The bionomic equilibrium points of the system are attained with their feasibility conditions. To get the optimal amount of harvesting effect with the Pontryagin's maximum principle, the harvesting parameter is considered as the control parameter.

Introduction

The ecosystem has a wide diversity in the population-community. Such a population-community in an ecosystem is unlikely to survive without dependence on other communities. In that case, there must be a significant relationship between different population-communities. Interactions between different populations are generally characterized by mutualism, competition, parasitism, and predation. In the above types of relationships, predatory activity is widespread in real-time. The predatory relationship between predator and prey is widely seen in real-life such as Lotka and Volterra (LV) system (see [1]). This LV system is employed to characterize the dynamical behavior of the species by their predatory relationship (see [2, 3]). Recently, mathematicians and ecologists have studied a wide variety of prey–predator systems (see [4-8]).

Many researchers studied a large number of fractional-order aspects of the biological process (see [9-13]). Those aspects should not be shown and explained by the integer-order calculus, but the fractional calculus can do it. The fractional-order derivative is related to the whole period of the biological process when compared with the integer-order derivative. Fractional-order population dynamics have advantages over the classical integer-order population dynamics. Fractional derivative is the generalization of the integer-order derivative. By the interests of fractional calculus have brought a comprehensive study from researchers. There are a lot of definitions of fractional derivatives with different properties. Riemann–Liouville’s definition and Caputo’s definition are the most popular ones. But, both of them have incredible difficulties and do not satisfy some significant properties that have met the integer-order calculus [14]. For example, both the Riemann–Liouville definition and Caputo definition derivatives do not fulfill the product rule, quotient rule, chain rule and have some challenges in the calculation. Many fractional-order physical/biological dynamics satisfy these properties. In 2014, Khalil et al. introduced a new fractional calculus called conformable fractional calculus in [15]. This conformable fractional calculus satisfies the properties of the integer-order calculus, such as derivative of the product and quotient of two functions and chain rule (see [16] and [17]). By these advantages, we are motivated to analyze the conformable fractional-order LV system.

It is known that almost all population dynamics have a long-term memory. It means that the dynamical behavior of the LV system also depends upon the fractional orders. In this study, we have proposed the new problem of the four species conformable-fractional-order LV system with the various fractional orders, In this work, we have provided evidence for the small changes of fractional orders that can provide a significant impact on the dynamic behavior of the LV system. In the literature, several valuable results such as stability analysis, persistence, permanence, and extinction are studied about the LV system (see [18-29]).

In recent days, we have to see the two unfortunate pandemic situations are made by the impact of immigration. The first one is, the migration of that asymptomatic infected individuals, as well as population movements, can play a crucial role in spreading the COVID-19 virus. And the next one is a large number of locusts migrate into the region and cause damages to the crops together with the local locusts. From these two cases, we can observe that immigration is playing an important role in population dynamics. In an ecosystem, all the predator individuals are depending on their prey individuals for their food. In this situation, they need to move or migrate from one region to another for their food where the highest density level of the prey species is surviving. According to the migration on predator, the stability of the ecosystem will be changed. So that the immigration on predator has a vital role in the stability of prey–predator interactions (see [18-24]). By this importance of predator’s immigration, we developed the fractional-order LV system with the predator immigration effect. Also, the impact of immigration on the stability behavior of the LV system is analyzed. Our obtained results generalize the existing integer-order LV system with the immigration factor (see [18-24]).

If the dynamical result of the LV system has stable with the permanence of all species, the system is regarded as a healthy ecosystem. Unfortunately, we can see that the existence of stability in the environment has contradictory issues. Appropriate action by the majority level of individuals of an area is the cause of the unstable LV system. If that action of particular individuals cannot be controllable, then we should prepare a supplement strategy that would prevent the destruction of the ecosystem. One of the best strategies is to apply the harvesting force as a control parameter on the particular individuals that are contributing to the damage (see [25]). In the literature (see [26-29]), many authors studied the LV system with the harvesting effect. From this analysis, we can observe that the immigration and harvesting effects having much influence on the stability of the LV system. To the best of the authors’ knowledge, the four species LV system on the fractional order has not been investigated with both immigration and harvesting effects. With this motive, this article investigates the new problem of a conformable fractional-order LV system with the immigration and harvesting effects. The main contributions of this work are given as follows:This paper is organized as follows: Sect. 2 gives the preliminary results based on the conformable fractional-order derivative and constructs the four species conformable fractional-order LV system with immigration and harvesting effects. Section 3 is dedicated to analyzing the existence and uniqueness of the solution of the system and described the solution behavior, uniform boundedness, persistence, permanence, and extinction results. Section 4 verifying the impact of fractional orders and the immigration on stability behavior of the system with the numerical example. In Sect. 5, the optimal harvesting policy of the LV system is obtained by the Pontryagin’s maximum principle at the non-trivial bionomic equilibrium points.

The proposed conformable fractional-order four species LV system includes prey, secondary-level predator, third-level predator, and top-level predator with the various fractional-orders

In this model, each species having intercommunication with another three species. Moreover, the immigration effect is applied to the predator species and harvesting effects are implemented in all species.

The impact of the fractional orders and immigration effect on the dynamical behavior of the LV system is investigated. From this investigation, we can see how the small changes in fractional orders and immigration effect can change the dynamical behavior of the system.

In this study, we have newly achieved the optimal harvesting of the fractional-order LV system by the Pontryagin’s maximum principle at the bionomic equilibrium points.

Preliminaries and system description

In this section, the required basic definition and theorem are defined and a description of the proposed conformable fractional-order LV system is presented.

Conformable fractional calculus

Definition 1

(See [15]) Given a differentiable function Then, the conformable fractional derivative of f of order is defined by

Theorem 1

(See [15]) Let and f be differentiable at a point then

Conformable fractional-order LV system

Let , and denote the densities of prey, secondary-level predator, third-level predator, and top-level predator, respectively, at time For the notation purpose, is simply notated as Assume that is a prey of its predators and ; is a prey of its predators and ; and is a prey of its predator . Species grow logistically with the carrying capacity respectively. The immigration effect is applied to the predator species respectively. The harvesting effect is applied to the species respectively. With the above assumptions, the following four-dimensional conformable fractional-order LV system is configured,with the initial condition at for all The non-negative parameters and are intrinsic growth rate and carrying capacity of each , respectively, for simplification is denoted as and denotes the interaction coefficient between the each species , , The parameter denotes the immigration rate of the each predator species Also, the parameters and denote the death rate, catch-ability coefficient, and total effect applied for harvesting of the each species , respectively.

For all the LV system (1) of equations can be defined in the following general form:Here, defined by and

According to Theorem 1 in [15], for the fractional-order LV system (2) can be written as,Here, The integer-order LV system (3) (corresponding to the fractional-order LV system (1)) can be written in the general form with the initial condition Here,

Remark 1

Comparing the LV systems (2) and (4), we can understand that the integer-order LV system (4) is the equivalent form of the fractional-order LV system (2).

Solution behavior of the system

This section verifying the existence and uniqueness of the system solution, uniform boundedness, persistence, permanence, and extinction of solution of the LV system (2).

Existence and uniqueness of the system solution

Theorem 2

Assume that the function defined in the LV system (4) is Lipschitz continuous and there exists a constant sufficient for the existence and uniqueness of solution X(t) of the LV system (2) with the initial condition in the region Here,

Proof

Let X(t) and . The solution X(t) of the LV system (2) with the initial condition and can be written as,In Eq. (5), put then It can be easily obtain that,Here, and T denote the transpose of the matrix. For all Choosing supremum norm, for continuous function , , and for the matrix Then, we can easily obtained that,Here, and From the inequality (6), the mapping is a contraction if By contraction principle, has a unique fixed point X(t) which satisfies (5). Hence, X(t) is unique solution of the LV system (2).

Uniformly bounded, persistent, permanent and extinction

Theorem 3

(Uniformly bounded) If the conditions and for hold, then the LV system (2) is uniformly bounded.

For , define the function, and The derivative If we choose and for then we obtain,Here, By applying the Gronwall’s inequality, then . As we have Hence, the solution X(t) of the LV system (2) in are restricted in the region . Hence, all the species are uniformly bounded in for any initial .

Remark 2

Biologically, it is confirmed that the LV system (2) is uniformly bounded when the following two conditions are satisfied.

Sum of the growth rate and immigration rate is equal to the sum of death rate and harvesting effect.

The prey consumption rates are equal to the corresponding predation rates.

Definition 2

The LV system (2) is said to be weakly persistent if every solution X(t) satisfies, the following two conditions are satisfied. The LV system (2) is said to be strongly persistent if every solution X(t) satisfies the following two conditions.

for all

for all

Definition 3

The LV system (2) is said to be permanent for the solution X(t) with if there exist constants such that, and .

Theorem 4

(Persistent and permanent) If and hold, then the LV system (2) is persistent and permanent, i.e., any positive solution X(t) of the LV system (2) satisfies

As the solution X(t) of the LV system (2) is positive, it follows that,By the Lemma 3 in [27], and for arbitrary there exist positive constants According to the assumptions, and , the LV system (2) is equivalently written as,Here, for all Then, By the Lemma 3 in [27], and for arbitrary , there exist positive constants From (7) and (8), it is clear that the LV system (2) is persistent and permanent.

Remark 3

From this analysis, we have observed that, if the predator species’ level should not exceed 1, then the LV system (2) is both persistent and permanent.

Definition 4

The LV system (2) is said to be extinction if there is a positive solution X(t) that satisfies

Theorem 5

(Extinction) If the following two conditions,

hold for all then the species will be extinct as

As the solution X(t) of the LV system (2) is positive, it follows that:Since, , we get It means that the species is extinct as . Using the assumption, , then the LV system (2) can be written as,Using the assumption, , for then Hence, the species are extinct as .

Remark 4

As the prey species will be extinct if sum of the death rate and harvesting rate of is greater the growth rate of . As , the other three predator species are extinct if their corresponding prey levels should not exceed 1 and the sum of the death rate and harvesting rate of is greater than the sum of the growth, immigration and consumption rates of

Stability analysis

This section analyzes the stability properties of the LV system (2) at the positive interior equilibrium point . Here, the stability analysis of the LV system (2) is obtained via the Routh–Hurwitz stability test.

Routh–Hurwitz stability test

First to construct the Jacobian matrix for all and of the LV system (2) at the positive interior equilibrium point for analyzing the stability.For, The characteristic equation of the Jacobian matrix at the positive interior equilibrium point ,Here,Thus, the above discussion gives the following stability theorem.

Theorem 6

[30]

The positive interior equilibrium point is stable for the LV system (2) if it satisfies (Routh–Hurwitz (RH) stability test) and with values are defined in Eq. (10).

Remark 5

From Theorem 6, it can be observed that the fractional-orders (here, viewed in terms of ) also perform an essential role in stability analysis. It is proved that the stability of the LV system (2) depends on fractional-orders of the LV system (2). In Sect. 4.2, stability behavior of the LV system (2) with the various fractional orders is graphically illustrated.

Remark 6

If the condition hold for the LV system (2), then the positive equilibrium point of the LV system (2) exists.

The impact of fractional-orders in stability analysis

To demonstrate the impact of fractional-orders in stability of the LV system (2) with the parameter set 1 and a different set of fractional orders.

Parameter set-1

Consider the following non-negative parameters with the initial conditions

Remark 7

The parameter set-1 with the fractional orders does not satisfy both the RH stability test in Theorem 6 (see Table 1) and the condition of the positive interior equilibrium point . So that the LV system (2) is unstable (see Fig. 1). The parameter set-1 with the fractional-orders satisfies both the RH stability test in Theorem 6 (see Table 1) and the condition of the positive interior equilibrium point . So that the LV system (2) is stable (see Fig. 2). The parameter set-1 with the integer-orders satisfies both the RH stability test in Theorem 6 (see Table 1) and the condition of the positive interior equilibrium point . So that the LV system (2) is stable (see Fig. 3). From this analysis, the small changes in the fractional orders of the system can stabilize and destabilize the LV system. Hence, the fractional orders having much influence on the dynamical behavior of the LV system (2) (see Table 1). The dynamical and phase-space representations of the unstable and stable behaviors of the LV system (2) are graphically represented in Figs. 1, 2, 3.

Table 1

Impact of fractional orders in stability analysis—RH stability test

α_1,α_2,α_3,α_4	E(x1∗,x2∗,x3∗,x4∗)	G1′	G4′	u_1	u_2	Stability
0.8, 0.95, 0.85, 0.5	x1∗=0.0823, x2∗=0.0002,	−74.4946	2.9102×10^5	-8.0896×10^4	-2.7694×10^8	Unstable
	x3∗=0.0000, x4∗=0.1551
0.95, 0.79, 0.9, 0.8	x1∗=0.0989, x2∗=0.2273,	705.8360	3.4405×10^7	1.6386×10^14	1.8348×10^10	Stable
	x3∗=0.0680, x4∗=0.1598
1, 1, 1, 1	x1∗=0.2690, x2∗=0.5188,	2.3781×10^3	2.9999×10^11	6.3612×10^9	9.0996×10^18	Stable
	x3∗=0.1527, x4∗=0.4366

Fig. 1

Unstable system (2) with

Fig. 2

Stable system (2) with

Fig. 3

Stable system (2) with

The impact of immigration on a stability analysis

To analyze the impact of immigration on stability of the LV system (2) with the parameter set-2.

Parameter set-2

Consider the following non-negative parameters with the initial conditions

Remark 8

From this case study of the immigration effect on each species of the LV system (2), in the following cases 2, 3, 5, 6, and 7, the LV system (2) is unstable. Since the parameter set-2 with the respective various set of immigration parameters does not satisfy conditions of both the RH stability test in Theorem 6 (see Table 2) and positive interior equilibrium point. In cases 1 and 4, the LV system (2) is stable (see Table 3). Since the parameter set-2 with the respective set of immigration parameters satisfying both the RH stability test in Theorem 6 (see Table 2) and the condition of the positive interior equilibrium point. From this case study on the immigration effect, we can conclude that the small changes in the immigration effect have a significant impact on the dynamic behavior of the LV system (2). The dynamical and phase-space representations of the stability (unstable or stable) behaviors of the LV system (2) are graphically represented in figures from Figs. 4, 5, 6, 7, 8, 9, 10.

Table 2

Impact of immigration in stability analysis—RH stability test

Case	Immigration on xj′s	E(x1∗,x2∗,x3∗,x4∗)	G1′	G4′	u_1	u_2	Stable/
	(I_2,I_3,I_4)						Unstable
1	x_2	x1∗=0.0000, x2∗=2.3934,	5.3492×10^3	8.0895×10^10	1.2191×10^21	9.5081×10^11	Stable
	(I_2,0,0)	x3∗=1.4523, x4∗=0.5047
2	x_3	x1∗=0.0000, x2∗=0.0403,	0.8180×10^4	-1.7612×10^8	-2.7317×10^16	-1.6076×10^10	Unstable
	(0,I_3,0)	x3∗=0.0006, x4∗=0.4325
3	x_4	x1∗=0.7005, x2∗=0.0195,	-2.0632×10^3	5.6309×10^10	9.0737×10^8	-1.2600×10^18	Unstable
	(0,0,I_4)	x3∗=0.0000, x4∗=1.0706
4	x_2,x_3	x1∗=0.0000, x2∗=2.3934,	4.5611×10^3	2.3042×10^11	6.1983×10^10	2.8316×10^20	Stable
	(I_2, I_3,0)	x3∗=1.4522, x4∗=0.5047
5	x_2,x_4	x1∗=0.1512, x2∗=0.6920,	-0.1876×10^4	2.6597×10^10	3.1555×10^8	-1.4773×10^17	Unstable
	(I_2,0,I_4)	x3∗=0.0000, x4∗=1.1076
6	x_3,x_4	x1∗=0.2628, x2∗=0.1075,	0.8204×10^4	-6.6101×10^10	-4.3843×10^8	-4.5115×10^17	Unstable
	(0,I_3,I_4)	x3∗=0.0008, x4∗=1.0778
7	x_2,x_3,x_4	x1∗=0.0113, x2∗=0.6476,	0.6829×10^4	-7.5470×10^10	-4.6353×10^8	-5.2534×10^17	Unstable
	(I_2,I_3,I_4)	x3∗=0.0000, x4∗=1.1272

Table 3

Bionomic equilibrium points

Case	Feasibility conditions	Bioinomic equilibrium points
1	(1) 3.0100 > 0.0662	x_1∞=0.1110, e_1∞=3.2709
	(2) 0.1882 > 0.0088	x_2∞=0.1110, e_2∞=0.1994
	(3) 0.2152 > 0.0457	x_3∞=0.1250, e_3∞=0.2119
	(4) 0.0669 > 0.0163	x_4∞=0.1429, e_4∞=0.0723

Fig. 4

Case-1 Stable LV system (2) with the immigration

Fig. 5

Case-2 Unstable LV system (2) with the immigration

Fig. 6

Case-3 Unstable LV system (2) with the immigration

Fig. 7

Case-4 Stable LV system (2) with the immigration

Fig. 8

Case-5 Unstable LV system (2) with the immigration

Fig. 9

Case-6 Unstable LV system (2) with the immigration

Fig. 10

Case-10 Unstable LV system (2) with the immigration

Bionomic equilibrium and optimal harvesting

In this section, the trivial and non-trivial bionomic equilibrium points of the LV system (2) are achieved. The optimal harvesting policy of the LV system (2) is gained by the Pontryagin’s maximum principle.

Bionomic equilibrium

In this subsection, we analyze the bionomic equilibrium points of the LV system (2) with the harvesting effect applying to all the species. Let and denote the harvesting cost per unit effect and price per unit of , respectively. Therefore, is the net revenue function defined by Here, represent the net revenues for respective The bionomical equilibrium point is obtained from the following simultaneous equations:Here, we discuss the existence of the bionomic equilibrium, which can be done in the following two cases (trivial and non-trivial). If , then the harvesting cost is greater than the revenue of each This means harvesting effect is not applied on the species , . If , then the harvesting cost is less than the revenue of each This means harvesting is applied on the species .

Case 1: .

If the harvesting cost of the all species does not exceed the revenue cost, then the harvesting effect is applied to all species , Hence, From (11), we have If the following four feasibility conditions are hold, then Case 2:

If the harvesting cost of the all species exceeds the revenue cost, then the LV system (2) will be closed.

Thus, the above discussion gives the following result.

Theorem 7

The LV system (1) has a positive non-trivial bionomic equilibrium point if the conditions 12 to 15 are satisfied.

Example 5.1 Here, we analyzed the existence of the positive non-trivial bionomic equilibrium points of the LV system (2) with the parameter set-3.

Parameter set-3

Consider the following non-negative parameters , , , , , , , , , , and with the initial conditions

Remark 9

Using the Theorem 7, the existence of the non-trivial interior equilibrium point and its feasibility conditions (12) to (15) of the LV system (2) is presented in Table 3.

Optimal harvesting policy

Here, we discuss the optimal harvesting policy by using the Pontryagin’s maximum principle. Define the cost function Here, represents the instantaneous annual discount rate. To maximize the cost function subject to the LV system (2) and the control variables satisfying Here, is the feasible upper limit of harvesting effect Now, construct the Hamiltonian function H, Here, are the adjoint variables. From (16), for we get,Here, is called the switching function. The optimal control should be a combination of the extreme control and the singular control as Hamiltonian H is a linear in the control variables Clearly, the optimal control that maximizes H satisfies,Consequently, the optimal harvesting policy satisfiesHere, is called singular control of the optimal control , which satisfying According to (17), we obtain the following equations, for To determine the optimal harvesting effect, introduce the corresponding adjoint equations, From Eq. (16), we get, for With the help of Eqs. (11), and substituting Eqs. (18) in (19), and integrating, for Here, we omit the integration constants to ensure that are bounded for as Substituting (18) in (20), for we get the following control parameter ,Here, The optimal harvesting effect can be determined by Eqs. (21).

Remark 10

We have analyzed the optimal harvesting policy by applying Pontryagin’s maximum principle. From the economic and biological aspects of renewable resource management, the harvesting effect is used for all four species For environment management, planning of harvest strategies, and keeping sustainable growth are needed. In the ecosystem, we have been concerned with the harvesting effect as a control parameter and its optimal point is . If the harvesting effect is applied less than the optimal point , then all species will coexist. It can control the optimal level and environmental balance. If the harvesting effect is used greater than the optimal point , then the ecosystem will go to dangerous destruction.

Conclusion

In this work, the new problem of the conformable fractional-order four species LV system is built with both immigration and harvesting effects. The existence and uniqueness of the system solution, uniform boundedness, persistence, permanence, and extinction of the LV system (2) are obtained. The stability of the system is obtained with the RH stability criterion. This work verified the small changes in fractional-orders, and the immigration effect gives a significant impact on the dynamic results of the LV system (2). The numerical results are gained for the obtained theoretical results. The optimal harvesting policy of the LV system (2) is obtained at the bionomic equilibrium point by Pontryagin’s maximum principle. Nowadays, researchers are concentrating their study on the mathematical modeling of infectious diseases known as the SIR (Susceptible–Infected–Recovered) epidemic model. In the future, it is motivated that the fractional-order SIR infection model as an essential research. Also, still it is open for stabilization of conformable fractional-order time-delay population system in this direction.