Next Article in Journal
Evaluation of the Protection Effectiveness of Natural Protected Areas on the Qinghai–Tibet Plateau Based on Ecosystem Services
Next Article in Special Issue
The Role of Big Data in Promoting Green Development: Based on the Quasi-Natural Experiment of the Big Data Experimental Zone
Previous Article in Journal
Implementation and Strategies of Community Music Activities for Well-Being: A Scoping Review of the Literature
Previous Article in Special Issue
The Relationship between Environmental Regulation, Green-Technology Innovation and Green Total-Factor Productivity—Evidence from 279 Cities in China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Identification of Key Carbon Emission Industries and Emission Reduction Control Based on Complex Network of Embodied Carbon Emission Transfers: The Case of Hei-Ji-Liao, China

1
School of Business, Shenyang University, Shenyang 110064, China
2
School of Information Management, Nanjing University, Nanjing 210023, China
3
Institute of Industrial and Economic Policy, Beijing Economic and Technological Development Zone (BDA), Beijing 100070, China
4
School of Public Finance and Taxation, Capital University of Economics and Business, Beijing 100070, China
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2023, 20(3), 2603; https://doi.org/10.3390/ijerph20032603
Submission received: 16 December 2022 / Revised: 20 January 2023 / Accepted: 25 January 2023 / Published: 31 January 2023
(This article belongs to the Special Issue The Recent Development of Environmental Management in Asia)

Abstract

:
Similar to the problems surrounding carbon transfers that exist in international trade, there are severe carbon emission headaches in regional industrial systems within countries. It is essential for emission reduction control and regional industrial restructuring to clarify the relationship of carbon emissions flows between industrial sectors and identify key carbon-emitting industrial sectors. Supported by the input–output model (I-O model) and social network analysis (SNA), this research adopts input–output tables (2017), energy balance sheets (2021) and the energy statistics yearbooks (2021) of the three Chinese provinces of Hei-Ji-Liao to construct an Embodied carbon emission transfer network (ECETN) and determine key carbon-emitting industrial sectors with a series of complex network measurement indicators and analysis methods. The key abatement control pathways are obtained based on the flow relationships between the chains in the industrial system. The results demonstrate that the ECETNs in all three provinces of Hei-Ji-Liao are small-world in nature with scale-free characteristics (varying according to the power function). The key carbon emission industry sectors in the three provinces are identified through centrality, influence, aggregation and diffusion, comprising coal mining, the chemical industry, metal products industry, machinery manufacturing and transportation in Liaoning Province; coal mining, non-metal mining, non-metal products, metal processing and the electricity industry in Jilin Province; and agriculture, metal processing and machinery manufacturing in Heilongjiang. Additionally, key emission reduction control pathways in the three provinces are also identified based on embodied carbon emission flow relationships between industry sectors. Following the above findings, corresponding policy recommendations are proposed to tackle the responsibility of carbon reduction among industrial sectors in the province. Moreover, these findings provide some theoretical support and policy considerations for policymakers.

1. Introduction

Global warming, which has become a common challenge for the whole world, is not only causing ecological damage but also threatening the survival and development of mankind. From the United Nations Framework Convention on Climate Change (UNFCCC) in 1992 to the Paris Agreement in 2015, the international community has reached a high level of consensus on reducing CO2 and other greenhouse gas emissions and implementing green and low-carbon development [1]. China’s economic growth has attracted worldwide attention since its reform and opening up, while it has consumed a large amount of energy. In the face of this situation, China proposed the “2030 China Emission Reduction Target” in 2015 to reduce CO2 emissions per unit of GDP by 60% to 65% in 2030 compared to 2005 and reach the peak of CO2 emissions as early as possible. In 2020, China went a step further and claimed that carbon emissions will peak by 2030 and the economy will become carbon neutral by 2060. China has been working towards this goal [2].
However, it is not easy for China to reduce emissions [3]. From one perspective, over 90% of China’s total carbon emissions come from industrial production and the Chinese government has introduced a series of policies to control the reduction of direct carbon emissions from industries where energy is the main driver, so as to achieve its emissions reduction targets as soon as possible [4]. However, it is challenging to optimize and upgrade the industrial structure to green energy and low-carbon in a short period of time and all industrial sectors in the industrial system are interdependent and integrated [5]. There are complex production and consumption activities between industrial sectors, leading to large amounts of direct carbon emissions that are not only used to supply their own demand for final products [6], but also flow into other industrial sectors in the form of Embodied carbon emissions along the industrial chain and the flow of intermediate products. In this way, a complex network of carbon emission flows is formed. Therefore, a single abatement policy is not enough, as illustrated in Figure 1. From another perspective, China is a vast country and the diversity of its geographical location and resource endowments has resulted in a wide range of different levels of development and importance of industrial sectors in different regions [7]. It is only through a diversified approach that low carbon emission reduction policies can be better suited to the industrial development of different regions [8].
Therefore, the relationship between Embodied carbon emission flows in different regions and different industrial sectors should be researched to realize China’s economic transformation toward green and low-carbon development [9]. Identifying the key carbon-emitting industrial sectors in the Embodied carbon flow network and analyzing the key Embodied carbon emission pathways will assist policymakers in formulating industrial restructuring plans [10], clarifying the implementation steps of emission reduction targets, improving energy utilization [3,11], promoting cleaner production in enterprises and mitigating overall environmental pollution emissions [12].
In international research, the concept of “Embodied carbon emissions” came to prominence in the 1970s. The first to introduce the concept was the International Federation of Institutes of Advanced Studies (IFIAS) Working Group on Energy Analysis. This concept was employed to measure the total amount of a resource consumed directly and indirectly throughout the industrial process of a product or service [13,14]. In principle, any resource name can be attached after “Embodied” (e.g., Embodied carbon and Embodied energy) [15,16,17]. As Embodied accounting is relatively similar to input–output economics, many of the concepts in input–output models can be applied to the Embodied analysis. Hence, the concept of “Embodied carbon emissions” is the sum of the direct and indirect carbon emissions emitted by goods and services throughout the production process [18]. The nature of economic activity suggests that Embodied carbon emissions are a crucial element of CO2 emissions in the economic system [19]. Therefore, it is necessary to analyze the structural characteristics of the Embodied carbon flow network and identify the key carbon-emitting industrial sectors in the industrial system and the key Embodied carbon flow pathways [20]. This will help design scientifically sound inter-industry synergistic emissions reduction policies and lay the foundation for better and cleaner production in industrial sectors and improved responsibility towards carbon reduction.
The current research on Embodied carbon emissions has been conducted from two main perspectives: the trade level and the industry level. In the case of international trade, trade liberalization has a vital impact on greenhouse gas emissions. International trade has experienced unprecedented growth since the Second World War, with global trade increasing 33-fold between 1950 and 2012 and the share of global trade in GDP increasing from 5.5% to 21%. Simultaneously, greenhouse gas (GHG) emissions have risen significantly. This issue is increasingly being explored in depth by scholars. Antweiller et al. [21] used a general equilibrium model of trade and environment to test the effect of international trade opening on SO2 concentrations. Meanwhile, Grether et al. [22] performed a qualitative analysis to demonstrate the impact of international trade on global SO2 emissions. These studies reveal that trade openness has led to some environmental improvements in some areas. Nonetheless, they ignore the impact of trade on greenhouse gas emissions. Based on this study, Cole and Elliott [23] analyzed the impact of trade on four environmental indicators, such as CO2. Using data on CO2 from 32 developed and developing countries in 1975–1995, they concluded that trade openness leads to an increase in CO2 emissions and the scale effect is much larger than the technology effect. Frankel and Rose [24] examined open relationships for seven environmental quality indicators including CO2 for a given level of per capita income in several countries, with CO2 data from 150 countries. The results implied that trade growth has a negative impact on CO2 emissions. Since CO2 emissions are global and some of the costs will be borne by foreigners, there will be no strict CO2 emission limits in their own countries. In other words, most studies support the conclusion that international trade increases CO2 emissions.
Regarding domestic inter-regional trade, numerous studies have revealed that, as international economic development becomes uneven, less developed regions and high-carbon polluting industries from developed regions are transferred. In the case of China, there is a spatial shift in carbon emissions from the east to the central and western regions. Wang et al. [25] measured carbon emissions and inter-provincial carbon transfers in multiple provinces in China with a multi-regional input–output model. The results unveil that the provinces in China with the largest net carbon transfer out are located in the eastern coastal, southern coastal and Beijing-Tianjin-Hebei regions and the provinces with the largest net carbon transfer in are located in the northwest region. Wang et al. [14] argued that differences in regional net carbon transfers are related to regional economic development strategies, regional division of labor and regional industrial structure. Chen et al. [26] investigated the issue of inter-regional carbon equity and concluded that there is carbon inequity in regional trade, with Embodied net carbon transfer in regions mostly located in less developed regions and Embodied carbon transfer out regions mostly located in economically developed regions. The above studies demonstrate the patterns of Embodied carbon emissions transfer between different regions and identify the main contributors to Embodied carbon emissions, contributing to the formulation of regional carbon reduction policies.
In addition to the trade aspect, there are more complex relationships between industries in terms of carbon emissions flows. Wang et al. [27] explored the Embodied carbon emission flows in regional industrial chains by combining input–output methods and DEA models to identify the key industrial sectors and key pathways of the Embodied carbon emission flows between regions. This study is a critical breakthrough in regional inter-industry Embodied carbon emission flows, while overlooking the relationship of Embodied carbon emission flows between industrial sectors from an industrial system perspective. A more effective research method that can reveal the dynamic interactions between sectors holistically should be selected to analyze the Embodied carbon flows between industrial sectors from an industrial system perspective. The complex network model is a preferable solution to these problems. It is an effective research method for analyzing the internal structural characteristics of complex systems and can be utilized to uncover the relationships between the internal parts of the system and determine the functions, roles and values of the network nodes. Complex network models are currently adopted in a wide range of fields. Ma et al. [28] investigated some deep features of the network of ore-energy flows in global trade based on complex network theory. Jiang et al. [29] explored the network topology of Embodied energy flows between global industry sectors and derived important evolutionary patterns. Jiang et al. [30] combined a complex network model with the EEIO model to characterize the macroscopic nature of the implicit carbon transfer network in China’s industrial sectors.
As suggested in the above research on Embodied carbon flows, much of the current research focuses on the relationship between Embodied carbon flows in the trade process, involving international trade and domestic trade. Research on Embodied carbon emissions flows from the industrial sector stresses the relationship between industrial sectors at the national level. There is little research on the relationship between industrial sectors and carbon emissions among provinces and regions. On the choice of research method, Moran et al. [31] constructed a multi-regional input–output model covering 41 countries and 35 industries using MRIO and WIOD to calculate the Embodied carbon flows in China’s foreign trade from 1995 to 2011. Lin et al. [10] examined the factors influencing the dynamics of Embodied carbon emissions from the industrial sector in China by combining the I–O model with the LMDI analytical model. These studies demonstrate direct links between carbon emissions from different industrial sectors, but do not effectively capture the indirect links between the various industrial sectors in the industrial system. Moreover, they fail to recognize the position, role and value of the specific industrial sector within the industrial system. The paper aims to answer the following key questions.
(1) How can key carbon-emitting industrial sectors be identified at the provincial level?
(2) How can key Embodied carbon emission flow paths be captured?
(3) How can effective inter-provincial carbon emission reduction, industrial structure optimization and clean production capacity be provided?
With the purpose of filling these gaps, the input–output model and social network analysis method are combined in this research to explore the important position of the industrial sector in China’s industrial development in the three provinces of Heilongjiang, Jilin and Liaoning. Meanwhile, the Embodied carbon emission flow relationships among industrial sectors in the three provinces of Hei-Ji-Liao are analyzed deeply at the provincial level, key carbon-emitting industrial sectors are identified and key Embodied carbon emission flow paths are derived. This paper provides theoretical support and policy formulation reference for reducing carbon emissions, contributing to optimizing industrial structure and improving clean production capacity between provinces.
The remainder of the paper is organized as follows. In Section 2, the model and the data sources are presented. In Section 3, the empirical results are analyzed. In Section 4, the results are discussed and policy implications are provided. Finally, conclusions are drawn in Section 5.

2. Research Methods and Data

2.1. Research Framework and Model Selection

In this research, an input–output model and a complex network model were chosen to analyze the relationship between embodied carbon emissions flows among industrial sectors in Heilongjiang, Jilin and Liaoning provinces of China, so as to identify the key carbon-emitting industries in each province. Figure 2 illustrates the research framework and model selection for this paper, with specific model construction details described in Section 2.1 and Section 2.2, respectively. Section 2.4 provides the data sources and calculation process.

2.2. Data Sources

In this paper, four main types of data sets are used, consisting of input–output (IO) tables, energy balance sheets, energy statistics yearbooks and Guidelines for National Greenhouse Gas Inventories (IPCC). IO tables are published by the statistical bureaus of Liaoning, Jilin and Heilongjiang provinces. The energy balance sheet and the Energy Statistics Yearbook (for Liaoning, Jilin and Heilongjiang provinces) cover eight major fossil energy sources: coal, coking coal, crude oil, gasoline, paraffin, fuel oil and natural gas [32,33]. Based on the relevant parameters in the IPCC Guidelines for National Greenhouse Gas Inventories, the CO2 emission factor λ k for energy source k can be calculated as:
λ k = N C V k × C C k × 44 12
where N C V k represents the average low-level heat generation of the energy type k ; C C k indicates the carbon emission factor of the energy type k ; 44 12 denotes the ratio of the relative molecular mass of CO2 to C. The CO2 emission factors for each energy source in each province calculated according to Equation (12) are presented in Appendix A. Energy balance sheets and Energy Statistics Yearbook employ the 2021 data set and IO tables for each province use the 2017 data set. Additionally, data processing is performed to merge the 42 industrial sectors in the IO table into 30, so as to ensure that the industrial sector classification in the IO table is consistent with the industrial energy consumption classification in the Energy Statistics Yearbook. The merged industrial sector classification is detailed in Table 1.

2.3. Analysis of Input–Output Model

2.3.1. Input–Output Analysis

The input–output analysis is a method based on the Leontief inverse matrix. It is intermediate use plus end use equals total output [34,35], expressed in mathematical notation as:
j = 1 n x i j + Y i = X i    ( i ,   j = 1 ,   2 ,   ,   n )
where x i j represents the value of the production of industry sector i to demand for the product of industry sector j and the value of intermediate inputs in industry sector i is obtained by summing up with j ; Y i denotes the amount of end-use value produced by industry sector i to satisfy final needs; X i signifies the total output of industry sector i .

2.3.2. Direct Consumption Coefficient

The industry sector is represented by A as a matrix of direct consumption factors. In matrix A , the element in row i and column j is a i j , which indicates the amount of value directly consumed by the production of industry sector i on industry sector j . a i j is calculated as:
a i j = x i j X j
Shifting the term in Equation (3) gives:
x i j = a i j X j
Substituting Equation (4) into Equation (2) yields:
j = 1 n a i j X j + Y i = X i
In Equation (5), the total output of the n industrial sectors can be obtained as a system of n linear equations. The direct consumption factor matrix A , the end-use matrix Y , the unit matrix I and the total output matrix X are expressed as:
A = [ a 11 a 12 a 1 ,   n 1 a 1 n a 21 a 22 a 2 , n 1 a 2 n a n 1 , 1 a n 1 , 2 a n 1 , n 1 a n 1 , n a n 1 a n 2 a n , n 1 a n n ]   Y = [ Y 1 Y 2 Y n 1 Y n ]
I = [ 1 0 0 0 0 1 0 0 0 0 1 0 0 0 0 1 ]   X = [ X 1 X 2 X n 1 X n ]
Then Equation (5) can be written as:
[ a 11 a 12 a 1 ,   n 1 a 1 n a 21 a 22 a 2 , n 1 a 2 n a n 1 , 1 a n 1 , 2 a n 1 , n 1 a n 1 , n a n 1 a n 2 a n , n 1 a n n ] × [ X 1 X 2 X n 1 X n ] + [ Y 1 Y 2 Y n 1 Y n ] = [ X 1 X 2 X n 1 X n ]
Then,
A X + Y = X
Transforming Equation (6) by shifting the terms gives:
X = ( I A ) 1 Y
where ( I A ) 1 indicates the Leontief inverse matrix, which creates a bridge between end-use and total input. Specifically, either end-use and total input is a constant and the other can be derived.

2.3.3. Competitive Embodied Carbon Emissions Model

CO2 emissions from energy consumption in the industrial sector [13,36] are expressed as:
C i = k = 1 n C i k = k = 1 n ( θ i k × k )                               ( i , k = 1 , 2 , , n )
where C i denotes the emissions of CO2 from the sum of energy consumption in industry sector i ; C i k represents the CO2 emissions resulting from the consumption of energy type k by industry sector i ; θ i k signifies the physical amount of energy of type k consumed by industry sector i ; k stands for the CO2 emission coefficient for energy of type k .
The CO2 emissions from the direct input energy requirement per unit of output in industry sector i are defined as the direct emission factor [37,38], denoted as E i ( i , j = 1 , 2 , , n ) , which is calculated by:
E i = C i X i = k = 1 n C i k = k = 1 n ( θ i k × k ) X i
The formula for calculating the implied carbon emissions C due to end use is Equation (10):
C = E X = E ( I A ) 1 Y             ( E   i s   a   d e t e r m i n a n t   c o n s i s t i n g   o f   E i . )

2.4. Analysis of ECETN Model Structure

In this section, the ECETN model is constructed and metrics from complex networks are utilized for overall network analysis.

2.4.1. Construction of the ECETN Model

Nodes, edges and weights are key factors in the construction of complex networks [39]. Through the calculations in Section 2.1 and Section 2.2, a complex network for industry sectors is constructed as nodes, the relationship of embodied carbon emission flows between two sectors is constructed as directed edges and the weights of the edges are constructed as embodied carbon emissions, which is ECETN. Furthermore, the structural features of the ECETN model are captured by analyzing the complex network metrics of the ECETN model. In the ECETN model, the number of nodes and edges is based on economic relationships. There is a flow of products or services between nodes and this generates a flow of embodied carbon emissions, as exhibited in Figure 3.
A threshold that can effectively evaluate the relationship of Embodied carbon emission flows is set in this paper to further explore the relationship of Embodied carbon emission flows between industry sectors. When the flow is less than or equal to the threshold, the Embodied carbon emissions flow between industry sectors is negligibly small and then the edge between nodes does not exist [40]. When the flow is greater than the threshold, the Embodied carbon emissions flow between industry sectors is large and is written as an edge:
f i j { 1   ( e i j > ε i ) 0     ( e i j ε i )               ε i = 1 n j = 1 n e i j
where f i j represents the Embodied carbon emission flow relationship between industry sectors; e i j indicates the implied carbon flows between industry sectors i and j ; ε i denotes the threshold value, which is the average Embodied carbon emissions between the two nodes.

2.4.2. Analytical Indicators for the ECETN Model

In this section, the ECETN model is analyzed using complex network metrics, as detailed below.
(1) Centrality Analysis
Centrality, as one of the main focuses of research in complex network analysis, describes the central position that nodes have in the network. In this paper, degree centrality, intermediate centrality and betweenness centrality are used for analysis [6].
① Degree Centrality
Degree centrality refers to the number of other industry nodes directly connected to the industry node. The more other nodes are connected, the more important and influential the node is in the ECETN. The formula is:
C R D ( x ) = O d x + I d x ( 2 n 2 )
where n denotes the network size, O d x indicates the point-out degree of x , I d x refers to the point-in degree of x and C R D signifies the degree centrality of x .
② Betweenness Centrality
According to the American sociologist Freeman (1979), if a node is between several nodes in a complex network, it must have a low degree and this relatively low degree point may play a vital intermediary role. Thus, it is at the center of the network. In ECETN, betweenness centrality is adopted to describe the degree of control that industry nodes have over the flow of Embodied carbon emissions in a complex network. The formula is:
  b j k ( i ) = g j k ( i ) g j k
C A B i = j n k n b j k ( i ) , j k i ,       j < k
where g j k indicates the number of edges between point j and point k ; g j k ( i ) denotes the number of edges passing through the third point i .
③ Closeness centrality
Closeness centrality describes the distance between an industry node and other nodes in the ECETN. The shorter the distance, the easier it is for the node to Embodied carbon emission flows in the ECETN. Hence, the more likely it is to be at the center of the ECETN, the higher the centrality. It is calculated as:
C A P i 1 = j = 1 n d i j
C R P i 1 = C A P i 1 n 1
where d i j denotes the distance between point i and point j connecting the edges, as well as the number of edges.
(2) Influence analysis
Influence describes the ability of the Embodied carbon emission flows of the nodes in the ECETN to act on other industries. In this paper, the direct consumption factor matrix is employed to assign weights to the adjacency matrix. The KATZ index is utilized to measure the influence of industries, calculated as:
T = a C + a 2 C 2 + + a k C k + = k = 1 a k C k = ( I a C ) 1 I
where a represents the decay factor, C k denotes the number of indirect paths of i , j through k   (after weighting by the direct consumption factor matrix), t i j refers to the elements of the matrix T , t j = i t i j signifies the column sums and t i = j t i j indicates the row sums.
(3) Clustering Analysis
In ECETN, K-cores, Main Core Networks and clustering coefficients are chosen to complete the screening of key carbon emission industry nodes for the aggregation analysis of Embodied carbon emission flows [29,41,42].
① K-core diagrams are mainly used to reflect the importance and agglomeration of nodes. They are one of the methods to analyze cohesive subgroups of complex networks based on the degree of nodes. K-cores are defined as follows. Let the graph G = ( V ,   E ) , V be the set of points in graph G and E be the set of edges in graph G (if a is G directed graph, E denotes the set of arcs), defined as W V . Then, the subgraph H k = ( W , E | W ) in graph G is a K-cores and also needs to satisfy when and only when v W , the degree d e g r e e ( v ) k of the point v and H k is the largest subgraph with this feature.
② Main Core Network
Generally, the K-cores subgraph with the largest kernel value is called the Main Core Network ( M n ). In the Main Core Network ( M n ), the industry nodes exhibit strong connectivity and high aggregation. In this paper, the K-cores of the ECETN are first determined. Afterward, the Main Core Network is identified by the maximum kernel value. Finally, the industry nodes in the Main Core Network are the most qualified as key carbon-emitting industries [43].
③ Clustering coefficients
Apart from these two indicators, the Clustering Coefficient is employed to quantitatively describe the degree of aggregation of ECETN. Specifically, a larger Clustering Coefficient indicates a higher degree of aggregation of ECETN [44]. The higher the Clustering Coefficient, the higher the degree of aggregation of the ECETN. The specific calculation formula is:
C i = E i k i ( k i 1 ) 2
C = 1 N i = 1 N C i
where E i denotes the actual number of edges that exist between the k i neighboring nodes of node i , C i represents the clustering coefficient of an industry node i of degree k i in the network and C refers to the clustering coefficient of the ECETN.
(4) Diffusivity analysis
In ECETN, key carbon-emitting industries also need to measure the strength of their diffusion ability in the network, that is, the diffusivity of the network nodes, in addition to having a high degree of aggregation. In this research, the strength of association in key carbon-emitting industries is evaluated using the two main indicators of the extent of association and diffusion effects.
  • ① Extent of association
The extent of association is made up of the forward and backward association extents, describing the number of nodes directly connected to other nodes, namely, the number of industries in the ECETN directly supplied/demanded by the Embodied carbon emission flows. In this research, the node’s out-degree ( O d i ) and in-degree ( I d i ) are calculated to quantify the number of other industries in the network directly driven by the industry’s Embodied carbon emissions, so as to measure the importance of the node in the ECETN.
  • ② Diffusion effects
The diffusion effect is a combination of nodes in the Main Core Network examining the primary key carbon-emitting industries, with the alternative nodes pulling or driving other nodes in the periphery through forward and backward linkages. The main two indicators are the forward and backward diffusion effects. The specific calculation formulas are:
D E i M o u t p u t i = d i O d i = 1 O d i j V M O D i j
D E j M i u t p u t j = d j I d j = 1 I d j i V M I D i j
where D E i o u t p u t indicates the forward diffusion effect of the alternative nodes in the Main Core Network; D E j i n p u t denotes the backward diffusion effect of the alternative nodes in the Main Core Network; V M represents the non-critical node at the periphery of the Main Core Network.

3. Empirical Results

3.1. Calculation of Embodied Carbon Emission Flows between Industrial Sectors

3.1.1. Direct Carbon Emissions from the Industrial Sector

In this research, the direct carbon emissions and direct carbon emission intensity due to energy consumption by industry sectors in Heilongjiang Province, Jilin Province and Liaoning Province for 2017 are estimated. As suggested in Figure 4a, direct carbon emissions from fossil energy consumption by sector in the three provinces are lowest in Jilin Province and highest in Heilongjiang Province. Electricity & heat production and supply (C24) is the largest direct carbon-emitting sector in Liaoning and Heilongjiang provinces. Petroleum processing and coking (C11) is the largest direct carbon-emitting sector in Jilin province, accounting for 41.2% (Liaoning), 43.6% (Heilongjiang) and 52.7% (Jilin) of total direct carbon emissions in each province. This is followed by metal smelting and rolling processing (C14), which account for 22.2% (Liaoning), 10.9% (Jilin) and 14.6% (Heilongjiang) of total direct carbon emissions by province [45].
In Figure 4b, the industrial sectors with high direct carbon emission intensity in the three provinces are petroleum processing and coking (C11), metal smelting and rolling processing (C14) and production and supply of heat and electricity (C24). These industrial sectors are all basic energy-based sectors. Since the main economic activity in the three provinces of Hei-Ji-Liao is industry, industrial development brings a large consumption of energy such as coal. However, the key industrial sectors that contribute to carbon emissions in each province should be captured, energy use should be effectively restructured and new energy industries should be vigorously developed, so as to lessen the increase in carbon emissions from energy consumption in the three provinces. As observed in Figure 4b, the direct carbon emission intensity of the industrial sectors in the three provinces is stable, except for those with high carbon emissions in C11, C14 and C24. This indicates that the sectors with high carbon intensity are initially controlled to assure the reduction of total carbon emissions in each province.
The direct carbon emissions caused by the consumption of energy in the production process of a particular industrial sector also have an impact on the carbon emissions of upstream and downstream industrial sectors along the chain and intermediate products, forming the flow of Embodied carbon emissions. Total carbon emissions from the industrial sector are composed of direct carbon emissions and Embodied carbon emissions. Considering the dynamics and complexity of the Embodied carbon emission flow process between industry sectors, the research perspective of complex networks is taken.

3.1.2. Embodied Carbon Emissions from the Industrial Sector

The Embodied carbon emissions and the Embodied carbon intensity of the industrial sector in Hei-Ji-Liao provinces were derived from the input–output model. As revealed in Figure 5, Embodied carbon emissions and Embodied carbon intensity fluctuate considerably in three provinces, suggesting that Embodied carbon emission flows between industrial sectors exert a significant impact on their carbon emissions. Particularly, Heilongjiang Province presents the most pronounced performance. In Figure 5a, the Embodied carbon emissions of all eight industrial sectors in Heilongjiang Province exceed their corresponding direct carbon emissions. Figure 5b also reflects the sharpest fluctuations in the fold of its Embodied carbon emissions intensity. From this perspective, the Embodied carbon emission flows between industry sectors should be emphasized to more effectively control the total amount of carbon emissions.

3.2. ECETN Model Construction and Topology Analysis

A network diagram and a ring diagram of the Embodied carbon flow network between industry sectors in Hei-Ji-Liao provinces are drawn using Gephi-0.9.2 and Circos-v19 software (Figure 6, Figure 7 and Figure 8) to fully demonstrate the topological structure of the Embodied carbon flow network between industry sectors. Theoretically, there should be 900 input–output relationships between the 30 industrial sectors. Nevertheless, the actual number of relationships existing in the 2017 input–output tables of the three provinces of Hei-Ji-Liao is less than 900, at 605 (Liaoning), 615 (Jilin) and 470 (Heilongjiang). Therefore, the Embodied carbon emission flows between the industrial sectors in each province are measured using the input–output relationships of the three provinces in a ring diagram (Figure 6b, Figure 7b and Figure 8b). Additionally, Embodied carbon flows between industrial sectors are analyzed, reflecting that the production and consumption of goods and services in industrial sectors are separated. The sector’s input and output to itself does not fall into this category. As a result, the sides of the Embodied carbon flow network are 575 (Liaoning), 585 (Jilin) and 440 (Heilongjiang). In Figure 6a, Figure 7a and Figure 8a, the circles represent the industry sectors and the numbers denote the industry sector numbers. The size of the circles represents the point intensity. The larger the circle, the higher the point intensity. The thickness of the edge demonstrates the weight of the implied carbon emission flow. The thicker the edge, the larger the implied carbon emission flow. The arrows signify the flow direction. The structural characteristic indicators of the corresponding ECETNs in Hei-Ji-Liao provinces are calculated using complex network theory. The results obtained are listed in Table 2.
The small-world character of complex networks is measured using clustering coefficient and average shortest path. Specifically, the value of the characteristic path length between points in the network is small, close to a random network, while the aggregation coefficient of the network is high, close to a regular network. The clustering coefficients of the three provincial ECETNs are 0.873 (Liaoning), 0.902 (Jilin) and 0.84 (Heilongjiang), as calculated by Equations (17) and (18). This indicates that the effective association number of the network is high and most industrial sectors are highly clustered. The average shortest paths for ECETN in the three provinces are 1.26 (Liaoning), 1.052 (Jilin) and 1.151 (Heilongjiang), implying that the network has better connectivity, with each industrial sector requiring only 1.26 steps/1.052 steps/1.151 steps to other sectors. Consequently, the three ECETNs have large clustering factors and small average shortest paths, all of which are small-world networks. The small-world character of the ECETN suggests that carbon emission flows between the various industrial sectors in the network are relatively tightly related. The changes in one sector of the network can quickly spill over to other sectors, leading to changes in the entire industrial system. From another perspective, the ECETN is extremely sensitive and vulnerable in certain circumstances, bringing crucial opportunities to collaborative emission reductions between industry sectors in the three provinces of Hei-Ji-Liao.

3.3. Analysis of Key Carbon-Emitting Industry Sectors in ECETN

As unveiled by measuring the correlation structure and correlation intensity of the nodes in the ECETN, the key carbon emission industry sector should have the following two characteristics. First, concerning correlation structure, the key carbon emission industry occupies a vital position in the ECETN and can influence the carbon emissions of the industrial system, especially for the carbon emissions of upstream and downstream industries. Second, regarding the intensity of the correlation, the key carbon-emitting industry has a strong ability to gather resources, guaranteeing that a small change in it can lead to changes in the surrounding industries or even the whole ECETN [46]. Additionally, key carbon-emitting sectors have a strong diffusion effect, revealing a diffusion effect to other industries. This assures that emission reduction measures in key sectors can spread to the ECETN and lead to emission reductions in all sectors. In summary, industries that meet the above criteria and have strong potential for development are considered key carbon-emitting industry sectors. In this section, the key carbon-emitting industry sectors in the ECETN are identified through analysis of the mediating role, impact intensity analysis and structural analysis of cohesive subgroups [25,47].

3.3.1. Analysis of the Mediating Role of the Industrial Sector

The degree centrality, intermediate centrality and near centrality are measured based on the concept of node centrality of ECETNs in Hei-Ji-Liao provinces to effectively evaluate the value of industrial nodes in ECETNs, with threshold values of each centrality as Liaoning province (17.073, 1.606, 27.891), Jilin province (14.634, 2.207, 25.786) and Heilongjiang Province (14.634, 1.188, 38.679). This research lists the series of industries with industry sector centrality greater than the threshold value in the ECETNs of the three provinces and considers the industries in the table as strong centrality industry sectors. The details are provided in Table 3, Table 4 and Table 5.
The strong centrality industries in Hei-Ji-Liao provinces in terms of degree centrality, intermediary centrality and proximity centrality are counted.
As suggested in Table 3, the industrial sectors with high Embodied carbon emission spillover capacity (i.e., the capacity of the CO2 emitted in the network) in the ECETN of Liaoning province are mostly concentrated in the metal and non-metal products (C15, C13), chemical manufacturing (C12, C25, C19), energy and mining (C2, C3, C4, C5, C24), manufacturing (C6, C17, C19, C20) and public services (C28, C29) sectors. Particularly, the three centrality indicators of C15 and C17 are both high, reflecting the important position of the two industries in ECETN, their strong control over other industries and their Embodied carbon emission spillover ability. They are essential industries worthy of attention.
In the ECETN of Jilin province (Table 4), sectors with high Embodied carbon emission spillover capacity are concentrated in metal and non-metal processing (C13, C14, C15), construction and product manufacturing (C16, C17, C22, C27), chemical manufacturing and others (C12, C23, C26), agriculture, forestry and fishery (C1), energy and mining (C2, C4, C5, C11, C24) and transportation (C18, C28) sectors. Specifically, the two centrality indicators, C16 and C17, are both high, demonstrating the importance of the two sectors in the ECETN and their strong control over other sectors, as well as the carbon emission spillover capacity. They are vital sectors worthy of attention.
In the ECETN of Heilongjiang province (Table 5), sectors with high Embodied carbon emission spillover capacity are concentrated in the agriculture, forestry and fishery (C1), construction (C27), wood products (C9), metal and non-metal products (C13, C14, C15), manufacturing (C16, C22) and energy processing (C12, C24) sectors. Particularly, the two centrality indicators, C1 and C15, are both high, implying the importance of these two sectors in the ECETN and their strong control over other sectors, as well as the ability to spill over carbon emissions. They are critical areas of concern.
The above analysis yields the key industrial sectors in the three provinces of Hei-Ji-Liao with strong Embodied carbon emission spillover capacity and mediating role. Then, the effective control of carbon emissions from these key industrial sectors can influence the flow of Embodied carbon emissions between their upstream and downstream sectors, reduce the flow of implied carbon emissions in the supply chain and thus curtail the carbon emissions of the whole ECETN. Concurrently, the clean development of these key industrial sectors can enhance the efficiency of energy use along the chain, contributing to the lessened carbon emissions of the industrial system.

3.3.2. Impact Intensity Analysis of Industry Sectors

In ECETN, the effect and influence ability of nodes on other nodes is called the point intensity of the industrial sector. In this paper, the 0-1 table is obtained by processing the input–output tables from three provinces of Hei-Ji-Liao. The KATZ index is calculated for the assignment adjacency matrix of the 0-1 table to obtain the point intensity value of each node. As suggested by calculating the impact intensity thresholds of the three provinces as 0.359 (Liaoning), 0.312 (Jilin) and 0.344 (Heilongjiang), the industries greater than this threshold in each province are considered high-impact industries and these industrial sectors carry a considerably large amount of Embodied carbon emissions, as listed in Table 6.
The results of KATZ index statistics in ECETNs (Table 6) demonstrate 23 industrial sectors with KATZ indexes greater than the threshold in Liaoning Province, except for oil and gas extraction, science and technology and services and public services. Additionally, the KATZ indexes of all industries are higher than the threshold (0.359), among which the top 15 industries are construction (C27), transportation and traffic (C18, C28), mining and mining products (C2, C4, C14, C15), energy (C11, C12, C24) and manufacturing (C17, C18, C19). The above industrial sectors as high-impact industrial sectors of implied carbon emissions directly influence the total amount of Embodied carbon emission flows in ECETN. Notably, the impact intensity of agriculture, forestry, animal husbandry & fishery and metal products is much higher than that of other sectors, indicating that these two sectors have a more significant pulling effect on the generation of carbon emissions in other sectors and act as aggregation centers of Embodied carbon emission flows in the network.
In Jilin province, there are 21 industries above the threshold, among which the public service industry (C29, C30) has a weak influence and the top 15 industries are: agriculture, forestry, animal husbandry & fishery (C1), processing and manufacturing (C16), energy and chemical industry (C2, C12, C24, C25, C23, C26), metallurgical manufacturing (C4, C5, C13, C14, C15) and Machinery manufacturing (C16). Among them, the mining industry is at the center of the Embodied carbon emission flow. In other words, coal and non-metal mining are both the supply and consumption centers of Embodied carbon emissions, are in the midstream of the industrial system and serve as crucial producers and suppliers of intermediate products.
In Heilongjiang Province, there are 18 industries above the threshold and the high-impact industries are more concentrated in agriculture (C1), construction (C27), metal products (C13, C15), energy and chemical industries (C3, C10, C12, C14) and manufacturing (C16, C17, C22). Agriculture, construction and manufacturing have always been vital pillar industries in Heilongjiang province, revealing that they exert a strong influence on the industrial system. They carry a large amount of Embodied carbon emissions while impacting the carbon generation and emissions of other industrial sectors.
After the screening results of centrality and influence in ECETNs of Hei-Ji-Liao provinces are compared, industries with strong centrality and high influence were selected as the primary key carbon emission industries. The results are presented in Table 7.
The distribution of the primary key carbon emission industries in the three provinces of Hei-Ji-Liao suggests that agriculture, coal mining and metal products are the most critical Embodied carbon emission industries in the industrial system, coinciding with the industrial structure of northeast China. Since these three industrial sectors are also critical supports for Northeast China’s industries, effective control of carbon emissions from these three industries can facilitate the achievement of carbon emission reduction targets in the regional economic systems of the three provinces.

3.3.3. Correlation and Cohesion Analysis of ECETNs

  • ① K-core and Main Core Network
The direct consumption coefficient matrix of Hei-Ji-Liao provinces was transformed into an adjacency matrix. Then, the K-cores of the adjacency matrix were visualized by Netdraw to obtain K-cores plots, as rendered in Figure 9.
M n is the submap of the ECETN with the largest K-cores, as well as the area with the highest density of correlations between the Embodied carbon emissions flows of the industrial sectors, which are strongly clustered. It is also the area where small changes can significantly lead to changes in the Embodied carbon emissions of neighboring industrial sectors or even industrial systems. Accordingly, the provincial main core networks were extracted and screened for high aggregation industrial sectors with Embodied carbon emission flows. As observed in Figure 9, the maximum core value is 9 for Liaoning Province and 7 for both Jilin and Heilongjiang Provinces, uncovering that the aggregation of Embodied carbon emission flows between industrial sectors in the main core network is slightly stronger in Liaoning than in Jilin and Heilongjiang.
There are 20 industrial sectors in Liaoning’s main core network, mostly concentrated in agriculture and animal husbandry (C1), mineral extraction (C2, C3, C5,), processing and machinery manufacturing (C7, C10, C13, C14, C15, C16, C17, C18, C19, C20), energy and chemicals (C11, C12, C25), distribution services (C28) and public services (C29, C30); 21 industrial sectors in the Jilin main core network, mostly concentrated in agriculture and animal husbandry (C1), mineral extraction (C2, C3, C4, C5), processing and manufacturing and energy (C6, C7, 13, C14, C15, C16, C17, C18, C19), construction and distribution services (C25, C28) and public services (C29, C30); 18 industrial sectors in Heilongjiang’s main core network, mostly concentrated in agriculture and animal husbandry (C1), mineral extraction (C2, C3, C4), processing and machinery manufacturing (C7, C13, C14, C15, C16, C17, C18, C20, C25), energy and chemicals (C12, C23, C25,), distribution services (C28) and public services (C29, C30).
  • ② Clustering coefficients
In addition to analyzing the industry sectors with strong aggregation of implied carbon emissions by determining the main core network using K-cores, clustering coefficients are adopted to determine the degree of aggregation of the industry sectors in the ECETN, with larger clustering coefficients indicating a higher degree of aggregation in the ECETN. The average clustering coefficients of ECETNs in the three provinces were measured to be 0.727 (Liaoning), 0.697 (Jilin) and 0.639 (Heilongjiang). The clustering coefficient thresholds for industry sectors in the three provinces were 0.614, 0.621 and 0.465 for Liaoning, Jilin and Heilongjiang, respectively. These values were for strongly clustered industry sectors. The industrial sectors whose clustering coefficient values were all greater than the threshold in the main core network were those with high aggregation of Embodied carbon emission flows.
The industry sectors with clustering coefficients greater than the threshold were selected as alternative key carbon-emitting industry sectors according to the comparison between the primary key carbon-emitting industry sectors identified above and the highly aggregated industry sectors in the main core network. The final alternative key carbon emission industry sectors for Liaoning are coal mining (C2), chemicals (C12), processing and manufacturing (C15, C17, C19) and distribution services (C28, C29). The alternative key carbon sectors in Jilin are mining and processing (C2, C5, C13, C14), processing and manufacturing (C15, C16), energy and chemicals (C12, C24) and distribution services (C28). The alternative key carbon sectors in Heilongjiang are agriculture and animal husbandry (C1) and metallurgical manufacturing (C13, C14, C15, C16). These are listed in Table 8.
  • ③ ECETNs’ cohesive structure analysis
Cohesive subgroup partitioning of ECETNs is performed based on modular algorithms using Gephi software. The maximum modularity of Liaoning, Jilin and Heilongjiang is 0.426, 0.301 and 0.214, respectively. This result is ideal for the division of cohesive subgroups, as illuminated in Figure 10. The size of the circles in the diagram represents the size of the point intensity of the industrial sector. The larger the circle, the greater the point intensity. The thickness of the edges represents the weight of the Embodied carbon flow. The thicker the edge, the greater the Embodied carbon flow. Figure 10 demonstrates that the nodes in the Liaoning ECETN are divided into two cohesive subgroups, three in Jilin and four in Heilongjiang.
The composition of industrial sectors and Embodied carbon flows within each cohesive subgroup of the three provinces are illustrated in Table 9. According to the analysis of Embodied carbon emission flows within the cohesive subgroups, Liaoning Province has the largest number of flows within subgroup I. Subgroup II has a larger number of industrial sectors, including manufacturing, construction and services and other related sectors; while these sectors are more interrelated and have more complex flow relationships, the overall mobility is not the largest. Although there are not many industrial sectors in Subgroup I, agriculture and coal mining are the pillar industries of the province and are closely related to the upstream and downstream industrial sectors, driving the overall carbon emission transfer volume. There are three cohesive subgroups in Jilin Province, in descending order: III, II and I. Subgroup III contains the largest number of alternative key carbon-emitting industries and their high-centeredness and strong influence result in large inter-sectoral carbon transfers within subgroup III. Meanwhile, subgroup I comprises no alternative key carbon-emitting industries and has the weakest carbon transfer capacity. The Embodied carbon emission flows of the cohesive subgroups in Heilongjiang Province are II, III, I and IV in descending order. The agriculture and livestock sector (C1) in the II cohesion subgroup has a reasonably strong capacity to move carbon emissions. It is a linked industry in several industrial sectors and, as the province’s mainstay industry, should receive focused attention.
Regarding the flow of Embodied carbon emissions between cohesive subgroups, the inflows and outflows of the subgroups in Liaoning province are balanced. Nevertheless, the inflows of subgroup II in Jilin province and subgroup II in Heilongjiang province are much larger than the outflows and are the absorption sites of Embodied carbon emissions in ECETNs. In contrast, the outflows of subgroup I in Jilin province and subgroup III in Heilongjiang province are much larger than the inflows. They have a stronger spillover effect of Embodied carbon emissions and are the dispersal sites of Embodied carbon emissions in ECETNs. The industrial sectors in some cohesive subgroups are in the middle and lower reaches of the chain and their development and growth require large amounts of energy and raw materials from the upstream industrial sectors, resulting in large carbon inflows. The upstream sectors of the chain are the suppliers of raw materials and essential support for the development of the downstream industries, contributing to significant carbon outflows.

3.3.4. Diffusivity Tests for Alternative Key Carbon-Emitting Industry Sectors

In this section, the diffusion of alternative key carbon-emitting industry sectors in ECETNs is examined by province, involving two main components: the breadth of association and the diffusion effect. Particularly, the breadth of association is calculated through the outward and inward degrees of the industry sector. The diffusion effect suggests that the forward diffusion effect, forward diffusion effect, backward diffusion effect and backward diffusion effect of the industry sectors are measured to determine the diffusion of Embodied carbon emissions in ECETNs for alternative key carbon-emitting industry sectors, as illustrated in Figure 11.
Overall, the diffusion effect in Liaoning (0.34) is better than in Heilongjiang (0.15) and Jilin (0.12). The analysis of alternative key carbon-emitting industry sectors demonstrates that, although some alternative key carbon-emitting industry sectors have strong centrality, high influence and high aggregation, Embodied carbon emission flows are only within the industry chain, compared to their weak contribution to the low carbon emission reduction of industry sectors outside the industry chain, and cannot promote cleaner production in the industry system; so they should not be considered as key carbon-emitting industry sectors and, for example, Liaoning (C17, C29), Jilin (C12, C16, C28) and Heilongjiang (C13) should be removed.
In Figure 11, the key carbon-emitting industrial sectors in Liaoning and Heilongjiang exhibit a better diffusion effect. The diffusion effect ratio of key carbon-emitting industry sectors in Liaoning is larger, implying that these sectors promote the rapid flow of Embodied carbon emissions in the industrial system. On the contrary, Jilin’s key carbon-emitting industrial sector is not satisfactory since it only has a unidirectional diffusion effect. Moreover, the diffusion effect is weak and does not significantly contribute to the Embodied carbon emission flow in the industrial system, revealing the difficulty of emission reduction control through the industrial sector in Jilin.
With respect to industry sectors, the Embodied carbon diffusion effects of C2, C15 and C19 in Liaoning are better than those of C12 and C28. In other words, the mining and manufacturing industries play a dominant role in carbon emissions in the industry system and promoting cleaner development of these industries improves the energy efficiency of the industry chain, leading to curtailed carbon emissions in the industry system as a whole. The values of C14 and C15 association breadth and diffusion effects in Jilin are both large, suggesting that metal manufacturing occupies a crucial position in the ECETN. Both its own Embodied carbon emissions and the Embodied carbon emission flows to the industrial chain and even to the industrial system have a huge impact. The key carbon-emitting industrial sectors in Heilongjiang are concentrated in agriculture and livestock (C1) and manufacturing (C14), both of which have strong diffusion capacity and the highest degree of association. Hence, effective control of carbon emissions from C1 and C14 can significantly lessen the overall emissions of the industrial system.

4. Discussion

In recent years, China’s rapid economic development has led to an increase in carbon emissions and thus increased pressure on China to reduce these. However, China, as a responsible power, has been advocating a new, fair and reasonable climate governance system at the international level. While developing its economy, China has fulfilled its commitment to energy conservation and emission reduction, especially after the incorporation of “carbon peaking and carbon neutral” into the overall layout of ecological civilization. Concurrently, China has placed more emphasis on environmental protection, green and low-carbon development and clean production to wrestle with the global climate change situation. In the social and economic system, the production process of industrial sectors is accompanied by a large number of carbon emissions and production and consumption between industrial sectors are closely related. In the production process, the industrial sectors constitute a complex network of carbon emissions between them. Energy saving and emission reduction can be achieved by effectively identifying the key nodes in the network and grasping the carbon emission relationship between industrial sectors. In this paper, ECETNs are established through the input–output model and the SNA method. Specifically, the structural characteristics of ECETNs in Hei-Ji-Liao provinces are demonstrated and the key carbon emission industry sectors in each province’s industrial system are identified through indicators such as centrality, influence, agglomeration and diffusion.
Firstly, the results of the overall network characteristics analysis suggest that the ECETNs in the three provinces have significant small-world characteristics. Thus, a close Embodied carbon emission transmission relationship is observed between most industrial sectors. Moreover, the changes in carbon emissions in any one industrial sector bring about a chain reaction of carbon emissions in upstream and downstream related industrial sectors in the industrial chain, which in turn will impact the changes in total carbon emissions in the whole industrial system.
Secondly, the individual structural characteristics of the network reflect that some industrial sectors are characterized by strong centrality and high influence in the ECETNs. In other words, a few industrial sectors play a crucial role in the network. Among them, agriculture, coal mining and metal products are the most critical carbon-emitting industries in the three provinces, coinciding with the industrial structure of the Northeast. Therefore, differentiated low-carbon measures should be adopted during the process of formulating industrial emissions reduction policies according to the positioning of individual industries in ECETNs.
Finally, an analysis of the clustering of the key industrial sectors of interest unveils clear clustering characteristics in the ECETNs of the three provinces. Coal mining (C2), chemicals (C12), processing & manufacturing (C15, C17, C19) and distribution services (C28, C29) in Liaoning; mining & processing (C2, C5, C13, C14), processing & manufacturing (C15, C16), energy & chemicals (C12, C24) and distribution services (C28) in Jilin; agriculture and livestock (C1) and metallurgical manufacturing (C13, C14, C15, C16) in Heilongjiang all have high agglomeration. The diffusivity tests of the high-concentration sectors demonstrate that mining and manufacturing in Liaoning dominate the carbon emissions of the industrial system and the energy efficiency of the industrial chain can be reinforced by promoting the clean development of these sectors, so as to curtail the carbon emissions of the whole industrial system. Jilin’s metal manufacturing industry occupies a key position in the ECETN. Both its own Embodied carbon emissions and the Embodied carbon emissions flowing into the chain or even the industrial system can exert a significant impact. The main carbon-emitting industrial sectors in Heilongjiang Province are concentrated in agriculture and manufacturing. The total emissions of the industrial system can be significantly reduced by effectively controlling their carbon emissions.
Based on the above study, the Embodied carbon emission flow focus pathways are captured following the maximum marginal power relationships of key carbon-emitting industrial sectors in the ECETNs of the three provinces to better accomplish regional emission reduction control. The specific results are provided in Table 10.
There are three critical paths in Liaoning’s ECETN: two starting with coal mining and metal products, one starting with the chemical industry and ultimately all ending with the construction and service sectors. There are four critical paths in Jilin’s ECETN: two paths starting with the energy sectors end up in transport equipment manufacturing; two paths starting with non-metal beneficiation and metal products end up in services. There are two critical paths in Heilongjiang’s ECETN: both paths, starting with metalworking and general equipment manufacturing, end up in the service sector. The key pathways of Embodied carbon emissions in Hei-Ji-Liao are similar, starting from traditional sectors such as mining & energy and ending in the construction and service sectors. Therefore, the key sectors of Embodied carbon emissions and the flow path of Embodied carbon emissions should be revealed when formulating carbon reduction policies. Besides, the role played by each industry sector along the key path should be demonstrated to discover the optimal carbon reduction measures. This will enable the effect of carbon reduction to move from ‘point’ to ‘line’ and from ‘line’ to ‘surface’ and to form regional linkages. On this basis, a carbon tax policy should also be considered. Many countries already have more mature carbon tax policies. China has been enriching and improving its own carbon tax policy since 2012. The ultimate goal is the sustainable development of a low-carbon economy.

5. Conclusions

This study is based on data from the input–output table of Hei-Ji-Liao in 2017, the energy balance table in 2021 and the energy statistics yearbook in 2021. The Embodied carbon emission flow relationships in the Hei-Ji-Liao industrial system were measured using input–output models and SNA. The results demonstrated that the Embodied carbon emissions of the three provinces were a key factor in increasing total carbon emissions, while the amount of Embodied carbon emissions was determined by several key carbon-emitting industry sectors. Through the above analysis, an Embodied carbon emission flow network (ECETN) between industry sectors of Hei-Ji-Liao was constructed. An analysis of the overall characteristics of ECETN revealed its small-world character. The individual characteristics of the ECETN were analyzed using complex network technical indicators (centrality, influence, agglomeration, diffusion) to identify the role and position of key carbon-emitting industrial sectors in the network. Following the Embodied carbon emission flow paths between key carbon-emitting sectors, relevant recommendations were provided for the formulation of carbon emission reduction and industrial restructuring policies.
Meanwhile, there are some limitations to this research. Firstly, the Embodied carbon emission flows between industry sectors in the province were only considered in constructing the ECETN, while the Embodied carbon emission flows between regions were neglected. In future research, the Embodied carbon emission flows between industrial sectors in the region will be further explored to ensure regional economic development while promoting regional low carbon emission reduction. Secondly, a large number of service industry sectors were combined in this study for the purpose of industry sector classification due to the limitations of energy consumption data. This failed to provide an understanding of the specific carbon emissions of the sub-sectors. In future research, the service sector will be refined and its sector-specific carbon emissions will be analyzed. Finally, this paper focuses on the analysis of Embodied carbon emissions and emission reduction control in Northeast China. The findings are not generalizable owing to the differences in geographical location and resource endowment. Therefore, multi-regional comparative studies will be conducted in the future to provide more precise policy recommendations for policymakers.

Author Contributions

All authors were involved in preparing the manuscript. Conceptualization, S.S., Y.L. and Z.Z.; methodology, S.S.; software, T.Z.; validation, Z.Z., W.Z. and Y.L.; formal analysis, S.S. and T.Z.; investigation, W.Z. and Z.Z.; resources, Y.L. and S.S.; data curation, Y.L.; writing—original draft preparation, T.Z. and S.S.; writing—review and editing, W.Z., Z.Z., Y.L. and S.S.; visualization, W.Z.; supervision, Z.Z. and S.S.; project administration, S.S. and Z.Z.; funding acquisition, S.S. All authors have read and agreed to the published version of the manuscript.

Funding

This paper was supported by grants from the Economic Development Research Foundation of Liaoning province (Grant No. 2023lslqnkt-056) and the Social Science Fund Project of Liaoning Province (Grant No. L20BZZ006).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

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

Conflicts of Interest

The author declares no conflict of interest.

Nomenclature

Parameters
i ,   j industry sector
Aa matrix of direct consumption factors
a i j the amount of value directly consumed by the production of industry sector i on industry sector j
I the unit matrix
X the total output matrix
Y the end-use matrix
C i the emissions of CO2 from the sum of energy consumption in industry sector i
k energy type
C i k the CO2 emissions resulting from the consumption of energy type k by industry sector i
θ i k the physical amount of energy of type k consumed by industry sector i
k When an enterprise adopts environmental the CO2 emission coefficient for energy of type k
E i the CO2 emissions from the direct input energy requirement per unit of output in industry sector i as the direct emission factor
f i j the Embodied carbon emissions flow relationship between industry sectors
e i j the Embodied carbon flows between industry sectors i and j
ε i the threshold value, which is the average Embodied carbon emissions between the two nodes
n the network size
O d x the point-out degree of x
I d x The risk loss brought by NIMBY facilities to surrounding residents
C R D the degree centrality of x
g j k the number of edges between point j and point k
d i j the distance between point i and point j connecting the edges
C k the number of indirect paths of i , j through k
D E i o u t p u t the forward diffusion effect of the alternative nodes in the Main Core Network
D E j i n p u t the backward diffusion effect of the alternative nodes in the Main Core Network
N C V k the average low-level heat generation of the energy type k
C C k the carbon emission factor of the energy type k
Variables
x i j the value of the production of industry sector i to demand for the product of industry sector j ,
Y i the amount of end-use value produced by industry sector i to satisfy final needs
Y i the amount of end-use value produced by industry sector i to satisfy final needs
X i the total output of industry sector i .
Acronyms
ECETNEmbodied carbon emission transfer network
I-O modelInput–output model
DEAData envelopment analysis
MRIOMulti-Regional input–output Model
SNAsocial network analysis
WIODWorld input–output Data

Appendix A

Table A1. CO2 emission factors for various energy sources.
Table A1. CO2 emission factors for various energy sources.
EnergyAverage Low-Level Heat Output
NCV (kJ/kg)
Carbon Emission Factor
(kg/106K\kJ)
CO2 Emission Factor
(kg-CO2/kg)
Crude Coal20,90826.322.02
Coke28,43531.383.27
Crude Oil41,81620.083.08
Petrol43,07018.902.98
Paraffin43,07019.603.10
Diesel oil42,65220.203.16
Fuel oil41,81621.103.24
Gas38,93115.322.19
Tabulated data from IPCC “Guidelines for National Greenhouse Gas Inventories (2006)”.

References

  1. Chen, Y.S.; Liu, S.; Wu, H. How can Belt and Road countries contribute to glocal low-carbon development? J. Clean. Prod. 2020, 256, 120717. [Google Scholar] [CrossRef]
  2. Liu, Q.J.; Qi, R.H.; Zhao, Y. Comparative analysis of the marginal abatement cost modeling for coal-fired power plants in China. J. Clean. Prod. 2022, 356, 131883. [Google Scholar] [CrossRef]
  3. Zhang, H. Effects of electricity consumption on carbon intensity across Chinese manufacturing sectors. Environ. Sci. Pollut. Res. 2019, 26, 27414–27434. [Google Scholar] [CrossRef] [PubMed]
  4. An, Y.F.; Zhou, D.Q.; Yu, J.; Shi, X.; Wang, Q. Carbon emission reduction characteristics for China’s manufacturing firms: Implications for formulating carbon policies. J. Environ. Manag. 2021, 284, 112055. [Google Scholar] [CrossRef]
  5. Lin, B.Q.; Huang, C.C. Analysis of emission reduction effects of carbon trading: Market mechanism or government intervention? Sustain. Prod. Consum. 2022, 22, 28–37. [Google Scholar] [CrossRef]
  6. Shi, J.; Li, H.; Guan, J.; Sun, X.; Guan, Q.; Liu, X. Evolutionary features of global embodied energy flow between sectors: A complex network approach. Energy 2017, 140, 395–405. [Google Scholar] [CrossRef]
  7. Huang, W.Q.; Wang, Q.F. Review of recent progress of emission trading policy in China. J. Clean. Prod. 2022, 349, 131480. [Google Scholar] [CrossRef]
  8. Xu, T.T.; Kang, C.Y.; Zhang, H. China’s efforts towards carbon neutrality: Does energy-saving and emission-reduction policy mitigate carbon emissions? J. Environ. Manag. 2022, 316, 115286. [Google Scholar] [CrossRef]
  9. Guo, F.; Meng, S.Q.; Sun, R.J. The evolution characteristics and influence factors of carbon productivity in China’s industrial sector: From the perspective of embodied carbon emissions. Environ. Sci. Pollut. Res. 2021, 28, 50611–50622. [Google Scholar] [CrossRef]
  10. Lin, B.Q.; Raza, M.Y. Coal and economic development in Pakistan: A necessity of energy source. Energy 2020, 207, 118244. [Google Scholar] [CrossRef]
  11. Tian, Y.X.; Yang, X.M. Asymmetric effects of industrial energy prices on carbon productivity. Environ. Sci. Pollut. Res. 2020, 27, 42133–42149. [Google Scholar] [CrossRef] [PubMed]
  12. Wang, C.; Zhan, J.Y.; Bai, Y.P.; Chu, X.; Zhang, F. Measuring carbon emission performance of industrial sectors in the Beijing-Tianjin-Hebei region, China: A stochastic frontier approach. Sci. Total Environ. 2019, 685, 786–794. [Google Scholar] [CrossRef] [PubMed]
  13. Gao, P.; Yue, S.J.; Chen, H.T. Carbon emission efficiency of China’s industry sectors: From the perspective of embodied carbon emissions. J. Clean. Prod. 2021, 283, 124655. [Google Scholar] [CrossRef]
  14. Wang, X.; Yu, J.; Song, J.; Di, X.; Wang, R. Structural evolution of China’s intersectoral embodied carbon emission flow network. Environ. Sci. Pollut. Res. 2021, 28, 21145–21158. [Google Scholar] [CrossRef] [PubMed]
  15. Dixit, M.K. Life cycle recurrent embodied energy calculation of buildings: A review. J. Clean. Prod. 2019, 209, 731–754. [Google Scholar] [CrossRef]
  16. Yang, W.X.; Gao, H.; Yang, Y.P. Embodied Carbon in China’s Export Trade: A Multi Region Input-Output Analysis. Int. J. Environ. Res. Public Health 2022, 19, 3894. [Google Scholar] [CrossRef]
  17. Yan, Y.F.; Wang, R.; Zheng, X.X.; Zhao, Z. Carbon endowment and trade-embodied carbon emissions in global value chains: Evidence from China. Appl. Energy 2020, 227, 115592. [Google Scholar] [CrossRef]
  18. Deng, G.Y.; Xu, Y. Accounting and structure decomposition analysis of embodied carbon trade: A global perspective. Energy 2017, 137, 140–151. [Google Scholar] [CrossRef]
  19. Fang, D.L.; Chen, B. Information-based ecological network analysis for carbon emissions. Appl. Energy 2019, 238, 45–53. [Google Scholar] [CrossRef]
  20. Li, X.Z.; Wang, X.J.; Zhang, Y. Spatial Differences in Emission Reduction Effect of Sterilization of Manufacturing Industry Export in China. Emerg. Mark. Financ. Trade 2021, 57, 2331–2355. [Google Scholar] [CrossRef]
  21. Antweiller, W.; Copeland, B.R.; Taylor, M.S. Is Free Trade Good for the Environment? Am. Econ. Rev. 2001, 91, 877–908. [Google Scholar] [CrossRef] [Green Version]
  22. Grether, J.M.; Mathys, N.A.; de Melo, J. Trade, Technique and Composition Effects: What is Behind the Fall in World-Wide SO2 Emissions 1990–2000? Fond. Eni Enrico Mattei Note Di Lav. 2007, 93. [Google Scholar] [CrossRef] [Green Version]
  23. Cole, M.A.; Elliott, R.J.R.; Okubo, T.; Zhang, L.Y. Importing, outsourcing and pollution offshoring. Energy Econ. 2021, 103, 105562. [Google Scholar] [CrossRef]
  24. Frankel, J.A.; Rose, A.K. Is trade good or bad for the environment? Sorting out the causality. Rev. Econ. Stat. 2005, 87, 85–91. [Google Scholar] [CrossRef] [Green Version]
  25. Wang, F.; Feng, G.F.; Wu, L.H. Regional contribution to the decline of national carbon intensity in China’s economic growth. Econ. Res. J. 2013, 48, 143–155. [Google Scholar]
  26. Chen, B.; Li, J.S.; Wu, X.F.; Han, M.Y.; Zeng, L.; Li, Z.; Chen, G.Q. Global energy flows embodied in international trade: A combination of environmentally extended input–output analysis and complex network analysis. Appl. Energy 2018, 10, 98–107. [Google Scholar] [CrossRef]
  27. Wang, R.Z.; Hao, J.X.; Wang, C.; Tang, X.; Yuan, X. Embodied CO2 emissions and efficiency of the service sector: Evidence from China. J. Clean. Prod. 2020, 247, 119116. [Google Scholar] [CrossRef]
  28. Ma, N.; Li, H.J.; Tang, R.W.; Dong, D.; Shi, J.; Wang, Z. Structural analysis of indirect carbon emissions embodied in intermediate input between Chinese sectors: A complex network approach. Environ. Sci. Pollut. Res. 2021, 26, 17591–17607. [Google Scholar] [CrossRef]
  29. Jiang, M.; Gao, X.; Guan, Q.; Hao, X.; An, F. The structural roles of sectors and their contributions to global carbon emissions: A complex network perspective. J. Clean. Prodn. 2019, 208, 426–435. [Google Scholar] [CrossRef]
  30. Jiang, M.H.; An, H.Z.; Gao, X.Y.; Liu, S.; Xi, X. Factors driving global carbon emissions: A complex network perspective. Resour. Conserv. Recycl. 2019, 146, 431–440. [Google Scholar] [CrossRef]
  31. Moran, D.; Wood, R. Convergence between the EORA, WIOD, EXIOBASE, and opener’s consumption-based carbon accounts. Econ. Syst. Res. 2014, 26, 245–261. [Google Scholar] [CrossRef]
  32. Zheng, H.; Bai, Y.; Wei, W.; Meng, J.; Zhang, Z.; Song, M.; Guan, D. Chinese provincial multi-regional input–output database for 2012, 2015, and 2017. Sci. Data 2021, 8, 224. [Google Scholar] [CrossRef] [PubMed]
  33. Zheng, H.; Zhang, Z.; Wei, W.; Song, M.; Dietzenbacher, E.; Wang, X.; Meng, J.; Shan, Y.; Qu, J.; Guan, D. Regional determinants of China’s consumption-based emissions in the economic transition. Environ. Res. Lett 2020, 15, 074001. [Google Scholar] [CrossRef]
  34. Hou, H.M.; Bai, H.T.; Ji, Y.J.; Wang, Y.; Xu, H. A historical time series for interindustrial embodied carbon transfers within China. J. Clean. Prod. 2020, 264, 121738. [Google Scholar] [CrossRef]
  35. Leontief, W. Quantitative input and output relations in the economic system of the United States. Rev. Econ. Stat. 1936, 18, 105–125. [Google Scholar] [CrossRef] [Green Version]
  36. Lv, K.J.; Feng, X.; Kelly, S.; Zhu, L.; Beng, M. A study on embodied carbon transfer at the provincial level of China from a social network perspective. J. Clean. Prod. 2019, 225, 1089–1104. [Google Scholar] [CrossRef]
  37. Chen, W.; Wu, F.Y.; Geng, W.X.; Yu, G. Carbon emissions in China’s industrial sectors. Resour. Conservat. Recycl. 2016, 117, 264–273. [Google Scholar] [CrossRef]
  38. Li, D.Z.; Huang, G.Y.; Zhang, G.M.; Wang, J. Driving factors of total carbon emissions from the construction industry in Jiangsu Province, China. J. Clean. Prod. 2020, 276, 123179. [Google Scholar] [CrossRef]
  39. An, Q.E.; An, H.Z.; Wang, L. Analysis of the embodied energy flow network between Chinese industries. J. Syst. Eng. 2014, 29, 754–762. [Google Scholar]
  40. Du, Q.; Xu, Y.D.; Wu, M.; Sun, Q.; Bai, L.; Yu, M. A network analysis of indirect carbon emission flows among different industries in China. Environ. Sci. Pollut. Res. 2018, 25, 24469–24487. [Google Scholar] [CrossRef]
  41. Barrat, A.; Barthelemy, M.; Pastor-Satorras, R.; Vespignani, A. The architecture of complex weighted networks. Proc. Natl. Acad. Sci. USA 2004, 101, 3747–3752. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Wasserman, S.; Faust, K. Social Network Analysis: Methods and Applications. Riv. Ital. Sci. Politica 1994, 25, 582–584. [Google Scholar]
  43. White, H.C.; Boorman, S.A.; Breiger, R.L. Social structure from multiple networks. I. Block models of roles and positions. Am. J. Soc. 1976, 81, 730–780. [Google Scholar] [CrossRef]
  44. Doreian, P.; Batagelj, V.; Ferligoj, A. Generalized blockmodeling of two-mode network data. Soc. Netw. 2004, 26, 29–53. [Google Scholar] [CrossRef] [Green Version]
  45. Shen-Orr, S.S.; Milo, R.; Mangan, S.; Alon, U. Network motifs in the transcriptional regulation network of Escherichia coli. Nat. Genet. 2002, 31, 64–68. [Google Scholar] [CrossRef]
  46. Wu, C.Y.; Huang, X.J.; Chai, X.W. Analysis of industrial structure adjustment and carbon reduction potential in Jiangsu province: Based on EIO-LCA model. China Populat. Resour. Environ. 2015, 25, 43–51. [Google Scholar]
  47. Zhang, Y.G. Economic development pattern change impact on China’s carbon intensity. Econ. Res. J. 2010, 4, 120–133. [Google Scholar]
Figure 1. Carbon emission trends.
Figure 1. Carbon emission trends.
Ijerph 20 02603 g001
Figure 2. Technical framework and model selection for research.
Figure 2. Technical framework and model selection for research.
Ijerph 20 02603 g002
Figure 3. Relationship between nodes and edges in the ECETN model.
Figure 3. Relationship between nodes and edges in the ECETN model.
Ijerph 20 02603 g003
Figure 4. Direct carbon emissions and direct carbon emission intensity due to fossil energy consumption of various industrial sectors in Hei-Ji-Liao. (a) Direct carbon emissions of various industrial sectors in Hei-Ji-Liao. (b) Direct carbon emission intensity of various industrial sectors in Hei-Ji-Liao.
Figure 4. Direct carbon emissions and direct carbon emission intensity due to fossil energy consumption of various industrial sectors in Hei-Ji-Liao. (a) Direct carbon emissions of various industrial sectors in Hei-Ji-Liao. (b) Direct carbon emission intensity of various industrial sectors in Hei-Ji-Liao.
Ijerph 20 02603 g004
Figure 5. Embodied carbon emissions and Embodied carbon emission intensity due to fossil energy consumption of various industrial sectors in Hei-Ji-Liao. (a) Embodied carbon emissions of various industrial sectors in Hei-Ji-Liao. (b) Embodied carbon emission intensity of various industrial sectors in Hei-Ji-Liao.
Figure 5. Embodied carbon emissions and Embodied carbon emission intensity due to fossil energy consumption of various industrial sectors in Hei-Ji-Liao. (a) Embodied carbon emissions of various industrial sectors in Hei-Ji-Liao. (b) Embodied carbon emission intensity of various industrial sectors in Hei-Ji-Liao.
Ijerph 20 02603 g005
Figure 6. Graph of Embodied carbon emissions flows between industry sectors in Liaoning. (a) Embodied carbon emissions network; (b) Embodied carbon emissions ring.
Figure 6. Graph of Embodied carbon emissions flows between industry sectors in Liaoning. (a) Embodied carbon emissions network; (b) Embodied carbon emissions ring.
Ijerph 20 02603 g006
Figure 7. Graph of Embodied carbon emissions flows between industry sectors in Jilin. (a) Embodied carbon emissions network; (b) Embodied carbon emissions ring.
Figure 7. Graph of Embodied carbon emissions flows between industry sectors in Jilin. (a) Embodied carbon emissions network; (b) Embodied carbon emissions ring.
Ijerph 20 02603 g007
Figure 8. Graph of Embodied carbon emissions flows between industry sectors in Heilongjiang. (a) Embodied carbon emissions network; (b) Embodied carbon emissions ring.
Figure 8. Graph of Embodied carbon emissions flows between industry sectors in Heilongjiang. (a) Embodied carbon emissions network; (b) Embodied carbon emissions ring.
Ijerph 20 02603 g008
Figure 9. K-cores of the ECETNs in the three provinces of Hei-Ji-Liao. (a) K-core for ECETN in Liaoning; (b) K-core for ECETN in Jilin; (c) K-core for ECETN in Heilongjiang.
Figure 9. K-cores of the ECETNs in the three provinces of Hei-Ji-Liao. (a) K-core for ECETN in Liaoning; (b) K-core for ECETN in Jilin; (c) K-core for ECETN in Heilongjiang.
Ijerph 20 02603 g009
Figure 10. Coalescent subgroups of the ECETNs in Hei-Ji-Liao. (a) Coalescent subgroups for ECETN in Liaoning; (b) Coalescent subgroups for ECETN in Jilin; (c) Coalescent subgroups for ECETN in Heilongjiang.
Figure 10. Coalescent subgroups of the ECETNs in Hei-Ji-Liao. (a) Coalescent subgroups for ECETN in Liaoning; (b) Coalescent subgroups for ECETN in Jilin; (c) Coalescent subgroups for ECETN in Heilongjiang.
Ijerph 20 02603 g010
Figure 11. Diffusivity test for alternative key carbon-emitting industry sectors in Hei-Ji-Liao. (a) Diffusivity test in Liaoning; (b) Diffusivity test in Jilin; (c) Diffusivity test in Heilongjiang.
Figure 11. Diffusivity test for alternative key carbon-emitting industry sectors in Hei-Ji-Liao. (a) Diffusivity test in Liaoning; (b) Diffusivity test in Jilin; (c) Diffusivity test in Heilongjiang.
Ijerph 20 02603 g011
Table 1. Industry Sector Classification.
Table 1. Industry Sector Classification.
CodeIndustrial SectorCodeIndustrial SectorCodeIndustrial Sector
C1Agriculture, forestry and fish farmingC11Petroleum Processing & CokingC21Instrument and cultural office machinery manufacturing
C2Coal Mining & WashingC12Chemical industryC22Craft & Other Manufacturing
C3Oil and gas extractionC13Non-metallic mineral productsC23Scrap and waste
C4Metal miningC14Metal Smelting & RollingC24Production and supply of electricity and heat
C5Non-metallic and other minerals miningC15Metalwork industryC25Gas production and supply industry
C6Food manufacture and tobacco processingC16General Equipment ManufactureC26Water production and supply industry
C7Textile industryC17Specialized equipment manufacturingC27Building industry
C8Textile, clothing, footwear, hats, leather and down and their productsC18Transport equipment manufacturingC28Transport, storage and postal
C9Wood processing and furniture manufacturingC19Electrical machinery and equipment manufacturingC29Wholesale, retail, accommodation and catering
C10Paper, printing and sporting goods manufacturingC20Communications equipment and other electronic equipment manufacturingC30Other services
Table 2. Results for structural characteristic indicators of ECETNs.
Table 2. Results for structural characteristic indicators of ECETNs.
ProvinceNo. of NodeNo. of EdgeAverage DegreeNetwork DensityClustering CoefficientAverage Path Length
Liaoning3057519.170.6610.8731.36
Jilin3058519.50.6720.9021.352
Heilongjiang3044016.6670.5060.841.411
Table 3. Distribution of the values of centrality of 15 industrial sectors in Liaoning Province in 2017.
Table 3. Distribution of the values of centrality of 15 industrial sectors in Liaoning Province in 2017.
Industry SectorsDegree CentralityBetweenness CentralityCloseness Centrality
Gas production and supply industry (C25)36.595.7537.96
Chemical industry (C12)26.831.6136.28
Metalwork industry (C15)60.9816.5142.27
Non-metallic mineral products (C13)36.591.7936.94
Wholesale, retail, accommodation and catering (C29)41.462.4438.68
Transport, storage and postal (C28)46.347.4239.42
Specialized equipment manufacturing (C17)73.1710.2544.57
Metal mining (C4)43.902.7838.68
Production and supply of electricity and heat (C24)31.712.5136.61
Oil and gas extraction (C3)46.343.2839.05
Food manufacture and tobacco processing (C6)34.151.9037.62
Electrical machinery and equipment manufacturing (C19)31.711.9136.94
Coal Mining & Washing (C2)39.022.3538.32
Communications equipment and other electronic equipment manufacturing (C20)34.152.2737.27
Non-metallic and other minerals mining (C5)36.592.2337.62
Table 4. Distribution of the values of centrality of 18 industrial sectors in Jilin Province in 2017.
Table 4. Distribution of the values of centrality of 18 industrial sectors in Jilin Province in 2017.
Industry SectorsDegree CentralityBetweenness CentralityCloseness Centrality
Chemical industry (C12)17.073.2733.07
Non-metallic mineral products (C13)31.716.4336.94
Building industry (C27)26.835.9334.17
Coal Mining & Washing (C2)39.029.1837.27
Metal mining (C4)31.718.7935.97
Non-metallic and other minerals mining (C5)24.393.6535.35
Petroleum Processing & Coking (C11)17.072.2432.28
Metal Smelting & Rolling (C14)14.633.2134.17
Metalwork industry (C15)39.028.1137.62
General Equipment Manufacture (C16)41.464.8937.27
Specialized equipment manufacturing (C17)46.347.3039.42
Craft & Other Manufacturing (C22)31.713.4635.97
Scrap and waste (C23)34.154.6736.61
Production and supply of electricity and heat (C24)31.712.3235.65
Water production and supply industry (C26)31.715.2735.65
Agriculture, forestry and fish farming (C1)36.592.2137.27
Transport equipment manufacturing (C18)29.273.4635.35
Transport, storage and postal (C28)34.155.2336.94
Table 5. Distribution of the values of centrality of 10 industrial sectors in Heilongjiang Province in 2017.
Table 5. Distribution of the values of centrality of 10 industrial sectors in Heilongjiang Province in 2017.
Industry SectorsDegree CentralityBetweenness CentralityCloseness Centrality
Production and supply of electricity and heat (C24)26.831.6956.94
Building industry (C27)41.462.9663.08
Petroleum Processing & Coking (C11)51.226.0367.21
Non-metallic mineral products (C13)31.711.7758.57
Metal Smelting & Rolling (C14)39.022.5161.19
Metalwork industry (C15)95.1235.4295.35
General Equipment Manufacture (C16)58.547.2770.69
Wood processing and furniture manufacturing (C9)31.711.8659.42
Craft & Other Manufacturing (C22)48.784.6666.13
Agriculture, forestry and fish farming (C1)39.027.0762.12
Table 6. Top 15 industry sectors in KATZ index distribution.
Table 6. Top 15 industry sectors in KATZ index distribution.
Liaoning ProvinceJilin ProvinceHeilongjiang Province
CodeValueCodeValueCodeValue
C270.84C121.03C241.15
C280.91C130.64C270.61
C120.65C270.66C140.65
C240.76C240.72C10.71
C140.66C140.56C151.78
C111.49C51.21C281.24
C151.61C21.65C291.84
C40.86C250.56C130.66
C170.65C40.60C160.71
C190.83C280.81C120.69
C301.07C230.56C170.64
C11.36C150.97C30.94
C290.65C261.19C221.47
C180.80C10.70C191.11
C20.73C160.58C100.87
Note: Calculated from the data of the three-province assignment adjacency matrix.
Table 7. Primary selection of key carbon-emitting industry sectors in Hei-Ji-Liao.
Table 7. Primary selection of key carbon-emitting industry sectors in Hei-Ji-Liao.
ProvinceStrong Center Degree IndustriesHigh Impact IndustriesPriming Key Carbon-Emitting Industries
Liaoning(C2, C3, C4, C5, C6, C12, C13, C15, C17, C19, C20, C24, C25, C28, C29), 15(C1, C2, C4, C11, C12, C14, C15, C17, C18, C19, C24, C27, C28, C29, C30), 15(C2, C4, C12, C15, C17, C19, C24, C28, C29), 9
Jilin(C1, C2, C4, C5, C11, C12, C13, C14, C15, C16, C17, C18, C22, C23, C24, C26, C27, C28), 18(C1, C2, C4, C5, C12, C13, C14, C15, C16, C23, C24, C25, C26, C27, C28), 15(C1, C2, C4, C5, C12, C13, C14, C15, C16, C24, C26, C28), 12
Heilongjiang(C1, C9, C11, C13, C14, C15, C16, C22, C24, C27), 10(C1, C3, C10, C12, C13, C14, C15, C16, C17, C19, C22, C24, C27, C28, C29), 15(C1, C13, C14, C15, C16, C22, C24, C27), 9
Table 8. Alternative key carbon-emitting industry sectors’ clustering coefficients in Hei-Ji-Liao provinces.
Table 8. Alternative key carbon-emitting industry sectors’ clustering coefficients in Hei-Ji-Liao provinces.
ProvinceLiaoning
CodeC2C12C15C17C19C28C29
Clustering coefficient0.640.620.640.850.640.620.64
ProvinceJilin
CodeC2C5C12C13C14C15C16C24C28
Clustering coefficient0.640.640.630.620.640.650640.630.65
ProvinceHeilongjiang
CodeC1C13C14C15C16
Clustering coefficient0.770.480.480.480.48
Table 9. Description of the structure of the cohesive subgroup in Hei-Ji-Liao provinces.
Table 9. Description of the structure of the cohesive subgroup in Hei-Ji-Liao provinces.
Cohesive Subgroup No.Industrial Sector No.No. of Industrial SectorsInternal Embodied Carbon Flow (Million Tons)Embodied Carbon Emissions Outflow (Million Tons)Embodied Carbon Emissions Inflow (Million Tons)
Liaoning
IC1, C2, C6, C7, C11, C12, C14, C248152,223.32322,586.00250,336.40
IIC3, C4, C5, C8, C9, C10, C13, C15, C16, C17, C18, C19, C20, C21, C22, C23, C25, C26, C27, C28, C29, C3022122,943.29467,400.10539,210.30
Jilin
IC3, C7, C8, C9, C10, C19, C20, C21, C25919,791.0198,883.7026,735.80
IIC1, C6, C18, C28, C29, C30638,268.0345,572.60244,693.00
IIIC2, C4, C5, C11, C12, C13, C14, C15, C16, C17, C22, C23, C24, C26, C271558,576.96319,424.80169,809.60
Heilongjiang
IC2, C3, C4, C5, C7, C8, C9, C10, C18, C20, C21, C23, C25, C261432,802.30468,315.70145,062.50
IIC1, C6, C12, C28, C29, C306147,400.10136,213.401,150,714.80
IIIC11, C13, C14, C15, C16, C19, C22, C249127,540.40986,130.80122,174.40
IVC1716140.2075,557.1014,362.20
Note: red markers are alternative key carbon-emitting industrial sectors.
Table 10. Key Embodied carbon emission pathways in Hei-Ji-Liao.
Table 10. Key Embodied carbon emission pathways in Hei-Ji-Liao.
ProvinceCritical PathsFlow Volume (104 t)Proportion
LiaoningC2-C14-C11-C13-C30297,1813.2
C12-C27-C19-C18-C28-C30249,3042.68
C15-C24-C11-C27380,1534.09
JilinC2-C24-C29-C6-C1864,8342.76
C5-C17-C27-C2926,4981.13
C24-C13-C27-C6-C1874,3413.17
C15-C14-C12-C29-C3042,7531.82
HeilongjiangC14-C15-C12-C27-C28-C1-C30244,1552.63
C16-C19-C22-C29-C3078,6441,04
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Shan, S.; Li, Y.; Zhang, Z.; Zhu, W.; Zhang, T. Identification of Key Carbon Emission Industries and Emission Reduction Control Based on Complex Network of Embodied Carbon Emission Transfers: The Case of Hei-Ji-Liao, China. Int. J. Environ. Res. Public Health 2023, 20, 2603. https://doi.org/10.3390/ijerph20032603

AMA Style

Shan S, Li Y, Zhang Z, Zhu W, Zhang T. Identification of Key Carbon Emission Industries and Emission Reduction Control Based on Complex Network of Embodied Carbon Emission Transfers: The Case of Hei-Ji-Liao, China. International Journal of Environmental Research and Public Health. 2023; 20(3):2603. https://doi.org/10.3390/ijerph20032603

Chicago/Turabian Style

Shan, Shaonan, Yulong Li, Zicheng Zhang, Wei Zhu, and Tingting Zhang. 2023. "Identification of Key Carbon Emission Industries and Emission Reduction Control Based on Complex Network of Embodied Carbon Emission Transfers: The Case of Hei-Ji-Liao, China" International Journal of Environmental Research and Public Health 20, no. 3: 2603. https://doi.org/10.3390/ijerph20032603

APA Style

Shan, S., Li, Y., Zhang, Z., Zhu, W., & Zhang, T. (2023). Identification of Key Carbon Emission Industries and Emission Reduction Control Based on Complex Network of Embodied Carbon Emission Transfers: The Case of Hei-Ji-Liao, China. International Journal of Environmental Research and Public Health, 20(3), 2603. https://doi.org/10.3390/ijerph20032603

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