Abstract

The present paper deals with a dynamic reaction model of a fishery. The dynamics of a fishery resource system in an aquatic environment consists of two zones: a free fishing zone and a reserve zone. To protect fish population from over exploitation, a control instrument tax is imposed. The existence of its steady states and their stability are studied. The optimal harvest policy is discussed next with the help of Pontryagin's maximum principle. Our theoretical results are confirmed by numerical simulation.

1. Introduction

With the growing need of human for more food and energy, several resources have been increasingly exploited. It has caused wide public concern to protect the ecosystem. A scientific management of commercial exploitation of the biological resource like fisheries and forestry is necessitated. The first attempt of mathematical modelling of resource management problems was made in the article by Hotelling [1], who dealt with the economies of exhaustible resources using the calculus of variations, which was not a very familiar topic for the researchers as Hotelling’s article was largely ignored until the 1973 energy crisis which led to a surge of interest in the field [2]. The techniques and issues associated with bioeconomic exploitation of these resources have been discussed in detail by Clark [3, 4]. During the last few decades, several investigations regarding fishing resource have been conducted [516]. Dubey et al. [17] proposed and analyzed a mathematical model as follows: ????(??)=????(??)1-??(??)???-??1??(??)+??2????(??)-??????(??),??(??)=????(??)1-??(??)???+??1??(??)-??2??(??),(1.1) where ??(??) and ??(??) represent biomass densities of the same fish population inside the free fishing zone, and reserve area, respectively, at a time ??. Let the fish subpopulation of the free fishing zone migrate into reserve area at a rate ??1 and the fish subpopulation of the reserve area migrate into free fishing zone at a rate ??2. ?? is the catch ability coefficient of the fish population in free fishing zone, and ?? denotes the effort devoted to the harvesting. Kar and Misra [18] modified the model proposed in [17] by Dubey et al.

On the other hand, regulation of exploitation of biological resources has become a problem of major concern nowadays in view of the dwindling resource stocks and the deteriorating environment. Exploitation of marine fisheries naturally involves the problems of law enforcement. Several governing instruments are suggested for the choice of a regulatory control variable. These are imposition of taxes and license fees, leasing of property rights, seasonal harvesting, direct control, and so forth. Various issues associated with the choice of an optimal governing instrument and its enforcement in fishery were discussed by Anderson and Lee [19]. Pradhan and Chaudhuri [20] studied a dynamic reaction model of a fishery. A regulatory agency controls exploitation of the fishery by imposing a tax per unit biomass of the landed fish. Ganguli and Chaudhuri [21] studied the bionomic exploitation of single species fishery using taxation as a control variable. Pradhan and Chaudhuri [22] studied a dynamic reaction model of two-species fishery with taxation as a control instrument. It deals with a dynamic reaction model of a fishery consisting of two competing species, each of which obeys the logistic law of growth. A regulatory agency controls exploitation of the fishery by imposing a tax per unit biomass of the landed fish. Kar [23] discussed a problem of selective harvesting in a ratio-dependent predator-prey fishery in which both the predator and prey obey the logistic law of growth. Kar [24] also studied a ratio-dependent prey-predator model with selective harvesting of prey species. A regulatory agency controls exploitation by imposing a tax per unit biomass of the prey species.

In order to keep a sustainable fishing resource. We will take some actions in fishing areas to protect certain fish stocks by restricting the fishermen’s fishing action. Such restriction would be implemented in the form of taxation. Following [22], we take ?? as dynamic (i.e., time-dependent) governed variable. We assume that the fish population in the free fishing zone is subject to a harvesting effort governed by the differential equations. To conserve fish populations, the regulatory agency imposes a tax ??>0 per unit biomass of the landed fish. (??<0 denotes the subsidies given to the fishermen). Keeping these aspects in view, the dynamics of the system may be governed by the following system of equations:????(??)=????(??)1-??(??)???-??1??(??)+??2????(??)-????(??)??(??),??(??)=????(??)1-??(??)???+??1??(??)-??2????(??),?[](??)={??????(??-??)??(??)-??-??}??(??),??(0)>0,??(0)>0,??(0)>0,0<??=1,0=??<1,(1.2) where ?? is constant rate of depreciation of capital. ?? and ?? are coefficients of proportionality, ???? is called the stiffness parameter measuring the strength of reaction of effort to the perceived rent; for further details about the biological sense of ?? and ??, we refer to [23, 24]. ?? is the constant price per unit biomass of the fish species, and ?? is the constant cost per unit of harvesting effort.

From [17], we know that there is no migration of fish population from reserve area to free fishing zone (i.e., ??2=0) and ??-??1<0, then ??(??)<0. Similarly, if there is no migration of fish population from free fishing zone to reserve area (i.e., ??1=0) and ??-??2<0, then ??(??)<0. Hence throughout our analysis, we assume that??-??1>0,??-??2>0.(1.3)

The aim of this paper is to find a proper taxation policy which would give the best possible benefit through harvesting to community while preventing the extinction of the fishing species. The structure of this paper is as follows. In the next section, we study the steady-state existence of positive equilibrium. In Section 3, by analyzing the corresponding characteristic equations, we discuss the local stability of equilibria. In Section 4, the global stability of the system is discussed by constructing a suitable Lyapunov function. In Section 5, optimal tax policy is discussed using Pontryagin’s maximum principle [25]. In Section 6, we try to interpret our results by numerical simulation. At last, we give some discussions.

2. Existence of Equilibria

We find the steady-states of (1.2) by equating the derivatives on the left-hand sides to zero and solving the resulting algebraic equations. This gives three possible steady states, namely: ??0(0,0,0), ??1(??,??,0), and ??2(??*,??*,??*). Existence of ??0(0,0,0) is obviously. We will consider the existence of ??1(??,??,0) in the following. We know that ?? and ?? are the positive solutions of the following equations:??2??(??)=????2(??)??-???-??1?????(??),1?????(??)=2?-????(??)+????2(??)??.(2.1)

From [17], we get a cubic equation in ??(??) as????3(??)+????2(??)+????(??)+??=0,(2.2) where??=????2????22??2,???=-2??????-??1?????22??,?????=??-??1?2????22-???-??2???????2,???=??-??2???2???-??1?-??1.(2.3)

Please note that there may be many possibilities that the above equation has positive solutions. We assume that the following inequalities hold:?????-??1?2????2<???-??2?????,???-??2????-??1?<??1??2,(2.4) then above equation has a unique positive solution ??(??)=??.

Knowing the value of ??, the value of ?? can then be computed from (2.1). It may be noted here that for ?? to be positive, we must have??????>??-??1.(2.5) Then we have following theorem.

Theorem 2.1. If inequalities (2.4) and (2.5) are satisfied, then (1.2) has a nonnegative equilibrium ??1(??,??,0).

Next, we will consider the existence of the positive equilibrium ??2(??*,??*,??*), where ??*, ??*, and ??* are positive solutions of?????(??)1-??(??)???-??1??(??)+??2???(??)-????(??)??(??)=0,????(??)1-??(??)???+??1??(??)-??2??(??)=0,????{??(??-??)??(??)-??}-??=0.(2.6)

From (2.6), we obtain??*=??????+??????(??-??)??>0(2.7) for????>??,(2.8)*=???-??2?+????-??2?2+4(??/??)??1??*2??/??>0.(2.9) Again, if???-??1?????????(??-??)??(??????+??)>1,(2.10) then??*=1????*????-????*??-??1???*+??2??*?>0.(2.11)

Thus we have following theorem.

Theorem 2.2. If (2.8) and (2.10) are satisfied, then (1.2) has a unique interior equilibrium point ??2(??*,??*,??*).

3. Local Stability Analysis

We first consider the local stability of equilibria. The variational matrix of the system (1.2) is?????????(??,??,??)=??-2????(??)??-??1-????(??)??2??-????(??)1??-2????(??)??-??20???????????????(??)(??-??)0????{??(??-??)??(??)-??}-??.(3.1)

At ??0(0,0,0), the characteristic equation of ??0(0,0,0) is???1??-????2+??2??+??3?=0,(3.2) where??1??=-(??????+??),2?=-??-??1+??-??2?,??3=???-??1????-??2?-??1??2.(3.3)

In this cubic equation, one root is ??1<0, the sum of other two roots is -??2>0. So (3.2) at least has one positive root. Therefore, ??0 is unstable. We get the following theorem.

Theorem 3.1. The equilibrium ??0 of (1.2) is unstable.

At ??1(??,??,0), the characteristic equation of ??1(??,??,0) is???1??-????2+??2??+??3?=0,(3.4) where??1=??????(??-??)????-(??+??????),2???=-1????+??????+??2????+?????????<0,3=????2??2+??????2??1??????+????????????>0.(3.5) Similarly, one of the eigenvalues of the variational matrix ??1(??,??,0) is ??????(??-??)??-(??????+??). This eigenvalue is positive or negative according to whether ??<(1/????)(??????-(??????+??)/????) or ??>(1/????)(??????-(??????+??)/????). The sum of other two roots is -??2>0; the product of other two roots is ??3>0, so (3.4) has two positive roots. Therefore, ??1 is unstable. Then we have following theorem.

Theorem 3.2. If inequalities (2.4) and (2.5) are satisfied, then the equilibrium ??1 of (1.2) is unstable.

To determine the local stability character of the interior equilibrium ??2(??*,??*,??*), we compute the variational ??2(??*,??*,??*) about (??*,??*,??*)??2???*,??*,??*?=?????????-2????*??-??1-????*??2-????*??1??-2????*??-??20????????*???????(??-??)00.(3.6) The characteristic equation of the variational matrix ??2(??*,??*,??*) is given by??3+??1??2+??2??+??3=0,(3.7) where??1=????*??+??2??*??*+????*??+??1??*??*??>0,2=?????*??+??2??*??*??????*??+??1??*??*?+??2??*??????*(??-??)-??1??2,??3=?????*??+??1??*??*?(??-??)??2??????*??*>0.(3.8) By the Routh-Hurwitz criterion, it follows that all eigenvalues of (3.7) have negative real parts if and only if??1>0,??3>0,??1??2>??3.(3.9) Here, ??1>0, ??3>0, and??1??2-??3=??????*??+??2??*??*??????*??+??1??*??*?-??1??2??????*??+??1??*??*+????*??+??2??*??*?+?????*??+??2??*??*???2??*??????*=?(??-??)??????*??*????+??1????*2????*+??2????*2????*??????*??+??1??*??*+????*??+??2??*??*?+?????*??+??2??*??*???2??*??????*(??-??)>0.(3.10) Hence ??2(??*,??*,??*) is locally asymptotically stable. We get the following theorem.

Theorem 3.3. If (2.8) and (2.10) are satisfied, then the unique interior equilibrium point ??2(??*,??*,??*) of (1.2) is locally asymptotically stable.

4. Global Stability

In this section, we will consider the global stability of the unique interior equilibrium of system (1.2) by constructing a suitable Lyapunov function. We have the following theorem.

Theorem 4.1. If (2.8) and (2.10) are satisfied, then the unique interior equilibrium point ??2(??*,??*,??*) of (1.2) is globally asymptotically stable.

Proof. Define a Lyapunov function ????(??(??),??(??),??(??))=??(??)-??*?-??*ln??(??)??*?+??1????(??)-??*?-??*ln??(??)??*?+??2????(??)-??*?-??*ln??(??)??*?,(4.1) where ??1 and ??2 are positive constants to be chosen suitably in the subsequent steps. It can be easily verified that ??(??(??),??(??),??(??)) is zero at the equilibrium point and positive for all other positive values of ??(??), ??(??), and ??(??).
Differentiating ?? with respect to ?? along the solutions of (1.2), a little algebraic manipulation yields????=???????(??)-??*???(??)????(??)????+??1???(??)-??*???(??)????(??)????+??2???(??)-??*???(??)????(??).????(4.2) Choosing ??1=??2??*/??1??*, ??2=1/????(??-??), a little algebraic manipulation yields ??????????=-?????(??)-??*?2-????2??*????1??*???(??)-??*?2-??2??*?????(??)??(??)*??(??)-??*???(??)2<0.(4.3) So ??2(??*,??*,??*) is globally asymptotically stable. The proof is complete.

5. Optimal Harvest Policy

The objective of the regulatory agency is to maximize the total discounted net revenues that the society drives due to the harvesting activity. Symbolically, this objective amounts to maximizing the present value ?? of a continuous time stream of revenues given by???=80??-????(??????(??)-??)??(??)????,(5.1) where ?? denotes the instantaneous annual rate of discount.

Our objective is to determine a tax policy ??=??(??) to maximize ?? subject to the state (1.2) and the control constraint??min=??(??)=??max(5.2) on the control variable ??(??).

We apply Pontryagin’s maximum principle in Burghes and Graham [26] to obtain the optimal equilibrium solution to this control problem. The Hamiltonian of this control problem is??=??-????(??????(??)-??)??(??)+??1??(??)????(??)1-??(??)???-??1??(??)+??2???(??)-????(??)??(??)+??2??(??)????(??)1-??(??)???+??1??(??)-??2???(??)+??3([]??){??????(??-??)??(??)-????(??)-????(??)},(5.3) where ??1(??), ??2(??), and ??3(??) are adjoint variables. Hamiltonian ?? must be maximized for ??(??)?[??min,??max]. Assuming that the control constraints are not binding (i.e., the optimal solution does not occur at ??(??)=??min or ??max), we have singular control [27] given by????????=??3(??)(-????????(??)??(??))=0???3(??)=0.(5.4)

Now the adjoint equations are????1(??)????=-????????=-??-??????????(??)-??1?(??)??-2????(??)??-??1?-????(??)-??2(??)??1,(5.5)????2(??)????=-????????=-??1(??)??2-??2?(??)??-2????(??)??-??2?,(5.6)????3(??)????=-???????????=--????(??????(??)-??)-??1?(??)????(??)=0.(5.7) From (5.7), we obtain??1(??)=??-?????????-?????(??).(5.8) To obtain an optimal equilibrium solution, we shall rewrite (5.6) by considering the interior equilibrium as????2(??)????=??2??-????+??1??2(??),(5.9) where??1=2????*??+??2-??,??2=???????*???-??2.(5.10) The solution of this linear equation is??2(??)=??0????1??-??2??1??+??-????.(5.11) The shadow price (the term “shadow price” refers to the fact that the asset’s value is not its direct sale value but the value imputed from its future productivity [4]) ??2(??)?????? is bounded as ???8 if ??0=0. Then??2??(??)=-2??1??+??-????.(5.12) Similarly,??1??(??)=-2??1??+??-????,(5.13) where??1=2????*??+??1+????*-??,??2=??1??2??1+??-??????*.(5.14)

Substituting the value of ??1(??) from (5.8) into (5.13), we get????-????*??=-2??1.+??(5.15)

Now using the values of ??*, ??*, and ??* from Section 3 into (5.15), we get an equation for ??(??). Let ???? be a solution (if it exists) of this equation. Using this value of ??(??)=????, we get the optimal equilibrium point (????,????,????).

The existence of an optimal equilibrium solution has been created, which satisfies the necessary conditions of the maximum principle. As stated by Clark [3], an optimal approach path, which is composed of the combination of bang-bang control and nonequilibrium singular controls, is fairly difficult to find. Faced with the same difficulty, Clark [3] researched a simple model of two ecologically independent fish populations. The present model is far more complex than the model mentioned by Clark. So we only consider an optimal equilibrium.

From the above analysis carried out in this section, we observe the following.(1)From (5.4), (5.12), and (5.13), we note that ????(??)??????(??=1,2) is independent of time in an optimum equilibrium. Hence they satisfy the transversality condition at 8, that is, they remain bounded as ???8.(2)Considering the interior equilibrium, (5.7) can be written as ??1????*=???????*???-??-????=??-????????.????(5.16) This implies that the total users cost of harvest per unit effort is equal to the discounted values of the future price at the steady-state effort level.(3)From (5.15), we get ??????*??-??=-2????*??1+???0as???8.(5.17) Thus, the net economic revenue ??(??8,??8,??8,??)=0.

This shows that an infinite discount rate results in the complete dissipation of economic revenue. For zero discount rate, it is indicated that the present value of continuous time stream gains its maximum value.

6. Numerical Example

In this section, we use Matlab 7.0 to simulate a numerical example to illustrate our results.

Let ??=1.5, ??=2.7, ??=200, ??=250, ??1=0.5, ??2=0.4, ??=0.2, ??=10, ??=10, ??=0.8, ??=0.4, ??=19, and ??=0.03, in appropriate units.

Then for the above values of the parameter, optimal tax becomes ????=15.1 and corresponding stable optimal equilibrium is (52.88, 223.90, 11.48). The time-paths of the free fishing zone ??(??), reserve area ??(??), and effort ??(??) are shown in Figure 1, and the three dimensional phase space portrait is depicted in Figure 2. From the solution curves, we infer that the system is globally stable about the interior equilibrium point. Biological meaning is that this optimal tax system achieves relatively balance and all the fish population may be maintained at an appropriate equilibrium level.

In Figures 35, variation of free fishing zone, reserve area, and harvesting effort against time are plotted for different tax levels. From these plots, we observe that as the rate of tax increases, free fishing zone populations and reserve area populations increase while harvesting effort decreases as expected.

We take ??=0 and ??=15.1. Now taking the tax ??=0 in (2.6), we have ??*=10.86, ??*=215.29, and ??*=44.26, respectively. But for the optimal tax ??=15.1, we have ??*=52.88, ??*=223.90, and ??*=11.48.

From these results, it is clear that if the fishermen have to pay no tax, then they use a large amount of effort compared to the case when the fishermen have to pay the optimal tax. As a result, the steady-state values of the two species for the case of ??=0 are much less than those for the case of ??=15.1.

The computer analyzed results for the time course display of the two species ?? and ?? and the phase space trajectory for ??=0 and ??=15.1 using these parameter values is also shown in Figures 35.

By comparing Figures 1 and 6, we can see that, in the model without reserve area, there exists oscillation of fish population in free fishing zone for a long time when they tend to positive equilibrium. However, in our model, fish population tend to positive equilibrium quickly in reserved zone, which is beneficial for population protection.

7. Discussion

In this paper, we study an optimal harvesting problem for fishery resource with prey dispersal in a two-patch environment: one is a free fishing zone and the other is a reserve zone, focusing attention on the use of taxation as an optimal governing instrument to control exploitation of the fishery. In Sections 3 and 4, we have discussed the local and global stability of the system. It has been observed that, in the case of no taxation, even under continuous harvesting in the free fishing zone, the fish population may be maintained at an appropriate equilibrium level. On the other hand, in the case of taxation, population may be also sustained at an appropriate equilibrium level. But from Figures 3 and 4, we know that as the rate of tax increases free fishing zone populations and reserved area populations increase while harvesting effort decreases. This situation is in accordance with reality. In the model without reserve area, there exists oscillation of fish population in free fishing zone for a long time when they tend to positive equilibrium. However, in our model, fish population tend to positive equilibrium quickly in reserved zone, which is beneficial for population protection.

Acknowledgments

This work was partially supported by the NNSF of China (10961018), the NSF of Gansu Province of China (1107RJZA088), the NSF for Distinguished Young Scholars of Gansu Province of China (2011GS03829), the Special Fund for the Basic Requirements in the Research of University of Gansu Province of China and the Development Program for HongLiu Distinguished Young Scholars in Lanzhou University of Technology.