Next Article in Journal
Are Renewables as Friendly to Humans as to the Environment?: A Social Life Cycle Assessment of Renewable Electricity
Previous Article in Journal
Searching for New Directions for Energy Policy: Testing the Cross-Effect of Risk Perception and Cyberspace Factors on Online/Offline Opposition to Nuclear Energy in South Korea
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Carbon Emission Allocation in a Chinese Province-Level Region Based on Two-Stage Network Structures

1
Faculty of Mathematics and Statistics, Hubei University, Wuhan 430062, China
2
School of Economics and Management, Beijing University of Technology, Beijing 100000, China
3
Research Center for Eco-Environment Sciences Chinese Academy of Sciences, Beijing 100085, China
*
Author to whom correspondence should be addressed.
Sustainability 2019, 11(5), 1369; https://doi.org/10.3390/su11051369
Submission received: 2 February 2019 / Revised: 27 February 2019 / Accepted: 28 February 2019 / Published: 5 March 2019

Abstract

:
With the increasingly severe global environment and climate change, the growing social attention toward the environmental problems has prompted local governments to make policy adjustments. The formulation of the carbon emission right allocation scheme is important for policy-makers. Many researchers have studied the problem of carbon emission right allocation by using data envelopment analysis (DEA) models. However, the existing literature using traditional models consider each Decision-Making Unit (DMU) as a “black box” without taking the internal structure into account, but in fact, it is more accurate for formulating the scheme when considering the inner operation of DMUs. This paper investigates the allocation plan of carbon emission right among each province in China from 2007–2016 based on a two-stage DEA model. The results indicate that, first, there is no space for carbon emission in the north, northeast, and northwest from 2007–2016, while in the southern regions, it always exists. In addition, the carbon emission permits of the southern and eastern regions are increasing, but in the southwestern regions, the carbon emission space barely fluctuated during this decade. Second, the potential of carbon emission reduction of each region tends to be stable after 2014, and in the north and northwest, it fluctuated greatly from 2007–2016. Besides, the northwest region has had the potential of emission reduction since 2010, while it also exists in the northern region after 2014.

1. Introduction

Growing social concerns about the environment with the increasingly severe global environment and climate change are pushing the local government to make policy adjustments [1]. As one of the largest energy consumption and CO 2 emission countries [2,3], China has shown its determination to actively deal with the problem of climate change. At the 2009 Copenhagen climate conference, China made the commitment to reduce the carbon emission per unit of Gross Domestic Product (GDP) by 40–45% compared with 2005 by 2020 [4,5], then in the 2014 Sino-U.S. joint statement on climate to reduce carbon emission intensity by 60–65% compared with 2005 by 2030 [6]. At the end of 2017, the Chinese government announced the launch of a nationwide carbon trading market to achieve energy conservation and emission reduction targets. However, in the process of policy implementation, policymakers not only need to encourage provinces to achieve emission reduction targets, but also have to avoid disturbing market order [7]. Therefore, under the circumstances of the fixed amount of carbon emission right, how to formulate a reasonable and effective scheme of carbon emission right allocation is urgently needed based on the development and the amount of actual carbon emission of each province.
There are many international studies on the distribution of carbon emissions. In 1990, three major initial distribution schemes were proposed in the Amendment for the first time, namely public auction, fixed price sale, and free distribution. Some literature, such as Cramton et al. [8] and Kampas et al. [9], compared the characteristics of these schemes and concluded that each method has its own advantages in different situations. In 2005, the European Union Emissions Trading Scheme (EU ETS), which is the biggest emission trading market, was launched with the aim to achieve a given reduction target for aggregate CO 2 emissions at minimal cost [10,11]. Many scholars have investigated the carbon emission allocation rule in the EU ETS. Martin et al. [12] concluded that EU ETS had a robust negative impact on emission reduction, but no evidence supported the view that the EU ETS had strong detrimental effects on economic performance in exploring whether EU ETS had an impact on the economy. Oestreich et al. [11] and Daskalakis et al. [13] found that there exists a statistically-significant carbon premium in the study of the impact of EU ETS on the financial performance of enterprises. As the world’s largest carbon emitter, China faces greater pressure to reduce emissions. Therefore, it is meaningful to investigate the Chinese provincial carbon emission rights allocation. In addition, the experience of Chinese provincial carbon allocation research has significance to other countries in the world because the measures adopted by the Chinese government to allocate the carbon emissions right to each province are similar to those of the European Union Emissions Trading Scheme (EU ETS) to participating countries.
In recent years, Data Envelopment Analysis (DEA), which is an axiomatic non-parametric mathematical programming technique, has been increasingly used to allocate the fixed cost among homogeneous Decision-Making Units (DMUs) [14,15]. The rationale for its popularity in solving the fixed cost allocation problem is that it has several advantages. First, subjectivity can be avoided because there is no need to impose weights on inputs and outputs in advance [16,17]. Second, the DEA model has a high level of computational tractability and practicality since it is formulated and solved by linear programming [18,19]. Third, DEA provides decision-makers the possibility of considering the effect of feasible allocation plans on performance evaluation [20].
In China, there is a large regional imbalance between provinces. These imbalances include natural resources and the structure of economic and energy consumption [21]. Therefore, provinces should shoulder the burdens of carbon emission reduction jointly, but differentially [22,23]. The allocation of carbon emission reduction responsibility is actually the allocation of carbon emission right of each province since each province wants to take on less responsibility for emission reduction. Therefore, the allocation of carbon emission right is equivalent to the allocation of fixed cost.
In terms of allocation principles, DEA-based fixed cost allocation models in the existing literature can be classified into three categories. The first one is the principle of efficiency-invariance proposed by Cook and Kress [24], which assumes that the efficiency score of each DMU should remain unchanged after allocation. However, the method of cost allocation based on the efficiency invariance principle may be determined entirely by the input side of DMUs, namely the amount of cost allocation is the same for DMUs with the same inputs, but different outputs. The second one is the efficiency-maximization approach, first proposed by Beasley [25], which assumes that all DMUs would be efficient after the allocation under a set of common weights. However, the model is non-linear, and the gap between the maximum value and the minimum value in the final allocation scheme is large, which leads to difficulty in implementation. The last one is the game-based approach, first proposed by Nakabayashi and Tone [26]. This principle combines game theory with DMUs, that is the more one side gets, the less others get. In a recent study, Li et al. [27] developed a DEA-game cross-efficiency approach by integrating cooperative game theory and the cross-efficiency to generate a unique and fair allocation plan.
In previous literature, almost all literature dealing with carbon emission allocation with DEA was based on the above three principles. Gomes and Lins [28] proposed a Zero Sum Gains DEA (ZSG-DEA) model based on the input-oriented radial CCR-DEA model to allocate CO 2 emissions allowance among the Annex I parties and Non-Annex I countries of the Kyoto Protocol. Lozano et al. [29] proposed a DEA approach with three phases, i.e., maximizing aggregated desirable production, minimizing undesirable total emissions, and minimizing the consumption of input resources, to reallocate the emission permits. Wei et al. [30] estimated the CO 2 reduction potential and marginal abatement costs for 29 provinces over the period of 1995–2007 by using an extended Slacks-Based Measure (SBM) model. Wang et al. [31] proposed a non-radial ZSG-DEA model for emission allowance allocation on the provincial level for China by 2020. Miao et al. [23] also used a non-radial ZSG-DEA model to allocate CO 2 emissions between different Chinese provinces, but treating CO 2 as the undesirable output variable. Feng et al. [32] proposed a novel method that combined centralized models and compensation schemes for carbon emission abatement allocation. Zhou et al. [33] presented a DEA approach with multiple abatement factors based on the principles of equity and efficiency to allocate CO 2 emission quotas in 71 Chinese cities over 2005–2012. Yu et al. [34] considered the potential collaboration between industrial decision-making units and proposed a nonlinear DEA approach to allocate the regional industrial carbon abatement tasks.
Summarizing the related research, we find that all of the above studies considered each DMU as a “black box” without taking into account the internal structure. However, it is more accurate for formulating the scheme when considering the inner operation of DMUs. As a result, a new DEA model with a network structure was introduced by Färe [35]. There are some studies in the existing literature evaluating the efficiency of DMUs and allocating the fixed cost based on two-stage network DEA approaches, while little work has been done on the fixed cost allocation of practical problems by using the DEA-network model. Wanke et al. [36] adopted a network-DEA centralized efficiency model to measure the efficiency in Brazilian banking. Li et al. [37] used a three-stage DEA model to measure the efficiency of China’s manufacturing sector measures on green productivity growth during the 11th Five-Year Period. In addition to efficiency evaluation, some scholars have studied fixed cost allocation based on the two-stage DEA model in recent years. Yu et al. [38] proposed a fixed cost allocation method based on the two-stage network DEA and the cross-efficiency to allocate the fixed cost adequately. Zhu et al. [39] discussed the fixed cost allocation problem by using a two-stage network DEA method and proposed three procedures to obtain a fair cost allocation plan based on different objectives in reality. Ding et al. [40] first presented additive two-stage models to evaluate the performance for each DMU when allocating the fixed cost, and via introducing the concepts of satisfaction degree and fairness degree, an approach was proposed to obtain an optimal allocation plan. However, their approach cannot guarantee a unique allocation plan. Li et al. [41] first used the DEA model to measure the relative efficiency while taking the internal structure and possible allocated costs into account and proposed a unique allocation plan for all DMUs and sub-stages in view of the operation size of each DMU. Although the efficiency values of all DMUs are maximized and the allocation plan is unique, the algorithm of Li et al. [41] is not rigorous, since they neglected the use of the dual solution.
In this paper, we adopt the two-stage DEA-network model based on Li et al. [41] to investigate the allocation plan of carbon emission right among each province in China, and the algorithm for solving the unique allocation scheme is modified to make the result more accurate. To the best of our knowledge, none of the existing literature deals with carbon allocation based on the two-stage DEA-network model. Due to the complexity of the social economic structure in China, it is obviously inaccurate for the allocation plan to use the traditional DEA models that only consider the initial inputs and the final outputs of DMUs to deal with carbon emission right allocation. Nevertheless, a two-stage DEA-network model combines two stages via the intermediate variables, which means that the output of the first stage is the input of the second stage. The total fixed carbon emission will be allocated to each phase to participate in the whole production process. To sum up, the contributions of this paper are three-fold. Firstly, this paper applies the two-stage DEA-network model to investigate the allocation plan of carbon emission right. The two-stage DEA-network model is adopted to investigate the allocation plan of carbon emission right for the first time. Secondly, we point out that the algorithms solving the unique allocation scheme of the existing literature are not rigorous since they neglect the information of the dual solution, so the algorithms are modified in this paper to make the result more accurate. Finally, a fair, effective, and unique allocation plan is generated by combining the two-stage DEA-network model with the practical problem of carbon emission allocation right in China. By means of analyzing the allocation amount in the first and second stage, as well as the carbon emission space and emission reduction potential of each province, some guidance and suggestions can be provided for policy-makers.
The remainder of this paper is organized as follows. In Section 2, the methodology for carbon emission allocation is presented. In Section 3, an empirical analysis of the carbon emission allocation of each province is analyzed based on the two-stage DEA-network models proposed above, as well as the space and potential of carbon emission reduction of different regions in China. Then, some suggestions are given to policy-makers based on the results. Finally, Section 4 presents the main conclusions of the work.

2. Materials and Methods

2.1. Preliminary

DEA was first introduced by Charnes et al. [42] in 1978 and was called the CCR model, which assumed Constant Return to Scale (CRS). Subsequently, it was extended to the BCC model by Banker et al. [43], which assumed Variable Return to Scale (VRS), i.e., inputs and outputs cannot be increased or decreased proportionately. We will briefly describe the classical BCC model in this section.
Suppose there is a set of homogeneous DMUs, and each D M U j ( j = 1 , , n ) consumes m inputs x i j ( i = 1 , , m ) to produce s outputs y r j ( r = 1 , , s ) . The relative efficiency score for any given D M U d is calculated by solving the following BCC model:
θ d = M a x r = 1 s u r y r d + u 0 i = 1 m v i x i d s . t . r = 1 s u r y rj + u 0 i = 1 m v i x ij 1 , j = 1 , 2 , , n u r , v i 0 , r = 1 , 2 , , s ; i = 1 , 2 , , m
where u r and v i are unknown weights attached to the r th output and i th input, respectively, and the optimal objective function θ d is defined as the BCC efficiency score of D M U d ; where u 0 > 0 represents increasing return to scale, u 0 < 0 represents decreasing return to scale, and u 0 = 0 represents constant return to scale. Given u 0 = 0 , then the above model will be transformed to the CCR model.
Model (1) is nonlinear in a fractional programming form, and it can be transformed to a linear model by the Cooper–Charnes transformation. Let t = 1 / i = 1 m v i x i d , μ r = t u r , υ i = t v i , μ 0 = t u 0 . Then, the model becomes:
θ d = M a x r = 1 s μ r y r d + μ 0 s . t . r = 1 s μ r y r j i = 1 m υ i x i j + μ 0 0 , j = 1 , 2 , , n i = 1 m υ i x i d = 1 μ r , υ i 0 , r = 1 , 2 , , s ; i = 1 , 2 , , m
Solving Model (2), a series of weights ( μ r , υ i , μ 0 ) and the efficiency value of each DMU are obtained. The efficiency of D M U d is θ d = r = 1 s μ r y r d + μ 0 , which ranges from zero to one, and D M U d is identified as DEA efficient if θ d = 1 , otherwise DEA inefficient.

2.2. Carbon Emission Right Allocation Model Based on Two-Stage Structures

2.2.1. Two-Stage Network Structure

Consider the two-stage network structure shown in Figure 1.
In the first stage, the inputs of D M U j are x i j ( i = 1 , , m ) , and the outputs are the intermediate variables z p j ( p = 1 , , q ) . In the second stage, the intermediate variables z p j are taken as inputs to produce the final outputs y r j ( r = 1 , , s ) .
Suppose that there exists a total fixed cost R to be allocated, and the amount of allocation to each D M U j is R j , i.e.,
j = 1 n R j = R , R j 0 , j = 1 , 2 , , n
The allocation amount R j of each D M U j is allocated to the first stage and the second stage to participate in production activities with R 1 j and R 2 j , respectively, such that:
R j = R 1 j + R 2 j , R 1 j 0 , R 2 j 0 , j = 1 , 2 , , n
Based on the two-stage DEA model [41], we obtain the following model to calculate the overall efficiency when D M U d ( d = 1 , , n ) is under evaluation:
θ d = M a x p = 1 q φ p z p d + φ 0 + r = 1 s μ r y r d + μ 0 i = 1 m υ i x i d + υ m + 1 R 1 d + p = 1 q φ p z p d + υ m + 1 R 2 d s . t . θ 1 j = p = 1 q φ p z p j + φ 0 i = 1 m υ i x i j + υ m + 1 R 1 j 1 , j = 1 , 2 , , n θ 2 j = r = 1 s μ r y r d + μ 0 p = 1 q φ p z p j + υ m + 1 R 2 j 1 , j = 1 , 2 , , n j = 1 n R 1 j + R 2 j = R , R 1 j , R 2 j 0 , j = 1 , 2 , , n μ r , φ p , υ i 0 , r = 1 , , s ; p = 1 , , q ; i = 1 , , m
where θ 1 j and θ 2 j are the efficiency values of D M U j in the first stage and the second stage, respectively, and θ d is the overall efficiency of the D M U d to be evaluated. It is noteworthy that the overall efficiency θ d can be represented by a combination of θ 1 d and θ 2 d , namely θ d = λ d 1 θ 1 d + λ d 2 θ 2 d . λ d 1 and λ d 2 represent the relative importance of Sub-stages 1 and 2, respectively. It can be measured by the proportion of inputs in each stage on the total inputs, that is,
λ d 1 = i = 1 m υ i x i d + υ m + 1 R 1 d i = 1 m υ i x i d + υ m + 1 R 1 d + p = 1 q φ p z p d + υ m + 1 R 2 d , λ d 2 = p = 1 q φ p z p d + υ m + 1 R 2 d i = 1 m υ i x i d + υ m + 1 R 1 d + p = 1 q φ p z p d + υ m + 1 R 2 d
Model (5) is fractional programming, so we need to transform it to a linear model (6) (see Appendix A for details).
θ d = M a x ( p = 1 q ϕ p z p d + ϕ 0 + r = 1 s u r y r d + u 0 ) s . t . p = 1 q ϕ p z p j + ϕ 0 i = 1 m v i x i j r 1 j 0 , j = 1 , 2 , , n r = 1 s u r y r j + u 0 p = 1 q ϕ p z p j r 2 j 0 , j = 1 , 2 , , n i = 1 m v i x i d + r 1 d + p = 1 q ϕ p z p d + r 2 d = 1 j = 1 n ( r 1 j + r 2 j ) = v m + 1 R , r 1 j , r 2 j 0 , j = 1 , 2 , , n u r , ϕ p , v i 0 , v m + 1 > 0 , r = 1 , , s ; p = 1 , , q ; i = 1 , , m
For D M U d , the optimal weights of the model are ( u r d , ϕ p d , v i d , v m + 1 d , r 1 j d , r 2 j d , ϕ 0 d , u 0 d ) , and the optimal value can be calculated by ( θ d = p = 1 q ϕ p d z p d + ϕ 0 d + r = 1 s u r d y r d + u 0 d ) . In addition, the fixed cost allocation plan from the perspective of D M U d is R j d = R 1 j d + R 2 j d ( j = 1 , , n ) , where R 1 j d = r 1 j d / v m + 1 d and R 2 j d = r 2 j d / v m + 1 d .
Li et al. [41] proved that there exists at least a set of weights to make all DMUs and their sub-stage be simultaneously efficient. Therefore, an efficient allocation plan can be obtained as follows:
p = 1 q ϕ p z p j i = 1 m v i x i j R 1 j + ϕ 0 = 0 , j = 1 , 2 , , n r = 1 s u r y r j p = 1 q ϕ p z p j R 2 j + u 0 , j = 1 , 2 , , n j = 1 n ( R 1 j + R 2 j ) = R , R 1 j , R 2 j 0 , j = 1 , 2 , , n u r , ϕ p , v i 0 , r = 1 , 2 , , s ; p = 1 , 2 , , q ; i = 1 , 2 , , m
It can be seen from the efficient allocation plan (7) that the weights attached to the allocated cost are meaningless; for simplicity, the weight is usually set to one [25]. In the efficient allocation plan (7), there are ( m + q + s + 2 n + 2 ) variables, but only ( 2 n + 1 ) equality constraints. Therefore, the solution of (7) may not be unique.

2.2.2. Carbon Emission Right Allocation Based on the Two-Stage DEA Model

In this subsection, we obtain a unique allocation plan by taking the operation size of each DMU into consideration. Miao et al. [23] pointed out in the study of regional carbon emission allocation in China that the carbon emission allocation right based on per capita carbon emission would expand the gap between provinces. Specifically, for some developed regions, such as Beijing and Shanghai, their population is relatively less than that of underdeveloped provinces, but they produce more carbon emission. If the allocation is carried out according to per capita carbon emission, each province will not be efficient, and the disadvantage of allocating based on emission per unit of GDP is that it will restrict the development of some relatively underdeveloped provinces. Similar literature include Zhang et al. [44], Zhou et al. [7], and Kong et al. [45]. Therefore, various factors need to be taken into account to allocate the carbon emission rights of all provinces in China equitably and effectively.
This paper formulates a unique allocation plan based on the size of DMUs; it is reasonable and acceptable to allocate a large amount of carbon emission right to DMUs with large inputs and outputs and a small amount to DMUs with small inputs and outputs. Suppose the size parameters of the first stage and second stage are α j and β j , respectively. Then, the total amount of carbon emission allocation right of D M U j will be allocated to each stage by α j R , β j R , and j = 1 n ( α j + β j ) = 1 . The data processing and the operation size parameters are presented in Appendix B.
In order to obtain a unique allocation plan, we need to introduce a series of deviation variables that measure the difference between the efficient allocation R 1 j and R 2 j given in (7) and the proportionate allocation α j R and β j R .
Let α j R R 1 j = C j and β j R R 2 j = D j , that is C j represents the deviation between R i j and α j R for the first stage and D j represents the deviation between R 2 j and β j R for the second stage, and both C j and D j should be minimized. To this end, we minimize the maximum combined deviation ( C j + D j ) . The model (8) is as follows:
M i n M a x ( C k + D k ) s . t . α j R R 1 j = C j , j = 1 , 2 , , n β j R R 2 j = D j , j = 1 , 2 , , n p = 1 q ϕ p z p j i = 1 m v i x i j R 1 j + ϕ 0 = 0 , j = 1 , 2 , , n r = 1 s u r y r j p = 1 q ϕ p z p j R 2 j + u 0 , j = 1 , 2 , , n j = 1 n ( R 1 j + R 2 j ) = R , R 1 j , R 2 j 0 , j = 1 , 2 , , n u r , ϕ p , v i 0 , r = 1 , 2 , , s ; p = 1 , 2 , , q ; i = 1 , 2 , , m
In Model (8), the objective function is the deviation, and the value should be as small as possible. The first and the second constraints represent the deviation of each stage, and the subsequent constraints are the efficient allocation set. Since the model (8) is nonlinear, we need to change it into a linear model (9) (see Appendix C for details).
M i n M a x ρ s . t . c 1 j + c 2 j + d 1 j + d 2 j ρ j = 0 , j = 1 , 2 , , n ρ j ρ 0 , j = 1 , 2 , , n α j R R 1 j = c 1 j + c 2 j j = 1 , 2 , , n β j R R 2 j = d 1 j d 2 j , j = 1 , 2 , , n p = 1 q ϕ p z p j i = 1 m v i x i j R 1 j + ϕ 0 = 0 , j = 1 , 2 , , n r = 1 s u r y r j p = 1 q ϕ p z p j R 2 j + u 0 , j = 1 , 2 , , n j = 1 n ( R 1 j + R 2 j ) = R , R 1 j , R 2 j 0 , j = 1 , 2 , , n c 1 j , c 2 j , d 1 j , d 2 j 0 , j = 1 , 2 , , n u r , ϕ p , v i 0 , r = 1 , 2 , , s ; p = 1 , 2 , , q ; i = 1 , 2 , , m

2.2.3. Algorithm for Solving the Above Model

Suppose the optimal weights of the model (9) are ( ρ 1 , c 1 j 1 , c 2 j 1 , d 1 j 1 , d 2 j 1 , ρ j 1 , u r 1 , ϕ p 1 , v i 1 , R 1 j 1 , R 2 j 1 , ϕ 0 1 , u 0 1 , j , r , p , i ) . The allocation plans can be affected if we neglect the fact that the model (9) may have multiple optimal solutions and ignoring the information provided by the dual problem. Therefore, we need to use the dual programming of the model (9). The duality problem is formulated below.
m a x j = 1 n h j α j R + j = 1 n k j β j R + w R s . t . j = 1 n g j 1 f j + h j 0 , j = 1 , 2 , , n f j h j 0 , j = 1 , 2 , , n f j + k j 0 , j = 1 , 2 , , n f j k j 0 , j = 1 , 2 , , n f j + g j 0 , j = 1 , 2 , , n j = 1 n t j y r j 0 , r = 1 , 2 , , s j = 1 n l j z p j j = 1 n t j z p j 0 , p = 1 , 2 , , q j = 1 n l j x i j 0 , i = 1 , 2 , , m h j l j + w 0 , j = 1 , 2 , , n k j t j + w 0 , j = 1 , 2 , , n j = 1 n l j 0 j = 1 n t j 0 f j , g j , h j , k j , l j , t j 0 , j = 1 , 2 , , n w 0
Based on duality theory, when the optimal value of a dual variable is positive, the inequality constraints associated with this variable in the original problem will hold with equality [46]. By solving the dual model, the D M U j that satisfies the equation ρ j = ρ 1 will be determined based on the sign of the corresponding variable c j being positive. Then, all DMUs can be divided into two subsets:
Γ 1 = j c 1 j + c 2 j + d 1 j + d 2 j = ρ 1 , Γ 2 = j c 1 j + c 2 j + d 1 j + d 2 j < ρ 1
then minimizing the deviations of DMUs in Γ 2 . The process will not be terminated until the minimum deviation has been determined for all DMUs. The algorithm flow is described as follows:
Step 1: Let t = 1 , solving the model (9) to obtain the optimal weights ( ρ 1 , c 1 j 1 , c 2 j 1 , d 1 j 1 , d 2 j 1 , ρ j 1 , u r 1 , ϕ p 1 v i 1 , R 1 j 1 , R 2 j 1 , ϕ 0 1 , u 0 1 , j , r , p , i ) . If the deviation ρ j = c 1 j + c 2 j + d 1 j + d 2 j of D M U j is equal to ρ 1 , then we solve the dual model to determine whether the corresponding dual variable c j of D M U j is positive. If c j > 0 , then the D M U j is a member of set Γ 1 = j c 1 j + c 2 j + d 1 j + d 2 j = ρ 1 , and the allocation plan of D M U j can be determined as R 1 j and R 2 j ; otherwise, the D M U j with other D M U j whose deviations are not equal to ρ 1 will be regarded as members of set Γ 2 , i.e., Γ 2 = j c 1 j + c 2 j + d 1 j + d 2 j < ρ 1 . Denote the number of DMUs in Γ 1 as n 1 .
Step 2: If n 1 = m + q + s + 1 (the number of flexible variables in the efficient allocation set (7) is ( m + q + s + 1 ) ), then the algorithm terminates, and the final allocation plan is uniquely given by ( R 1 j , R 2 j , j ) . Otherwise, go to Step 3 if n 1 < m + q + s + 1 .
Step 3: Let t = t + 1 , solving the model (12) to obtain the optimal weight ( ρ ( t + 1 ) , c 1 j ( t + 1 ) , c 2 j ( t + 1 ) , d 1 j ( t + 1 ) , d 2 j ( t + 1 ) , ρ j ( t + 1 ) , u r ( t + 1 ) , ϕ p ( t + 1 ) , v i ( t + 1 ) , R 1 j ( t + 1 ) , R 2 j ( t + 1 ) , ϕ 0 ( t + 1 ) , u 0 ( t + 1 ) , j , r , p , i ) .
M i n ρ s . t . c 1 j + c 2 j + d 1 j + d 2 j ρ j = 0 , j = 1 , 2 , , n ρ j = ρ 1 , j Γ 1 ρ j = ρ ( t 1 ) , j Γ 2 t 1 ρ j ρ 0 , j Γ 2 t α j R R 1 j = c 1 j + c 2 j j = 1 , 2 , , n β j R R 2 j = d 1 j d 2 j , j = 1 , 2 , , n p = 1 q ϕ p z p j i = 1 m v i x i j R 1 j + ϕ 0 = 0 , j = 1 , 2 , , n r = 1 s u r y r j p = 1 q ϕ p z p j R 2 j + u 0 , j = 1 , 2 , , n j = 1 n ( R 1 j + R 2 j ) = R , R 1 j , R 2 j 0 , j = 1 , 2 , , n c 1 j , c 2 j , d 1 j , d 2 j 0 , j = 1 , 2 , , n u r , ϕ p , v i 0 , r = 1 , 2 , , s ; p = 1 , 2 , , q ; i = 1 , 2 , , m
If the deviation of D M U j is equal to ρ ( t + 1 ) , then we solve the dual programming of the model (12) in the same way to determine whether the corresponding dual variable c j of D M U j is positive. If c j > 0 , then the D M U j is a member of set Γ 2 t + 1 , and the allocation plan of D M U j can be determined as R 1 j ( t + 1 ) and R 2 j ( t + 1 ) ; otherwise, the D M U j with other D M U j whose deviations are not equal to ρ ( t + 1 ) will be regarded as members of set Γ 2 t + 2 . Then, the set Γ 2 t is divided into two subsets as follows: Γ 2 t + 1 = j c 1 j + c 2 j + d 1 j + d 2 j = ρ ( t + 1 ) , Γ 2 t + 2 = j c 1 j + c 2 j + d 1 j + d 2 j < ρ ( t + 1 ) . Denote the number of DMUs in Γ 2 t + 1 as n t + 1 .
Step 4: If n 1 + n 2 + + n t + 1 = m + q + s + 1 , then the algorithm terminates, and the final allocation plan is determined by ( R 1 j ( t + 1 ) , ( R 2 j ( t + 1 ) , j ) . Otherwise, go to Step 3 to continue.

3. Results and Discussion

The data source, allocation results, and some policy suggestions are presented in this section.

3.1. Variables and Data

By reviewing the principle of CO 2 emissions quota allocation in present studies, it can be seen that allocation principles of CO 2 emissions right are often based on economic ability, historical carbon emissions levels, and energy consumption, as well as some other indicators. For example, Wang et al. [31] considered population, total energy consumption, and CO 2 emissions as inputs and GDP as the only output. Feng et al. [32], Miao et al. [23], and Kong et al. [45] regarded labor, capital stock, and total energy consumption as inputs, GDP as the desirable output, and CO 2 emissions as the undesirable output. Jiang et al. [47] concluded that human economic activities in China were the dominant effect of carbon emission, while energy intensity and population growth were the most significant driving force. Wang et al. [48] and Ma et al. [49] indicated that the increase of carbon emission was basically promoted by aggressive economic output and increased energy consumption since the Chinese economy highly depends on energy consumption, and CO 2 emissions from energy consumption mainly came from the residential consumption sector and transportation industry. Therefore, it is reasonable to take the improvement of resident living standards in cities and transportation into account. In this paper, two indicators of urbanization level and private car ownership were selected as intermediate variables to measure the impact of residential consumption and transportation on carbon emission and the economy, respectively. In addition, like most literature, we use capital stock (K), labor force (L), and energy consumption (E) as inputs in the first stage, while GDP as the final desirable output. The following Figure 2 shows the detailed inputs and outputs of the system.
Dealing with undesirable outputs is also a problem that many scholars have studied in the DEA field. For example, the Directional Distance Function (DDF model, Chung et al. [50]) assumes that undesirable outputs are imposed with weak disposability and DMUs follow a predetermined direction to approach the effective frontier. The limitations of this model are that a directional vector needs to be specified beforehand, and a different final value results from different choices of directional vectors. The hyperbolic model (Färe et al. [51]) also assumes that the undesirable outputs are imposed with weak disposability. Unlike the DDF model, the locus of DMUs moving to the effective frontier in this model is hyperbolic, i.e., the inefficiency is measured by expanding the desirable outputs and reducing undesirable outputs by the same proportion. The limitation of this model is that it is difficult to solve since the model is nonlinear. The Seiford and Zhu model (SZ model, Seiford and Zhu [52]) assumes that undesirable outputs are imposed with strong disposability and transformed by adding a positive scalar on original undesirable outputs after multiplying them by −1. The limitation of this model is that the positive scalar is hard to specify since different transfer vectors will result in different efficiency values. The Undesirable Output as the Input model (UINP model, Hailu and Veeman. [53]) treats undesirable outputs as inputs, but this model cannot accurately represent the production process. In this paper, we adopted the approach consistent with most of the literature on fixed cost allocation dealing with undesirable outputs for simplicity, namely regarding undesirable outputs as inputs.
As an undesirable output, we regarded CO 2 emission as an input in the first stage since DMUs are expected to minimize their input consumption, as well as undesirable outputs.
The data covered 30 provinces during the period 2007–2016. However, Tibet was excluded due to the absence of energy consumption data. Data were collected from the China Statistical Yearbook 2007–2016. The description and processing of variables are as follows:
In terms of capital stock, this refers to all existing resources in provinces, and it is a comprehensive index reflecting the scale, speed, and structure of social investment (Wang et al. [54]). The capital stock is unavailable in any Statistical Yearbook. Therefore, we have to estimate this indicator via using the following method:
K i , t = I i , t + ( 1 δ t ) K i , t 1
where I i , t , δ i , and K i , t represent gross investment, depreciation rate, and capital stock for province i at time t, respectively. This paper adopted the method of Zhang et al. [55] to set δ = 9.6 % . Furthermore, we extended the capital stock of each province to 2016 based on 2005.
In terms of labor force, this is one of the indicators that must be considered in the production process. This paper selected the number of employees in the urban unit as the labor force indicator, since they are the main contributors of GDP and carbon emission.
As for energy consumption, this is the main contribution of China’s economic development. Besides, carbon emission is mainly generated by the burning of fossil energy. Therefore, energy consumption was taken as an input.
For carbon emission, it needed to be estimated since it was not available in existing data sources. At present, there are three international principles for carbon emission accounting [56]. The first one is producer responsibility, i.e., the responsibility of carbon emission belongs to the producers; the second principle is the consumer responsibility, which blames CO 2 emission on consumption; and the third one is a compromise between the first and the second principles, i.e., we need to consider the economic structure of the provinces when calculating carbon emission. Although the third one is more equitable, it was difficult for us to collect data because the China regional input-output table is investigated every five years. Therefore, this paper attributed carbon emissions to producers and believes that the carbon emissions of each province was mainly generated by the consumption of fossil energy in the province [57]. These can be estimated by multiplying the total amounts of energy consumptions by their corresponding carbon emission coefficients [54]. The carbon emission coefficients were obtained from the Intergovernmental Panel on Climate Change 2006. The calculation formula is as follows [30]:
C O 2 = i = 1 6 E i × C F i × C C i × C O F i × 44 12
where i represents different fossil energy indicators, which includes coal, gasoline, kerosene, diesel, fuel oil, and natural gas. E i represents the total energy consumption of type i. C F i represents the carbon emission factor of the ith energy. C C i is the net calorific value of the ith energy source. C O F i is the carbon oxidation factor, and 44 12 is the ratio of the mass of one carbon atom combined with two oxygen atoms to the mass of one oxygen atom. The net calorific value and carbon emission factor can be obtained from IPCC 2006, as shown in Table 1. The carbon oxidation factor of each energy defaults to one.
The urbanization level of each province can reflect people’s consumption patterns, which are closely related to carbon emissions and GDP. The index is measured by the ratio of urban population to total population.
With the upgrading of residents’ consumption patterns, the purchasing of private cars showed a trend of rapid growth. The ownership of private cars not only reflects economic development and consumption level, but also relates to carbon emission due to the change of travel mode.
This paper took GDP as the only desirable output. In case of the same input from each province, it is believed that the higher the GDP is, the higher the efficiency value the province obtains. Descriptive statistics of the input and output indicators for 30 provinces appear in Table 2, which contains the average level, the standard variance, the minimum level, and the maximum level. Note in Table 2 that the large standard deviation of each input index indicates that there is a great gap between provinces. Therefore, the VRS DEA model was adopted in this paper.
To facilitate the analysis, 30 provinces in China were divided into six regions including north, northeast, east, south, southwest, and northwest according to their geographical division. The north contained five provinces: Beijing, Tianjin, Hebei, Shanxi, and Inner Mongolia. The northeast included three provinces: Liaoning, Jilin, and Heilongjiang. The east included seven provinces: Shanghai, Jiangsu, Zhejiang, Anhui, Fujian, Jiangxi, and Shandong. The south included six provinces: Henan, Hubei, Hunan, Guangdong, Guangxi, and Hainan. The southwest involved five provinces: Chongqing, Sichuan, Guizhou, Yunnan, and Tibet. The northwest was comprised of five provinces: Shaanxi, Gansu, Qinghai, Ningxia, and Xinjiang. In order to put the different indicators of different regions together to compare the differences, this paper simultaneously reduced the asset stock, energy consumption, CO 2 emission, and GDP by 100-times because of the great difference between the indicators. The differences in input and output indicators among different regions are shown in Figure 3.
As shown in Figure 3, GDP was relatively large in the eastern and southern regions, while it varied little in the northern, northeastern, and southwestern regions, and it was smallest in the northwest. At the same time, other indicators such as labor force, private car ownership, and capital stock showed the same trend as GDP in each region. This phenomenon was caused by the relatively developed economy and technology in the east and south of China, while the relatively undeveloped situation in the northwest. Energy consumption and CO 2 emission in the north and east were higher, but in the southwest and northwest, they were less since economic development in the north and northeast mainly depends on heavy industry and energy utilization. As for the southwest and northwest, such as Yunnan, Gansu, and Ningxia, there were less energy consumption and CO 2 emissions because of the sparse population and tourism. The urbanization rate was higher in the north, northeast, and east, while it was lower in the south, southwest, and northwest, but they all exceeded 0.5. According to the analysis above, the efficiency evaluation of each province needs to take all the above indicators into account.

3.2. Results Analysis

In order to analyze the carbon emission right allocation of each province, this paper took data into the model (2) and the model (6), respectively, to obtain the efficiency before allocation by setting the fixed cost allocation R = 0 , and the results are shown in Table 3.
Experimental results show that when using the BCC model for efficiency evaluation, the efficiency of 19 provinces reached one, and all provinces exceeded 0.9 except for seven provinces, including Shanxi, Liaoning, Heilongjiang, Shaanxi, etc. While using the two-stage network model to evaluate the efficiency of provinces, there only Hainan and Qinghai provinces reached the efficiency value of one, followed by Ningxia, Beijing, Shanghai, Tianjin, Guangdong, Zhejiang, etc. This is because each DMU will chose a weight that was most advantageous to itself to achieve maximum efficiency when using the BCC model for efficiency evaluation [58], videlicet the traditional BCC model usually overestimated the efficiency of DMUs and neglected the internal structure of DMUs via only considering the inputs and outputs. By contrast, with some inputs as internal variables, the two-stage network model not only reflected the efficiency of each phase, but also the actual production situation.
This paper first used the method mentioned above to analyze the carbon distribution in each province in 2016. In 2016, the total carbon emission, namely fixed cost allocation, was 978,536.52 tons. Firstly, the model (6) was used to calculate the efficiency of the two subsystems and the overall efficiency after allocation. The result shows that the sub-stages and overall efficiency values of each province were one. Then, we used the model (9) to get the amount of allocation of each subsystem, and finally, the overall allocation was obtained, as shown in Table 4.
The second column of Table 4 shows the actual carbon emissions of the provinces in 2016, among which the top five provinces for carbon emission were Shandong, Inner Mongolia, Shanxi, Jiangsu, and Hebei. In other words, carbon emission in the northern and eastern regions was relatively higher due to more energy consumption in the northern area and more labor force in the eastern. Based on the two-stage network DEA model proposed in this paper, the top five provinces with the highest carbon allocation were Jiangsu, Shandong, Guangdong, Zhejiang, and Henan, which indicates that the allocation plan was not only related to carbon emissions, but also related to factors such as labor force, asset stock, urbanization rate, and private car ownership.
The last column in the Table 4 is the space of carbon emission that was attained based on the difference between the carbon allocation and actual emission of each province. The space of carbon emission reflects the maximum amount of CO 2 that each province can emit. If the value is negative, this indicates that the corresponding province should improve energy efficiency and reduce CO 2 emission. Otherwise, it demonstrates that the province has enough carbon emission right. The results of carbon emission permits of each province in 2016 are shown in Figure 4.
As shown in Figure 4, provinces with negative carbon emission permit value included Hebei, Shanxi, Inner Mongolia, Liaoning, Jilin, Heilongjiang, Shandong, Hainan, Guizhou, Yunnan, Shaanxi, Gansu, Qinghai, Ningxia, and Xinjiang. In particular, some provinces located in the northern and western regions such as Shanxi, Inner Mongolia, Guizhou, Shaanxi, Ningxia, and Xinjiang had large negative corresponding values.
In order to compare the carbon emission permits of each region horizontally, this paper calculated the average space of carbon emission in six regions from 2007–2016, and the results are displayed in Figure 5.
In Figure 5, the values of carbon emission permits in north and northeast China were negative from 2007–2016, and the phenomenon of excess emissions was most serious in 2011, 2012, and 2013. Although the situation has improved since 2014, it is still necessary to develop the production technology level and energy utilization to reduce CO 2 emissions. As for the northwest, carbon emission permit values continued to be negative after 2010. In addition, the carbon emission permits of the southwest were almost zero, which indicates that there was almost no difference between the actual consumption and carbon emissions. The space of carbon emission in southern and eastern areas was not only large, but also gradually increasing because of the economy and technology.
The indicator of carbon emission reduction potential was used to measure the emission reduction potential of a region, which is calculated as 1-(allocation/actually emission). If the value is below zero, i.e., the carbon allocation is greater than the actual consumption in the province, which indicates that there is no emission reduction potential, on the contrary, the province has a large potential of emission reduction. The carbon emission reduction potential of provinces in 2016 is shown in Figure 4.
As can be seen from Figure 4, the northeast and northwest regions, especially Shanxi, Inner Mongolia, Qinghai, Ningxia, and Xinjiang provinces, have huge potential for emission reduction, while there are few in the eastern and southern regions. In order to compare the carbon emission reduction potential of each region over time, we calculated the average potential of carbon emission reduction of the six regions from 2007–2016, and the results are shown in Figure 6 below.
In Figure 6, the potential of carbon emission reduction in the north, northeast, southwest, and east was basically stable from 2007–2016. Furthermore, the carbon emission reduction potential of the northeast and southwest remained less, while that of northern and eastern China was below zero. The potential in regions such as the south and the northwest changed greatly, and the trend was basically the same. In addition, the potential of carbon emission reduction in the northwest increased significantly and gradually became stable after 2012.
To investigate the relationship between the carbon emission permits and the initial efficiency value of provinces from 2007–2016, Figure 7 is drawn as follows.
The abscissa of Figure 7 represents the average space of carbon emission of each province from 2007–2016; the ordinate represents the average efficiency value of each province; and the color bar on the right represents the average potential of emission reduction. If the value of average emission reduction potential is below zero, then the color is dark, and the potential of emission reduction is less. As can be seen from Figure 7, 30 provinces were divided into four categories according to the average efficiency value (i.e., 0.7) and the space of carbon emission (i.e., zero). The provinces with an efficiency value greater than 0.7 and carbon emission permits greater than zero included Qinghai, Hainan, Beijing, Tianjin, Guangdong, Shanghai, and Jiangsu. The provinces with an efficiency value greater than 0.7 and the carbon emission permit less than zero included Ningxia and Inner Mongolia. The provinces with an efficiency value less than 0.7 and carbon emission space less than zero involved Shanxi, Hebei, Guizhou, Shandong, Yunnan, Gansu, Jilin, Xinjiang, Heilongjiang, Henan., and Shaanxi, while provinces with an efficiency value less than 0.7 and carbon emission permit greater than 0 were Liaoning, Anhui, Chongqing, Jiangxi, Hubei, Guangxi, Sichuan, Hunan, Fujian, and Zhejiang. Shanxi, Inner Mongolia, Ningxia, Guizhou, and other provinces have great potential for emission reduction. In addition, coastal areas had high efficiency values and carbon emission permits, but the efficiency value of the eastern, central, and southern regions were low, while northern regions had high efficiency values.
From Figure 8, Beijing and Hainan have no potential for emissions reduction, followed by Tianjin, Guangdong, Fujian, Shanghai, and other coastal regions. Additionally, some provinces such as Qinghai, Zhejiang, Guangxi, Jiangsu, Jiangxi, Hunan, Chongqing, Sichuan, and Hubei almost have no potential to reduce emissions, while Anhui, Liaoning, Jilin, Shandong, Henan, Yunnan, Shaanxi, Gansu, Hebei, and Heilongjiang have emissions reduction potential, and the potential of emission reduction in Xinjiang, Ningxia, Guizhou, Inner Mongolia, and Shanxi is the largest.

3.3. Policy Suggestions

Based on the above analysis, we draw the following conclusions. Firstly, there was no space for carbon emission in the north, northeast, and northwest from 2007–2016, particularly in the north and northwest where the carbon emission space was seriously insufficient. In southern regions, the space for carbon emission always existed during this period. In addition, the carbon emission permits of southern and eastern regions were increasing year by year. The actual carbon emission in southwestern regions was equivalent to the carbon emission allocation right, i.e., the carbon emission space barely fluctuated during this decade. Secondly, the potential of carbon emission reduction of each region tended to be stable after 2014, and in the north and northwest, it fluctuated greatly from 2007–2016. Besides, the northwest region had the potential of emission reduction since 2010, while this also existed in the north region after 2014.
Based on the above conclusions, this paper puts forward some policy suggestions of energy saving and emission reduction in provinces of China under the background of a low-carbon economy.
First, we need to define the scope of key regions and evaluate their emission reduction potential. During 2007–2016, the emission reduction potential of provinces such as Shanxi, Inner Mongolia, Guizhou, Heilongjiang, Hebei, Ningxia, Xinjiang, Gansu, Shaanxi, Yunnan, and Henan was higher. Shandong, Jilin, Liaoning, and Anhui had a certain potential to reduce emissions, and the emission reduction potential of other developed areas such as Guangdong, Beijing, and Tianjin was very low. It is worth noting that Hainan and Qinghai had significantly increased their emission reduction potential in recent years, as well as Shanxi, Ningxia, Xinjiang, and other places have great potential for emission reduction. In order to achieve the goal of emission reduction, we should focus on these areas.
Second, we should adjust and optimize the industrial structure and promote the green development of industry. In terms of the northeast, we should focus on these key areas and optimize industrial layout. At the same time, the prohibited or restricted industries, production technology, and industrial catalogs should be clearly defined. We should revise and improve the entry conditions for high energy consumption, high pollution, and resource-oriented industries. In addition, we should accelerate adjusting the distribution of industries in different regions and relocate, transform, or upgrade the enterprises with much pollution. Besides, we should strictly control the industrial production capacity with large energy consumption and high pollution, thoroughly deal with the industrial pollution, and strengthen comprehensive treatment of enterprises with a large amount of pollution. In the end, we should expand the scale of green industries such as industries of energy conservation, environmental protection, clean production, and clean energy, and new energy sources also need to be developed.
Finally, we should accelerate energy restructuring and build a low-carbon energy system. Especially, we should pay attention to the northeast of China as the focus. The specific measures are as follows: We should optimize the energy structure and develop new energy sources to reduce the ratio of coal and other fossil energy sources in total consumption. Moreover, we should improve energy efficiency and accelerate the development of clean and new energy sources. As for some regions with more private cars such as the eastern and southern, we should promote the use of vehicles with new energy, for example the number of buses, taxis, private cars, and other vehicles with new energy and clean energy should be increased and updated. In addition, we should take measures such as economic compensation to restrict the use of old vehicles and supervise the excess emissions in key areas. With respect to the northwestern region, we should optimize the energy structure and accelerate the upgrading of rural areas from coal to electricity and from coal to natural gas. Besides, some new energy such as nuclear energy, wind energy, solar energy, and other new energy sources should be developed to replace fossil fuels in a safe and effective manner.

4. Conclusions

Under the background of advocating energy saving and emission reduction, how to allocate carbon emission rights has become a hot research topic among global scholars. Especially in China, it is urgent to formulate a reasonable and effective scheme of carbon emission right allocation. In recent years, many researchers have applied DEA to the fixed cost allocation problem, among which carbon emission right allocation has become one of the most important applications, but all of this literature using traditional models considers each DMUs as a “black box” without taking the internal structure into account. However, it is more accurate for formulating the scheme when considering the inner operation of DMUs. This paper firstly adopts the two-stage DEA-network model to investigate the allocation plan of carbon emission right, and the algorithm in solving the unique allocation scheme is modified to make the result more accurate. Based on the result analysis, corresponding policy suggestions are given. First, we need to define the scope of key regions and evaluate their emission reduction potential. Second, we should adjust and optimize the industrial structure and promote the green development of industry. Finally, we should accelerate energy restructuring and build a low-carbon energy system.
Although it is reasonable for this paper to firstly apply the two-stage DEA network model to the research of carbon emission right allocation in China and the model established can be extended to other countries in the world, it has some limitations at the same time. On the one hand, for the treatment of undesirable output, we regard carbon emission as input; although policymakers hope the undesirable output to be as small as the input, this is actually inconsistent with the actual production process. On the other hand, this article does not consider the problems of the sharing of resources in the two stages and the addition of new inputs in the second stage, while the process of production is more realistic. Therefore, the above two points are our future research directions. Besides, the model adopted in this paper takes the size of the operational units into account to formulate a unique allocation plan, so it is also an improvement point to find a better method for dealing with the problem that allocation schemes are not unique.

Author Contributions

Conceptualization, X.J. and X.T.; methodology, X.J., B.Z., and X.T.; software, X.J.; validation, X.T. and K.R.; formal analysis, X.J.; investigation, X.J.; resources, X.J.; data curation, X.J.; writing, original draft preparation, X.J.; writing, review and editing, X.T. and C.W.; visualization, X.J.; supervision, B.Z. and X.T.; project administration, X.T.; funding acquisition, X.T.

Funding

This research was funded by the Natural Science Foundation of China Grant No. 61772011, the Open Project Foundation of Intelligent Information Processing Key Laboratory of Shanxi Province Grant No. CICIP2018002, the Social Science Foundation of Beijing Grant No. 16JDGLC005, the Project of Great Wall Scholar, Beijing Municipal Commission of Education, Grant No. CIT & TCD20170304 and No. CIT & TCD20180305, the Frontier Science Key Program of the Chinese Academy of Sciences Grant No. QYZDY-SSW-DQC004, the Beijing Municipal Science and Technology Project Grant No. Z171100004417025, and the Science and Technology Planning Project of Guangdong Province, China, Grant No. 2016B02024007.

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A

Model (5) is a fractional model, so we need to transform it into a linear programming model by C-C transformation. Let τ = 1 i = 1 m υ i x i d + υ m + 1 R 1 d + p = 1 q φ p z p d + υ m + 1 R 2 d , u r = τ · μ r , ϕ p = τ · φ p , v i = τ · υ i , u 0 = τ · μ 0 , ϕ 0 = τ · φ 0 , v m + 1 = τ · υ m + 1 . Then, we have:
θ d = M a x ( p = 1 q ϕ p z p d + ϕ 0 + r = 1 s u r y r d + u 0 ) s . t . p = 1 q ϕ p z p j + ϕ 0 i = 1 m v i x i j v m + 1 R 1 j 0 , j = 1 , 2 , , n r = 1 s u r y r j + u 0 p = 1 q ϕ p z p j v m + 1 R 2 j 0 , j = 1 , 2 , , n i = 1 m v i x i d + v m + 1 R 1 d + p = 1 q ϕ p z p d + v m + 1 R 2 d = 1 j = 1 n ( R 1 j + R 2 j ) = R , R 1 j , R 2 j 0 , j = 1 , 2 , , n u r , ϕ p , v i 0 , v m + 1 > 0 , r = 1 , , s ; p = 1 , , q ; i = 1 , , m
Model (A1) is still nonlinear due to v m + 1 R 1 j and v m + 1 R 2 j ; in order to transform the model (A1) into a linear model, let v m + 1 R 1 j = r 1 j , v m + 1 R 2 j = r 2 j , then Model (A1) can be rewritten as Model (6).

Appendix B

Before calculating the operation size parameters, we use formula (A2) to normalize the original data.
x ^ i j = x i j j = 1 n x i j , z ^ p j = z p j j = 1 n z p j , y ^ r j = y r j j = 1 n y r j , i , p , r , j
Suppose the size parameters of the first stage and second stage are α j and β j , respectively. The size parameters are as follows:
α j = p = 1 q z ^ p j · i = 1 m x ^ i j j = 1 n ( p = 1 q z ^ p j · i = 1 m x ^ i j + r = 1 s y ^ r j · p = 1 q z ^ p j ) , β j = r = 1 s y ^ r j · p = 1 q z ^ p j j = 1 n ( p = 1 q z ^ p j · i = 1 m x ^ i j + r = 1 s y ^ r j · p = 1 q z ^ p j )
where j = 1 n ( α j + β j ) = 1 .

Appendix C

Since the model (8) is nonlinear, we need to do the following processing to change it into a linear model.
Let α j R R 1 j + α j R R 1 j = 2 c 1 j , α j R R 1 j α j R + R 1 j = 2 c 2 j , β j R R 2 j + β j R R 2 j = 2 d 1 j , β j R R 2 j β j R + R 2 j = 2 d 2 j , then the model (8) is transformed to Model (A4).
M i n M a x ( c 1 k + c 2 k + d 1 k + d 2 k ) s . t . α j R R 1 j = c 1 j + c 2 j j = 1 , 2 , , n β j R R 2 j = d 1 j d 2 j , j = 1 , 2 , , n p = 1 q ϕ p z p j i = 1 m v i x i j R 1 j + ϕ 0 = 0 , j = 1 , 2 , , n r = 1 s u r y r j p = 1 q ϕ p z p j R 2 j + u 0 , j = 1 , 2 , , n j = 1 n ( R 1 j + R 2 j ) = R , R 1 j , R 2 j 0 , j = 1 , 2 , , n c 1 j , c 2 j , d 1 j , d 2 j 0 , j = 1 , 2 , , n u r , ϕ p , v i 0 , r = 1 , 2 , , s ; p = 1 , 2 , , q ; i = 1 , 2 , , m
Although the model (A4) is linear, it cannot be directly solved since this model involves multiple objectives. Therefore, we let M a x ( c 1 k + c 2 k + d 1 k + d 2 k ) = ρ , then the model (A4) can be rewritten as the model (A5).
M i n M a x ρ s . t . c 1 j + c 2 j + d 1 j + d 2 j ρ j = 0 , j = 1 , 2 , , n ρ j ρ 0 , j = 1 , 2 , , n α j R R 1 j = c 1 j + c 2 j j = 1 , 2 , , n β j R R 2 j = d 1 j d 2 j , j = 1 , 2 , , n p = 1 q ϕ p z p j i = 1 m v i x i j R 1 j + ϕ 0 = 0 , j = 1 , 2 , , n r = 1 s u r y r j p = 1 q ϕ p z p j R 2 j + u 0 , j = 1 , 2 , , n j = 1 n ( R 1 j + R 2 j ) = R , R 1 j , R 2 j 0 , j = 1 , 2 , , n c 1 j , c 2 j , d 1 j , d 2 j 0 , j = 1 , 2 , , n u r , ϕ p , v i 0 , r = 1 , 2 , , s ; p = 1 , 2 , , q ; i = 1 , 2 , , m

References

  1. Zhang, K.; Chen, H.H. Environmental Performance and Financing Decisions Impact on Sustainable Financial Development of Chinese Environmental Protection Enterprises. Sustainability 2017, 9, 2260. [Google Scholar] [CrossRef]
  2. Shi, X. Inequality of opportunity in energy consumption in China. Energy Policy 2019, 124, 371–382. [Google Scholar] [CrossRef]
  3. Bing, G.; Xiaochen, Z.; Qing, G.; Ordieres-Meré, J. Discovering the patterns of energy consumption, GDP, and CO2 emissions in China using the cluster method. Energy 2019, 166, 1149–1167. [Google Scholar]
  4. Yi, W.J.; Zou, L.L.; Guo, J.; Wang, K.; Wei, Y.M. How can China reach its CO intensity reduction targets by 2020? A regional allocation based on equity and development. Energy Policy 2011, 39, 2407–2415. [Google Scholar] [CrossRef]
  5. Li, A.; Lin, B. Comparing climate policies to reduce carbon emissions in China. Energy Policy 2013, 60, 667–674. [Google Scholar] [CrossRef]
  6. Dong, F.; Long, R.; Yu, B.; Wang, Y.; Li, J.; Wang, Y.; Dai, Y.; Yang, Q.; Chen, H. How can China allocate CO2 reduction targets at the provincial level considering both equity and efficiency? Evidence from its Copenhagen Accord pledge. Resour. Conserv. Recycl. 2018, 130, 31–43. [Google Scholar] [CrossRef]
  7. Zhou, P.; Wang, M. Carbon dioxide emissions allocation: A review. Ecol. Econ. 2016, 125, 47–59. [Google Scholar] [CrossRef]
  8. Cramton, P.; Kerr, S. Tradeable carbon permit auctions: How and why to auction not grandfather. Energy Policy 2002, 30, 333–345. [Google Scholar] [CrossRef]
  9. Kampas, A.; White, B. Selecting permit allocation rules for agricultural pollution control: A bargaining solution. Ecol. Econ. 2003, 47, 135–147. [Google Scholar] [CrossRef]
  10. Daskalakis, G.; Psychoyios, D.; Markellos, R.N. Modeling CO2 emission allowance prices and derivatives: Evidence from the European trading scheme. J. Bank. Financ. 2009, 33, 1230–1241. [Google Scholar] [CrossRef]
  11. Oestreich, A.M.; Tsiakas, I. Carbon emissions and stock returns: Evidence from the EU Emissions Trading Scheme. J. Bank. Financ. 2015, 58, 294–308. [Google Scholar] [CrossRef]
  12. Martin, R.; Muûls, M.; Wagner, U.J. The impact of the European Union Emissions Trading Scheme on regulated firms: What is the evidence after ten years? Rev. Environ. Econ. Policy 2015, 10, 129–148. [Google Scholar] [CrossRef]
  13. Daskalakis, G. Temporal restrictions on emissions trading and the implications for the carbon futures market: Lessons from the EU emissions trading scheme. Energy Policy 2018, 115, 88–91. [Google Scholar] [CrossRef]
  14. Mardani, A.; Zavadskas, E.K.; Streimikiene, D.; Jusoh, A.; Khoshnoudi, M. A comprehensive review of data envelopment analysis (DEA) approach in energy efficiency. Renew. Sustain. Energy Rev. 2017, 70, 1298–1322. [Google Scholar] [CrossRef]
  15. Li, Y.; Li, F.; Emrouznejad, A.; Liang, L.; Xie, Q. Allocating the fixed cost: an approach based on data envelopment analysis and cooperative game. Ann. Oper. Res. 2019, 274, 373–394. [Google Scholar] [CrossRef]
  16. Xie, Q.; Li, Y.; Wang, L.; Chao, L. Improving discrimination in data envelopment analysis without losing information based on Renyi’s entropy. Cent. Eur. J. Oper. Res. 2018, 26, 1053–1068. [Google Scholar] [CrossRef]
  17. Feng, L.; Zhu, Q.; Liang, L. A new data envelopment analysis based approach for fixed cost allocation. Ann. Oper. Res. 2019, 274, 347–372. [Google Scholar]
  18. Sueyoshi, T.; Yuan, Y.; Goto, M. A literature study for DEA applied to energy and environment. Energy Econ. 2017, 62, 104–124. [Google Scholar] [CrossRef]
  19. Xie, Q.; Chen, X.; Li, L.; Rao, K.; Tao, L.; Ma, C. Image Fusion Based on Kernel Estimation and Data Envelopment Analysis. Int. J. Inf. Technol. Decis. Mak. 2019, 17. [Google Scholar] [CrossRef]
  20. Li, Y.; Yang, F.; Liang, L.; Hua, Z. Allocating the fixed cost as a complement of other cost inputs: A DEA approach. Eur. J. Oper. Res. 2009, 197, 389–401. [Google Scholar] [CrossRef]
  21. Zhou, P.; Zhang, L.; Zhou, D.; Xia, W. Modeling economic performance of interprovincial CO2 emission reduction quota trading in China. Appl. Energy 2013, 112, 1518–1528. [Google Scholar] [CrossRef]
  22. Yu, S.; Wei, Y.M.; Wang, K. Provincial allocation of carbon emission reduction targets in China: An approach based on improved fuzzy cluster and Shapley value decomposition. Energy Policy 2014, 66, 630–644. [Google Scholar] [CrossRef]
  23. Miao, Z.; Geng, Y.; Sheng, J. Efficient allocation of CO2 emissions in China: A zero sum gains data envelopment model. J. Clean. Prod. 2016, 112, 4144–4150. [Google Scholar] [CrossRef]
  24. Cook, W.D.; Kress, M. Characterizing an equitable allocation of shared costs: A DEA approach1. Eur. J. Oper. Res. 1999, 119, 652–661. [Google Scholar] [CrossRef]
  25. Beasley, J. Allocating fixed costs and resources via data envelopment analysis. Eur. J. Oper. Res. 2003, 147, 198–216. [Google Scholar] [CrossRef] [Green Version]
  26. Nakabayashi, K.; Tone, K. Egoist’s dilemma: A DEA game. Omega 2006, 34, 135–148. [Google Scholar] [CrossRef]
  27. Li, F.; Zhu, Q.; Liang, L. Allocating a fixed cost based on a DEA-game cross efficiency approach. Expert Syst. Appl. 2018, 96, 196–207. [Google Scholar] [CrossRef]
  28. Gomes, E.; Lins, M.E. Modelling undesirable outputs with zero sum gains data envelopment analysis models. J. Oper. Res. Soc. 2008, 59, 616–623. [Google Scholar] [CrossRef]
  29. Lozano, S.; Villa, G.; Brännlund, R. Centralised reallocation of emission permits using DEA. Eur. J. Oper. Res. 2009, 193, 752–760. [Google Scholar] [CrossRef]
  30. Wei, C.; Ni, J.; Du, L. Regional allocation of carbon dioxide abatement in China. China Econ. Rev. 2012, 23, 552–565. [Google Scholar] [CrossRef]
  31. Wang, K.; Zhang, X.; Wei, Y.M.; Yu, S. Regional allocation of CO2 emissions allowance over provinces in China by 2020. Energy Policy 2013, 54, 214–229. [Google Scholar] [CrossRef]
  32. Feng, C.; Chu, F.; Ding, J.; Bi, G.; Liang, L. Carbon emissions abatement (CEA) allocation and compensation schemes based on DEA. Omega 2015, 53, 78–89. [Google Scholar] [CrossRef]
  33. Zhou, Z.; Liu, C.; Zeng, X.; Jiang, Y.; Liu, W. Carbon emission performance evaluation and allocation in Chinese cities. J. Clean. Prod. 2018, 172, 1254–1272. [Google Scholar] [CrossRef]
  34. Yu, A.; You, J.; Rudkin, S.; Zhang, H. Industrial carbon abatement allocations and regional collaboration: Re-evaluating China through a modified data envelopment analysis. Appl. Energy 2019, 233, 232–243. [Google Scholar] [CrossRef]
  35. Färe, R. Measuring Farrell efficiency for a firm with intermediate inputs. Acad. Econ. Papers 1991, 19, 329–340. [Google Scholar]
  36. Wanke, P.; Barros, C. Two-stage DEA: An application to major Brazilian banks. Expert Syst. Appl. 2014, 41, 2337–2344. [Google Scholar] [CrossRef]
  37. Li, K.; Lin, B. Impact of energy conservation policies on the green productivity in China’s manufacturing sector: Evidence from a three-stage DEA model. Appl. Energy 2016, 168, 351–363. [Google Scholar] [CrossRef]
  38. Yu, M.M.; Chen, L.H.; Hsiao, B. A fixed cost allocation based on the two-stage network data envelopment approach. J. Bus. Res. 2016, 69, 1817–1822. [Google Scholar] [CrossRef]
  39. Zhu, W.; Zhang, Q.; Wang, H. Fixed costs and shared resources allocation in two-stage network DEA. Ann. Oper. Res. 2017, 1–18. [Google Scholar] [CrossRef]
  40. Ding, T.; Zhu, Q.; Zhang, B.; Liang, L. Centralized fixed cost allocation for generalized two-stage network DEA. INFOR Inf. Syst. Oper. Res. 2017, 1–18. [Google Scholar] [CrossRef]
  41. Li, F.; Zhu, Q.; Chen, Z. Allocating a fixed cost across the decision making units with two-stage network structures. Omega 2019, 83, 139–154. [Google Scholar] [CrossRef]
  42. Charnes, A.; Cooper, W.W.; Rhodes, E. Measuring the efficiency of decision making units. Eur. J. Oper. Res. 1978, 2, 429–444. [Google Scholar] [CrossRef]
  43. Banker, R.D.; Charnes, A.; Cooper, W.W. Some models for estimating technical and scale inefficiencies in data envelopment analysis. Manag. Sci. 1984, 30, 1078–1092. [Google Scholar] [CrossRef]
  44. Zhang, N.; Wei, X. Dynamic total factor carbon emissions performance changes in the Chinese transportation industry. Appl. Energy 2015, 146, 409–420. [Google Scholar] [CrossRef]
  45. Kong, Y.; Zhao, T.; Yuan, R.; Chen, C. Allocation of carbon emission quotas in Chinese provinces based on equality and efficiency principles. J. Clean. Prod. 2019, 211, 222–232. [Google Scholar] [CrossRef]
  46. Guajardo, M.; Jörnsten, K. Common mistakes in computing the nucleolus. Eur. J. Oper. Res. 2015, 241, 931–935. [Google Scholar] [CrossRef] [Green Version]
  47. Jiang, X.T.; Wang, Q.; Li, R. Investigating factors affecting carbon emission in China and the USA: A perspective of stratified heterogeneity. J. Clean. Prod. 2018, 199, 85–92. [Google Scholar] [CrossRef]
  48. Wang, Q.; Zeng, Y.E.; Wu, B.W. Exploring the relationship between urbanization, energy consumption, and CO2 emissions in different provinces of China. Renew. Sustain. Energy Rev. 2016, 54, 1563–1579. [Google Scholar] [CrossRef]
  49. Ma, X.; Wang, C.; Dong, B.; Gu, G.; Chen, R.; Li, Y.; Zou, H.; Zhang, W.; Li, Q. Carbon emissions from energy consumption in China: Its measurement and driving factors. Sci. Total Environ. 2019, 648, 1411–1420. [Google Scholar] [CrossRef] [PubMed]
  50. Chung, Y.H.; Färe, R.; Grosskopf, S. Productivity and undesirable outputs: A directional distance function approach. J. Environ. Manag. 1997, 51, 229–240. [Google Scholar] [CrossRef]
  51. Färe, R.; Grosskopf, S.; Lovell, C.K.; Pasurka, C. Multilateral productivity comparisons when some outputs are undesirable: A nonparametric approach. Rev. Econ. Stat. 1989, 71, 90–98. [Google Scholar] [CrossRef]
  52. Seiford, L.M.; Zhu, J. Modeling undesirable factors in efficiency evaluation. Eur. J. Oper. Res. 2002, 142, 16–20. [Google Scholar] [CrossRef]
  53. Hailu, A.; Veeman, T.S. Non-parametric productivity analysis with undesirable outputs: An application to the Canadian pulp and paper industry. Am. J. Agric. Econ. 2001, 83, 605–616. [Google Scholar] [CrossRef]
  54. Wang, S.; Ma, Y. Influencing factors and regional discrepancies of the efficiency of carbon dioxide emissions in Jiangsu, China. Ecol. Indic. 2018, 90, 460–468. [Google Scholar] [CrossRef]
  55. Jun, Z.; Guiying, W.; Jipeng, Z. The Estimation of China’s provincial capital stock: 1952–2000. Econ. Res. J. 2004, 10, 35–44. [Google Scholar]
  56. Chang, N. Sharing responsibility for carbon dioxide emissions: A perspective on border tax adjustments. Energy Policy 2013, 59, 850–856. [Google Scholar] [CrossRef]
  57. Chung, W.S.; Tohno, S.; Shim, S.Y. An estimation of energy and GHG emission intensity caused by energy consumption in Korea: An energy IO approach. Appl. Energy 2009, 86, 1902–1914. [Google Scholar] [CrossRef]
  58. Tang, X.; Yongjun, L.I.; Wang, M.; Liang, L.; Management, S.O. A DEA efficiency probability dominance method. Syst. Eng. Theory Pract. 2016, 36, 2641–2647. [Google Scholar]
Figure 1. Two-stage network process of D M U j .
Figure 1. Two-stage network process of D M U j .
Sustainability 11 01369 g001
Figure 2. Two-stage network process for carbon emission right allocation.
Figure 2. Two-stage network process for carbon emission right allocation.
Sustainability 11 01369 g002
Figure 3. The differences between each region in 2016.
Figure 3. The differences between each region in 2016.
Sustainability 11 01369 g003
Figure 4. Carbon emission permits and emission reduction potential of provinces in 2016.
Figure 4. Carbon emission permits and emission reduction potential of provinces in 2016.
Sustainability 11 01369 g004
Figure 5. Carbon emission space of each region.
Figure 5. Carbon emission space of each region.
Sustainability 11 01369 g005
Figure 6. Carbon emission reduction potential of each region.
Figure 6. Carbon emission reduction potential of each region.
Sustainability 11 01369 g006
Figure 7. The average results of provinces.
Figure 7. The average results of provinces.
Sustainability 11 01369 g007
Figure 8. The average potential of emission reduction.
Figure 8. The average potential of emission reduction.
Sustainability 11 01369 g008
Table 1. Carbon emission coefficients.
Table 1. Carbon emission coefficients.
Energy TypesThe Net Calorific ValueCarbon Emission Factors
UnitValue
coalkJ/kg20,90825.8
gasolinekJ/kg43,07018.9
kerosenekJ/kg43,07019.5
dieselkJ/kg42,65220.2
fuel oilkJ/kg41,81621.1
natural gaskJ/m 3 38,93115.3
Table 2. Descriptive statistics of the input and output levels for 30 provinces in 2016.
Table 2. Descriptive statistics of the input and output levels for 30 provinces in 2016.
Descriptive StatisticsAverageSDMaximumMinimum
Capital stock91,033.1656,517.92237,203.2114,362.83
Labor force595.22423.181957.5763.09
Energy consumption15,665.7311,388.6847,772.591393.82
Carbon emission32,617.8823,409.29102,669.793139.03
Urbanization rate0.590.110.880.44
Private car ownership543.30401.751550.6572.96
GDP25,963.9519,937.9980,854.912572.49
Table 3. The efficiency values of the provinces before the allocation.
Table 3. The efficiency values of the provinces before the allocation.
ProvinceBCC.Eff a BCC.RankStage1.Eff b Stage2.Eff c Overall.Eff d Overall.Rank
Beijing1110.76920.90544
Tianjin110.81430.76570.79256
Hebei0.95572210.30160.650821
Shanxi0.7910290.89320.34780.635923
Inner Mongolia110.76130.51360.654220
Liaoning0.7551300.66010.43580.570930
Jilin0.9067260.72020.55130.649522
Heilongjiang0.8761270.73740.54740.656819
Shanghai110.738010.84935
Jiangsu110.67860.65430.668817
Zhejiang0.9399240.98350.44870.71838
Anhui0.9463230.68250.48400.602027
Fujian110.72580.61530.679314
Jiangxi0.9705200.67460.57530.634624
Shandong110.80430.47480.657418
Henan110.71190.43120.595228
Hubei110.63060.60680.621426
Hunan110.77380.56380.682213
Guangdong110.97440.58530.78247
Guangxi110.82630.53720.695510
Hainan1110.971811
Chongqing0.9596210.67630.66710.672615
Sichuan110.99370.39620.69599
Guizhou110.80060.53870.684112
Yunnan1110.38710.693611
Shaanxi0.8729280.62970.50720.582329
Gansu110.70450.62000.669616
Qinghai111111
Ningxia1110.89540.95103
Xinjiang0.9394250.70670.53020.633625
a BCC.Eff: BCC Efficiency; b Stage 1.Eff: Stage 1 Efficiency; c Stage2.Eff: Stage 2 Efficiency; d Overall.Eff: Overall.Efficiency.
Table 4. Carbon emission right allocation in 2016.
Table 4. Carbon emission right allocation in 2016.
ProvinceActual Carbon EmissionsStage 1 AllocationStage 2 AllocationOverall AllocationCarbon Emission Space
Beijing5825.549220.1711,844.0721,064.2415,238.69
Tianjin10,918.0810,014.6410,111.0220,125.669207.58
Hebei60,150.4047,892.501088.8148,981.31−11,169.09
Shanxi73,066.1623,885.3479.1623,964.50−49,101.66
Inner Mongolia75,158.7125,141.857356.6432,498.49−42,660.23
Liaoning40,263.7726,521.466550.2233,071.68−7192.09
Jilin20,477.6613,092.275993.1519,085.42−1392.23
Heilongjiang30,376.1813,555.965683.9419,239.90−11,136.28
Shanghai16,686.646331.6919,833.4226,165.109478.46
Jiangsu62,239.0354,698.2042,451.3197,149.5134,910.48
Zhejiang34,570.8938,748.0516,548.6755,296.7120,725.82
Anhui34,802.0724,804.3710,418.3235,222.69420.61
Fujian17,400.9516,955.2516,313.3733,268.6215,867.68
Jiangxi17,773.4114,345.578871.9223,217.495444.08
Shandong102,669.7971,214.8025,463.8796,678.67−5991.12
Henan51,235.0441,066.0213,592.9354,658.963423.92
Hubei28,889.4523,236.9218,089.6841,326.6012,437.15
Hunan27,141.8123,832.3916,556.7840,389.1713,247.37
Guangdong44,437.5544,489.4439,439.0883,928.5239,490.97
Guangxi15,978.4914,217.728184.2522,401.966423.48
Hainan3139.03206.161568.581774.75−1364.29
Chongqing13,991.0210,921.459613.7320,535.176544.15
Sichuan24,729.4329,454.7711,626.9741,081.7416,352.31
Guizhou29,718.1911,652.603631.3815,283.99−14,434.20
Yunnan18,010.1016,016.661576.5417,593.20−416.90
Shaanxi41,315.7320,707.807241.7327,949.53−13,366.19
Gansu14,272.066183.351847.258030.61−6241.46
Qinghai4557.2761.23584.19645.42−3911.84
Ningxia17,931.763150.05165.033315.07−14,616.68
Xinjiang40,810.3112,048.272543.5914,591.85−26,218.45

Share and Cite

MDPI and ACS Style

Jin, X.; Zou, B.; Wang, C.; Rao, K.; Tang, X. Carbon Emission Allocation in a Chinese Province-Level Region Based on Two-Stage Network Structures. Sustainability 2019, 11, 1369. https://doi.org/10.3390/su11051369

AMA Style

Jin X, Zou B, Wang C, Rao K, Tang X. Carbon Emission Allocation in a Chinese Province-Level Region Based on Two-Stage Network Structures. Sustainability. 2019; 11(5):1369. https://doi.org/10.3390/su11051369

Chicago/Turabian Style

Jin, Xi, Bin Zou, Chan Wang, Kaifeng Rao, and Xiaowen Tang. 2019. "Carbon Emission Allocation in a Chinese Province-Level Region Based on Two-Stage Network Structures" Sustainability 11, no. 5: 1369. https://doi.org/10.3390/su11051369

APA Style

Jin, X., Zou, B., Wang, C., Rao, K., & Tang, X. (2019). Carbon Emission Allocation in a Chinese Province-Level Region Based on Two-Stage Network Structures. Sustainability, 11(5), 1369. https://doi.org/10.3390/su11051369

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