Journal of Applied Mathematics
Volume 2019 (2019), Article ID 2601648, 13 pages
https://doi.org/10.1155/2019/2601648

Review Article

Modeling Dynamics of Prey-Predator Fishery Model with Harvesting: A Bioeconomic Model

Charles Raymond, Alfred Hugo, and Monica Kung’aro

Department of Mathematics, University of Dodoma, P.O. Box 338, Dodoma, Tanzania

Correspondence should be addressed to Alfred Hugo; alfredkhugo@gmail.com

Received 7 March 2019; Accepted 15 May 2019; Published 2 June 2019

Academic Editor: Turgut Öziş

Copyright © 2019 Charles Raymond et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

A mathematical model is proposed and analysed to study the dynamics of two-prey one predator system of fishery model with Holling type II function response. The effect of harvesting was incorporated to both populations and thoroughly analysed. We study the ecological dynamics of the Nile perch, cichlid, and tilapia fishes as prey-predator system of lake Victoria fishery in Tanzania. In both cases, by nondimensionalization of the system, the equilibrium points are computed and conditions for local and global stability of the system are obtained. Condition for local stability was obtained by eigenvalue approach and Routh-Hurwitz Criterion. Moreover, the global stability of the coexistence equilibrium point is proved by defining appropriate Lyapunov function. Bioeconomic equilibrium is analysed and numerical simulations are also carried out to verify the analytical results. The numerical results indicate that the three species would coexist if cichlid and tilapia fishes will not be overharvested as these populations contribute to the growth rates of Nile perch population. The fishery control management should be exercised to avoid overharvesting of cichlid and tilapia fishes.

1. Introduction

In today's life, the relationship between predator and prey became an important aspect to discuss in ecology. The prey-predator system has attracted many researchers to study the interaction between the species [1]. Thus, we use mathematical ecology aspect to study the interacting species. The current study considers lake Victoria fishery found in Tanzania as a case study. However we have not gone to the field, but the current trends on the dynamics of species in the lake have been obtained from reading different literature such as NPFMP [2], FAO [3], LVFO [4], Barack [5], and Barilwa [6]. This literature explains the current trend of the lake and the fishery management of the lake Victoria. The lake is comprised of a lot of species such as stocked Nile perch, Lates niloticus, tilapia fish, Oreochromis niloticus, the cyprinid, Rastrineobola argentea, catfishes, insects, cichlids, crocodiles, and many zooplanktons and phytoplanktons NPFMP [2].

The particular study focuses on Nile perch as the predator while cichlid fishes and tilapia fishes are considered as prey populations. All three species are encountering the harvesting aspect. However harvesting without limitations may have detrimental effects on fish population because it decreases the population and sometimes leads a certain species to extinction, Ganguli [7]. The lake Victoria fishery today is either overexploited or in a state of full exploitation because of greater fishing effort and increased competition between fishers, vessels, or nations over the resource. The particular study intends to apply mathematical techniques to ensure the sustainability of the species in lake Victoria without compromising the biological, economic, and social objectives for the benefit of present and future generations. Prey-predator model in fishery was also studied by Kar [8], Chakraborty [9], and Yunfei and Yongzhen [10], while studies by Tapas et al. [11], Ganguli at el [7], Kar [8], Gian [12], Chaudhuri and Kar [13], and Kar [8] analysed the bioeconomic aspect of prey- predator system and observed that increasing harvesting efforts result in population decreases. The particular study intends to analyse the bioeconomic impact for the lake Victoria fishery activities in Tanzania.

2. Materials and Methods

2.1. Model Description, Formulation, and Analysis

It will be assumed that the Nile perch depends completely on cichlid and tilapia fishes as their favorite food (because of easy to capture and their taste) where cichlid and tilapia fishes have unlimited sources of food. The dynamics therefore follow the Holling type II function response. In this case, , , and represent the population of cichlid fish, tilapia fish, and Nile perch, respectively, at any time , where all species involved in modeling are also encountering harvesting aspect. The growth rate of cichlid fishes and tilapia fishes follows the logistic law and the birth rate should always be positive. Terms representing interspecific competition among the prey species are included in the model and the model is then divided into three nonlinear autonomous ordinary differential equations describing how the population densities of the three species would vary with time. The following assumptions are made in order to construct the model:

(i)Cichlid and tilapia fishes have an unlimited food supply in the lake Victoria.(ii)The Nile perch is completely dependent on the cichlid and tilapia fishes as the only favorite food source.(iii)Interspecific competition among tilapia and cichlid fishes is exploitative.(iv)In absence of the predator, prey species grow logistically. That is, the population of the cichlid and tilapia fishes would increase exponentially until it reaches the maximum density of the Lake, which is its environmental carrying capacity .(v)The predation functional response of the Nile perch towards both cichlid fishes and tilapia fishes is assumed to follow Michaelis-Menten kinetics and is modeled using a Holling type II functional form with predation coefficients and and the half saturation constants and .

By considering the underlying assumptions of the incorporated populations, we formulate the system of model equations as:

with initial data values .

All parameters in the model are assumed to be positive and and are per capita intrinsic growth rates of cichlid and tilapia fishes, respectively, while and are environmental carrying capacities of cichlid and tilapia fishes, respectively, and are coefficients for interspecific competition, and are predation coefficients for cichlid fishes and tilapia fishes, respectively, , , and are effort harvesting rates, is natural mortality rate of Nile perch, and and are conversion parameters for cichlid fishes and tilapia fishes by Nile perch, while , , and are catchability coefficients.

For ease of computation, we are rescaling model (1) to reduce the number of parameters as follows: take , then the system of model (1) becomes

with initial values , and

2.2. Equilibrium Points of System (2)

The equilibrium states of the model are obtained by setting and we assume that the predator has positive mortality rate . The following are the possible equilibrium points of the system , , , , , and . Therefore,

(i) The equilibrium point with  From system (2), in the absence of tilapia fish and Nile perch we have and this exists when (ii) The equilibrium point with  From system (2), in the absence of cichlid fishes and Nile perch ( and ) we have and this exists if (iii) The equilibrium point with and  From system (2), in the absence of Nile perch () we have

 This exists if, for and ,  Also exists if, for and , (iv) The equilibrium point with and  From system (2), in the absence of tilapia fishes () we have

 This exists if

  which is possible when and  In terms of original parameter it implies and (v) The equilibrium point with and . From system (2), in the absence of cichlid fishes () we have

 This exists if and this is possible when and In terms of original parameter it means and . The condition implies that and (vi) Coexistence equilibrium point

Following the procedure by Dubey [14], the endemic equilibrium point is obtained as follows:

From (7) we have

From (8) we have

From (8) and (9) we have

From (10) and (11) we have

and from (10) and (12) we get

Equations (13) and (14) are two functions of and . To prove the existence of , the conditions under which and meet in the interior of the positive plane at the point are found. Now the values of , and can be obtained from (7), then from (9) we observe that, as , tends to . is the value of at which the function would cut the axis in the plane. So is given by

We notice that is the same as of . From (8), as , y tends to given by

where

and are the points at which the functions and would cut the y-axis in the plane, respectively. Also from (13), where

and

We note that if and and this requires and . Similarly from (14) where

and

We also note that if and and this requires and . Since for , we have and for , we have , then and will meet if . We therefore state the existence of the positive equilibrium point in the following theorem.

Theorem 1. The positive equilibrium point will exist if the following conditions are satisfied:

where and are as defined in (15) and (16), respectively.

In terms of original parameter, implies that ; i.e., the growth rate of cichlid fishes must be greater than the harvesting effort imparted. Condition gives . That is, the harvesting rate of Nile perch should be greater than the harvesting rate of tilapia fishes. For economic purpose the condition needs to be satisfied which implies that the rate at which the Nile perch is harvested must be greater than its death rate.

2.3. Local Stability of Equilibrium Points

To analyse the local stability of the equilibrium point we consider the Jacobian matrix;

where

(i), the Jacobian matrix evaluated at gives the eigenvalues, , , and . We see that and are always positive and so is unstable.(ii). The Jacobian matrix (23) is evaluated at with the following eigenvalues:

 The eigenvalues are negative if and . Hence, the equilibrium point is locally asymptotically stable if the following conditions hold: , , and . The condition implies . For the local stability of , , the harvesting rate of cichlid fishes must be less than their intrinsic growth rate. Other inequalities show parameters that are vital for the local stability of .(iii). The Jacobian matrix (23) is evaluated at and the following eigenvalues obtained:

 The eigenvalues above are negative if . Hence, the equilibrium point is locally asymptotically stable if conditions hold. The condition implies . For the local stability of , , the harvesting rate of tilapia fishes must be less than their intrinsic growth rate.(iv) = . The eigenvalues of are obtained by solving the characteristic equation;

 where , , , , , , and . This can be expressed in form of

 By Routh-Hurwitz criteria (Murray, 1989), the are negative if .(v) = . The eigenvalues are obtained by solving the characteristics equation;

 where

 The characteristic equation is in the form

 where , , and . By Rouths stability criterion, the equilibrium point is stable if (i) , , and ; (ii) . Otherwise it is unstable.(vi). The eigenvalues evaluated at are obtained by solving the characteristic equation

  where , , and . The characteristic equation is in the form ;  where , , and . By Rouths stability criterion, the equilibrium point is stable if (i) , , and and (ii) . Otherwise it is unstable.

2.4. Global Stability of the Coexistence Equilibrium Point

To analyse global stability of coexistence equilibrium point, a suitable Lyapunov function is chosen, from which conditions for the global asymptotic stability of the coexistence point are derived. The approach is based on work by Chaudhuri [13] and Dubey [14].

Theorem 2. The coexistence equilibrium point is globally asymptotically stable if
(i) , (ii) , and (ii) .

Proof. Consider the following Lyapunov function:

Simplifying the above Lyapunov function,

where , ,and .
Differentiating with respect to time, i.e.,

which simplifies to

Thus, is a quadratic form which can be expressed as , where and A is symmetric matrix given by

with , , , , , and . We note that the point is globally asymptotically stable if ; that is, the matrix A is positive definite Chaudhuri [13]. Now the matrix A is positive if , , , , , and . gives , gives , gives , gives , and gives . This completes the proof.

2.5. Bioeconomic Equilibrium

The term bionomic equilibrium is an amalgamation of the concepts of biological equilibrium as well as economic equilibrium Kar [8]. From system (2), a biological equilibrium is given by . The economic equilibrium is said to be achieved when the total revenue obtained by selling the harvested biomass (TR) equals the total cost for the effort devoted to harvesting (TC).

Let be the fishing cost per unit effort for cichlid fishes, the fishing cost per unit effort for tilapia fishes, the fishing cost per unit effort for Nile perch, the price per unit biomass of cichlid fishes, the price per unit biomass of tilapia fishes, and the price per unit biomass of Nile perch; , , and are catchability coefficients of cichlid fishes, tilapia fishes, and Nile perch, respectively.

Then we have , , and .

, , and are the economic rent (net revenue) of cichlid fishes, tilapia fishes, and Nile perch, respectively. Therefore, the economic rent (net revenue) at any time is given by , which is

Then is the bioeconomic equilibrium where , , , , , and are the bioeconomic values of cichlid fishes, tilapia fishes, Nile perch, harvesting effort of cichlid fishes, harvesting effort of tilapia fishes, and harvesting effort of Nile perch, respectively, and it is given by the simultaneous equation:

In order to determine the bioeconomic equilibrium, we now consider the following cases.

Case I. If , that is, the fishing cost per unit effort for cichlid fishes is greater than the revenue in the cichlid fish fishery, then fishermen will be in loss and naturally they would withdraw their participation from cichlid fish and the fishery will be closed . Only tilapia fish and Nile perch fishery remain operational (i.e., and ). Now we have and . From (39) when , we have

Since and , then (44) can be written in quadratic form as

Therefore

where , , .
From (46),

Therefore, if

Also,

Thus, if

Therefore, the bioeconomic equilibrium exists if conditions (48) and (50) hold.

Case II. If , that is, the fishing cost per unit effort for tilapia fishes is greater than the revenue in the tilapia fish fishery, then fishermen will be in loss and naturally, they would withdraw their participation from tilapia fish fishery and the fishery will be closed . Only cichlid fishes and Nile perch fishery remain operational (i.e., and ). Now we have and . From (46) when , we have

Since and , then (52) can be written in quadratic form as

Therefore

where

From (39),

if

Also,

Thus, if

Therefore, the bioeconomic equilibrium exists if conditions (57) and (59) hold.

Case III. If , that is, the fishing cost per unit effort for Nile perch is greater than the revenue in the Nile perch fishery, then fishermen will be in loss and naturally, they would withdraw their participation from Nile perch fishery and the fishery will be closed. Only cichlid fish and tilapia fish fishery remain operational (i.e., and ). Now we have and . Substituting and in (44), we get

Thus, if

Also,

Thus, if

Therefore, the bioeconomic equilibrium exists if conditions (61) and (63) hold.

Case IV. If and , then the fishing cost is greater than revenues for all three species and the whole fishery will be closed.

Case V. If and , then the fishing cost is less than revenues for all three species; that is, the fishery is more profitable and hence it would attract more fishermen and the whole fishery will be in operation. Now we have , and , then

Thus, , if

And,

Thus, ,

Also,

Thus, if

Thus, the nontrivial bioeconomic equilibrium point exists if and only if conditions (65), (67), and (69) hold together.

3. Results and Discussion

3.1. Numerical Results

The model system is simulated using the inbuilt ODE solvers coded in Matlab programming language and figures are plotted using parameter values presented in Table 1.

Table 1: Parameter Values of the model.

Parameter
Value
source

2.07
[13]

2.09
[13]

200
[13]

100
[13]

0.001
[15]

0.02
[15]

0.002
[16]

0.03
[16]

1.5
[17]

0.1
[17]

0.2
[17]

0.14
[17]

0.13
[18]

0.125
[18]

1
[17]

Figure 1 presents the dynamics of the population interaction with respect to time. The figure indicates that cichlid fishes population grows faster compared to other species. The dynamics was also visualized in as presented in Figure 2. Figure 3 shows the impact in harvesting cichlid fish population; it is observed that increasing harvesting effort leads to decrease Nile perch population as well. In Figure 4 the same scenario happens when tilapia fish population harvested tends also to decrease Nile perch population while increasing harvesting effort in Nile perch population tends to increase cichlid and tilapia fishes as presented in Figure 5. Figure 6 presents the exponential growth rate of cichlid and tilapia fishes in absence of Nile perch population.

Figure 1: Graph of prey and and predator against time (in years). It shows the variation of population density with time.

Figure 2: Graph of prey and and predator against time (in years). The limit cycle shows that population density will slightly change but not oscillate the boundary.

Figure 3: Variation of the fishes population with different harvesting efforts in cichlid with fixed values of and .

Figure 4: The effect of Tilapia harvest with different values of with fixed values of and .

Figure 5: The effect of Nile perch harvest with different values of with fixed values of and .

Figure 6: The interaction between cichlid and tilapia fishes in absence of Nile perch.

3.2. Discussion

A mathematical model was proposed and analysed to study the dynamics of a two-prey-one predator system with harvesting aspects. The model was used to study the ecological dynamics of the Nile perch-cichlid-tilapia fishes prey-predator system of the lake Victoria fishery. The harvesting rate was found to play a crucial role in stabilizing the system. Figures 3 and 4 show that cichlid and tilapia fishes tend to extinction when the harvesting rates and exceed their intrinsic growth rates and . For sustainability harvesting of cichlid and tilapia fishes their intrinsic growth rates should be kept smaller. However, harvesting of the three species at a rate much lower than their intrinsic growth rate would not lead to collapse of the system. This result is similar to one obtained by Chaudhuri and Kar [13] on the existence conditions for the system they studied. In absence of the predator, the two prey species can coexist and are stable if the interspecific competition among them is maintained at minimum level or negligible also if both prey species are not harvested beyond its intrinsic growth rate. Numerical analysis results indeed confirmed this as shown in Figure 6.

Theorem 1 showed that the three species would coexist if cichlid and tilapia fishes were not harvested beyond their intrinsic growth rates; the Nile perch converted the biomass of cichlid and tilapia fishes into fertility at a rate greater than the Nile perch’s natural mortality rate and the time it took to handle the cichlid or tilapia fish. The findings of Dubey and Upadhyay [14] indicated that the predator’s mortality rate and food conversion coefficients played a crucial role in determining the stability behaviour of the equilibrium points.

In order for the bioeconomic equilibrium to exist, the fishing cost per unity effort for all species, price per unity biomass, catchability coefficient, harvesting effort, intrinsic growth rate of prey species, and the mortality rate of predator play a vital role. Hence for the bioeconomic equilibrium to exist, the fishing cost per unity effort for all three species should be less than the revenue in their fishery. Similar results were also obtained by Ganguli and Kar [7] when they studied the optimal harvesting of a prey-predator model with variable carrying capacity.

One of the major observations from numerical simulation results is that the predator population density increased significantly when the harvesting rate of both prey species decreased. This implies that a gradual increase in the number of cichlid and tilapia fishes would result in significant increase in the number of Nile perch and vice versa, which would in the long term lead the population density of the cichlid and tilapia fishes to fall to a level lower than the original one. Kar and Chaudhuri [16] also discovered that predator population density can increase significantly when the harvesting rate of prey species is decreased.

4. Conclusion

We deduce from both analytical and numerical results that if the harvesting rate of the cichlid and tilapia fishes exceeds their intrinsic growth rate, the population of the cichlid, tilapia fishes, and Nile perch would become extinct with time. However, analytical and numerical results also show that harvesting of the three species at a rate much lower than their intrinsic growth rate would not lead to collapse of the system. Thus, in order to use fish as a resource and produce maximum economic benefit while maintaining sustainable fishery species, the harvesting rate of species should never be allowed to exceed their growth rate.

Conflicts of Interest

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

References

M. Sambath and K. Balachandran, “Spatiotemporal dynamics of a predator-prey model incorporating a prey refuge,” Journal of Applied Analysis and Computation, vol. 3, no. 1, pp. 71–80, 2013.

NPFMP, “Nile Perch Fishery Management Plan for Lake Victoria 2015 - 2019 Ebene, Mauritius,” 2015.

FAO, “Poverty in coastal fishing communities, Advisory commettee on fishery research,” Rome, Italy, 2015.

LVFO, Regional Status Report on The Lake Victoria Biennial Frame Surveys between 2000 And 2012 Lake Victoria Fisheries Organization, 2013.

A. Barack, C. Emily, and F. Joanne, A Mathematical Ecohistory of Lake Victoria, Dartmouth College, Dartmouth, 1998.

J. S. Barilwa, “The Present status of the fishery for tilapia in Lake Victoria(Uganda), Mwanza,” in Proceeding of the workshop on the Prevailing Activities on the Lake Victoria Basin with Particular Reference to the fisheries of Lake, 1990.

C. Ganguli and T. K. Kar, optimal harvesting of prey predator model with variable carrying capacity, Indian Institute of Engineering Science and Technology, India, 2017.

T. K. Kar, Global Dynamics and Controllability of a Harvested Prey-Predator System with a Holling Type III Functional Responce, Yokohama National university, Japan, 2011.

K. Chakraborty, K. Das, and T. K. Kar, “An ecological perspective on marine reserves in prey-predator dynamics,” Journal of Biological Physics, vol. 39, no. 4, pp. 749–776, 2013.

L. Yunfei, Y. Rong, and P. Yongzhen, “A prey-predator model with harvesting for fishery resources with reserve area,” Journal of Applied Mathematical Modelling, vol. 37, pp. 3048–3062, 2012.

D. Tapas, R. Mukherjee, and K. Chaudhuri, “Harvesting a prey- predator fishery in the presence of toxicity,” Journal of Biological Dynamics, vol. 3, no. 5, pp. 447–462, 2008.

G.-I. Bischi, F. Lamantia, and D. Radi, “A prey-predator fishery model with endogenous switching of harvesting strategy,” Applied Mathematics and Computation, vol. 219, no. 20, pp. 10123–10142, 2013.

T. K. Kar and K. S. Chaudhuri, “Harvesting in a two-prey one-predator fishery: a bioeconomic model,” The ANZIAM Journal, vol. 45, no. 3, pp. 443–456, 2004.

B. Dubey and R. K. Upadhyay, “Persistence and extinction of one prey and two Predator system,” Non Linear Modeling and Control Journal, vol. 9, no. 4, pp. 307–329, 2004.

E. Akukizibwe, Modeling and Analysis of a Two Prey-One Predator System with Harvesting, Holling Type II and Ratio-Dependent Responses, Makerere University, Uganda, 2010.

T. K. Kar and K. S. Chaudhuri, “On nonselective harvesting of a multi-species fishery,” International Journal of Mathematical Education in Science and Technology, vol. 33, no. 4, pp. 543–556, 2002.

T. K. Kar, S. K. Chattopadhyay, and K. R. P. Chandan, “A Bio economic model of two prey one predator system,” International Journal of Applied Mathematics and Informatics, vol. 27, pp. 1411–1427, 2009.

T. K. Kar, K. Chakraborty, and U. K. Pahari, “A prey-predator model with alternative prey: mathematical model and analysis,” Canadian Applied Mathematics Quarterly, vol. 18, no. 2, pp. 137–168, 2010.