1. Introduction
Carbon dioxide emissions are one of the most important greenhouse gas emissions caused by human economic activities [
1,
2]. Many scholars have studied carbon emissions, aiming to find the essential drivers of carbon emissions and reasonably allocate carbon emission reduction responsibilities [
3,
4,
5]. Along with the fact that economic globalization is deepening, the world economy has become a complicated and integrated system. There is still a big challenge to locate the vital targets for carbon emission reduction in such a complicated economic system, especially based on the global supply chain perspective.
In the existing literature, when assessing the impacts of human economic activities on the environment, scholars mostly focus on accounting for embodied carbon emissions rather than territorial carbon emissions [
6,
7,
8]. This is because embodied carbon emissions include all information along the entire supply chain, including production, transportation and consumption, which can be helpful for locating the essential drivers of the global carbon emission growth [
9,
10]. Input-output analysis is the major method of accounting for embodied carbon emissions [
11,
12]. There are also scholars who have applied life cycle assessment (LCA) to study the flow of carbon emissions along the microcosmic production chains [
13,
14,
15]. These studies have provided a strong foundation for accounting for embodied carbon emissions and finding the essential carbon emission drivers. However, the research above mostly focused on the direct relationships when determining the significant carbon emission contributors. The major measurements for locating the important drivers of the carbon emissions were commonly the total embodied carbon emissions of industrial sectors (sectors) or countries, and the exported carbon emissions or imported carbon emissions embodied in the bilateral trades between sectors or countries [
16,
17,
18,
19]. The question of how the influence will be spread among all sectors if the production activities of these sectors are changed to reduce carbon emissions still remains unanswered.
In general, the production activities of one sector are not only decided by that sector but also depend on the inputs from other sectors. The flows of the goods and services construct the economic relationships between sectors and provide the possibilities of the spreading of the influence brought by the change of some sectors. Recently, some scholars have proposed that the industrial chain should be regarded as an interactive and complicated system [
20,
21]. The analysis of the carbon emissions embodied in economic activities should not only focus on two ends of the supply chain but also include the interaction within the intermediate production process between economies [
22,
23,
24]. In most existing research, this intermediate production process was still treated as a black box. Within intermediate production processes, there are complicated linkages between sectors. The change of the production activities of sectors might influence each other. Due to the existence of the supply chain, this kind of influence might further spread to the downstream sectors [
25,
26,
27]. This means that when adjusting the production activities of sectors for reducing carbon emissions, the production and carbon emissions of the related sectors are also influenced. Hence, in the complicated supply chain system, it is important to find the key sectors which have strong capabilities to spread influence.
To describe and analyze the structure of the complicated supply chain, many scholars have combined the traditional input-output analysis with the complex network approach and have proposed the concept of the supply chain network [
27,
28,
29,
30]. Many scholars have studied the complicated characteristics of the supply chain from different perspectives [
31,
32,
33]. By accounting for the embodied carbon emission flows between sectors, scholars found the different roles of sectors in the supply chain networks based on a systemic perspective [
34]. Furthermore, scholars found that some roles of sectors in the supply chain networks would significantly contribute to carbon emissions [
35]. However, few studies focused the potential mechanism on how the structure of the supply chain network influences the carbon emissions. In previous studies, the susceptible–infective (SI) model and its derived models have been widely applied to simulate propagation behavior in the complex network, including the propagation of rumor, virus and traffic congestion [
36,
37,
38]. It has been proven by many scholars that the SI model has the ability to study dynamic mechanisms in the networks. As mentioned above, the economic relationships between sectors make it possible to spread the change of embodied carbon emissions of one sector to others. It is consistent with the concept of the propagation of a virus. Thus, this paper proposes the influence spreading model based on the typical susceptible and infective (SI) model to study the possible cascading influence brought by the change of some individual sectors in the global supply chain system.
The major aim of this paper is to locate the principal sectors for carbon emission reduction in the global supply chain network. Different from existing studies, this paper identifies the key sectors as the policy targets according to their spreading capabilities in global carbon emission flows. This approach can improve the effects of production activity adjustments of sectors and save economic costs for reducing carbon emissions. First, embodied carbon emissions between sectors are accounted for based on EXIOBASE data. Then, this paper applies the complex network theory to build the global embodied carbon emission flow network and quantify the different roles of sectors in the network. Finally, the spreading capabilities of sectors with different roles in the global supply chain are simulated in a proposed influence spreading model.
2. Methods and Materials
2.1. Embodied Carbon Emission Accounting
According to previous research [
39,
40], the carbon emissions embodied in the output of a production process are defined as the sum of the carbon emissions embodied in the intermediate production process and the direct carbon emissions, shown as:
denotes the direct carbon emissions produced by sector . represents the intermediate input from sector to sector . denotes the total output of sector . is the carbon emissions embodied in one unit of inputs from sector , while is the carbon emissions embodied in one unit of outputs of sector .
The equation can be rewritten as the matrix form as follows [
41]:
represents the column vector of the direct carbon emissions. is the intermediate matrix. denotes the diagonal matrix of the total output. is the column vector of the embodied carbon emission coefficient.
Right-multiplying this equation with
and replacing
with
, the embodied carbon emission coefficient could be obtained [
41]:
The embodied carbon emission flows across the production activities are thus obtained as [
42]:
The unit of embodied carbon emissions is kilogram (kg).
2.2. Global Embodied Carbon Emission Flow Network Construction
In this paper, the sectors are taken as the nodes, and embodied carbon emissions between sectors are taken as the edges. The weights of edges are the volume of the embodied carbon emissions. The aim of this paper is to study the influence spreading between different sectors, thus, the self-loops within the sectors are not included. Then, a directed and weighted network are obtained, which is called the global embodied carbon emission flow network.
In the traditional studies regarding embodied carbon emissions, scholars usually use the domestic consumption, imports and exports of embodied carbon emissions to quantify the responsibilities of sectors or countries. In this paper, based on a network approach, four network indicators are applied to measure the importance of sectors. The strength-out is chosen to measure the influence degree of one sector on others in the traditional way. Three other special network indicators, including closeness-out, betweenness and clustering coefficient, are also chosen to quantify the spreading capabilities of sectors from the different network structure perspectives, so that comparisons could be made with traditional research.
2.2.1. Strength-Out
In a network, the strength-out of one node denotes the sums of the weights of all edges from this node to others. In this paper, the strength-out of one sector measures how many embodied carbon emissions would flow to others from this sector. Hence, the strength-out is defined as the direct spreading strength of the sectors. It is calculated as follows [
35]:
where
denotes the set of nodes which connect to node
.
represents the weights of edges from node
to node
.
2.2.2. Closeness-Out
In a network, the closeness-out of one node is defined as the inverse of the average shortest path length to all other nodes in the network. In this paper, the closeness-out of a sector measures how many sectors will be passed at least to access every other sector from this sector. Thus, the closeness-out of one sector is defined as its global spreading capability [
43]. It is calculated as follows:
where
denotes the number of nodes.
represents the shortest path from node
to node
.
2.2.3. Betweenness
In a network, the betweenness of one node denotes the number of shortest paths which pass through it. In this paper, the betweenness of one sector measures how many shortest paths between two other sectors would pass through these sectors. Thus, the sectors with high betweenness play important roles as bridges between sectors. This is defined as the media capability. The betweenness is obtained as follows [
44]:
where
denotes the number of shortest paths from node
to node
that pass through node
.
denotes the number of the shortest paths from node
to node
.
2.2.4. Clustering Coefficient
In a network, the clustering coefficient of one node denotes the connected degree between its adjacent nodes. In this paper, the clustering coefficient measures the connection degree of one sector’s neighbor sectors. If one sector has a high clustering coefficient, it is more possible that its collaboration partners form an industrial cluster around it. Thus, the clustering coefficient is defined as the clustering capability. For a weighted network, the formula for calculating the clustering coefficient is [
45]:
where
is the strength of node
.
represents the element of the adjacency matrix.
denotes the degree of node
.
is the weight of the edge from node
to node
.
2.3. The Influence Spreading Model
To simulate the influence spreading performance of different sectors in the global embodied carbon emission flow network, in this paper, an influence spreading model is proposed based on the typical susceptible–infected (SI) model. In an SI model, each individual has two states: susceptible and infected. Let
be the number of individuals. At time
, the proportion of these two states are
and
, respectively, and there is always
. Let
be the infection rate. At each time, the probability that susceptible individuals will be infected by each infected individual will be
. Because the number of the infected individuals is
at time
, the number of new infected individuals will be
. Thus, the SI model is represented as follows:
In the typical SI model, there are usually not physical linkages between individuals, for example, virus propagation between persons. Infection depends on the probability to a great extent. Therefore, the infection rate is usually assumed as a constant. However, there are actual economic and technical relationships between sectors in this paper. This represents three differences from the typical SI model. First, only the sectors that have linkages with the infected sectors will be susceptible sectors. Other sectors that do not have linkages with the infected sectors will stay in the third status which is called “untouchable”. Second, whether the sectors will be influenced by other sectors depends on some reasonable factors, not only the probability. Thus, the infection rates will be different between sectors. Third, due to the change of the infected sectors, the infection rate of the susceptible sectors also changes over time. Hence, the influence spreading model is proposed as follows:
is a vector containing all infection rates of susceptible sectors at time .
In the real economy, many factors determine whether one sector will be influenced by another. Here, two major factors are assumed to determine the influence spreading process: the strength of linkages and geographic distance. The flows of goods and services between sectors in the multi-regional input-output table are essentially based on international trade. It is easy to understand that trade would not happen if two countries have a very weak trade relationship and are very far from each other, because the benefit and cost of the trade would not be balanced. Hence, this paper here assumes that if there is a relatively strong connection and close distance between two sectors, they will easily have impact on each other. Hence, for each susceptible sector, the infection rate will be:
For sector , the probability that it would not be influenced by the change in sector is . Sector might be impacted by multiple sectors at time , so the probability that it would not be influenced by all infected sectors is . Thus, the probability that sector would be infected at time is . is the volume of the embodied carbon emissions flowing from the infected sector to the susceptible sector . denotes the geographic distance between the countries of sector and sector . If sector and sector belong to the same country, =0. This paper assumes that the strength of linkages and geographic distance are two equal factors on influence spreading, thus the normalized procedures are used to eliminate the difference of the dimensions. () is a constant. When not considering the impacts of linkage strength and distance, the infection rate will be constant, which is consistent with the typical SI model. In general, due to the existence of the linkages between sectors, there are always possibilities for one sector to spread its influence to another. The difference between different sectors of infection could be assumed as the strength of linkages and geographic distance. Hence, here this paper set = 1.
2.4. Data Source
In this paper, the multi-regional input-output table in 2015 from the EXIOBASE database was selected for empirical study (
https://www.exiobase.eu/, version 3.6, accessed on 20 January 2022). It included 200 products produced and consumed in 44 countries (areas) and 5 aggregate regions. In this paper, 200 products were aggregated into 26 products. Because the geographic distance between aggregate regions is difficult to quantify, 44 countries (areas) were selected as the boundary of the world economy system in this paper. The sectors and countries (areas) were listed in
Table A1 and
Table A2 in
Appendix A respectively. The carbon emission data were also retrieved from the environment account of the EXIOBASE database. The geographic distance data were from the CEPII (
http://www.cepii.fr/, accessed on 18 December 2021).
3. Results
3.1. The Global Embodied Carbon Emission Flow Network
The original global embodied carbon emission flow network is a 1144 × 1144 matrix. It has 1,024,296 edges, and it is nearly a full matrix. As this paper already mentioned in our previous works [
35,
46], in complex network theory, it is difficult to find the essential structure characteristics in a fully connected network. The common and effective solution for dealing with this problem is to set a threshold to filter the relatively unimportant connections in the network. Thus, this paper retained the edges which occupied 99% of the total weights in the final networks after ranking the edges by the weights.
In the final global embodied carbon emission flow network, there are 1093 nodes and 34,176 edges, meaning 51 sectors are insular. The network density is 0.029, indicating the network has a not very well-connected structure. There is a total of 2.687 × 10
13 kg of embodied carbon emissions within the network (3.951 × 10
13 kg total in 2015). According to
Figure 1, the degree, degree-in and degree-out of sectors follow power law distribution
. This indicated the global embodied carbon emission flow network had the significant scale-free characteristic, and the sectors were connected with preferential attachment. This meant there were a few sectors which played more important roles than other sectors in the structure of the global embodied carbon emission flow network.
3.2. The Principal Sectors with Different Spreading Capabilities
The top 10 sectors with different spreading capabilities were chosen as the principal sectors, shown in
Table 1. The results shown there are different important sectors with different roles in the global embodied carbon emission flow network. Many sectors in China had very significant export volumes of embodied carbon emissions, including mining and quarrying (CN11), petroleum, chemicals and non-metallic mineral products (CN15), metal and metal products (CN16), electrical and machinery (CN17), electricity, gas and water (CN20), transport (CN23) and financial intermediation and business activity (CN25). Other sectors, including metal and metal products (IN16) and electricity, gas and water (IN20) in India, and electricity, gas and water (US20) in the US, also exported many embodied carbon emissions. In general, if carbon emissions embodied in the production process of these sectors are reduced, it will have significant influence on the importers of intermediate products.
Similar to the strength-out, some typical heavy industrial sectors also had high closeness-out, including petroleum, chemicals and non-metallic mineral products, and electrical and machinery in China. Although other sectors, such as forestry, etc., in the Netherlands (NL13) and Spain (13), mining and quarrying in Portugal (PT11), petroleum, chemicals and non-metallic mineral products in Russia (RU15), Germany (DE15) and the US (US15), and transport in Russia (RU13) did not export as much embodied carbon emissions as China, the changes in these locations can quickly influence many sectors in the global supply chain. This is because most sectors are located in Europe. The countries in Europe usually have the closer connections with each other than other countries.
The sectors with high clustering coefficients include fish in Lithuania (LT05), dairy and eggs in Croatia (HR06), meat n.e.c. in China (CH08), textiles, leather and wearing apparel in Sweden (SE12) and Finland (FI12), forestry, etc., in Romania (RO13), South Korea (KR13) and Mexico (MX13), and post and telecommunications in Bulgaria (BG24). It can be seen that most sectors were food and light industrial sectors. The products of these sectors were usually produced and consumed in the domestic markets. Hence, these sectors and their partners usually form relatively tight clusters in the global economy system.
There are also some sectors which had high betweenness, including petroleum, chemicals and non-metallic mineral products in China (CH15), Russia (RU15), the US (US15), Germany (DE15) and India (IN15), metal and metal products in China (CH16), electrical and machinery in China (CH17) and Germany (GE17), and transport equipment in Germany (DE18). These sectors usually purchase the raw materials or intermediate products for production, then sell the products to other sectors. Additionally, China, Russia, the US, India and Germany all play very important roles in the international trade. Thus, these sectors had significant bridge effects in the global supply chain.
3.3. The Simulation of the Influence Spreading Process
3.3.1. The Sensitivity Analysis of the Infection Thresholds
To compare the difference between different spreading capabilities of sectors, the top 10 sectors were chosen as infection sources. Then, the influence spreading model was applied to simulate the influence spreading process in the global embodied carbon emission flow network. As this paper mentioned before, there are always possibilities of infection if one sector connects to another by setting = 1. To compare the differences of the spreading capabilities of the selected sectors under different scenarios, this paper set 1001 infection thresholds from 0 to 1. If the infection rate of one susceptible sector is larger or equal to the infection threshold, it will be infected.
According to
Figure 2, the total infection proportion had a significant decrease when the infection threshold increased from 0 to 0.005, then was kept stable when the infection threshold was between 0.005 and 1. This indicated that the influence spreading process of embodied carbon emissions in the network was sensitive when the infection threshold was between 0 and 0.005. When the infection threshold was in this range, the influence spreading process would have significant differences in different scenarios. Therefore, six scenarios were chosen to analyze the influence spreading process of embodied carbon emissions with infection thresholds of 0, 0.001, 0.002, 0.003, 0.004 and 0.005. The reason for the small values of the infection thresholds is because the global embodied carbon emission flow network had significant scale-free characteristics. The majority of carbon emissions were embodied in only a few goods and service flows between sectors. Thus, the proportions of the embodied carbon emissions between most sectors were relatively low. Due to the purpose of the infection threshold, its specific value would not influence the results.
According to
Figure 3, the selection of infection thresholds had significant influence on the influence spreading process of carbon emissions between sectors. The proportion of the infected sectors significantly decreased by about 0.95 to 0.10 on average with the increase of the infection threshold. In addition, the spread speed also decreased with the change of the infection threshold. When the infection threshold was 0, the spread process could be completed after three rounds of infection. When the infection threshold was 0.005, it required about eight rounds of infection.
3.3.2. The Performance of Sectors with Different Spread Capabilities
Compared to the randomly selected sectors, the sectors with high strength-out, closeness-out and betweenness generally had better performance in the selected six scenarios. Especially in the high infection threshold scenarios, the influence spreading process might not happen if randomly selecting some sectors as the policy targets (
Figure 3d–f). This indicated that it is reasonable and necessary to select some important sectors as the policy targets. The special positions of these sectors in the supply chain would be helpful for the influence spreading of the related policies. Among the selected four capabilities of sectors in the global embodied carbon emission flow network,
Figure 3a–f shows that the sectors with high betweenness and closeness-out had good spreading performance. With the increase of the infection threshold, the advantage of these sectors became more and more significant. However, the sectors with high strength-out did not perform as well as the sectors with high betweenness and closeness-out. Although these sectors had important direct relationships with neighbor sectors, the influence might be limited based on the full industrial chain perspective. Otherwise, the spreading process did not happen if the sources had a high clustering effect in the network, according to
Figure 3. This might be because the influence is limited within the clusters rather than spreading to the outside of clusters.
4. Discussion and Conclusions
To locate the policy targets for carbon emission reduction, this paper applied the influence spreading model to compare the spreading capabilities of different sectors with different roles in the global embodied carbon emission flow network. First, the embodied carbon emission flows between sectors were accounted for using multi-regional input-output analysis based on the EXIOBASE data. Then, this paper built the global embodied carbon emission flow network and quantified four-dimensional roles of sectors, including strength-out, closeness-out, betweenness and clustering coefficient, using a complex network approach. Finally, the influence spreading model was proposed by improving the typical susceptible–infected (SI) model to simulate the influence spreading process of different sectors. The main results are as follows.
First, the network analysis showed that the global embodied carbon emission flow network had a significant scale-free characteristic. It indicated that the connections between sectors were not even and random. There were a few important sectors in the network. The results also showed that the roles of different sectors in the network were different. For example, electricity, gas and water in China had the greatest export volume of embodied carbon emissions, but it had a very limited global influence in the network. This means this sector might have a significant influence on its neighbor sectors.
Second, the results of the sensitivity analysis showed that the effective infection thresholds were relatively low, which were between 0 and 0.005. This might be caused by the significant scale-free characteristic of the global embodied carbon emission flow network. The majority of carbon emissions were embodied in only a few goods and service flows between sectors. With the increase of the infection threshold, the proportion of the infected sectors significantly decreased from about 0.95 to 0.10, and the spread speed decreased from about three rounds to eight rounds of infection.
Third, the simulation of the influence spreading process showed that the sectors with special roles in the network surely had better performance when compared with randomly selected sectors. This meant that the structure of the global supply chain could be used for the spreading of policy influence. Although previous studies have proven that some sectors with high energy consumption are responsible for climate change, our results indicated that sectors with important indirect relationships in the supply chain have a significant role in the influence spreading of related policies, including the global spreading capability (closeness-out) and the media capability (betweenness). Change in the sectors with high closeness-out and betweenness might have a better spreading speed and scope, which could be helpful for obtaining the maximized effects with minimal cost.
In summary, although some sectors, such as petroleum, chemicals and non-metallic mineral products in China, electricity, gas and water in India, and electricity, gas and water in the US exported significant embodied carbon emissions, the influences on the total carbon emissions of the entire global supply chain might be relatively limited. If the production activities of these sectors change, only the neighbor sectors would be influenced. Although the sectors including petroleum, chemicals and non-metallic mineral products in Germany, electrical and machinery in China, and transport equipment in Germany did not export significant embodied carbon emissions, changes in production activities would widely influence most sectors in the global supply chain. Hence, from a global perspective, the contribution of the sectors with high global spreading capability and media capability to carbon emissions should be recognized and given attention.