# A Differential Game of Ecological Compensation Criterion for Transboundary Pollution Abatement under Learning by Doing

Authors: Zhigang Chen; Rongwei Xu; Yongxi Yi
Journal: Discrete Dynamics in Nature and Society
 (2020)

---

## Abstract

This paper studies a stochastic differential game of transboundary pollution abatement between two kinds of ecological compensation and the abatement policy, in which the learning by doing is taken into account. Emission and pollution abatement between upstream and downstream region in the same basin is a Stackelberg game, and the downstream regions provide economic compensation for pollution abatement in the upstream region. We discuss the feedback Nash equilibrium strategies of proportional compensation and investment compensation, and it is found that an appropriate ecological compensation ratio can improve the investment level of pollution abatement in the two regions by accumulating experience in the process of learning by doing. In the long term, the investment compensation mechanism is an effective transboundary pollution abatement measure that can continuously reduce the water pollution stock in the upstream and downstream.


---

## Body

## 1. Introduction

Environmental degradation and its relationship with economic activities have been a subject of intense debate in the academic community and among policy makers. One of the most important environmental problems is water pollution. Pollution in one region can be transported across hundreds and even thousands of kilometers according to Li [1]. Hence, we call it “transboundary pollution”. The natural and transregional mobility of the water flow make the regional division and management of water resources very complicated, which leads to the increasingly serious and deteriorating problem of transboundary water pollution in the river basin. Dam construction and random discharge of harmful pollutants happen frequently, and regional contradictions are becoming more and more serious, even triggering tremendous panic in the whole society. Environmental protection of river basin transboundary water pollution has become the focus and difficulty of local government departments.Differential games provide an effective tool to study pollution control problems and to analyze the interactions between the participants’ strategic behaviors and dynamic evolution of pollution. Some researchers have paid their attention to the transboundary industrial pollution problems in recent years. For example, Yeung [2] first derived the time consistent solutions in a cooperative differential game and studied the pollution management in a stochastic differential game framework. In paper [3], a cooperative stochastic differential game of transboundary industrial pollution is presented, and a payment distribution mechanism is derived to maintain the subgame consistency. Additionally, there are several published studies of transboundary pollution problems from other views, such as renewable resource, clean technologies, harmonization of international and domestic law, abatement cost, and R&D spillovers (for instance, [4–7]). Papers [1, 3, 8], taking emission permits trading into account, study the problem of transboundary pollution. It has been shown in this literature that coordination of countries’ emission strategies leads to a lower total level of pollution and to higher total welfare than when countries use noncooperative emission strategies.This diffusion of pollutants is common. In each region, those who suffer from pollution wish that the polluter in neighboring regions would either reduce polluting or compensate for the damage, and there is generally a game between the sufferers and the polluter around pollution abatement and benefit compensation [9]. Ecological compensation is an important mechanism in environmental economics and can be used to promote and coordinate harmonious development between humans and nature, and it is an effective way to reduce water environment pollution and limit the discharge of pollutants. In recent years, in watershed studies of ecological compensation, researchers have calculated standard values of ecological compensation that account for water, land, waterways, and society and have specifically considered aspects such as water volume, water ecology, and regional economy [10]. Evaluation of the ecological service function [11] and the willingness-to-pay approach [12, 13] have been used in most of these studies. A river basin ecological compensation model reflects multiple factors. A compensation model that only includes a single aspect will fail to meet the fundamental requirements of river basin ecological compensation. This means that the traditional method of focusing solely on water quality should be abandoned, and objective factors that reflect regional heterogeneity should be included in the methods for ecological compensation standards in river basins [14, 15]. Eco-compensation criterion is considered as one of the appropriate instruments for internalizing the spillover of public goods, which encourages the establishment of a longitudinal compensation relation between polluted downstream and upstream regions through capital supply, industrial transfer, water rights, and carbon trading [16–18]. Heitzig et al. [19] believe that, through reasonable ecological and economic compensation methods, pollution abatement losses in some areas of the basin’s water environment can be compensated, and failed environmental compensation cooperation can be further corrected.Articles using differential game methods to study ecological compensation mechanisms are relatively scarce [20, 21], especially combining ecological compensation standards accumulated by transboundary pollution abatement technologies. The interaction between ecological, social, and economic systems, as well as pollution abatement between administrative regions, from upstream to downstream, has been considered to ensure the foundation of ecological compensation for the water environment in the basin. This work studies the transboundary water pollution control problem between two regions, characterizing and contrasting feedback Nash strategies under two kinds of ecological compensation mechanisms. The main purpose is to characterize the parameter spaces in which governments or organizations in two neighboring regions can formulate a consistent ecological compensation mechanism and abatement policy under the influence of regional impacts. Besides, we attempt to explore a differential game model of the interregional ecological compensation mechanism covering the upstream and downstream regions in a continuous time. At the same time, the accumulation of pollution investment technology—learning by doing—is taken into account in this work. The experience will be accumulated continuously in the process of learning by doing, which means that the abatement cost will decrease with the improvement of abatement technology.The rest of this paper is organized as follows. In Section2, we will establish our basic dynamic general equilibrium model between two kinds of ecological compensation and pollution abatement. In addition, the optimal abatement levels and optimal pollution stock paths under proportional compensation and investment compensation are presented in Section 3 and Section 4, respectively. Some discussions are provided with several numerical examples in Section 5. Finally, Section 6 concludes the paper.


## 2. The Basic Model

There are two adjacent regionsn=1,2 in a river basin, which we call upstream region 1 and downstream region 2, in our transboundary pollution model. It is assumed that both regions discharge organic pollutants into the river basin. For region nn=1,2, production always leads to a quantity of by-products, namely, emissions Entn=1,2. We assume that Qnt is the industrial output of upstream and downstream region, which indicates that, at time t, the industrial production of region 1 and region 2 is Q1t and Q2t. The instantaneous emissions of pollutants from industrial output are Entn=1,2; that is, at time t, the pollution emissions of region 1 and region 2 are E1t and E2t, respectively. It is assumed that pollution emissions in various regions are positively related to industrial production, and the instantaneous linear production function can be expressed as(1)Qnt=QnEnt.According to literature [22–24], it is assumed that the regional industrial income function is RnQnt, which can be expressed by the following quadratic concave functional form in terms of emissions:(2)RnEnt=EntAn−12Ent,0≤Ent≤2An.Environmental pollution damage cost is a linear function of pollution stockPnt, following [23], the cost function is Dnt=DnPntDn≥0. Among them, Dn indicates the degree of environmental damage per unit pollution stock to region n.The difference between the cooperative solution and the noncooperative solution depends on the spillover effects of benefits and costs, which are likely asymmetric due to geographic and economic differences between the regions. The upstream region imposes a negative unidirectional water burden upon the downstream region by preventing the latter from reaching an unconstrained optimal water quality. Compensation through side payments is a viable way to create incentive for cooperation. Even where resource flows occur unidirectionally, precluding mutual control over negative externalities, there is the possibility of Pareto optimal or “win-win” solutions through side payments.According to the payment principle of the ecological compensation mechanism, the downstream polluted regions provide certain economic compensation for the investment of the upstream pollution control. This compensation is used to promote pollution control in the upstream regions, thereby reducing the transfer of pollutants to the downstream basins and improving the quality of the water environment. It is assumed that there are two types of compensation: (i) Proportional compensation means to compensate for a certain proportion of investment in pollution abatement. We assume that the proportion of ecological compensation in the downstream region to the upstream region isBt, and the proportion is determined by the downstream region. (ii) Investment compensation means that the downstream region invests in the upstream region to reduce pollution so as to ensure that fewer pollutants are transferred to the downstream region. The investment intensity function of pollution abatement is Intn=1,2. It is known that pollution abatement can be realized only when technique and labor are invested. Therefore, we should face the abatement cost which could decrease the net revenue. Following [22], we assume that the abatement cost can be described by the following quadratic form:(3)CnIn=12μnIn2t,μn≥0.where μnn=1,2 represents positive constants. It measures the difference between the two regions’ ability in mastering the abatement technology. Equation (3) means that the marginal cost is increasing with respect to the level of pollution abatement. By means of [25, 26], the experience of applying pollution abatement technology Gnt is measured by the cumulative abatement from time 0 to t; that is,(4)Gnt=Gn0+gn∫00Insds,Gn0=Gn0.We defineGn·t>0, where Gn0n=1,2 denotes the initial experience level of applying pollution abatement technology. Similar to the above, gnn=1,2 is a positive parameter and it represents the differences between the two regions’ ability in accumulating experience. According to the learning-by-doing theory, the amount of cumulative experience will lead to a decline in the unit cost.As we all know, with the change of time and temperature, the water in nature has a certain self-purification ability. Without the loss of generality, it is assumed that the basin has the purification rate of water pollution, namely, the coefficient of self-purification ability of waterθnn=1,2, θ1,θ2≥0. The pollution in the upstream of the basin mainly comes from the discharge of local polluting industries, while the pollution in the downstream comes from the transfer of pollution from the upstream of the basin, besides the discharge of polluting enterprises in the downstream region. The amount of pollution transferred from upstream to downstream is dP1t, where d is assumed to be the transfer coefficient and 0≤d≤1.The number of initial emission permits in region 1 isE10 and in region 2 is E20. It is assumed that the emission trading market is a fully competitive market and the price of emission trading right is constant at w. If the amount of emission exceeds the initial allocation, the emission right can be purchased in the permit market. On the contrary, if there is a surplus of the emission right, it can be reserved for the next year or sold in the emission permit trading market.Given the above assumptions, we can get the concrete expression of the income functionWnt of the two regions:(5)W1t=E1tA1−12E1t−wE1t−I1t−E10−12μ1I12t+Gnt−Gn0,W2t=E2tA2−12E2t−wE2t−E10−Gnt−Gn0μ1BtI12t.


## 3. Proportional Compensation

The amount of compensation provided by the downstream region in the same basin will affect the investment enthusiasm of the upstream region in pollution control and the amount of pollutant discharge in the region. In turn, the amount of pollutant discharge in the upstream region will affect the amount of pollutant discharge in the downstream region, which constitutes a dynamic game relationship between the two sides. Considering the decision-making in continuous time, this constitutes a dynamic differential game relationship.Under this model, both regions aim to maximize the net present value of their own long-term income. The pollution discharge in the upstream region will affect the income of the downstream region by affecting the pollution stock in the same water basin. The decision-making problem of independent discharge from the two regions constitutes a differential game problem withE1t, E2t, I1t as control variables and P1t, P2t as state variables, aiming at maximizing the net present value of their respective income.The stocks of pollution in the upstream and downstream regions of the basin areP1t and P2t at time t. Also, the pollution stock in the two basins can be expressed by the following two differential equations:(6)P1t•=E1t−I1t−θ1+dP1t,P10=P10,P1t≥0,P2t•=E2t−θ2P2t+dP1t,P20=P20,P2t≥0,where θ1P1t and θ2P2t represent the stocks of self-purification pollution. The investment intensity of pollution abatement in upstream region is I1t. dP1t is the amount of pollution transferred from upstream to downstream.The current goal of region 1 is to maximize the expected present flow of instantaneous net revenue in terms of the emission path and the abatement level. Hence, the objective functional and the constraint conditions of region 1 can be given as follows:(7)J1Bt=maxE1≥0,I1≥0∫0∞e−rtA1E1t−12E12t−wE1t−I1t−E10−12μ1−BtI12t+g1∫0tI1sds−D1P1tdt,(8)s.t.P1t•=E1t−I1t−θ1+dP1t,P10=P10,P1t≥0.As above, the current goal of downstream region 2 is to maximize the expected present flow of instantaneous net revenue in terms of the emission path and the abatement level. Hence, the objective functional and the constraint conditions of region 2 can be given as follows:(9)J2Bt=maxE2≥0,I1>0∫0∞e−rtA2E2t−12E22t−wE2t−E20−12μBtI12t−D2P2tdt,(10)s.t.P1t•=E1t−I1t−θ1+dP1t,P10=P10,P1t≥0,P2t•=E2t−θ2P2t+dP1t,P20=P20,P2t≥0.In a transboundary pollution control problem, the neighboring regions can be seen as players and they aim at maximizing the joint or their own net present profits under noncooperative games. Under the model of proportional compensation mechanism, the problem is a Stackelberg differential game that determines the amount of pollution emission and the compensation ratio from the downstream of the basin and determines the amount of pollution emission and abatement in the upstream region. Each player tries his best to maximize his own net revenue by choosing the optimal emission path and the optimal abatement level under a noncooperative game. In order to obtain the optimality conditions for the optimal control problems, we use Pontryagin’s maximum principle.The current-value Hamiltonian function of (7) is(11)H1B=A1E1t−12E12t−wE1t−I1t−E10−12μ1−BtI12t+g1∫0tI1sds−D1P1t+λ1E1t−I1t−θ1+dP1t.The current-value Hamiltonian function of (9) is(12)H2B=A2E2t−12E22t−wE2t−E20−12μBtI12t−D2P2t+λ2E2t−θ2P2t+dP1t+λ3E1t−I1t−θ1+dP1t,where λnn=1,2,3 are the dynamic adjoint variables associated with the state equation about Pnt. Here, the dual variables λ1, λ2, and λ3, also called shadow prices or common state variables, are Lagrange multipliers, which are the derivatives of the two players’ value functions, i.e., revenues, with respect to the pollution stock Pnt.To maximize (11) and (12), the first-order conditions of the current Hamiltonian function are the following:(13)∂H1B∂E1t=A1−w−E1t+λ1B=0,(14)∂H1B∂I1t=w−1−BtμI1t+g1I1tI1t·−λ1B=0,(15)∂H2B∂E2t=A2−w−E2t+λ2B=0.Along with the current value costate equations:(16)λ1B•=rλ1−∂H1B∂P1t=r+θ1+dλ1+D1,(17)λ2B•=rλ2−∂H2B∂P2t=r+θ2λ2+D2,(18)λ3B•=rλ3−∂H2B∂P1t=r+θ1+dλ3−λ2d.The initial state variable values and cross-sectional conditions are as follows:(19)P10=P10,P20=P20,PnT≥0,λ1T=λ2T=λ3T=0.Solving (13), (14), and (15), we have(20)E1B∗t=A1−w+λ1B,(21)E2B∗t=A2−w+λ2B,(22)I1t·=g1I1tλ1B−w+1−BtμI1t.Thus, we can obtain the following results from (16), (17), and (18):(23)λ1B=−D1r+θ1+d,(24)λ2B=−D2r+θ2,(25)λ3B=−D2dr+θ2r+θ1+d.Proposition 1.
The conditions given above are necessary. Next, we use the sufficiency theorem (O. L. Mangasarian’s theorem) to prove that the necessary conditions are also sufficient.The functionsW⋅ and P˙t are differentiable in this problem and all functions about P,E,I are concave. In addition, the function P˙t is a linear function about P,E,I, so, in the optimal solution, the positive and negative values of λt can be ignored (no sign restriction is required). Therefore, we know from Mangasarian’s theorem that, for the optimal control problem in this paper, the necessary condition of the Pontryagin maximum principle is also a sufficient condition for the global maximization of J⋅.For the proof of the theorem validity, see Appendix.Converting (22) into (12), we get the present-value Hamiltonian function of downstream region 2:(26)H2B=A2E2t−12E22t−wE2t−E20−12Btμw−λ1Bμ1−Bt−g1/I1t·2−D2P2t+λ2E2t−θ2P2t+dP1t+λ3E1t−w−λ1Bμ1−Bt−g1/I1t·−θ1+dP1t.To maximize (26), the first-order conditions of the current Hamiltonian function are the following:(27)0=∂H2B∂Bt=−12μw−λ1Bμ1−Bt−g1/I1t·2−λ3Bμw−λ1Bμ1−Bt−g1/I1t·2−Btμw−λ1Bμ1−Bt−g1/I1t·μw−λ1Bμ1−Bt−g1/I1t·2.Solving (27), we get the optimal investment compensation ratio function of downstream region 2 to upstream region 1 as follows:(28)B∗t=w−λ1B+2λ3Bg1w−λ1B−2λ3BμI1t·−1.Numerically, solving an ordinary differential equation system means computing a sequence of points which are close to the graph of the real solution. The numerical solution algorithm is based on the numerical integration of the system using the Runge–Kutta method. In numerical analysis, the Runge–Kutta methods are an important family of implicit and explicit iterative methods for the approximation of solutions of ordinary differential equations. Our aim in this section is to study the numerical solution of the differential equation (22) to determine an optimal control strategy, and its time evolution is discussed. The parameters of (22) which are used in the numerical solution are presented in Table1 [27–29].However, it is not difficult to find from (22) and (28) that abatement levels I1t in region 1 and eco-compensation ratio Bt are interconnected dynamic systems. In order to obtain a numerical solution of abatement levels, we assume that the eco-compensation ratio Bt is 0.2, 0.4, 0.6, 0.7, 0.8, and 0.9, respectively. The optimal numerical solution of abatement levels under different Bt is obtained using the Runge–Kutta algorithm and is displayed in Figure 1.From Figure1, we find that when the proportion of ecological compensation continues to increase, that is, when downstream region 2 gives more and more economic compensation to upstream region 1, the abatement level in upstream region 1 will continue to increase. When the compensation ratio exceeds a certain value, about 0.6 in our simulation, the abatement level will be drastically reduced. In other words, a suitable ecological compensation ratio will increase the enthusiasm of pollution control investment in the upstream region. However, excessive economic compensation from the downstream regions will also greatly reduce this enthusiasm, thereby increasing the reliance on such economic compensation.Next, we simulated the system equation (6) of the pollution stock in the upstream and downstream and obtained the optimal numerical solutions P1 and P2 of the pollution stock in two regions. The results are shown in Figure 2.It can be seen from Figure2 that, with the passage of time, under the dual effects of pollution abatement investment in region 1 and economic compensation from downstream region 2, the pollution stock in the upstream will gradually decrease as the investment of pollution abatement increases. In the short term, the pollution stock in the downstream region will increase sharply because of the transfer of pollutants from the upstream region. However, in the long run, with the increasing investment in pollution abatement in the upstream region, the pollution stock transfer will decrease relatively, which has a great mitigation effect on the pollution stock in the downstream region.In order to evaluate the impact of parameter change on the prediction results of the model, we performed sensitivity analysis by changing the parametergn; meanwhile, other parameters were fixed to study the effect of the change of system parameters values on objective function value. The analysis results are shown in Figures 3 and 4. From both figures, it can be concluded that the pollution abatement level is the most sensitive to the change of the parameter gn, and pollution stock is the least sensitive to the change of the parameter gn.Table 1
The parameters used in the numerical example.
A1A2wrD1D2E10E208.57.650.20.080.20.2456P10P20θ1θ2dμ1μ2μ31001100.10.110.50.40.50.3Figure 1
Abatement levels under differentBt.Figure 2
Pollution stock of the upstreamP1 and downstream P2.Figure 3
Sensitivity analysis of investment in pollution abatement.Figure 4
Sensitivity analysis of pollution stock.


## 4. Investment Compensation

In the investment compensation model, the downstream regions in the basin invest in pollution abatement in the upstream to reduce the degree of water pollution in the upstream region and the amount of pollution transfer downstream. In the long run, the cost of investing in pollution control in the downstream region is relatively lower, which is conducive to the long-term benefits of water environment improvement in the downstream region. The compensation model of investment in different places can better exert the effect of interregional cooperation and integrate various forms of capital to solve the dilemma of pollution abatement in river basins.Different from the proportional compensation model, the assumption of the pollution abatement cost learning-by-doing model, the pollution stock function in the upstream and downstream will change accordingly.The investment cost of water pollution abatement in downstream region 2 and the investment cost of water pollution control in upstream region 1 can be expressed as follows:(29)C22I22=12μ2I222t,μ2≥0,(30)C21I21=12μ3I22t+I21t2−I222t=μ3I22tI21t+12μ3I212tμ3≥0,where I22t indicates the investment intensity of downstream region 2 in local area, and I21t indicates the investment intensity of downstream in upstream region 1. Correspondingly, the learning-by-doing function has also changed due to different investment intensity:(31)Gnnt=Gnn0+gnn∫0tInnsds,Gnn0=Gnn0.In this model, the stocks of pollution in the upstream and downstream regions of the basin areP1t and P2t at time t, and the pollution stock in the two basins can be expressed by the following two differential equations:(32)P1t•=E1t−I11t−I21t−θ1+dP1t,P10=P10,P1t≥0,P2t•=E2t−I22t−θ2P2t+dP1t,P20=P20,P2t≥0,where θ1P1t and θ2P2t represent the stocks of self-purification pollution. The investment intensity of pollution abatement in upstream region is I11t. dP1t is the amount of pollution transferred from upstream to downstream.The current goal of region 1 is to maximize the expected present flow of instantaneous net revenue in terms of the emission path and the abatement level. Hence, the objective functional and the constraint conditions of region 1 can be given as follows:(33)J1St=maxE1,I21,I11≥0∫0∞e−rtA1E1t−12E12t−wE1t−I11t−I21t−E10−12μ1I112t+g11∫0tI11sds−D1P1tdt,(34)s.t.P1t•=E1t−I11t−I21t−θ1+dP1t,P10=P10,P1t≥0.As above, the current goal of downstream region 2 is to maximize the expected present flow of instantaneous net revenue in terms of the emission path and the abatement level. Hence, the objective functional and the constraint conditions of region 2 can be given as follows:(35)J2St=maxE2,I21,I22>0∫o∞e−rtA2E2t−12E22t−wE2t−I22t−E20−12μ2I222t−12μ3I212t−μ3I21tI22t+g21∫0tI21sds+g22∫0tI22sds−D2P2tdt,(36)s.t.P1t•=E1t−I11t−I21t−θ1+dP1t,P10=P10,P1t≥0,P2t•=E2t−I22t−θ2P2t+dP1t,P20=P20,P2t≥0.The current-value Hamiltonian function of (33) is(37)H1S=A1E1t−12E12t−wE1t−I11t−I21t−E10−12μ1I112t+g11∫0tI11sds−D1P1t+λ1Et−I11t−I21t−θ1+dP1t.The current-value Hamiltonian function of (35) is(38)H2S=A2E2t−12E22t−wE2t−I22t−E20−12μ2I222t−12μ3I212t−μ3I21tI22t+g21∫0tI21sds+g22∫0tI22sds−D2P2t+λ2E2t−I22t−θ2P2t+dP1t+λ3E1t−I11t−I21t−θ1+dP1t,where λnn=1,2,3 are the dynamic adjoint variables associated with the state equation about Pnt. Here, the dual variables λ1, λ2, and λ3, also called shadow prices or common state variables, are Lagrange multipliers, which are the derivatives of the two players’ value functions, i.e., revenues, with respect to the pollution stock Pnt.To maximize (37) and (38), the first-order conditions of the current Hamiltonian function are the following:(39)∂H1S∂E1t=A1−w−E1t+λ1S=0,(40)∂H1S∂I11t=w−μ1I11t+g11I11tI11t·−λ1S=0,(41)∂H2S∂E2t=A2−w−E2t+λ2S=0,(42)∂H2S∂I22t=w−μ2I22t−μ3I21t+g22I22tI22t·−λ2S=0,(43)∂H2S∂I21t=−μ3I22t−μ3I21t+g21I21tI21t·−λ3S=0.Along with the current value costate equations:(44)λ1S•=rλ1S−∂H1S∂P1t=r+θ1+dλ1S+D1,(45)λ2S•=rλ2S−∂H2S∂P2t=r+θ2λ2S+D2,(46)λ3S•=rλ3S−∂H2S∂P1t=r+θ1+dλ3S−dλ2S.Solving (39)–(43), we have(47)E1S∗t=A1−w+λ1S,(48)E2S∗t=A2−w+λ2S,(49)I11t·=g11I11tλ1S−w+μ1I11t,(50)I22t·=g22I22tλ2S−w+μ2I22t+μ3I21t,(51)I21t·=g21I21tλ3S+μ3I21t+μ3I22t.Proposition 2.
The conditions given above are necessary. Next, we use the sufficiency theorem (O. L. Mangasarian’s theorem) to prove that the necessary conditions are also sufficient.The functionW⋅ and P˙t are differentiable in this problem and all functions about P,E,I are concave. In addition, the function P˙t is a linear function about P,E,I, so, in the optimal solution, the positive and negative values of λt can be ignored (no sign restriction is required). Therefore, we know from Mangasarian’s theorem that, for the optimal control problem in this paper, the necessary condition of the Pontryagin maximum principle is also a sufficient condition for the global maximization of J⋅.For the proof of the theorem validity, see the appendix.From (49), (50), and (51), we can find that abatement levels I11t, I21t, and I22t are interconnected dynamic systems. In order to obtain a numerical solution of the abatement levels, we use the Runge–Kutta algorithm to solve them. The numerical solutions are plotted in Figure 5, in which the parameters are from Table 1.In the initial stage, the investment level of pollution abatement in the downstream region is higher than that in the upstream region, so as to achieve the overall goal of water pollution control. Due to the long-term effect of investment compensation mechanism, with the increase of investment compensation level in the downstream region, the quality of water environment will be improved in a certain time. Because of this, the investment pressure of downstream pollution control will be reduced, and the investment level of pollution abatement will tend to grow slowly. This can be clearly seen from theI22t curve in Figure 5.Next, we simulated the system equation (32) of the pollution stock in the upstream and downstream and obtained the optimal numerical solutions P1 and P2 of the pollution stock in two regions. The results are shown in Figure 6.From Figure6, we can see that, under the dual effects of pollution abatement investment in region 1 and economic compensation from downstream region 2, the pollution stock in the upstream will gradually decrease as the investment of pollution abatement increases. In the short term, the pollution stock in the downstream region will increase sharply because of the transfer of pollutants from the upstream region. However, in the long run, with the increasing investment in pollution abatement in the upstream region, the pollution stock transfer will decrease relatively, which has a great mitigation effect on the pollution stock in the downstream region.In order to evaluate the impact of parameter change on the prediction results of the model, we performed sensitivity analysis by changing the parametergn; meanwhile, other parameters were fixed to study the effect of the change of system parameters values on objective function value. The analysis results are shown in Figure 7.Figure 5
Abatement levels of upstream and downstream.Figure 6
Pollution stock of the upstreamP1 and downstream P2.Figure 7
Sensitivity analysis of pollution abatement levels and pollution stock.
(a)(b)(c)(d)


## 5. Discussion

The biggest difference between the proportional compensation mechanism and the investment compensation mechanism is that the proportional compensation directly compensates the pollution abatement behavior in the upstream region in the form of funds. Therefore, this kind of compensation has the characteristics of timely compensation, high capital availability rate, and rapid effect on pollution abatement enthusiasm and can produce better environmental treatment effects in a short period of time. This can be seen from the path of the pollution stock change in Figure8, where the dotted line indicates the trace of the pollution stock under the proportional compensation mechanism.Figure 8
Comparison between two kinds of compensation in pollution stock.However, the investment compensation mechanism needs a higher investment cost of environmental governance and longer time to improve the mechanism. In the short term, the effect of investment compensation mechanism on water pollution control in river basins is relatively poor, but, in the long term, the investment in pollution abatement in the upstream region can play a strong sustained role and can continuously reduce the pollution level in the waters for a longer period of time. The full line track in Figure8 reflects this process. Therefore, the investment compensation mechanism is a sustainable investment mechanism for water transboundary pollution abatement.


## 6. Conclusion

This paper discusses the dynamic optimal strategies in the transboundary pollution game under two kinds of ecological compensation mechanism and the abatement policy, in which the learning by doing is taken into account. The experience will be accumulated continuously in the process of learning by doing, which means that the abatement cost will decrease with the improvement of abatement technology. By solving dynamic equations and numerical simulation, we find that all the results show that the two regions will promote the abatement levels and reduce the pollution stocks when the learning by doing is considered. At the same time, we also found the following:(i)
With the continuous improvement of the proportion of ecological compensation, the abatement level in upstream region 1 will continue to increase. When the compensation ratio exceeds a certain value, about 0.6 in our simulation, the abatement level will be drastically reduced. However, under the investment compensation mechanism, in the initial stage, the investment level of pollution abatement in the downstream region is higher than that in the upstream region. Due to the long-term effect of investment compensation mechanism, with the increase of investment compensation level in the downstream region, the quality of water environment will be improved in a certain time.(ii)
On the other hand, whether in proportional compensation or in investment compensation, the pollution stock in the upstream will gradually decrease as the investment of pollution abatement increases. In the short term, the pollution stock in the downstream region will increase sharply because of the transfer of pollutants from the upstream region. However, in the long run, the pollution stock transfer will decrease relatively, which has a great mitigation effect on the pollution stock in the downstream region. However, the decreasing process and speed of pollution stocks are different in both cases.(iii)
Generally speaking, a suitable ecological compensation ratio will increase the enthusiasm of pollution control investment in the upstream region. However, excessive economic compensation from the downstream regions will also greatly reduce this enthusiasm, thereby increasing the reliance on such economic compensation.(iv)
In the long term, the investment in pollution abatement in the upstream region can play a strong sustained role and can continuously reduce the pollution level in the waters for a longer period of time. Therefore, for the environmental regulation department, the investment compensation mechanism is a sustainable one for water transboundary pollution abatement.



---

*Source: 7932845-2020-07-27.xml*