Next Article in Journal
Hydrodynamics of Moonpool-Type Floaters: A Theoretical and a CFD Formulation
Next Article in Special Issue
A Concept for Solving the Sustainability of Cities Worldwide
Previous Article in Journal
A Proportional Digital Controller to Monitor Load Variation in Wind Turbine Systems
Previous Article in Special Issue
A Strategy for Planned Product Aging in View of Sustainable Development Challenges
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Examining the Relationship between Energy Consumption and Unfavorable CO2 Emissions on Sustainable Development by Going through Various Violated Factors and Stochastic Disturbance–Based on a Three-Stage SBM-DEA Model

1
School of Business, Macau University of Science and Technology, Taipa, Macau 999078, China
2
Faculty of Hospitality and Tourism Management, Macau University of Science and Technology, Taipa, Macau 999078, China
*
Author to whom correspondence should be addressed.
Energies 2022, 15(2), 569; https://doi.org/10.3390/en15020569
Submission received: 10 December 2021 / Revised: 5 January 2022 / Accepted: 11 January 2022 / Published: 13 January 2022
(This article belongs to the Special Issue Green Economy and Sustainable Development)

Abstract

:
The article applies a three-stage Slacks-Based Measure-Data Envelopment Analysis (SBM-DEA) pattern to examine the relationship between energy consumption and unfavorable CO2 emissions on green sustainable development, for the 11 cities of the Guangdong-Hong Kong-Macao Greater Bay Area (GBA) during 2010–2016, by going through various violated factors and stochastic disturbance. Labor, capital and energy resource are chosen as input variables, while GDP and CO2 emission as output variables. During the three phases consisting of the SBM-DEA model (first stage and third stage) and SFA analysis (second stage), CO2 emission is considered as an unfavorable outcome, while stochastic statistical disturbances and external environmental influences are identified. The results show that the average efficiency of the GBA cities is 0.708, with only Shenzhen, Macao SAR and Hong Kong SAR having an efficiency of 1 during the whole study period. Based on the findings, suggestions are made for the GBA cities’ sustainable development aspects.

Graphical Abstract

1. Introduction

The World Commission on Environment and Development defined sustainable development as “development that meets current demands without jeopardizing future generations’ ability to meet their needs” in 1987 [1]. In June 1992, in Rio de Janeiro, Brazil, at the United Nations Conference on Environment and Development, Agenda 21 first proposed sustainable development as the common development strategy of mankind toward the 21st century, transforming sustainable development strategy from a concept to a global scale of action and making sure that the most central part is economic sustainable development. Liu (1997) defines economic sustainable development as, “What we call sustainable development economy can be expressed as sustainable economic development, which should be the economy with the lowest ecological cost and social cost of economic development” [2]. Such economic sustainable development is economic development under the interests of future generations and the protection of natural resources. In addition to economic sustainability being at the core of sustainable development, energy consumption and environmental protection are also major focus of sustainable development strategies. Both show a one-way Granger causality and an inverted U-shaped relationship with sustainable development, respectively. In the early days, labor and energy-intensive industries were used for high economic development. But now, with sustainable development as a core value, cities are improving their energy consumption efficiency and increasing investment in environmental protection. At the same time, it is expected that they would attain sustainable development in the future. Here, the Greater Bay Area (GBA) in China is used as an example of sustainable development.
Figure 1 shows the total energy consumption, total CO2 emission and total GDP of the GBA cities from 2010–2016. Although the total GDP kept on increasing since 2010, the total energy consumed has not changed tremendously. The amount of total CO2 emission also demonstrated a decline to some extent.
Improving energy efficiency and considering undesirable output concurrently is one of the mainstreams in achieving the goal of sustainable development. However, most of the existing research literature is divided by provincial level or only includes 9 cities in Guangdong Province, China (Guangzhou, Shenzhen, Zhuhai, Foshan, Jiangmen, Zhaoqing, Huizhou, Dongguan, and Zhongshan), and few studies include Hong Kong SAR and Macau SAR, but GBA is a complete regional concept, so this study hopes to include data from Hong Kong SAR and Macau SAR by using the three-stage SBM-DEA model to study the complete DBA energy efficiency data for a total of 11 cities in GBA, the better usage of energy consumption toward green sustainable development can be investigated and use the data from these studies to provide better policy recommendations for those inefficient GBA cities, and eventually build an efficient and green GBA.
The followings are the research’s main contributions: (1) This paper, with the best of my knowledge, is the first three-stage SBM-DEA model to be applied to discuss energy consumption and unfavorable CO2 emissions on sustainable development from the various violated aspects, like industrial structure, energy structure etc. (2) When applying the model, previous researches focus on provincial or industrial analysis, while this study focuses on city-level evaluation, especially on GBA. (3) According to the UNEP Governing Council’s definition of sustainable development, this research applies Hong Kong SAR and Macao SAR as decision-making units (DMUs) and then compare with other mainland cities. It is very important view of this research.

2. Literature Review

With the growing concern about adverse climate conditions and extreme weather due to global warming, a lot of efforts have been made to maintain the environmental sustainability. The GBA is considered one of the most economically dynamic regions in China since the reform and opening up. The development of GBA is of strategic importance in the development of China. In the early stages of development, most cities in the GBA achieved economic growth through labor- and energy-intensive industries, leading to over-exploitation of resources and environmental pollution. However, with the introduction of green and energy reforms in the Guangdong-Hong Kong-Macao Greater Bay Area in recent years, especially since the area consists of the cities of Guangdong as well as Hong Kong and Macao SARs, the study of the region has become representative. In February 2019, the state promulgated the Outline of the Development Plan of the Guangdong-Hong Kong-Macao Greater Bay Area, which states that priority is given to resource conservation, environmental protection, and natural restoration [3]. Therefore, how the government can improve the ecological environment and protect existing natural resources through various policies has become an important issue. Among them, controlling energy consumption is one of the keys to sustainable development, and this paper hopes to improve the energy efficiency of cities in the Guangdong-Hong Kong-Macao Greater Bay Area through policy recommendations, so as to achieve green and low-carbon development in the Guangdong-Hong Kong-Macao Greater Bay Area.

2.1. Sustainable Development

The IUCN initially established the notion of sustainable development in the World Conservation Strategy in 1980. This document states, “Sustainable development emphasizes the human use of the biosphere to manage it in order to meet the maximum sustainable interests of current generation while preserving its potential for the needs and desires of future generations” [4]. In 1989, the Governing Council of the United Nations Environment Program provided a more specific definition of sustainable development in its Statement on Sustainable Development: “It is defined as development that meets current demands without jeopardizing future generations’ ability to meet their own needs, and it does not imply any violation of national sovereignty. According to the UNEP Governing Council, achieving sustainable development involves both domestic and international cooperation, include providing aid to developing nations in line with their national development programs’ aims and development objectives. Furthermore, sustainable development entails a favorable international economic environment that promotes long-term economic growth and development in all countries, particularly in developing countries, which is critical for good environmental management. The maintenance, rational use, and enhancement of the natural resource base that promotes ecological stress tolerance and economic growth is also part of sustainable development. Furthermore, sustainable development refers to the integration of environmental concerns and considerations into development plans and policies, rather than a new kind of aid or development financing conditionality” [5].
In the above specific definition of sustainable development, it can be found that it encompasses various elements such as natural resources, economic growth, and the environment, which can be summarized as ecological development, economic development, and social development. Because of its central position in the sustainable development system, economic sustainable development has its unique definition. The British economist Barbier defines sustainable development from an economic perspective as “the maximization of net economic benefits while preserving natural resource quality and service provision” [6]. The British economist Pierce defines it as “economic development on the premise that natural capital remains unchanged, or that the use of resources today should not reduce real income in the future” [7]. The difference between the two is whether economic development is at the expense of resources and the environment. Barbier believes in maximizing economic benefits while preserving natural resource quality and using the money gained from environmental pollution and ecological damage as compensation for environmental and ecological construction. Pierce, on the other hand, believes that economic development should be done without destroying the world’s natural resource base, and should not destroy before repairing. Yang (2002) defines sustainable economic development as “the continuous improvement of the economic welfare of the present generation on the basis of certain resources and environment, while being able to ensure that the economic welfare received by future generations is not less than the economic welfare enjoyed by the present generation” [8]. In addition to emphasizing the protection of natural resources, this definition also proposes another provision: emphasizing intergenerational equality.
The Greater Bay Area is an excellent object study for studying China’s sustainable growth. The GBA includes the nine cities of Guangzhou, Shenzhen, Foshan, Zhuhai, Dongguan, Huizhou, Zhongshan, Jiangmen and Zhaoqing of the Guangdong Province, the Hong Kong Special Administrative Region (Hong Kong SAR) and the Macao Special Administrative Region (Macao SAR). GBA is considered as one of the most economically vibrant regions in China since the reformation. With a population of around 86 million, a total area of 56,000 km2 and a GDP of USD 1669 billion in 2020 [9], the development of GBA has also its strategic importance in the development of China. In the early stage of development, most GBA cities achieved economic growth through labor and energy-intensive industries, resulting in over-exploitation of resources and pollution of the environment. However, with the introduction of green and energy reform in the GBA recently, especially with the GBA comprising of Guangdong municipalities as well as the Macao SAR and the Hong Kong SAR, the study on this area becomes representative. According to the February 2019 “Guangdong-Hong Kong-Macao Greater Bay Area Outline Development Plan” [3], priority is given to resource conservation, environmental protection and restoration of nature. Thus, efforts should be made to improve ecological and environmental quality, and then a greener and low-carbon development of GBA can be achieved.
In China, which is in the middle of industrialization, the most important factors hindering the implementation of sustainable development strategy are the significant share of high-energy-consumption industries, the unbalanced industrial structure, the backward technology and equipment, the large population base, the large number of poor people, and the generally low living standard of the people. Increasing industrial restructuring efforts, accelerating technology and equipment replacement, and actively advocating energy-saving production and lifestyle are the primary solutions to promote sustainable development in China.
Then, it is essential to do research on energy consumption efficiency for sustainable development, for the 11 cities of the Guangdong-Hong Kong-Macao Greater Bay Area (GBA) during 2010–2016, by going through various violated factors including investment, human capital, and environmental protection. Especially, compared with the high economic growth rate at the current stage, the low growth rate might lower the consumption of energy resources and other environmental aspects of nature, which play a catalytic role in the global process of sustainable development [10].

2.2. The Relationship between Energy Consumption, Environmental Protection, and Sustainable Development

Early Western mainstream economists fully recognized the role of capital, technological progress, labor, and institutions on economic growth, but did not pay attention to the impact of natural resources and environmental factors on economic growth. During the oil crisis in the 1970s, when economies suffered from economic growth difficulties of different degrees due to the scarcity of energy resources, natural resources and environmental factors began to be noticed by economists. The concept of sustainable development was introduced in the 1980s, and natural resources and environmental factors became important indicators for the analysis of sustainable development [11].

2.2.1. The Relationship between Sustainable Development and Energy Consumption

Energy is the material basis of social production and human survival, and is inextricably linked to the development of a country’s economy. Sustainable development emphasizes the integrated and coordinated development of economy, resources, environment and population. Therefore, sustainable development of the energy economy means that the use and development of energy meet the needs of economic development without causing serious or irreversible damage to the ecosystem. It addresses the current generations’ energy needs without posing a danger to the future generations’ ability to satisfy their demands.
The original concept of a sustainable energy economy dates back to 1972, when the book Limits to Growth was published, which introduced the concept that social development and economic growth cannot be supported without energy, causing a heated controversy at that time [12]. Rashe and Tatom (1977) first introduced energy into the Douglas production function, hoping to uncover a fundamental pattern between economic growth and energy consumption [13]. Stem (1993) used a vector autoregressive (VAR) pattern with four variants (capital, energy consumption, labor, and GDP) and applied Granger multivariate causality tests to find a one-way Granger causality from energy consumption to GDP in the United States for 1947–1990 [14]. Willem et al., (2010) presented an explicit model of energy economic growth using a historical deductive approach in 2010, concluding that the current socio-economic model is unsustainable and that energy security is critical [15]. Zhang et al., (2011) examined the relationship between economic growth and energy consumption using the error correction model, cointegration test, and Granger causality test. They found that there is a significant one-way Granger causality of economic growth on energy consumption and that the two are in a long-term equilibrium relationship. [16]. Li et al., (2018) constructed a VAR model of the energy-economy-environment in Henan Province based on the statistical data of total energy consumption, GDP and industrial emissions from 2000–2014. Through impulse response function and variance decomposition, a long-term stable cointegration relationship between energy consumption, economic growth and environmental pollution in Henan Province is obtained, while energy consumption plays a certain positive impact on economic growth [17]. He et al., (2018) analyzed the relationship between energy consumption and economic growth in China since the 1950s by means of the elastic decoupling index and the generalized LMDI method and concluded that both total energy consumption and GDP growth exhibit an exponential growth curve [18]. Zhang and Wang (2018) examined the spatial variation of total factor energy efficiency in 30 provincial and municipal regions of China from 2006 to 2015 based on the haze constraint and concluded that economic development has significantly contributed to the improvement of energy efficiency [19]. Wang and Li (2019) studied Chinese coal cities and found a significant positive correlation between coal resources and economic growth in the sample cities [20]. Previous literature suggests a close relationship between energy consumption and sustainable economic development.

2.2.2. The Relationship between Sustainable Development and Environmental Protection

Environmental protection is the general term for human actions to ensure sustainable social and economic development and to solve real or potential environmental problems. Environmental sustainability is the goal of sustainable development that satisfies future generations’ demands by protecting the environment and reducing ecological burdens through measures such as reducing pollutant emissions.
The most representative early study of the relationship between economic sustainability and environmental protection was the empirical verification of the environmental Kuznets curve by American economists Grossman and Kruger in 1990, which demonstrated the relationship between environmental pollution and economic growth is inverted U-shaped [21]. Subsequently, more studies demonstrated this curve. Shafik et al., (1992) found that the water quality of rivers in cities was deteriorating because of the fast growth of the economy; the concentration of suspended solids and sulfur dioxide in the atmosphere in cities showed an inverted U-shaped relationship with per capita income level; and sulfur dioxide emissions in cities increased with per capita income climbed [22]. Sherry and David (2008) studied the inverted U-shaped relationship between environmental pollution and economic growth from the perspective of theoretical analysis [23]. Salvador et al., (2010) used differential dynamics models to simulate the association between economic growth, CO2 emissions and population, etc. [24]. Liang, et al., (2011) calculated the ecological environment index of 31 Chinese provinces and cities using entropy weight method and fuzzy comprehensive evaluation method, and analyzed the ecological environment status of different locations in the Yangtze River region, and carried out a study on the ecological environment sustainability with the Yangtze River basin as an example [25]. Feng et al., (2018) established a VAR model to analyze qualitatively and quantitatively the factors influencing carbon emissions in Beijing during 1996–2016, and the study found that carbon emissions per capita showed a positive relationship with GDP per capita [26]. Wang et al., (2018) studied the relationship between carbon emissions and GDP per capita in 27 provinces in China based on EKC theory. It was found that the relationship between carbon dioxide emissions and economic development of the 27 provinces as a whole, or geographically divided into eastern, central, and western provinces, met the inverted “U” curve of the EKC hypothesis [27]. According to Xu et al., (2019), green development is an expression of sustainable development, which is essentially a balanced relationship between the environment and the economy to promote economic development and ecological protection, so that the two complement each other [28]. Previous literature shows a close link between environmental protection and sustainable economic development, and it is mostly presented in an inverted U-shape.

2.2.3. Other Impacts on Sustainable Development

In addition to energy consumption and environmental protection factors, there are many external environmental factors that affect sustainable development, such as: industrial structure, energy structure, government financial expenditure, etc.
Tian (2007) proposed that the restructuring of industry forms a structural consumption of resources and thus leads to environmental changes, which largely affect the sustainable development of the economy [29]. Based on data from Shanxi Province, Wang and Gao (2018) used a double difference model to conclude that the increase in the degree of fiscal decentralization is conducive to the upgrading of industrial structure and thus has a positive impact on the ecological environment [30]. Yang (2000) analyzed the causal link between consumption of various energy sources and GDP (electricity, natural gas, oil, and coal) using the cointegration technique using sample data from 1954–1997 in Taiwan and found that there is a two-way Granger causality between GDP and coal and electricity and a one-way Granger causality from GDP to oil and natural gas to GDP [31]. Tian (2014) found that the fiscal expenditure policy has an indirect effect through influencing the economy and thus the ecological environment is significantly based on the analysis of fiscal revenue and expenditure data of Chinese provinces from 2008 to 2011 [32]. Jiang (2018) based on the analysis of data from 2007–2015 concluded that environmental spending has both economic and environmental nature and has a significant positive correlation with economic growth [33]. Hoff and Stiglitz (2001) argue that fiscal expenditure policies should focus on the updating and application of pollution treatment technologies to reduce pollution emission levels through technological innovation in order to enhance ecological protection and promote economic development [34]. Shi (2002) in “Recommendations on China’s energy consumption decreased” proposed that opening up to the outside world has a significant role in improving the efficiency of energy use [35]. Xu (2004) selected Shandong Province as a sample to analyze the supply and demand of human capital in Shandong Province and made a qualitative and quantitative analysis study respectively calculated that the support capacity of human capital to sustainable development reached 66.65% [36].

2.3. SBM-DEA Model

Data Envelopment Analysis (DEA) is a non-parametric approach applied to obtain the efficiency of multiple Decision-Making Units (DMUs). It is like a peer group comparison using a frontier to determine efficient and inefficient units relatively. DEA has been widely used as a method in DMUs’ efficiency analysis since the production function is not necessary for evaluation.
DEA was first introduced by Charnes, Cooper and Rhodes in 1978, and numerous researches have been done on methodological developments, practical applications, the status of variables and data variation, etc [37]. Since the first issued paper on using DEA in energy efficiency by Färe, Grosskopf and Logan (1983) on the relative efficiency of Illinois electric utilities, there are more and more researches on energy efficiency evaluation using DEA thereafter [38]. Bian and Yang (2010) used DEA models for estimating the aggregated efficiency of resources and the environment of the 30 provinces in China [39]. Guo et al., (2017) studied the dynamic DEA model to evaluate efficiencies based on fossil-fuel CO2 emissions in OECD countries and China [40]. Iftikhar et al., (2018) applied the network DEA model under the free disability assumption for all undesirable outputs in 19 major economies to assess the energy and CO2 emissions efficiency [41]. Mardani et al., (2017) also conducted an extensive review which indicated that DEA has shown to be a great evaluation tool for future analysis on energy efficiency issues, indicating that DEA is a common tool for energy and environmental efficiency researches [42].
Traditional DEA models include the CCR model and BBC model that deal with multiple inputs and multiple outputs to calculate efficiency. However, with the consideration of undesirable outputs such as atmospheric pollutants such as CO2 emission, traditional DEA models do not perform well in increasing desirable outputs while decreasing undesirable outputs at the same time. Thus, Tone (2001) proposed a non-radial, non-angled slacks-based measure (SBM) model, which can improve radial models like the CCR model and BBC model that do not consider the slacks of inputs and outputs [43]. Du et al., (2016) constructed an SBM-DEA model to evaluate the total factor energy efficiency of 29 provincial-administrative regions of China during 1997–2011 and the influential factors [44]. Huang and Wang (2017) studied the total-factor energy efficiency of 276 cities in China during 2000–2012 by using a three-stage SBM model, with consideration of influential factors and undesirable output [45]. Shang et al., (2020) conducted a study using the SBM-DEA model to calculate the total factor energy efficiency including 30 provinces and municipalities of China from 2005 to 2016 [46]. Through previous literature, it is shown that SBM-DEA Model is a commonly used tool for assessing energy efficiency and can bring effective calculations for the energy efficiency assessment of DBA in this paper.

3. Methodological Framework

Various methods can be used to examine the energy efficiency regarding sustainability. One of the most direct ways is to evaluate the economic output in terms of energy input, i.e., GDP/energy consumption. In this study, however, a production possibility frontier theory is chosen as the method to examine energy efficiency. Instead of using only one input and one output to calculate efficiency, a nonparametric approach (DEA) is used, which has the advantage of measuring the relative energy efficiency of DMUs with multiple inputs and outputs. In addition, by application of a DEA method for efficiency analysis, the influences of external factors and stochastic disturbances can be removed.
Charles, Cooper and Rhodes are famous scholars in operation research who have developed a new systematic analysis approach known as Data Envelopment Analysis (DEA). During these years, different elaborations on DEA have been carried out by numerous researchers to tackle different situations and scenarios.
In a DEA model, each decision-making unit (DMU) is assumed to contain inputs xRn, ygRS1 desirable inputs and ybRS2 undesirable outputs. The production possibility sets are thus defined as follows:
P = x , y g ,   y b           | x X λ ,   y g Y g λ ,   y b Y b λ , λ 0
which λ R n is the vector of intensity.
Based on the above production possibility set, a classic DEA model, known as a CCR model, can be obtained as follows:
θ * = min θ ε i = 1 m S i + r = 1 s S r + s . t .   j = 1 n X i j λ j + S i = θ X i j   ,         i = 1 , 2 , , m j = 1 n Y r j λ j S r + = Y r j   ,     r = 1 , 2 , , s         θ , λ , S i , S r + 0   ;             j = 1 , 2 , ,   n
where S i is the excess variable and S r + is the insufficient variable respectively, making effective frontiers to expand horizontally or vertically to form the envelope. The variable θ represents the efficiency of the DMU.
Besides CCR, BCC model is another classic DEA model and is described as follow:
min θ 0   s . t .   j = 1 n λ j x r j θ 0 x i 0 ,             i = 1 , 2 , , m j = 1 n λ j y r j y r j ,                                   r = 1 , 2 , , s                               j = 1 n j = 1   j 0 ,                                             j = 1 , 2 , , n
BBC model refers to pure technical efficiency, whereas CCR model refers to both scale efficiency and technical efficiency.
In this research, a three-stage DEA method is proposed to calculate energy consumption efficiency on sustainable development. The idea is to determine the comprehensive energy efficiency in phase 1, then eliminate the influence of stochastic disturbances and external environmental influences in the second stage, so as to enable the real efficiency to be evaluated in the third stage. The three-stage DEA methodology framework for the efficiency in the GBA cities is described in Figure 2.

3.1. The Initial Phase DEA: The Undesirable-SBM Pattern with Original Inputs

In comparison to CCR and BCC models, the undesirable-SBM model deals with input excess and output shortfall, according to Tone (2001) [43]. The formula is written as:
ρ * = min 1 1 m i = 1 m S i s i 0 1 + 1 S 1 + S 2 r = 1 S 1 S r g y r 0 g + r = 1 S 2 S r b y r 0 b
s . t .   x 0 = X λ + S y 0 g = Y g λ S g y 0 b = Y b λ S b λ 0 ,   S 0 ,   S g 0 ,   S b 0
From the above definitions, each DMU has m inputs, S1 desirable outputs and S2 undesirable outputs. The slacks of inputs, desirable outputs and undesirable outputs are S, Sg and Sb respectively. When a DMU is efficient, S = 0, Sg = 0 and Sb = 0. If either S, Sg and Sb is not equal to 0, then the objective function ρ* is not 1, 0 ≤ ρ* <1, and this DMU is described as inefficient.

3.2. The Second Phase: Frontier Analysis with a Random Component

DEA is an ideal way to assess efficiency without having to know the production functions. However, as Fried et al., (2002) suggested, stochastic disturbances, environmental impacts, and managerial efficiency can all affect the efficiency calculated by DEA [47]. In order to deal with this disadvantage and to evaluate the real efficiency of each DMU, a SFA model is established to decompose the slacks in the first stage due to measurement errors in the input and output variables. To describe a SFA model, the environmental variables can be used as explanatory variables and the relaxation values in the first stage as interpreted variables:
S n i = f n Z i ; β n + ν n i + μ n i                                       n = 1 , 2 , N ,         i = 1 , 2 I
In Equation (4), Sni stands for the relaxation value of the n-th input of the i-th DMU and f n Z i ; β n indicates the impact of the environmental variables on the relaxation value of the input and Z i = z 1 i , z 2 i , z k i ,       i = 1 , 2 , I , and k are the external environmental variables. βn is the coefficient of environmental variables. The term ν n i + μ n i   is the mixed error, in which ν n i is a random error term and ν n i ~ N 0 , σ ν 2 ; μ ni is the management efficiency and μ n i ~ N + 0 , σ μ 2 . The term ν n i and μ n i are independent and does not relate to each other. When γ = σ μ i 2 σ μ i 2 + σ ν i 2 is close to 1, it means that the difference in efficiency is affected by the managerial efficiency. However, when γ = σ μ i 2 σ μ i 2 + σ ν i 2 is near to 0, it suggests that efficiency remains primarily untouched by stochastic disruption.
The mixed error term can be further decomposed as follow [28]:
E ν i k / ν i k + μ i k = s i k f i z k ; β i E μ i k i k + μ i k
Then, by applying the Dengyue’s approach, E μ i k | ν i k + μ i k can be solved as following [48]:
E μ i k / ν i k + μ i k = σ λ 1 + λ 2 φ ε k λ σ ε k λ σ + ε k λ σ
with λ = σ u σ v ,   ε k = ν i k + μ i k ,   σ 2 = σ u 2 + σ v 2 ,   φ ,   are the standard normal distribution’s density and distribution functions respectively.
In order to adjust the input variables for real efficiency calculation under the same external environment, the SFA model is used to identify and exclude the associated environmental variables and stochastic disturbances. The formula for adjusted input variables X n i * is then described as:
X n i * = X n i + max i Z i β n Z i β n + max i ν n i ν n i
where n = 1 , 2 , N   i = 1 , 2 , I ; X ni * denotes the adjusted input after the SFA model while X ni indicates the initial input variable of the first phase; max i Z i β n Z i β n describes the external environmental variable adjustment and max i ν n i ν n i represents the adjusted stochastic disturbance value.

3.3. The Third Phase: The Undesirable-SBM Pattern with Adjusted Input Variables

During phase 3, the efficiency is once again evaluated using undesirable-SBM model but replacing the original inputs with the adjusted input variables X n i * instead. By comparing with the results of the first phase, the impact of stochastic disruption and external environmental variables have been removed as the outcomes of the third phase, indicating a real efficiency of the GBA cities.

4. Variables and Data Sources

The SBM-DEA Model’s Input and Output Variables

During the examination of energy efficiency by using the SBM-DEA model, production assumptions are not necessary and multiple input and output variables can be handled at the same time. This is beneficial in cases when examining energy efficiency when details are uncertain. Table 1 shows a selection of input and output variables commonly used in DEA models.
Based on previous studies, the input variables can be justified by the total number of employees, the total investment in fixed assets and total energy consumption in this study.
Regarding output variables, both desirable output and undesirable output are included, especially with the global agenda of low carbon development. Therefore, both the environmental influence and the economic output should be considered. Hence, GDP is chosen as the desired output, while the undesirable output is CO2 emissions based on Table 1. Table 2 presents the input and output variables that are chosen, while Table 3 demonstrates the relevant statistical descriptions.
For differentiation, if the chosen input and output variables are appropriate, the Spearman’s Rank correlation coefficients are calculated to demonstrate the strength of the link between inputs and outputs. The relevant outcomes in Table 4 show that the coefficients between input and output variables possess a positive correlation with significant values of 1% and 5%. Thus, the input and output variables have been chosen appropriately.
External environmental variables do also exert influence on efficiency as mentioned by Liu et al. [49]. In terms of selecting these environmental variables, based on existing literatures, five major factors are analyzed in this research, which include industrial structure, energy structure, technology, government intervention and openness of economy [44,50,51]. In the SFA regression, Table 5 lists the factors chosen as explanatory variables.
Secondary industries are known to be more energy intensive industries than tertiary industries. Thus, the industrial composition of a city exerts huge difference in energy demand and its energy efficiency. Besides, the types of energy resources used also affect the performance of outputs too. Coal is known to be an inexpensive but unclean source of energy. CO2 emission is one of the environmental concerns when coal is being combusted tremendously during development. Thus, advancement in technology is necessary to clean up air pollution or even to utilize clean energy as a replacement. In this regard, the government can intervene with fiscal incentives to incur research and development in technology. Besides, the frequency of import and export often remarks the openness of an economy, hopefully with more technology transfer of goods and services taking place at the same time.
The data used in this research is mainly comprised of variables of 11 cities, including 9 Guangdong cities (Guangzhou, Shenzhen, Foshan, Zhuhai, Dongguan, Huizhou, Zhongshan, Jiangmen and Zhaoqing), and 2 special administrative regions, namely the Hong Kong Special Administrative Region (Hong Kong SAR) and Macao Special Administrative Region (Macao SAR) from 2010 to 2016. In most of the previous literatures on sustainable development of the cities of China, data of HK SAR and Macao SAR are often excluded due to statistical caliber. Therefore, in order to fulfill the consistency of input and output variables as well as for all cities in GBA, all data were derived from the official Statistical Yearbook of each Guangdong municipality, the Census and Statistics Department of the Hong Kong SAR and the Statistics and Census Service of the Macao SAR from 2010 to 2016. Since there is no official CO2 emission data of the GBA cities available, relevant CO2 data from 2010 to 2016 was directly obtained from research by Zhou et al. [52]. According to the study of Zhou et al., (2018), the total CO2 emissions is calculated based on territorial emission accounting approach of the Intergovernmental Panel on Climate Change (IPCC), with inventories consisting of 17 kinds of fossil fuels, 47 socio-economic sectors and 7 industrial processes.

5. Results and Discussion

5.1. Results

5.1.1. The First Phase Undesirable-SBM Model: The Comprehensive Efficiency Calculation Results

During the first stage, with consideration of the undesirable output CO2 emission, the undesirable-SBM pattern stated in article 3.1 was chosen to calculate the comprehensive efficiency of the 11 GBA cities from 2010–2016. The calculation was done with MaxDEA (Beijing Realword Software Company Ltd, Beijing, China) and the outcomes were displayed in Table 6.

5.1.2. The Second Phase: The Analysis of Influence of External Environmental Factors on Efficiency

During phase 2, the stochastic frontier analysis (SFA) approach is applied to analyze the influence of exterior environmental elements on the relaxation variables of inputs, namely the relaxations of the total number of employees, the slacks of total investment of fixed assets and the slacks of the total energy consumption. External environmental variables are elements that will affect energy efficiency, but not in a direct and controllable way. Five factors include the ratio of secondary industry GDP to total GDP, the ratio of coal consumption to the total energy consumption, the ratio of R&D investment to total GDP, the ratio of local fiscal expenditure to total GDP and the ratio of import and export to total GDP. The explanatory factors are these external environmental variables, while the slack inputs are considered as explained variables. Frontier 4.1 software is being used to create the SFA model. Table 7 shows the empirical findings.
After conducting the second stage SFA regression model, the input data in the first stage is adjusted with the coefficient values as shown in Table 7. With this adjustment, the effects from stochastic disturbances and external environmental elements can be eliminated. Adjusted input and original output variables can be run by the undesirable SBM-model once again to obtain the real energy efficiencies. The adjusted results of energy efficiency are compared with that of the original shown in Figure 3.

5.2. Discussion

As seen in Table 5, Shenzhen, Macao SAR and Hong Kong SAR are operating at the efficient frontier from 2010–2016, with energy efficiencies equal to 1. This implies that Shenzhen, Macao SAR and Hong Kong SAR are efficient in terms of energy consumption. However, other cities such as Guangzhou, Zhuhai, Foshan, Jiangmen, Zhaoqing, Huizhao, Dongguan, and Zhongshan, their efficiencies are all below the average value of 0.494 from 2010–2016. Among inefficient cities, only Guangzhou and Dongguan are close to the average value, but still, the tendency of comprehensive efficiencies is not promising and have not been increasing obviously from 2010–2016. On the contrary, the comprehensive efficiencies of Shenzhen, Hong Kong SAR and Macao SAR demonstrate better performance in efficiencies and these cities have similar economic activities, industrial structure, technology and market openness.
According to previous literature demonstrating that there is a clear one-way Granger causality between sustainable economic development and energy consumption, the energy efficiency of the eight cities Guangzhou, Zhuhai, Foshan, Jiangmen, Zhaoqing, Huizhou, Dongguan and Zhongshan has been below the average from 2010–2016, and as the cities develop year by year, the low energy efficiency will have an impact on future sustainable economic development, and an efficient and green Energy consumption policy is an important issue at this stage to ensure sustainable development in the future.
From the SFA results, certain evaluated coefficients listed in Table 6 are significant at the 1%, 5%, and 10% significance levels respectively, indicating that certain external environmental conditions do have an impact on the GBA cities’ energy efficiency. Besides, the correlation coefficients calculated in the SFA regression also represent the relationship between input relaxations and external environmental conditions.
For the GDP of secondary industries as a percentage of overall GDP, it shows a positive and 1% significant correlation with the slack of investment in fixed assets. An increase in the ratio of secondary industry GDP to the total GDP will increase the slack of the investment in fixed assets, which in turn bringing unfavorable impacts on energy efficiencies. Secondary industry is known to be an energy intensive industrial category and an adjustment policy in industrial structure to less energy intensive industry may have a positive effect on energy efficiencies of the GBA cities.
Coal combustion in energy generation is the known source of CO2 emission. Thus, in Table 6, it shows that when the ratio of coal depletion to overall energy consumption rises, the slacks of energy consumption rise as well. However, the result is not significant.
R&D investment as a percentage of GDP shows positive and significant correlation coefficients with the slacks of number of employees and energy consumption. This finding is not similar to Zhao et al. [53], as investment of R&D usually encourages technology improvement and innovation that can contribute to efficiency in production process, causing a decrease in labor and investment demand, as well as more energy saving technologies. However, Liu (2020) mentioned that technological progress affects negatively on energy efficiency in her research [54]. One of the possible reasons maybe due to rebound effect, indicating that an increase in energy efficiency triggered by technology advancement is being offset by the increase in energy consumption demand.
The ratio of local fiscal expenditure to total GDP is negatively but not significantly correlated with the input slacks. An increase in government involvement in economic activities may be interpreted as incentives or drivers in energy efficiency improvement such as energy savings technology or implementation of relevant efficiency policies. The significant relationship between the ratio of local fiscal expenditure to GDP and investment in fixed assets signals the government’s part in driving energy efficiency performance.
For the empirical study, the intensity of import and export activities shows a negative relationship with all the three input variables. Any increase in import and export as a percentage of GDP will decrease the slacks of inputs. Among all other environmental factors, import and export as a percentage of GDP appears to be the only factor that is significant to all the input variables of the GBA cities. It is worthwhile to note that the average import and export to GDP value of the GBA cities has been declining from 1.38 in 2010 to 0.97 in 2016 and the average energy efficiency has also dropped since 2010.
The average energy efficiency of the third stage during 2010–2016 has increased from 0.494 to 0.708, compared with the first phase. This explains that the overall energy efficiencies may have been underestimated in the first stage and are affected by the external environmental factors and stochastic disturbances. During phase 3 calculation, the average energy efficiencies of the GBA cities range between 0.604 and 0.851 throughout the study period from 2010–2016. From the first phase calculation shown in Figure 4, the values of energy efficiencies have increased but the average energy efficiency of the GBA cities does not improve during the period 2010–2016. In fact, the average energy efficiency had reached a peak of 0.85 in 2011 and since then, it started to drop to 0.605 in 2016. By eliminating the external environmental factors and stochastic disturbances after the third stage, a decrease in energy efficiency implies a real decline in managerial efficiency of the GBA cities since 2011.
From Figure 5, it is shown that Shenzhen, Hong Kong SAR and Macao SAR are energy efficient cities, with an energy efficiency value of 1. When stochastic disturbances and external environmental elements are excluded after the SFA regression, the number of energy efficient cities varied from 4 to 7 in the studied years, including Shenzhen, Zhuhai, Zhaoqing, Zhongshan, Dongguan, Macao SAR and Hong Kong SAR. However, the number of cities with an average energy efficiency of 1 both in the first and the third stage of the SBM-DEA throughout 2010–2016 remain only 3, namely Shenzhen, Hong Kong SAR and Macao SAR.

6. Conclusions and Policy Recommendations

6.1. Conclusions

The following conclusion can be obtained using a three-stage SBM-DEA model to determine energy efficiencies and investigate the influential elements of the GBA cities:
(1)
Among the 11 GBA cities, only Shenzhen, Hong Kong SAR and Macao SAR have an energy efficiency of 1 from 2010–2016, both in the initial phase as well as in the third phase. This means that all these three cities operated at the efficient frontier during the study period. Other cities such as Guangzhou, Zhuhai, Foshan, Jiangmen, Zhaoqing, Huizhao, Dongguan, and Zhongshan, their energy efficiencies were all below the average value of 0.494 (the first stage) from 2010–2016;
(2)
By eliminating the external environmental factors and stochastic disturbances by the stochastic frontier analysis during phase 2, the GBA cities’ average energy efficiency during 2010–2016 has increased from 0.494 to 0.708 during phase 3 SBM-DEA model. This explains that energy efficiency has been underestimated at the first stage. Besides, Shenzhen, Hong Kong SAR and Macao SAR, Dongguan was the other GBA city that had an average energy efficiency of 0.87, which is above the average of all GBA cities. At the same time, a decreasing trend of energy efficiency after the third phase SBM-DEA pattern since 2010 had been observed. This implies a real decline in managerial energy efficiency of the GBA cities since 2010;
(3)
Through the second stage stochastic frontier analysis, the influence of external environmental factors is investigated and among all, the ratio of import and export to GDP shows a negative and significant relationship to all three input variables, meaning an increase in the ratio will cause a decline in the input slacks, which favors energy efficiency.

6.2. Policy Recommendations

Based on the empirical analysis and conclusion, several policy recommendations are proposed for energy efficiency and sustainable growth:
(1)
Acceleration of economic structure transformation. Instead of achieving the sole target of GDP growth, more emphasis should be made on environmental protection such as increasing the ratio of clean fuels and minimizing the emission of atmospheric pollutants such as CO2 emission, especially under the agenda of sustainable development.
(2)
Financial support from the government. Government intervention is always an important driver in encouraging energy-efficient and low-carbon production, especially during the initial stage in fixed assets investment. Government can also promote public awareness and enforce tougher environmental protection standards.
(3)
Promotion of imports and exports. When products are less energy-intensive and environmentally friendly, an increase in imports and exports will not only facilitate a higher GDP, but also indirectly boost technological innovation and living standards.
This study, however, has a few limitations. First, this is the first research of this kind with the Guangdong-Hong Kong-Macao Bay Area cities as the DMUs and using the SBM-DEA model for energy efficiency evaluation, as far as we know. Thus, it is recommended that other DEA methods should be studied to examine the energy efficiency for comparison of results. Second, the structure of statistical data of mainland China and that of HKSAR and Macao SAR are quite different and may cause uncertainties in results. Further development of a homogeneous GBA statistical data collection system would be helpful to minimize the associated uncertainties. Last, different GBA cities exhibit different energy efficiency performances and are affected by different influential factors. Investigation on other social or economic influential factors for a wider aspect of policy formulation would be helpful.

Author Contributions

Conceptualization, W.F. and Y.S.; methodology, W.F. and Y.S.; software, W.F. and Y.S.; validation, W.F.; formal analysis, W.F. and Y.C.; investigation, W.F.; resources, W.F.; data curation, W.F. and Y.C.; writing—original draft preparation, W.F. and Y.S.; writing—review and editing, W.F. and Y.S.; visualization, W.F. and Y.S.; supervision, Y.C.; project administration, Y.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Not available.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. World Commission on Environment and Development. Our Common Future; Jilin People’s Publishing House: Jilin, China, 1997; p. 52. [Google Scholar]
  2. Sihua, L. Some questions about sustainable development and sustainable economy. Contemp. Financ. 1997, 6. [Google Scholar]
  3. Outline Development Plan for the Guangdong-Hong Kong-Macao Greater Bay Area. 2018. Available online: https://www.bayarea.gov.hk/en/outline/plan.html (accessed on 18 February 2021).
  4. Guangyu, L.; Ziliang, S. Economy, Environment and Law; Science Publishing House: Karachi, Pakistan, 2000; p. 19. [Google Scholar]
  5. Qigui, Z. Sustainability Assessment; Shanghai University of Finance and Economics Publishing House: Shanghai, China, 1999; pp. 19–20. [Google Scholar]
  6. Barbier, E.B. Economics. Natural Resource Scarcity and Development; Earth Science: Hubei, China, 1989. [Google Scholar]
  7. Pearce, D.; Wofford, J. No End of the World—Economics—Environment and Sustainable Development; China Financial and Economic Publishing House: Beijing, China, 1996; pp. 49–69. [Google Scholar]
  8. Wenjin, Y. Economic Sustainability Theory; China Environmental Science Publishing House: Beijing, China, 2002; pp. 7–8. [Google Scholar]
  9. Constitutional and Mainland Affairs Bureau, H.K.S. Guangdong-Hong Kong-Macao Greater Bay Area. 2018. Available online: https://www.bayarea.gov.hk/en/about/overview.html (accessed on 18 February 2021).
  10. Lei, S.; Lituo, L.; Limao, W.; Fengnan, C.; Chao, Z.; Ming, S.; Shuai, Z. Scenario projection of China’s energy consumption in 2050. Nat. Resour. J. 2015, 30, 361–373. [Google Scholar]
  11. Lifeng, Z. Research on China’s Energy Supply and Demand Forecasting Model and Development Countermeasures. Ph.D. Thesis, School of Economics, Capital University of Economics and Trade, Beijing, China, 2006. [Google Scholar]
  12. Meadows, D.; Randers, J.; Meadows, D. Limits to Growth; Machinery Industry Publishing House: Shanghai, China, 2013. [Google Scholar]
  13. Rashe, R.; Tatom, J. Energy resources and potential GNP. Fed. Resevre Bank StLouis Rev. 1977, 59, 68–76. [Google Scholar] [CrossRef] [Green Version]
  14. Stern, D.I. Energy use and economic growth in the USA: A multivariate approach. Energy Econ. 1993, 15, 137–150. [Google Scholar] [CrossRef]
  15. Nel, W.P.; van Zyl, G. Defining limits: Energy constrained economic growth. Appl. Energy 2010, 87, 168–177. [Google Scholar] [CrossRef]
  16. Xinxin, Z.; Guangbin, L.; Lu, C. A study on the relationship between energy consumption and economic growth in China based on Granger test. J. Shanxi Univ. Financ. Econ. 2011, 33, 26–27. [Google Scholar]
  17. Xiaofei, L.; Lichen, Z.; Kewen, L. Research on the dynamic relationship between energy—Economy—Environment in Henan Province based on VAR model. J. Henan Univ. Technol. 2018, 1, 48–54. [Google Scholar]
  18. He, Z.; Yang, Y.; Song, Z.; Liu, Y. Mutual evolutionary trends and drivers of energy consumption and economic growth in China. Geogr. Res. 2018, 8, 1528–1540. [Google Scholar]
  19. Zhang, Z.W.; Wang, Z.L. A study on spatially divergent measures of energy efficiency under technological heterogeneity and haze constraints. East China Econ. Manag. 2018, 32, 65–74. [Google Scholar]
  20. Baoqian, W.; Jingya, L. An empirical study on the effect of “resource curse” in Chinese coal cities. Stat. Decis. Mak. 2019, 35, 121–125. [Google Scholar]
  21. Grossman, G.M.; Krueger, A. Economic growth and the environment. Q. J. Econ. 1995, 110, 353–377. [Google Scholar] [CrossRef] [Green Version]
  22. Shafik, N.; Bandayopadhyay, S. Economic Growth and Environmental Quality: Time Series and Cross Country Evidence; World Bank Publications: Washington, DC, USA, 1992. [Google Scholar]
  23. Bartz, S.; Kelly, D.L. Economic growth and the environment: Theory and facts. Resour. Energy Econ. 2008, 30, 115–149. [Google Scholar] [CrossRef] [Green Version]
  24. Puliafito, S.E.; Puliafito, J.; Grand, M. Modeling population dynamics and economic growth as competing species: An application to CO2 global emissions. Ecol. Econ. 2008, 65, 602–615. [Google Scholar] [CrossRef] [Green Version]
  25. Wei, L.; Konglai, Z. A study of ecological sustainability in the Yangtze river basin. Explor. Econ. Issues 2011, 8, 159–165. [Google Scholar]
  26. Feng, M.; Yang, S.; Zifu, Z. VAR model analysis of carbon emission influencing factors—Based on Beijing city data. Sci. Manag. Res. 2018, 36, 78–81. [Google Scholar]
  27. Wang, F.; Yang, X.; Yang, T. An empirical study on the relationship between carbon emissions and economic growth based on EKC hypothesis. Ecol. Econ. 2018, 34, 19–23. [Google Scholar]
  28. Xianchun, X.; Xue, R.; Zihao, C. Big data and green development. China Ind. Econ. 2019, 4, 5–22. [Google Scholar]
  29. Tao, T. A study on the relationship between industrial restructuring and the implementation of sustainable development strategy: The case of Hubei province. Mod. Shopp. Mall 2007, 02, 237. [Google Scholar]
  30. Wang, L.Y.; Gao, Y.R. Fiscal decentralization and industrial structure upgrading: Empirical evidence from a quasi-natural experiment of “directly administered counties”. Financ. Trade Econ. 2018, 39, 145–159. [Google Scholar]
  31. Yang, H.Y. A note of the causal relationship between energy and GDP in Taiwan. Energy Econ. 2000, 22, 309–317. [Google Scholar] [CrossRef]
  32. Dan, T. An Empirical Study of the Environmental Effects of Fiscal Balance Policy in China. Ph D Thesis, Wuhan University, Wuhan, China, 2014. [Google Scholar]
  33. Nan, J. Can fiscal spending on environmental protection help achieve a win-win situation for both the economy and the environment? J. Zhongnan Univ. Econ. Law 2018, 1, 95–103. [Google Scholar]
  34. Hoff, K.; Stiglitz, J.E. Modern Economic Theory and Development; World Bank: Washington, DC, USA, 2000. [Google Scholar]
  35. Dan, S. Recommendations on China’s energy consumption decreased. China Econ. Trade Her. 2002, 18, 26. [Google Scholar]
  36. Yinliang, X. Analysis of Human Capital Supply and Demand in Shandong Province and its Ability to Support Sustainable Development. Master ‘sThesis, Shandong Normal University, Jinan, China, April 2004. [Google Scholar]
  37. Cook, W.D.; Seiford, L.M. Data envelopment analysis (DEA)—Thirty years on. Eur. J. Oper. Res. 2009, 192, 1–17. [Google Scholar] [CrossRef]
  38. Färe, R.; Grosskopf, S.; Logan, J. The relative efficiency of Illinois electric utilities. Resour. Energy 1983, 5, 349–367. [Google Scholar] [CrossRef]
  39. Bian, Y.; Yang, F. Resource and environment efficiency analysis of provinces in China: A DEA approach based on Shannon’s entropy. Energy Policy 2010, 38, 1909–1917. [Google Scholar] [CrossRef]
  40. Guo, X.; Lu, C.C.; Lee, J.H.; Chiu, Y.H. Applying the dynamic DEA model to evaluate the energy efficiency of OECD countries and China. Energy 2017, 134, 392–399. [Google Scholar] [CrossRef]
  41. Iftikhar, Y.; Wang, Z.; Zhang, B.; Wang, B. Energy and CO2 emissions efficiency of major economies: A network DEA approach. Energy 2018, 147, 197–207. [Google Scholar] [CrossRef]
  42. 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]
  43. Tone, K. A slacks-based measure of efficiency in data envelopment analysis. Eur. J. Oper. Res. 2001, 130, 498–509. [Google Scholar] [CrossRef] [Green Version]
  44. Du, H.; Matisoff, D.C.; Wang, Y.; Liu, X. Understanding drivers of energy efficiency changes in China. Appl. Energy 2016, 184, 1196–1206. [Google Scholar] [CrossRef]
  45. Huang, H.; Wang, T. The total-factor energy efficiency of regions in China: Based on three-stage SBM model. Sustainability 2017, 9, 1664. [Google Scholar] [CrossRef] [Green Version]
  46. Shang, Y.; Liu, H.; Lv, Y. Total factor energy efficiency in regions of China: An empirical analysis on SBM-DEA model with undesired generation. J. King Saud Univ. Sci. 2020, 32, 1925–1931. [Google Scholar] [CrossRef]
  47. Fried, H.O.; Lovell, C.K.; Schmidt, S.S.; Yaisawarng, S. Accounting for environmental effects and statistical noise in data envelopment analysis. J. Product. Anal. 2002, 17, 157–174. [Google Scholar] [CrossRef]
  48. Denyue, R. Three stage DEA model management inefficiency estimation note. Stat. Res. 2012, 29, 104–107. [Google Scholar]
  49. Liu, H.; Zhang, Z.; Zhang, T.; Wang, L. Revisiting China’s provincial energy efficiency and its influencing factors. Energy 2020, 208, 118361. [Google Scholar] [CrossRef]
  50. Jie, H. Spatial network structure of energy and environmental efficiency in China and its influencing factors. Resour. Sci. 2018, 40, 14. [Google Scholar]
  51. Li, H.; Fang, K.; Yang, W.; Wang, D.; Hong, X. Regional environmental efficiency evaluation in China: Analysis based on the Super-SBM model with undesirable outputs. Math. Comput. Model. 2013, 58, 1018–1031. [Google Scholar] [CrossRef]
  52. Zhou, Y.; Shan, Y.; Liu, G.; Guan, D. Emissions and low-carbon development in Guangdong-Hong Kong-Macao Greater Bay Area cities and their surroundings. Appl. Energy 2018, 228, 1683–1692. [Google Scholar] [CrossRef] [Green Version]
  53. Zhao, H.; Guo, S.; Zhao, H. Provincial energy efficiency of China quantified by three-stage data envelopment analysis. Energy 2019, 166, 96–107. [Google Scholar] [CrossRef]
  54. Qingqing, L. Green Total Factor Energy Efficiency Measurements and Impact Factors; Jilin University: Jilin, China, 2020. [Google Scholar]
Figure 1. Total energy consumption, total CO2 emission and total GDP of the GBA cities from 2010–2016. (Data Source: International Energy Agency, 2021).
Figure 1. Total energy consumption, total CO2 emission and total GDP of the GBA cities from 2010–2016. (Data Source: International Energy Agency, 2021).
Energies 15 00569 g001
Figure 2. The three-stage DEA’s methodological framework.
Figure 2. The three-stage DEA’s methodological framework.
Energies 15 00569 g002
Figure 3. In the period 2010–2016, the average energy efficiency of the GBA cities’ first stage (original inputs) and third stage (adjusted inputs) was compared.
Figure 3. In the period 2010–2016, the average energy efficiency of the GBA cities’ first stage (original inputs) and third stage (adjusted inputs) was compared.
Energies 15 00569 g003
Figure 4. Energy efficiencies of the GBA cities in 2010–2016 (the first stage).
Figure 4. Energy efficiencies of the GBA cities in 2010–2016 (the first stage).
Energies 15 00569 g004
Figure 5. Energy efficiencies of the GBA cities in 2010–2016 (the third stage).
Figure 5. Energy efficiencies of the GBA cities in 2010–2016 (the third stage).
Energies 15 00569 g005
Table 1. List of selected input and output variables commonly used in energy efficiency DEA models.
Table 1. List of selected input and output variables commonly used in energy efficiency DEA models.
AuthorTitleInputOutput
Khalili-Damghani, K., Tavana, M., & Haji-Saami, E. (2015)A data envelopment analysis model with interval data and undesirable output for combined cycle power plant performance assessmentFossil fuelEnergy, CO2, SO2, SO3, NOx
Zha, Y., Zhao, L., & Bian, Y. (2016)Measuring regional efficiency of energy and carbon dioxide emissions in China: A chance constrained DEA approachLabor, capital, coal, oil, natural gasGDP, CO2
Zhao, L., Zha, Y., Liang, N., & Liang, L. (2016)Data envelopment analysis for unified efficiency evaluation: an assessment of regional industries in ChinaCapital input (Investment of fixed assets), labor, energyIndustrial production value, waste gas, wastewater
Chen, L., & Jia, G. (2017)Environmental efficiency analysis of China’s regional industry: a DEA based approachTotal number of employees of industry, total energy consumption, total investment in fixed assets of industryGDP, total emissions of SO2, industrial solid waste
Guo, X., Lu, C., Lee, J., & Chiu, Y. (2017)Applying the dynamic DEA model to evaluate the energy efficiency of OECD countries and ChinaLand area, population, and energy use, carry-over variable, energy stockCO2 emission and GDP
Huang, J., Du, D., & Hao, Yu. (2017)The driving forces of the change in China’s energy intensity: An empirical research using DEA-Malmquist and spatial panel estimationsEnergy, capital stock, labor, economic structure and GDPGDP
Sueyoshi, T., & Yuan, Y. (2017)Social sustainability measured by intermediate approach for DEA environmental assessment: Chinese regional planning for economic development and pollution preventionCapital, labor and energyGross regional product, CO2, SO2, soot (dust), wastewater, COD, ammonia nitrogen
Sueyoshi, T., Yuan, Y., Li, A., & Wang, D. (2017)Methodological comparison among radial, non-radial and intermediate approaches for DEA environmental assessmentCapital, labor, energyGross Regional Product, CO2, SO2, Smoke and dust, wastewater, COD, ammonia nitrogen
Zhao, L., Zha, Y., Wei, K., & Liang, L. (2017)A target-based method for energy savings and carbon emission reduction in China based on environmental data envelopment analysisLabor, capital stock, total energy consumptionGDP, CO2
Cayir Ervural, B., Zaim, S., & Delen D. (2018)A two-stage analytical approach to assess sustainable energy efficiencyTotal renewable energy potential, network length, total installed power of renewable energy, transformer capacityGross energy generation from renewable sources, number of consumers, total exports, GDP per capita, human development index, total energy production, population, area
Iftikhar, Y., Wang, Z., Zhang, B., & Wang, B. (2018)Energy and CO2 emissions efficiency of major economies: A network DEA approachLabor, capital, energy, populationCO2 emission, middle income class, high income class, low income class
Khoshroo, A., Izadikhah, M., & Emrouznejad, A. (2018)Improving energy efficiency considering reduction of CO2 emission of turnip production: A novel data envelopment analysis model with undesirable output approachAmount of energy for labor, machinery, diesel fuel, chemical fertilizers, seed and waterTurnip yield, GHG
Nadimi, R., & Tokimatsu, K. (2019)Evaluation of the energy system through data envelopment analysis: Assessment tool for Paris AgreementHeat, FEC, oil products, IMR, LEB, MYS, P_NRE, R_NRE, RE, IWACO2, electricity, GDP, GNI per capita, QoL
Piao, S., Li, J., & Ting, C. (2019)Assessing regional environmental efficiency in China with distinguishing weak and strong disposability of undesirable outputsLabor (number of employees), water (total domestic and industrial water consumption), energy and capitalGDP, CO2 (previous research), SO2, wastewater and waste
Wu, J., Li, M., Zhu, Q., Zhou, Z., & Liang, L. (2019)Energy and environmental efficiency measurement of China’s industrial sectors: A DEA model with non-homogeneous inputs and outputsCapital, labor, coal, oil, natural gasGross industrial output value, volatile hydroxyl-benzene, cyanide, COD, petroleum, ammonia-nitrogen, demand, petroleum, and ammonia-nitrogen
Yang, Z., & Wei, X. (2019)The measurement and influences of China’s urban total factor energy efficiency under environmental pollution: Based on the game cross-efficiency DEAProductive capital stockGDP, wastewater, SO2, smoke and dust, pollutants (comprehensive index)
Zhao, H., Guo, S., & Zhao, H. (2019)Provincial energy efficiency of China quantified by three-stage data envelopment analysisLabor force (total amount of employees at year end), capital (total investment on fixed assets), energy consumption (total energy consumption)GDP divided by amount of SO2 emission
Zhou, Z., Xu, G., Wang, C., & Wu, J. (2019)Modeling undesirable outputs with a DEA approach based on an exponential transformation: An application to measure the energy efficiency of Chinese industryEnergy, total assetTotal profit, industrial wastewater, industrial waste gas
Table 2. The choice of input and output variables.
Table 2. The choice of input and output variables.
VariableDefinition of VariablesUnits
InputsTotal number of employees10,000 people
Total investment in fixed assets100 million yuan
Total energy consumption10,000 tons of coal equivalent
Desirable outputGross domestic product (GDP)100 million yuan
Undesirable outputCarbon dioxide (CO2) emissionmillion tons
Table 3. The input and output variables, as well as the environmental elements, are statistically described.
Table 3. The input and output variables, as well as the environmental elements, are statistically described.
YearVariablesNo. of Employees (10,000 ppl)Investment on Fixed Assets (100 Million Yuan)Energy Consumption (10,000 Tons of Coal Equivalent)GDP
(100 Million Yuan)
CO2 Emission (Million Tons)Secondary Industry GDP to Total GDPCoal Consumption to Total Energy ConsumptionR&D Investment (10,000 Yuan) to GDPLocal Fiscal Expenditure (10,000 Yuan) to GDPImport and Export to GDP
2010Average359.171348.831641.304979.1630.450.390.310.010.111.38
Variance60,080.561,152,070.332,032,778.5723,002,949.16715.360.050.060.000.001.13
Maximum758.143263.574775.6014,928.9299.320.630.760.030.173.60
Minimum31.48230.1098.35965.122.130.050.000.000.060.23
2011Average373.251499.221716.275691.8632.850.380.360.020.111.33
Variance68,132.361,390,138.242,205,935.3428,291,300.48774.880.050.070.000.001.13
Maximum828.863826.455013.4016,257.63103.900.630.940.040.153.67
Minimum32.76297.0998.001169.412.300.040.000.000.060.23
2012Average383.951700.311760.436223.6732.670.370.350.020.111.26
Variance75,512.841,698,756.092,312,417.0633,425,972.92697.180.050.070.000.001.05
Maximum898.544348.505163.4517,120.1699.730.630.880.040.163.61
Minimum34.32380.63105.031279.642.340.040.000.000.070.23
2013Average400.431892.621799.006806.9633.020.420.360.020.111.23
Variance88,584.021,900,441.682,425,523.2438,804,566.21700.950.040.090.000.001.04
Maximum967.144454.555333.5717,971.06100.180.620.980.040.153.56
Minimum36.10448.2098.111437.042.560.040.000.000.070.22
2014Average416.912062.031791.557317.6930.410.420.350.020.121.17
Variance101,984.852,096,958.272,452,465.9644,691,526.73359.550.040.080.000.000.92
Maximum1034.584889.505496.4618,993.8766.020.630.830.040.163.49
Minimum38.81678.07107.661580.502.710.050.000.000.070.23
2015Average431.512299.131887.087758.3929.700.410.310.020.151.11
Variance114,433.112,394,742.572,722,659.3152,982,707.12325.650.030.060.000.000.75
Maximum1100.805405.955688.8920,155.9867.330.620.740.040.223.19
Minimum39.65726.85117.571691.852.900.070.000.000.090.27
2016Average446.182496.351891.418378.4630.220.400.310.020.140.97
Variance128,117.702,761,952.372,839,642.3160,900,156.05312.080.030.060.000.000.68
Maximum1165.735703.595852.6020,930.5165.550.610.680.040.233.05
Minimum38.97640.43121.831810.673.110.070.000.000.080.23
Note: No. of employees, investment in fixed assets and energy consumption are input variables; GDP and CO2 emissions are output variables; the ratio of secondary industry GDP to total GDP, the ratio of coal consumption to total energy consumption, the ratio of R&D investment to GDP, the ratio of local fiscal expenditure to GDP and the ratio of import and export to GDP are environmental factors.
Table 4. Between input and output variables, the Spearman’s rank correlation coefficients.
Table 4. Between input and output variables, the Spearman’s rank correlation coefficients.
InputsLaborCapitalEnergy
Outputs
GDP0.782 ***
(0.004)
0.827 ***
(0.0001)
0.627 **
(0.04)
CO2 emission0.873 ***
(0.0004)
0.927 ***
(0.00004)
0.709 ***
(0.01)
Note: The significance level of *** and ** are 1% and 5% respectively. The relevant p-value is indicated by numbers in brackets.
Table 5. The selection of external environmental variables.
Table 5. The selection of external environmental variables.
FactorsExplanatory Variables
Industrial structure ISRatio of secondary industry GDP to total GDP
Energy structure ESRatio of coal consumption to the total energy consumption
Technology TRatio of R&D investment to total GDP
Government intervention GIRatio of local fiscal expenditure to total GDP
Openness of economy OERatio of import and export to total GDP
Table 6. Comprehensive efficiencies of the GBA cities from 2010 to 2016.
Table 6. Comprehensive efficiencies of the GBA cities from 2010 to 2016.
Cities2010201120122013201420152016Average
Guangzhou0.3930.4380.4780.4390.4250.4050.3960.425
Shenzhen11111111
Zhuhai0.3070.3020.3090.3000.3480.3330.2940.313
Foshan0.3740.3710.3780.3480.3590.3730.3540.365
Jiangmen0.1940.2070.2280.2280.2880.2690.2230.234
Zhaoqing0.2650.2640.2690.2590.3150.3000.2760.278
Huizhao0.1500.1590.1680.1680.2040.2060.1940.178
Dongguan0.3180.3560.3590.3220.3760.4520.4220.372
Zhongshan0.2030.2150.2310.2430.3440.3280.3150.269
Hong Kong SAR11111111
Macao SAR11111111
Average0.4730.4830.4930.4820.5140.5150.4980.494
Table 7. SFA model parameters and estimation results in the second stage.
Table 7. SFA model parameters and estimation results in the second stage.
Slacks
Explanatory VariableNumber of EmployeesInvestment in Fixed AssetsEnergy Consumption
Constant term6.66
(42.56) ***
6.78
(3.72) ***
9.09
(24.11) ***
The ratio of secondary industry GDP to total GDP−0.07
(−0.90)
1.97
(2.99) ***
0.23
(1.48)
The ratio of coal consumption to total energy consumption−0.02
(−0.46)
−0.25
(−0.73)
0.05
(0.66)
The ratio of R&D investment to GDP0.09
(2.24) **
0.25
(0.64)
0.15
(1.70) *
The ratio of local fiscal expenditure to GDP−0.04
(−0.69)
−1.12
(−1.66) *
−0.07
(−0.54)
Import and export to GDP−0.20
(−3.26) ***
−1.52
(−3.39) ***
−0.30
(−2.22) **
sigma-squared12.08
(2.34) **
6.25
(2.12) **
18.32
(3.07) ***
gamma0.99
(11414.37) ***
0.93
(26.08) ***
0.99
(4651.11) ***
LR test of the one-sided error411.04 *102.42 ***359.24 ***
Note: The significance level of ***, ** and * are 1%, 5% and 10% respectively. The matching t-statistics of the computed parameters are indicated by numbers in brackets 5.1.3. The third phase: the real energy efficiency calculation results.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Fong, W.; Sun, Y.; Chen, Y. Examining the Relationship between Energy Consumption and Unfavorable CO2 Emissions on Sustainable Development by Going through Various Violated Factors and Stochastic Disturbance–Based on a Three-Stage SBM-DEA Model. Energies 2022, 15, 569. https://doi.org/10.3390/en15020569

AMA Style

Fong W, Sun Y, Chen Y. Examining the Relationship between Energy Consumption and Unfavorable CO2 Emissions on Sustainable Development by Going through Various Violated Factors and Stochastic Disturbance–Based on a Three-Stage SBM-DEA Model. Energies. 2022; 15(2):569. https://doi.org/10.3390/en15020569

Chicago/Turabian Style

Fong, Wengchin, Yao Sun, and Yujie Chen. 2022. "Examining the Relationship between Energy Consumption and Unfavorable CO2 Emissions on Sustainable Development by Going through Various Violated Factors and Stochastic Disturbance–Based on a Three-Stage SBM-DEA Model" Energies 15, no. 2: 569. https://doi.org/10.3390/en15020569

APA Style

Fong, W., Sun, Y., & Chen, Y. (2022). Examining the Relationship between Energy Consumption and Unfavorable CO2 Emissions on Sustainable Development by Going through Various Violated Factors and Stochastic Disturbance–Based on a Three-Stage SBM-DEA Model. Energies, 15(2), 569. https://doi.org/10.3390/en15020569

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