Next Article in Journal
The Condorcet Jury Theorem with Information Acquisition
Previous Article in Journal
Screening Teams of Moral and Altruistic Agents
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Double-Weighted Bankruptcy Method to Allocate CO2 Emissions Permits

1
LAMSADE, CNRS, Université Paris-Dauphine, Université PSL, 75016 Paris, France
2
LARODEC, Institut Superieur de Gestion de Tunis, University of Tunis, Tunis 2000, Tunisia
*
Author to whom correspondence should be addressed.
Submission received: 6 September 2021 / Revised: 12 October 2021 / Accepted: 14 October 2021 / Published: 23 October 2021
(This article belongs to the Section Applied Game Theory)

Abstract

:
Global warming, as a result of greenhouse gases, is exceeding the planet’s temperature stabilization capacities. Thus, greenhouse gas emissions must be reduced. We analyse a bankruptcy situation aimed at allocating emissions permits of CO2, the predominant greenhouse gas emitted by human activities. Inspired by the Constrained Equal Awards (CEA) solution for bankruptcy situations, we introduce a new allocation protocol based on the extension of the CEA solution over double-weighted bankruptcy situations, including two exogenous parameters aimed at providing a balance, in the request of emissions permits, between economic activities and the production of renewable energy. In these bi-criteria allocation problems, we focus on a computational approach to find an allocation protocol that does not prioritize any particular parameter. As an application of our method, we first consider CO2 permit allocation problems in European Union (EU) countries, using real data about the gross domestic product (GDP), the production rate of renewable energies, and countries’ ‘demands’ of CO2 emissions from 2010 to 2014. Then, we compare our approach with the CEA solution and its single-weighted extension to show the impact of using two weights over the distribution of CO2 emissions permits; we analyse the correlation between allocations of CO2 emission permits and the distribution of power within the EU Council to study the acceptability of alternative allocations.

1. Introduction

Climate change and global warming is a major challenge for the international community. Greenhouse gases, in particular carbon dioxide (CO2), accumulated in the atmosphere, cause irreversible damages, and the time required for their repairs could last hundreds of years. Emissions must be reduced by moving to clean and renewable energy resources, and the study of responsible sharing policies, in particular CO2 emissions, is one of the main goals of international climate policies nowadays. The literature on international cooperation on climate change has been organized around various international cooperation forms [1]. The existing methods proposed for emission allowances in CO2 can be classified into four groups: indicator approach based on one or more indicators, optimization approach, game theory approach, and hybrid approaches [2]. Certain allocation criteria are based upon equity principles, such as the grandfathering, consisting of giving an equal right to emit for all countries, the egalitarianism, giving to each person the same right to emit, or historical responsibility, according to which long-standing polluters should face a greater reduction rate, or other principles of distributive justice [3]. Other allocation criteria, where countries are compensated for reducing their emissions, follow a principle of efficiency, such as the approach proposed in [4], where the authors argue that efficiency can be defined as a type of “fairness”. In practice, alternative methods have been proposed to allocate CO2 emissions among countries and taking into account criteria, such as energy, the gross domestic product (GDP), population, energy production, etc. (see [2], for a recent survey), but no consensus has been reached on any particular allocation protocol.
However, the increasing amount of gas emissions and the limited absorption rate of the planet, on the one hand1, and the fact that countries, in 2020, have already consumed in 8 months2, all the resources that ecosystems can produce in a single year, it is obvious that we are facing a claims problem [5,6]. Claims problems, also called bankruptcy situations, are situations where a given quantity of a divisible resource must be allocated among a group of agents. Each agent has his own demand, and the total amount of divisible resources cannot satisfy the demand of all agents. Thus, the problem is to formulate a fair allocation rule or solution [7] (for a survey on bankruptcy situations, solutions, and their properties see, for instance, [8,9]). Bankruptcy situations can model different real ecological problems where a scarce resource, such as food, fishing quotas, medical suppliers, carbon budget, or water has to be shared among countries, or among the individuals of a population [10].
In this paper, we focus on the allocation of the most important greenhouse gas, carbon dioxide, whose emissions permits must be shared among European Union (EU) countries, and where each country has its own claim of CO2 emissions permits. This is not the first time that bankruptcy situations have been applied to study this problem at the EU level. For instance, two bankruptcy allocation solutions from the literature, namely, the Constrained Equal Awards (CEA) [7,8,11] and the Weighted Constrained Equal Awards (WCEA) [12,13], have been used in [14] to allocate CO2 emissions permits using EU countries’ CO2 claims and an exogenous weight based on the GDP. For other applications of bankruptcy situations to the allocation of emissions permits, see also [6]. As far as we know, however, this is the first study where the CO2 allocation problem has been investigated as a claims problem from a multi-criteria point of view. Specifically, in this paper, we introduce a preliminary application of a new allocation method for double-weighted bankruptcy situations based on the combination of two distinct criteria: the ability of countries to efficiently use their CO2 emissions permits, and the capacity of countries to produce energy in a sustainable manner via renewable sources. More precisely, using an extended version of the CEA solution for bankruptcy situations, we leverage the request of CO2 emissions permits of EU countries, taking into account both the economic growth of EU countries (in terms of their GDP) and their sustainable policies, measured as the amount of renewable energy produced by each country. We also compare the allocation of CO2 emissions permits provided by the methods considered in this paper with the distribution of power among the members of the EU Council, which is the main decision making body in the EU. As a main result of our analysis, we show that the double-weighted allocation protocol well represents the effective power of countries within the EU Council, and for this reason it could benefit from a more general consensus within the EU.
The remainder of the paper is organized as follows. The next section is devoted to the introduction of bankruptcy situations, and to the presentation of related basic notions and notation. Then, Section 3 introduces a new family of bi-criteria allocation methods for double-weighted bankruptcy situations, namely, the Double-Weighted Constrained Equal Awards (DWCEA) methods, and presents an algorithm to compute a specific allocation protocol. Section 4 is devoted to the application of the proposed protocol introduced in Section 3 to allocate CO2 emissions permits among EU countries over the years 2010–2014, and to the comparison of the results provided by the alternative solutions considered in this study. Section 5 concludes.

2. Preliminary

A bankruptcy situation or bankruptcy problem occurs when there is an infinitely divisible resource, the estate, to be divided among several agents having different claims and the same preferences, but there are natural upper limits on the allocation: none should be awarded more than his demand. Formally, a bankruptcy situation [15], or a claims problem, is defined by a tuple ( N , E , c ) (or, simply, ( E , c ) if the set N is already clearly identified), where N = { 1 , 2 , , n } is a set of agents, the estate E R > 0 is such that 0 < E < i N c i , and c R 0 N is a vector of agents’ claims. So, the problem is to divide E among the agents of N, each agent i N having her/his own claim c i 0 .
To allocate the estate E among the agents in N, we define an allocations vector as a real valued vector x R N , respecting the following properties:
  • rationality: x i 0 , for all i N ;
  • claim boundedness: x i c i , for all i N ;
  • efficiency: i N x i = E .
So, according to an allocation vector, every agent should receive a non-negative allocation, smaller or equal than her/his claim, and the entire estate should be divided completely among the agents. We are interested in defining a general allocation method, also called rule or solution, which is a function that associates to each bankruptcy situation an allocation vector. In the literature on bankruptcy situations, several rules have been proposed to allocate the estate (see, for instance, [7,8,9,11]).
In this work, we will focus on the Constrained Equal Awards (CEA) rule, and on some weighted version of the CEA rule. Formally, for each bankruptcy situation ( E , c ) with N as the set of players, the C E A rule yields an allocation vector C E A R 0 N , such that:
C E A i ( E , c ) = min { λ , c i } ,
for each agent i N and where λ R 0 and
i N min { λ , c i } = E .
This method can be defined as an iterative procedure. In the first iteration, all agents receive the same amount of estate, which must be smaller or equal than the smallest claim among the agents. Then, the agents with the smallest claim leave the game. A new iteration starts by equally allocating the remaining part of the estate, if any, among the remaining agents, and paying attention to allocate to each agent no more than the second smallest claim. Then, the agents claiming the second smallest claim leave the game, and the procedure is repeated among the remaining agents until there is no more estate to be shared.
A weighted bankruptcy situation [12,13], is a triple ( N , E , c , w ) where ( N , E , c ) is a bankruptcy situation and w = ( w 1 , w 2 , , w n ) is a vector of positive weights. The Weighted Constrained Equal Awards (WCEA) solution [12,13] is inspired by the CEA rule, and takes into account not only players’ claims, but also their weights, which modify the impact of the claims over the final allocation. The allocation vector provided by the WCEA rule is the unique vector W C E A R 0 N such that:
W C E A i ( E , c , w ) = min { c i , λ * w i } ,
for each agent i N and where the parameter λ * is such that
i N min { c i , λ * w i } = E
In the last section of this paper, we also make use of some concepts from cooperative game theory. Specifically, we compare allocation vectors for bankruptcy situations with the vector provided by the Shapley and Shubik (Sh–Sh) power index [16] for a simple game representing the voting rule at the EU Council [17]. A simple game is a pair ( N , v ) where N is a set of agents, or voters, and v is the characteristic function of the game that associates to each coalition S N a value v ( S ) { 0 , 1 } such that v ( N ) = 1 and v ( ) = 0 . The standard interpretation of a simple game is that v ( S ) = 1 means that coalition S is winning, according to some voting rule, while v ( S ) = 0 means that S is losing, according to the same voting rule. The Shapley and Shubik power index [16] is defined as the map ϕ associating to any simple game ( N , v ) a vector of numbers representing the P-power of voters [18], which is interpreted as the voters expected shares of a fixed prize for winning the elections, and is computed as the expected marginal contribution of each voter i N over all possible permutations of voters, i.e., ϕ i ( v ) = S N \ { i } s ! ( n s 1 ) ! n ! v ( S { i } ) v ( S ) , where s is the cardinality of coalition S.

3. Double-Weighted Constrained Equal Awards Rule (DWCEA)

Previous studies on bankruptcy situations and their solutions can be seen as an allocation approach where only one criterion, agent’s claim, is considered. Weighted bankruptcy situations [12] keep into account one more parameter to determine allocation vectors. In this section, we introduce a richer framework for bankruptcy situations, considering two vectors of weights at once, in addition to the vector of claims.
We denote by ( N , E , c , w 1 , w 2 ) a double-weighted bankruptcy situation on a set of agents N and where E R > 0 is a strictly positive estate, c = ( c i ) i N is the vector of non-negative claims such that i N c i > E , and w 1 = ( w i 1 ) i N and w 2 = ( w i 2 ) i N are two vectors of non-negative weights. We denote by DW the family of all double-weighted bankruptcy problems. An allocation rule ϕ : DW R 0 N is a map that associates to any double-weighted bankruptcy situation ( E , c , w 1 , w 2 ) DW an allocation vector in R N such that i N ϕ i ( E , c , w 1 , w 2 ) = E . Some interesting properties for an allocation rule for double-weighted bankruptcy situations are the following. Equal treatment: an allocation rule ϕ satisfies equal treatment if for all ( E , c , w 1 , w 2 ) DW , if i , j N are such that c i = c j , w i 1 = w j 1 and w i 2 = w j 2 , then ϕ i ( E , c , w 1 , w 2 ) = ϕ j ( E , c , w 1 , w 2 ) . Composition: An allocation rule ϕ satisfies composition if for all ( E , c , w 1 , w 2 ) DW such that w i 1 = w i 2 = 1 for all i N and for all 0 E E we have ϕ ( E , c , w 1 , w 2 ) = ϕ ( E , c , w 1 , w 2 ) + ϕ ( E E , c ϕ ( E , c , w 1 , w 2 ) , w 1 , w 2 ) .  Invariance under claims truncation: an allocation rule ϕ satisfies invariance under claims truncation if for all ( E , c , w 1 , w 2 ) DW we have ϕ ( E , c , w 1 , w 2 ) = ϕ ( E , c E , w 1 , w 2 ) , where c E R N is such that c i E = min { c i , E } for all i N .
In the attempt to further extend the WCEA rule to situations with two vectors of weights, in this paper, we introduce a Double-Weighted Constrained Equal Awards (DWCEA) solution, which associates to any double-weighted bankruptcy situation a particular allocation vector x R 0 N satisfying the following conditions:
x i ( E , c , w 1 , w 2 ) = min ( c i , λ 1 * w 1 , λ 2 * w 2 )
for each agent i N and with λ 1 * , λ 2 * R 0 such that
i N min ( c i , λ 1 * w 1 , λ 2 * w 2 ) = E
As it often happens in multi-criteria decision problems, in general, we face a wide choice of allocation vectors that satisfy the constraints specified by relations (5) and (6) for some λ 1 * and λ 2 * . Therefore, in order to select a specific allocation vector, in the next section we introduce a computational approach to find parameters λ 1 * and λ 2 * that do not require to assume any specific priority on the criteria represented by the two weight vectors.

A DWCEA Solution Algorithm with No Priority over Criteria

The algorithm studied in this section is based on the general principle that a DWCEA solution should grant the same importance to the different criteria represented by the two weight vectors. Based on this assumption, we introduce Algorithm 1 providing a solution for problem (5) under the constraint (6) with the objective to find feasible λ 1 * and λ 2 * without arbitrarily promoting the use of one of the two criteria to drive the choice of these parameters.
Algorithm 1. Double-Weighted method’s algorithm
Input: Estate E, Set of player N, Claims vector c, weight vectors w 1 , w 2
Output: an allocation vector x R 0 N
  ForEach i N : λ i 1 = c i w i 1 , λ i 2 = c i w i 2
   λ 1 S S o r t I n c r λ 1 , λ 2 S S o r t I n c r λ 2
   m o y = m e d i a n λ 1 S + m e d i a n λ 2 S 2
   λ 1 P r e v m o y , λ 2 P r e v m o y
  ForEach i N : x i = min ( c i , w i 1 λ 1 P r e v , w i 2 λ 2 P r e v )
   S = i N x i
  Begin Bloc 1:
  if S = E then
    λ 1 * = λ 1 P r e v
    λ 2 * = λ 2 P r e v
   ForEach i N :
    x i = m i n ( c i , w i 1 λ 1 * , w i 2 λ 2 * )
  end if
  End Bloc 1
   I t 0
  while I t | N | + 1 2 AND S E do
   if S > E then
    Algorithm A1
   end if
   if S < E then
    Algorithm A2
   end if
  end while
  return x
Algorithm 1 starts by calculating, for every player i N , a value λ 1 , i = c i w i 1 and a value λ 2 , i = c i w i 2 , then it sorts these values in incremental order, yielding vectors λ 1 S and λ 2 S . In the initial iteration, the same value m o y = m e d i a n λ 1 S + m e d i a n λ 2 S 2 is affected to parameters λ 1 P r e v and λ 2 P r e v , which are the temporary parameters used to compute the allocation. Starting the search for parameters λ 1 * and λ 2 * from the median m o y ensures a kind of neutrality for the importance of both criteria in the allocation computation, and it guarantees a minimal worst-case number of iterations before reaching the final λ 1 * and λ 2 * , since this procedure will be iterated at most N 2 + 1 times. At each iteration of the algorithm, the sum of the allocations, S = i N min ( c i , w i 1 λ 1 P r e v , w i 2 λ 2 P r e v ) , can take three possible values:
(1)
S = E , i.e., the sum of allocations is equal to the estate. In this case, we immediately have that λ 1 * = λ 1 P r e v and λ 2 * = λ 2 P r e v .
(2)
S > E , i.e., the sum of allocations is strictly larger than the estate. So, λ 1 and/or λ 2 must be decreased. In this situation, Algorithm 1 selects the first sorted λ 1 S < λ 1 P r e v as X and the first sorted λ 2 S < λ 2 P r e v as Y. Then, we choose the biggest value among X and Y to change the corresponding λ j P r e v value, j { 1 , 2 } , and we keep the same value for the other λ k j P r e v ; then, we recompute S and we reiterate the procedure until we obtain S = E , as shown in Algorithm A1 (see Appendix A).
(3)
S < E , i.e., the sum of allocations is strictly smaller than the estate, we have to increase λ 1 and/or λ 2 . So, λ j P r e v , j { 1 , 2 } , will be updated to the smallest value among X, the first λ 1 S > λ 1 P r e v , and Y, the first λ 2 S > λ 2 P r e v , and we keep the same value for the other λ k j P r e v ; then, we recompute S and we re-iterate the procedure until we obtain S = E , as shown in Algorithm A2 (see Appendix A).
After a number of iterations I t > N 2 + 1 , half of the agents’ demands are considered, and if the sum of allocations does not satisfy the efficiency constraint, E = i N min ( c i , w i 1 λ 1 P r e v , w i 2 λ 2 P r e v ) , the algorithm starts the “LastComputing” procedure. The same “LastComputing” procedure is also called at an iteration I t such that S < E and such that S > E at iteration I t 1 (or, inversely, at iteration I t such that S > E and such that S < E at iteration I t 1 ). Basically, as shown in the corresponding pseudo-code in Appendix A, the Last “LastComputing” procedure transforms the double-weighted bankruptcy situation with two weight vectors in a single weighted bankruptcy situation, using the parameter λ j P r e v , j { 1 , 2 } , which did not change since the previous iteration.
It is obvious that the solution provided by Algorithm 1 satisfies the equal treatment property. It also satisfies the composition property, as any solution satisfying conditions (5) and (6) when applied to a double-weighted bankruptcy problem ( E , c , w 1 , w 2 ) with w 1 = w 2 yields the allocation provided by the WCEA solution applied to the weighted bankruptcy problem ( E , c , w 1 ) , and the WCEA solution satisfies the composition property on the class of weighted bankruptcy problems (see [12,13]). Instead, we cannot guarantee that the solution provided by Algorithm 1 satisfies the property of invariance under claims truncation, as the use of vector c E instead of c may affect the procedure to compute λ 1 P r e v and λ 2 P r e v along the iterations of Algorithm 1.
We now introduce an example of calculations provided by Algorithm 1.
Example 1.
Consider the double-weighted bankruptcy situations ( E , c , w 1 , w 2 ) with N = { 1 , 2 , 3 } as the set of agents, E = 20 and the other parameters as shown in Table 1.
At the first iteration of Algorithm 1, the vectors λ 1 S = ( 2 , 5 2 , 5 ) and λ 2 S = ( 1 2 , 4 , 5 ) are defined to sort the ratios ( c i w i 1 ) i N and ( c i w i 2 ) i N , respectively, and the value m o y = 5 2 + 4 2 = 13 4 is calculated as the average of the two medians of λ 1 S and λ 2 S . This value m o y is assigned to the variables λ 1 P r e v and λ 2 P r e v . The provisional allocation for the three agents is computed as x 1 = min { 6 , 3 × 13 4 , 12 × 13 4 } = 6 for agent 1, x 2 = min { 10 , 2 × 13 4 , 2 × 13 4 } = 6.5 for agent 2 and x 3 = min { 20 , 8 × 13 4 , 5 × 13 4 } = 16.25 for agent 3. For the sum of allocations S = 28.75 is larger than the estate E = 20 , the procedure calls Algorithm A1 (see Appendix A). Now, at the new iteration ( I t = 1 in Table 2), the procedure selects the largest value between λ 1 S and λ 2 S which is also smaller than m o y = 13 4 , and such a value is used to set the new λ j P r e v for the corresponding weight j { 1 , 2 } , while the value λ k j P r e v is maintained equal to its value at the previous iteration. So, at I t = 1 , λ 1 P r e v = 5 2 and λ 2 P r e v = m o y = 13 4 . A new allocation for λ 1 P r e v and λ 2 P r e v is computed, and the procedure continues as detailed in in Table 2, till the sum of allocations becomes strictly smaller than the estate E = 20 (this happens at the iteration I t = 3 of Table 2. At this point, the largest (and closest to m o y ) value between the λ 1 P r e v and λ 2 P r e v is selected (so, λ 1 P r e v = 2 in the specific case), and the L a s t C o m p u t i n g process is called to calculate the final allocation corresponding to the WCEA allocation defined by relations (3) with claims vector ( min { c i , w i 1 λ 1 P r e v ) ) i N . The relevant parameters calculated at each iteration of Algorithm 1 are shown in Table 2.
Example 2.
Consider the double-weighted bankruptcy situations ( E , c , w 1 , w 2 ) with N = { 1 , 2 , 3 } as the set of agents, E = 30 and the other parameters as shown in Table 1. Similar to Example 2, we show all the iterations of Algorithm 1 in Table A1, Appendix B. Notice that at iterations I t = 0 and I t = 1 , it is called Algorithm A2, instead of Algorithm A1, as in Example 1, for S < E at I t = 0 and S > E at I t = 1 .

4. DWCEA Applied to CO2 Emissions Permits

In the context of global climate negotiations, an imperative step is to find an agreement or a strategy to be adopted for the reduction of CO2 emissions. However, CO2 emissions permits are limited and countries have to find a common consensus over methods to allocate emissions permits. In this paper, we analyse EU country claims on CO2 emissions permits, taking into account the limits in emission of CO2 recommended by the Kyoto protocol [19].
In order to retrospectively determine an allocation method that can be considered both efficient and equitable by EU countries, we focus on data provided by the World Bank Open Data project3. Precisely, for each of the five years from 2010 to 2014, we consider a double-weighted bankruptcy situation ( N , E y , c y , w y , 1 , w y , 2 ) , y = 2010 , , 2014 , with the 27 EU countries as the set N of agents, and the following features as vectors of claims and weights: the actual CO2 emission data from 2010 to 2014, as the vector of claims c y ; the quantity of GDP over the same time interval, as the first weight vector w y , 1 reflecting the economic growth rate of a country; finally, the production of renewable energy, as the second weight vector w y , 2 quantifying the sustainability of policies adopted by each country. Figure 1 summarizes the total fraction of these three parameters for the 27 EU countries during the period 2010–2014. According to the Kyoto protocol [19], which imposed by 2010 a reduction of 22 % of the CO2 total emissions in 1990, we set the estate E of each of the five bankruptcy problems considered over the years from 2010 to 2014, equal to the 78 % of the total amount of CO2 emitted by all EU countries in 1990.
We applied Algorithm 1 to each bankruptcy situation ( N , E y , c y , w y , 1 , w y , 2 ) to obtain the DWCEA allocations for years y = 2010 , , 2014 , and we computed the CEA allocations over the bankruptcy situations ( E y , c y ) as well as the WCEA allocations over the weighted bankruptcy problems ( E y , c y , w y , 1 ) and ( E , c , w y , 2 ) , respectively, for each year from 2010 to 2014. All of those allocations of CO2 emissions permits are reported in Table A2, Table A3, Table A4 and Table A5 in Appendix B.
As expected, the allocation yielded by the CEA solution completely satisfies small claims of emission permits. Instead, countries with high demands (i.e., Germany, UK, Italy and France), are drastically limited in their claims and receive the same amount of emission permits. The allocation provided by the WCEA solution based on GDP as the unique weight, favours countries with high GDP by giving more than half of the estate to the four countries with the highest GDP (i.e., Germany, UK, France, and Italy), while using the renewable energy as the unique weight, it completely satisfies some countries with intermediate emissions claims, (e.g., Denmark, Sweden, Finland, and Portugal), as well as countries having high claims and high renewable energy production (e.g., Germany and Spain).
Considering both the GDP and renewable energy as weights in Algorithm 1 to provide a DWCEA method, the countries respecting a specific threshold of both weights receive their total claims of emission permits, as it happens for Denmark, Sweden, Austria, and Spain in 2010, and for France and UK in 2014. Compared to the allocation generated by the WCEA solution based on GDP only as the weight, we can notice that the DWCEA allocation increases the amount of emission permits for 13 countries in 2010, 16 countries in 2011 and 2012, and 19 countries in 2013 and 2014. Instead, compared to the results obtained by the WCEA based on renewable energy as unique weight, we can observe an improvement in allocations for 17 countries for 2010 and 2011, and at least 14 countries from 2012 to 2014. Furthermore, in 2010, Denmark and Sweden were completely satisfied with whatever allocation method is used; Malta and Austria, in 2010, obtained an equal amount by DWCEA and WCEA based on renewable energy only, and Ireland, Finland, and France, in 2014, received the same amount of emissions permits according to DWCEA and WCEA based on GDP only.
We also notice similarities among allocations of some countries over the five studied years. To be more specific, we investigated those similarities by means of an unsupervised clustering technique, namely, the K-means method, based on the 1-distance notion to measure the similarity of each country’s allocation distribution over the years from each cluster’s reference centre [20]. This unsupervised clustering technique aims at grouping records (countries) of the data set into K distinct clusters according to their similarities, and each record can only be found in one cluster at a time. Note that countries’ emissions are naturally divided in small, medium, and high emission levels. For this reason, we applied the K-means method with K = 3 , to define three distinct groups reflecting the impact on CO2 emissions of different EU countries. The clustering on CO2 emission records over the five years is reported in Table 3(1) and shows a group G1 of countries with low emissions, a group G2 formed by countries with medium emissions and, finally, a group G3 containing high emissions countries.
The results of the clustering procedure using the CO2 emissions and the allocation distributions over the five years provided by the considered solutions CEA, WCEA with GDP only as weight, and WCEA with renewable energy only as weight, and DWCEA are shown in Table 3(2–5), respectively. Table 3(2) shows the clusters in accordance with the allocation distribution over the five years yielded by the CEA solution: notice that countries in groups G1 and G2 have low and medium demands and are totally satisfied with respect to their claims from 2010 to 2014. However, in Group 3, only the Netherlands, which has a medium emission distribution, is totally satisfied over 5 years. Other members of cluster G3 receive the same amount despite their different claims (except Spain, who has been totally satisfied since 2013).
Groups formed according to the allocation distribution over the five years by the WCEA solution based on GDP only, are presented in Table 3(3). We observe that group G1 is defined by countries with a lower GDP, group G2 by countries with intermediate GDP, and the last group, G3, is formed by countries with high GDP. Two countries from group G2 (i.e., Denmark and Sweden), and only one country from group G3 (i.e., France), are totally satisfied throughout the five years. So, these groups show a high association between emissions permits claims and GDP, whose effect plays in favour of countries with a large internal economic production.
In Table 3(4), concerning the clusters on allocation distributions over the five years provided by the WCEA solution based on the renewable energy only over the five years, group G3 contains countries with high emissions, which are also characterized by high levels of renewable energy production. In this group, Spain and Germany obtain their claims for all studied period. Even if the effect of the activities producing CO2 is still predominant, here, countries with a high production of renewable energy show a similar allocation distribution over time, and those that are fully satisfied for at least three years are in groups with intermediate and high claims and also high renewable energy production at the same time.
Finally, Table 3(5) presents the clusters generated on the allocation distribution over the years by the DWCEA allocation method. In this case, clusters are more homogeneous because they are generated by a balance among their claims, their economic levels, and their attitudes to employ green energy. According to the DWCEA allocation distribution over time, only Ireland, Denmark, Sweden, and Austria receive their full claims for at least three years due to their weight on renewable energy production. So, the trade-off between GDP and renewable energy production in the DWCEA allocation method seems to play in favour of countries characterized by intermediate levels of CO2 claims and GDP and, at the same time, good levels of renewable energy production.
The clusters provided by the K-means analysis can be interpreted as groups of countries equally affected by each proposed allocation method over the years. So, it seems reasonable to assume that countries in a cluster are guided by common interests and similar goals when faced with the possibility to accept or not a proposed allocation method. Therefore, we analysed the ability of these groups to influence the decision making process in the EU, and in particular on the voting process in the EU Council, which is the main collegiate body defining the overall political priorities of the EU, over the five year interval 2010–2014. For this reason, we focused on the ability of the different clusters of countries to impose a decision according to the voting system adopted during the same period by the EU Council. According to the EU rules (following the Treaty of Lisbon, effective since 2009, and operative for the voting rule of the EU Council since 2014), a decision is approved by the EU Council if it is supported by a coalition of at least 55 % of countries representing at least 65 % of the EU population (and keeping into account that a coalition may block a decision if it contains at least four countries globally representing at least 35 % of the EU population). Using the calculator provided by the EU Council4, we found the alternative combinations of clusters that may lead to the approval of a decision according to the voting system of the EU Council, which are summarized by Table 4.
Considering the clusters generated by K-means just on the emissions distributions over the five years, and assuming that all countries within a cluster may only cooperate within the same cluster, or together with all the countries in another cluster, a decision can be approved when all countries in cluster G1 and G3 cooperate, or when all countries cooperate together (the grand coalition). Concerning the clusters obtained on the CEA allocation distributions over the years, the approval of a new agreement is reached by the collaboration of group G3 with groups G2 or G1, respectively, or when all countries of the three groups cooperate to form a grand coalition. Looking at clusters formed with the allocation distributions yielded by the WCEA solution based only on GDP, we observe the same winning coalitions as for the CEA solution, whereas for clusters obtained over the WCEA solution based on renewable energy, an agreement could be reached either by the coalition of countries in clusters G1 and G2, or by all countries together. A similar configuration of winning coalitions of clusters is obtained with the allocation distribution provided by the DWCEA solution, where again the cluster G3 plays a less relevant role to form winning coalitions in the EU Council. So, even if countries in cluster G3 are never fully satisfied for at least three years according to the DWCEA solution, their possible common interest to increase their own allocation would be prevented by the voting mechanism of the EU Council, which makes cluster G3 non-decisive in forming a winning coalition. A similar argument holds for countries in cluster G1 that are never totally satisfied, as well as for the group formed by clusters G1 and G3 together. On the contrary, the group formed by G1 and G2 together could form a blocking coalition against the use of the DWCEA solution, but the fact that group G2 contains many countries that are totally satisfied using the DWCEA method, makes the formation of such a coalition less likely.
As we already observed, the CEA rule fully satisfies countries with a low claim of CO2 emissions permits, but it might encounter an objection to its application from countries with larger claims, who see their demands of CO2 emissions permits strongly reduced. In order to mitigate this effect, and improve the acceptability of an allocation by countries with larger claims, the DWCEA solution can be seen by EU governments and populations as a fair compromise, keeping into account both the efficiency in production of countries, as represented by the GDP, and their environmental impact, measured by the rate of renewable energy production. To measure the level of acceptability of the different allocations of CO2 emission permits, we compared the allocation vectors over the years 2010–2014 with the Sh–Sh power index [16] computed on the simple game representing the EU Council voting rule. This power index is used as a benchmark to represent the actual shares of power of the EU countries according to the voting rule established by the Lisbona treatment (the data of the Sh–Sh power index used in this paper are from the paper [17], Table A1 on page 7, and refer to the population data of 2015). It is in fact well established that the Sh–Sh power index in a simple game is an appropriate measure of the P-power of voters (see, for instance, [17,18]), which is a measure for evaluating the outcome that each voter in a simple game can expect before playing the game.
Precisely, to assess the association between the Sh–Sh power index distribution and the CO2 emission permits allocation vectors, we computed the Pearson correlation between the vector yielded by the Sh–Sh power index of the EU Council and the sum of CO2 allocations over the interval 2010–2014 for the most powerful EU countries having Sh–Sh power index larger or equal to 3 % of the total power, as well as for the remaining countries (see Table 5; the choice of the cut-off follows from the fact that countries with a Sh–Sh index larger or equal than 3 % form the smallest set of countries having two-thirds of the total power).
Figure 2 and Figure 3 show the distribution (in percentage) of total CO2 allocation vectors over the five years for the most powerful EU countries and the least powerful ones, respectively. In Figure 2, we observe the systematic cut operated by the CEA rule (dotted line) on the demand of CO2 permits for the six most powerful countries (precisely, for Poland, Spain, Italy, UK, France, and Germany), while the the adoption of the DWCEA rule shows an increase of the association between the allocation of CO2 and the Sh–Sh power index of more than 20 % in terms of the Pearson correlation (see Table 5 for the correlation coefficients of the different CO2 allocation vectors versus the Sh–Sh index vector). On the other hand, for the group of countries with Sh–Sh index smaller than 3 % , the loss of correlation with the Sh–Sh power index faced by the DWCEA allocation with respect to the CEA allocation is much less important, and is in the order of 10 % (see Table 5). So, we argue that the loss of acceptability of the DWCEA versus the CEA for the groups of less powerful countries, is largely compensated by the double gain of acceptability for the most powerful ones, and further justified by the use of criteria based on production efficiency and environmental preservation, which may further encourage small countries to adopt better technologies to improve their energetic policies.

5. Conclusions

The most critical objective in climate change study is to encourage countries worldwide to limit their emissions by finding a common strategy. In our work, we propose a new bi-criteria method in order to allocate CO2 emissions permits in EU countries. The DWCEA method, by giving the same importance to two parameters (GDP and production of renewable energy), ensures a balance between efficiency and green policies. Using basic clustering techniques on allocation records over the years 2010–2014, we showed that the DWCEA allocation method generates allocation distributions over time, having similar behaviours within clusters of countries characterized by specific levels of decision influence within the EU Council voting system, and prevent the cluster of larger countries from claiming extra shares of emissions permits.
We also showed that the allocation suggested by the DWCEA solution could encounter a larger level of acceptability, with respect to another solution, by the most influential or powerful EU countries (according to the well established Shapley and Shubik measure of power [16]), and without drastically affecting the acceptability by less influential countries. The justification of the trade-off in favour of the larger and most powerful countries is based on the use of criteria aimed at taking into account both production efficiency and environmental preservation.
We want to stress the fact that the GDP parameter, which we use extensively in this preliminary application to estimate the production of a country, is not necessarily the most appropriate attribute to represent the sustainable growth of a country. We leave for future research the use of alternative criteria, more oriented to reward an eco-friendly production.
Finally, future research devoted to a more detailed analysis of the axiomatic foundation of the allocation rule defined according to Algorithm 1 would be relevant. In particular, considering the variability of the parameters involved in each production process, it would be interesting to study the monotonicity property of the solution with respect to changes in countries’ claims, populations, and in the two weights describing the evolving environmental policies of countries.

Author Contributions

R.T. carried out algorithm conception and implementation, data acquisition, and data analysis. S.M. drafted a first version of the article. All authors contributed to the interpretation of the results. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data supporting reported results were essentially collected from the World Bank Open Data project: https://data.worldbank.org/indicator/ (accessed on 6 September 2021).

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

In this Appendix section, we provide the pseudo-code of all relevant routines and functions that are used by Algorithm 1 to compute an allocation according to the DWCEA solution introduced in Section 3.
Algorithm A1. Allocations’ sum is higher than Estate
Input: Estate E, Set of players N, Claims vector c, weight vectors w 1 , w 2
Output: an allocation vector x R 0 N
  while S > E and I t | N | + 1 2  do
   It ++;
   Bloc 2:
    X F i r s t λ 1 S
    Y F i r s t λ 2 S
   if X > Y then
     λ 1 P r e v X
   else
     λ 2 P r e v Y
   end if
   End Bloc 2
   ForEach i N : x i = min ( c i , w i 1 λ 1 P r e v , w i 2 λ 2 P r e v )
    S = i N x i
   if S = E then
     λ 1 * = λ 1 P r e v
     λ 2 * = λ 2 P r e v
    ForEach i N :
     x i = m i n ( c i , w i 1 λ 1 * , w i 2 λ 2 * )
   end if
   if S < E then
    LastComputing
   end if
  end while
  return x
Algorithm A2. Estate is higher than allocations’ sum
Input: Estate E, Set of players N, Claims vector c, weight vectors w 1 , w 2
Output: an allocation vector x R 0 N
  while S < E and I t | N | + 1 2  do
   It ++;
   Bloc 3:
    X F i r s t λ 1 S
    Y F i r s t λ 2 S
   if X < Y then
     λ 1 P r e v X
   else
     λ 2 P r e v Y
   end if
   End Bloc 3
   ForEach i N : x i = min ( c i , w i 1 λ 1 P r e v , w i 2 λ 2 P r e v )
    S = i N x i
   if S = E then
     λ 1 * = λ 1 P r e v
     λ 2 * = λ 2 P r e v
    ForEach i N :
     x i = m i n ( c i , w i 1 λ 1 * , w i 2 λ 2 * )
    end if
    if S > E then
    LastComputing
   end if
  end while
  return x
Algorithm A3. LastComputing Algorithm DWCEA
Input: Estate E, Set of players N, Claims vector c, weight vectors w 1 , w 2
Output: an allocation vector x R 0 N
  if ( λ 1 P r e v = L o w e r B o u n d λ 1 S ) OR ( λ 1 P r e v = U p p e r B o u n d λ 1 S ) then
    M m i n ( c , λ 1 p r e v w 1 )
    x W C E A ( E , M , w 2 )
  else
    M m i n ( c , λ 2 p r e v w 2 )
    x W C E A ( E , M , w 1 )
  end if
  return x
Algorithm A4. Function First λ j S for j = 1 or j = 2
Input: λ P r e v , λ j S : a vector of sorted λ j for j = 1 or j = 2
Output: X
  if S > E then
   %We are looking for the first λ < λ P r e v
    i 0
   while λ j S [ i ] < λ j P r e v do
     i + +
   end while
    X λ j S [ i 1 ]
  else
   %We are looking for the first λ > λ P r e v
    i N
   while λ j S [ i ] > λ j P r e v do
     i
   end while
    X λ j S [ i + 1 ]
  end if
  return X

Appendix B

Table A1. Parameters computed at each the iteration I t of Algorithm 1 on the double-weighted bankruptcy situation of Example 2, with λ 1 S { 2 , 5 2 , 5 } and λ 2 S { 1 2 , 4 , 5 } .
Table A1. Parameters computed at each the iteration I t of Algorithm 1 on the double-weighted bankruptcy situation of Example 2, with λ 1 S { 2 , 5 2 , 5 } and λ 2 S { 1 2 , 4 , 5 } .
I t = 0 : λ 1 P r e v = λ 2 P r e v = m o y = 5 2 + 4 2 = 13 4
Agent c i w i 1 × λ 1 P r e v w i 2 × λ 2 P r e v x i
16 3 × 13 4 = 9.75 12 × 13 4 = 39 6
210 2 × 13 4 = 6.5 2 × 13 4 = 6.5 6.5
320 8 × 13 4 = 26 5 × 13 4 = 16.25 16.25
S = 28.75 < E
I t = 1 : X 5 , Y 4 , λ 1 P r e v = m o y , λ 2 P r e v = 4
Agent c i w i 1 × λ 1 P r e v w i 2 × λ 2 P r e v x i
16 3 × 13 4 = 9.75 12 × 4 = 48 6
210 2 × 13 4 = 6.5 2 × 4 = 8 6.5
320 8 × 13 4 = 26 5 × 4 = 20 20
S = 32.5 > E
LastComputing Process: λ 1 * λ 1 P r e v = 13 4 , λ 2 * 7 2
Agent c i w i 1 × λ 1 * w i 2 × λ 2 P r e v x i
16 3 × 13 4 = 9.75 12 × 7 2 = 42 6
210 2 × 13 4 = 6.5 2 × 7 2 = 7 6.5
320 8 × 13 4 = 26 5 × 7 2 = 17.5 17.5
S = 20 = E
Table A2. Allocation vectors provided by CEA solution for each year in the interval 2010–2014 (values are in millions of tons of CO2).
Table A2. Allocation vectors provided by CEA solution for each year in the interval 2010–2014 (values are in millions of tons of CO2).
CountryCEA2010CEA2011CEA2012CEA2013CEA2014
Austria67.50265.0262.27362.48658.712
Belgium110.82499.94495.10796.9793.351
Bulgaria44.11449.34744.70839.642.416
Croatia20.17219.80917.99417.5516.843
Cyprus7.7087.4266.925.9486.062
Czech republic111.579106.908101.0398.67596.475
Denmark46.64140.64536.42838.53333.498
Estonia18.10818.60617.62419.89319.519
Finland62.08256.81649.13447.2247.301
France220.42227.373236.057242.088248.389
Germany220.42227.373236.057242.088248.389
Greece83.85779.84280.04369.48267.319
Hungary50.23347.84344.58342.14142.086
Ireland40.05535.63235.59234.85534.066
Italy220.42227.373236.057242.088248.389
Lativia8.0757.2947.0637.0816.975
Lithuania13.46913.78813.83212.6412.838
luxembourg10.96810.93910.66410.0519.659
Malta2.562.5412.6812.342.347
Netherlands183.053174.168170.31173.255167.303
Poland220.42227.373236.057242.088248.389
Portugal48.13747.62346.01445.42745.053
Romania79.41384.8881.72370.94570.003
Slovakia36.24134.52532.76533.09130.678
Slovenia15.33515.0914.78214.15112.812
Spain220.42227.373236.057237.035233.977
Sweden52.02451.73447.04844.84743.421
UK220.42227.373236.057242.088248.389
Table A3. Allocation vectors provided by the WCEA solution with GDP as a unique weight (denoted by WCEA1) for each year in the interval 2010–2014 (values are in millions of tons of CO2).
Table A3. Allocation vectors provided by the WCEA solution with GDP as a unique weight (denoted by WCEA1) for each year in the interval 2010–2014 (values are in millions of tons of CO2).
CountryWCEA 1 2010WCEA 1 2011WCEA 1 2012WCEA 1 2013WCEA 1 2014
Austria57.41759.4460.1462.48658.712
Belgium70.84672.6673.13469.48273.843
Bulgaria7.4157.9167.91811.2587.893
Croatia8.7668.68.32511.7288.018
Cyprus3.7453.7814.1315.9484.359
Czech republic30.39831.42830.46142.2828.912
Denmark46.64140.64536.42838.53333.498
Estonia2.8563.1953.3855.0753.25
Finland36.30637.73237.70747.2237.926
France353.033331.805333.228302.278303.276
Germany500.647518.086520.573485.25541.275
Greece43.8639.6836.08648.4332.976
Hungary19.18219.4118.78127.30119.494
Ireland32.51932.95433.13434.85534.066
Italy311.348313.84304.476334.097299.357
Lativia3.4813.8913.6784.8633.648
Lithuania5.4395.9946.2949.3726.75
luxembourg7.7968.2738.30910.519.228
Malta1.2811.3111.3532.0481.563
Netherlands122.542123.225121.76398.675122.378
Poland70.22772.91273.498105.84775.847
Portugal34.91433.76431.78227.30231.947
Romania24.41825.41925.21638.67527.754
Slovakia13.11313.53713.72219.8836.943
Slovenia7.0357.0726.8099.7156.943
Spain209.75205.165196.247173.255191.561
Sweden52.02451.73447.04844.84743.421
UK357.657361.187391.03363.443419.818
Table A4. Allocation vectors provided by the WCEA solution with renewable energy production as a unique weight (denoted by WCEA2) for each year in the interval 2010–2014 (values are in millions of tons of CO2).
Table A4. Allocation vectors provided by the WCEA solution with renewable energy production as a unique weight (denoted by WCEA2) for each year in the interval 2010–2014 (values are in millions of tons of CO2).
Country NameWCEA 2 2010WCEA 2 2011WCEA 2 2012WCEA 2 2013WCEA 2 2014
Austria67.50258.969034150.324247648.689983149.2450121
Belgium66.504410572.796992368.655486966.084277765.7096249
Bulgaria7.8639152522.121173814.280286916.64539815.3401465
Croatia1.8503332.27796482.88868253.820543025.1262437
Cyprus0.785315751.58389741.74680331.913196891.74733253
Czech republic33.488875847.027515540.346398438.468714140.0343097
Denmark46.6419.058469436.42838.53333.498
Estonia10.940631810.22414679.75355156.985801477.50746658
Finland62.08256.81649.13447.2247.301
France166.551486174.477866165.300608151.862197159.723832
Germany758.86732.498739.861756.900473719.883
Greece32.940230514.121602139.068581246.823592342.4486051
Hungary30.46594811.189612316.536857715.065694115.6488235
Ireland33.660999835.63230.273392629.423446931.3141833
Italy278.07708330.482862342.210321345.318320.411
Latvia1.237141251.6906772.73235383.639169614.4647929
Lithuania3.991125255.62372565.1656445.880008786.12944408
Luxembourg1.710482251.63728721.42055211.468539511.60401819
Malta0.010757750.0889830.17671940.204776430.41340675
Netherlands119.324963109.1376584.34952970.612762163.9181956
Poland85.71775296.1550298100.8795985.584844197.3435094
Portugal48.13747.62346.01445.42745.053
Romania4.4859817536.767775619.43913430.3771245.871613
Slovakia7.3690587510.86482439.31854998.8053862811.1344218
Slovenia2.485040252.82965942.9226672.814213162.86077471
Spain270.911270.548264.779237.035233.977
Sweden52.02415.683253847.04844.84743.421
UK239.037472256.728998243.602043284.205863323.526244
Table A5. Allocation vectors provided by the DWCEA solution for each year in the interval 2010–2014 (values are in millions of tons of CO2).
Table A5. Allocation vectors provided by the DWCEA solution for each year in the interval 2010–2014 (values are in millions of tons of CO2).
CountryDWCEA2010DWCEA2011DWCEA2012DWCEA2013DWCEA2014
Austria67.565.0259.286789644.809523158.712
Belgium72.6182.10442180.882747431.792056373.5489375
Bulgaria8.5910.011983736.42838.5337.86135324
Croatia2.022.569213033.4031454.172075657.9857891
Cyprus0.861.786405932.05790182.398283744.34209095
Czech republic36.5639.747320647.531926482.839800928.7973403
Denmark46.6410.216636219.48200429.2545158933.498
Estonia3.94.040155734.107462934.789233293.23685663
Finland49.5647.720568848.780845887.008947237.7754291
France181.84196.78566194.739968190.366523303.276
Germany683.42655.230259631.695144622.817849539.121412
Greece35.9615.927113622.90116439.810733532.8452471
Hungary26.1812.620255438.566374744.84719.4161707
Ireland36.7535.63235.59234.85534.066
Italy303.6372.736609369.469345.318298.165642
Latvia1.351.906837793.21897483.997471533.63399825
Lithuania4.366.342744666.0856247.370872066.72286212
luxembourg1.871.846621861.67354661.840884469.19093239
Malta0.010.10.20819240.256697041.55669538
Netherlands130.28123.09139899.37183434.755181121.891022
Poland93.5892.212500489.186524226.03166375.5455926
Portugal47.6642.7023719.607905518.885567731.8198291
Romania4.932.148096343.789469161.035221227.643745
Slovakia8.0512.253941810.978145411.037972613.9883644
Slovenia2.713.191444313.4431823.527750716.91533585
Spain270.91259.474424238.13784688.5164725190.798557
Sweden52.0217.688429547.04837.522125143.421
UK260.966289.548588286.982282356.26558418.879797

Notes

1
2
3
https://data.worldbank.org/indicator/ (accessed on 6 September 2021).
4

References

  1. Jin, Z.-G.; Cai, W.-J.; Wang, C. Simulation of climate negotiation strategies between China and the US based on game theory. Adv. Clim. Chang. Res. 2014, 5, 34–40. [Google Scholar] [CrossRef]
  2. Zhou, P.; Wang, M. Carbon dioxide emissions allocation: A review. Ecol. Econ. 2016, 125, 47–59. [Google Scholar] [CrossRef]
  3. Rose, A.; Zhang, Z.X. Interregional burden-sharing of greenhouse gas mitigation in the United States. Mitig. Adapt. Strateg. Glob. Chang. 2004, 9, 477–500. [Google Scholar] [CrossRef] [Green Version]
  4. Zhou, P.; Sun, Z.R.; Zhou, D.Q. Optimal path for controlling CO2 emissions in China: A perspective of efficiency analysis. Energy Econ. 2014, 45, 99–110. [Google Scholar] [CrossRef]
  5. Duro, J.A.; Giménez-Gómez, J.M.; Vilella, C. The allocation of CO2 emissions as a claims problem. Energy Econ. 2020, 86, 104652. [Google Scholar] [CrossRef] [Green Version]
  6. Giménez-Gómez, J.M.; Teixidó-Figueras, J.; Vilella, C. The global carbon budget: A conflicting claims problem. Clim. Chang. 2016, 136, 693–703. [Google Scholar] [CrossRef] [Green Version]
  7. O’Neill, B. A problem of rights arbitration from the Talmud. Math. Soc. Sci. 1982, 2, 345–371. [Google Scholar] [CrossRef] [Green Version]
  8. Thomson, W. Axiomatic and game-theoretic analysis of bankruptcy and taxation problems: A survey. Math. Soc. Sci. 2003, 45, 249–297. [Google Scholar] [CrossRef]
  9. Thomson, W. Axiomatic and game-theoretic analysis of bankruptcy and taxation problems: An update. Math. Soc. Sci. 2015, 74, 41–59. [Google Scholar] [CrossRef] [Green Version]
  10. Madani, K.; Zarezadeh, M. Bankruptcy methods for resolving water resources conflicts. In Proceedings of the World Environmental and Water Resources Congress: Crossing Boundaries, Albuquerque, Mexico, 20–24 May 2012. [Google Scholar]
  11. Aumann, R.J.; Maschler, M. Game theoretic analysis of a bankruptcy problem from the Talmud. J. Econ. Theory 1985, 36, 195–213. [Google Scholar] [CrossRef]
  12. Casas-Méndez, B.; Fragnelli, V.; García-Jurado, I. Weighted bankruptcy rules and the museum pass problem. Eur. J. Oper. Res. 2011, 215, 161–168. [Google Scholar] [CrossRef]
  13. Moulin, H. Priority rules and other asymmetric rationing methods. Econometrica 2000, 68, 643–684. [Google Scholar] [CrossRef]
  14. Trabelsi, R.; Moretti, S.; Krichen, S. Using Bankruptcy Rules to Allocate CO2 Emission Permits. In International Conference on Game Theory for Networks; Springer: Cham, Switzerland, 2019. [Google Scholar]
  15. Thomson, W. How to Divide When There Isn’t Enough. No. 62; Cambridge University Press: Cambridge, UK, 2019. [Google Scholar]
  16. Shapley, L.S.; Shubik, M. A method for evaluating the distribution of power in a committee system. Am. Political Sci. Rev. 1954, 48, 787–792. [Google Scholar] [CrossRef]
  17. Kóczy, L.Á. Brexit and Power in the Council of the European Union. Games 2021, 12, 51. [Google Scholar] [CrossRef]
  18. Felsenthal, D.S. A well-behaved index of a priori p-power for simple n-person games. Homo Oeconomicus 2016, 33, 367–381. [Google Scholar] [CrossRef]
  19. Vesterdal, M.; Svendsen, G.T. How should greenhouse gas permits be allocated in the EU? Energy Policy 2004, 32, 961–968. [Google Scholar] [CrossRef]
  20. Sinwar, D.; Kaushik, R. Study of Euclidean and Manhattan distance metrics using simple k-means clustering. Int. J. Res. Appl. Sci. Eng. Technol. 2014, 2, 270–274. [Google Scholar]
Figure 1. Percentage distribution of total EU CO2 emissions, GPD and renewable energy production of each EU country over the interval time of five years from 2010 to 2014 (all data were collected from the World Bank Open Data project in 2020, see https://data.worldbank.org/indicator/EN.ATM.CO2E.KT?view=chart (accessed on 6 September 2021)).
Figure 1. Percentage distribution of total EU CO2 emissions, GPD and renewable energy production of each EU country over the interval time of five years from 2010 to 2014 (all data were collected from the World Bank Open Data project in 2020, see https://data.worldbank.org/indicator/EN.ATM.CO2E.KT?view=chart (accessed on 6 September 2021)).
Games 12 00078 g001
Figure 2. Percentages of total CO2 emissions permits for the eight most powerful countries (Sh–Sh power index 3 % ) for the four allocation rules considered in this paper. Countries are ordered according to the Sh–Sh power index (solid line).
Figure 2. Percentages of total CO2 emissions permits for the eight most powerful countries (Sh–Sh power index 3 % ) for the four allocation rules considered in this paper. Countries are ordered according to the Sh–Sh power index (solid line).
Games 12 00078 g002
Figure 3. Percentages of total CO2 emission permits for EU countries with Sh–Sh power index strictly smaller than 3 % for the four allocation rules considered in this paper. Countries are ordered according to the Sh–Sh power index (solid line).
Figure 3. Percentages of total CO2 emission permits for EU countries with Sh–Sh power index strictly smaller than 3 % for the four allocation rules considered in this paper. Countries are ordered according to the Sh–Sh power index (solid line).
Games 12 00078 g003
Table 1. A double-weighted bankruptcy situation ( E , c , w 1 , w 2 ) with N = { 1 , 2 , 3 } as set of agents and estate E = 20 .
Table 1. A double-weighted bankruptcy situation ( E , c , w 1 , w 2 ) with N = { 1 , 2 , 3 } as set of agents and estate E = 20 .
Agent i c i w i 1 λ i 1 = c i w i 1 w i 2 λ i 2 = c i w i 2
163212 1 2
2102525
3208 5 2 54
Table 2. Parameters computed at each the iteration I t of Algorithm 1 on the double-weighted bankruptcy situation of Example 1, with λ 1 S { 2 , 5 2 , 5 } and λ 2 S { 1 2 , 4 , 5 } .
Table 2. Parameters computed at each the iteration I t of Algorithm 1 on the double-weighted bankruptcy situation of Example 1, with λ 1 S { 2 , 5 2 , 5 } and λ 2 S { 1 2 , 4 , 5 } .
I t = 0 : λ 1 P r e v = λ 2 P r e v = m o y = 5 2 + 4 2 = 13 4
Agent c i w i 1 × λ 1 P r e v w i 2 × λ 2 P r e v x i
16 3 × 13 4 = 9.75 12 × 13 4 = 39 6
210 2 × 13 4 = 6.5 2 × 13 4 = 6.5 6.5
320 8 × 13 4 = 26 5 × 13 4 = 16.25 16.25
S = 28.75 > E
I t = 1 : X 5 2 , Y 1 2 , λ 1 P r e v = 5 2 , λ 2 P r e v = m o y
Agent c i w i 1 × λ 1 P r e v w i 2 × λ 2 P r e v x i
16 3 × 5 2 = 7.5 12 × 13 4 = 39 6
210 2 × 5 2 = 5 2 × 13 4 = 6.5 5
320 8 × 5 2 = 20 5 × 13 4 = 16.25 16.25
S = 27.25 > E
I t = 2 : X 2 , Y 1 2 , λ 1 P r e v = 2 , λ 2 P r e v = m o y
Agent c i w i 1 × λ 1 P r e v w i 2 × λ 2 P r e v x i
16 3 × 2 = 6 12 × 13 4 = 39 6
210 2 × 2 = 4 2 × 13 4 = 6.5 4
320 8 × 2 = 16 5 × 13 4 = 16.25 16
S = 26 > E
I t = 3 : X 2 , Y 1 2 , λ 1 P r e v = 2 , λ 2 P r e v = 1 2
Agent c i w i 1 × λ 1 P r e v w i 2 × λ 2 P r e v x i
16 3 × 2 = 6 12 × 1 2 = 6 6
210 2 × 2 = 4 2 × 1 2 = 1 1
320 8 × 2 = 16 5 × 1 2 = 2.5 2.5
S = 9 , 5 < E
LastComputing Process: λ 1 * λ 1 P r e v , λ 2 * 2
Agent c i w i 1 × λ 1 * w i 2 × λ 2 * x i
16 3 × 2 = 6 12 × 2 = 24 6
210 2 × 2 = 4 2 × 2 = 4 4
320 8 × 2 = 16 5 × 2 = 10 10
S = 20 = E
Table 3. Outcomes of the K-means clustering application over allocation records provided by alternative allocation methods. Names in bold indicate countries that are totally satisfied at least for 3 years.
Table 3. Outcomes of the K-means clustering application over allocation records provided by alternative allocation methods. Names in bold indicate countries that are totally satisfied at least for 3 years.
Table 3(1): Clustering EU country emissionsMaps: EU country emissions
Group 1Group 2Group 3 Games 12 00078 i001
MaltaRomaniaSpain
CyprusGreecePoland
LatviaBelgiumFrance
LuxembourgCzech republicItaly
LithuaniaNetherlandsGermany
Slovenia
Estonia
Croatia
Slovakia
Ireland
Bulgaria
Denmark
Portugal
Hungary
Sweden
Finland
Austria
Table 3(2): Clustering CEA allocationsMaps: CEA allocations
Group 1Group 2Group 3 Games 12 00078 i002
MaltaRomaniaNetherlands
CyprusGreecePoland
LatviaBelgiumFrance
LuxembourgCzech republicItaly
LithuaniaDenmarkGermany
SloveniaBulgariaSpain
EstoniaPortugal
CroatiaHungary
SlovakiaSweden
IrelandFinland
Austria
Table 3(3): Clustering WCEA GDP allocationsMaps: WCEA GDP allocations
Group 1Group 2Group 3 Games 12 00078 i003
MaltaAustriaSpain
CyprusGreeceFrance
LatviaBelgiumItaly
LuxembourgCzech republicGermany
LithuaniaNetherlands
SloveniaPoland
EstoniaIreland
CroatiaDenmark
SlovakiaPortugal
BulgariaSweden
RomaniaFinland
Hungary
Table 3(4): Clustering WCEA renewable energyMaps: WCEA renewable energy
Group 1Group 2Group 3 Games 12 00078 i004
MaltaAustriaItaly
CyprusGreeceSpain
LatviaBelgiumGermany
LuxembourgCzech republic
LithuaniaNetherlands
SloveniaPoland
EstoniaFinland
CroatiaDenmark
SlovakiaPortugal
BulgariaSweden
HungaryFrance
IrelandRomania
Table 3(5): Clustering DWCEA allocationsMaps: DWCEA allocations
Group 1Group 2Group 3 Games 12 00078 i005
MaltaAustriaItaly
CyprusGreeceSpain
LatviaBelgiumGermany
LuxembourgCzech republicFrance
LithuaniaNetherlands
SloveniaPoland
EstoniaFinland
CroatiaDenmark
SlovakiaPortugal
BulgariaSweden
Hungary
Ireland
Romania
Table 4. Cluster combinations that are able to approve a decision according to the EU Council voting systems.
Table 4. Cluster combinations that are able to approve a decision according to the EU Council voting systems.
EmissionsCEAWCEA GDP WCEA Renewable DWCEA
G1 + G3G1 + G3G1 + G3G1 + G2G1 + G2
G1 + G2 + G3G2 + G3G2 + G3G1 + G2 + G3G1 + G2 + G3
G1 + G2 + G3G1 + G2 + G3
Table 5. Pearson correlation coefficients for CO2 allocation vectors provided by the four allocation rules considered in this paper over the interval 2010–2014, and for two groups of countries (groups are based on the Sh–Sh power index, cut-off 3 % ).
Table 5. Pearson correlation coefficients for CO2 allocation vectors provided by the four allocation rules considered in this paper over the interval 2010–2014, and for two groups of countries (groups are based on the Sh–Sh power index, cut-off 3 % ).
CEAWCEA GDP WCEA Ren DWCEA
8 EU countries (Sh–Sh 3 % )0.690.950.820.90
17 EU countries (Sh–Sh < 3 % )0.890.770.790.78
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Moretti, S.; Trabelsi, R. A Double-Weighted Bankruptcy Method to Allocate CO2 Emissions Permits. Games 2021, 12, 78. https://0-doi-org.brum.beds.ac.uk/10.3390/g12040078

AMA Style

Moretti S, Trabelsi R. A Double-Weighted Bankruptcy Method to Allocate CO2 Emissions Permits. Games. 2021; 12(4):78. https://0-doi-org.brum.beds.ac.uk/10.3390/g12040078

Chicago/Turabian Style

Moretti, Stefano, and Raja Trabelsi. 2021. "A Double-Weighted Bankruptcy Method to Allocate CO2 Emissions Permits" Games 12, no. 4: 78. https://0-doi-org.brum.beds.ac.uk/10.3390/g12040078

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