How to Effectively Control Energy Consumption Growth in China’s 29 Provinces: A Paradigm of Multi-Regional Analysis Based on EAALMDI Method

: Controlling energy consumption to reduce greenhouse gas emissions has become a global consensus in response to the challenge of climate change. Most studies have focused on energy consumption control in a single region; however, high-resolution analysis of energy consumption and personalized energy policy-making, for multiple regions with differentiated development, have become a complicated challenge. Using the logarithmic mean Divisia index I (LMDI) decomposition method based on energy allocation analysis (EAA), this paper aims to establish a standard paradigm for a high-resolution analysis of multi-regional energy consumption and provide suggestions for energy policy-making, taking 29 provinces of China as the sample. The process involved three steps: (1) determination of regional priorities of energy consumption control by EAA, (2) revealing regional disparity among the driving forces of energy consumption growth by LMDI, and (3) deriving policy implications by comparing the obtained results with existing policies. The results indicated that 29 provinces can be divided into four groups, with different priorities of energy consumption control according to the patterns of coal ﬂows. Most provinces have increasing levels of energy consumption, driven by increasing per capita GDP and improving living standards, while its growth is restrained by decreasing end-use energy intensity, improving energy supply efﬁciency, and optimization of industrial structures. However, some provinces are not following these trends to the same degree. This indicates that policy-makers must pay more attention to the different driving mechanisms of energy consumption growth among provinces.


Introduction
To reduce the possibility of serious, universal, and irreversible climate change and the effects it will have on human beings and ecosystems, a substantial and sustained reduction in greenhouse gas (GHG) emissions is necessary [1]. However, global energy consumption is still growing rapidly [2]. According to the latest forecast, global energy consumption will increase by more than half by 2060 [3]. Fossil fuel combustion is the main source of GHG emissions and still dominates global energy consumption, so controlling energy consumption growth is one of the means to effectively reduce GHG emissions globally and promote low-carbon energy transition [4,5]. This is especially significant for developing countries, whose energy consumption is rapidly increasing.
As each region or country is composed of several heterogeneous subregions, it is important to make differentiated policies for them. The formulation of differentiated (1) As enabled by the new paradigm, we have achieved the first application of EAA-LMDI to a multi-regional analysis. This paradigm provides a high-resolution way to formulate differentiated policies to control the energy consumption growth of multiple regions. (2) With high technical resolution, this paper reveals the similarities and differences in energy consumption, the growth of this consumption in China's provinces, and the reasons behind this growth. For the first time, energy allocation Sankey diagrams of 29 provinces in China are mapped and the LMDI decomposition results that consider the K PECR are presented. (3) The case study of China indicates a new policy-making model. When making regional energy consumption control targets, this model for determining the targets based on the characteristics of the main energy flow and the driving factors of energy consumption growth will be more effective than the commonly used current model, which directly determines the targets involving energy consumption and its growth.
The paper contents are as follows: a literature review of existing methods to analyze regional energy consumption is provided in Section 2; an introduction of the methodology and data can be found in Section 3; results and discussion are provided in Section 4; a brief discussion of limitations and uncertainties can be found in Section 5; and conclusions and recommendations for future research are available in Section 6.

Literature Review
Previous studies of multi-regional energy consumption were mainly based on MRIO. As for the studies of energy consumption in a single region or a few regions, econometric and decomposition analysis methods are common. In particular, the EAA-LMDI method can have high technical resolution, but was only used to study a single region.

Multi-Regional Input-Output Methods
MRIO is a common method for studies on multi-regional energy consumption. Through MRIO, scholars have calculated the energy consumption, carbon emissions, and other environmental impacts embodied in final consumption and cross-regional energy trade.
Compared with direct energy consumption, the research based on MRIO pays more attention to the indirect energy consumption embodied in final consumption and cross regional trade, as is shown in Table 1. Through the construction of cross-regional embodied energy and a carbon flow network, the multi-regional energy footprint and carbon leakage can be better explained. The analysis of embodied energy in final consumption can explain the change in regional energy consumption from the consumption side. However, these studies rarely consider the status and trends of energy systems in various regions, so their implications for making energy policy are limited. Table 1. Studies concerning multi-regional energy consumption based on multi-regional input-output (MRIO) methods.

No.
Author Research Object Main Conclusion 1 Lan et al. [23] Energy footprint of 186 countries When a country's per capita GDP is increasing, that country's total energy footprint becomes increasingly concentrated on imports or consumption. 2 Hong et al. [24] Embodied energy of the Chinese construction industry Energy flows are from resource-abundant areas in the central part of China to resource-deficient areas in the eastern coast. 3 Mi et al. [25] China's embodied carbon emissions In the "new normal", the strongest factors that offset CO 2 emissions have shifted from efficiency gains to structural upgrading. 4 White et al. [26] Water-energy-food nexus in East Asia China's current national export-oriented economic growth strategy in East Asia is not sustainable, as China is a net virtual exporter of nexus resources to Japan and South Korea. 5 Duan et al. [27] Embodied energy in final consumption of 121 countries The total energy consumption conditions across different countries throughout the world became more and more equal from 2000 to 2013. 6 Li et al. [28] Embodied carbon in international trade At a global level, the network density increases, indicating the widely expanding carbon leakage among economies. 7 Xia et al. [29] Embodied coal in the global economy 64.99% of the world's coal direct supply is ultimately embodied in international trade, thereby indicating the crucial roles of trade links in global coal utilization. 8 Guo et al. [30] Embodied energy of the service industry in six Asian cities Service industries consume 17.02-46.40% of the total embodied energy in the six cities. 9 Chen et al. [31] Embodied energy in international trade At a global level, a small-world scenario has been found, implying that economies are highly connected through embodied energy transfers. 10 Zhang et al. [32] China's provincial-level embodied energy Over a half of the national interregional transfers of embodied energy via domestic trade were induced by the three economic circles' final demand.

Econometric Methods
In the published literature using econometric methods, most focused on the driving factors of energy consumption, fitting the historical data of energy indicators and economic indicators. Only a few studies further predicted future energy consumption based on historical energy and economic data [11,33].
As illustrated by the studies listed in Table 2, econometric methods have an advantage in revealing the interactions among various indicators by testing their causal relationships. However, little consideration is given to the indicators related to the details of energy systems. For example, only a few studies included energy structures [11,12], and studies seldom considered energy supply efficiency. In addition, due to the free choice of driving factors, the factors considered cannot cover the whole process. Therefore, these studies have limited implications for making energy policy, which normally requires a high-resolution understanding of the energy system to determine key energy indicators, including the structure of primary energy, energy supply efficiency, the structure of energy end-use, and their driving factors. Urbanization intensifies the growth of energy consumption, but has a differential impact on the carbon emissions in different regions. 5 Wang et al. [9] GDP, carbon emission from fossil energy, and cement industry The carbon emissions of each province are increasing rapidly, but the law of growth operates differently in each province. 6 Liu et al. [35] Urbanization Urbanization helps to reduce the local energy consumption level, but its spillover effects increase the energy consumption level in neighboring areas.

Sheng et al. [36] Urbanization
There is a causal relationship between urbanization and the increase in energy consumption, but it cannot be proved that there is a causal relationship between urbanization and energy efficiency. 8 Si et al. [37] Energy policy Some energy policies have a positive effect on curbing energy consumption, but other energy policies have the opposite effect. 9 Nan and Gao [10] GDP China's current economic growth comes at the cost of huge energy consumption and environmental pollution.

Decomposition Analysis Method
In the published literature using the decomposition analysis methods, the methods can be divided into Index Decomposition Analysis [38] and Structural Decomposition Analysis (SDA) [39] methods. The fundamental difference between these two methods is a difference in the research objective [40]. IDA focuses on the relationship between commodity production and energy consumption, while SDA focuses on the relationship between final commodity consumption and energy consumption. Moreover, SDA needs to rely on the Input-Output model framework and IDA does not. SDA can better reflect the impact of the final demand structure on energy consumption during economic development. Meanwhile, IDA can better describe the process of physical energy supply and analyze it with technical drivers, which is more useful for energy policy-making. Therefore, IDA is more suitable for the purpose of studying energy policies in this paper.
In early studies, IDA decomposition results usually contained unexplained residuals. To solve this problem, Ang et al. proposed a new decomposition method named LMDI, which proved that the decomposition result does not contain unexplained residuals, and recommended LMDI as the preferred method for using IDA [41,42]. Since then, LMDI has been widely used for studying national and regional energy consumption and carbon emissions. Many scholars have used LMDI to analyze China's energy consumption at the national and regional levels. In Elsevier ScienceDirect, there are 103 papers on China's energy consumption and LMDI, most of which involve a LMDI decomposition analysis of energyrelated carbon emissions. In the literature on energy consumption, the focus is on energy consumption in the economic sector, high energy-consuming industries, and the building sector. Table 3 surveys the literature on energy consumption in the economic sector. Investment effect, energy intensity, economic structure, energy structure, and labor effect 1991-2011; China. Energy intensity is the main factor to reduce energy consumption. Investment and labor promote the growth of energy consumption.

2.
Sun et al. [44] GDP, economic structure, and energy intensity 2000-2010; Shenyang. The rapid growth of energy is mainly due to economic development, and energy intensity is the main restraining factor.

3.
Wang et al. [45] GDP, economic structure, and energy intensity 2003-2012; Tianjin. Industrial production scale is the main factor of energy consumption growth, and energy intensity has a restraining effect on growth. 4 Wang et al. [46] Population, GDP per capita, and energy intensity 1970-2012; China. The growth in per capita income and population promotes the growth of energy consumption. The energy intensity has a restraining effect, but it weakens after 2000. 5 Lima et al. [47] GDP, economic structure, and energy intensity 1990-2012; 6 countries. The decrease in energy intensity is not enough to offset the increase in energy consumption brought about by economic development. 6 Wang et al. [48] GDP, economic structure, and energy intensity 2006-2015; Hunan. The scale effect promotes the rapid growth of energy consumption; the efficiency effect can reduce energy consumption, while the structure has no significant inhibition on the growth of energy consumption. 7 Lin et al. [49] Energy structure, energy intensity, economic structure, spatial structure, investment efficiency, fixed asset investment per capita, and population proportion 2000-2016; 30 provinces. The expansion of fixed asset investment per capita is the main factor leading to the growth in energy consumption, energy intensity and investment efficiency are the main energy-saving factors, and the decoupling state of economy and energy consumption is mainly determined by secondary industry. 8 Feng et al. [50] GDP, economic structure, and energy intensity 2005-2013; 21 cities in Guangdong. The 21 cities are in different stages of economic development and energy utilization; Guangdong Province still has a great potential to improve its energy utilization efficiency. 9 Li et al. [51] Population, economic activity, and energy intensity 2011-2015; China. Population growth and per capita GDP growth are the leading factors that stimulate energy consumption, while the increase in consumption intensity is the main factor restricting consumption growth. 10 Wang et al. [52] Labor input, investment, energy intensity, energy structure, and technical level 1990-2015; China and India. The main influence of decoupling of energy and the economy is similar in China and India. Investment effect is the biggest driving force, while energy intensity is the biggest restraining factor, and China's energy-saving effect is better than India's. 11 Wang et al. [53] Population, GDP per capita, economic structure, energy intensity, and energy structure 2006-2015; 30 Chinese provinces. Output scale effect and production structure effect play a leading role in the growth of energy consumption in China, while the energy intensity effect can effectively inhibit the growth of energy consumption in China.
According to the above literature, the advantage of LMDI is that it can comprehensively analyze the impact of various driving factors on energy consumption growth, including economic factors and energy factors. Common driving factors include population, GDP per capita, economic structure, energy structure, and energy intensity [54]. Based on the abovementioned driving factors, the research on energy consumption (and energy-related carbon emissions) in specific sectors further considers factors such as the internal structure of the industry [55][56][57][58][59] (and carbon emission factors [60][61][62][63][64]). As these driving factors are usually related to policy-making, LMDI has better policy implications than econometric methods. Therefore, in general, LMDI is suitable for this study, but the technical resolution of energy systems needs to be further improved, and the impact of energy supply efficiency needs to be considered further.

EAA-LMDI Method
EAA-LMDI is a new advancement of LMDI to improve the technical resolution of energy systems; it originated from the development of the EAA method. Cullen and Allwood [65] are earlier scholars who proposed a high-resolution EAA method to compensate for the energy loss in the energy conversion process to each end-use department, and described the whole process from primary energy exploitation to end-use consumption by a Sankey diagram. Ma et al. [66] further applied this EAA method to China. After that, Chong et al. [13][14][15] proposed the EAA-LMDI method by combining EAA with the LMDI method, which could further improve the technical resolution of the LMDI method to include the impact of energy supply efficiency on energy consumption growth.
The main advantage of the EAA-LMDI method is that it can improve the understanding on the physical change in complex energy systems by mapping energy allocation Sankey diagrams in a time series, which provides whole maps of energy flows to illustrate the mechanism of how energy consumption is generated by the relationships among various sectors in terms of energy end-use, energy conversion, and energy supply. Secondly, by introducing K PECR , which is a reciprocal of the network efficiency from primary energy to end-use energy, the corresponding LMDI method can further disassemble the energy intensity in traditional LMDI research into end-use energy intensity and energy supply efficiency. Thus, the technical resolution of LMDI decomposition analysis is improved, and the impact of energy supply efficiency on energy consumption is reflected.
However, the main disadvantage of EAA-LMDI is the large data volume, as it is based on a time series of energy balance tables containing hundreds of energy data points. When it is applied to a multi-regional study, the volume and diversity of data will increase dramatically, which makes it difficult to pick out useful implications from a huge volume of data results. Therefore, if we want to reap the high-resolution advantages of EAA-LMDI to study multi-regional problems, we must first invent a standard paradigm that can effectively use the advantages of EAA-LMDI to help us enhance the understanding of the growth of multi-regional energy consumption, and easily develop policy recommendations based on it.

Methodology and Data
Although the main purpose of this paper was to build an analysis paradigm, this was based on the data treatment and analysis process of EAA-LMDI. Therefore, the basic principles and processes of EAA-LMDI are first introduced in brief, followed by a detailed introduction of the paradigm and a brief introduction of data sources.

Energy Allocation Analysis and Primary Energy Consumption Responsibility
The core idea of energy allocation analysis is to convert the end-use consumption of all secondary energy into the amount of primary energy needed to produce these secondary energies [65,66], which is named the primary energy consumption responsibility. The EAA-LMDI method calculates the primary energy consumption responsibility of the end-use sectors by introducing the primary energy consumption responsibility conversion factor (K PECR ), which is defined as the quantity of primary energy needed to generate a unit heat value of secondary energy j in the energy conversion process. K PECR can reflect the impact of secondary energy supply efficiency on regional energy consumption, and ensure that the primary energy consumption responsibility of the end-use sector is strictly equal to the amount of primary energy input into the regional energy system.
The calculation of K PECR and primary energy consumption responsibility can be performed as follows: (1) Obtain energy data from the regional energy balance table and reclassify according to the different energy types. (2) Construct a regional energy input-output table, and calculate the Leontief inverse matrix (L ij ) of the input-output table.
L 1,J -L 4,J represent the quantity of raw coal, crude oil, natural gas, and other energy, respectively, consumed in order to produce one unit of secondary energy j for final demand.
K PECR,j refers to the primary energy quantity conversion factor of the regional secondary energy j after the correction of the import impact, and K PECR,j,import refers to the primary energy quantity conversion factor of the external secondary energy j.
After calculating the primary energy consumption responsibility of each end-use sector, a Sankey diagram can be used to describe the whole process of energy utilization. Energy allocation Sankey diagrams of each province are divided into three stages: primary energy supply, energy conversion, and end-use sector. As there are many areas involved, the end-use sectors were simplified here. Figure 1 shows the structure of an energy allocation Sankey diagram.

EAA-LMDI Calculation
After introducing KPECR, the standard quantity of energy was transformed into primary energy consumption responsibility. The primary energy consumption responsibility of the economic sector and the residential sector can be expressed as follows: where i represents different economic or living sectors (primary, secondary and tertiary industries/town, rural) and j represents different energy types (raw coal, crude oil, natural gas, cleaned coal, coke, coal products, gas, gasoline, diesel, kerosene, oil products, heat, and electricity). KPECR,j represents the conversion efficiency of secondary energy j in the regional energy system. Accordingly, , in Equation (4) represents the efficiency from the secondary energy used by the end-use sector to the final economic output. That is to say, by introducing KPECR,j into the formula, the energy intensity in the traditional energy consumption decomposition formula can be divided into two factors, end-use energy inten-

EAA-LMDI Calculation
After introducing K PECR , the standard quantity of energy was transformed into primary energy consumption responsibility. The primary energy consumption responsibility of the economic sector and the residential sector can be expressed as follows: where i represents different economic or living sectors (primary, secondary and tertiary industries/town, rural) and j represents different energy types (raw coal, crude oil, natural gas, cleaned coal, coke, coal products, gas, gasoline, diesel, kerosene, oil products, heat, and electricity).
K PECR,j represents the conversion efficiency of secondary energy j in the regional energy system. Accordingly, in Equation (4) represents the efficiency from the secondary energy used by the end-use sector to the final economic output. That is to say, by introducing K PECR,j into the formula, the energy intensity in the traditional energy consumption decomposition formula can be divided into two factors, end-use energy intensity and energy supply efficiency, so as to improve the technical resolution of the LMDI decomposition analysis.
The decomposition results of LMDI can be expressed in the form of addition decomposition or multiplication decomposition. Addition decomposition is in the form of an absolute value, and multiplication decomposition is in the form of a relative value. In order to make a better multi-regional comparison and avoid the impact of regional energy consumption base, this study took LMDI multiplication decomposition as the main form of results display.

Standard Paradigm for Comparative Analysis of Multi-Regional Energy Consumption
The purpose of establishing the paradigm is to quickly process a large amount of data of EAA-LMDI, and to obtain knowledge about the status quo and historical trends of energy consumption in multiple regions so as to avoid the discussion of each region one by one. The whole paradigm consists of three parts, including a multi-regional energy allocation analysis, the LMDI decomposition analysis of multi-regional energy consumption, and policy implications and limitations of the paradigm, as shown in Figure 2 and explained in detail below. Flowchart of the multi-regional energy consumption analysis paradigm.
Step 1: Multi-regional EAA The purpose of this step was to complete energy allocation analysis in multiple regions, and compare the similarities and differences between the status quo of energy systems in different regions with the help of the energy allocation Sankey diagram. First, referring to the steps in Section 3.1.1, an energy input-output table was built for each region, and primary energy consumption responsibility was calculated by deriving KPECR. Second, based on the calculated primary energy consumption responsibilities, an energy allocation Sankey diagram for each region is drawn. Finally, by quickly scanning and comparing the main energy flows in the energy Sankey diagram, the main characteristics of energy systems in different regions were compared.
In this step, to quickly compare multi-regional energy systems, we classify the regions according to the main energy flows in the energy allocation Sankey diagrams, and then dis- Step 1: Multi-regional EAA The purpose of this step was to complete energy allocation analysis in multiple regions, and compare the similarities and differences between the status quo of energy systems in different regions with the help of the energy allocation Sankey diagram. First, referring to the steps in Section 3.1.1, an energy input-output table was built for each region, and primary energy consumption responsibility was calculated by deriving K PECR . Second, based on the calculated primary energy consumption responsibilities, an energy allocation Sankey diagram for each region is drawn. Finally, by quickly scanning and comparing the main energy flows in the energy Sankey diagram, the main characteristics of energy systems in different regions were compared.
In this step, to quickly compare multi-regional energy systems, we classify the regions according to the main energy flows in the energy allocation Sankey diagrams, and then discuss the main characteristics of energy systems in different types of regions. For example, in China, coal, as the main energy source, has a far greater impact than other energy sources. At the same time, coal consumption is also a key issue in China's environmental protection and low-carbon energy development. Consideration of the similar structure of end-use energy consumption in different regions (except Beijing and Shanghai), and the energy conversion stage is relatively complicated, so it is difficult to describe the characteristics of coal consumption in only one or two stages. Therefore, the production, import, and export of coal can be used as the criteria for multi-regional division in China's case. We named this criterion the "taxonomy of main energy flow characteristics". This criterion can be changed for various regions. For example, for some countries whose energy consumption is mainly based on oil, the main flow should be changed into oil flows.
Step 2: Multi-regional LMDI decomposition The purpose of this step was to complete a LMDI decomposition calculation of the energy consumption growth of the economic sector and residential sector in each region, and, through the comparison and discussion of the decomposition results, summarize the general impacts and differences of the multi-regional energy consumption growth mechanism. First, referring to the steps in Section 3.1.2, a LMDI decomposition calculation of energy consumption growth in each region was carried out. Secondly, combining the changes of various driving factors, through the comparison of the LMDI decomposition results of all regions, the impact of driving factors on energy consumption in multiple regions was summarized. Due to the relatively uniform changes in energy consumption in the residential sector, we omitted the discussion of each driver, respectively.
In order to avoid the impact of different energy consumption bases on the multiregional comparison, this paper mainly discusses the multiplication results of the LMDI decomposition. In the analysis of the economic sector, first of all, the historical trend of total energy consumption and the general impact of energy consumption change in multiple regions were compared. Subsequently, special cases that did not conform to the general impact under each driving factor are discussed one by one. In each discussion, the energy consumption growth rate and the growth are ranked from high to low.
Step 3: Policy implication analysis Combined with the research results and the existing multi-regional energy policies, the policy implications were further explored. First of all, we reviewed the formulation and actual implementation of the multi-regional energy policy objectives in the selected time period. Secondly, we compared the actual policy-making and completion with the calculation results to find out which factors cause the actual implementation to be inconsistent with the target. Finally, combined with the research in this paper and our understanding of multi-regional energy consumption, a new model of multi-regional energy policy objectives could be proposed.
In the case of China, we reviewed the formulation and actual implementation of China's multi regional energy policy objectives during the 12th Five-Year Plan period, and compared them with our research results. The actual results show that the current multi-regional energy policy-making in China does not fully consider regional differences and the actual energy situation. Therefore, we tried to propose a more scientific and feasible policy-making model based on the EAA-LMDI method.

Data Sources
This study selects the energy consumption of China's provinces in 2011-2016 as the research object. The energy consumption data of each province come from the energy balance table of each province in the China Energy Statistics Yearbook [67]. Population and economic data are from the National Bureau of Statistics [68]. It should be noted that, in addition to a lack of Tibetan statistical data, we found that there were problems with the power data of Guangdong Province in 2016 (the calculated power production efficiency of Guangdong Province in 2016 was 6.5%), so Guangdong was removed from the research.

Energy Allocation Sankey Diagram of Each Province in 2016
This study has drawn an energy allocation Sankey diagram for China's provinces in 2016, as shown in Figure 3. See Supplementary Materials for detailed energy allocation Sankey diagrams. The energy allocation Sankey diagram contains three stages: energy supply, energy conversion, and end-use. The first stage (energy supply) shows the quantity of all primary energy available for local use, including local production and imports. The second stage (energy conversion) shows the quantity of primary energy flowing into the energy conversion and then the derived quantity of energy carriers. The third stage (enduse) shows the quantity of all energy carriers flowing into the end-use sector, and also the distribution of energy consumption responsibilities among different end-use sectors. In all stages, all energy quantities were converted into corresponding primary energy consumption responsibilities, thereby ensuring a balanced energy quantity was maintained throughout the process.
While most regions in China have similar structures of end-use sectors, the energy conversion process is more complicated. In this study, we mainly classified and discuss multi-regional energy systems based on the characteristics of local energy production, import, and export, which is named the "taxonomy of main energy flow characteristics".
According to the "taxonomy of main energy flow characteristics", energy systems in 29 provinces can be divided into four categories. In the first category, the shape of coal flow is similar to "π"-that is, coal supply depends on imports from other regions, and is mainly used in power and heat generation. Most of the provinces' energy systems can be classified into this category, especially the eastern coastal regions. In the second category, the shape of the coal flow is similar to "⊥"-that is, it has a large amount of local coal production and export to other regions. The energy systems of several major coal-producing provinces conform to this pattern, including Shanxi, Shaanxi, and Inner Mongolia. In the third category, the shape of coal flow is similar to "Z"-that is to say, self-sufficiency is basically realized through the local production of coal, and the quantity of coal import and export is relatively small compared with local production and consumption. The energy systems of Anhui, Guizhou, Xinjiang, Ningxia, and Gansu conform to this pattern. In the fourth category, some provinces have only a small quantity of coal production and consumption, and other energy types occupy the main position. For example, Beijing has a larger proportion of natural gas and electricity, Shanghai and Hainan have a larger proportion of oil products, and Sichuan, Yunnan, and Qinghai have a larger proportion of hydropower.
According to this classification, considering that coal is the main source of China's carbon emissions, energy consumption in regions with a high coal dependence should be given priority. Eight regions in categories (2) and (3) are the main export areas of China's coal, and also consume a large amount of coal, so they need to be controlled first. Some regions in category (1), such as Hebei, Jilin, Shandong, and Henan, also have high coal dependence, and should also be key areas of energy consumption control. The regions in category (4) have a lower carbon energy structure, which can temporarily reduce the priority of energy consumption control in these areas.

LMDI Decomposition of Energy Consumption Growth in the Economic Sector
In this section, we calculated the LMDI decomposition results of energy consumption in China's 29 provinces' economic sectors from 2011 to 2016. Figures 4 and 5 show the additive decomposition results and multiplicative decomposition results of LMDI, respectively, with different colors indicating the quantitative influence of various driving factors to energy consumption growth. In Figure 4, all provinces are ranked from left to right following the total amount of energy consumption growth, as indicated by the red points. In Figure 5, all provinces are ranked from left to right following the growth rate of energy consumption.
carbon emissions, energy consumption in regions with a high coal dependence should be given priority. Eight regions in categories (2) and (3) are the main export areas of China's coal, and also consume a large amount of coal, so they need to be controlled first. Some regions in category (1), such as Hebei, Jilin, Shandong, and Henan, also have high coal dependence, and should also be key areas of energy consumption control. The regions in category (4) have a lower carbon energy structure, which can temporarily reduce the priority of energy consumption control in these areas.

LMDI Decomposition of Energy Consumption Growth in the Economic Sector
In this section, we calculated the LMDI decomposition results of energy consumption in China's 29 provinces' economic sectors from 2011 to 2016. Figures 4 and 5 show the additive decomposition results and multiplicative decomposition results of LMDI, respectively, with different colors indicating the quantitative influence of various driving factors to energy consumption growth. In Figure 4, all provinces are ranked from left to right following the total amount of energy consumption growth, as indicated by the red points. In Figure 5, all provinces are ranked from left to right following the growth rate of energy consumption.   Table 4 shows the general impact of each driving factor on energy consumption growth in the economic sector. The results revealed that, of the 29 regions, there are 20 that have increased their energy consumption in the economic sector. Of all the driving  Table 4 shows the general impact of each driving factor on energy consumption growth in the economic sector. The results revealed that, of the 29 regions, there are 20 that have increased their energy consumption in the economic sector. Of all the driving factors, GDP per capita has played the most significant role. Over the past five years, the growth rate of China's GDP per capita has remained above 6.2% every year. Rapid economic growth has led to a greater energy demand, which is the most important reason for the energy consumption growth of the economic sector in most provinces. Population and end-use energy structure are two other important driving factors of energy consumption growth. In the past five years, China's population growth rate has remained around 0.5%, which has led to increased energy demand. The adjustment of end-use energy structure mainly refers to the use of electricity and heat in the terminal to replace coal-based primary energy. As the electricity needs to bear more primary energy consumption responsibility, some provinces eventually consume more primary energy. The growth of the population and the economy is inevitable for China, which is still developing, but the adjustment of an end-use energy structure can lead to cleaner energy utilization and less environmental pollution. However, the increase in energy consumption caused by these three driving factors is inevitable. On the other hand, most provinces are curbing the growth of energy consumption by improving the economic structure, end-use energy intensity, and energy supply efficiency. China is reducing the proportion of the secondary industry in the economic structure and providing more support for the development of tertiary industry, such as the Internet and the financial industry. As the end-use energy intensity of the tertiary industry is significantly lower than that of secondary industry, this adjustment can effectively reduce energy consumption. At the same time, across both secondary and tertiary industries, China is trying to promote structural adjustment and the upgrading of industry, reduce the production of products with high energy consumption and low added value, and turn to the production of goods with higher economic value, so as to effectively improve end-use energy intensity. Each region is also continuously upgrading technology and equipment to improve the energy supply efficiency of the whole energy system, which also makes a contribution to reducing energy consumption.

The Special Cases of Each Driving Factor
In the vast majority of regions, China's population is still growing, resulting in the growth of energy consumption, as is shown in Figure 6. However, the population of Northeast China (Heilongjiang, Jilin, and Liaoning) has declined over the past five years, resulting in a decrease in energy consumption, which is not consistent with the overall trend of population growth in China. Due to the sustained low economic growth, the loss of population in Northeast China has become a trend that cannot be ignored, though it actually reduced energy consumption.
The vast majority of regions have maintained rapid economic growth, resulting in a significant increase in energy consumption as is shown in Figure 7. However, similar to the population, the growth rate of GDP per capita in Liaoning Province over these five years was negative, resulting in a reduction in energy consumption. The decline in the region's population and GDP per capita also reflects the relatively poor development of the region.
Energy supply efficiency Inhibition 28

The Special Cases of Each Driving Factor
In the vast majority of regions, China's population is still growing, resulting in the growth of energy consumption, as is shown in Figure 6. However, the population of Northeast China (Heilongjiang, Jilin, and Liaoning) has declined over the past five years, resulting in a decrease in energy consumption, which is not consistent with the overall trend of population growth in China. Due to the sustained low economic growth, the loss of population in Northeast China has become a trend that cannot be ignored, though it actually reduced energy consumption. The vast majority of regions have maintained rapid economic growth, resulting in a significant increase in energy consumption as is shown in Figure 7. However, similar to the population, the growth rate of GDP per capita in Liaoning Province over these five years was negative, resulting in a reduction in energy consumption. The decline in the region's population and GDP per capita also reflects the relatively poor development of the region. As is shown in Figure 8, the adjustment of economic structure has reduced energy consumption in most provinces. China is transforming its economic structure from the secondary industry to the tertiary industry. As the tertiary industry has lower energy intensity, this adjustment has greatly reduced energy consumption.
The only one that does not conform to the general impact of economic structure is Guizhou. On the one hand, the proportion of the secondary industry in Guizhou has increased, which may be related to industrial transfers; on the other hand, the energy intensity of the secondary industry in Guizhou has declined significantly, which means the increase in the proportion of secondary industry has not caused a major energy consumption increase. Therefore, the comprehensive impact of Guizhou's economic structure is a slight decline in total energy consumption. As is shown in Figure 8, the adjustment of economic structure has reduced energy consumption in most provinces. China is transforming its economic structure from the secondary industry to the tertiary industry. As the tertiary industry has lower energy intensity, this adjustment has greatly reduced energy consumption.
The only one that does not conform to the general impact of economic structure is Guizhou. On the one hand, the proportion of the secondary industry in Guizhou has increased, which may be related to industrial transfers; on the other hand, the energy intensity of the secondary industry in Guizhou has declined significantly, which means the increase in the proportion of secondary industry has not caused a major energy consumption increase. Therefore, the comprehensive impact of Guizhou's economic structure is a slight decline in total energy consumption.
The only one that does not conform to the general impact of economic structure is Guizhou. On the one hand, the proportion of the secondary industry in Guizhou has increased, which may be related to industrial transfers; on the other hand, the energy intensity of the secondary industry in Guizhou has declined significantly, which means the increase in the proportion of secondary industry has not caused a major energy consumption increase. Therefore, the comprehensive impact of Guizhou's economic structure is a slight decline in total energy consumption. With the continuous industrial restructuring, technological improvement and policy support, most regions have achieved remarkable results in reducing energy intensity, as is shown in Figure 9. The reduction in energy intensity plays a great role in restraining the growth of energy consumption. With the continuous industrial restructuring, technological improvement and policy support, most regions have achieved remarkable results in reducing energy intensity, as is shown in Figure 9. The reduction in energy intensity plays a great role in restraining the growth of energy consumption. However, the adjustment of end-use energy intensity in seven provinces has restrained the growth of energy consumption. Through analysis and discussion, we think these provinces can be divided into two situations. For the provinces such as Liaoning and Heilongjiang, the decrease in economic benefits leads to an increase in terminal energy intensity (the reduction in product added value per unit end-use energy consumption). In provinces such as Shanxi and Xinjiang, which are relatively rich in minerals and resources, there are some high-energy-consumption industries, resulting in the enhancement of the end-use energy intensity. In other words, the improvement of energy intensity in Shanxi and Xinjiang is mainly due to the influence of national industrial transfer policy. In most regions, replacing coal with electricity is a major trend due to the country's clean energy strategy, as is shown in Figure 10. Due to the energy loss in the process of power generation, electricity needs to bear more responsibility for primary energy consumption, which leads to increases in energy consumption. However, the adjustment of end-use energy structure in six provinces has restrained the growth of energy consumption. The decline in the proportion of electric consumption is the most important factor in the decrease in energy consumption, because electric power consumption needs to bear more primary energy consumption responsibility, followed by heat and oil products. However, the adjustment of end-use energy intensity in seven provinces has restrained the growth of energy consumption. Through analysis and discussion, we think these provinces can be divided into two situations. For the provinces such as Liaoning and Heilongjiang, the decrease in economic benefits leads to an increase in terminal energy intensity (the reduction in product added value per unit end-use energy consumption). In provinces such as Shanxi and Xinjiang, which are relatively rich in minerals and resources, there are some high-energy-consumption industries, resulting in the enhancement of the end-use energy intensity. In other words, the improvement of energy intensity in Shanxi and Xinjiang is mainly due to the influence of national industrial transfer policy.
In most regions, replacing coal with electricity is a major trend due to the country's clean energy strategy, as is shown in Figure 10. Due to the energy loss in the process of power generation, electricity needs to bear more responsibility for primary energy consumption, which leads to increases in energy consumption. However, the adjustment of end-use energy structure in six provinces has restrained the growth of energy consumption. The decline in the proportion of electric consumption is the most important factor in the decrease in energy consumption, because electric power consumption needs to bear more primary energy consumption responsibility, followed by heat and oil products. In most regions, replacing coal with electricity is a major trend due to the country's clean energy strategy, as is shown in Figure 10. Due to the energy loss in the process of power generation, electricity needs to bear more responsibility for primary energy consumption, which leads to increases in energy consumption. However, the adjustment of end-use energy structure in six provinces has restrained the growth of energy consumption. The decline in the proportion of electric consumption is the most important factor in the decrease in energy consumption, because electric power consumption needs to bear more primary energy consumption responsibility, followed by heat and oil products. As is shown in Figure 11, the improvement in energy supply efficiency mainly comes from electric power production, followed by heat. Shanghai is the only region where electric supply efficiency is reduced. For Shanghai, we do not think that the decrease in efficiency is due to a change of technology or equipment, because the technology and equipment in Shanghai remain at a high level. We believe that the decrease in electric conversion efficiency in Shanghai is due to the shortage of generation hours. The shortage of generation hours means that the high-efficiency capacity cannot be fully utilized and the power generation equipment will start and stop frequently, which leads to a reduction in efficiency. As is shown in Figure 11, the improvement in energy supply efficiency mainly comes from electric power production, followed by heat. Shanghai is the only region where electric supply efficiency is reduced. For Shanghai, we do not think that the decrease in efficiency is due to a change of technology or equipment, because the technology and equipment in Shanghai remain at a high level. We believe that the decrease in electric conversion efficiency in Shanghai is due to the shortage of generation hours. The shortage of generation hours means that the high-efficiency capacity cannot be fully utilized and the power generation equipment will start and stop frequently, which leads to a reduction in efficiency.

LMDI Decomposition of Energy Consumption Growth in the Residential Sector
Compared with the economic sector, the impact of energy consumption change in the residential sector is relatively simple. Figures 12 and 13 show the addition decompo-

LMDI Decomposition of Energy Consumption Growth in the Residential Sector
Compared with the economic sector, the impact of energy consumption change in the residential sector is relatively simple. Figures 12 and 13 show the addition decomposition and multiplication decomposition results, respectively, of LMDI. Figure 11. LMDI multiplication decomposition results of energy consumption growth caused by energy supply efficiency in the economic sector in 2011-2016.

LMDI Decomposition of Energy Consumption Growth in the Residential Sector
Compared with the economic sector, the impact of energy consumption change in the residential sector is relatively simple. Figures 12 and 13 show the addition decomposition and multiplication decomposition results, respectively, of LMDI.  In 26 out of 29 provinces, energy consumption in the residential sector has increased, as is shown in Table 5. After the LMDI addition decomposition, the three provinces with the highest energy consumption growth were found to be Hebei, Zhejiang, and Jiangsu. After the LMDI multiplication decomposition, the three provinces with the fastest growth in energy consumption were found to be Hainan, Chongqing, and Anhui. Among all the driving factors, the impact of per capita energy consumption is the most obvious, which also came from the improvement in living standards over the last five years. In addition, the increase in population, the promotion of urbanization, and the promotion of clean energy utilization also lad to a growth in energy consumption in the residential sector. Only the energy supply efficiency can restrain energy consumption growth. In 26 out of 29 provinces, energy consumption in the residential sector has increased, as is shown in Table 5. After the LMDI addition decomposition, the three provinces with the highest energy consumption growth were found to be Hebei, Zhejiang, and Jiangsu. After the LMDI multiplication decomposition, the three provinces with the fastest growth in energy consumption were found to be Hainan, Chongqing, and Anhui. Among all the driving factors, the impact of per capita energy consumption is the most obvious, which also came from the improvement in living standards over the last five years. In addition, the increase in population, the promotion of urbanization, and the promotion of clean energy utilization also lad to a growth in energy consumption in the residential sector. Only the energy supply efficiency can restrain energy consumption growth. Table 5. General impact of each driving factor on energy consumption growth in the residential sector.

Driving Factors
General Impact on Energy Consumption Growth Special cases that do not conform to the general impact include the following: (1) The decrease in population in Heilongjiang, Jilin, and Liaoning restrained the growth of energy consumption. (2) The decrease in urbanization rates in Zhejiang and Fujian led to lower energy consumption. The urbanization rate of these two provinces was mainly affected by the national household registration policy. Many migrant workers chose to move their household registration back to the countryside. (3) The per capita energy consumption of Heilongjiang, Jilin, Inner Mongolia, and Shanxi decreased. These areas reduced their per capita energy consumption by promoting central heating instead of bulk coal combustion. (4) The adjustment of end-use energy structure in Xinjiang, Hebei, and Shanghai resulted in a decrease in energy consumption. The share of electricity and heat consumption in these three provinces declined, while the share of fossil energy consumption (Xinjiang's oil products and natural gas, Hebei's raw coal, and Shanghai's oil products) increased. (5) The decrease in electric supply efficiency in Shanghai resulted in more energy consumption.

Formulation and Completion of Existing Policy Objectives
In order to corroborate the results of this study, this article investigates the achievement of energy policy targets in each region during the 12th Five-Year Plan period (2011)(2012)(2013)(2014)(2015). The 12th Five-Year Policy targets of each region were primarily concerned with GDP growth rate and energy intensity, as shown in Table 6.
Compared with the economic target, the energy intensity target did not fully take into account the actual situation in the region, but it only slightly fluctuated according to the overall national target. Energy intensity reduction targets in 26 provinces were in the range of 15-17%. At the same time, the differences in terms of the actual achievement of the energy intensity targets were significantly greater than the differences in the targets themselves.
According to the actual achievement of economic and energy intensity targets, the provinces can be divided into three categories for discussion: (1) Simultaneously achieved economic growth and energy intensity reduction targets-15 provinces: Tianjin, Jilin, Zhejiang, Anhui, Fujian, Shandong, Henan, Hubei, Hunan, Guangdong, Guangxi, Chongqing, Guizhou, Yunnan, and Qinghai. According to the decomposition results of LMDI, the decrease in energy intensity in these regions largely offset the increase in energy consumption brought about by economic growth, thus inhibiting the growth of energy consumption. (2) Achieved energy intensity reduction targets but failed to meet economic targets-15 provinces: Beijing, Hebei, Inner Mongolia, Liaoning, Heilongjiang, Shanghai, Jiangsu, Jiangxi, Hainan, Sichuan, Shaanxi, Gansu, and Ningxia (Ningxia's energy intensity did not include the coal chemical project in eastern Ningxia). Judging from the decomposition results of LMDI, these regions put more effort into the adjustment of industrial structure, which may have affected the economic growth to a certain extent. In addition, in our calculations, Liaoning and Heilongjiang did not meet their energy intensity reduction targets. After inspection, we found that, from 2015 to 2016, Liaoning's GDP decreased by 22.4%, and Heilongjiang's GDP remained basically unchanged, which may be related to changes in the caliber of data statistics. The underlying data issues make our calculations somewhat different from the energy intensity data published by the country. (3) Unmet energy intensity reduction targets: Shanxi and Xinjiang. Both regions are important energy bases for China. While undertaking energy production tasks, many high-energy-consuming industries have also grown up in these two regions, resulting in high energy intensity. Combined with the results of LMDI decomposition, Xinjiang and Shanxi's economic sector's energy consumption growth ranks in the top two in the country, and changes in energy intensity resulted in more energy consumption, which is also consistent with the completion of the target. Compared with the 12th Five-Year Plan, the 13th Five-Year Plan of each region simultaneously gives the control targets of economic growth, energy intensity reduction, and total energy consumption growth. Among them, the policy targets of most regions are not in strict accordance with the general formula of "energy consumption = GDP * energy intensity". On the one hand, this shows that each region has a more positive attitude and higher expectations for the "dual control" policy, but on the other hand, it also reflects the lack of a scientific basis to the formulation of regional policy targets. In addition, there are 28 provinces whose energy intensity reduction targets were between 14% and 17%. Therefore, it can be considered that the regional differences in energy intensity reduction targets are still not fully considered.

Suggestions on Scientific
Formulation of China's Multi-Regional "Dual Control" Policy Targets In this study, we also provide a new idea for how to scientifically formulate China's multiregional energy "dual control" targets. In this study, population, GDP per capita, economic structure, end-use energy intensity, end-use energy structure, and energy supply efficiency are considered the main driving factors. We can express energy consumption by using the following formula: energy consumption = population * GDP per capita * economic structure * end-use energy intensity * end-use energy structure * energy supply efficiency Therefore, each region can take corresponding measures according to its own situation around the above driving factors. For example, in more developed provinces, the number of migrants needs to be controlled. Several important energy bases, such as Xinjiang and Shanxi, need to improve the supply efficiency of the energy system. Provinces with more heavy industry can pay more attention to end-use energy saving and improve their end-use energy intensity. For provinces rich in renewable energy, optimizing the energy consumption structure and increasing the proportion of no fossil energy are important development directions for the energy system.
In addition to the driving factors, another consideration is the taxonomy of main energy flow characteristics. For regions with a higher proportion of coal, such as the eight areas belonging to categories (2) and (3) in Section 4.1, reducing energy consumption will be more helpful for reducing carbon emissions. Therefore, these areas should adopt more stringent energy-saving policies. For regions that have only a small quantity of coal production and consumption, such as the six areas belonging to category (4) in Section 4.1, relatively loose energy consumption control policies can be adopted. In these regions, the proportion of renewable energy is relatively high, and the increase in energy consumption will not significantly increase the carbon emissions. These regions can even encourage energy consumption to stimulate the consumption of renewable energy.
Future energy policy targets can also be set based on these driving factors. Compared with energy consumption, driving factors are better regulated and their future targets are easier to achieve. According to the development of the regional economy and energy, policy targets can be set for each driving factor separately. The "dual control" targets are determined based on the targets of these driving factors. The characteristics of main energy flow can be used as a reference for judging the priority of controlling energy consumption in different regions. If it is necessary to adjust the established targets, we can take appropriate methods and policies for each driving factor and adjust the changes expected in them so as to achieve the adjustment of the "dual control" targets. Compared with the existing policy-making methods, this idea can better reflect the actual situation in each region, and has a more scientific basis.

Limitations and Uncertainties
There are still some limitations and uncertainties in this study, including the limitations of the paradigm itself and the particularity of the case of China.
(1) This paradigm is based on the EAA-LMDI method. Although the method has the advantages of high resolution, some factors that need to be considered in policy are not included (resource availability is not considered, policy is not quantitative enough, and no fossil energy proportion factor is not introduced in LMDI). (2) The method to decide "taxonomy of main energy flow characteristics" is not accurate enough. In China, the problem of coal is prominent, but other countries may not have a similar situation. Therefore, the determination and analysis of this characteristic value may need further improvement. (3) The resolution of end-use consumption is not high enough. Compared with the continuous growth of China's energy consumption, energy consumption has remained stable in many countries and regions. Policymakers are more concerned with the distribution of energy consumption in the end-use sector, and need to analyze and discuss a higher resolution of end-use energy consumption. (4) A few provinces are not included. For example, Guangdong is not included in the research because of data errors, and Tibet lacks sufficient energy data. In the future, we can consider how to deal with an incorrect energy balance table and strengthen the research on Tibet. (5) The policy analysis section is relatively rough and needs further improvement.

Conclusions and Recommendations for Future Research
Based on the EAA-LMDI method, this study established a standard paradigm to analyze the growth mechanism of multi-regional energy consumption, and adopted it to analyze the status quo and historical trend of China's multi-regional energy system and energy consumption in 2011-2016. By drawing energy allocation Sankey diagrams of 29 provinces in 2016, we showed the status quo of the energy system in each province, and highlighted the main characteristics of the energy system in each province based on the "taxonomy of main energy flow characteristics". Based on the LMDI method, we decomposed the influence of driving factors on energy consumption of 29 provinces in 2011-2016, and summarized the energy consumption growth mechanism of each driving factor in different regions. We also discussed the existing "dual control" targets, and put forward a new way to formulate the "dual control" objectives of each region in the future.
The results show that this paradigm can be effectively used in China's multiregional energy consumption research; the main findings include: (1) According to the shape of coal flow, China's provinces can be divided into four categories. The shape of coal flow in most provinces is similar to "π", and provinces with coal flow shape similar to "⊥" and "Z" should be priority regions for energy consumption control. (2) The most important driving factor to promote the growth of energy consumption is GDP per capita, followed by population and end-use energy structure. The main factors that restrain the growth of energy consumption are the economic structure and end-use energy intensity, followed by energy supply efficiency. Compared with the improvement of energy supply efficiency, the decrease in end-use energy intensity has a more obvious effect on reducing energy consumption. The adjustment of the enduse energy structure (mainly by replacing coal with electricity) and the improvement of residents' living standards lead to more energy consumption. (3) The energy intensity of the secondary industry in some resource-rich areas and old industrial bases has increased, resulting in a significant increase in energy consumption. Fewer generation hours in some provinces affects the efficiency of electricity conversion, resulting in more energy consumption. (4) The formulation of China's regional energy policy targets during the 12th and 13th Five-Year-Plan periods did not fully consider the differences in energy consumption among regions. For the formulation of the multiregional "dual control" target, the development targets of the main driving factors can be determined first, and then the "dual control" policy targets are determined by the targets of the driving factors.
The core value of our work is to establish a framework that can analyze multi-regional energy consumption with high resolution, and provide research ideas for relevant researchers. Although we take China as a case study, this framework can also be applied to other regions in the world. This work can also provide reference for energy policy makers. We put forward the following policy recommendations for the formulation of China's multi-regional energy policy targets: (1) Strengthen the linkage between energy policy-making and other policy-making. Energy consumption is deeply affected by population, economic development, industrial structure and other factors. When formulating policies on population, economy and industry, the impact on energy consumption should be fully considered. Meanwhile, it is also necessary to verify whether the formulated energy policy targets can meet the development needs of the main driving factors. Furthermore, the adjustment of the policy targets is necessary when faced with conflict. (2) The taxonomy of main energy flow characteristics can be used as a reference to determine the priority of controlling regional energy consumption. The regions with a higher proportion of coal production and consumption, such as the eight areas belonging to categories (2) and (3) in Section 4.1, should be the priority regions of energy consumption control and adopt more stringent energy-saving policies.
The regions with a small quantity of coal production and consumption, such as the six areas belonging to category (4) in Section 4.1, can, relative to the regions with higher coal use, loosen the control of total energy consumption and encourage the development of non-fossil energy.
There are still some limitations and uncertainties in this paradigm that need to be further improved, as mentioned in Section 5. Moreover, the impact of energy trade between different regions on energy consumption is not fully considered. Therefore, in future research we will try to build a regional energy trade network, and find appropriate methods to analyze the impact of energy trade on regional energy consumption.

Conflicts of Interest:
The authors declare no conflict of interest.

LMDI
The logarithmic mean Divisia index I; EAA Energy allocation analysis; GHG Greenhouse gas; MRIO Multi-regional input-output methods; IDA Index decomposition analysis; EAA-LMDI The logarithmic mean Divisia index I decomposition method based on energy allocation analysis; K PECR The primary energy consumption responsibility conversion factor; SDA Structural Decomposition Analysis; Subscript i Economic subsector i involved during the LMDI decomposition of energy consumption growth in economic sector, including primary industry, secondary industry and tertiary industry/area i involved during the LMDI decomposition of residential energy consumption growth, including urban areas or rural areas; Subscript j Energy type j; K PECR,j The primary energy quantity conversion factor of the regional secondary energy j after the correction of the import impact; K PECR,j,import The primary energy quantity conversion factor of the external secondary energy j; E SQ Energy expressed in the form of standard quantity; E PECR Energy expressed in the form of primary energy consumption responsibility; E SQ,j,domestic Quantity of energy produced locally in standard quantity; E SQ,j,import Quantity of imported energy in standard quantity; P Population; GDP Gross domestic product; GDP i Value added within economic subsector i; E SQ,i Total energy consumption of economic subsector i expressed in standard quantity form during the LMDI decomposition of energy consumption growth in economic sector/total energy consumption of area i expressed in standard quantity form during the LMDI decomposition of residential energy consumption growth; E SQ,ij Energy type j consumption expressed in standard quantity form in the economic subsector i during the LMDI decomposition of energy consumption growth in economic sector/energy type j consumption expressed in standard quantity form of area i during the LMDI decomposition of residential energy consumption growth.