Next Article in Journal
The Relationship between Scabies and Stroke: A Population-Based Nationwide Study
Previous Article in Journal
Perspectives on Disability and Non-Communicable Diseases in Low- and Middle-Income Countries, with a Focus on Stroke and Dementia
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optimal Coordination Strategy for International Production Planning and Pollution Abating under Cap-and-Trade Regulations

1
Nonlinear Science Center, College of Economics and Management, Shandong University of Science and Technology, Qingdao 266590, China
2
Department of Management Science and Statistics, College of Business, The University of Texas at San Antonio, San Antonio, TX 78249-0634, USA
*
Authors to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2019, 16(18), 3490; https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph16183490
Submission received: 27 August 2019 / Revised: 15 September 2019 / Accepted: 16 September 2019 / Published: 19 September 2019
(This article belongs to the Section Environmental Science and Engineering)

Abstract

:
Because both pollution emissions and production policies often are international in scope, it is necessary to find optimal coordination strategies for international production planning and pollution abating. Differential game models are developed for multiple neighboring countries to reach optimal decisions on their production planning and pollution abating under cap-and-trade regulations. Non-cooperative and cooperative differential games are presented to depict the optimal tradeoffs between production planning and pollution abating. Hamilton-Jacobi-Bellman (HJB) equations are then employed to analyze the asymmetric and symmetric feedback solutions. Numerical simulations are used to illustrate the results. Five different dividends are also discussed. With the proposed strategies, more improvement will be directed toward production supplies and environmental issues than ever before.

1. Introduction

Since the middle of the 20th century, the rapid population growth and economic development have led to serious environmental pollution that has critical impacts on crop yields, climate change, biodiversity, and ecosystems, amongst others. As a result, these impacts have significant social and economic consequences affecting values, customs, human health, economic development as well as social welfare [1]. In recent years, the growing concerns about industrial pollution, the strict regulations of international conventions on sustainability and environmental protection, and an increase in the number of pro-environmental consumers have led industries to dedicate significant efforts to developing green practices [2,3,4,5]. However, as of 2017, almost all advanced industrialized countries have had falling their emission rates, but still have not met their pledged emission reduction targets that they made in Paris [6]. Therefore, it is entirely necessary to conduct studies in pollution abatement.
The Paris climate agreement, a new global agreement to combat climate change within the United Nations Framework Convention on Climate Change (UNFCCC), was signed in 2016. The agreement aims at holding global warming to well below 2 °C and to “pursue efforts” to limit it to 1.5 °C [7]. As of March 2019, 195 countries have signed the agreement, and 186 of which have submitted emission abatement pledges as nationally determined contributions. To meet the pledges, some instruments are needed to regulate and ultimately reduce the amount of pollution emitted into the atmosphere. One such instrument is the cap-and-trade mechanism. In recent years, several countries, including China, have implemented cap-and-trade programs to set a national limit on carbon emissions by the heavy-polluting industries such as power generation, iron and steel, and chemicals, amongst others [8]. Firms must have enough permits to cover their emissions if they exceed the cap. Generally, firms can invest in clean technologies to become more efficient, switch to lower carbon fuels or purchase additional permits. Indeed, they may sell some permits if they do not use up the cap. By comparison with pollution tax, the cap-and-trade mechanism creates a commodity out of the right to emit pollutants and allows the commodity to be traded on the free market. However, there exist several dilemmas, such as a firm’s higher income depends on more output which may lead to more pollution emissions, and fewer permit purchases or more permit sales depend on greater abatement efforts which may lead to profit losses [9]. Thus this study deals with real demand for production planning and pollution abating under cap-and-trade regulations.
Pollution emissions often show the feature of being international, that is, there exists transboundary pollution which is a typical example of a negative externality and is also an essential instance of global environmental agreement failure [10]. When a negative externality occurs in neighboring countries, a country does not take responsibility for external costs but passes them on to others. Thus the producing country has lower marginal costs than it would otherwise have and may choose to produce more of the product than it would when it was required to pay all associated pollution costs. Hence, some governance mechanisms could be designed to “internalize” an externality. One such mechanism is to employ a differential game to optimize their production planning and pollution abating under cap-and-trade regulations.
Given the concerns mentioned above, the key contribution of this study is on developing differential game models for multiple neighboring countries to reach optimal decisions by optimizing their production planning and pollution abating under cap-and-trade regulations. By using Hamilton-Jacobi-Bellman (HJB) equations, the impacts of the production efficiency and the initial allocations of emission permits on the optimal production plans and the pollution abating decisions can be determined by analyzing the feedback equilibria of the proposed models.
The remainder of this paper is organized as follows. The relevant literature is reviewed in Section 2. In Section 3, differential game models are proposed for multiple neighboring countries. Utilizing the HJB equations, the non-cooperative differential games are studied in Section 4, and the cooperative differential games are discussed in Section 5. In Section 6, numerical simulations are performed to verify the results. In Section 7, the cooperative and symmetric dividends are defined, and the evolutions are analyzed by using numerical methods. In Section 8, the potential policies and managerial implications are discussed. Finally, a summary and future research directions are presented in Section 9.

2. Literature Review

As a powerful analytical tool, the differential game has been widely used in many scientific fields. You, Jiang and Li [11] employed the differential game to design optimal coordination strategies of regional emission abatement collaboration. Feng and Liu [12] used the differential game to consider a single supply chain structure concerning the online word-of-mouth effect. Schüller, Staňková and Thuijsman [13] put forward a differential game of pollution control while there are transboundary environmental effects of national policies. Di Liddo [14] set up a differential game to coordinate a pharmaceutical company controlling the drug price and a social planner determining the number of people to treat.
Many researchers have made significant efforts in searching for optimal strategies of pollution control, as shown in Table 1. Ploeg and Zeeuw [15] proposed a model of international pollution control. Yeung [16] employed a cooperative differential game to set up a model of transboundary industrial pollution control. He and Hua [17] studied transboundary pollution control strategies about production quantity, pollution taxes and abatement investment. Bertinelli, Camacho and Zou [18] proposed strategies to motivate business efforts to reduce CO2 emissions through a capture and storage mechanism. This study employs differential games to balance optimally business production planning and pollution abating to optimize income.
Transboundary/international pollution emission and abatement are now becoming one of the hottest topics in scientific research. A differential game is one of the powerful tools for studying transboundary pollution emission and abatement. Long [19] employed a differential game to depict the transnational pollution between two sovereigns, and analyzed its symmetric open-loop Nash equilibrium and Stackelberg leadership. Breton, Zaccour and Zaha [20] proposed a finite-horizon differential game model to analyze the joint implementation of environmental projects. List and Mason [21] developed an asymmetric differential game and determined the cooperative and non-cooperative equilibria. Yu and Xin [22] proposed a stochastic differential game to describe greenhouse gas (GHG) emission decision making of developed and developing countries. Yeung [16] presented a cooperative differential game of transboundary industrial pollution with industries and governments being separate entities and designed a payment distribution mechanism. Jørgensen [23] studied a differential game of waste management (disposal) between neighboring regions with strategic and stock externalities. Masoudi and Zaccour [24] considered a two-player differential game of international emissions to represent the interactions between developed and developing countries, and showed their cooperative and non-cooperative solutions. Wrzaczek, Shevkoplyas and Kostyunin [25] applied a differential game to formulate an overlapping generation model on optimal emissions with a continuous age structure. Bertinelli, Camacho and Zou [18] proposed a finite-horizon differential game to depict CO2 capture and storage and transboundary pollution and provided its explicit short-run dynamics with symmetric open-loop and a specific Markovian Nash strategy. Li [26] extended the model of Yeung [16] by taking into account emission permits trading, and found the optimal emission paths for the two regions. Gromova and Plekhanova [27] considered a differential game of transboundary pollution with developed and developing countries with a duration that is assumed to be exponentially distributed. Benchekroun and Martín-Herrán [28] studied a differential game of transboundary pollution between countries with myopic or non-myopic foresight. Huang, He and Hua [17] studied a cooperative differential game of transboundary industrial pollution between two asymmetric regions. Moreaux and Withagen [29] studied a differential game of carbon capture and storage from point sources with CO2 pollution and fossil fuel exhaustibility. Li [30] studied the behavior of the government about how to invest in emission abatement, how to regulate emission caps, and how to sell emission permits in a region divided into two economic sectors. Chang, Qin and Wang [31] involve the learning by doing in abatement in a transboundary pollution game. From a mathematical point of view, a transboundary pollution emission model can be used to approximately describe the transboundary GHG and other kinds of transboundary pollution emissions.
Although at first glance there are many similarities in the literature mentioned above, there are always some subtle differences from an economic point of view, mainly reflected in their objective functions and stock dynamics equations of carbon or pollution. In this study, to depict the complex scenarios of coordination strategies for international production planning and pollution abating under cap-and-trade regulations, we try to integrate the pollutant stock level, the abatement effort, the production quantity, the initial allocation of pollution emission permit and the emission permit trading price into the objective function and stock dynamics equation.

3. Modelling

All notations used through the paper are summarized in Table 2.
Although usually only two countries are involved in transboundary pollution emissions, a more general case is considered in this study. Assume n countries, labelled by i = 1 , 2 , , n , compete as pollution emission oligopolies with production activities creating pollutant as an undesirable by-product, and share the same pollution state. They also optimize tradeoffs between production planning and pollution abating under cap-and-trade regulations. The pollutant emitted by the n countries is assumed to be identical. Furthermore, e i ( t ) is assumed to be proportional to q i ( t ) , i.e.,
e i ( q i ( t ) ) = μ q i ( t ) .
The production quantity q i ( t ) can bring country i a certain level of net revenue, so that R i ( t ) can be written as
R i ( q i ( t ) ) = α i e i ( t ) 1 2 e i ( t ) 2 = α i μ q i ( t ) 1 2 μ 2 q i ( t ) 2 .
The production activity of each country generates not only net revenue but also a level of pollution with negative effects on its welfare. The environmental damage cost D i ( s ( t ) ) is supposed to satisfy
D i ( s ( t ) ) = 1 2 β i s 2 ( t ) .
As we know, the higher the pollutant concentration is, the more costly the damage is to the country. The welfare of one country depends not only on its own policies but also on those of other countries [32]. However their excessive production activities can cause increasing environmental pollution, and even have negative impacts on their welfare levels. While one country bears the pollution abatement costs, all countries can share the benefits from any reduction in emissions. It is reasonable for each country to do its best to opt for a free ride. Hence, one of the critical things for country i is to decide how much effort it should make on pollution abatement. Suppose the pollution abatement cost function G i ( u i ( t ) ) of country i is quadratic and depends on the effort on pollution abatement of the country
G i ( u i ( t ) ) = 1 2 σ u i 2 ( t ) ,
which means that the more pollution abatement effort the country makes, the more costly the pollution abatement is to the country.
Analogous to human beings, a country’s abatement action also has egoistic and altruistic sides, i.e., the former acts for the country’s self-interest and directly benefits the country’s revenue, and the latter is helpful for the others. The total pollution abatement amount of country i , also divided into egoistic and altruistic parts, is supposed to be proportional to its abatement effort, denoted by δ u i ( t ) . The egoistic part of pollution abatement amount B i ( u i ( t ) , s ( t ) ) of country i is simultaneously determined by its abatement effort and the pollutant stock level, as shown in (5)
B i ( u i ( t ) , s ( t ) ) = γ i u i ( t ) s ( t ) .
The trading revenue R i ( ) of country i is written as
R i ( q i ( t ) , s ( t ) , u i ( t ) ) = k ( F i + B i ( u i ( t ) , s ( t ) ) e i ( q i ( t ) ) ) = k ( F i + γ u i ( t ) s ( t ) μ q i ( t ) ) .
Apparently, the following three cases holds: R i ( ) > 0 if country i has permit surplus, and R i ( ) = 0 if the permit of country i is just the right amount, and R i ( ) < 0 if country i has permit shortage.
From (2)–(6), the instantaneous revenue of country i is given by (7) in the following
π i ( t ) =   R i ( q i ( t ) ) + R i ( q i ( t ) , s ( t ) , u i ( t ) ) G i ( u i ( t ) ) D i ( s ( t ) ) = μ α i q i ( t ) +   k ( F i + γ u i ( t ) s ( t ) μ q i ( t ) ) 1 2 ( μ 2 q i ( t ) 2 +   β i s 2 ( t ) +   σ u i 2 ( t ) )
The more goods the n countries produce, the more pollutants are emitted and, eventually, the higher the pollutant stock level is. However, each country can make its effort to decrease the pollutant stock level, and natural absorption or decomposition also reduces some pollutants. Suppose the decay function of pollutant stock is linear given by ε s ( t ) . Then the evolvement dynamics of the pollutant stock level is given by (8) in the following
s ˙ ( t ) = μ i = 1 n q i ( t ) δ i = 1 n u i ( t ) ε s ( t ) , s ( 0 ) = s 0 , s ( t ) 0 ,
In a non-cooperative differential game, each country independently chooses its production plan and pollution abatement to maximize its own discounted infinite-horizon utility as stated in (9) in the following
max q i ( t ) , u i ( t ) J i = 0 e r t π i ( t ) d t , i = 1 , 2 , , n ,
subject to (8).
In a cooperative differential game, n countries make joint decisions of the production plan and pollution abatement controlled by a central authority, such as the European Union, to maximize the sum of the discounted infinite-horizon utilities of all countries as stated in (10) in the following
max q i ( t ) , u i ( t ) J = 0 e r t i = 1 n π i ( t ) d t ,
subject to (8).
In the following, N C A , N C S , C A and C S in the superscripts indicate the asymmetric non-cooperative, symmetric non-cooperative, asymmetric cooperative and symmetric cooperative differential games, respectively. For expositional convenience, time-dependence and superscripts will be omitted from the notations if no confusion is caused.

4. Non-Cooperative Differential Games

4.1. Asymmetric Case (NCA)

4.1.1. Solutions of the Model

From (8) and (9), the HJB equation for country i is given in (11) in the following
r V i N C A = max q i , u i { π i + V i N C A s d s ( t ) d t } ,
where V i N C A = V i N C A ( s ) is the value function of country i .
The following optimal feedback production and abatement strategies of country i are obtained by taking first-order conditions of the HJB equation in (11) w.r.t. q i and u i , respectively
q i N C A   =   1 μ ( α i k + V i N C A s ) ,
u i N C A   =   1 σ ( γ i k s δ V i N C A s ) .
Although the coefficients of the value functions in the following propositions can be explicitly characterized, for simplicity and regularity, they are characterized implicitly for all, including the symmetric, asymmetric, cooperative and non-cooperative, scenarios.
Proposition 1.
In the asymmetric non-cooperative differential game, the Nash equilibrium solutions of the value functions and instantaneous levels of production and abatement, represented by V i N C A ( s ) , q i N C A ( t ) and u i N C A ( t ) , respectively, are
V i N C A ( s ) = η i 0 + η i 1 s + η i 2 s 2 ,
q i N C A   =   1 μ ( α i k + η i 1 + 2 η i 2 s ) ,
u i N C A   =   1 σ ( k s γ i δ ( η i 1 + 2 η i 2 s ) ) ,
where
η i 0 = 1 r ( η i 1 ( j = 1 n α j k n + σ + δ 2 σ ( j = 1 n η j 1 η i 1 2 ) ) + k F i + 1 2 ( k α i ) 2 ) ,
η i 1 = 1 r σ ( 2 η i 2 ( ( σ + δ 2 ) ( j = 1 n η j 1 η i 1 ) + σ ( j = 1 n α j k n ) ) η i 1 ( ε σ + k δ j = 1 n γ j 2 ( δ 2 + σ ) j = 1 n η j 2 ) ) ,
η i 2 = 1 2 σ r ( k 2 γ i 2 σ β i 4 η i 2 ( k δ j = 1 n γ j + ε σ + ( δ 2 + σ ) ( η i 2 2 j = 1 n η j 2 ) ) ) .
Proof of Proposition 1.
Substituting the optimal solutions from (12) and (13) into the HJB equation in (11) and simplifying yield
r V i N C A ( s ) = 1 2 ( ( k α i ) 2 β i s 2 ) + k F i + k 2 γ i 2 s 2 2 σ σ + δ 2 2 σ ( V i N C A s ) 2 + V i N C A s ( j = 1 n α j + j = 1 n V j N C A s ε s k n δ σ ( k s j = 1 n γ j δ j = 1 n V j N C A s ) )
Differentiating (14) w.r.t. s , the following is obtained
V i N C A ( s ) s = η i 1 s + 2 η i 2 s .
Substituting (14) and (21) into (20), then equating the coefficients of η i 0 + η i 1 s + η i 2 s 2 in (20), the results in (15) and (16) are obtained. □

4.1.2. Optimal Trajectory of the Pollutant Stock

The ordinary differential equation in (22) in the following is obtained by directly substituting (15) and (16) into (8)
s ˙ N C A ( t ) = C 1 N C A + θ N C A s N C A ( t ) ,
where C 1 N C A = σ + δ 2 σ j = 1 n η j 1 + j = 1 n α j k n and θ N C A = 2 ( σ + δ 2 ) σ j = 1 n η j 2 k δ σ j = 1 n γ j ε .
The optimal trajectory of the pollutant stock in (23) in the following is obtained by solving (22)
s N C A ( t ) = C 1 N C A θ N C A + C 2 N C A e t θ N C A ,
where C 2 N C A = s 0 + C 1 N C A θ N C A , and θ N C A < 0 ensures the stability of the steady state of the pollutant stock.

4.2. Symmetric Case (NCS)

4.2.1. Solutions of the Model

Remark 1.
For the symmetric case, let α i   =   α , β i   =   β , q i   =   q , u i   =   u , γ i = γ and F i   =   F for i = 1 , 2 , , n .
From (8) and (10), the HJB equation for any country can be written as
r V N C S = max q , u { π + V N C S s d s ( t ) d t } ,
where V N C S = V N C S ( s ) is the value function of the country.
The following optimal feedback production and abatement strategies of the country are obtained by taking the first-order conditions of the HJB equation in (24) w.r.t. q and u , respectively
q N C S = 1 μ ( α k + n V N C S s ) ,
u N C S = 1 σ ( k γ s n δ V N C S s ) .
Proposition 2.
In the symmetric non-cooperative differential game, the Nash equilibrium solutions of the value functions and instantaneous levels of production and abatement, represented by V N C S ( s ) , q N C S ( t ) and u N C S ( t ) , respectively, are
V N C S ( s ) = η 0 + η 1 s + η 2 s 2 ,
q N C S   =   1 μ ( α k + n η 1 + 2 n η 2 s ) ,
u N C S   =   1 σ ( k γ s n δ ( η 1 + 2 η 2 s ) ) ,
where
η 0 = 1 r ( k F + 1 2 ( α k ) 2 + n η 1 ( α k ) + n 2 η 1 2 2 σ ( σ + δ 2 ) ) ,
η 1 = 1 r ( 2 n η 2 ( ( α k ) + n η 1 σ ( σ + δ 2 ) ) η 1 ( n k γ δ σ + ε ) ) ,
η 2 = 1 2 σ r ( k 2 γ 2 + 4 n 2 η 2 2 ( σ + δ 2 ) 4 η 2 ( ε σ + n k γ δ ) β σ ) .
Proof of Proposition 2.
Substituting the optimal solutions from (25) and (26) into the HJB equation in (24) and simplifying yield
r V N C S ( s ) = k F + 1 2 ( ( α k ) 2 β s 2 ) + k 2 γ 2 s 2 2 σ + ( n α ε s k n n k γ δ s σ + n 2 ( σ + δ 2 ) 2 σ V N C S s ) V N C S s
Differentiating (27) w.r.t. s , the following is obtained
V N C S s = η 1 + 2 η 2 s .
Substituting (27) and (34) into (33) and equating the coefficients of η 0 + η 1 s + η 2 s 2 in (33), the results in (28) and (29) are obtained. □

4.2.2. Optimal Trajectory of the Pollutant Stock

The ordinary differential equation in (35) in the following is obtained by substituting (28) and (29) into (8)
V N C S s = η 1   +   2 η 2 s ,
where C 1 N C S = n ( α k ) + n 2 η 1 ( 1 + δ 2 σ ) and θ N C S = 1 σ ( 2 n 2 δ 2 η 2 k n γ δ ) + 2 n 2 η 2 ε .
The optimal trajectory of the pollutant stock in (36) in the following can be obtained by solving (35)
s N C S ( t ) =   C 1 N C S θ N C S + C 2 N C S e t θ N C S ,
where C 2 N C S   =   s 0 + C 1 N C S θ N C S , and θ N C S < 0 ensures the stability of the steady state of the pollutant stock.

5. Cooperative Differential Games

5.1. Asymmetric Case (CA)

5.1.1. Solutions of the Model

The following HJB equation for country i is obtained from (8) and (10),
r V C A = max q i , u i { i = 1 n π i + V C A s d s ( t ) d t } ,
where V C A = V C A ( s ) is the value function of country i .
The following optimal feedback production and abatement strategies for country i are obtained by taking the first-order conditions of the HJB equation in (37) w.r.t. q i and u i , respectively
q i C A = 1 μ ( α i k + V C A s ) ,
u i C A = 1 σ ( k s γ i δ V C A s ) .
Proposition 3.
In the asymmetric cooperative differential game, the optimal solutions of the value functions and instantaneous levels of production and abatement, represented by V C A ( s ) , q i C A ( t ) and u i C A ( t ) , respectively, are
V C A ( s ) = η 0 + η 1 s + η 2 s 2 ,
q i C A   =   1 μ ( α i k + η 1 + η 2 s ) ,
u i C A   =   1 σ ( k s γ i δ ( η 1 + η 2 s ) ) ,
where
η 0 = 1 r ( k j = 1 n F j ( k 1 ) j = 1 n α j + 1 2 j = 1 n α j 2 + k n ( k 2 ) ( 1 + δ 2 ) j = 1 n η j 1 2 + ( 1 + δ 2 σ ) j = 1 n η j 1 ) ,
η 1 = 1 r ( 2 σ ( σ + δ 2 ) j = 1 n η j 2 2 ( 1 + δ 2 ) j = 1 n η j 1 η j 2 ε k δ σ ( j = 1 n γ j + ( 1 σ ) j = 1 n γ j η j 1 ) ) ,
η 2 = 1 2 σ r ( 4 k δ ( σ 1 ) j = 1 n γ j η j 2 k 2 ( σ 2 ) j = 1 n γ j 2 σ ( j = 1 n β j + 4 ( 1 + δ 2 ) j = 1 n η j 2 2 ) ) .
Proof of Proposition 3.
Substituting the optimal solutions in (41) and (42) into the HJB equation in (37) and simplifying yield
r V C A ( s ) = j = 1 n α j 2 + k ( j = 1 n F j j = 1 n α j + k n 2 ) s 2 2 j = 1 n β j k 2 s 2 2 σ ( σ 2 ) j = 1 n γ j 2 + ( j = 1 n α j k n ε s + k δ s ( 1 2 σ ) j = 1 n γ j   +   n 2 ( 1 δ 2 σ ( σ 2 ) ) V C A s ) V C A s
Differentiating (40) w.r.t. s , the following is obtained
V C A s   =   η 1   +   2 η 2 s .
Substituting (40) and (47) into (46) and equating the coefficients of η 0 + η 1 s + η 2 s 2 in (46), the results in (41) and (42) are obtained. □

5.1.2. Optimal Trajectory of the Pollutant Stock

The ordinary differential equation in (48) in the following is obtained by substituting (41) and (42) into (8)
s ˙ C A ( t ) = C 1 C A + θ C A s ( t ) ,
where C 1 C A = j = 1 n α j k n + n η 1 ( 1 + δ 2 σ ) and θ C A = 2 n η 2 ( 1 + δ 2 σ ) k δ σ j = 1 n γ j ε .
The optimal trajectory of the pollutant stock in (49) in the following is obtained by solving (48)
s C A ( t ) = C 1 C A θ C A + C 2 C A e t θ C A ,
where C 2 C A = s 0 + C 1 C A θ C A , and θ C A < 0 ensures the stability of the steady state of the pollutant stock.

5.2. Symmetric Case (CS)

5.2.1. Solutions of the Model

The HJB equation for the countries in (50) in the following is from (8) and (9) together with Remark 1
r V C S = max q , u { n π + V C S s d s ( t ) d t } ,
where V C S = V C S ( s ) is the value function of the countries.
The following optimal feedback strategies of production and abatement of the countries are obtained by taking the first-order conditions of the HJB equation in (50) w.r.t. q and u , respectively
q C S = 1 μ ( α k + V C S s ) ,
u C S = 1 σ ( k γ s δ V C S s ) .
Proposition 4.
In the symmetric cooperative differential game, the optimal solutions of the value functions and instantaneous levels of production and abatement, represented by V C S ( s ) , q C S ( t ) and u C S ( t ) , respectively, are
V C S ( s ) = η 0 + η 1 s + η 2 s 2 ,
q C S = 1 μ ( α k + η 0 + η 1 s + η 2 s 2 ) ,
u C S   =   1 σ ( k γ s δ ( η 0   +   η 1 s +   η 2 s 2 ) ) ,
where
η 0 = n r ( k F + 1 2 ( k α ) 2 η 1 ( k α ) + η 1 2 2 ( 1 + δ 2 σ ) ) ,
η 1 = 1 σ r ( 2 n η 2 ( η 1 ( σ + δ 2 ) σ ( k α ) ) η 1 ( k n γ δ + ε σ ) ) ,
η 2 = n 2 σ r ( k 2 γ 2 β σ ) 2 η 2 σ r ( k n γ δ + ε σ + n η 2 ( σ + δ 2 ) ) .
Proof of Proposition 4.
Substituting the optimal solutions of (51) and (52) into the HJB equation in (50) and simplifying yield
r V C S ( s ) = 1 2 σ n k 2 γ 2 s 2 + n 2 ( 2 k F + ( k α ) 2 s 2 β ) + ( n α ε s k n σ ( γ δ s + σ ) + n 2 σ ( σ + δ 2 ) V C S s ) V C S s
Differentiating (53) w.r.t. s , the following is obtained
V C S s = η 1 + 2 η 2 s .
Substituting (53) and (60) into (59) and equating the coefficients of η 0 + η 1 s + η 2 s 2 in (59), the results in (54) and (55) are obtained. □

5.2.2. Optimal Trajectory of the Pollutant Stock

The ordinary differential equation in (61) in the following is obtained by substituting (54) and (55) into (8)
s ˙ C S ( t ) = C 1 C S + θ C S s C S ( t ) ,
where C 1 C S = 2 μ ( α k + n η 1 ( 1 + δ 2 σ ) ) and θ C S = 1 σ k n γ δ + ε 2 n 2 η 2 ( 1 + δ 2 σ ) .
The optimal trajectory of the pollutant stock in (62) in the following is obtained by solving (61)
s C S ( t ) = C 1 C S θ C S + C 2 C S e t θ C S ,
where C 2 C S = s 0 + C 1 C S θ C S , and θ C S < 0 ensures the stability of the steady state of the pollutant stock.

6. Numerical Simulation

In the following, numerical results of the case of a duopoly game, i.e., n =   2 , are demonstrated. The effects of the initial allocations of pollution emission permits F 1 and F 2 ( F ) and the pollution emission coefficient μ on the value functions V 1 and V 2 ( V ), on the production levels q 1 * and q 2 * , on the abatement efforts u 1 * and u 2 * ( u * ), and on the pollutant stock levels s ( t ) are studied through numerical simulation.
The selection of important parameter values in the models is similar to that in List and Mason [21]. The following example comes from global warming and U.S. carbon emissions (List and Mason, 2001). Nordhaus [33] derived that the carbon decay rate in the atmosphere per year is about 1%. Therefore, ε = 0.01 is used in the numerical study. He also estimated that the discount rate ranges from 1% to 7%. Hence, r = 0.02 is used in the numerical study. According to the ICE End of Day Reports (https://www.theice.com/market-data/end-of-day-reports), k = 21 dollars per ton is used. Based on empirical evidences [21,34,35,36,37], β 1 = 0.0051 , β 2 = 0.005 , α 1 = 27 , α 2 = 28 , γ 1 = 0.000001 , γ 2 = 0.0000011 , δ = 2 , σ = 0.1 , μ = 0.01 , s ( 0 ) = 167 (billion tons), F 1 = 56 (billion tons) and F 2 = 51 (billion tons) are used. In the NCS and CS scenarios, β = 0.0051 , α = 27 , γ = 0.000001 , and F = 56 are used. In all the plots in the following, t varies from 0 to 200 with an increment of 1.
The optimal solutions of the value functions, instantaneous levels of production and abatement, and the optimal trajectories of the pollutant stock levels are obtained by substituting the above parameter values into (14)–(16) and (23) for NCA, (27)–(29) and (36) for NCS, (40)–(42) and (49) for CA and (53)–(55) and (62) for CS, respectively. After simplification, these optimal solutions are presented in the following.
CaseNCA
  • s N C A ( t ) = 8.932 + 158.068 e 0.522583 t , q 1 N C A ( t ) = 589.56 0.633 s ( t ) ,
  • q 2 N C A ( t ) = 690.117 0.615 s ( t ) , u 1 N C A ( t ) = 2.088 + 0 . 127 s ( t ) ,
  • u 2 N C A ( t ) = 1.977 + 0 . 123 s ( t ) , V 1 N C A ( s ) = 1193.29 0.104 s ( t ) 0.00317 s 2 ( t ) ,
  • V 2 N C A ( s ) = 1094.84 0.099 s ( t ) 0.00307 s 2 ( t ) .
CaseNCS
  • s N C S ( t ) = 8.29039 + 158.71 e 0.519944 t , q NCS ( t ) = 589.56 0.633 s ( t ) ,
  • u NCS ( t ) = 1.876 + 0 . 124 s ( t ) , V NCS ( s ) = 1193.42 0.094 s ( t ) 0.0031 s 2 ( t ) .
CaseCA
  • s C A ( t ) = 0.485 + 166.515 e 0.900294 t , q 1 C A ( t ) = 584.678 1.08465 s ( t ) ,
  • q 2 C A ( t ) = 684.678 1.085 s ( t ) , u 1 C A ( t ) = 2.088 + 0 . 127 s ( t ) ,
  • u 2 C A ( t ) = 3 . 064 + 0 . 217 s ( t ) , V C A ( s ) = 2288.47 0.153 s ( t ) 0.005423 s 2 ( t ) .
CaseCS
  • s C S ( t ) = 0.451 + 166.549 e 0.895778 t , q C S ( t ) = 585.859 1.079 s ( t ) ,
  • u C S ( t ) = 2.828 + 0 . 216 s ( t ) , V C S ( s ) = 2387.12 0.141 s ( t ) 0.005396 s 2 ( t ) .

6.1. The Optimal Pollutant Stock Levels

Figure 1 shows the optimal evolution trajectories of pollutant stock levels of NCA, NCS, CA and CS as time t changes. It can be seen from Figure 1 that all pollutant stock levels decrease quickly at the beginning and finally reach low steady levels s NCA , s NCS , s CA and s CS , respectively, with s NCA = 8.932 , s NCS = 8.29039 , s CA = 0.485 and s CS = 0.451 .
Obviously, s NCA > s NCS > s CA > s CS holds, which means that the pollutant stock levels in symmetric cases are better than those in asymmetric cases, and those in cooperative cases are better than those in non-cooperative cases. It further shows that both symmetry and cooperation are helpful in reducing pollutant stock levels.

6.2. The Optimal Production Levels

Figure 2 shows the evolutionary trajectories of the optimal production levels of NCA, NCS, CA and CS as time t changes. It can be seen from Figure 2 that all optimal production levels increase quickly in the beginning and finally reach high steady levels q 1 NCA , q 2 NCA , q NCS q 1 CA , q 2 CA and q CS , respectively, with q 1 NCA = 583.906 , q 2 NCA = 684.624 , q NCS = 585.475 , q 1 CA = 584.153 , q 2 CA = 684.153 and q CS = 585.372 .
Apparently, q 1 NCA + q 2 NCA > 2 q NCS and q 1 CA + q 2 CA > 2 q CS mean that both the optimal production levels and the pollution emission levels in asymmetric games are more significant than those in symmetric games. Furthermore, q 1 NCA + q 2 NCA > q 1 CA + q 2 CA and q NCS > q CS show that both the optimal production levels and the pollution emission levels in non-cooperative games are higher than those in cooperative games. It further shows that both symmetry and cooperation are helpful in reducing pollution emission levels.

6.3. The Optimal Pollution Abatement Levels

Figure 3 shows the evolutionary trajectories of the optimal pollution abatement levels of NCA, NCS, CA and CS as time t changes. It can be seen from Figure 3 that all pollution abatement levels decrease quickly at the beginning and finally reach low steady levels u 1 NCA , u 2 NCA , u NCS , u 1 CA , u 2 CA and u CS , respectively, with u 1 NCA = 3.2207 , u 2 NCA = 3.07729 , u NCS = 2.90665 , u 1 CA = 3.16955 , u 2 CA = 3.16956 and u CS = 2.92573 .
Apparently, u 1 NCA + u 2 NCA > 2 u NCS and u 1 CA + u 2 CA > 2 u CS mean that both the optimal pollution abatement levels in asymmetric games are more significant than those in symmetric games. u NCS < u CS and u 1 NCA + u 2 NCA < u 1 CA + u 2 CA show that both the optimal pollution abatement levels in non-cooperative games are smaller than those in cooperative games. It further shows that symmetry goes against pollution abatement effort levels and cooperation is helpful to pollution abatement effort levels.

6.4. The Optimal Value Functions

Figure 4 shows the evolutionary trajectories of the optimal value functions of NCA, NCS, CA and CS as time t changes. It can be seen from Figure 4 that all value functions increase quickly in the beginning and finally reach high steady levels V 1 NCA , V 2 NCA , V NCS V 1 CA , V 2 CA , and V CS , respectively, with V 1 NCA = 1193.29 , V 2 NCA = 1094.84 , V NCS = 1193.42 , V CA = 2288.47 and V CS = 2387.12 .
Apparently, V 1 NCA + V 2 NCA < 2 V NCS and V CA < V CS mean that both the optimal value functions in asymmetric games are smaller than those in symmetric games. Furthermore, V 1 NCA + V 2 NCA < V CA and V NCS < V CS show that both the optimal value functions in non-cooperative games are smaller than those in cooperative games. It further shows that both symmetry and cooperation are helpful to the value functions.

7. Dividends Analysis

7.1. Cooperative Dividends

Definition 1.
The cooperative dividend is the difference between optimal payoffs of the cooperative game and the non-cooperative game.
Definition 2.
A cooperative agreement is feasible if and only if the total optimal cooperative payoffs are greater than the sum of the individual non-cooperative payoffs.
The so-called agreement feasibility refers to the collective rational choice of many parties. Definition 2 shows that a cooperative agreement is feasible if sum of individual payoffs of the cooperation strategy is more than that of the current non-cooperative strategy.
Figure 5 shows the evolutionary trajectories of the different optimal dividend levels as time t changes, in which CPsc denotes the cooperative dividend of a symmetric game, CPca denotes the cooperative dividend of an asymmetric game, CPcs denotes the symmetric dividend of a cooperative game, CPncs denotes the symmetric dividend of a non-cooperative game, and SCPcs denotes the cooperative and symmetric dividend of a game.

7.1.1. The Cooperative Dividend of an Asymmetric Game

Definition 3.
The cooperative dividend of an asymmetric game is the difference between optimal payoffs of its asymmetric cooperative scenario and its asymmetric non-cooperative scenario.
According to Definition 3, the cooperative dividend of an asymmetric game can be written as
C P c a = V C A ( V 1 N C A + V 2 N C A ) = 2.5804 150.372 e 1.80059 t + 155.917 e 1.04517 t 26.3881 e 0.9200294 t + 49.7445 e 0.522583 t .
Figure 5 also shows that the cooperative dividend of an asymmetric game decreases quickly at the beginning and finally reaches a low steady level C P c a = lim t C P c a = 2.5804 .

7.1.2. The Cooperative Dividend of a Symmetric Game

Definition 4.
The cooperative dividend of a symmetric game is the difference between optimal payoffs of its symmetric cooperative scenario and its symmetric non-cooperative scenario.
According to Definition 4, the cooperative dividend of a symmetric game can be written as
C P s c = V C S 2 V N C S = 2.20858 149.676 e 1.79156 t + 156.386 e 1.03989 t 24.3631 e 0.895778 t + 46.1037 e 0.519944 t .
Figure 5 also shows that the cooperative dividend of a symmetric game decreases quickly at the beginning and finally reaches a low steady level C P s c = lim t C P c a = 2.20858 .

7.2. Symmetric Dividends

Definition 5.
The symmetric dividend is the difference between optimal payoffs of the symmetric game and the asymmetric game.
Definition 6.
A symmetric agreement is stable if and only if the total optimal symmetric payoffs are greater than the total optimal asymmetric payoffs.

7.2.1. The Symmetric Dividend of a Non-Cooperative Game

Definition 7.
The symmetric dividend of a non-cooperative game is the difference between optimal payoffs of its symmetric non-cooperative scenario and its asymmetric non-cooperative scenario.
According to Definition 7, the symmetric dividend of a non-cooperative game can be written as
C P n s c = 2 V N C S ( V 1 N C A + V 2 N C A ) = 99.0347 + 155.917 e 1.04517 t 156.386 e 1.03989 t + 49.7445 e 0.522583 t 46.1037 e 0.519944 t .
Figure 5 also shows that the symmetric dividend of a non-cooperative game decreases quickly at the beginning and finally reaches a low steady level C P n s c = lim t C P n s c = 99.0347 .

7.2.2. The Symmetric Dividend of a Cooperative Game

Definition 8.
The symmetric dividend of a cooperative game is the difference between optimal payoffs of its symmetric cooperative scenario and its asymmetric cooperative scenario.
According to Definition 8, the symmetric dividend of a cooperative game can be written as
C P c s = V C S V C A = 98.6629 + 150.372 e 1.80059 t 149.676 e 1.79156 t + 26.3881 e 0.900294 t 24.3631 e 0.895778 t .
Figure 5 also shows that the symmetric dividend of a cooperative game decreases quickly at the beginning and finally reaches a low steady level C P c s = lim t C P c s = 98.6629 .

7.3. The Cooperative Symmetric Dividend of a Game

Definition 9.
The cooperative symmetric dividend of a game is the difference between optimal payoffs of its symmetric non-cooperative scenario and its asymmetric cooperative scenario.
Definition 10.
A cooperative symmetric agreement is stable if and only if the sum of the optimal cooperative symmetric payoffs is greater than the sum of the optimal asymmetric non-cooperative payoffs.
According to Definition 9, the cooperative symmetric dividend of a game can be written as
S C P c s = V C S ( V 1 N C A + V 2 N C A ) = 101.243 149.676 e 1.79156 t + 155.917 e 1.04517 t 24.3631 e 0.895778 t + 49.7445 e 0.522583 t .
Figure 5 also shows that the cooperative symmetric dividend of a game decreases quickly at the beginning and finally reaches a low steady level S C P c s = lim t S C P c s = 101.243 .
Obviously, S C P c s > C P n s c > C P c s > C P c a > C P s c holds, which means that the cooperative symmetric dividend of a game is the largest and the cooperative dividend of a symmetric game is the smallest.

8. Policies and Managerial Implications

Results mentioned above show that symmetry and cooperation are not only helpful in reducing both pollutant stock levels and pollutant emission levels but also helpful to payoffs of all countries. These results also prove the Paris agreement’s idea which is that as each country implements its pledge, others can learn what is feasible, and that collaborative global climate protection will emerge. So, we can preliminarily reach the following five necessary judgments.
First, the competition between countries is integral and an essential strategic factor. A country should consciously put competition factors into consideration in its global strategy.
Second, in global competition, a country needs to build its comparative advantage, which is the basis of cooperation. Only by taking the initiative in the competition can it better carry out international cooperation.
Third, from the perspective of global governance, the compromise between countries in international competition does not necessarily mean clandestine collusion against the value of other countries. On the contrary, it is likely to lead to a more effective allocation of global resources.
Fourth, the formation of regional alliances similar to the European Union will help in improving the symmetry of participants, and also be more conducive to production efficiency and pollution reduction.
At last, the future global climate agreements following the Kyoto Protocol and the Paris Climate Agreement should strive to enhance cooperation, symmetry, and transparency in international production and pollution reduction.

9. Conclusions

In this study, infinite-horizon cooperative and non-cooperative differential games are proposed to describe strategies of production planning and pollution abating of n countries under cap-and-trade regulations. Their asymmetric and symmetric feedback solutions are obtained using HJB equations.
These countries in the game simultaneously make two decisions—how much product to produce and how much pollution to abate under cap-and-trade regulations. In other words, each country makes decisions about how much pollution to abate or how much pollution emission right to buy or sell when deciding how much product to produce. Furthermore, the equilibrium solutions of the value functions, instantaneous levels of production and abatement, and the optimal trajectories of the pollutant stock levels vary at the beginning and finally reach their steady levels.
Both symmetry and cooperation are shown to have significant impacts on the production planning and pollution abating decisions. These impacts are:
  • Both symmetry and cooperation are helpful in reducing pollutant stock levels, abating pollution emission levels, and improving value functions.
  • Dividends increase in the following order: The cooperative dividend of a symmetric game, the cooperative dividend of an asymmetric game, the symmetric dividend of a cooperative game, the symmetric dividend of a non-cooperative game, and the cooperative and symmetric dividend of a game.
As concluded above, the cooperation of production plan and pollution emission significantly affect the decisions of the countries. However, how to share their information and how to improve their production cooperation among the n countries are interesting research topics. Moreover, equilibria mentioned above may vibrate under uncertain circumstances.
The following extensions are of interest for future research:
  • The deterministic differential games are based on assumptions of specific environmental damages which make the problem tractable. Further development is to consider uncertain environmental damages for which a stochastic differential game will be appropriate.
  • The ordinary differential games are based on the assumption that the environmental damage cost depends on the total pollutant stock level without considering its spatial diffusion. Partial differential games may be employed to consider the spatial diffusion.
  • The differential games are only constrained by the dynamical pollutant stock levels. If the pollution emission permits can be traded with dynamical prices, two or more dynamical constraints may be considered as shown in Xin and Sun [38].

Author Contributions

Conceptualization: B.X. and W.P.; formal analysis: B.X. and M.S.; visualization: B.X. and W.P.; writing—original draft: B.X.; validation: W.P. and M.S.; writing—review and edit: B.X. and M.S.; final manuscript preparation: B.X.; and all the authors have read and approved the final manuscript.

Funding

The work was funded by the National Planning Office of Philosophy and Social Science of China [No. 16FJY008]; the National Natural Science Foundation of China [No. 11801060]; and the Natural Science Foundation of Shandong Province [No. ZR2016FM26].

Acknowledgments

The authors would like to express sincere gratitude to the referees for their valuable suggestions and comments.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wang, Y.; Wang, Q.; Hang, Y.; Zhao, Z.; Ge, S. CO2 emission abatement cost and its decomposition: A directional distance function approach. J. Clean. Prod. 2018, 170, 205–215. [Google Scholar] [CrossRef]
  2. Santibanez-Gonzalez, E.D.R. A modelling approach that combines pricing policies with a carbon capture and storage supply chain network. J. Clean. Prod. 2017, 167, 1354–1369. [Google Scholar] [CrossRef]
  3. Pourhejazy, P.; Kwon, O.K. A Practical Review of Green Supply Chain Management: Disciplines and Best Practices. J. Int. Logist. Trade 2016, 14, 156–164. [Google Scholar] [CrossRef]
  4. Lee, P.T.-W.; Kwon, O.K.; Ruan, X. Sustainability Challenges in Maritime Transport and Logistics Industry and Its Way Ahead. Sustainability 2019, 11, 1331. [Google Scholar] [CrossRef] [Green Version]
  5. Pourhejazy, P.; Kwon, O.K.; Lim, H. Integrating Sustainability into the Optimization of Fuel Logistics Networks. KSCE J. Civ. Eng. 2019, 23, 1369–1383. [Google Scholar] [CrossRef]
  6. Victor, D.G.; Akimoto, K.; Kaya, Y.; Yamaguchi, M.; Cullenward, D.; Hepburn, C. Prove Paris was more than paper promises. Nature 2017, 548, 25. [Google Scholar] [CrossRef]
  7. Rogelj, J.; den Elzen, M.; Höhne, N.; Fransen, T.; Fekete, H.; Winkler, H.; Schaeffer, R.; Sha, F.; Riahi, K.; Meinshausen, M. Paris Agreement climate proposals need a boost to keep warming well below 2  °C. Nature 2016, 534, 631. [Google Scholar] [CrossRef] [PubMed]
  8. Zhang, B.; Lai, K.; Wang, B.; Wang, Z. The clean development mechanism and corporate financial performance: Empirical evidence from China. Resour. Conserv. Recycl. 2018, 129, 278–289. [Google Scholar] [CrossRef]
  9. Zhang, B.; Du, Z.; Wang, Z. Carbon reduction from sustainable consumption of waste resources: An optimal model for collaboration in an industrial symbiotic network. J. Clean. Prod. 2018, 196, 821–828. [Google Scholar] [CrossRef]
  10. Calvo, E.; Rubio, S.J. Dynamic Models of International Environmental Agreements: A Differential Game Approach. Int. Rev. Environ. Resour. Econ. 2013, 6, 289–339. [Google Scholar] [CrossRef] [Green Version]
  11. You, D.; Jiang, K.; Li, Z. Optimal Coordination Strategy of Regional Vertical Emission Abatement Collaboration in a Low-Carbon Environment. Sustainability 2018, 10, 571. [Google Scholar] [CrossRef]
  12. Feng, J.; Liu, B. Dynamic Impact of Online Word-of-Mouth and Advertising on Supply Chain Performance. Int. J. Environ. Res. Public Health 2018, 15, 69. [Google Scholar] [CrossRef] [PubMed]
  13. Schüller, K.; Staňková, K.; Thuijsman, F. Game Theory of Pollution: National Policies and Their International Effects. Games 2017, 8, 30. [Google Scholar] [CrossRef]
  14. Di Liddo, A. Price and Treatment Decisions in Epidemics: A Differential Game Approach. Mathematics 2018, 6, 190. [Google Scholar] [CrossRef]
  15. Ploeg, F.V.D.; Zeeuw, A.J.D. International aspects of pollution control. Environ. Resour. Econ. 1992, 2, 117–139. [Google Scholar] [CrossRef] [Green Version]
  16. Yeung, D. Dynamically Consistent Cooperative Solution in a Differential Game of Transboundary Industrial Pollution. J. Optim. Theory Appl. 2007, 134, 143–160. [Google Scholar] [CrossRef]
  17. Huang, X.; He, P.; Hua, Z. A cooperative differential game of transboundary industrial pollution between two regions. J. Clean. Prod. 2015, 120, 43–52. [Google Scholar] [CrossRef]
  18. Bertinelli, L.; Camacho, C.; Zou, B. Carbon capture and storage and transboundary pollution: A differential game approach. Eur. J. Oper. Res. 2014, 237, 721–728. [Google Scholar] [CrossRef]
  19. Long, N.V. Pollution control: A differential game approach. Ann. Oper. Res. 1992, 37, 283–296. [Google Scholar] [CrossRef]
  20. Breton, M.; Zaccour, G.; Zahaf, M. A differential game of joint implementation of environmental projects. Automatica 2005, 41, 1737–1749. [Google Scholar] [CrossRef]
  21. List, J.A.; Mason, C.F. Optimal Institutional Arrangements for Transboundary Pollutants: Evidence from a Differential Game with Asymmetric Players. J. Environ. Econ. Manag. 2001, 42, 277–296. [Google Scholar] [CrossRef]
  22. Yu, W.; Xin, B. Governance Mechanism for Global Greenhouse Gas Emissions: A Stochastic Differential Game Approach. Math. Probl. Eng. 2013, 312513. [Google Scholar] [CrossRef]
  23. Jørgensen, S. A dynamic game of waste management. J. Econ. Dyn. Control 2010, 34, 258–265. [Google Scholar] [CrossRef]
  24. Masoudi, N.; Zaccour, G. A differential game of international pollution control with evolving environmental costs. Environ. Dev. Econ. 2013, 18, 680–700. [Google Scholar] [CrossRef]
  25. Wrzaczek, S.; Shevkoplyas, E.; Kostyunin, S. A differential game of pollution control with overlapping generations. Int. Game Theory Rev. 2014, 16, 1450005–1450014. [Google Scholar] [CrossRef]
  26. Li, S. A Differential Game of Transboundary Industrial Pollution with Emission Permits Trading. J. Optim. Theory Appl. 2014, 163, 642–659. [Google Scholar] [CrossRef]
  27. Gromova, E.; Plekhanova, K. A differential game of pollution control with participation of developed and developing countries. Contrib. Game Theory Manag. 2015, 8, 64–83. [Google Scholar]
  28. Benchekroun, H.; Martín-Herrán, G. The impact of foresight in a transboundary pollution game. Eur. J. Oper. Res. 2015, 251, 300–309. [Google Scholar] [CrossRef]
  29. Moreaux, M.; Withagen, C. Optimal abatement of carbon emission flows. J. Environ. Econ. Manag. 2015, 74, 55–70. [Google Scholar] [CrossRef] [Green Version]
  30. Li, S. Dynamic optimal control of pollution abatement investment under emission permits. Oper. Res. Lett. 2016, 44, 348–353. [Google Scholar] [CrossRef]
  31. Chang, S.; Qin, W.; Wang, X. Dynamic optimal strategies in transboundary pollution game under learning by doing. Physica A Stat. Mech. Appl. 2017, 490, 139–147. [Google Scholar] [CrossRef]
  32. Breton, M.; Sbragia, L.; Zaccour, G. A Dynamic Model for International Environmental Agreements. Environ. Resour. Econ. 2010, 45, 25–48. [Google Scholar] [CrossRef]
  33. Nordhaus, W.D. To slow or not to slow: The economics of the greenhouse effect. Econ. J. 1991, 101, 920–937. [Google Scholar] [CrossRef]
  34. Nordhaus, W.D.; Yohe, G.W. Future paths of energy and carbon dioxide emissions. In Changing Climate: Report of the Carbon Dioxide Assessment Committee; National Academy Press: Washington, DC, USA, 1983; p. 87. [Google Scholar]
  35. Pearce, D.W.; Cline, W.R.; Achanta, A.N.; Fankhauser, S.; Pachauri, R.K.; Tol, R.S.J.; Vellinga, P. The social costs of climate change: Greenhouse damage and the benefits of control. In Climate Change 1995: Economic and Social Dimensions of Climate Change; Cambridge University Press: Cambridge, UK, 1996; pp. 179–224. [Google Scholar]
  36. Shogren, J. The Benefits and Costs of the Kyoto Protocol; American Enterprise Institute: Washington, DC, USA, 1999. [Google Scholar]
  37. Smith, J.; Tirpak, D. The Potential Effects of Global Climate Change on the United States: Report to Congress; US Environmental Protection Agency, Office of Policy, Planning, and Evaluation, Office of Research and Development: Washington, DC, USA, 1989.
  38. Xin, B.; Sun, M. A differential oligopoly game for optimal production planning and water savings. Eur. J. Oper. Res. 2018, 269, 206–217. [Google Scholar] [CrossRef]
Figure 1. Evolutions of optimal trajectories of the pollutant stock levels.
Figure 1. Evolutions of optimal trajectories of the pollutant stock levels.
Ijerph 16 03490 g001
Figure 2. Evolutions of the optimal production levels.
Figure 2. Evolutions of the optimal production levels.
Ijerph 16 03490 g002
Figure 3. Evolutions of the optimal pollution abatement levels.
Figure 3. Evolutions of the optimal pollution abatement levels.
Ijerph 16 03490 g003
Figure 4. Evolutions of the optimal value functions.
Figure 4. Evolutions of the optimal value functions.
Ijerph 16 03490 g004
Figure 5. Evolutions of different optimal dividends.
Figure 5. Evolutions of different optimal dividends.
Ijerph 16 03490 g005
Table 1. Some game models of transboundary emissions of CO2/pollution.
Table 1. Some game models of transboundary emissions of CO2/pollution.
StudyObjective FunctionStock DynamicsGame Type
Long (1992) 0 e r t ( R i ( e i ( t ) ) D i ( s ( t ) ) ) d t d s ( t ) d t = i = 1 2 e i ε s ( t ) Noncooperative
List and Mason (2001) 0 e r t ( R i ( e i ( t ) ) D i ( s ( t ) ) ) d t d s ( t ) d t = i = 1 2 e i ε s ( t ) Cooperative,
Noncooperative
Yeung (2007) 0 T e r t ( R i ( q i ( t ) ) G i ( u i ( t ) ) D i ( q i ( t ) , s ( t ) ) ) d t d s ( t ) d t = i = 1 2 μ i q i ( t ) j = 1 2 γ j u j ( t ) s ( t ) ε s ( t ) Cooperative
Masoudi and Zaccour (2013) 0 e r t ( R i ( e i ( t ) ) D i ( s ( t ) ) ) d t d s ( t ) d t = μ i = 1 2 e i ε s ( t ) Cooperative,
Noncooperative
Bertinelli, Camacho and Zou (2014) 0 T e r t ( G i ( s i ( t ) , u i ( t ) ) D i ( s ( t ) ) ) d t + e r T S ( s ( T ) ) d s ( t ) d t = i = 1 2 e i + μ i = 1 2 u i ε s ( t ) Non-cooperative
Li (2014) 0 e r t ( R i ( e i ( t ) ) + M i ( e i ( t ) ) D i ( s ( t ) ) ) d t d s ( t ) d t = i = 1 2 e i ε s ( t ) Cooperative,
Noncooperative
Gromova and Plekhanova (2015) 0 e r t ( R i ( e i ( t ) ) D i ( s ( t ) ) ) d t d s ( t ) d t = μ i = 1 2 e i ε s ( t ) Cooperative,
Noncooperative
Benchekroun and Martín-Herrán (2015) 0 e r t ( R i ( e i ( t ) ) D i ( s ( t ) ) ) d t d s ( t ) d t = i = 1 2 e i ε s ( t ) Non-cooperative
Huang, He and Hua (2015) 0 T e r t ( R i ( q i ( t ) ) G i ( u i ( t ) ) D i ( s ( t ) ) ) d t e r T S ( s ( T ) ) d s ( t ) d t = i = 1 2 q i ( t ) j = 1 2 γ j u j ( t ) s ( t ) ε s ( t ) Cooperative,
Noncooperative
This study 0 e r t ( R i ( q i ( t ) , s ( t ) , u i ( t ) ) G i ( u i ( t ) ) D i ( s ( t ) ) ) d t d s ( t ) d t = μ i = 1 n q i ( t ) δ i = 1 n u i ( t ) ε s ( t ) Cooperative,
Noncooperative
Note that notations in Table 1 are defined in Section 3.
Table 2. Notations.
Table 2. Notations.
SymbolsDescriptions
q i ( t ) the quantity of goods produced by country i in time t
e i ( t ) the amount of pollutant emitted by country i in time t
μ a pollution emission coefficient which represents the current technology level of clean production, μ > 0
R i ( t ) the net revenue function of country i in time t , which is strictly concave and quadratic, and depends on its pollution emission amount (List and Mason, 2001; Breton et al., 2005; Yeung, 2007; Li, 2014)
α i a positive constant for country i
D i ( s ( t ) ) the environmental damage cost to country i in time t
β i a positive scaling parameter for country i
s ( t ) a pollutant stock level in time t contributed by the emissions of the n countries
u i ( t ) the abatement effort of country i in time t
G i ( u i ( t ) ) the pollution abatement cost function of country i in time t
σ a positive scaling parameter
δ a positive scaling parameter, δ > 0
B i ( u i ( t ) , s ( t ) ) the egoistic part of pollution abatement amount of country i in time t
γ i a positive constant.
F i the initial allocation of pollution emission permit of country i
k ( t ) the emission permit trading price in time t
R i ( ) the trading revenue of country i
ε a natural absorption rate of pollutants, ε > 0

Share and Cite

MDPI and ACS Style

Xin, B.; Peng, W.; Sun, M. Optimal Coordination Strategy for International Production Planning and Pollution Abating under Cap-and-Trade Regulations. Int. J. Environ. Res. Public Health 2019, 16, 3490. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph16183490

AMA Style

Xin B, Peng W, Sun M. Optimal Coordination Strategy for International Production Planning and Pollution Abating under Cap-and-Trade Regulations. International Journal of Environmental Research and Public Health. 2019; 16(18):3490. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph16183490

Chicago/Turabian Style

Xin, Baogui, Wei Peng, and Minghe Sun. 2019. "Optimal Coordination Strategy for International Production Planning and Pollution Abating under Cap-and-Trade Regulations" International Journal of Environmental Research and Public Health 16, no. 18: 3490. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph16183490

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop