Spatial Heterogeneity of Energy-Related CO 2 Emission Growth Rates around the World and Their Determinants during 1990 – 2014

Understanding the spatial heterogeneity and driving force identification of energy-related CO2 emissions (ECEs) can help build consensus for mitigating CO2 emissions and designing appropriate policies. However, previous studies on ECEs that focus on both the global-regional scale and the interaction of factors have been seldom conducted. In this paper, ECE data from 143 countries from 1990 to 2014 were selected to analyze regional differences in ECE growth rates by using the coefficient of variation. Then a geographical detector was used to analyze the key determinant factors on ECE growth rates around the world and in eight types of regions. The results show that: (1) the ECE growth rate in the Organization for Economic Cooperation and Development (OECD) region is low and tended to decrease, while in the non-OECD region it is high and tended to increase; (2) the coefficient of variation and detection factor of ECE growth rates at a regional scale are higher than those at a global scale; (3) in terms of the key determinant factors, population growth rate, growth rate of per capita GDP, and energy intensity growth rate are the three key determinant factors of ECE growth rates in the OECD region and most of the non-OECD regions such as non-OECD European and Eurasian (NO-EE), Asia (NO-AS), non-OECD Americas (NO-AM). The key determinant factors in the African (NO-AF) region are population growth rates and natural gas carbon intensity growth rates. The key determinant factors of the Middle East (NO-ME) are population growth rate, coal carbon intensity growth rate and per capita GDP growth rate; (4) the determinant power of the detection factor, the population growth rate at the global scale and regional scale is the strongest, showing a significant spatial consistency. The determinant power of per capita GDP growth rate and energy intensity growth rate in the OECD region, respectively, rank second and third, also showing a spatial consistency. However, the carbon intensity growth rates of the three fossil fuels contribute little to the growth rate of ECEs, and their spatial coherence is weak; (5) from the perspective of the interaction of detection factors, six detection factors showed bilinear or non-linear enhancement at a global and a regional scale, and the determinant power of the interaction of factors was significantly enhanced; and (6) from the perspective of ecological detection, the growth rate of carbon intensity and the growth rate of natural gas carbon intensity at the global scale and NO-ME region are significantly stronger than other factors, with a significant difference in the spatial distribution of its incidence. Therefore, the OECD region should continue to reduce the growth of energy intensity, and develop alternative energy resources in the future, while those that are plagued by carbon emissions in non-OECD regions should pay more attention to the positive influence of lower population growth Energies 2017, 10, 367; doi:10.3390/en10030367 www.mdpi.com/journal/energies Energies 2017, 10, 367 2 of 17 rates on reducing the growth rate of energy-related CO2 emissions. Reducing energy intensity growth rates and reducing, fossil energy consumption carbon intensity.


Introduction
In 2014, fossil energy accounted for 82% of the global energy use as the main contributor to global energy consumption [1].Energy-related CO 2 emissions are one of the research hotspots in the field of climate change [2].According to the International Energy Agency (IEA), the global annual average growth rate of energy-related CO 2 emissions fell from 1.96% in 1990-2010 to 1.60% in 2010-2014, but continued to grow at an average annual rate of 1.9% during the 1990-2014 period [1].There is a significant difference in the average annual growth rate of energy-related CO 2 emissions in each region, such as Asia (NO-AS) (5.83%) with the highest average annual growth rate from 1990 to 2014, while the lowest is non-OECD European and Eurasian (NO-EE) (−1.97%).
Energy-related CO 2 emissions are currently being studied in two major ways: (1) identifying the driving forces for energy-related CO 2 emissions, which is critical for regions and countries to design carbon emission reduction policies [3].Economic development is considered to be one of the key determinant factors affecting carbon emissions, but decoupling studies show that their relationship needs to consider them in different regions and different stages of development [4].Later, the introduction of Structural Decomposition Analysis (SDA) and Index Decomposition Analysis (IDA) decomposition methods overcame the limitations of previous studies [5][6][7].The Kaya formula has become one of the main analytical frameworks for studying the changing characteristics of CO 2 emissions and their influencing factors [8].Population, energy intensity, affluence and energy carbon intensity are the key determinant factors [9].On the basis on the Kaya formula, subsequent research on drivers in structure, size, technology and other directions continues to expand and enrich.However, literature on whether these driving forces of CO 2 emissions (ECEs) play an independent role interact is still scarce; (2) spatial heterogeneity analysis on ECEs mostly in China, the United States, the former Soviet Union, the European Union, ASEAN (Association of Southeast Asian Nations) and other typical countries and regions [2,[10][11][12][13].Studies in 33 countries have shown that reducing energy intensity and improving energy efficiency are critical to reducing CO 2 emissions in fast developing countries such as China and India [11,14,15].Increasing energy productivity in the industrial sector contributed most to carbon reduction in China (72%) in 2013, followed by households [3].However, some studies have suggested that technological progress has less contribution to CO 2 emissions in developing countries, such as China, compared to population growth and income levels [16].Technological innovation and energy structure transition are important factors in reducing carbon emissions in Europe [17].In addition, studies have also focused on changes in the main drivers of the same region or country over time.Such as the seasonal differences in carbon emissions from fossil fuels in North America [12,13].The increase in carbon emissions in the United States in 1997-2007 was mainly driven by economic growth, with lower emissions after 2007 mainly due to economic recession and substitution of natural gas for coal [12].During the period of economic growth in the former Soviet Union, the increase in wealth was partially offset by a reduction in energy intensity, which was mainly due to a decline in wealth and a reduction in the share of fossil energy use [18].Growth of CO 2 emissions in Beijing during the period from 1997 to 2010 is mainly due to changes in industrial structure and population growth [19].The short-term causal relationship of CO 2 emissions in the G20 region is discussed in [20].In addition, the study also extends to specific sectors such as lifestyles and consumption patterns [16], urbanization [20], international trade [21], industrial sector, power generation [22], and the transportation industry [20,23], which has greatly enriched the content of Energies 2017, 10, 367 3 of 17 ECE research through extensive empirical and comparative studies within and among countries and broadened the research horizons.In addition, literatures have considered the following aspects: lifestyles and consumption patterns [16], urbanization [20], international trade [21] or industrial sectors [23], electricity production [22], transportation [20,23] and so on, which have greatly enriched ECE studies.However, the literature has mainly focused on a few typical countries and regions, and there is a lack of comparisons of spatial heterogeneity in many regions and countries at the global scale over a long period of time, and a lack of cognition from macro-view and comparisons of the contribution of driving forces to ECEs in the same time and space.In recent years, some studies have paid more and more attention to the study of ECE estimation, prediction and policy evaluation in some regions and countries [24][25][26][27][28].
There are mainly two decomposition approaches for exploring the impact factors on CO 2 emissions, one is the structural decomposition based on input-output methods, such as SDA and Data Envelopment Analysis (DEA) [11,29], and the other is based on index decomposition approaches, such as IDA and Logarithmic Mean Divisia Index (LMDI) [15,30].These approaches have contributed significantly to revealing the impact factors on CO 2 emissions in some regions and countries.However, previous studies on ECEs that focus on both global-regional scale and the interaction of factors have been seldom conducted.Geophysical detectors are widely used to reveal the causes of multi-scale spatial heterogeneity and the interactions between factors, without having to accept a hypothesis of test for homoskedasticity and normal distribution.Therefore, this paper, from the growth rate perspective, based on 143 countries as the basic geographic unit, attempts to reveal the spatial heterogeneity of ECEs growth rate at the global-region multi-scale by geophysical detectors and their determinant factors.

Research Framework
The literature shows that countries' global carbon emissions vary due to the combined effect of population, per capita GDP, energy intensity, CI and other factors, but there is a lack of analysis of multi-scale global regional contrasts.Fossil energy consumption is a large contributor to ECEs, however, nuclear, hydro and other renewable energy sources and other non-fossil energy sources are not considered CO 2 emission sources [1], and the term the energy in this paper only refers to fossil fuel energy, and CI can be further decomposed into CI coal , CI oil , CI gas .From the perspective of growth, combined with the availability of data, this paper selects the growth of the above factors as major indicators, and analyses the spatial heterogeneity of global ECE growth rates by the coefficient of variation.Secondly, the growth rate of ECEs and the impact of the index were clustered into 5 levels and then matched each other.Thirdly, geophysical methods are used to explore the interaction between the key determinant factors and the factors in the multi-scale space from the global-local multi-scale perspective, ECE factor detection, interactive detection and ecological detection analysis.
Finally, the paper discusses the future energy structure transformation and emission reduction of various regions, and puts forward policy recommendations (Figure 1).variation.Secondly, the growth rate of ECEs and the impact of the index were clustered into 5 levels and then matched each other.Thirdly, geophysical methods are used to explore the interaction between the key determinant factors and the factors in the multi-scale space from the global-local multi-scale perspective, ECE factor detection, interactive detection and ecological detection analysis.
Finally, the paper discusses the future energy structure transformation and emission reduction of various regions, and puts forward policy recommendations (Figure 1).

Data Source and Processing
Energy and emissions data are mainly from the International Energy Agency (IEA), while population(P) and GDP data are mainly from the World Bank (WB) (Table 1).The study area includes 143 countries.According to the global regional division by BP and IEA, the world is divided into eight types of regions: OECD Americas (O-AM, including four countries), non-OECD Americas (NO-AM IEA estimates include use of fuels, but exclude emissions from non-energy use of fuels: • CO 2 emissions from fuel combustion, including international marine and aviation bunkers but excluding process emissions.The coefficient of variation has been widely used in the study of geographical spatial heterogeneity, and its advantage is that it can eliminate the influence of units and mean differences on the results.In this paper, we use C v to measure the difference between ECE growth rates in different countries and their degree, and describe the temporal and spatial heterogeneity of carbon emission growth rates.The formula as follows: where C v is coefficient of variation of one regional ECEs growth rate; S is an unbiased estimation for a regional standard for ECEs growth rate; G is the mean of a region ECEs growth rate; G i is the growth rate of energy-related CO 2 emissions in country i in a region; n is the number of countries, which is a ECE region.

The Geographical Detector
Its core idea is that geographical things always exist in a specific spatial position, while the environmental factors affecting the development and changes in the space are different, if the environmental factors and geographical things in space have significant consistency, then this kind of environmental factors play a decisive role in the occurrence and development of geographical matter [32].In recent years, as an important method to detect the cause and determinants of the spatial heterogeneity of a specified element, it has been gradually applied to the study of social, economic and natural problems [33][34][35].Factor detection, interactive detection and ecological detection are three functions of the geographical detector: (1) Factor detector identifies the risk factors that contribute to the risk and is used to reflect the degree to which a determinant explains the prevalence of the growth rate of energy-related CO 2 emissions (G).We rename it the P x,G , defined as follows [32]: where G denotes the growth rate of ECEs; P x,G denotes the explanatory power of the impact factor x on G.A study area is composed of n units and is stratified in h = 1, 2, . . ., L stratum; n and σ 2 denote the area and the variance of G prevalence of the study area, respectively.n h and σ 2 h denote the h layer sample size and variance of G, respectively.The value of P x,G is required to be within [0, 1]: 1 if the determinant completely controls G, 0 if the determinant is completely unrelated to G. The greater the value, the stronger the explanatory power of the impact factor x to the ECE growth rate.A value of zero indicates that the impact factor is completely unrelated to the ECE growth rate, and a value of 1 indicates that the impact factor can fully explain the ECE growth rate distribution differences.(2) The interaction detector identifies the interaction between the ECE growth rate factor x i and the factor x j , explaining the impact factor as an independent function or having interactions [32].
If P(x i ∩ x j ) < min(P(x i ), P(x j )), which shows that the interaction between factor x i and x j is weak and univariate, where the symbol '∩' denotes the intersection between x i and x j .If min(P(x i ), P(x j )) < P(x i ∩ x j ) < max(P(x i ), P(x j )), which indicates that the interaction between factor x i and x j is the unidirectionally weakened.If P(x i ∩ x j ) > max(P(x i ), P(x j )), which indicates that the interaction between factor x i and x j is the bilinear strengthening.If P(x i ∩ x j ) > P(x i ) + P(x j ), which indicates that the interaction between factor x i and x j is nonlinear and Energies 2017, 10, 367 6 of 17 enhance.If P(x i ∩ x j ) = P(x i ) + P(x j ), which indicates that the interaction between factor x i and x j is independent.
(3) The ecological detector mainly identifies the impact differences of two factors.By comparing the differences of total variance of ECE growth indexes among different influencing factors, this paper explores the influence of different influencing factors on the ECE growth rate distribution and the role of whether there are significant differences, measured by F test [32]: where F denotes the F-test value in statistical test, n x i ,p and n x j ,p denote the sample size of impact factor x i and x j in the unit p, respectively, σ 2 x i ,p and σ 2 x j ,p denote the variance of the impact factor x i and x j , respectively, and the statistical expression must followthe distribution x j ,p .The impact factor x i controlling the ECE growth rate is significantly higher than the impact factor x j if the null hypothesis is rejected when p < 0.05.

Regional Growth Characteristics of Energy-Related CO 2 Emissions (ECEs)
The growth rates of global ECEs are classified into five levels by the natural breaks classification method.The difference of ECE growth rates between 1990 and 2014 is significant.The negative and positive increase in ECEs occurred mainly in the "heartland" of Eurasia, where the largest increase in emissions happened in China (7011.0 million tons), the lowest in Russia (−695.7 million tons) (Figure 2a).Higher growth rates are mainly in the low latitudes, such as in Asia, the Middle East, East Africa and West Africa, Latin America and other regions in the low latitudes of the developing countries (Figure 2b).In terms of specific countries, the average annual growth rate of Benin's emissions is the highest (13.83%), and the one of Georgia is the lowest (−5.93%).The reason for the differences in emission growth rates is to a certain extent connected with the development stage and level of countries.x controlling the ECE growth rate is significantly higher than the impact factor j x if the null hypothesis is rejected when p < 0.05.

Global Distribution of Energy-Related CO2 Emission Growth Rates
3.1.1.Regional Growth Characteristics of Energy-Related CO2 Emissions (ECEs) The growth rates of global ECEs are classified into five levels by the natural breaks classification method.The difference of ECE growth rates between 1990 and 2014 is significant.The negative and positive increase in ECEs occurred mainly in the "heartland" of Eurasia, where the largest increase in emissions happened in China (7011.0 million tons), the lowest in Russia (−695.7 million tons) (Figure 2a).Higher growth rates are mainly in the low latitudes, such as in Asia, the Middle East, East Africa and West Africa, Latin America and other regions in the low latitudes of the developing countries (Figure 2b).In terms of specific countries, the average annual growth rate of Benin's emissions is the highest (13.83%), and the one of Georgia is the lowest (−5.93%).The reason for the differences in emission growth rates is to a certain extent connected with the development stage and level of countries. (a)

Regional Characteristics of Impact Factors on the Growth Rate of ECEs
From the global viewpoint, the coal carbon intensity growth rate (2.47%) is the fastest growth of all factors, and EI shows the slowest growth (−0.88%), which indicates that the high demand for coal is mainly for thermal power generation, and construction in some industries, especially in a number of developing countries.From the regional view, the average annual population growth rate in the OECD region is generally lower, and Non-OECD countries, except NO-EE with negative growth, showed more than 1% growth (Table 2).The average annual per capita GDP growth rate in NO-AS is the highest, which the other Non-OECD areas are slightly higher than OECD; the growth rate of EI in NO-AS and NO-EE is the lowest, while on the contrary NO-ME is the highest.The highest average annual growth rate of CI coal appears in NO-ME (up 10.26%), far more than other regions, however, that in O-AM and O-EU region showed negative growth.The average growth rate of CI oil and CI gas in non-OECD countries is higher except, that in NO-EE which is negative, but OECD countries had low or negative growth.NO-AS has become the fastest growing region in CI coal , while NO-EE has become the slowest growing region in CI coal .NO-ME has become the fastest growing region in CI gas .

Regional Differentiation of Energy-related CO 2 Emission Growth Rates
The regional variation coefficients can reflect the degree of regional differences.According to coefficients of regional variation from 1990, 2000, 2010 and 2014, the following main characteristics can be seen (Table 3): (1) from a global perspective, the coefficient of variation of global ECEs is larger than that of the sub-regions, and the coefficient of variation is increasing, which indicates that the differences of emission growth in sub-regions are greater and growing over time; (2) from a regional perspective.Firstly, the coefficient of variation constantly tend to increase in non-OECD regions, especially NO-AS.In contrast, that of the OECD regions, especially O-AO, tend to decrease.Secondly, the coefficients of variation (less than 2) are smaller and smaller, mainly in the OECD regions, such as O-AM, O-EU, O-AO.The main reason is that there is a high degree of consistency and strong stability for ECE growth rates in the OECD countries.Furthermore, there is a high similarity in the influencing factors, development level and development stage in those countries.Thirdly, the coefficient of variation (more than 2) tends to be widened, mainly in non-OECD regions such as NO-EE, NO-AF, NO-AS.This is mainly due to the large differences in the level of development, stage of development and speed of development among the Non-OECD countries, especially with the rapid development of the BRICs and other emerging economies such as China and India in recent years.Finally, the coefficient of variation of NO-ME is lower than that of other non-OECD regions.The main reason is that there is greater similarity in the energy consumption structure, industrial structure, development level and development stage among NO-ME countries.

Factor Detection Analysis
(1) Factor Detection at a Global Scale Factor detection is mainly used to detect the degree of interpretation of the growth rate of the factors on ECE growth rates.The order of contribution to ECE growth rate on a global scale is: population growth rate > the growth rate of per capita GDP > energy intensity growth rate > natural gas carbon intensity growth rate > coal carbon intensity growth rate > petroleum carbon intensity growth rate (Table 4).In this paper, we determine whether the factor is the key determinant factors with P value (=0) and q statistic value (≥0.10), then a key determinant factor in the world includes the population growth rate, The growth rate of per capita GDP, energy intensity growth rate (Table 5).The factor of population growth rate (0.240) has the largest determinative power, indicating that global ECE growth is most strongly controlled by the population growth rate factor, which has the strongest consistency with the ECE growth rates.The reason is that the increase in population leads to a significant increase in energy demand, both in developed countries and developing countries.and the increase in urbanization, population aging, and economic levels has resulted in an increase in ECEs due to the increase in population growth, which explains why although OECD countries have low population growth rates, the population growth rate has become the key determinant factors in the eight types of regions.In addition, it also shows that the current global energy efficiency and the proportion of renewable energy consumption growth may lag behind the increase in population, and it is sometimes difficult to offset the impact on population growth.The sub-important decisive factor is the growth rate of per capita GDP whose determinative power is 0.178, which indicates that the growth rate of per capita GDP is one of the important factors affecting global ECE growth rates.The main reason is that the increase in the growth rate of per capita GDP results in an increase in per capita income, lifestyle, consumption patterns will be further enhanced, with a corresponding increase in energy demand.When the energy consumption structure is still dominated by fossil fuels, it will directly lead to growth of the ECE rates.The contribution of energy intensity growth to global ECE growth ranks third among all factors.The main reason is that there are great differences in development stages among different countries.The values of X 4 , X 5 , and X 6 are generally smaller (less than 0.1) on the global scale, due in part to differences in energy carbon intensity across the world, such as X 4 , X 5 and X 6 , of the decision-making power is partially offset and weakened.
(2) Factor Detection at a Regional Scale The order of contribution to ECE growth rate in O-AM is: population growth rate > the growth rate of per capita GDP > energy intensity growth rate > natural gas carbon intensity growth rate > coal carbon intensity growth rate > petroleum carbon intensity growth rate (Table 4).The key determinant factors on O-AM are population growth rate, the growth rate of per capita GDP and energy intensity growth rate for the reason that with the rapid development of alternative energy, there are differences in the energy consumption structures and fossil energy growth rates among countries (Table 5).The order of contribution to ECE growth rates in NO-AS is: population growth rate > the growth rate of per capita GDP > energy intensity growth rate > coal carbon intensity growth rate > natural gas carbon intensity growth rate > petroleum carbon intensity growth rate.The key determinant factors in NO-AS are: population growth rate, the growth rate of per capita GDP and energy intensity growth rate for the rapid development of renewable energies such as nuclear energy, hydropower, and the difference of energy consumption structure and growth rate of energy consumption.The order of contribution to ECE growth rates in O-EU is: population growth rate > the growth rate of per capita GDP > energy intensity growth rate > petroleum carbon intensity growth rate > coal carbon intensity growth rate > natural gas carbon intensity growth rate.The key determinant factors in the O-EU region are: population growth rate, the growth rate of per capita GDP, energy intensity growth rate because of the rapid development of the nuclear energy, hydropower, wind energy and other renewable energies, and the larger difference of natural gas growth rate among those countries.The order of contribution to ECE growth rates in the O-EE region is: population growth rate > energy intensity growth rate > the growth rate of per capita GDP > petroleum carbon intensity growth rate > natural gas carbon intensity growth rate > coal carbon intensity growth rate.The key determinant factors in O-EE are: population growth rate, energy intensity growth rate, the growth rate of per capita GDP, because of the inconsistency of factors such as the fossil energy consumption structure and growth rate in O-EE countries.The order of contribution to ECE growth rates in NO-AF is: population growth rate > GDP growth rate > energy intensity growth rate > natural gas carbon intensity growth rate > coal carbon intensity growth rate > petroleum carbon intensity growth rate.The key determinant factors in NO-AF are: population growth rate, the growth rate of per capita GDP, energy intensity growth rate, natural gas carbon intensity growth rate, for the main reason that the stability and consistency of the natural gas carbon intensity increase is slightly better than that of the petroleum carbon intensity and coal carbon intensity in NO-AF countries.In addition, compared to the proportion of annual consumption, the proportion of annual carbon emissions and the average annual carbon intensity of coal and oil, those of natural gas in the non-OECD region are usually much higher (Table 2).The order of contribution to ECE growth rates in NO-AS is: population growth rate > the growth rate of per capita GDP > energy intensity growth rate > natural gas carbon intensity growth rate > petroleum carbon intensity growth rate > coal carbon intensity growth rate.The key determinant factors in the NO-AS region are: population growth rate, the growth rate of per capita GDP, energy intensity growth rate, for the great difference of the fossil energy consumption of carbon intensity and structure of energy consumption.The order of contribution to ECEs growth rate in NO-AM is: population growth rate > the growth rate of per capita GDP > energy intensity growth rate > natural gas carbon intensity growth rate > coal carbon intensity growth rate > petroleum carbon intensity growth rate.
The key determinant factors on NO-AM are: population growth rate and per capital GDP growth rate, for the reason of the larger difference of the growth rate of petroleum carbon intensity in NO-AM.The order of contribution to ECE growth rates in NO-ME is: population growth rate > growth rate of carbon intensity of coal > growth rate of per capita GDP > growth rate of energy intensity > growth rate of natural gas carbon intensity > growth rate of petroleum carbon intensity.The key determinant factors in NO-ME is population growth rate, coal carbon intensity growth rate and per capital GDP growth rate, which have a high similarity and consistency in NO-ME.In addition, although the carbon intensity of the Middle East is smaller than that of other factors, the growth rate of the carbon intensity is the largest because coal is used to generate electricity instead of oil in some countries.
(3) Factor Detection from the Factor Itself The decisive force of the regional scale detection factor is generally high, which indicates that the decisive factors contribute more to the smaller spatial scale.From the perspective of detection factors, the factors show a certain consistency and different characteristics in the change of ECE growth rates in different regions.Population growth rate (X 1 ) is more consistent on the global scale and in all regions, ranking the first among all factors (Table 5).This suggests that the increase in population growth over the past 25 years, regardless of the degree of regional economic development, results in an increase in energy demand and an increase in CO 2 emissions.X 1 , X 2 , and X 3 have been shown to be more prominent in the OECD and in some non-OECD regions, especially in the OECD.The detection factors of X 1 , X 2 and X 3 played a more pronounced decisive role in some areas, especially showing a high consistency in OECD.The main reason is just strong similarity for OECD countries in socio-economic development speed, development level, industrialization process, urbanization level, population growth.With the smaller regional differences and the rapid development of renewable energy sources, the difference in the growth rate of energy emissions is relatively small.The determinant power of X 2 and X 3 maintain a certain consistency in the non-OECD region, due to a large difference between X 2 and X 3 in NO-EE and NO-ME.With the larger difference of the energy intensity growth of X 4 , X 5 and X 6 , their contribution to the growth rate of ECEs is weak, which indicates that the differences of ECE growth rates are more affected by factors such as population growth rate, the growth rate of per capita GDP and energy intensity growth rate, which reflect socio-economic, policy and technical factors.The determination factors of X 4 , X 5 and X 6 in non-OECD region are slightly larger than those in the OECD.Note: "-"denotes that the value is less than 0.5 or less.

Interaction Detector Analysis
The interaction detector is used to check whether two factors impacting on the ECEs growth rate work independently or not.The result is as follows: (1) Global scale: the factors X 4 and X 5 were found to enhance each other to increase the growth rate of ECEs, which is a nonlinear relationship.The interactions of factors X 4 and X 6 , and interactions of factors X 5 and X 6 .had a similar effect.The relationship between the other detector factors is bilinear, when the factors enhance each other to increase the decisive power of factors to ECEs growth rate.X 5 and X 6 were found to enhance each other (X 5 ∩X 6 (0.077) > 0.060 = X 5 (0.009) + X 6 (0.051)), which indicates synergistic effect on ECEs is generated by power of determinant (P x,G ) of X 5 ∩X 6 .(2) Regional scale: In the O-AM region, all factors were found to enhance each other, and the relationship of the power of determinant (P x,G ) of X 2 ∩X 6 is nonlinear, which is the same as the power of determinant (P x,G ) of X 5 ∩X 6 , while the relationship of the other factors is bilinear.X 2 and X 6 were found to enhance each other to (X 2 ∩X 6 (0.399) > 0.379 = X 2 (0.342) + X 6 (0.037)), which indicates synergistic effect on ECEs is generated by power of determinant (P x,G ) of X 2 ∩X 6 .
In the O-AO region, all factors were found to enhance each other, and the relationship of power of determinant (P x,G ) of X 1 ∩X 6 is nonlinear, which is the same as power of determinant (P x,G ) of X 5 ∩X 6 .While the relationship of the other factors is bilinear.X 1 and X 6 were found to enhance each other to (X 1 ∩X 6 (0.537) > 0.523 = X 1 (0.506) + X 6 (0.017)), which indicates synergistic effect on ECEs is generated by power of determinant (P x,G ) of X 1 ∩X 6 .
In the O-EU region, all factors were found to enhance each other, and the relationship of power of determinant (P x,G ) of X 1 ∩X 6 is nonlinear, which is the same as the power of determinant (P x,G ) of X 2 ∩X 6 , X 3 ∩X 6 , and X 5 ∩X 6 .The relationship of the other factors is bilinear.X 1 and X 6 were found to enhance each other to (X 1 ∩X 6 (0.541) > 0.537 = X 1 (0.527) + X 6 (0.008)), which indicates synergistic effect on ECEs is generated by power of determinant (P x,G ) of X 1 ∩X 6 .
In the NO-EE region, all factors were found to enhance each other, and the relationship of power of determinant (P x,G ) of X 1 ∩X 4 is nonlinear, which is the same as the power of determinant (P x,G ) of X 1 ∩X 5 , X 2 ∩X 3 , X 2 ∩X 4 , and X 5 ∩X 6 , while the relationship of the other factors is bilinear.X 1 and X 5 were found to enhance each other to (X 1 ∩X 5 (0.486) > 0.478 = X 1 (0.428) + X 5 (0.050)), which indicates synergistic effect on ECEs is generated by power of determinant (P x,G ) of X 1 ∩X 5 .
In the NO-AF region, all factors were found to enhance each other, and the relationship of power of determinant (P x,G ) of X 1 ∩X 5 is nonlinear, which is the same as power of determinant (P x,G ) of X 5 ∩X 2 , X 5 ∩X 3 , X 5 ∩X 4 , and X 5 ∩X 6 .The relationship of the other factors is bilinear.X 1 and X 5 were found to enhance each other to (X 1 ∩X 5 (0.225) > 0.220 = X 1 (0.219) + X 5 (0.001)), which indicates synergistic effect on ECEs is generated by power of determinant (P x,G ) of X 1 ∩X 5 .
In the NO-AS region, all factors were found to enhance each other, and the relationship of power of determinant (P x,G ) of X 1 ∩X 4 is nonlinear, which is the same as power of determinant (P x,G ) of X 1 ∩X 5 , X 2 ∩X 4 , X 2 ∩X 5 , X 3 ∩X 5 , X 4 ∩X 5 , X 4 ∩X 6 , and X 5 ∩X 6 , while the relationship of the other factors is bilinear.X 1 and X 5 were found to enhance each other to (X 1 ∩X 5 (0.430) > 0.425 = X 1 (0.421) + X 5 (0.004)), which indicates synergistic effect on ECEs is generated by power of determinant (P x,G ) of X 1 ∩X 5 .
In the NO-AM region, all factors were found to enhance each other, and the relationship of power of determinant (P x,G ) of X 4 ∩X 5 is nonlinear, which is the same as power of determinant (P x,G ) of X 5 ∩X 6 .The relationship of the other factors is bilinear.X 2 and X 3 were found to enhance each other to (X 2 ∩X 3 (0.142) > 0.177 = X 2 (0.110) + X 3 (0.067)), which indicates synergistic effect on ECEs is generated by power of determinant (P x,G ) of X 2 ∩X 3 .
In the NO-ME region, all factors were found to enhance each other, and the relationship of power of determinant (P x,G ) of X 1 ∩X 5 is nonlinear, which is the same as power of determinant (P x,G ) of X 2 ∩X 3 and X 5 ∩X 6 , while the relationship of the other factors is bilinear.X 1 and X 5 were found to enhance each other to (X 1 ∩X 5 (0.184) > 0.179 = X 1 (0.176) + X 5 (0.002)), which indicates synergistic effect on ECEs is generated by power of determinant (P x,G ) of X 1 ∩X 5 .
We draw the conclusion that there was no weakening effect, single linear weakening effect and independent effect between the six factors.The interaction detection results show that the determinant power of any two factors show bilinear or non-linear enhancements in the study area.That is, under the control of any two detectors, the difference of ECE growth rates will decrease and the determinant power of factors was significantly enhanced.

Ecological Detection Analysis
The results of the ecological detection focused on whether there were significant differences in the effects of detection factors on the spatial distribution of ECEs growth rates.The results show that: (1) on a global scale, the incidence of ECE growth rates between X 5 and X 6 is significantly different at the 95% confidence level; (2) on a regional scale, the incidence of ECE growth rates between X 5 and X 6 is significantly different at the 95% confidence level in the Africa region.There is no statistically significance in the statistics for explanatory power between the detection factors in other regions and no significant difference in spatial distribution, because the explanatory power of other factors are smaller.Compared to other factors, the explanatory power of the growth rate of oil carbon intensity and the growth rate of natural gas carbon intensity factor is much bigger.Therefore, reducing the growth rate of oil and carbon intensity and natural gas carbon intensity is a common problem faced by all regions, especially in Africa.

Conclusions
This paper examines the changes of the global energy-related CO 2 emission growth rates during 1990-2014.We analyze the regional differences in the growth rate of global energy-related CO 2 emissions by the coefficient of variation, and explain the determinant factors of energy-related CO 2 emission growth rates by means of geographical detector methods, and provide a new perspective and methods for relevant carbon dioxide emission research.The results can provide a guide for those regions on how to reduce energy-related CO 2 emissions.Results show that: (1) Overall, the highest and lowest growth rate of global energy-related CO 2 emissions during 1990-2014 were mainly in the Non-OECD regions such as NO-AS, NO-EE and NO-ME.
The growth rate of energy-related CO 2 emissions in Non-OECD regions are higher than in the OECD regions due to the development stage and development level of the countries, especially in the NO-AS.(2) The coefficient of variation of energy-related CO 2 emissions growth rate on the global scale is larger than on the regional scale.As far as the eight types of regions is concerned, the spatial heterogeneity is smaller and tends to decrease in the OECD, while that in non-OECD regions tends to expand, with the spatial heterogeneity being the largest in NO-AS and the smallest in O-AO.(3) The determination factor of the detection factor on the global scale is not as good as that on the sub-scale scale.The main reason is that the differences of the detection factors are large at the global scale, which weakens the decisive force.( 4 (5) From the perspective of detection factors, the determinant power of the factors in the different regions of the energy-related CO 2 emissions growth rate of decisive force showed a specific consistency and difference.The population growth rate is the strongest of the determinant power of the factors in the world and in all regions, showing a significant degree of consistency.
In general, the determinant power of the growth rates of coal carbon intensity, petroleum carbon intensity and natural gas carbon intensity, with great differences in each area, have little effect on the energy-related CO 2 emission growth rates.However, the impact of the determinant power on non-OECD is slightly increased.For example, the growth rate of coal carbon intensity is the second decisive factor in non-OECD Middle East.(6) In the world and in each region, it is significantly enhanced with bilinearity or non-linearity for the determinant power of factors that any two factors interact with each other at the global scale or at the regional scale.(7) The high conversion between oil and natural gas in the process of production and consumption, in addition, the demand for oil and gas growth in the world is more synchronized.Thus, the interaction between the growth rate of oil and gas is more significant than that of other factors.As a result, at the global scale or in the non-OECD Middle East region, the determinant power of X 5 ∩X 6 are significantly stronger than other impact factors on energy-related CO 2 emissions growth rates, with significant differences in their spatial distribution.

Discussion
The future, the proportion and growth rate of fossil energy will be further reduced due to the rapid of renewable energy development, after the Paris Agreement.Based on the above analysis, the determinant power of the factors of energy intensity and energy carbon intensity has less influence on the OECD countries than that of non-OECD ones.Despite the relatively low population growth rate and the slow GDP growth rate in the OECD, these factors contribute most to the growth rate of energy-related CO 2 emissions.Therefore, in improving energy efficiency and increasing the development and utilization of alternative energy sources at the same time, more attention should be paid for the impact of factors population and affluence, for example, raising awareness of low-carbon consumption and low-carbon living habits in OECD countries.For non-OECD regions, the future policy focus should be to reduce the growth rate of energy intensity and to decrease the intensity of fossil energy consumption due to their relatively high determinant power while controlling the growth rate of per capita GDP.
However, this study is still preliminary with some limitations.Firstly, this paper only uses selected six important indicators from viewpoint of growth rate, yet other factors reflecting structure and technology, such as population, industry, renewable energy, are excluded from this study.These factors may also affect the reason of spatial heterogeneity of energy-related CO 2 emission growth rates, but, considering the difficulties of data acquisition for all of the countries, further research in consideration of these issues would contribute to understanding the spatial heterogeneity.Secondly, it is noteworthy that, although the coal consumption in non-OECD Middle East is relatively small, the growth rate of coal carbon intensity is the highest growth rate of all fossil energy carbon intensities, and the coal carbon intensity growth rate has become the area of growth of energy-related CO 2 emissions, second only to the population growth rate, and thus is also worthy of close research attention.Furthermore, this paper only classifies the world into eight types of region according to the IEA, further study on energy-related CO 2 could focus on types of global climate negotiations, types of national development stage, and types of national energy utilization by using the geographical detector method.

Figure 1 .
Figure 1.Conceptual framework of impact factors on energy-related CO2 emissions growth rate.Figure 1. Conceptual framework of impact factors on energy-related CO 2 emissions growth rate.

Figure 1 .
Figure 1.Conceptual framework of impact factors on energy-related CO2 emissions growth rate.Figure 1. Conceptual framework of impact factors on energy-related CO 2 emissions growth rate.

Figure 2 .
Figure 2. Change of global energy-related CO 2 emissions (1990-2014).(a) The average annual change; (b) the average annual growth rate.

Figure 3 .
Figure 3. Energy-related CO2 emissions of change (1990-2014) from energy.(a) The amount of change; (b) Average annual growth rate.Figure 3. Energy-related CO 2 emissions of change (1990-2014) from energy.(a) The amount of change; (b) Average annual growth rate.

Figure 3 .
Figure 3. Energy-related CO2 emissions of change (1990-2014) from energy.(a) The amount of change; (b) Average annual growth rate.Figure 3. Energy-related CO 2 emissions of change (1990-2014) from energy.(a) The amount of change; (b) Average annual growth rate.

Table 1 .
Index system for energy-related CO 2 emissions (ECEs) growth rate.

Table 2 .
Volume of variation and average annual growth rate of energy CO 2 emissions and their factors in the world from 1990 to 2014.

Table 3 .
Coefficient of variation of energy-related CO 2 emissions growth rate in the eight types of regions.

Table 4 .
The determinant power of impact factors on the energy-related CO 2 emissions growth rate in the eight types of regions and the world from 1990 to 2014. the stronger the evidence against the null hypothesis.X 1 denotes Population growth rate; X 2 denotes per-capital Gross Domestic Product growth rate; X 3 denotes energy intensity of the economy growth rate; X 4 denotes carbon intensity of coal consumption growth rate; X 5 denotes carbon intensity of oil consumption growth rate; X 6 denotes carbon intensity of gas consumption growth rate.

Table 5 .
Proportion and rank of decision force of detector factors in individual regions.
) At the global scale, the key determinant factors on energy-related CO 2 emission growth rates are population growth rate, the growth rate of per capita GDP, and energy intensity growth rate.At the regional scale, determinant power of the key determinant factors on the growth rate of energy-related CO 2 emissions are different in the eight types of regions.The determinant power of the key determinant factors on O-AM, O-EU and NO-AS are: population growth rate > the growth rate of per capita GDP > energy intensity growth rate.The determinant power of the key determinant factors on NO-EE are: population growth rate > energy intensity growth rate > the growth rate of per capita GDP.The determinant power of the key determinant factors in NO-AF are: population growth rate > the growth rate of per capita GDP > energy intensity growth rate > natural gas carbon intensity growth rate.The determinant power of the key determinant factors in NO-AM are: population growth rate > the growth rate of per capita GDP.The determinant power of the key determinant factors in NO-ME are: population growth rate > coal carbon intensity growth rate > the growth rate of per capita GDP.