Abstract

In China, the transportation sector contributes about 18% of the total carbon emissions. This research contributes to measuring the energy and carbon emission efficiency (ECEE) of regional transportation systems (RTS) in China considering uncertain carbon emissions. A radial chance-constrained data envelopment analysis (DEA) model is developed to estimate the overall efficiency, and a nonradial chance-constrained DEA model is presented to evaluate the pure energy efficiency (PEE) and the pure carbon emission efficiency (PCEE). We prove that the proposed chance-constrained DEA models can effectively address the uncertain carbon emissions when measuring efficiency. We find that most of China’s RTS have low ECEE and the inefficiencies are mainly due to the lower gasoline utilization efficiency and the lower kerosene utilization efficiency. In addition, east China performs better than central China, and central China performs better than west China. In China, the unbalanced regional development of the ECEE in transportation corresponds with the unbalanced regional economic development. We provide some valuable suggestions based on the evaluation of the potential cuts in each kind of energy and the potential decreases in carbon emissions.

1. Introduction

Since China started economic reforms in 1978, the country has been experiencing rapid economic development [1, 2]. According to the International Monetary Fund (IMF), China had a nominal Gross Domestic Product (GDP) of about 13.37 trillion U.S dollars in 2018. The nominal GDP value of China accounted for approximately 15.7 percent of the global economy [3]. IMF rated China as the largest economy by purchasing power parity (PPP) in the world [4].

China’s rapid economic growth has significantly expanded its energy consumption [5]. Therefore, China should make every effort to shift to a green economy. Chinese government has been implementing many policies and measures to increase energy efficiency and reduce carbon emissions. By 2020, China intends to contract carbon dioxide (CO2) emissions per unit of GDP by 18% compared with the level in 2015 [6].

The energy consumption and carbon emissions in the transportation sector are very high [7, 8]. The total energy consumption in the transportation sector of China was approximately 402 million tons of standard coal equivalent in 2017, an increase of 4.6% on the previous year. The transportation energy consumption in east China, central China, and west China (the detailed information of the three areas is shown in Section 3) occupied about 49% (197 million tons), 32% (127 million tons), and 19% (78 million tons) of the total energy consumption, respectively. In China, the transportation sector contributes about 18% of the total carbon emissions [9]. Hence, it is necessary to research the energy and carbon emission efficiency (ECEE) of the transportation sector in China.

Several mathematical approaches can be used to evaluate the ECEE. However, data envelopment analysis (DEA) is most widely used (see the following paragraphs) [10]. DEA is a nonparametric mathematical method to measure the efficiency of decision-making units (DMUs) [11, 12]. DEA has been applied to evaluate the performance of hospitals, supply chains, universities, and other entities [13]. The CCR model (proposed by Charnes et al.) and the BCC model (proposed by Banker et al.) are two basic DEA models [14, 15]. The former assumes that technology exhibits constant returns to scale (CRS), and the latter is developed under the condition of variable returns to scale (VRS). In the past forty years, many DEA models, e.g., integer-valued DEA models [16], super-efficiency DEA models [17], generalized DEA models [18], slack-based DEA models [19], network DEA models [20], and stochastic DEA models [21], have been proposed based on the two models.

Scholars have been studying the energy efficiency measurement for a long time. Ramanathan applied the CCR model and the BCC model to evaluate the energy efficiency of railways and roadways in India [22]. Hu and Wang used the DEA method to measure the total-factor energy efficiency of 29 regions in China [23]. Wu et al. proposed a DEA technique with nonhomogeneous variables to estimate the energy efficiency of industrial sectors in China [24]. Geng et al. proposed a DEA approach with the affinity propagation clustering algorithm (AP-DEA) to evaluate the energy efficiency of petrochemical industries [25]. However, undesirable outputs, e.g., carbon emissions, are not considered in the aforementioned DEA models. These energy efficiency evaluation models are not suitable in reality because the fossil fuel energy consumption will definitely generate carbon emissions [26]. Therefore, the efficiency scores obtained from these models may be distorted.

Zhou et al. proposed slacks-based measures (SBM) to evaluate the CO2 emission efficiency of 30 OECD countries (the Organization for Economic Cooperation and Development) [27]. Zhou et al. applied DEA approaches to evaluate the carbon emission efficiency of eight world regions [28]. Wang et al. developed DEA approaches to compare the environmental efficiency of several provinces in China by taking CO2 emissions into account [29]. Wei et al. evaluated the CO2 emission efficiency of China’s 29 provinces by an extended SBM-DEA model [30]. Wang et al. applied a DEA window analysis approach to evaluate the regional environmental efficiency in China considering CO2 emissions [31]. Guo et al. developed a dynamic DEA model to measure the energy efficiency of OECD countries by taking CO2 emissions into consideration [32]. Cheng et al. used a time substitution DEA approach to estimate the carbon emission efficiency of China [33]. Chu et al. proposed an SBM-DEA approach to evaluate the environmental efficiency of China’s transportation system considering CO2 emissions [34]. Cui developed a network DEA approach to measure the greenhouse gas emission efficiency of 28 international airlines [35]. Djordjevic and Krmac applied a nonradial DEA and TOPSIS (the Technique for Order of Preference by Similarity to Ideal Solution) approach to measure the environmental efficiency of European transport systems by taking greenhouse gas emissions into account [36]. Zhou et al. proposed a DEA method according to an exponential transformation to evaluate the energy efficiency of China’s industry considering greenhouse gas emissions [37]. Omrani et al. developed a DEA-cooperative game method to estimate the energy efficiency of the transportation system in Iran by taking greenhouse gas emissions into account [38]. The disadvantage of the abovementioned literature is that they assume deterministic carbon emissions. However, it is impossible to get precise data related to carbon emissions, which could be caused by measurement errors, natural uncertainty of carbon emissions, and so on [39]. The results would be erroneous if decision-makers do not consider the uncertainty of carbon emissions when measuring the ECEE [3943]. As a matter of fact, the deterministic DEA model is only a particular case of the uncertain DEA model [39].

Wu et al. [40], Jin et al. [41], Zha et al. [39], Chen et al. [42], and Zhou et al. [43] used the chance-constrained DEA method to evaluate the environmental efficiency considering uncertain undesirable outputs. However, as far as we know, there is no literature on measuring the efficiency of regional transportation systems (RTS) by taking the uncertainty of carbon emissions into account, especially China’s RTS. It is necessary to research the RTS because of the unbalanced regional development in China [44]. The Chinese government should set up different policies for different RTS to improve their ECEE.

The innovation of this paper is that, to the best of our knowledge, this is the first paper to measure the efficiency of RTS by taking uncertain carbon emissions into consideration. The contributions of this paper are that (1) a radial chance-constrained DEA model is presented to evaluate the overall efficiency; (2) a nonradial chance-constrained DEA model is developed to accurately estimate the pure energy efficiency (PEE) and the pure carbon emission efficiency (PCEE); (3) the proposed approaches are applied in the RTS in China and some valuable suggestions are provided.

This paper is structured as follows. The chance-constrained DEA models are developed by taking uncertain carbon emissions into consideration in Section 2. In Section 3, these models are applied to estimate the ECEE of thirty RTS in China, and some valuable suggestions are provided. Some conclusions are proposed in Section 4.

2. Methodology

A regional transportation system can be regarded as a DMU. Suppose that has nonenergy inputs , energy inputs , desirable outputs , and an undesirable output (carbon emissions).

To measure the environmental efficiency, Färe et al. handled undesirable outputs under the assumption of the weak disposability of undesirable outputs [45]. Nowadays, this method has been widely used. For given levels of inputs, it is impossible to simultaneously reduce carbon emissions and expand desirable outputs, e.g., passenger turnover, freight turnover, and transportation value-added. Therefore, to estimate the overall efficiency of RTS, we propose a radial DEA model (model 1) under the assumption of the weak disposability of the undesirable output (the amount of carbon emissions) [45, 46]. Model (1) is a radial DEA model as this model aims to proportionally reduce inputs and the undesirable output of [22, 45].where is the weight of ; is a non-Archimedean infinitesimal; and , , and are the slack variables for the nonenergy input, the energy input, and the desirable output, respectively; and respectively represent the assumption of CRS and VRS.

The optimal value is the technical efficiency score (or pure technical efficiency score) of when (or ). According to the definition of DEA efficiency, we have the following propositions: (i) is efficient when ( and ); (ii) is inefficient when .

Model (1) is called the radial deterministic DEA model because it can only be used when all the needed data are precise. However, as discussed above, the value of carbon emissions is uncertain. Therefore, we have to develop a stochastic DEA model. To do so, we first propose model (2) that is the dual of model (1). All slack variables are deleted in the following transformation processes for the sake of clarity.where , , , and are the weights for the desirable output, the nonenergy input, the th energy input, and the undesirable out, respectively.

As stated in Section 1, the stochastic chance-constrained DEA is an effective method to evaluate the environmental efficiency considering uncertain undesirable outputs. Thus, we also propose a radial stochastic chance-constrained DEA model (model 3) based on model (2) to measure the overall efficiency of RTS by taking uncertain carbon emissions into consideration.where indicates stochastic carbon emissions; represents the expected value of ; is the “probability” measurement of the corresponding constraint; and is the probability of satisfying the inequality . In other words, indicates a confidence level while represents a risk criterion [47, 48].

As model (3) is a nonlinear model, we should transform it into a linear model. Obviously, inequality (4) is only a transformation of the stochastic chance constraint in model (3).

Let us define that denotes the expected value of and indicates the variance of . The inequality (4) then can be transformed into the following inequality:

Assume that obeys (standard normal distribution) and denotes its cumulative distribution function. The inequality (5) then can be transformed into the following inequality:where represents the inverse function of . We assume and , where represents ’s expected value and denotes the standard deviation. Assume that follows . Then, we have , and model (3) can be transformed into the following model:

The dual of model (7) can be denoted by model (8) that aims to proportionally reduce inputs and uncertain carbon emissions of . It is worth noting that slack variables are added into model (8). We have transformed nonlinear model (3) into linear model (8) that is called the radial chance-constrained DEA model. The optimal objective function value of model (8) is the technical efficiency score (or pure technical efficiency score) of when (or ).

According to the definition of DEA efficiency, we have the following propositions: (i) is superefficient when ; (ii) is efficient when ( and ); (iii) is inefficient when .

The advantages of the radial chance-constrained DEA model are that (1) it can be used to deal with undesirable outputs that cannot be disposed of; (2) it can be used to handle stochastic variables. However, model (8) cannot be used to evaluate the PEE or PCEE. Therefore, we have to extend it to a nonradial chance-constrained DEA model (model 9). Model (9) optimizes the carbon emissions and energy consumption without changing the desirable outputs and nonenergy inputs. We evaluate the ECEE by assigning weights to the PEE and PCEE [29, 39].where indicates the PEE (pure efficiency of energy ); is the PCEE; and respectively represent the weights of the PEE and PCEE; and is the weight of the th energy input. There must be and . The PEE indicates the utilization performance of transportation fuels, and the PCEE reflects the performance of carbon emission cuts.

For an inefficient regional transportation system, its potential energy savings and carbon emission cuts can be calculated using the following equations, respectively:

3. Empirical Study

The proposed models are used to estimate the ECEE of China’s RTS in 2017.

3.1. Data

Mainland China includes 31 provinces, and the regional disparity is wide. These provinces can be classified into three areas, i.e., the eastern area, the central area, and the western area. The eastern area, which is the most developed and densely populated area, comprises of eleven provinces; The central area includes ten provinces; The western area, which is the least developed and sparsely populated area, is composed of ten provinces [44]. The data about transportation energy consumption in Tibet are not available so that we have to exclude Tibet. Therefore, thirty provinces are involved in this research.

We select employed persons (labor), fixed-asset investment (capital), and civil motor vehicles as three nonenergy input variables. Gasoline consumption, kerosene consumption, and diesel consumption are regarded as three energy input variables. Passenger turnover, freight turnover, and transportation value-added are selected as three desirable output variables. The amount of carbon emissions from transportation is defined as a single undesirable output.

The chosen nonenergy input, desirable outputs, and undesirable outputs in this paper are vital indicators for transport systems and have been used widely in previous research [36, 38]. Gasoline, kerosene, and diesel are selected as three energy inputs because they are the dominant transportation fuels in China. As shown in Figure 1, diesel makes up more than 40% of the total transportation energy consumption, gasoline occupies more than 20%, and kerosene accounts for more than 10% [49]. We cannot select all transportation fuels as energy inputs according to the rule of thumb [50].

The data related to employed persons (unit: 10 thousand persons), civil motor vehicles (unit: 10 thousand vehicles), fixed-asset investment (unit: 1 billion CNY), passenger turnover (unit: 1 billion passenger/km), freight turnover (unit: 1 billion ton/km), and transportation value-added (unit: 1 billion CNY) are obtained from [51]. The data related to energy consumption (unit: 10 thousand tons of standard coal equivalent) are calculated using the following equation:where (unit: 10 thousand tons of standard coal equivalent) represents the consumption of energy by ; (unit: 10 thousand tons) also indicates the consumption of energy by ; and denotes the standard coal coefficient of the energy (gasoline, 1.4714; kerosene, 1.4714; diesel, 1.4571) [52]. The values of are obtained from [49].

The amount of carbon emissions from transportation (units: 10 thousand tons) can be calculated using the following equation:where represents the carbon emissions of ; indicates the consumption of the energy by ; and denotes the carbon emission coefficient for the energy (gasoline, 0.5538; kerosene, 0.5714; diesel, 0.5921; coal, 0.7559; fuel oil, 0.6185; natural gas, 0.4483) [53]. Although we only select three transportation fuels as energy inputs, we want to calculate carbon emissions generated from all transportation fuels. Note that the electric power, which can be generated by wind power and hydropower, is not included in our research because the carbon emission coefficients of both wind power and hydropower are zero.

We suppose that the expected carbon emissions from transportation are equal to the estimated carbon emissions. Therefore, we have . The standard deviation of carbon emissions from transportation is calculated using the dataset of the estimated carbon emissions from 2013 to 2017. The descriptive statistics of all selected variables of RTS in mainland China are shown in Table 1. Table 2 shows the detailed data related to carbon emissions from transportation.

3.2. Results
3.2.1. Overall Efficiency

The overall efficiency of the thirty RTS in mainland China is calculated using the proposed radial chance-constrained DEA model (model 8). We also compare the efficiency scores resulting from model (8) with those resulting from the radial deterministic DEA model (model 1).

Tables 3 and 4 respectively show the efficiency scores under the assumption of CRS and VRS . We use four groups of values (, , , and ) to study the effects of on the results of efficiency evaluation, and we select as the risk criterion in this study as it is set 0.05 in the related literature [3943].

As shown in Table 3, the efficiency scores resulting from model (1) are the same as those resulting from model (8) with , which implies that the radial deterministic DEA model (model 1) is only a particular case of the radial chance-constrained DEA model (model 8). Furthermore, only four RTS are regarded as inefficient DMUs based on their overall technical efficiency scores resulting from model (1). However, based on their overall technical efficiency scores resulting from model (8), we have thirteen inefficient RTS when , , or . The results prove that decision-makers should take the uncertainty of carbon emissions into account when they measure the efficiency of RTS. Otherwise, the results may be doubtful. As shown in Table 4, we can get the same results under the assumption of VRS .

The proposed radial chance-constrained DEA model (model 8) is more effective than the radial deterministic DEA model (model 1) because (i) it has higher discriminatory ability than model (1) as stated above; (ii) we cannot rank efficient RTS because all of their efficiency scores obtained from model (1) are equal to one (for example, as shown in Table 3, Beijing, Tianjin, Hebei, Liaoning, Shanghai, Jiangsu, Zhejiang, Shandong, Guangdong, Hainan, Shanxi, Nei Mongol, Jilin, Heilongjiang, Anhui, Jiangxi, Henan, Hubei, Hunan, Chongqing, Guizhou, Shaanxi, Gansu, Qinghai, Ningxia, and Xinjiang.), but their efficiency scores resulting from the radial chance-constrained DEA model can be greater than one, which makes it possible to differentiate the efficient RTS.

As shown in Tables 3 and 4, the efficiency scores (technical efficiency score and pure technical efficiency score) of Qinghai are the highest, while Yunnan has the lowest scores. Therefore, Yunnan should make every effort to improve the performance of its transportation system.

As shown in Table 5, based on the overall efficiency scores of the three areas under both CRS and VRS, the performance of east China is the best while west China is the worst. It indicates that the unbalanced regional development of overall efficiency in transportation corresponds with the unbalanced regional economic development.

3.2.2. PEE and PCEE

The proposed nonradial chance-constrained DEA model (model 9) is used to evaluate the PEE and PCEE of China’s RTS. We set because carbon emissions are generated by combusting energy fuels. Moreover, the weights for the three energy inputs can be calculated using the following equation [39]:where is the carbon emission coefficient of energy (as shown in Section 3.1). We finally obtain , , and .

Tables 6 and 7 show the efficiency scores of these RTS resulting from the nonradial chance-constrained DEA model (model 9) under the assumption of CRS and VRS , respectively. It can be observed that the transportation system of Ningxia does not have feasible solutions for its kerosene utilization efficiency. The reason is that the consumption of kerosene in Ningxia’s transportation system is zero.

Table 6 shows that the technical efficiency score of energy utilization and carbon emissions in Qinghai is the highest while Yunnan has the lowest score. The results are consistent with those obtained from model (8) under the CRS condition. The technical efficiency score of gasoline utilization in Henan is the highest while Heilongjiang has the lowest score. The technical efficiency score of kerosene utilization in Shanxi is the highest while Hainan has the lowest score. The technical efficiency score of diesel utilization in Heilongjiang is the highest while Yunnan has the lowest score. The technical efficiency score of carbon emissions in Qinghai is the highest while Yunnan has the lowest score.

The average , , , , and under CRS are 0.745, 0.673, 0.618, 0.877, and 0.772, respectively, which implies that the technical inefficiencies of energy utilization and carbon emissions in transportation are mainly due to the lower technical efficiency of gasoline utilization and the lower technical efficiency of kerosene utilization .

Table 7 shows that the pure technical efficiency score of energy utilization and carbon emissions in Qinghai is the highest while Yunnan has the lowest score. The results are consistent with those obtained from model (8) under the assumption of VRS. There is no pure technical superefficient province of gasoline utilization or kerosene utilization. The pure technical efficiency score of gasoline utilization in Heilongjiang is the lowest, and the pure technical efficiency score of kerosene utilization in Yunnan is the lowest. The pure technical efficiency score of diesel utilization in Heilongjiang is the highest while Yunnan has the lowest score. The pure technical efficiency score of carbon emissions in Qinghai is the highest while Yunnan has the lowest score.

The average , , , , and under VRS are 0.833, 0.741, 0.697, 0.966, 0.869, respectively, which implies that the pure technical inefficiencies of energy utilization and carbon emissions in transportation are mainly due to the lower pure technical efficiency of gasoline utilization and the lower pure technical efficiency of kerosene utilization .

As shown in Table 8, under both CRS and VRS, the ECEE of the eastern area is the best while the western area is the worst. The results are consistent with those obtained from model (8). The PCEE of the eastern area is also the best under both CRS and VRS. However, the rankings of the central area and the western area are different under different assumptions of production technology. The PEE (, , or ) of the eastern area is not always the highest in mainland China. For instance, the eastern area has the lowest technical efficiency score of gasoline utilization.

3.3. Discussion and Suggestion

The uncertainty of carbon emissions can significantly affect the results of efficiency evaluation. Therefore, decision-makers should take uncertain carbon emissions into account when they measure the performance of RTS [54, 55]. The risk criterion also has effects on the results of efficiency evaluation. Hence, a scientific risk criterion is required when decision-makers measure the performance of DMUs using the proposed chance-constrained DEA models.

The performance of Yunnan’s transportation system is the worst in mainland China. It has the lowest overall efficiency scores (the technical efficiency score and the pure technical efficiency score), the lowest pure technical efficiency score of kerosene utilization, the lowest efficiency scores (the technical efficiency score and the pure technical efficiency score) of diesel utilization, and the lowest efficiency scores (the technical efficiency score and the pure technical efficiency score) of carbon emissions. Therefore, Yunnan should put forth the greatest possible effort to increase the efficiency of its transportation system. The transportation system of Heilongjiang should make every effort to increase its efficiency of gasoline utilization. The transportation system of Hainan should put forth its best effort to increase its technical efficiency of kerosene utilization.

In China, the unbalanced regional development of the ECEE in transportation corresponds with the unbalanced regional economic development. East China performs better than central China, and central China performs better than west China. Hence, China should try its best to improve the performance of RTS in the western area.

Figure 2 shows the potential cuts in each kind of energy and the potential decreases in carbon emissions of each inefficient regional transportation system under the assumption of CRS. It is noteworthy that there are twelve efficient DMUs, i.e., Tianjin, Hebei, Liaoning, Shanghai, Shandong, Anhui, Jiangxi, Hunan, Guizhou, Gansu, Qinghai, and Ningxia (as shown in Table 6). It can be observed that the eastern area (especially, Jiangsu, Zhejiang, and Guangdong) can save many energy fuels and reduce lots of carbon emissions as the eastern area is the largest transportation fuel consumer and the largest carbon emitter in the three areas. Hence, China should also improve the efficiency of RTS in the eastern area.

4. Conclusion

This study proposes chance-constrained DEA models to estimate the energy and carbon emission efficiency (ECEE) of regional transportation systems (RTS) considering uncertain carbon emissions. The overall efficiency of thirty regional transportation systems in mainland China is calculated using the radial chance-constrained DEA model. We prove that the radial deterministic DEA model is only a particular case of the radial chance-constrained DEA model. The pure energy efficiency (PEE) and the pure carbon emission efficiency (PCEE) of the thirty regional transportation systems are estimated using the proposed nonradial chance-constrained DEA model. We prove that this model can provide more valuable information for decision-makers, i.e., how to save each transportation fuel and reduce carbon emissions. We find that most of China’s regional transportation systems have low energy and carbon emission efficiency and the inefficiencies are mainly due to the lower gasoline utilization efficiency and the lower kerosene utilization efficiency.

We suggest that (1) China should try its best to improve the performance of the regional transportation systems in the western area because its efficiency is the lowest; (2) China should also improve the performance of the regional transportation systems in the eastern area because the eastern area can save many energy fuels and reduce lots of carbon emissions (especially, Jiangsu, Zhejiang, and Guangdong).

Measuring the energy and carbon emission efficiency (ECEE) of regional transportation systems is an important and hot topic. The results of this study are helpful for the Chinese government to set up different policies for different RTS to improve their ECEE. In fact, the proposed approaches can be also applied to measure the RTS in other countries. In addition, we can measure the ECEE of RTS considering other uncertain inputs and outputs.

Data Availability

The research data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Authors’ Contributions

Jianwei Ren and Chunhua Chen were responsible for the conceptualization of the study; Jianwei Ren and Jiewei Zhang curated the data; Jianwei Ren was responsible for funding acquisition; , Jianwei Ren and Chunhua Chen was responsible for methodology; Software, Bin Gao and Chunhua Chen provided software; Jianwei Ren, Jiewei Zhang, and Chunhua Chen wrote the original draft; Bin Gao and Chunhua Chen performed reviewing and editing.

Acknowledgments

This research was funded by the National Natural Science Foundation of China (Grant no. 71862026), the China Postdoctoral Science Foundation (Grant no. 2018T110209), the Natural Science Foundation of Inner Mongolia (Grant no. 2018MS07006), and the State Scholarship Fund of China Scholarship Council (Grant no. 20180815502).