Cooperative Game for Fish Harvesting and Pollution Control

This paper studies fishery strategies in lakes, seas, and shallow rivers subject to agricultural and industrial pollution. The flowing pollutants are modeled by a nonlinear differential equation in a general manner. The logistic growth model for the fish population is modified to cover the pollution impact on the fish growth rate. We start by presenting the stability analysis of the dynamical system to discern the different types of the evolution of the fish population according to human actions. A cooperative game is formulated to design strategies for preserving the fish population by controlling the pollution as well as the fish stock for harvesting. The sufficient conditions for implementing the cooperative strategy are investigated through an incentive design approach with an adaptive taxation policy for the players. Numerical results are presented to illustrate the benefit of the cooperative for fish population preservation but also for the players’ rewards.


Introduction
Prior to the 2019 coronavirus disease (COVID-19) pandemic, towns and cities had grown and attracted more than fifty percent of the world population. Their sustainability and competitive efficiency have become increasingly important issues to decision-makers. A city, more simply, is a geographic area where live people, with hard infrastructures and changing systems that interconnect the city with itself. Because of pandemics, economic drivers, energy, climate change, and resource supply, the concept of biodiversity preservation has gained more attention for the next generation of cities, called smart cities. The "smart city" can include a large and extensive range of systems, networks, and infrastructure elements. However, in this paper we use the term "smart city" to refer to more intelligent, more efficient, and more sustainable cities, and we focus on pollution control and economic supply. In support of the environment, risk management and sustainable development; natural resource management including systems for the reduction of pollutants; increasing energy efficiency; managing the human response to environmental stresses, and sustaining biodiversity are key elements of smart cities.
Seas and lakes provide a variety of economic uses: they are sources of drinking water, fishing, and recreational sports, and they provide pleasant locations for smart homes and drainage for agriculture. Runoff from fields flows into the streams and rivers that feed the lake. Much of the agricultural runoff includes phosphorus from fertilizers and animal wastes. Phosphorus is the primary nutrient of algae and weeds in the water. When it grows excessively from an infusion of phosphorus, algae blooms reduce the oxygen content of the lake and release toxins into the water. In addition, large amounts of waste produced by some industries and municipal sewage enter natural water bodies, either without treatment or through inappropriate and inadequate treatment processes [1]. In fact, the sources of water pollution are various [2]. Seas and lakes are polluted with high concentrations of phosphorus, Pb, Zn, Cu, Cd, and other toxins in their sediments. This pollution severely affects aquatic life and leads to a serious environmental, climatic, and public health problem. Fish and other marine species die due to high levels of pesticides and toxins in water and other types of water pollution, such as plastic and litter [3]. On the other hand, fish reproduction is seriously affected because of genetic mutations when toxins are absorbed. Fish are therefore mainly affected by human nuisance and it is necessary to give sufficient attention to this issue and to implement the necessary corrective measures [3,4].
Research activities are focused on finding suitable solutions to these problems. Anupam et al. [5] proposed a new mathematical model that uses fuzzy inferences to study the impacts of global warming, water pollution, and juvenile fish harvesting on the production of mature Hilsa fish. They used the Mamdani inference method and applied it to a model based on fuzzy rules. Most recent mathematical models take into account the spatialization of fisheries to study the control of a multi-site fishery [6]. For this, bio-economic models can be used. These models take into account economic aspects and, in particular, changes in investment and prices, in fisheries management models [7]. It is possible to eliminate overfishing and move towards a sustainable fishery without risk of extinction for the exploited species, by varying the number of fishing sites and the costs of operating the fishery [8]. For a more recent work on economic models, see [9], where the authors studied the dynamics of marine industrial interactions through a dynamical systems theory approach.
For pollution control, the authors in [10] considered a dynamic partial equilibrium model that combines the optimum exploitation of renewable resources with optimum pollution control. For them, the pollution accumulates as a slowly decomposing stock and is supposed to affect the growth and quality of the stock of renewable resources. They maximized a social welfare function that gives the present value of the difference between the benefits of natural resources and the costs of controlling pollution. Their analysis also gives a general result concerning the state of equilibrium of two problems of optimal control. In [11], a mathematical model for solving the dispersion of pollutants in a river is considered. The concentration of the pollutant is obtained using a finite element method. This is subject to optimal control of the water treatment plants in order to achieve minimum costs.
In this paper, we introduce a mathematical model describing the impact of pollution on the fish population. The pollution dynamic covers the time unit input of chemical residuals and the eutrophication process; see [12]. The dynamic of the fish population is obtained by discharging the harvested fish stock from the natural logistic growth model. The inwater part of the allowed portion, according to the regulations, of the fishable population is referred to as the "fish stock" [13]. To gain clarity in the mathematical analysis, we first study the properties of the eutrophication dynamics. Moreover, considering fish preservation, we investigate the gain optimization for polluters and fishermen in noncooperative and cooperative regimes. Two separate payoff functions are defined in the non-cooperative regime, while a single aggregated payoff for all the players is considered in the cooperative regime under the authority of a third entity. It is revealed that the cooperative regime is more supportive of fish preservation.
The outline of the paper is as follows: in Section 2, we discuss the management problem to develop the mathematical model. Section 3 is dedicated to the mathematical analysis, which is based on the basic reproduction number of the eutrophication dynamics. The non-cooperative optimization between the polluters and the fishermen with the constraint of preserving the fish population is presented in Section 4. In Section 5, we address the cooperative optimization and present sufficient conditions for its implementation.

Water Pollution and Fishing Industry
In this section, we describe the management problem of ecosystem preservation, present the mathematical models for the pollution and the fish population, and identify the control variables.

Management Problema
We identify the first class of players: farmers and mining industries surrounding rivers, lakes, and seas are polluters. Industry companies need to clean up material or deposit chemical waste into the water. Farmers use artificial fertilizers to grow crops; these fertilizers contain chemical elements such as phosphorus and zinc, which are ultimately washed out into the water. The dynamic of the chemical waste (such as the amount of phosphorus in the water) is given by [12,14]: Here, n is the total number of polluters (farming and industrial companies); a j is the input of chemical waste (due to industrial production and farming) produced by the player j; c ě 0 is proportional to the rate of loss of chemical elements due to the sedimentation, the outflow, and the sequestration in other biomasses. The term b x p 1 x p 1`m p (p ě 2, b ą 0) models the production of chemicals in the water because the dissolution process generates other types of chemicals through reactions between water components and the spilled chemical.
Fishes represent the second class of players. To illustrate the dynamics and the impact of the fishing industry, we consider that the fish population size evolves according to: where the first term h 1 is the inflow and the second term h 2 is the outflow due to harvesting activities. We agree that the function h 2 is essentially the harvesting activity as the natural death rate is neglected and we interpret the term h 1 as the drift of the logistic growth model [6], i.e., The carrying capacity K is non-increasing with the pollution level x 1 , while the natural growth rate r is taken to be constant. If h 1 ą h 2 , then the fish population growth continues, and for h 1 ă h 2 , the harvesting activities are at over-fishing levels that lead to the extinction of the fish population.

Mathematical Model
The resulting ecosystem from the above problem consists of two predators (polluters and fishermen) and a single prey (fish population). The polluters are indirect predators because they do not have a direct interest in the fish population, but their activities are harmful for the fish population and can cause the fish population's extinction. The fishermen are the direct predators for the fish population but they are naturally concerned with its preservation. If the watercourse is polluted, the fish population drops to the lowest level through Kpx 1 q. The fish industry will lose considerable profit because the harvesting cost becomes higher and the utility (the price of fish in the market) smaller due to the low quality. Thus, a smarter solution, addressed in this paper, is to design a policy for the fish industry that may obtain a better outcome by cooperating with the farmers to control the pollution level. In this way, the farmers may be compensated by a tax reduction.
We consider h 2 " x 2 2 u 2 , where u 2 measures the fish stock to be harvested. Accordingly, the dynamics of the fishery subject to industrial and agricultural pollution is: where u 1 ptq " ř n j"1 a j ptq globalizes the actions of the n polluters to lighten the notation. Then, the dynamical system (1) can be recast in the condensed form: 9 xptq " f pxptq, uptqq, starting at xp0q " x 0 . ( where the vector u " pu 1 , u 2 q describes the players' activities. Short-term suspension of the players' activities appears to be insufficient to restore the fish population because, for u 1 " u 2 " 0, the residual chemicals continue to jeopardize the natural growth of the fish population through eutrophication.

Analysis of the Model
This section is devoted to the stability analysis of the system (1). We start by addressing the stability of the eutrophication dynamics in a more general setting, which facilitates the stability analysis of the system, which is considered without chemical input and fishery activities, i.e., u 1 " u 2 " 0, x 2 ptq " rx 2 ptq´r x 2 2 ptq Kpx 1 ptqq , For this, we primarily focus on the analysis of the eutrophication process, which is the first equation of (3), and secondly, we end with the analysis of the complete dynamics.

Stability Analysis of the Eutrophication Process
Here, we consider the dynamics of the pollutant, described by the following nonlinear equation: x 1 ptq " kpx 1 ptqqx 1 ptq, with m ą 0 and the integer p ě 2. Remark that the mapping x 1 Þ Ñ kpx 1 qx 1`c x 1 is nonnegative and increasing for x 1 P R`, which implies that the semi-flow generated by (4) is monotone with respect to the initial conditions in R`. In this sequel, we will denote by t Þ Ñ Xpt, x 0 1 q the unique solution of the Cauchy problem (4) with data x 0 1 . The monotonicity of (4), for two initial conditions, x a and x b , is given by: 0 ď x a ď x b ùñ Xpt, x a q ď Xpt, x b q, @t ě 0.
To establish the causality of the function k on the existence of the equilibrium points and their asymptotic behaviors, we look at the variation of k and find that k is increasing on r0, x ‹ 1 s and decreasing on rx while kp0q "´c and lim We now introduce the basic reproduction number R 0 of the chemical eutrophication process: In contrast to the epidemiological approach of defining R 0 as an amplification factorsee [15,16]-we assign to R 0 the rate of generating new chemical components as a result of the reactions between the components during eutrophication. We note that kpx ‹ 1 q " cpR 0´1 q; thus, for R 0 ă 1, the total amount of generated chemicals during eutrophication is less than the number of involved input chemicals. If R 0 ą 1, the eutrophication tends to increase the amount and the number of chemicals, while R 0 " 1 means that the toxicity level, given by the amount of residual chemicals, stays constant because it is not affected by the eutrophication. We enumerate the equilibrium points of (4) according to R 0 . Proposition 1. The number of equilibrium points of the eutrophication dynamics (4) is given by R 0 as follows: If R 0 ą 1, additional to the free-pollution equilibrium x1 , there are two equilibrium points x ‚ 1 and x § 1 , which are positive. For p " 2, they are given by: It is clear that the pollution-free equilibrium x1 " 0 is always a stationary solution of (4). Furthermore, the other equilibrium of (4) is obtained by solving kpx 1 q " 0 in R. Indeed, using the variation of k stated earlier in this section, we observe that the graph of k intersects the abscissa x 1 " 0 once at x ‹ 1 when R 0 " 1 and twice at x ‚ 1 and x § 1 when R 0 ą 1. The following lemmas concern the behavior of the equilibrium points given in Proposition 1. In fact, Lemma 2 states that when R 0 ă 1, the pollution-free equilibrium x1 is globally asymptotically stable. The case R 0 " 1 is given in Lemma 3. In this latter case, the dynamics converges to x1 or x ‹ 1 whether the initial condition is in p0, x ‹ 1 q or in px ‹ 1 ,`8q. When R 0 ą 1, we have in Lemma 4 that x ‚ 1 is a uniform repeller. In addition, the dynamics converges to x1 or x § 1 whether the initial condition is in p0, x ‚ 1 q or in px ‚ 1 ,`8q.

Lemma 2.
If R 0 ă 1, then the pollution-free equilibrium point is globally asymptotically stable. Moreover, for initial state x 0 1 ą 0, the mapping t Þ Ñ Xpt, x 0 q is exponentially decreasing at the rate of cp1´R 0 q, i.e., Proof. Let us consider a positive initial condition x 0 1 ą 0; we use the monotonicity property (5) with x a " 0 and x b " x 0 1 to show that the solution t Þ Ñ Xpt, x 0 1 q satisfies: Xpt, 0q ď Xpt, x 0 1 q ùñ 0 ď Xpt, x 0 1 q, @t ě 0, since 0 is an equilibrium point.
Lemma 3. If R 0 " 1, the notion of stability is weakened because there is no basin of attraction covering the neighborhood of any two equilibrium points. However, the following properties hold: Proof. We recall that when R 0 " 1, the global maximum of k is attained at x ‹ 1 , which is also an equilibrium of the dynamical system (4), i.e., For an initial condition x 0 1 Ps0, x ‹ 1 s, the monotonicity property (5) applied to x a " x1 , which, because the function k is increasing in R`, implies that: Multiplying by 0 ď Xpt, x 0 1 q ď Xpt, x ‹ 1 q, we find that: which implies that: This shows that the mapping t Þ Ñ Xpt, x 0 1 q is decreasing, and since it is bounded from below by the pollution-free equilibrium x1 , we have Xpt, x 0 1 q Ñ x1 " 0 as t Ñ`8. For an initial condition greater than and since k is decreasing in px ‹ 1 ,`8q, we obtain: Multiplying by 0 ă x ‹ 1 ă Xpt, x 0 1 q, we obtain: which means t Þ Ñ Xpt, x 0 1 q is decreasing since it is bounded from below by the equilibrium x ‹ 1 . We conclude that Xpt, x 0 1 q Ñ x ‹ 1 as t Ñ`8.
It is worth mentioning that the fact that x ‚ 1 is a repeller, for R 0 ą 1, is a consequence of (1) and (2). It is therefore sufficient to prove (1) and (2). However, the proof of (1) is similar to that of Lemma 3; therefore, we elaborate only the proof of (2). We will do so stepwise, by considering first the initial conditions in rx ‚ 1 , x § 1 s and then those in px § 1 ,`8q.
Proof. Note that the global maximum of k, given by Using the fact that the function k is increasing in px ‚ 1 , x ‹ 1 q, we obtain: By arguing similarly for x 0 1 P rx ‹ 1 , x § 1 s and using the fact that the function k is decreasing Multiplying the terms of the previous inequality by 0 ă Xpt, x 0 1 q, we obtain in both cases: Xpt, x 0 1 q, @t ě 0, meaning that t Þ Ñ Xpt, x 0 1 q is increasing since it is bounded from above by x § 1 . This implies that Xpt, x 0 1 q Ñ x § 1 as t Ñ`8. For an initial condition greater than x § 1 , x 0 1 ą x § 1 , it yields from the monotonicity condition (5) with x a " x § 1 , x b " x 0 1 that: Since k is decreasing in px § 1 ,`8q Ă px ‹ 1 ,`8q, we obtain: Hence, t Þ Ñ Xpt, x 0 1 q is decreasing and bounded from below by the equilibrium

Stability Analysis of the Fish Pollution Dynamic
Here, we focus on the analysis of the dynamical system (3), i.e., where K is a non-increasing function of x 1 ptq. The following results give the equilibrium points related to the fish pollution dynamic. These points are obtained by combining the equilibrium points of the eutrophication process (Section 3.1) with the equilibrium points of the second equation of Equation (8), which are x2 " 0 andx 2 " Proposition 5. The number of equilibrium points of the fish pollution dynamics (8) depends on R 0 as follows: • If R 0 ă 1, the only equilibrium points are the trivial equilibrium and the best-case scenario equilibrium, which are given, respectively, by: E 0 " px1 , 0q " p0, 0q and E 1 " px1 , Kpx1 qq " p0, 1q.
• If R 0 " 1, additional to the trivial equilibrium and the best-case scenario equilibrium, there are two equilibrium points, which are: • If R 0 ą 1, additional to the trivial equilibrium and the best-case scenario equilibrium, there are four equilibrium points, which are: Here, x1 , x ‹ 1 , x ‚ 1 and x § 1 are given in Proposition (1).
Since the second equation of (8) is a Bernoulli differential equation, we consider the following Lemma, which will be used together with the results of Section 3.1 to analyze the fish pollution dynamics.

Lemma 6.
Let v P Cpr0,`8q, R`q be a given continuous function and w P C 1 pr0,`8q, R`q the solution of the following Bernoulli differential equation: wptq " rwptq´r Kpvptqq w 2 ptq, t ą 0 wp0q " w 0 P R`.
If v is a monotone function with vptq Ñ v˚as t Ñ`8, then, for each w 0 ą 0, we have: Proof. The solution to the Bernoulli differential equation is: We recall that K is a non-increasing function of vptq. As vptq Ñ v˚, Kpvptqq Ñ Kpv˚q. Let ą 0 be given and fixed. Then, there exists t 0 ě 0 such that Kpv˚q ď Kpvptqq ă Kpv˚q` for all t ě t 0 .
Hence, for all t ě t 0 , we have: Since is arbitrary, we obtain that wptq Ñ Kpv˚q as t Ñ`8.
We have now all the technical elements in order to completely describe the asymptotic behavior of the solutions of (8). In the following, we refer to E 2 , E 4 and E 6 as extinction-fish equilibrium and E 1 , E 3 , E 5 and E 7 as persistent-fish equilibrium. Furthermore, E 7 will be called the worst-case scenario equilibrium.
Next, we list in a series of lemmas the results concerning the asymptotic behavior of the system (8) in terms of the threshold R 0 . Lemma 8 states that when R 0 ă 1, the best-case scenario equilibrium E 1 is globally asymptotically stable. When R 0 " 1 (Lemma 8), the extinction-fish equilibrium E 2 is unstable while the dynamics of the system converges either to the best-case scenario equilibrium E 1 or to the persistent-fish equilibrium E 3 depending on the initial condition. For the case R 0 ą 1 (Lemma 9), the extinction-fish equilibrium E 4 and E 6 are unstable. In this case, the dynamics of the system converges either to the best-case scenario equilibrium E 1 or to the worst-case scenario equilibrium E 7 depending also on the initial condition.

Lemma 7.
If px 0 1 , x 0 2 q P R`ˆR˚, then the fish is uniformly persistent. More precisely, there exists δ ą 0 such that for each px 0 1 , x 0 2 q P R`ˆR˚, we have: Proof. Let px 0 1 , x 0 2 q P R`ˆR˚be given. Using the results of Section 3.1, we know that, independently of R 0 , t Ñ x 1 ptq is monotone and converges to some κ, where κ is either x1 , x ‹ 1 , x1 or x § 1 . Therefore, we deduce from Lemma 6 that x 2 ptq Ñ Kpκq as t Ñ`8. The proof is completed. Lemma 8. If R 0 ă 1, then the best-case scenario equilibrium E 1 is globally asymptotically stable in R`ˆR˚. More precisely, for each px 0 1 , x 0 2 q P R`ˆR˚, we have: Proof. Let px 0 1 , x 0 2 q P R`ˆR˚be given. If R 0 ă 1, then, from Lemma 2, t Þ Ñ x 1 ptq is monotone and converges to 0 as t Ñ`8. Thus, Lemma 6 implies that x 2 ptq Ñ Kp0q " 1 as t Ñ`8. Lemma 9. If R 0 " 1, then the extinction-fish equilibrium E 2 " px ‹ 1 , 0q is unstable and we have the following properties: (i) For each px 0 1 , x 0 2 q P R`ˆR˚with x 0 1 P p0, x ‹ 1 q, we have: i.e., t Þ Ñ px 1 ptq, x 2 ptqq converges to the best-case scenario equilibrium E 1 .
Proof. The fact that E 2 " px ‹ 1 , 0q is unstable is a consequence of (i) and (ii). The proof of (i) follows the same lines of the proof of Lemma 8 and thus we only prove (ii). Let px 0 1 , x 0 2 q P R`ˆR˚with x 0 1 P rx ‹ 1 ,`8q be given. Then, Lemma 3 ensures that t Þ Ñ x 1 ptq is monotone and converges to x ‹ 1 as t Ñ`8. Hence, we infer from Lemma 6 that: The proof of the next lemma uses the same arguments of the proofs of Lemmas 7-9.

Numerical Illustrations
For illustration, we consider the capacity function Kpx 1 q " 1 x 1`1 and the following values : p " 2, m " 1, c " 1, and r " 1. Therefore, the value of R 0 will depend on b. For b " 0.7, 2 and 3, we have R 0 " 0.35 ă 1, R 0 " 1 and R 0 " 1.5 ą 1, respectively. The dynamic of the model is given in Figure 1.  As stated in Proposition 5, the number of equilibrium points depends on R 0 . • R 0 ă 1. We have the trivial equilibrium E 0 " p0, 0q and the best-case scenario equilibrium E 1 " p0, 1q, which is a nodal sink and globally asymptotically stable. It characterizes non-pollution for a perfect natural growth of the fish population. • R 0 " 1. Here, additionally to the trivial equilibrium E 0 and the best-case scenario equilibrium E 1 , we have E 2 " px ‹ 1 , 0q " p1, 0q, which is unstable, and E 3 " px ‹ 1 , Kpx ‹ 1 qq " p1, 1{2q. In this case, we see that if the initial condition x 0 1 is in p0, 1q, the dynamic converges to E 1 as stated by Lemma 9. For initial condition x 0 1 starting at x 1 " 1, E 3 is a stable point and behaves as a nodal sink. The level of pollution is damping at half the level of fish population growth. • R 0 ą 1. Additionally to the trivial equilibrium E 0 and best-case scenario equilibrium E 1 , we have E 4 " px ‚ 1 " 0.38, 0q, E 5 " px ‚ 1 " 0.38, Kpx ‚ 1 qq, E 6 " px § 1 " 2.62, 0q and E 7 " px § 1 " 2.62, Kpx § 1 qq. The equilibrium point E 5 is unstable (as well as E 4 and E 6 ) and behaves as a source. This stands for a very low pollution state, where the pollutants are not yet damping or at least not at an alarming level regarding the fish population growth. The stationary point E 7 is a nodal sink and is locally asymptotically stable. It illustrates the impact of a high level of pollution on the fish population. This leads asymptotically to the survival of the fish population after the total disintegration of the pollutant.
The main issue that we address in this paper is to design a fishing/pollution policy that satisfies all parties. The outcome is to have the best positioning of the stable (asymptotically stable) stationary point for both fishermen and polluters and to be able to bring and/or maintain the system state around this equilibrium. Before designing this stability policy, we must determine the feasibility, i.e., the controllability issue.

Controllability
In the rest of the paper, we consider u 1 and u 2 as control variables and we intend to investigate their optimal dynamics to preserve the fish population and to generate a consistent gain for both of the players in terms of business outcomes. Here, we look at the feasibility of designing such dynamics for u 1 and u 2 by investigating the controllability of the model. The feasibility of a mutually beneficial strategy for all the players agrees with the controllability of system (2). The first-order controllability around an equilibrium point implies its short-term controllability. In this sense, we announce the condition for the fish population preservation policy as a controllability result. Proposition 11. The regulation of the total pollution u 1 is a necessary and sufficient condition on the feasibility of the stability policy.
Proof. The linear tangent system of (2) around the equilibrium point px˚, u˚q is given by: where A " ∇ x f | x"x˚,u"u˚a nd B " pB 1 , B 2 q " ∇ u f | x"x˚,u"u˚.
We denote by C 1 , C 2 and C the Kalman controllability matrices of (11) with respect to u 1 , u 2 and u, respectively. They are given by: Since the state variable x i is positive at any time, then rank(C 1 ) = rank(C) = 2 while rank(C 2 ) = 1. Consequently, the system (1) is, in the short term, controllable by u 1 while not being locally controllable by u 2 .
The non-controllability of (1) by u 2 is due to the fact that the effect of this control is limited only to the fish population, while the control u 1 acts on both dynamics, directly on the pollution level and via the limit growth rate K of the fish population. In spite of this evidence of regulating the pollution in the river, it will be politically difficult to attribute full responsibility to industries and farmers. Therefore, a cooperative strategy to yield them for the preservation of the fish population is required.
In the next sections, we investigate the evolution of the fish population as the consequence of unilateral/bilateral strategies adopted by the players. The unilateral strategies yield from a non-cooperative situation while the bilateral strategies are led by a third entity that decides the threshold of fish stock to be harvested, the tolerated level of pollution, and the tax to be paid by the players.

Non-Cooperative Regime: Quest for Profit
In the non-cooperative situation, the polluters as well as the fishermen maximize their payoff with the constraint of preserving a certain threshold amount of fish. It is worth noticing that the no-cooperation regime is not as bad as the lack of a third regulation entity; it is rather the failure to take into account the fishery actions' when defining a strategy for the polluters and vice versa. Let J 1 and J 2 be the payoffs of the polluters (industries + farmers) and of the fishermen, respectively, The goal of the polluter (resp. of the fisher) is to maximize the payoff J i for i " 1 (resp. for i " 2). The control u 1 represents the amount of pollutant deposited in the river. This chemical residual from the industrial/agricultural process generates the utility U 1 . Moreover, it governs the production upstream, so the cost process C 1 depends on u 1 . The control u 2 is the rate of harvesting fish so that x 2 u 2 is the exact amount of harvested fish. This yields the utility U 2 , which might depend also on the consumer demand. The cost C 2 is associated with the fish stock of x 2 u 2 . It depends on x 2 in the sense that the more fish there are in the river, the lower is the cost. It is therefore suitable to deal with the following functions: The cost function of the polluter depends on the industrial/agricultural residual as well as on the fish population in the sense that it covers also the environmentalist tax, and the utility function U 1 would take into account also the consumer demand for the fish. Since x 2 is normalized, the cost C 2 declines when x 2 grows, and it grows with u 2 . With the above choice of utilities and costs, the control system objective for each player is the optimal u i that generates the maximum gain from the business.

Optimal Control System
We seek the optimal control pu 1 , u 2 q by applying the maximum principle of Pontryaginsee [17]-to the problems: subject to p1q subject to p1q.
The set of admissible control values for the pollution is A 1 " u 1 : 0 ď u 1 ď Tu max ( , where u max is the maximum amount of chemical residual generated, by all the polluters, per unit of time. As the fish population is normalized 0 ď x 2 ď 1, the control set A 2 for the fish population is given by A 2 " r0, 1s.
It yields the following form of the optimal control dynamics for pollution and fishery management: u 1 ptq "´pγ 1 u 1 ptq´α 1 q˜´c`2 bx 1 ptq px 2 1 ptq`1q 2¸, 9 u 2 ptq " 2r´x 2 2 ptqu 2 ptq 2p1´x 2 ptqq u 2 ptq´r x 2 2 ptq Kpx 1 ptqqp1´x 2 ptqq pα 2 {γ 2`u2 ptqq, γ 1 u 1 pTq " g 1 1 px 1 pTqq`α 1 , One important observation from the optimality system (18) is that the scaling factors γ 1 for the pollution cost and α 2 for the fishery utility should be non-zero. From (18), γ 1 " 0 is an asymptotic situation with an explosive aquatic pollution characterized by u 1 ptq "`8. Small values of γ 1 indicate small costs to pay for the polluters: the farmers obtain the fertilizers almost free of charge, and industrial companies are not constrained by any regulation policy for their residual chemicals. A typical case promoting this scenario is when political decisions are biased by economical objectives that apply important subventions to obtain fertilizer for farmers or to boost the industrial production for companies, without supporting and controlling the implementation. The scaling factor α 2 " 0 describes a zero-utility for fishery activities and can happen when the demand is very low compared to the harvested fish quantity. Note that the scenario α 2 " 0 differs from the rarity of fishes, although both situations produce the same zero rewards for the fishermen.

Numerical Results
For simulation purposes, we consider a period of 15 months, during which, we suppose, that the fish population, with the carrying capacity K affected by the pollution level x 1 as Kpx 1 q " 1{px 1`1 q, grows according to (18), i.e., the seasonality is neglected. We look at the system of forward-backward equations at (1) and (18) as a boundary value problem and we deploy a fourth-order Runge-Kutta method at each step of the convergence in the shooting method to approximate the solution x 1 , x 2 , u 1 and u 2 . The selected values of the parameters are b " 0.1, the desintegration rate is c " 0.04, the natural growth rate of the fish population is r " 0.25, the coefficients of the utility functions are α 1 " 0.01, α 2 " 2, and the coefficients of the cost functions are γ 1 " γ 2 " 1.0.
The control u 2 is highly sensitive to the initial value of the state variables compared to the control u 1 , see Figure 2. When the pollution starts at a high level, the fish population decreases during the first year, before starting the monotonic growth. However, for lowlevel pollution, the fish population initially grows before becoming damped by the players' activities.
By combining one fixed initial condition state with several terminal conditions for the control variables (see Figure 3), we observe that the pollution control variable u 1 is a consistently increasing function. Moreover, we see that the pollution state grows with the terminal condition value of u 1 for the first two years but decays after. Subsequently, the fish population takes advantage of this and starts a slight growth.
The non-cooperative outcome, presented in the current Section 4, is a dynamic Nash equilibrium between the fishing industry and the other industries involved in pollution. From the perspective of global performance, this non-cooperative differs from the fully cooperative. The idea here can help to cooperate and enforce cooperation between the parties.

Cooperative Strategies and Biodiversity
Now, we examine the cooperative optimization of the players' profits. Instead of a direct bilateral collaboration between the players, the cooperation is coordinated by a third entity such as state representations or international organizations. In practice, countries and international regulations have designed different rules, including incentive mechanisms, taxation, control, risk, and surveillance. Below, we examine one particular disincentive mechanism for pollution control.

The Optimal Strategy
To preserve the fish population, the political entity defines the global payoff and, as strategy designer, is in charge of setting the values of the parameters in the players' utility and the cost functions. The global payoff function is given by: The utilities of the players are similar to the non-cooperative case ones but the cost function changes particularly for the fishermen. Although the primary objective of the cooperation is biodiversity preservation, the fishermen take substantial benefits from it. This is because the efficient running of the fisheries' activities is intrinsically associated with the level of fish population. Later on, the profit for the polluters will be discussed. We are then concerned with the optimal control problem: The pseudo-Hamiltonian is: Hpt, x, pq 1 ,q 2 q,ûq " pq 1 ,q 2 q¨f px,ûq`α 1û1 ptq´γ 1û 2 1 ptq 2`α 2 x 2 ptqû 2 ptq´γ 2 p1´x 2 ptqqû 2 2 ptq 2 .
A maximizer of H with respect to the controlû " pû 1 ,û 2 q provides an open-loop optimal controlû˚associated with the adjoint state pq 1 ,q 2 q given by: 9 q 1 ptq "´B H . We find that: running to the terminal condition:`q 1 pTq,q 2 pTq˘"´B x 1 gpxpTqq, B x 2 gpxpTqq¯. The pseudo-Hamiltonian H is concave in the control variableû; the first-order optimality condition is then sufficient as in the non-cooperative scenario. The optimal control u " pû 1 ,û 2 q of the cooperative strategy is given by BH Bû 1 " 0 and BH Bû 2 " 0, which are written in detail as follows: q 1 ptq`α 1´γ1û1 ptq " 0 andγ 2 p1´x 2 ptqqû 2 ptq "´q 2 ptqx 2 2 ptq`α 2 x 2 ptq.
In Figure 4, we illustrate the optimal dynamics of the state and control variables for a period of 10 months. Here, the shooting method convergence is hard to achieve; therefore, an optimization algorithm is deployed in advance in order to approximate the initial conditions u i p0q associated with the terminal conditions u i pTq. The main observation is that the dynamics of the state and control variables behave quite similarly to the noncooperative regime. In the next section, we investigate the benefits of the cooperative regime.

Taxation: Incentive Design
In this section, we investigate sufficient beneficial conditions for the implementation of the cooperative strategy. In order to achieve efficiency of the cooperative strategy, the decision-makers set, implement, and control the policy for the players, such that their rewards when cooperating are greater than the total sum of the value functions in the non-cooperative scenario. From the dynamics programming approach of Bellman in control theory, the value functions v 1 and v 2 obey the followings Hamilton-Jacobi-Bellman equations: Proposition 12. Given the same circumstances in the non-cooperative and cooperative scenarios (same business growth rate, same level of pollution, and same approved fish stock to be harvested), the sufficient conditions for the players to cooperate, without trying to cheat, are: Proof. Players should accept and align with the cooperative regime if, and only if, the aggregated rewards are greater than the sum of what they earn, as business outcomes, in the non-cooperative case, i.e., v 1`v2 ď w. The business growth rate is given by the time variation of the value function: in the non-cooperative scenario, it is given by Bt , while in the cooperative regime, we have the global growth rate Bw Bt . The same business growth rate in both regimes means that Bw Bt " Bt . Then, it remains to compare H 1`H2 and H with the same level of pollution and same fish stock to be harvested.
The utility function of the fishermen α 2 x 2 u 2 in the non-cooperative regime andα 2 x 2û2 in the cooperation scenario should remain the same to avoid rough variations of the fish price that amplify the life cost for the population. Similarly, the gross gain of the polluters is basically given by the level of the residual chemicals and might depend also on the market state, but not on whether there is cooperation or not. This yields the condition α 1 "α 1 and the analysis is reduced to compare: However, the polluters expect better business growth with regard to the amount of generated residual chemicals when cooperating, i.e., . Therefore, we find that the sufficient condition v 1`v2 ď w can be formulated as follows: It is obvious that (25) holds for´γ 1 ă´γ 1 and´γ 2 ă´γ 2 , which means that γ 1 ąγ 1 and γ 2 ąγ 2 .
The plotted running costs in Figure 5 are obtained by taking the arithmetic mean over 10 runs for each case. Our comparison analysis matches well with the theoretical statement, which claims that (24) is a sufficient but not a necessary condition to align all the players with the cooperative policy. In fact, we can see clearly that the plot γ i ąγ i for i " 1, 2 gives a consistent advantage to the cooperation. Comparison of the aggregated running cost in the cooperative regime (blue line) and the sum of the two running costs in the non-cooperative regime (red dashed line) starting at (x 0 1 " 0.5, x 0 2 " 0.2). Pollution control and fish population control after T " 10 months are u 1 pTq " 0.7 and u 2 pTq " 0.9, respectively, and α 2 " 2.01 andα 2 " 1.0.

Conclusions
We introduced a system of differential equations describing the impact of the pollution from industrial and agricultural companies in the fish population in watercourses subject to fishery activities. The stability analysis highlights the necessity of controlling the pollution input to preserve the fish population and, subsequently, the fishery industries. The noncooperative optimal control is not the worse-case scenario for fishermen as both classes of players incorporate the fish preservation into the payoff definition. However, the cooperative strategy under the authority of a third entity seems to be more advantageous for fish preservation. Implementing an incentive taxation approach for setting the parameters' values helps the third entity as it means that the players should not try to cheat.
One central aspect missed in this work is as follows: normally, it is not the pollution caused by an individual farmer or one single industry that harms people's health but it is the aggregate effect of all polluting entities. This is the starting point of our future work to introduce a mean-field approach [20,21] with the distribution of polluters and also to quantify the pollution generated by fishermen.

Conflicts of Interest:
The authors declares no conflict of interest.