#### Abstract

Environmental issues caused by energy consumption have attracted increasing attention recently. Shanxi Province, a typical energy-dominated region in China, has long-term dependency on coal industry generating extensive economic growth, which is detrimental to green development. Distinguished from previous studies ignoring driving factors of energy consumption, this paper establishes a vector autoregression (VAR) model to dynamically identify the drivers of energy consumption based on STIRPAT model in Shanxi Province from 1990 to 2015. It can be obtained from the impulse response analysis that a positive shock in population, GDP, and urbanization level, respectively, positively affect energy consumption, and a positive change in technology negatively affects energy consumption in the long run. The variance decomposition results indicate that fluctuation in energy consumption explained by the innovation of the urbanization level accounts for 23.18%, which plays a prevailing role in increasing energy consumption. Meanwhile, the forecasting results of GM (1,1) model manifest that energy consumption in Shanxi Province generally has an increasing trend from 2016 to 2025. Consequently, Shanxi can achieve green development through optimizing energy structure, promoting the transformation of resource-based cities, and promoting low-carbon technological innovation. This paper can be available for other resource-based regions analogous to Shanxi.

#### 1. Introduction

As the energy is the significant resource and motive power, its cheap and stable supply is the safeguard of the economic and social development [1]. With China’s enormously economic growth, energy consumption shows a rapid increase trend [2–4]. However, the combustion of fossil energy is the main source of CO_{2} emissions, which makes China’s government face great challenges of global climate change and domestic environment protection to guarantee energy supply and energy-related environment security [5]. Meanwhile, China has now become the largest CO_{2} emitter in the world and joined “The Paris Agreement” to work together for coping with “greenhouse effect” globally after 2020 [6]. This implies that it is essential for the decision-makers to implement appropriate energy policy for China’s sustainable development.

As a spatially large developing country, China not only has prominent regional differences in climate and geography, but also has great diversities in energy endowments, economic development patterns, and levels of household consumption. It is critical to identify the main driving forces of energy consumption growth at the regional level, so that appropriate regional-specific energy conservation policies can be adopted.

Shanxi Province in central China is rich in mineral resources, among which coal, coalbed methane (CBM), bauxite, and iron are the most affluent. Specifically, coal has reserves of 2767.85×10^{8} tons, accounting for 20.1% of the country’s total reserves. And the reserves of CBM are 1825.16×10^{8} m^{3}, accounting for 88.2% of the total. Shanxi Province is a traditional coalmining region and a significant energy and industrial base in China. Since 1949, Shanxi has yielded 15.3 billion tons of coal which account for one-fourth of the national total [7], with the net amount of coal transferred being 10.4 billion tons accounting for three quarters of the national total. And coke production and the amount of the transfer account for 40% and 60% of the national total, respectively [8], making remarkable contribution towards China's economic development.

However, the long-term extensive economic growth model in Shanxi has engendered excessive dependence on the coal and a single industrial structure, inducing the issues of resource exhaustion and the formation of many subsidence and goaf areas. Meanwhile, the primary energy consumed in Shanxi was produced by the burning of coal, which accounts for a large share of greenhouse gas (GHG) emissions. Moreover, energy demand continues to increase and so does the air pollution, causing severe health problems, crop damage, and acid rain, all of which are taking a social and economic toll. Hence, it is beneficial and necessary to study the energy consumption in Shanxi Province.

Literature based on energy consumption in China can be generally classified as national and provincial angles. Many scholars have researched energy consumption in China from a nationwide perspective, and the initial study focused on the dual relationship between energy consumption and economic growth. Soytas and Sari [9] examined relationship between energy consumption and economic development in China and found that there was no causal relationship between total energy consumption and GDP. Talha and Cakar [10] researched the causal relationship between energy consumption and real GDP based on asymptotic distribution theory. As for the driving factors of energy consumption, Zhang and Song [11] used Logarithmic Mean Divisia Index (LMDI) method to study the factors affecting the changes in energy consumption in China and uncovering the economic growth effect was the primary factor. Li et al. [12] identified the factors driving China's energy consumption and demonstrated demand elasticity of impact factors changed regularly. Xie et al. [13] studied autocorrelation and influenced factors of energy consumption change in China based on the spatial autocorrelation analysis and the autoregression model. There are also some studies aimed at energy consumption predictions. Pao et al. [14] forecasted carbon emissions, energy consumption in China based grey model. Xiong et al. [15] proposed a novel GM (1,1) model to predict China's energy consumption and production. Wang et al. [16] estimated primary energy consumption based on the hybrid AMVO-SVM model to predict the energy consumption of China.

Literature concerning the energy consumption at national level gave macro conclusions which were not applicable in all regions with a noticeable diversity across every region throughout China. Thus, scholars continued to study more concretely several representative provinces and cities. Most of the studies mainly analyzed the energy consumption of economically developed regions from a static point of view.

Wang and Feng [17] discussed energy resource, production, and consumption in Jiangsu province and put forward a long-term policy to ensure the sustainable development. Cai and Jiang [18] explored differences in energy consumption between rural and urban households and designed relative conservation implications in Shaanxi Province. Ren et al. [19] studied energy consumption and carbon emissions during the urbanization of Shandong Province and demonstrated urbanization had distinct positive impacts on energy consumption and carbon emissions. Chong et al. [20] analyzed the driving factors of energy consumption growth in Guangdong Province based on the LMDI decomposition model. Feng et al. [21] used system dynamics approach to probe into the influence of driving factors selected on the energy consumption and CO_{2} emission trends in Beijing. There are also a few studies that have comprehensively analyzed the energy consumption of Chinese provinces. Zhang and Lahr [22] applied the structural decomposition method to uncover the regional disparities in energy consumption and proposed regional-specific policies. Wang et al. [23] investigated the relationship between urbanization, energy consumption, and CO_{2} emissions based on the data from 30 Chinese provinces by a panel data model.

In terms of Shanxi Province, Wang [24] et al. conducted an empirical analysis of the relationship between economic growth and energy supply and demand in Shanxi Province. The results showed that energy consumption and energy supply together constituted a unidirectional driving engine for economic growth in Shanxi Province. Xu et al. [25] used the SVAR model to analyze the dynamic relationship between urbanization and energy consumption and concluded that urbanization in the medium and long term contributed more to the variance in energy consumption. Li and Yang [26] analyzed the energy consumption system and related subsystems in Shanxi Province by constructing a system dynamics model and provided corresponding policy recommendations to ensure sustainable development. Zhai [27] used the ARIMA model to predict the time series of energy consumption in Shanxi Province.

Shanxi is a national resource-based economy restructuring comprehensive reform pilot area in China. Reviewing existing literature, however, research conducted on driving energy consumption in Shanxi has been very scanty. Moreover, previous studies have mainly emphasized on the relationship between energy consumption and a few concerning variables but ignored its driving factors. Meanwhile, the ARIMA model frequently adopted does not have high prediction accuracy when the sample size is small.

In summary, to make a contribution to filling the gaps in existing studies, we conduct the paper from the following aspects. First, we use the Stochastic Impacts by Regression on Population, Affluence and Technology (STIRPAT) model to explore determinants of energy consumption of Shanxi Province at the regional level. Second, considering the dynamical research of energy systems lacked in Shanxi, a vector autoregression (VAR) model based STIRPAT model is established to quantitatively analyze the influencing factors of energy consumption growth, which is not subject to traditional economic theory and to avoid the estimated error caused by multicollinearity [26]. Third, with relatively little data, the grey forecasting model provides approximately realistic forecast values of energy consumption for corresponding policy establishment [28]. Finally, conclusions and policy recommendations to boost green sustainable development in Shanxi are presented based on results received, which is essential for controlling energy consumption and deserves further study for other resource-based regions analogous to Shanxi.

Therefore, this paper sheds light on the current studies in the following points. (1) It is the first try to combine the VAR model and the extended STIRPAT model to analyze the energy consumption issue in Shanxi Province, a vital energy-dominated region in China. (2) The influenced factors selected in the extended STIRPAT model are distinguished from previous research, to explore a different study perspective concerning energy consumption in Shanxi. (3) It is the first time that we use the GM (1,1) method, with relatively little data, to forecast energy consumption in Shanxi from 2016 to 2025.

#### 2. Method

##### 2.1. VAR Model

The VAR model, one of the easiest models to analyze and predict multiple related economic indicators, is proposed by Sims [29] based on statistic characteristic of data comprehensively making use of multiple variables information. Further, a VAR model can better cover enough information to simulate the system than that using merely one single time series, providing more realistic predicted values in forecasting when making the near future prediction.

The multivariable VAR (P) model can be obtained bywhere is a vector of endogenous variables, coefficient matrices, the one-step ahead prediction error with variance-covariance matrix , the intercept, and the lag length. As for the VAR model, the single coefficient can only reflect a local relationship and fails to capture the overall impacts of one variable on the other, which is meaningless to study the coefficients in this model. Thus, based on the VAR model established, the main purpose on this occasion is to analyze the dynamic relations between variables by impulse response analysis and variance decomposition.

Impulse response analysis can comprehensively reflect not only the dynamic relationship among various variables but also influence of the change in an endogenous variable on others in the system. Specifically, it demonstrates impacts on the current or future value when a random error term is shocked from inner or outer system.

Equation (1) corresponds with the stability condition if roots of the determinant are outside the unit circle, so that in this case the VAR model can be shown as an infinite vector VMA process, and (1) can be rewritten bywhere is the lag operator. For two conditioning random variables, if is given a unit of impulse in the base period: , , , , , the impulse response function of caused by can be expressed as

Variance decomposition will decompose the variance of one variable in the system into each disturbance, to analyze the innovation of the endogenous variables from the impact of each structure, so that the significance of various structural shocks can be further explored. In this case, (2) can be rewritten as

The content of each bracket denotes the sum of influence of the disturbance from the infinite past to the present, which can be evaluated by the variance. Supposing is no serial correlation, variances of can be decomposed into k types of unrelated effects. Hence, the relative variance contribution rate is defined to determine how much the variance of disturbance contributes to the variance of .

If the VAR model satisfies the stability condition, the approximate relative variance contribution rate is shown as follows:where refers to standard deviation of . The greater the , the larger the impact of the variable on the variable.

##### 2.2. GM (1,1) Model

The grey system was firstly proposed by Deng [30]. The grey forecasting model can be applied in cases with relatively fewer statistics, which is suitable for predicting the competitive environment where decision-makers can refer to mere limited historical data. The GM (1,1) procedure can be summarized as follows.

*Step 1. *The nonnegative raw time-sequence data is defined as

*Step 2. *Take accumulated generating operation (AGO) on into a new sequence :where

*Step 3. *The basic form of GM (1, 1) model is obtained by the following equation:where is a developing coefficient and represents the grey input. Parameter is represented bywhere and

*Step 4. *According to the estimated coefficients and , the response equation and the prediction value can be obtained byTo demonstrate the reliability of the GM (1, 1) model, the posterior variance test is utilized to discuss the results of fitting and prediction. Criteria of the posterior variance test contain the indicators, posterior variance ratio , and small error probability . The posterior variance ratio is defined aswhere and . represents the residual between the and

The small error probability is defined asThe range of indicator values and is such that , , where and indicate a highly accurate forecast, and indicate a reasonable forecast, and and indicate an unreasonable range

##### 2.3. STIRPAT Model

To study the driving factors of the environment, Ehrlich [31] first proposed IPAT model in 1971, i.e.,where denotes the environment pressure, represents population, indicates affluence, and is technology. The IPAT equation hypothesized that environmental impact changes proportionally or monotonically, which is not accordance with the reality condition.

To overcome this limitation, Rosa and Dietz [32] developed a model named STIRPAT. The specific model can be described aswhere , , , and are the same variables as those in (14); , , , and are coefficients representing stochastic impacts; denotes error term. The natural logarithm form is shown in

The STIRPAT model is entitled to be extended by incorporating extra factors for researching more possible causes for environmental pressure [33]. Therefore, how to select influenced factors is the key point of the whole study. In order to investigate the impacts of the driving forces of the energy consumption in Shanxi Province, considering the specific development situation and learning from relevant research experience, we rewrote (16) as follows:

Energy consumption can be regarded as a kind of environmental impact because it consumes natural resources, denoted by EC. denotes population size, is measured by per capita GDP, represents the total factor productivity (TFP), represents urbanization level, and , , , and are estimation parameters. So we rewrote (17) as follows:

##### 2.4. Data Illustration and Collection

The paper adopted the data interval 1990–2015; all data in this paper covering the period of 1990-2015 are derived from the Shanxi Statistical Yearbook [34]. More specifically, technology progress is measured by TFP, because this indicator accurately reflects the level of technology. This paper selects DEA-Malmquist productivity index method based on dynamic nonparametric frontier production surface. This method is based on balancing panel data containing multiple decision units, applying input-output data, identifying the best business decision-making unit for technical efficiency, and constructing the production frontier based on it. The technical efficiency level index is obtained by comparing the actual input or output value of the production process with the optimal value. And this method allows technology inefficiency to exist and is therefore widely used for total factor productivity estimation. The DEA 2.1 software is used in this paper to measure TFP of Shanxi Province in each year. The key of the DEA-Malmquist productivity index method is the selection of input and output variables. The inputs of this article are selected from the two indicators of physical capital and labor. The physical capital is represented by the capital stock, the labor index is selected by the number of employed people, and the output index is expressed by GDP. This paper draws on the research of Shan Haojie [35] to select capital stock data estimation methods and indicators. In summary, the definitions of all variables are shown in Table 1.

#### 3. Result

##### 3.1. The Changing Tendency of Energy Consumption

Figure 1 shows the trend of Shanxi's energy consumption during 1990 and 2015. We can obviously see that the total energy consumption changes by various stages. The first stage is from 1990 to 1999, during which energy consumption experienced a gentle growth with an average annual growth rate of 2.43%. The second stage is from 2000 to 2013, during which energy consumption of Shanxi increased rapidly with an annual growth rate of 8.14%. The third stage is from 2013 to 2015, when energy consumption of Shanxi suddenly decreased. This may be attributed to the reason that the downturn of coal markets and the rapid drop of energy prices since 2013 had caused a severely negative shock to the economy of Shanxi Province.

As is shown in Figure 2, from sectional perspective, the coal, electricity, and coke remained to be leading parts in energy structure in the past, accounting for approximately 90% of total energy consumption from 1994 to 2004. While this percentage declined with small fluctuations from 2005 to 2015, it has still remained over 80%. As for the oil, its consumption has always accounted for a comparatively tiny proportion of the whole; only 7.2% of total energy consumption in 2015. Last but not least, although the percentage of natural gas showed a general trend of growth, the increasing rate was very slow. This indicates that relevant departments should vigorously develop clean energy and improve energy structure.

The original VAR model requires that each variable should be stationary. For nonstationary variables, the differential processing needs to be performed to obtain a stationary time series to establish a VAR model, which will nevertheless lose the corresponding level information. With the help of cointegration theory, VAR models can be directly established as long as there exist cointegration relationships between variables [36].

It should be the first to conduct the stationarity test to determine the stationarity of variables and its order of integration, with the test methods including the augmented Dickey-Fuller (ADF) test and the nonparametric Phillips-Perron (PP) test [37]. In this paper, the ADF test is employed.

From Table 2, as for the logarithm form of EC, P, GDP, T, and URB, the correlative ADF statistics are greater than the critical value at the 10% significance level, indicating that the logarithm forms of these variables are all nonstationary. Tests performed on the first difference intensely demonstrate that the first difference of each time series is stationary at the 5% significance level; that is, these series are integrated of order first.

As for the cointegration test, variables are required to be integrated. The stationarity test result indicates that all variables are integrated of order first, which is reasonable to perform the cointegration test. There are two mainly used cointegration test methods: the Engle-Granger (E-G) test based on the regression residual applied to the bivariate and the Johansen-Juselius (JJ) test based on the regression coefficient for the multivariate. In our model, the JJ test method is applied.

Before the cointegration test, we should determine the optimal lag length, whose specific operation is expressed as follows: beginning from a larger lag order, the optimal lag length can be determined by comprehensive judgment through the LR value, FPE value, AIC value, SC value, and HQ value. From Table 3, the test results demonstrate that the optimum lag length is one.

According to characteristics of the statistics, we select the linear tendency of the sequence, but the equation only has the test form of intercept. The statistical results of the trace test shown in Table 4 illustrate that null assumption of no cointegration relationships is rejected at the 5% significance level. As a result, there exists a cointegration relationship between variables. It can also obtain the same conclusion from the maximum eigenvalue test results in Table 5.

##### 3.2. The Impulse Response Function Analysis

The stability of VAR model requires to be tested before the impulse response function analysis which is only suitable for the stable VAR model [38]. The model will be stable when estimated absolute values of all inverse characteristic roots of the characteristic equation are less than one; that is, inverse roots of AR characteristic polynomial are inside the unit circle.

From Figure 3, absolute values of all unit roots are inside the unit circle. Therefore, the VAR model meets the stability conditions so that impulse response function analysis can be explored.

From Figure 4, a shock to energy consumption positively effects itself, reaching the minimum in the fifth period, showing a U-shaped curve; a one standard deviation positive change to population positively affects energy consumption, climaxing in the second period, and this effect is stabilized after the seventh period; GDP responds positively to a one standard deviation positive change in energy consumption, and this response is relatively long-lived, reaching the maximum in the second period; a positive shock in technology negatively affects energy consumption after the third period, while this effect is very little previously; a positive change in the urbanization level exerts a positive effect on energy consumption, climaxing in the third period, showing an inverted U-shaped curve.

##### 3.3. Variance Decomposition

Table 6 indicates that the shock to energy consumption account is 100% and 44.47% fluctuation in itself in the first quarter horizon and tenth quarter horizon, respectively, and it has been a downward trend from 1 to 10 quarters. Fluctuation in energy consumption from 0 to 9.91% can be explained by the innovation to population, showing a rising trend from the first period to the tenth period and can also be attributed to innovations in GDP ranging from 0 to 17.97%, the innovation of which has been increasing over the ten quarters. The innovation to technology accounts for 0 and 4.47% fluctuation in energy consumption in the first horizon and tenth quarter horizon, respectively. The urbanization level can be attributed to innovations in energy consumption from 0 to 23.18%.

##### 3.4. Forecasting of GM(1,1) Model

According to the Shanxi-based energy consumption data, the prediction formula of GM(1,1) model for the energy consumption can be expressed as

The results of predicted values and actual ones from 1990 to 2015 are shown in Table 7. It can be calculated that the mean percentage error of GM(1,1) model is 7.264 %, the mean prediction accuracy is 92.736%, the posterior variance ratio is 0.132, and the P value of posterior variance test is 1.000. In accordance with the optimum range of indicator values C and P, as aforementioned, the model produced approximately perfect estimates.

The energy consumption from 2016 to 2025 can be predicted by (19). As can be seen from Table 8, energy consumption of Shanxi Province is continually increasing in the next few years. By 2025, the energy consumption will increase to 31103.82 ten thousand tons of standard coal.

#### 4. Discussion

The main drivers of energy consumption are based on the extended STIRPAT model. The results indicate the natural logarithm form of energy consumption, population, GDP, technology, and the urbanization level are all integrated of order first and subject to long-term equilibrium relationship. The VAR model can be built accordingly to probe into the dynamic connection between them.

Based on the impulse response analysis, a one standard deviation positive change to population positively affects energy consumption, showing that the population is a critical factor influencing the growth of energy consumption in Shanxi Province. The increase in population has augmented the demand for products and services, which has led to the growth of energy consumption demand.

GDP responds positively to a one standard deviation positive change in energy consumption, indicating that the economic growth has had a continuous effect on the growth of energy consumption in Shanxi Province. This is due to the reason that as the abundant coal resources in Shanxi Province, the local resource attributes determine that its long-term economic development will inevitably lead to the large-scale development and utilization of energy.

A positive shock in technology negatively affects energy consumption after the third period, which can be attributed that technological progress cannot curb energy consumption in the short term. As it will take a long time that the new technology diffusion after research and development, technology can only play a role in reducing energy consumption after a period of accumulation. This result shows that, in the long term, technological progress is an effective means of reducing energy consumption.

A positive change in the urbanization level exerts a positive effect on energy consumption. The reason for this trend is that cities are the main body of energy consumption in China. Along with the advancement of urbanization, the demand for energy in urban buildings, transportation, and residential buildings is also increasing. As for Shanxi, the promotion of urbanization process mainly depends on the development of industry, which has intensified energy consumption accordingly. The response curve shows an inverted “U” trend, indicating that there may exist an environmental Kuznets curve (EKC) between Shanxi’s coal consumption and urbanization rate.

From the variance decomposition results, fluctuation in energy consumption from 0 to 9.91% can be explained by the innovation to population. Due to the continuous improvement of living standards, the demand for energy production and consumption will also increase as a result, and the impact of population effects on energy consumption will gradually enhance. The innovation to GDP for fluctuation in energy consumption continues to be increasing, which accounts for 17.97% in the long run, meaning that in medium and long-term the economy of Shanxi Province is in its extensive growing stage, and hence the resource-consuming economic growth is accompanied by an increase in energy consumption. The innovation to technology accounts for 4.47% fluctuation in the long run, indicating that the technology level contributes a small part to energy consumption, but raising the level of coal utilization and energy efficiency through technological innovation is still a key factor in reducing energy consumption. A positive change in the urbanization level exerts a positive effect on energy consumption. This tendency mainly arises from the fact that the increase of urbanization rate has prominent contributions to medium and long-term energy consumption, because the burden of urban construction caused by agglomeration of rural people to cities has increased the use of fuel and cement and fossil energy consumption in daily life.

GM(1,1) model can be established to perfectly predict the energy consumption of Shanxi Province. According to the forecasting model, it is expected that the energy consumption of 2025 of Shanxi will reach 31103.82 ten thousand tons of standard coal. As a result, the government and relevant departments should effectively adopt reasonable policies to increase energy efficiency to avoid the crisis of energy exhaustion and realize the sustainable development.

#### 5. Conclusions and Policy Implications

This paper explores dynamic relationship between the driving forces determined by the extended STIRPAT model and energy consumption in Shanxi Province through VAR model. We found clean energy represented by the natural gas accounted for a comparatively tiny proportion of the whole. It can be obtained from the impulse response analysis that a positive shock in population, GDP, and the urbanization level, respectively, exerts a positive effect on energy consumption, with all of those responses being relatively long-lived. And a positive change in technology negatively affects energy consumption in the long run. The variance decomposition results indicate that fluctuation in energy consumption can be respectively explained by the innovation to population, GDP, technology, and the urbanization level accounting for 9.91%, 17.97%, 4.47%, and 23.18% in the long run. A large part of Shanxi's energy consumption growth can be explained by urbanization. Meanwhile, the forecast results of GM(1,1) model manifest that energy consumption in Shanxi Province generally has an increasing trend year by year.

These conclusions are illuminating to attract attention from local policy-makers. Specific implications are summarized as follows in aspects of energy structure, urbanization, and technology in order to ameliorate traditional coal-dominating situation, earnestly implement the new Five Development Concept, and accomplish in-depth reform of resource-based economy and green sustainable development in Shanxi Province.

*(**1) Reform the Energy Structure and Propel Diversified Energy Consumption.* Energy consumption of Shanxi is mainly dominated by the coal of primary energy source. Owing to its high-pollution and nonrenewable properties, the overuse of coal should be slashed dramatically and green renewable energy deserves to be boosted actively. Thus, relevant departments ought to expedite the implementation of “coal to gas” projects for civil and industrial use and other key alternative areas.

As a province with abundant CBM resource, Shanxi has a well-distributed CBM field with good conditions, so that it would be helpful to increase the development and utilization of CBM resources. Meanwhile, with rich resources of solar and wind energy, it is necessary and beneficial to steadily and smoothly push forward the construction of photovoltaic bases and northern Shanxi wind power bases, which will be of great help for the amelioration of conventional coal-dominating situation. The progress of other new energy sources should be strengthened persistently, with the development relationship between traditional and new energy sources being coordinated.

*(**2) Promote the Transformation of Resource-Based Cities and Clarify the Orientation of Urban Development Planning.* The economic growth mode of Shanxi Province is mainly driven by resources consumption. The promotion of urbanization mainly relies on the development of industry, and the urbanization process has obviously aggravated energy consumption. Therefore, in the process of urbanization, according to the local energy supply and demand situation, the urbanization rate should be appropriately slowed down, and new urbanization should be promoted to guarantee green energy consumption, and the transformation of resource-based cities should be promoted to target deep processing of resources, lengthen the coal industrial chain, and vigorously promote the novel coal chemical industry. At the same time, it is helpful to develop nonresource industries, establish high-quality cereals producing areas, vigorously develop the tourism industry, and promote the integration of culture and tourism.

Shanxi should correctly understand the development stage and characteristics of the city, define the future direction and path of the city, and give full play to cities’ agglomeration and radiation potential. It is of great significance to strengthen the cooperation in division of labor among big cities and its surrounding small and medium-sized ones, promote the integration process of the urban agglomerations and the integration of tourism resources, and build Shanxi Merchant Culture tourism boutique line. The local government should continue to expand the implementation scope of the “Hundred Town Construction Projects” and explore more potential small towns to give them a certain policy inclination and form a group of specialized towns with special characteristics.

*(**3) Increase Environmental Investment and Enhance Low-Carbon Technological Innovation.* The environmental technology should be given top priority, playing its critical role in technological progress in the energy conservation and development of low-carbon economy. Financial subsidies for “coal-to-electricity” of rural power grid reform and optimization project should be augmented to actualize a virtuous circle of energy production system.

Priorities for technology aspect should be vigorously driving the comprehensive utilization of coal to stimulate key areas such as coking and coal chemical industry to implement technology transformation for cleaner production. In addition, it is vital to continuously enhance the capability of independent innovation of the coal industry, extending the value chain of coal resources and increasing added value of coal industry. Furthermore, the implementation of ultra-low emissions and energy-saving transformation for coal-fired units and clean coal power generation technology will improve the overall quality of coal industry in Shanxi Province.

#### Data Availability

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

#### Conflicts of Interest

The authors confirm that the mentioned received funding in Acknowledgments did not lead to any conflicts of interest regarding the publication of this manuscript. And the authors declare there are no conflicts of interest in the manuscript.

#### Acknowledgments

This work was supported by the National Natural Science Foundation of China (no. 71471061) and the National Social Science Foundation of China (no. 15BGL145).