1. Introduction
Actively responding to global climate change has become an urgent issue for all countries [
1]. In December 2015, 196 countries signed the Paris Agreement to jointly promote a plan to curb global warming [
2]. China is the first developing country to propose energy-saving and emission reduction targets. It has committed itself to reducing carbon emission intensity in 2020 to 55–60% of the levels in 2005 [
3]. In November 2014, the Chinese government strived to realize the target of decreasing carbon emission as soon as possible [
4,
5] in accordance with the China–US Joint Statement on Climate Change. Many countries have also successively promulgated their own low-carbon action plans. At present, the situation of carbon emissions is still very serious, and is mainly driven by economic sectors with high energy consumption. Global greenhouse gas emission data show that in 2015, global carbon emissions were 4.91 billion tons. The contributions of various economic sectors were as follows: electricity production: 29%, transportation: 27%, industry: 21%, commercial and residential: 12%, and agriculture: 9% [
6]. The trend of high energy consumption and high emissions has caused tremendous ecological pressure in all countries. Economic activity is one of the most important factors for change in carbon emissions [
7]. Decomposing the contribution of various economic development factors to changes in carbon emissions can provide an effective basis for formulating emission reduction policies, and coordinating the balance between economic growth and ecological protection [
8].
The “Silk Road Economic Belt” and the “21st Century Maritime Silk Road” (One Belt and One Road) are open international initiatives proposed by the Chinese government. Their goal is to achieve shared benefits for all the countries along the route [
9,
10,
11]. Green development is an inevitable choice for the countries along the One Belt and One Road, and an objective requirement for the development of the world economy. At present, most of the countries along the One Belt and One Road are developing countries. The development stage and economic characteristics of the various countries are quite different. In 2016, the amount of carbon emissions in the region accounted for 56.1% of global carbon emissions [
12], and most of the countries had not established clear emission reduction targets and action paths. Therefore, selecting countries along the One Belt and One Road as targets to explore the issue of carbon emission pressures is of great importance for the global response to climate change issues. It is also an important means for the countries along the route to seek more feasible and effective emission reduction pathways and achieve sustainable development [
13].
In recent years, many researchers have used carbon footprint to assess the combined environmental impacts of carbon emissions from human economic activities [
14,
15]. In 1992, Rees [
16] and Wackernagel and Rees [
17,
18] proposed the concept of the ecological footprint. On this basis, Thomas Wiedmann and Jan Minx [
19] proposed the concept of the carbon footprint, which is used to measure the total amount of carbon dioxide emissions directly or indirectly caused by an activity, or the total amount of carbon dioxide accumulated in the product life cycle. At the same time, some scholars have defined the carbon footprint as the productive land area needed to absorb carbon emissions, indicating the ecological footprint of carbon emissions [
20]. Since then, others have analyzed carbon footprints from perspectives such as households [
21], products [
22,
23], industries, and sectors [
24,
25], cities [
26], countries [
27], and so forth. Most of the research focuses on the industrial and regional levels. Brown et al. [
28] took the example of 100 major metropolitan areas in the United States to analyze the carbon emissions from energy consumption in the transportation sector and the construction industry. Hertwich [
29] analyzed the carbon footprint of 73 countries and 14 clustered regions. Based on energy consumption data, Chuai and Li [
30] studied the temporal and spatial changes in the carbon footprint of six regions in Northeast China, North China, and East China, and concluded that the carbon emissions, carbon footprint, and carbon footprint intensity in East China were significantly higher than those in the other regions. Chen [
31] used the Stochastic Impacts by Regression on Population, Affluence and Technology (STIRPAT) model and the partial least squares model to assess Beijing’s carbon footprint of energy consumption from the perspectives of population, urbanization, and technology, and analyzed the impact of different factors on the carbon footprint.
Existing studies tend to focus on the carbon footprint itself and do not correlate the carbon footprint with other environmental factors. The introduction of ecological pressure into carbon footprint research can intuitively reflect the impact of carbon emissions on the ecological environment, which is an issue of notable interest. There is little relevant literature at present, and those papers that do exist are usually characterized by carbon footprint ecological pressure [
32,
33]. In order to seek the source of carbon footprint ecological pressure, it is necessary to decompose the contribution of its driving factors into the various ecological pressure affecting the carbon footprint [
34]. Based on studies of the ecological pressure of the carbon footprint of energy consumption, Zhang et al. used logarithmic mean divisia index (LMDI) to decompose the driving forces of different factors. They found that the promotion of economic factors outweighed the inhibitory effects of other factors, and that the ecological pressure brought about by economic development was the most prominent factor [
35].
At present, factor decomposition is mostly used in the study of carbon emission driving factors, and a flood of related research has been undertaken in this area since the mid-1980s. There are four main types of methods used. The first category comprises those that focus on identifying the factors affecting carbon emissions through structural decomposition techniques (SDA) under the input–output analysis framework, combined with input–output tables and models [
36]. Feng [
37] found that the economic growth of the United States during 1997–2007 was the key factor driving the increase of carbon emissions. The change in the energy structure over this period was an important factor in restraining carbon emission growth. Such studies often rely on input–output tables for non-continuous intertemporal analysis, and the types of impact factors are single, which makes it difficult to examine the impact of multiple types of socioeconomic factors. The second category uses the Computable General Equilibrium (CGE) model and scenario analysis to calculate the impact of different socioeconomic factors and their combinations on carbon emissions, in accordance with general equilibrium conditions. This type of research is also dependent on input–output tables, however it is not straightforward to uncover the underlying continuum of a time series [
38,
39]. The third category comprises the use of econometric models to study the impact of various factors on carbon emissions. There are many such methods, but all of them can reveal only statistical laws, not the intrinsic driving factors [
40,
41]. The fourth category uses the Kaya identity, along with extensions and improvements. They use factorial decomposition techniques (IDA) to calculate the contribution of different factors. This type of method has a variety of decomposition forms and can be selected based on research needs [
42,
43,
44,
45]. The logarithmic mean Divisia index (LMDI) method is such a method. It was proposed by Ang et al. in 1998, who deduced its mathematical properties, as well as summarizing a large number of applications of IDA models. Its advantages include using of diverse forms, simple calculations, and the ability to find intrinsic driving factors [
46,
47,
48,
49]. The method is used in many research fields, such as for the decomposition of energy consumption or carbon emission factors. For example, Vinuya et al. [
50,
51] used the LMDI decomposition model to find that economic and population growth is the largest factor driving carbon emissions. The increase in energy use efficiency and the decline in energy intensity can effectively restrain the increase in carbon emissions due to the impact of GDP and population increase. Ma and Stern improved the IDA model and used the LMDI method to decompose the energy intensity changes from 1980–2003, confirming that technological progress is the most important factor in reducing energy intensity [
52]. Zhang et al. used the LMDI method to find that the rapid growth of carbon emissions in China’s power generation industry is mainly due to the use of coal, and energy efficiency is also an important contributor [
53]. Inglesi-Lotz applied LMDI to study the main factors affecting the changes in CO
2 emissions in South Africa. They found that energy intensity had a negative impact on CO
2 emissions during 2008–2014 [
54]. The use of LMDI by Kopidou et al. showed that the two main drivers of industrial CO
2 emissions and employment are economic growth and resource intensity, and that the optimization of the energy structure is conducive to national emission reduction [
55].
The One Belt and One Road is a shared, win–win, open economic and trade cooperation initiative, covering many countries and regions. At present, most of the research on One Belt and One Road focuses on strategic significance, institutional design [
11,
56], and capacity cooperation [
57,
58,
59]. The primary focus has been how to promote the economic development of each country. However, there are several studies on the green and coordinated development of the international regional economy. The related literature focuses on the OECD [
60,
61], East Asia, ASEAN [
62,
63], the European Union and its sub-regions [
17,
64], and a few international regions, meaning less attention is paid to the coordinated development of the regional economy and environment in the One Belt and One Road area. In fact, most of the countries along the “One Belt and One Road” are developing or underdeveloped countries, many of which are not concerned with the problem of low carbon and emission reduction. In addition, the development stages and economic characteristics of the countries along the route are quite different. Therefore, selecting countries along the One Belt and One Road as the target and exploring the carbon footprint ecological pressures of those countries is of great significance to sustainable development; it can also contribute to global emission reduction targets.
In summary, previous research has mainly focused on the measurement and causes of carbon emissions and carbon footprints. There has been little research on ecological pressure caused by carbon footprints, and it does not directly reflect the impact of carbon emissions on the ecological environment.
The LMDI method has been widely used in carbon emission and carbon footprint research. In previous work looking at causal factors, the choice of factors was not varied, with the main focus always being the country. Explanations of changes to ecological pressure resulting from the global carbon footprint were insufficient. In addition, the scope of the countries included was relatively limited, and the differences in economic characteristics and resource endowments among countries are relatively small, which limits the applicability of the research. By selecting countries along the One Belt and One Road with large differences in economic characteristics and resource endowments—based on an ecological-pressure-based model of energy consumption and the carbon footprint—we adopted the LMDI decomposition method to introduce the driving factors that reflect the links between the global carbon footprint and ecological pressure. Factors include the per capita export goods and services trade volume, and the influence of international trade on GDP. Our approach makes the study of the ecological pressures arising from energy consumption as measured by the carbon footprint more straightforward.
The rest of the paper is organized as follows: In
Section 2, the concepts of the energy consumption carbon footprint (
CFEC) and carbon footprint ecological pressure
(EPcfec) are proposed, and the measurement and factor decomposition model of
EPcfec are explained and constructed.
Section 3 introduces the research scope, data sources, and regional divisions.
Section 4 analyzes the calculation results and spatial evolution of
EPcfec from 1994–2014, and uses the LMDI method to decompose
EPcfec. The conclusions and recommendations are given in
Section 5.
2. Methodology
2.1. Carbon Footprint Measurement Model of Energy Consumption
The carbon footprint is a quantified value that measures carbon emissions. It is based the concept of an ecological footprint and can directly measure the response of a natural system to carbon emissions from human activities [
65]. In this paper, we use a measure of carbon footprint based on energy consumption, mainly related to coal, oil, and natural gas. The carbon footprint is measured using the productive land area needed to store energy consumption (see Equation (1)). According to the data provided by the IPCC (2006), productive land includes forests, grasslands, arable land, gardens, and other agricultural lands, of which forest and grassland carbon reserves account for 93% of the total productive land [
66]:
where
CFEC refers to the carbon footprint of energy consumption (hm
2);
Qi is the energy consumption of the
ith energy (t); and
EFi is the energy carbon emission factor (t/tec), for which the study takes the average from several agencies, such as the EIA, as shown in
Table 1 [
67,
68,
69].
β is the conversion coefficient for converting carbon emissions into land area, for which the adopted WWF’s value is 6.49 t/hm
2 [
70].
2.2. The Econometric Model of EPcfec
The ecological pressure of the carbon footprint in energy consumption (
EPcfec) refers to the pressure of carbon emissions on the natural ecosystem and is the ratio of the carbon footprint to productive land area. The calculation is as follows:
where
CFEC is the carbon footprint as described above; and
Sf,
Sa,
Sc, and
Sp refer to the area of forest, cultivated land, permanent cropland, and permanent ranching land, respectively, with the unit hm
2.
When EPcfec ∈ (0, 1), the productive land can fully absorb the carbon emissions resulting from energy consumption, and the carbon emissions, therefore, exert little pressure on the ecological environment. The smaller the value of EPcfec, the lower the ecological pressure. This is termed the light pressure state.
When EPcfec = 1, the carbon emissions generated by energy consumption are equal to those of being absorbed by productive land; that is, the productive land has reached its capacity to absorb the carbon emissions from energy consumption. This is the balanced state.
When EPcfec ∈ (1, +∞), the productive land cannot fully absorb the carbon emissions, and the ecological environment faces a large amount of pressure from carbon storage. The larger the value of EPcfec, the greater the ecological pressure on the carbon footprint. At this point, it is in the high pressure state.
A schematic diagram of the three pressure states is shown in
Figure 1.
2.3. Decomposition Method of EPcfec
The paper applies the LMDI method to the decomposition of ecological pressure factors influencing the carbon footprint. The model introduces two factors related to the
EPcfec—the influence of international trade on GDP, per capita exports of goods and service trade—reflecting linked changes in the ecological pressure between countries. We decomposed the driving factors into five categories: energy structure, energy intensity, the influence of international trade on GDP, per capita export goods and services trade volume, and population intensity of productive land (see Equation (3)). The variables are defined in
Table 2.
where
C is the carbon emission of non-renewable energy as described above;
Q is the consumption of non-renewable energy;
E is the total energy consumption including renewable energy and non-renewable energy;
Y is the gross domestic product;
EXP is exports of goods and services;
POP is the total population of every country;
H is the productive land area;
ε means coefficient;
EF means carbon emission intensity, which is the
ith energy carbon emission coefficient;
ES means
the energy structure that is the proportion of non-renewable energy in the total energy consumption;
EI means the energy intensity that is the ratio of total energy consumption to GDP;
YE means the influence of international trade on GDP, which is the ratio of GDP to exports of goods and services;
EP means the per capita exports of goods and service trade, which is the ratio of exports of goods and services trade to the population;
PH means the population intensity of productive land, which is the ratio of the population to productive land; and
ε means the reciprocal of the conversion coefficient
β, as shown in
Table 2.
where:
where:
In Equations (4) and (5), ΔEPcfec is the change in EPcfec; ΔEPEF, ΔEPES, ΔEPEI, ΔEPYE, ΔEPEP, ΔEPPH, and ΔEPε denote the change in carbon intensity, energy structure, energy intensity, influence of international trade on GDP, per capita exports of goods and service trade, population intensity of productive land, and the conversion coefficient from the 0th to the Tth year, respectively—these are the contribution values; WT is the weight of influence; D represents the change in EPcfec, DEF, DES, DEI, DYE, DEP, DPH, and Dε are the rates of change in carbon intensity, energy structure, energy intensity, Influence of international trade on GDP, per capita import and export goods and service trade, population intensity of productive land, and the conversion factor from the 0th to the Tth year, respectively. ES0, EI0, YE0, EP0, PH0 and EST, EIT, YET, EPT, PHT denote the energy structure, energy intensity, influence of international trade on GDP, per capita export goods and service trade volume, and population intensity of productive land in the 0th year and Tth year, respectively.
3. Data Sources
According to the bilateral trade agreement of the One Belt and One Road Initiative, combined with the availability of data, 56 countries along the route were selected, as shown in
Table 3. Due to the unavailability of data eight countries were excluded: Armenia, Bosnia and Herzegovina, Hungary, Moldova, Montenegro, Nepal, Serbia, and East Timor. The time period used was 1994–2014.
The information analysis center of the Oak Ridge National Laboratory in the United States Carbon Dioxide Information Analysis Center (CDIAC) provides carbon emission data [
71]. The primary energy sources for the countries studied mainly comprise coal, oil, and natural gas. Compiling data on energy consumption and carbon emissions shows that the carbon emissions of countries along the One Belt and One Road are on the rise. The highest value of the carbon emissions in 2012 was 486.41 × 10
7 t after which a downward trend occurred, showing that various countries began to focus on green and low-carbon development models. Carbon emissions from the energy consumption of coal, oil, and natural gas are also decreasing. From 2000–2012, the carbon emissions from oil and natural gas production have changed by a small amount, and remained within the range 50 × 10
7–70 × 10
7 t. Of these, coal combustion is the main source of carbon emissions and is growing rapidly. In 2012, the carbon consumption of coal consumption reached 251.9 × 10
7 t, which was 1.5 times that of oil and natural gas consumption. The trends are shown in
Figure 2.
Data about the productive land area, population, and GDP was taken from the World Bank public data resource [
72]. Taking into account the integrity and continuity of the data, a small amount of data was processed using interpolation. In order to analyze regional differences, 56 countries along the route were divided into seven regions: Mongolia, Russia, Central Asia, West Asia, North Africa, Central and Eastern Europe, Southeast Asia, South Asia, and East Asia (
Table 3). The spatial layout is shown in
Figure 3.
5. Conclusions
Utilizing data from the countries along the One Belt and One Road from 1994–2014, we measured the ecological pressure of the carbon footprint in energy consumption (EPcefc), and revealed the temporal and spatial dynamic changes EPcefc. Furthermore, by using the LMDI method to decompose its driving factors, contributions to ecological pressure arising from the carbon footprint in energy consumption was analyzed for the entire region and in seven sub-regions of the One Belt and One Road.
In general, the value of the EPcefc along the One Belt and One Road increased annually as a whole. From 1994–2014, the value of the EPcefc of 56 countries along the One Belt and One Road increased annually. The growth rate increased at first and then decreased, which indicates that the overall EPcefc of the countries along One Belt and One Road is still high. However, the growth trend has gradually slowed down. In terms of sub-regions, the region with the most ecological stress on the carbon footprint was found to be Southeast Asia, and the lowest was Central Asia. The proportion of carbon footprint ecological pressure in Central and Eastern Europe and Mongolia was gradually decreasing over the study period. The proportion in Southeast Asia, South Asia, and China was gradually increasing, and Central Asia remains basically unchanged. The per capita export of goods and services, and the population density on productive land contribute to ecological pressure on the carbon footprint. Energy structure, the influence of international trade on GDP, and energy intensity exerted an inhibitory effect on the ecological pressure of the carbon footprint. The former had a greater effect than the latter, while the energy structure has less effect, and positive and negative fluctuations. The effect of overall per capita export of goods and services, is greater than the population intensity of productive land.
From the perspective of the seven sub-regions, the effects of various factors on the EPcefc are different, and the influence of each driving factor on different regions is also quite different. Energy intensity shows inhibition in all regions. During the period of 2010–2014, the inhibitory effects were as follows from small to large: Central Asia, Mongolia, Russia, South Asia, East Asia (China), Central and Eastern Europe, West Asia, North Africa and Southeast Asia. The energy structure had both inhibitory and promoting effects, but the effect was always small. Areas where it caused an inhibitory effect include Mongolia, Russia, Central Asia, West Asia, North Africa, and Central and Eastern Europe. Regions where it was a promoter include South Asia and East Asia. The impact of the Influence of international trade on GDP was generally inhibitory effect. The impact of the two factors of per capita exports of goods and services and the population intensity of productive land showed a promoting effect as a whole, and the effect of the former was greater than the latter. However, in the Mongolia–Russia region, the population intensity of productive land shows a slight increase.
This study mainly discusses issues regarding the ecological pressure of the carbon footprint in energy consumption, and other measurements of the ecological pressure of the carbon footprint from various angles, such as transportation and tourism. The conclusions of the paper can only be explained by the carbon emissions generated by energy consumption; however, this has its limitations. The regional division of land was mainly based on geography and not according to economic characteristics. The generality of the analysis conclusions is affected by this. Additionally, the productive land area conversion coefficient directly used the WWF’s value to measure the carbon footprint in energy consumption without considering the difference in storage capacity of different types of productive land for carbon emissions.