1. Introduction
With climate change becoming an increasingly serious issue, the reduction of carbon dioxide (CO
2) emissions has attracted extensive attention worldwide. As the greatest CO
2 emitter in the world [
1,
2], China has shown its determination for developing a low-carbon economy and promised to abate its CO
2 emissions per unit of gross domestic product (GDP) (i.e., carbon intensity) by 40–45% by 2020 compared with that in 2005 [
3]. Further, China set the latest target of abating carbon intensity by 18% by 2020, with 2015 as the reference year [
4]. In order to achieve the above international commitments for mitigating CO
2 emissions, China’s National Development and Reform Commission (NDRC) has launched seven pilot emission trading schemes (ETS) since 2013 [
5], which are specifically located in Shenzhen, Guangdong, Shanghai, Beijing, Tianjin, Chongqing, and Hubei. These regional carbon markets are considered as experimental explorations for the establishment of China’s nationwide emission trading scheme (CN-ETS), which is scheduled to be launched in 2017. It is reported the CN-ETS will cover seven emission-intensive industries, including paper making, electricity generation, metallurgy, non-ferrous metals, building materials, the chemical industry, and the aviation service industry [
6].
In the context of achieving the construction and operation of CN-ETS, it is of great urgency and necessity to obtain a good understanding of the participating sectors in terms of energy utilization and CO
2 emissions [
7]. In this regard, estimating the CO
2 emission performance (CEP), mitigation potential (MP), and marginal abatement costs (MAC) for these sectors can provide valuable information for the governments and participating enterprises. From the perspective of the government, a good knowledge of CEP, MP, and MAC could help design appropriate market mechanisms for the CN-ETS, e.g., the estimated MAC may be used as a reference for carbon pricing [
8,
9]. On the other hand, a comprehensive acquaintance of MAC among the participating sectors could help the participating enterprises to determine the best mitigation strategies [
10]. Moreover, to the best knowledge of the authors, there have been few studies on the CEP, MP, and MAC of the sectors covered in the CN-ETS, and this paper aims to fill this research gap.
Additionally, it is accepted that there exists significant heterogeneity in terms of the production technology among various sectors [
11], which is regarded as an obstacle to the objective evaluation of CEP, MP, and MAC [
12]. Therefore, taking the technology heterogeneity into consideration, we employ a joint framework consisting of the directional distance function (DDF) and meta-frontier analysis to estimate CEP, MP, and MAC under the meta-frontier and the group-frontier, respectively. Following this, we investigate the sectoral distributions of CEP, MP, and MAC under both frontiers, and analyze the differences between the two frontiers at sector levels. Finally, several policy recommendations for the CN-ETS are made based on the conclusions. The rest of the paper proceeds as follows:
Section 2 provides a literature review.
Section 3 introduces the methods and materials.
Section 4 presents the empirical results and discussion.
Section 5 draws conclusions and policy implications.
2. Literature Review
DDF, theorized and developed by Chung et al. [
13] and Chambers et al. [
14,
15], has been widely employed to study energy and environmental issues [
16,
17,
18,
19,
20]. The main advantage of the DDF is that it can achieve the expansion of desirable outputs and reduction of undesirable outputs (e.g., CO
2), simultaneously [
18]. Generally, two estimation techniques are often employed to estimate the DDF, namely the parametric and non-parametric methods. Compared with the former, the non-parametric technique does not need to pre-determine any functional and parametric forms, thereby avoiding the impacts of subjective factors on the results [
17]. In this regard, the data envelopment analysis (DEA), a well-developed nonparametric frontier tool, is often combined with DDF to evaluate CEP [
21,
22,
23], MP [
9,
24], and MAC [
25,
26,
27]. For instance, Watanabe and Tanaka [
21] employed the DEA-DDF to estimate the environmental performance for China’s industrial sector at province levels from 1994 to 2002. Liu et al. [
26] evaluated the carbon emission performance and marginal abatement cost for provinces in China by using a non-parametric DDF. Additionally, the same method was applied by Wei et al. [
9] to measure the reduction potential of CO
2 emissions for the thermal power plants in China’s Zhejiang province.
Considering the technology heterogeneity across decision making units (DMUs), Battese et al. [
28] and O’Donnell et al. [
29] incorporated the meta-frontier approach into DDF to formulate a joint framework. In the framework, DMUs with different production technologies are classified into several groups in which DMUs are deemed to be homogeneous, and then evaluated under the meta-frontier and the group-frontier, respectively. Recently, this combined methodology has been widely applied in the energy and environmental field [
30,
31,
32,
33,
34,
35,
36,
37,
38,
39]. For instance, using the combined method of DDF and the meta-frontier approach, Lin et al. [
33] evaluated the environmental performance of 63 countries during the period from 1981 to 2005. Further, Zhang et al. [
34] proposed a meta-frontier non-radial DDF by combining the meta-frontier approach with the non-radial DDF, and used it to assess the energy and CO
2 emission performance of electricity generation in Korea. Furthermore, the model was applied by Yao et al. [
36] to estimate China’s energy efficiency, carbon emission performance, and mitigation potential at regional levels. Additionally, based on the same model, Li and Song [
38] constructed a green development growth index to assess China’s green development at province levels.
The literature on energy and environmental issues is abundant at industry levels [
7,
10,
27,
40,
41]. For instance, Lee and Zhang [
27] measured the reduction potential and marginal abatement cost of CO
2 emissions for 30 of China’s manufacturing industries. Yuan et al. [
40] estimated the shadow prices of CO
2 emissions for China’s industrial sectors with the use of non-parametric DDF. Teng et al. [
41] employed multiple methods to derive the marginal abatement cost curves for China’s energy-intensive industries. Zhou et al. [
10] applied multiple distance function approaches to approximate the shadow prices of CO
2 emissions for Shanghai’s industrial sectors. Xiao et al. [
7] estimated the marginal abatement costs of CO
2 emissions for China’s industrial sectors during 2005–2011 by using a parametric DDF. Notably, it is found that the above studies mainly focus on the industrial sectors, while a comprehensive investigation on the participating sectors of CN-ETS has not been conducted. Furthermore, to the best of our knowledge, few studies take the industry heterogeneity into account, in addition to Xie et al. [
11] and Chung and Heshmati [
12]. In this context, we attempt to perform an empirical study on the participating sectors of CN-ETS in terms of CEP, MP, and MAC, taking into consideration the industry heterogeneity.
4. Results and Discussions
In this section, we first report the estimates of DDF values (), CO2 emission performance (CEP), mitigation potential (MP), and marginal abatement cost (MAC) under the meta-frontier and group-frontier at group levels. Then, we investigate the sectoral distributions of CEP, MP, and MAC under both frontiers, and analyze the differences between the two frontiers in terms of the CEP, MP, and MAC at sector levels.
4.1. Statistical Summary of Estimates
4.1.1. DDF Values
DDF values under the meta-frontier (
) and the group-frontier (
) are estimated by combining Equation (3) with Equations (6) and (7), respectively, and the results are listed in
Table 5. From the table, it can be observed that
and
vary from 0 to 0.888 and from 0 to 0.687, respectively. Moreover, the mean values of
are generally larger than those of
. From the perspective of total sectors, the average
is found to be 0.544, which is more than twice that of
(0.225), indicating that there is a significant difference. Additionally, the differences in the average DDF value between the two frontiers are found to be considerable at group levels, especially for the high-carbon and medium-carbon groups. Specifically, the average
(0.688) of the high-carbon group is found to be almost three times that of
(0.234). Further, the average DDF value of the medium-carbon group varies between the meta-frontier and group-frontier at greater levels, which is found to jump from 0.169 to 0.672. As for the low-carbon group, the difference between the two frontiers could be neglected. Notably, there are significant differences among the DDF values of the three groups when the meta-frontier is used as the basis of evaluation, while relatively moderate distinctions are observed when using the group frontiers.
4.1.2. CO2 Emission Performance
Table 6 reports the estimates of CO
2 emission performance under the meta-frontier (MCEP) and the group-frontier (GCEP). From the table, it is found that MCEP and GCEP vary from 0.112 to 1 and from 0.313 to 1, respectively. Moreover, the mean values of GCEP are generally higher than those of MCEP. From the perspective of total sectors, the average GCEP is found to be 0.775, while that under the meta-frontier is only 0.456, which indicates that there is a significant difference. Additionally, the differences in the average CEP scores between the two frontiers are found to be considerable at group levels, especially for the high-carbon and medium-carbon groups. Specifically, the average GCEP (0.766) of the high-carbon group is found to be more than twice that of MCEP (0.312). Similarly, the average CEP of the medium-carbon group is found to jump from 0.328 to 0.831, which is an increase of 153%. As for the low-carbon group, the difference between the two frontiers could be neglected. Notably, there are significant differences among the CEP scores of the three groups when the meta-frontier is used as the basis of evaluation, while relatively moderate distinctions are found when using the group frontiers.
4.1.3. Mitigation Potential
Table 7 reports the estimates of Mitigation Potential under the meta-frontier (MMP) and the group-frontier (GMP). From the table, it can be observed that MMP and GMP vary from 0% to 88.8% and from 0% to 68.7%, respectively. Moreover, the mean values of MMP are generally higher than those of GMP. From the perspective of total sectors, the average MMP is found to be 54.4%, which is more than twice the mean value of GMP (22.5%), indicating that there is a significant difference. Additionally, the differences in the average MP values between the two frontiers are found to be considerable at group levels, especially for the high-carbon and medium-carbon groups. Specifically, the average MMP of the high-carbon group (68.8%) is found to be almost three times that of GMP (23.4%). Further, the average MP value of the medium-carbon group varies between the meta-frontier and group-frontier at greater levels, which is found to jump from 16.9% to 67.2%. As for the low-carbon group, the difference between the two frontiers could be neglected. Notably, there are significant differences among the MP values of the three groups when the meta-frontier is used as the basis of evaluation, while relatively moderate distinctions are observed when using the group frontiers.
4.1.4. Marginal Abatement Cost
Table 8 reports the estimates of the marginal abatement cost under the meta-frontier (MMAC) and the group-frontier (GMAC). From the table, it is found that MMAC and GMAC vary from 240 to 12,970 Yuan/ton and from 160 to 16,160 Yuan/ton, respectively. From the perspective of the meta-frontier, the low-carbon group is found to have the highest weighted average MMAC of 4970 Yuan/ton, followed by the medium group with a value of 1310 Yuan/ton, and the high-carbon group is at the bottom, with a value of 500 Yuan/ton. On the other hand, the weighted average GMAC values of the three groups also follow the same order, which are found to be 10,320, 2080, and 400 Yuan/ton. The differences between the mean values of MMAC and GMAC at group levels are considerable, especially for the low-carbon group, which is found to have a twofold change in the values. Additionally, the average weighted MAC values of the total sectors under the group-frontier and the meta-frontier are found as 670 and 770 Yuan/ton, respectively. As a matter of fact, the sectors have to sacrifice relatively more IVA for reducing an additional unit of CO
2 emissions, if they are producing at the more efficient overall production frontier with relatively less mitigation potential. So, in general, marginal abatement costs are negatively correlated with the mitigation potential.
It is interesting to compare our MAC estimates with those of previous studies. A summary of MAC estimates is reported in
Table 9. From the table, it can be observed that the estimates of the previous studies lie in different ranges. The different samples and approaches in these studies are the main reasons for the differences in the results. The results of this study are found to be apparently lower than those of Xiao et al. [
7] in terms of the mean and weighted mean of MAC. This is reasonable because our study only covers the emission-intensive industries of the industrial sector, while Xiao et al. [
7] consider all industrial sectors as research objects, which is more expensive to reduce emissions. Additionally, the significant disparity of the average MAC between Yuan et al. [
40] and Xie et al. [
11] can also be attributed to the difference in research objects.
4.2. Sectoral Analysis
4.2.1. Distribution Analysis
Figure 1 indicates the estimates of CEP, MP, and MAC under the meta-frontier and the group-frontier at the sector level. From the figure, we can observe that there exist significant disparities in the CEP, MP, and MAC under both frontiers among the 39 sectors. Specifically, from the meta-frontier perspective, the sectoral distribution of MCEP presents an increasing trend from S1 to S39 (i.e., the decreasing order of carbon intensity), while that of MMP presents the opposite tendency. Similarly, the sectoral distributions of GCEP and GMP show the upward and downward trends along with the decrease of carbon intensity in their respective groups, respectively. For instance, in the high-carbon group, the sectoral distribution of GCEP presents an increasing tendency from S1 to S13, while a decreasing trend is observed for GMP. On the other hand, it is found that both of the sectoral distributions of MMAC and GMAC present an increasing tendency from S1 to S39.
In light of this interesting phenomenon, we conducted a correlation analysis between MMAC, GMAC, and carbon intensity, respectively, and the results are listed in
Table 10. According to the correlation results, there exists a significantly negative correlation (
p = 0.000 < 0.01, Pearson’s r = −0.537) between carbon intensity and MMAC, which has been previously reported in the literature, such as in Zhou et al. [
10]. Additionally, GMAC is also found to be significantly negatively correlated with carbon intensity (
p = 0.000 < 0.01, Pearson’s r = −0.566).
In order to provide an intuitive understanding of the distribution patterns of CEP, MP, and MAC, and make a comparison between the meta-frontier and the group-frontier, the kernel density curves of CEP, MP, and MAC under both frontiers are plotted based on the R-program. From
Figure 2, it can be found that there exist some differences in the distribution patterns of CEP, MP, and MAC between the two frontiers. Specifically, the kernel curve of CEP moves rightward and the dispersion range of points become narrower when the evaluation basis switches from the meta-frontier to the group-frontier, which means an increase in the mean value and a decrease in the variance of CEP, respectively. In the meanwhile, the kernel curve of MP moves leftward and the dispersion range of points become narrower, which means that the mean value and the variance of MP both decrease. On the contrary, the kernel curve of MAC holds the position unchanged, but exhibits a significant decrease in the peak value and an expansion in the dispersion range of points when the evaluation basis changes, which jointly means a considerable increase in the variance of MAC.
4.2.2. Difference Analysis
Figure 3 reports the differences in CEP, MP, and MAC between the meta-frontier and the group-frontier at sector levels. From the figure, we can observe that there exist significant disparities between the two frontiers in terms of CEP, MP, and MAC, and the differences exhibit unequal distributions among the 39 sectors. The difference in CEP (D_CEP) is found to be significant in the high-carbon and medium-carbon sectors. For instance, the D_CEP of Pressing of Steel (S2) is found to be 0.691. Similar results can be observed for sectors such as the Manufacture of Fireproof Materials (S13), Manufacture of Special Chemical Products (S19), and Manufacture of Explosives and Fireworks (S25). This means that there is significant technology heterogeneity between the meta-frontier and group-frontier for the two groups. Conversely, it is found that there is no difference in the low-carbon sectors, which may be attributed to the fact that low-carbon sectors play an important role in benchmarking efficiency levels under meta-frontier technologies, and thereby there are no technology gaps for the low-carbon sectors. Similarly, the difference in MP (D_MP) between the two frontiers is also found to be significant in the high-carbon and medium-carbon sectors, while the opposite is true for the low-carbon sectors.
On the other hand, the difference in MAC (D_MAC) exhibits a completely different situation, in which the D_MAC is found to be larger, equal to, and smaller than 0, depending on the observations. For instance, the Pressing of Non-ferrous Metals (S32) and Casting Pressing of Non-ferrous Metals (S37) are found to have the lowest and the highest D_MAC, with values of −5610 Yuan/ton and 8680 Yuan/ton, respectively. There is no definite relationship between the MMAC and GMAC of sectors; the MMAC of sectors can be larger, equal to, or smaller than the GMAC, depending on the relative slopes of the meta- and the group-frontiers, respectively [
19]. However, as mentioned earlier, both MMAC and GMAC are significantly correlated with carbon intensity. Therefore, it is interesting to explore the relationship between D_MAC and carbon intensity. In this regard, the correlation analysis between D_MAC and carbon intensity was conducted, and the results are listed in
Table 10. According to the correlation results, there is also a significantly negative correlation (
p = 0.047 < 0.05, Pearson’s r = −0.320) between D_MAC and carbon intensity.
4.3. Discussions
In general, whether from the perspective of group levels or from the perspective of sector levels, it is found that there exist significant disparities in the CEP, MP, and MAC under both frontiers among various sectors. Additionally, the differences between the two frontiers in terms of CEP, MP, and MAC are considerable, and exhibit unequal distributions among the 39 sectors. This can be attributed to the significant heterogeneity of production technology among various sectors. Considering that carbon intensity is widely considered as the measurement of carbon emission efficiency, the industries with a relatively low carbon intensity are more efficient than those with a high carbon intensity in terms of energy utilization and carbon emissions. Thus, low-carbon sectors are found to have a higher CEP, less MP, and larger MAC than medium-carbon and high-carbon sectors.
As the most efficient DMUs, low-carbon sectors are found to play a more important role in benchmarking efficiency levels under the meta-frontier. In this case, no difference between the meta-frontier and the group-frontier for the low-carbon group in terms of CEP and MP exists, while the opposite is true for the medium-carbon and high-carbon sectors. As a result, considerable differences in the CEP and MP between the two frontiers are observed for the medium-carbon and high-carbon sectors. As for MAC, theoretically, the MMAC can be larger, equal to, or smaller than the GMAC, depending on the relative slopes of the meta-frontier and the group-frontier. Nevertheless, it is found that D_MAC has a significantly negative correlation with carbon intensity.
5. Conclusions and Policy Implications
CN-ETS that covers seven emission-intensive industries is scheduled to be launched in 2017. In this context, it is of great urgency and necessity to obtain a good understanding of participating sectors in terms of energy utilization and carbon emissions. In this regard, estimating the CEP, MP, and MAC for these sectors can provide valuable information for the governments and participating enterprises. Therefore, taking the industry heterogeneity into consideration, we employed a joint framework consisting of the DDF and meta-frontier analysis to estimate CEP, MP, and MAC under the meta-frontier and the group-frontier, respectively. Following this, we investigated the sectoral distributions of CEP, MP, and MAC under both frontiers, and analyzed the differences between the two frontiers in terms of CEP, MP, and MAC at sector levels.
Based on the detailed analysis, the main conclusions are drawn as follows: First, there exist significant disparities in the CEP, MP, and MAC under both frontiers among various sectors. Specifically, high-carbon and medium-carbon sectors are found to display a low CEP, large MP, and high MAC, while the opposite situation is observed for low-carbon sectors which have a high CEP, small MP, and high MAC. Furthermore, the sectoral distributions of CEP, MP, and MAC are found to be different between the two frontiers. Additionally, the differences between the two frontiers in terms of CEP, MP, and MAC are considerable, and exhibit unequal distributions among the 39 sectors. Notably, the MAC values under both frontiers and the difference between them are all found to be significantly correlated with carbon intensity.
Based on the above conclusions, possible policy implications are provided for the government and participating enterprises, respectively. First of all, based on the estimates of CEP, MP, and MAC for sectors, our study allows policy makers to pinpoint where the greatest emissions cuts—at the least expense—can be made in China’s emission-intensive sectors. From the perspective of the government, different policies should be implemented for the critical emission reduction sectors based on their CEP, MP, and MAC. In this regard, high-carbon sectors with a low CEP, large MP, and low MAC should shoulder more responsibility for the reduction of emissions. In particular, necessary policy measures such as introducing carbon-emission-reduction technologies, eliminating backward production facilities, raising the industry entry threshold, and increasing the industry concentration should be implemented for high- and medium-carbon sectors. Furthermore, the CEP, MP, and MAC (especially GCEP, GMP, and GMAC) of participating sectors should be taken into consideration when formulating the criteria for the initial allocation of carbon allowances before transactions. Additionally, the weighted average MAC could offer a reference for the carbon price of CN-ETS, since the industries covered in CN-ETS are taken as the research objects. On the other hand, from the perspective of enterprises, by comparing the MAC of participating sectors with the carbon price of CN-ETS, participating enterprises could identify the least-costly emission reduction strategy from a list of policy alternatives such as abating carbon emissions, buying carbon allowances, and selling carbon allowances. In this context, high- and medium-carbon enterprises tend to cut their emissions and sell carbon allowances, while low-carbon enterprises may choose to buy carbon allowances.
Despite the contributions, this paper has a limitation in assuming the same growth rate for various sectors in light of their heterogeneity. Moreover, our study cannot consider dynamic effects since we are limited to using a one-year data cross-section. Additionally, due to the limitation of statistical data, this study selects 39 sectors as research objects rather than the specific enterprises covered in CN-ETS. However, this study could be easily extended to enterprises, so future research requires the collection of data at enterprise levels in order to obtain more accurate estimates.