Next Article in Journal
A Sustainable Home Energy Prosumer-Chain Methodology with Energy Tags over the Blockchain
Previous Article in Journal
Model Prediction of Secondary Soil Salinization in the Keriya Oasis, Northwest China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Impacts of GDP, Fossil Fuel Energy Consumption, Energy Consumption Intensity, and Economic Structure on SO2 Emissions: A Multi-Variate Panel Data Model Analysis on Selected Chinese Provinces

1
School of Economics and Management, North China Electric Power University, Beijing 102206, China
2
Beijing Key Laboratory of New Energy and Low-Carbon Development, North China Electric Power University, Beijing 102206, China
*
Author to whom correspondence should be addressed.
Sustainability 2018, 10(3), 657; https://doi.org/10.3390/su10030657
Submission received: 22 December 2017 / Revised: 16 February 2018 / Accepted: 21 February 2018 / Published: 1 March 2018

Abstract

:
Atmospheric pollution gradually become a focus of concern all over the world owing to its detrimental influence on human health as well as long range impact on global ecosystem. This paper investigated the relationship among SO2 emissions, GDP, fossil fuel energy consumption, energy consumption intensity, and economic structure of five provinces in China with the highest SO2 emissions spanning from 2002–2015 based on panel data model. Through comparatively analyzing the coefficients in the established panel data model for Hebei, Henan, Inner Mongolia, Shandong, and Shanxi, we can obtain that: (1) fossil fuel energy consumption made the most devotion to SO2 discharge compared with GDP, energy consumption intensity, and economic structure. And the more the fossil fuel energy consumption, the more the devotion made by it to SO2 discharge. (2) GDP devoted less to SO2 emissions than fossil fuel energy consumption, and the larger the scale of the economy, the greater the contribution made by it to SO2 emissions. (3) The higher the proportion of the secondary industry added value accounted in GDP, the more the devotion made by the economic structure and energy consumption intensity to SO2 emissions. Through analyzing the Granger causality examination results, it can be concluded that: (1) there existed a bi-directional causal relationship between fossil fuel energy consumption and SO2 emissions among five selected provinces. (2) There existed uni-directional causal nexus running from GDP to SO2 emissions, from energy consumption intensity to SO2 emissions, and from economic structure to SO2 emissions among five chosen provinces. Based on the empirical analysis, several policy implications were proposed to provide references for policy makers, which were (1) Giving full play to the guiding role of price signals, and improving the price policy for desulfurization. (2) Formulating a new comprehensive evaluation system to measure the regional development level considering economic development and environmental impacts. (3) Exploring renewable and sustainable energy sources to substitute for fossil fuel energy according to regional resources endowment. (4) Developing high value added and low pollution emissions industries and reducing the proportion of secondary industry.

Graphical Abstract

1. Introduction

As the second largest economy in the world, China is facing serious environmental issues owing to rapid development of economy. From the start of 2003, an uncommon heavy fog and haze weather swept over central and eastern region in China, which drew widespread attention from public [1,2]. Since then, the environment has gradually deteriorated. In 2014, the air quality of nearly 90% China’s cities have not met the standards [3]. Additionally, in accordance with the 2016 Environmental Performance Index, the air quality of China ranked 109th in 178 countries [4]. Such serious air pollution problems in China not only threaten the health of local people, but also pose great damage to the global air quality [5]. Thus, the reasons caused environmental problems, and the measures taken to control the pollution in China gain widespread concerns around the world [6]. Even though a series of measures have been implemented by Chinese government to enhance air quality, such as establishing the national air quality monitoring network, closing down small energy-intensive enterprises with heavy pollution, eliminating backward production capacity, and developing renewable energy, the effects have not been remarkable.
Due to China’s heavy dependence on fossil energy, including coal energy, and petroleum energy, as well as secondary industry, sulfur dioxide (SO2) is one of the significant air pollutants which emissions amount came to the top at 25.89 million tons in 2006 [7]. Although Chinese government set series of goals to cut down SO2 emissions, the volume of SO2 emissions is still a critical burden for the environment. Therefore, recent years, some literatures researched on the driving forces of SO2 emissions and environmental Kuznets curve for SO2 discharge to find solutions for reducing SO2 discharge. Literatures about factors affecting SO2 emission can be classified into two categories, direct aspects and indirect aspects. Since atmospheric pollutants are mainly generated from energy consuming, direct forces refer to energy-related causes, while indirect factors imply socio-economic elements which firstly act on energy-related factors and then influence air pollutant emissions. Sinha [8] researched on the influence of energy intensity and economic development on SO2 emissions. Yao et al. [9] employed index decomposition analysis method to analyze the overall industrial scale in the change of SO2 discharge and found that engineering emission reduction and supervision emission reduction contributed the most to reduce SO2 emissions, while structure emissions reduction has not shown a remarkable effect. Han et al. [10] established a new decomposition method combining Grossman decomposition model with logarithmic mean Divisia index (LMDI) method to analyze the factors affecting industrial SO2 emissions. They found that the scale effect was the significant factors influencing SO2 emissions. Literatures above all studied on indirect influencing factors on SO2 discharge, while some studies also researched on direct impact factors. Wang et al. [11] examined the impacts of treatment technology on SO2 discharge for industry employing LMDI method. Yang et al. [12] analyzed the influence of three critical direct elements: energy consumption, energy structure, and treatment technology on SO2 discharge in China employing LMDI. They found that energy consumption was the primary reason for SO2 increasing, while the improvement of treatment technology played a significant role in reducing SO2 emissions.
Empirical studies on environmental Kuznets curve can be retrospected to the work of Grossman and Krueger [13], who firstly found evidence for the existence of the inverted U-shaped nexus between some contaminant indicators and GDP for 42 countries. Since then, the existence of environmental Kuznets curve hypothesis was proved in carbon dioxide emissions by Holtz-Eakin and Selden [14], and Stern et al. [15], in sulfur discharge by Kaufmann et al. [16] as well as List and Gallet [17], in four critical atmospheric pollutants by Selden and Song [18], Harbaugh [19], and Miah [20], and in renewable energy by Angeliki and Konstantinos [21]. Stern and Common [22] examined the environmental Kuznets curve hypothesis for sulfur discharge employing a separate sample and a global sample of high-income Organization for Economic Cooperation and Development countries, and the findings validated the existence of this hypothesis among these countries. Shen [23] found a U-shape curve instead of an inverted U-shape curve between sulfur emissions and GDP based on provincial data of China during the period of 1993–2002. The research of Fodha and Zaghdoud [24] displayed an inverted U-shape curve nexus between sulfur discharge and economic development. Then, series of literatures queried that existing literatures based on income-pollutants environmental Kuznets curve model may provide biased evidence due to the failure to take relevant explanatory variables into consideration [25]. Therefore, some latest literatures try to study the nexus between income and pollutants by taking energy related variables, urbanization, and technique progress variables into the function. Kohler [26] found that there existed a long term relationship between environmental quality, per capita energy use and foreign trade in South Africa. Shahbaz et al. [27] investigated the co-integration relationship and causal relationship among industrialization, electricity consumption, and CO2 emissions of Bangladesh, and found that environmental Kuznets curve existed between industrial process and CO2 discharge. Wang et al. [28] investigated the relationship among economic development, urbanization, and sulfur dioxide discharge employing panel data model based on semi-parametric panel fixed effects regression. The results showed that an inverted U-shape curve existed between economic development and sulfur dioxide discharge. Zhou et al. [7] verified an inverse N-shape relationship between sulfur discharge and economic development, and proved the positive influence of technical progress on sulfur discharge reduction.
Although environmental Kuznets curve analysis model is widely applied in researching the relationship between economic growth and contaminants, many scholars and policymakers put forward their critiques on the concept of environmental Kuznets curve hypothesis and methodology in previous researches. Firstly, the inverse U-shape environmental Kuznets curve is not applicable for all kinds of pollutants, and only small amount of empirical analysis support an inverse U-shape environmental Kuznets curve for some primary atmospheric contaminants [19]. Shafik and Bandyopadhyay [29] examined a panel of 149 countries during the period of 1960–1990, and found that among ten indicators of environmental quality, only two indicators fitted for an environmental Kuznets curve path. Secondly, a large amount of studies on environmental Kuznets curve hypothesis researched cross-sectional data and summarized an exclusive development path for different countries or provinces, which was criticized for the invalidity of cross-sectional technique [30,31].
Above all, considering about the shortcomings of existing literatures on researching the relationship between sulfur emissions and socio-economic driving forces, this paper established a multi-variate panel data model for five provinces with the largest sulfur dioxide (SO2) emissions in China taking economic development, fossil fuel energy consumption, energy consuming intensity, and economic structure into consideration spanning the period of 2002–2015. The main contributions of this paper are as follows:
(1)
Combining economic development, energy consumption, technical advancement, and economic structure together to analyze the contribution of each variable to SO2 emissions. To the best of our knowledge, this paper is the first study in the field of investigating the relationship between sulfur discharge and socio-economic forces to simultaneously explore the contribution of economic growth, energy consumption, technical progress, and economic structure to SO2 emissions using panel data unit root test and panel co-integration theory. Additionally, granger causality test is also employed to investigate the causal relationship between these four data sequences and SO2 emissions.
(2)
The contribution of economic development, energy consumption, technical advancement, and economic structure to SO2 discharge can be quantitatively analyzed regarding to different provinces. Based on panel data model, cross-sectional technique can be fully exploited, thus, the contribution degree of four variables to SO2 discharge can be quantitatively measured according to five different provinces, and the causal relationship direction among these four variables and SO2 discharge can be obtained with respect to different provinces. Therefore, the policymakers of different provinces can formulate effective and practicable policies to reduce the discharge of SO2 according to the empirical analysis results.
The rest parts of this paper are conducted as below. Section 2 majors on introducing the methodology used in this paper and the framework of this research. Section 3 provides the data sources and pre-analysis. The empirical analysis will be carried out in Section 4. Section 5 draws the conclusions and provides policy implications.

2. Theoretical Framework and Methodology

2.1. Test for Cross-Sectional Dependence

The establishment of panel data model should start with cross-sectional dependence test which is necessary for choosing appropriate methods for testing unit root. This paper used the Pesaran cross-sectional dependence test for cross-sectional dependence examination [32]. The original panel data model can be written as:
y i t = α i + β i t x i t + μ i t
where i = 1, 2, …, N and t = 1, 2, …, T, βit demonstrates a vector of parameters for K × 1 to be evaluated, xit indicates a K × 1 vector of explanatory variables, αi implies parameters which do not change with provinces and μit is supposed to be independently and identically distributed. The null and the alternative hypotheses are represented as follow:
H 0 : ρ i j = ρ j i = c o r ( μ i t , μ j t ) = 0   f o r   i j
H 1 : ρ i j = ρ j i 0   f o r   s o m e   i j
where:
ρ i j = ρ j i = t = 1 T μ i t μ j t ( t = 1 T μ i t 2 ) 1 2 ( t = 1 T μ j t 2 ) 1 2
The statistic of the Breusch-Pagan [33] LM test is provided as Equation (5):
L M = i = 1 N 1 j = i + 1 N T i j ρ ^ i 2 χ N ( N 1 ) 2 2
where the ρ ^ i j means the coefficients estimated by the residuals of the model. Pesaran [32] improved the LM test by proposing an alternative test regarding to the average value of ρ ^ i j calculated by:
C D L M = 2 T N ( N 1 ) i = 1 N 1 j = i + 1 N ρ ^ i j N ( 0 , 1 )
Pesaran [32] certified the advancement of this test in small samples which is suitable for this study.

2.2. Test for Panel Unit Root

Panel unit root test methods employed in previous literatures can generally be classified into two categories. The first category takes cross-sectional independence into consideration, such as Hadrid [34], Choi [35], Levin, Lin and Chu (LLC test) [36], Im, Pesaran and Shin (IPS test) [37] and some others. The second category considers cross-sectional dependence, including the test methods proposed by Bai and Ng [38], Moon and Perron [39], Pesaran [40], Phillips and Sul [41], and Smith et al. [42]. The equation applied to test stationary is listed as Equation (7):
Δ y i t = ρ i y i t 1 + δ i X i t + ε i t
where i = 1, 2, …, N is used to represent province; t = 1, 2, …, T demonstrates time point; Xit indicates the explanatory variables containing fixed effects or individual time trend; ρi means the coefficient for auto-regression; and εit implies interference term of stable series.
As Equation (7) may be an autocorrelation, higher order differential delay terms were explored by Levin et al. [36]:
Δ y i t = ρ i y i t 1 + δ i X i t + j = 1 p i θ i j Δ y i t 1 + ε i t
where pi represents the lags amount in the regression. Im et al. [37] verified a t-bar statistic as the average value of the individual ADF statistic as Equation (9):
t b a r = 1 N × i = 1 N t ρ i
where tρi means the individual t-statistic to test for the original hypothesis.

2.3. Test for Panel Co-Integration

If the variables sequences are confirmed to be stable in the same order, the process should step to panel co-integration test employing Pedroni’s co-integration test method [43] of which the regression equation can be written as:
y i t = α i + δ i t + β 1 i x 1 i , t + β 2 i x 2 i , t + K + β M i x M i , t + e i t
where i implies various provinces, t indicates different time points, M represents the number of explanatory variables, β 1 i ,   β 2 i ,   β M i   demonstrates the coefficients of explanatory variables, α i means the intercept component, and e i t is the evaluated residual on behalf of the divergence from the long-run.
Pedroni co-integration examination method is made up of two groups. One is based on the within dimension approach including panel v-statistic, panel ρ-statistic, panel PP-statistic and panel ADF-statistic. The other one is based on the between dimension approach consisted of group ρ-statistic, group PP-statistic and group ADF statistic.

2.4. Test for Panel Data Model Form

The primary three forms for panel data model are random effects, fixed effects, and regression models for pooled. The Hausman examination method and the likelihood ratio (LR) examination method are usually employed to test for the proper form of panel data model [44].
The forms of panel data models contain the constant intercepts and coefficients model, the variable intercepts and invariable coefficients model, and the changed intercepts and coefficients model [45]. These three forms of panel data models are shown in Equations (11)–(13).
y i t = α + β x i t + μ i t
y i t = α i + β x i t + μ i t
y i t = α i + β i x i t + μ i t
where i represents the provinces, t demonstrates the time point, α i means the intercept, β i indicates coefficient, and μ it implies the error component.
In order to verify the panel data model form, F-test will be utilized to decide if the following two null hypotheses should be accepted or rejected by computing the residual sum of squares (RSS) of Equations (11)–(13).
H 1 : β 1 = β 2 = = β N F 1 = ( S 2 S 1 ) / [ ( N 1 ) k ] S 1 / [ N T N ( k + 1 ) ] F ( ( N 1 ) k , N ( T k 1 ) )
H 2 : α 1 = α 2 = = α N , β 1 = β 2 = = β N F 2 = ( S 3 S 1 ) / [ ( N 1 ) ( k + 1 ) ] S 1 / [ N T N ( k + 1 ) ] F ( ( N 1 ) ( k + 1 ) , N ( T k 1 ) )
where F1 is for the H1 hypothesis of which the slopes are invariable and the intercepts are different, F2 is for the H2 hypothesis of which the slopes and the intercepts are all unchanged, S1, S2, and S3 are the residual sum of squares of Equations (11)–(13). Additionally, N, T, and k indicate the amount of provinces, years, and explanatory variables.
Considering the significance level and T > k + 1, if F2 is less than the critical value, then we accept H2 and the panel data model form should be Equation (11), otherwise, it needs to examine hypothesis H1. If F1 is more than the critical value, H1 should be rejected and the panel data model form should be Equation (13), otherwise, it should be Equation (12).

2.5. Test for Causality

Engle and Granger [46] proposed that if two data sequences are co-integrated, then there exists Granger Causal relationship. At the aim of investigating causality among different data sequences, the Granger causality examination method [47] is utilized to testify if one data sequence has an impact on another. For Granger causality, if Y can be forecasted more precisely through using the data of both X and Y than using Y, then we can conclude that the variable X Granger causes Y. This method is employed in this paper to identify the causal relationships among different variables:
y t = α 0 + i = 1 m α i y t i + i = 1 m β i x t i + e t
x t = α 0 + j = 1 n α j y t j + j = 1 n β j x t j + e t
Equation (18) shows the null hypothesis of the Granger causality examination method, which demonstrates “X does not Granger-cause Y”. Equation (19) is used to test if Y Granger-causes X.
H 0 : β i = 0 , i = 1 , 2 , , m
H 0 : β j = 0 ,   j = 1 , 2 , , n

2.6. Theoretical Framework

The theoretical framework is demonstrated in Figure 1. The empirical analysis can be proceeded with following steps.
Step 1: Test for cross-sectional dependence
The empirical analysis will start with cross-sectional dependence test to identify the methods used in unit root test and panel co-integration test.
Step 2: Test for panel unit root
After confirming whether it is necessary to consider cross-sectional dependence, the stable of all variables need to be examined and the methods used in this stage depend on the results of step 1. Only if all variables are stable in the same order, the empirical analysis can proceed to co-integration test.
Step 3: Test for panel co-integration
After verifying all variables are stationary at the same order, Pedroni co-integration examination method will be used for testing whether there exists long-term relationship among SO2 emissions and all independent variables. If not, the panel data model cannot be established, otherwise, we need to determine the form of panel data model.
Step 4: Test for panel data model form
At this step, LR test and Hausman test are used to determine the fixed effect or random effect of panel data model. And then F-test will be employed to judge whether the panel data model is constant intercepts and coefficients model, the variable intercepts and invariable coefficients model, or the changed intercepts and coefficients model. After identifying the model form, the panel data model can be estimated.
Step 5: Test for Granger causal relationship
At the aim of further and better understanding the relationship among SO2 emissions, GDP, fossil fuel energy consumption, energy consumption intensity, and economic structure, Granger causality examination method will be used to explore the nexus among these variables. And the direction of causal nexus can provide policy making references for policy makers.

3. Data Sources and Pre-Analysis

3.1. Study Area and Data Sources

At the aim of exploring the contributions of various socio-economic factors to SO2 discharge and specifically carrying out policy recommendations to reduce SO2 emissions, this paper investigates the contributions of economic development, fossil fuel energy consumption, energy consuming intensity, and economic structure to SO2 discharge using the data of five provinces with the highest SO2 emissions in China during the period of 2002–2015. The provinces selected in this paper are Hebei, Henan, Inner Mongolia, Shandong, and Shanxi.
For economic development, gross domestic production (GDP) is selected to represent the growth of economy, and is converted to the constant price using 2002 as basic period. For fossil fuel energy consumption, it is made up of the use amount of coal, crude oil, and natural gas. The statistical unit is 104 ton for coal and crude oil consumption, and 108 cubic meter (m3) for natural gas consumption. On the purpose of maintaining consistency of statistical requirements, the consumption units are converted to 104 ton coal equivalent (tce) employing the coefficients suggested in the China Energy Statistic Yearbook, which are 0.7143 tce/t, 1.4286 tce/t, and 1.33 tce/103 m3 for the consuming of coal, crude oil, and natural gas [48,49,50,51,52,53,54,55,56,57,58,59,60,61,62,63,64,65,66,67]. For energy consuming intensity, it is represented by the value of energy consumption for the whole society divided by GDP (constant price taking 2002 as basic period). For economic structure, it is expressed by the value of the added value for secondary industry divided by GDP due to the high proportion of added value of secondary industry in GDP. All data of the variables mentioned above as well as the data of SO2 emissions are collected from the official website of National Statistics Bureau of China. This dataset on SO2 emissions in provincial level allowed this research from provincial aspects over a long period.
The representation forms of all independent variables and SO2 discharge used in the panel data model are listed in Table 1. The data of SO2 discharge, GDP, fossil fuel energy consumption, and energy consuming intensity are transformed into natural logarithmic form. The data of the economic structure is multiplied by 100 before being converted into the natural logarithmic form.

3.2. Pre-Analysis

Since the high level of SO2 discharge can cause acid rain weather which seriously damage to human health, it is necessary to research on the contributions of significant socio-economic factors on SO2 emissions to provide references for policy makers. Figure 2 displays the spatial distribution of SO2 discharge in 32 primary provinces of China in the year of 2010 and 2015. As can be seen from the figure, the SO2 emissions amount of Inner Mongolia, Shanxi, Hebei, Shandong, Henan, Sichuan, and Guizhou in 2010 were more than 1100 thousand tons, which contributed more than 40% of SO2 emissions in China. During the period of 2011–2015, the SO2 emissions amount of most provinces showed downward trend. However, the amount of SO2 discharge in Hebei, Henan, Inner Mongolia, Shandong, and Shanxi still maintained more than 1100 thousand tons, among which Shandong ranked the first with 1525.67 thousand tons and Hebei ranked the fifth with 1108.37 thousand tons. Therefore, these five provinces are selected to be the objects of this study aiming at better understanding the relationships between SO2 emissions and significant socio-economic indicators in these five provinces.
The development trend of GDP for the selected five provinces from 2002 to 2015 is demonstrated in Figure 3. As indicated in the figure, the economy scale of Shandong ranks the first, with 1027.55 billion yuan in 2002 and 4486.34 billion yuan in 2015. Before 2009, the scale of economic development in Hebei and Henan was close to each other, while after 2009, the growth rate of Hebei lagged behind Henan. The economic scales of Inner Mongolia and Shanxi are much smaller than other three provinces, with 194.09 billion yuan for Inner Mongolia and 232.48 billion yuan for Shanxi in 2002, and 1226.40 billion yuan and 856.93 billion yuan for Inner Mongolia and Shanxi in 2015, respectively.
The development tendency of fossil fuel energy consumption from 2002 to 2015 for the chosen five provinces is indicated in Figure 4. With the largest economic scale, the fossil fuel energy consumption of Shandong also ranked the first with 495.34 million tce in 2015 and 10.28% annual growth rate. Although the economic scale of Shanxi ranked the fifth, the fossil fuel energy consumption of it ranked the second with 378.68 million tce in 2015, which means the development of economy in Shanxi mainly depends on the consumption of fossil fuel energy. The amount of fossil fuel energy consumption in Inner Mongolia increased fast from 69.90 million tce in 2002 to 368.83 million tce in 2015 with 14.14% annual growth rate ranking the third. The fossil fuel energy consumption of Hebei and Henan showed a downtrend from 323.57 million tce for Hebei and 292.49 million tce for Henan in 2011 to 306.10 million tce for Hebei and 245.67 million tce for Henan in 2015.
The development trend of energy consuming intensity from 2002 to 2015 for five provinces is illustrated in Figure 5. As depicted in the figure, the energy consuming intensity of Hebei, Henan, Inner Mongolia, Shandong and Shanxi all showed decrease tendency from 2002 to 2015. Based on the analysis of GDP and fossil fuel energy consumption, we can conclude that the energy use efficiency of Shanxi was the lowest, thus the energy consuming intensity of Shanxi should be the highest. It can also be obtained from the figure that the energy consuming intensity of Shanxi ranked the first which is much higher than other four provinces.
The development tendency of economic structure for the selected five provinces from 2002 to 2015 is implied in Figure 6. As displayed in the figure, at the beginning of the analysis phase, the proportion of secondary industry added value accounted for GDP in Shandong was the highest with 50.46%, and then it reduced to 46.80% in 2015 ranking the fourth. The proportion of secondary industry added value accounted for GDP in Shanxi ranked the second in 2002, and then it decreased to 40.69% in 2015 ranking the fifth. The proportion of secondary industry added value accounted for GDP in Inner Mongolia, Henan and Hebei first increased to 55.97% in 2011, 57.28% in 2010, and 54.34% in 2008, and then dropped to 50.48%, 48.42%, and 48.27% in 2015.

4. Empirical Analysis

4.1. Results for Cross-Sectional Dependence Examination

The initial stage of empirical analysis should be cross-sectional dependence examination which decides the methods selected to test unit root. The results of Pesaran cross-sectional dependence examination method utilized in this study are shown in Table 2. As can be seen from p-value in Table 2, the null hypothesis can be rejected at 5% confidence level. Therefore, the methods employed to test data stationary and co-integration relationship should take cross-sectional dependence into consideration.

4.2. Results for Unit Root Examination

Since the following analysis should take cross-sectional dependence into consideration, at the stage of testing data series stationary, L.L&C, IPS, Augmented Dickey Fuller-Fisher (ADF-Fisher) and Phillips-Perron Fisher (PP-Fisher) examination approaches are selected. As shown in Table 3, regarding to the probability values in the brackets of different examination methods for variables, we can obtain that all variables are not stable in their level form. Then all variables need to be differenced. ΔlnSO2, ΔlnGDP, ΔlnFEC, ΔlnECI, and ΔlnES are utilized to represent the first difference form for SO2 emissions, GDP, fossil fuel energy consumption, energy consumption intensity, and economic structure in natural logarithmic form. Since all of the probability values in brackets are less than the confidence level, SO2 emissions and four independent variables are stable in the same order.

4.3. Results for Panel Co-Integration Examination

As all variables are first order differential stationary, the empirical analysis can step to panel co-integration test. Pedroni’s co-integration examination approach is utilized in this stage and the statistic values of it are illustrated in Table 4. Since all values of probability are less than the confidence level, it can be verified that the long-term co-integration relationship is existed among SO2 discharge, GDP, fossil fuel energy consumption, energy consumption intensity, and economic structure.

4.4. Model Form Determination

As there exists long-run co-integration relationship among all variables, the empirical analysis can proceed to model estimation. First, the effect of panel data model needs to be identified. Generally, LR test and Hausman test are used to determine whether the panel data model is fixed effect or random effect. The test results are shown in Table 5. With regard to LR examination approach, the probability value is less than the significance level, which implies that the panel data model should be fixed effect. As for Hausman examination results, based on the value of cross-section random and probability values of all independent variables, it is verified that the panel data model should be fixed effect.
Then the panel data model form of the unchanged intercepts and coefficients model, the changed intercepts and unchanged coefficients model, and the changed intercepts and coefficients model need to be determined employing F-test. Firstly, we need to obtain three sum square residual values represented by S1, S2, and S3 of Equations (11)–(13), respectively. Secondly, we can calculate the values of F1 and F2 statistics based on Equations (14) and (15). Thirdly, the panel data model form can be confirmed. If the value of F2 is smaller than the critical value F2,α((N − 1)(K + 1),(NTN(K + 1)), the panel data model is confirmed to be the model with constant intercepts and coefficients, if not, the hypothesis H1 will be examined. If the value of F1 is larger than the critical value F1,α((N − 1)K,(NTN(K + 1)), the panel data model is verified to be the model with changed intercepts and coefficients, if not, it is judged to be the model with variable intercepts and invariable coefficients. In accordance with the results for F-test, both of the values of F1 and F2 statistics are more than the critical values at the given significance level. Therefore, the established panel data model should be the changed intercepts and coefficients model as Table 6 shown.

4.5. Estimation for Panel Data Model and Provincial Comparative Analysis

Based on the above examination process, a fixed effect panel data model with changed intercepts and coefficients consists of SO2 discharge, GDP, fossil fuel energy consumption, energy consumption intensity, and economic structure can be established for Hebei, Henan, Inner Mongolia, Shandong, and Shanxi. The estimated coefficients of the panel data model and the test results for the effectiveness of the model are listed in Table 7. The value of R2 is 0.9981 which means the fitting effect of the model is pretty good. The value of F-statistic is 129.0515, which is much more than the critical value representing the coefficients of the model are all significant. Therefore, it can be confirmed that the established panel data model is valid and significant.
As all variables are converted to logarithmic form, the coefficients represent elasticities. Through analyzing the coefficients, we can conclude that:
(1)
Compared with GDP, energy consumption intensity, and economic structure, fossil fuel energy consumption makes the greatest contribution to SO2 discharge. Since fossil fuel energy contain large amount of sulfur element, the combustion of fossil fuel energy will release a great deal of sulfide which is the main source of sulfur dioxide. Through comparing the different contributions of fossil fuel energy consumption to SO2 emissions for different provinces, it can be summarized that the more the fossil fuel energy consumption, the greater the contribution made by it to SO2 emissions. Since the fossil fuel energy consumption amount of Shandong ranked the first, the contribution of fossil fuel energy consumption to SO2 emissions is the greatest, followed by Shanxi. The contribution of fossil fuel energy consumption to SO2 emissions in Henan is the smallest.
(2)
The contribution of GDP to SO2 emissions is less than fossil fuel energy consumption, but more than energy consumption intensity and economic structure. Compared with the contribution of fossil fuel energy consumption to SO2 emissions, GDP make much less contribution to SO2 emissions. Through comparing the different devotions of GDP to SO2 emissions among five selected provinces, it can be concluded that the larger the scale of the economy, the greater the devotion made by it to SO2 emissions. Since the economy scale of Shandong ranked the first, the devotion of GDP to SO2 emissions in Shandong is the greatest, followed by Henan. GDP in Shanxi contributes the smallest to SO2 emissions.
(3)
The more the proportion of the secondary industry added value accounted in GDP, the greater the contribution of the economic structure and energy consumption intensity to SO2 emissions. From the overall trend of the proportion for secondary industry added value accounted in GDP, the proportions of secondary industry added value accounted in GDP for Shandong and Inner Mongolia are higher than other three provinces, and the contributions of energy consumption intensity and economic structure to SO2 emissions of Shandong and Inner Mongolia are higher than Hebei, Henan, and Shanxi.

4.6. Analysis for Granger Causality Relationship

The results of Granger causality test for the established panel data model are listed in Table 8. If the probability values in brackets are less than the specific significance level, there exists Granger causality relationship between these two variables. As illustrated in Table 8, we can summarize that:
(1)
For the selected five provinces, there exists a bi-directional causality relationship between fossil fuel energy consumption and SO2 emissions. This indicates that the decrease of fossil fuel energy consumption can contribute to the reduction of SO2 emissions, while the decrease of SO2 emissions will have a negative impact on fossil fuel energy consumption.
(2)
There exist uni-directional causal relationships running from GDP to SO2 emissions, from energy consumption intensity to SO2 emissions, and from economic structure to SO2 emissions. That means the increase of GDP will lead to the raise of SO2 emissions, the decrease of energy consumption intensity will bring about the decline of SO2 emissions, and the decrease of the proportion for the secondary industry added value accounted in GDP will result in the reduction of SO2 emissions. But there is no feedback from SO2 emissions to GDP, energy consumption intensity, and economic structure.

5. Conclusions and Policy Implications

This paper investigated the relationship among SO2 emissions, GDP, fossil fuel energy consumption, energy consumption intensity, and economic structure of five provinces in China with the highest SO2 emissions employing panel data model approach spanning from 2002 to 2015. Based on panel unit root examination to confirm all variables are stable in the same order, panel co-integration examination to verify all variables are co-integrated in the long-term, and model form determination tests, the panel data model can be estimated and the Granger causality examination can be conducted, which provided strong evidence on the complicated relationships among these variables. Through analyzing the established panel data model, the main conclusions are as follows:
(1)
Fossil fuel energy consumption makes the greatest contribution to SO2 discharge compared with GDP, energy consumption intensity, and economic structure. And the more the fossil fuel energy consumption, the greater the devotion made by it to SO2 emissions.
(2)
GDP makes less contribution to SO2 emissions than fossil fuel energy consumption, and the larger the scale of the economy, the greater the contribution made by it to SO2 emissions.
(3)
The higher the proportion of the secondary industry added value accounted in GDP, the more the devotion made by the economic structure and energy consumption intensity to SO2 emissions.
Through investigating the causal relationship among SO2 emissions, GDP, fossil fuel energy consumption, energy consumption intensity, and economic structure of Hebei, Henan, Inner Mongolia, Shandong, and Shanxi using Granger causality examination method, we can obtain that:
(1)
A bi-directional causal relationship exists between fossil fuel energy consumption and SO2 emissions among five selected provinces.
(2)
Uni-directional causal nexus exist running from GDP to SO2 emissions, from energy consumption intensity to SO2 emissions, and from economic structure to SO2 emissions among five chosen provinces.
The empirical analysis of this research can make policy makers better understand the complex relationships between SO2 emissions, GDP, fossil fuel energy consumption, energy consumption intensity, and economic structure of Hebei, Henan, Inner Mongolia, Shandong, and Shanxi. In this way, the significant influencing factors of SO2 discharge can be found, and effective and practicable policies can be formulated according to the contributions of different factors and causality relationships. Therefore, regarding to the above econometric analysis, the following recommendations are put forward:
(1)
Giving full play to the guiding role of price signals, and improving the price policy for desulfurization. Based on the empirical analysis, it can be found that the consumption of fossil fuel energy makes the greatest contribution to SO2 emissions. However, in the long-run, energy consumption of various provinces in China will still depend on fossil fuel energy to a large extent. Therefore, in the process of using fossil fuel energy, desulfurization equipment is needed to reduce sulfur emissions. Considering about the increase of the cost for enterprises using desulfurization equipment, it is necessary to implement a certain price subsidy policy so that the increased cost due to the use of desulfurization technique can be covered, which can also encourage enterprises to develop desulfurization technology.
(2)
Formulating a new comprehensive evaluation indicator to measure the regional development level considering economic development and environmental impacts. Currently, GDP is deemed as the significant indicator to evaluate the regional development level, which neglects the environmental impact. Therefore, it is essential to establish a comprehensive evaluation system for measuring regional development level, which can not only consider the level of economic development, but also take the impact of pollutant emissions on the environment in the process of rapid economic development into consideration.
(3)
Exploring renewable and sustainable energy sources to substitute for fossil fuel energy. Considering the high pollution and limited nature of fossil fuel energy, we should actively exploit other renewable energy to take over the use of fossil resources. Based on regional resources endowment, people in inland areas should energetically develop the use of wind energy and solar energy, such as Hebei, Henan, Inner Mongolia, and Shanxi, while people in coastal areas should positively explore the use of hydropower and tidal energy, such as Shandong.
(4)
Developing high value added and low pollution emissions industries and reducing the proportion of secondary industry. According to the empirical analysis, we can obtain that the higher the proportion of the secondary industry accounted in GDP, the more the devotion made by the economic structure and energy consumption intensity to SO2 emissions. Therefore, policy makers should draft related policies for economic structure adjustment, which should aims at reducing the proportion of secondary industry and developing high value added and low pollution emissions industries. Additionally, policies related to improve energy using efficiency should be executed.
Based on the above analysis, although the results of this study are inspiring, we also have lots of work needed to be done in the future. The policy effect of reducing SO2 emissions should be quantified and added to the panel data model to analyze the function of policies in decreasing SO2 discharge. Additionally, it is also necessary to analyze the influence of renewable energy and its relevant policies on SO2 emissions. Therefore, policymakers can grasp the policy effect, and propose more effective policies and measures.

Acknowledgments

This research is sponsored by the National Natural Science Foundation of China under Grant No. 71373076, the Fundamental Research Funds for the Central Universities under Grant No. 2017MS060, the Beijing Social Science Foundation Project under Grant No. 17GLC042, and the Beijing Excellent Talents Training Program under Grant No. 2017000020124G189.

Author Contributions

Huiru Zhao proposed the theoretical framework of this investigation. Haoran Zhao completed the paper. Sen Guo modified the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhu, L.; Gan, Q.; Liu, Y.; Yan, Z. The impact of foreign direct investment on SO2 emissions in the Beijing-Tianjin-Hebei region: A spatial econometric analysis. J. Clean. Prod. 2017, 166, 189–196. [Google Scholar] [CrossRef]
  2. Wang, H.; An, J.; Shen, L.; Zhu, B.; Pan, C.; Liu, Z.; Liu, X.; Duan, Q.; Liu, X.; Wang, Y. Mechanism for the formation and microphysical characteristics of submicron aerosol during heavy haze pollution episode in the Yangtze River Delta, China. Sci. Total Environ. 2014, 490, 501–508. [Google Scholar] [CrossRef] [PubMed]
  3. Wang, Y.; Ying, Q.; Hu, J.; Zhang, H. Spatial and temporal variations of six criteria air pollutants in 31 provincial capital cities in China during 2013–2014. Environ. Int. 2014, 73, 413–422. [Google Scholar] [CrossRef] [PubMed]
  4. Hsu, A. 2016 Environmental Performance Index; Yale University: New Haven, CT, USA, 2016. [Google Scholar]
  5. Qu, Y.; An, J.; He, Y.; Zheng, J. An overview of emissions of SO2 and NOx and the long-range transport of oxidized sulfur and nitrogen pollutants in East Asia. J. Environ. Sci. 2016, 44, 13–25. [Google Scholar] [CrossRef] [PubMed]
  6. Zhang, H.; Wang, Z.; Zhang, W. Exploring spatiotemporal patterns of PM2.5 in China based on ground-level observations for 190 cities. Environ. Pollut. 2016, 216, 559–567. [Google Scholar] [CrossRef] [PubMed]
  7. Zhou, Z.; Ye, X.; Ge, X. The Impacts of Technical Progress on Sulfur Dioxide Kuznets Curve in China: A Spatial Panel Data Approach. Sustainability 2017, 9, 674. [Google Scholar] [CrossRef]
  8. Sinha, A. Trilateral association between SO2/NO2 emission, inequality in energy intensity, and economic growth: A case of Indian cities. Atmos. Pollut. Res. 2016, 7, 647–658. [Google Scholar] [CrossRef]
  9. Yao, L.; Liu, J.; Zhou, T.; Wang, R. An analysis of the driving forces behind pollutant emission reduction in Chinese industry. J. Clean. Prod. 2016, 112, 1395–1400. [Google Scholar] [CrossRef]
  10. Han, Y.; Xiong, X.; Liu, Y.; Pan, Y.R.; Zhang, Y.B. The Analysis of Factors Affecting SO2 Emission of Chinese Industry. In Proceedings of the 2011 International Conference on Computer and Management (CAMAN), Wuhan, China, 19–21 May 2011; pp. 1–11. [Google Scholar]
  11. Wang, Q.; Wang, Y.; Zhou, P.; Wei, H. Whole process decomposition of energy-related SO2 in Jiangsu Province, China. Appl. Energy 2017, 194, 679–687. [Google Scholar] [CrossRef]
  12. Yang, X.; Wang, S.; Zhang, W.; Li, J.; Zou, Y. Impacts of energy consumption, energy structure, and treatment technology on SO2 emissions: A multi-scale LMDI decomposition analysis in China. Appl. Energy 2016, 184, 714–726. [Google Scholar] [CrossRef]
  13. Grossman, G.M.; Krueger, A.B. Economic growth and the environment. Q. J. Econ. 1995, 110, 353–377. [Google Scholar] [CrossRef]
  14. Holtz-Eakin, D.; Selden, T.M. Stoking the fires? CO2 emissions and economic growth. J. Public Econ. 1995, 57, 85–101. [Google Scholar] [CrossRef]
  15. Stern, D.I.; Common, M.S.; Barbier, E.B. Economic growth and environmental degradation: The environmental Kuznets curve and sustainable development. World Dev. 1996, 24, 1151–1160. [Google Scholar] [CrossRef]
  16. Kaufmann, R.K.; Davidsdottir, B.; Garnham, S.; Pauly, P. The determinants of atmospheric SO2 concentrations: Reconsidering the environmental Kuznets curve. Ecol. Econ. 1998, 25, 209–220. [Google Scholar] [CrossRef]
  17. List, J.A.; Gallet, C.A. The environmental Kuznets curve: Does one size fit all? Ecol. Econ. 1999, 31, 409–423. [Google Scholar] [CrossRef]
  18. Selden, T.M.; Song, D. Environmental quality and development: Is there a Kuznets curve for air pollution emissions? J. Environ. Econ. Manag. 1994, 27, 147–162. [Google Scholar] [CrossRef]
  19. Harbaugh, W.T.; Levinson, A.; Wilson, D.M. Reexamining the empirical evidence for an environmental Kuznets curve. Rev. Econ. Stat. 2002, 84, 541–551. [Google Scholar] [CrossRef]
  20. Miah, M.D.; Masum, M.F.; Koike, M. Global observation of EKC hypothesis for CO2, SOx and NOx emission: A policy understanding for climate change mitigation in Bangladesh. Energy Policy 2010, 38, 4643–4651. [Google Scholar] [CrossRef]
  21. Menegaki, A.N.; Tsagarakis, K.P. Rich enough to go renewable, but too early to leave fossil energy? Renew. Sustain. Energy Rev. 2015, 41, 1465–1477. [Google Scholar] [CrossRef]
  22. Stern, D.I.; Common, M.S. Is there an environmental Kuznets curve for sulfur? J. Environ. Econ. Manag. 2001, 41, 162–178. [Google Scholar] [CrossRef]
  23. Junyi, S.H. A simultaneous estimation of environmental Kuznets curve: Evidence from China. China Econ. Rev. 2006, 17, 383–394. [Google Scholar]
  24. Fodha, M.; Zaghdoud, O. Economic growth and pollutant emissions in Tunisia: An empirical analysis of the environmental Kuznets curve. Energy Policy 2010, 38, 1150–1156. [Google Scholar] [CrossRef]
  25. Stern, D.I. The rise and fall of the environmental Kuznets curve. World Dev. 2004, 32, 1419–1439. [Google Scholar] [CrossRef]
  26. Kohler, M. CO2 emissions, energy consumption, income and foreign trade: A South African perspective. Energy Policy 2013, 63, 1042–1050. [Google Scholar] [CrossRef]
  27. Shahbaz, M.; Uddin, G.S.; Rehman, I.U.; Imran, K. Industrialization, electricity consumption and CO2 emissions in Bangladesh. Renew. Sustain. Energy Rev. 2014, 31, 575–586. [Google Scholar] [CrossRef]
  28. Wang, Y.; Han, R.; Kubota, J. Is there an environmental Kuznets curve for SO2 emissions? A semi-parametric panel data analysis for China. Renew. Sustain. Energy Rev. 2016, 54, 1182–1188. [Google Scholar] [CrossRef]
  29. Shafik, N.; Bandyopadhyay, S. Economic Growth and Environmental Quality: Time-Series and Cross-Country Evidence; World Bank: Washington, DC, USA, 1992. [Google Scholar]
  30. Unruh, G.C.; Moomaw, W.R. An alternative analysis of apparent EKC-type transitions. Ecol. Econ. 1998, 25, 221–229. [Google Scholar] [CrossRef]
  31. Stern, D.I. Progress on the environmental Kuznets curve? Environ. Dev. Econ. 1998, 3, 173–196. [Google Scholar] [CrossRef]
  32. Pesaran, M.H. General Diagnostic Tests for Cross Section Dependence in Panels; CESifo Working Papers No. 1233; CESifo: München, Germany, 2004; pp. 255–260. [Google Scholar]
  33. Breusch, T.S.; Pagan, A.R. The Lagrange multiplier test and its applications to model specification in econometrics. Rev. Econ. Stud. 1980, 47, 239–253. [Google Scholar] [CrossRef]
  34. Hadri, K. Testing for stationarity in heterogeneous panel data. Econom. J. 2000, 3, 148–161. [Google Scholar] [CrossRef]
  35. Choi, I. Unit root tests for panel data. J. Int. Money Financ. 2001, 20, 249–272. [Google Scholar] [CrossRef]
  36. Levin, A.; Lin, C.F.; Chu, C.S.J. Unit root tests in panel data: Asymptotic and finite-sample properties. J. Econom. 2002, 108, 1–24. [Google Scholar] [CrossRef]
  37. Im, K.S.; Pesaran, M.H.; Shin, Y. Testing for unit roots in heterogeneous panels. J. Econom. 2003, 115, 53–74. [Google Scholar] [CrossRef]
  38. Bai, J.; Ng, S. A PANIC attack on unit roots and cointegration. Econometrica 2004, 72, 1127–1177. [Google Scholar] [CrossRef]
  39. Moon, H.R.; Perron, B. Testing for a unit root in panels with dynamic factors. J. Econom. 2004, 122, 81–126. [Google Scholar] [CrossRef]
  40. Pesaran, M.H. A simple panel unit root test in the presence of cross-section dependence. J. Appl. Econom. 2007, 22, 265–312. [Google Scholar] [CrossRef]
  41. Phillips, P.C.B.; Sul, D. Dynamic panel estimation and homogeneity testing under cross section dependence. Econom. J. 2003, 6, 217–259. [Google Scholar] [CrossRef]
  42. Smith, L.V.; Leybourne, S.; Kim, T.H.; Newbold, P. More powerful panel data unit root tests with an application to mean reversion in real exchange rates. J. Appl. Econom. 2004, 19, 147–170. [Google Scholar] [CrossRef]
  43. Pedroni, P. Panel Cointegration: A Symptotic and Finite Sample Properties of Pooled Time Series Tests with an Application to the PPP Hypothesis. Econ. Theory 2004, 20, 597–625. [Google Scholar] [CrossRef]
  44. Ou, J.; Liu, X.; Li, X.; Chen, Y. Quantifying the relationship between urban forms and carbon emissions using panel data analysis. Landsc. Ecol. 2013, 28, 1889–1907. [Google Scholar] [CrossRef]
  45. Wang, S.; Fang, C.; Wang, Y.; Huang, Y.; Ma, H. Quantifying the relationship between urban development intensity and carbon dioxide emissions using a panel data analysis. Ecol. Indic. 2015, 49, 121–131. [Google Scholar] [CrossRef]
  46. Engle, R.F.; Granger, C.W.J. Co-integration and error correction: Representation, estimation, and testing. Econometrica 1987, 55, 251–276. [Google Scholar] [CrossRef]
  47. Granger, C.W.J. Investigating causal relations by econometric models and cross-spectral methods. Econom. Econom. 1969, 37, 424–438. [Google Scholar] [CrossRef]
  48. National Bureau of Statistics of China. China Energy Statistical Yearbook 1996; China Statistics Press: Beijing, China, 1996.
  49. National Bureau of Statistics of China. China Energy Statistical Yearbook 1997; China Statistics Press: Beijing, China, 1997.
  50. National Bureau of Statistics of China. China Energy Statistical Yearbook 1998; China Statistics Press: Beijing, China, 1998.
  51. National Bureau of Statistics of China. China Energy Statistical Yearbook 1999; China Statistics Press: Beijing, China, 1999.
  52. National Bureau of Statistics of China. China Energy Statistical Yearbook 2000; China Statistics Press: Beijing, China, 2000.
  53. National Bureau of Statistics of China. China Energy Statistical Yearbook 2001; China Statistics Press: Beijing, China, 2001.
  54. National Bureau of Statistics of China. China Energy Statistical Yearbook 2002; China Statistics Press: Beijing, China, 2002.
  55. National Bureau of Statistics of China. China Energy Statistical Yearbook 2003; China Statistics Press: Beijing, China, 2003.
  56. National Bureau of Statistics of China. China Energy Statistical Yearbook 2004; China Statistics Press: Beijing, China, 2004.
  57. National Bureau of Statistics of China. China Energy Statistical Yearbook 2005; China Statistics Press: Beijing, China, 2005.
  58. National Bureau of Statistics of China. China Energy Statistical Yearbook 2006; China Statistics Press: Beijing, China, 2006.
  59. National Bureau of Statistics of China. China Energy Statistical Yearbook 2007; China Statistics Press: Beijing, China, 2007.
  60. National Bureau of Statistics of China. China Energy Statistical Yearbook 2008; China Statistics Press: Beijing, China, 2008.
  61. National Bureau of Statistics of China. China Energy Statistical Yearbook 2009; China Statistics Press: Beijing, China, 2009.
  62. National Bureau of Statistics of China. China Energy Statistical Yearbook 2010; China Statistics Press: Beijing, China, 2010.
  63. National Bureau of Statistics of China. China Energy Statistical Yearbook 2011; China Statistics Press: Beijing, China, 2011.
  64. National Bureau of Statistics of China. China Energy Statistical Yearbook 2012; China Statistics Press: Beijing, China, 2012.
  65. National Bureau of Statistics of China. China Energy Statistical Yearbook 2013; China Statistics Press: Beijing, China, 2013.
  66. National Bureau of Statistics of China. China Energy Statistical Yearbook 2014; China Statistics Press: Beijing, China, 2014.
  67. National Bureau of Statistics of China. China Energy Statistical Yearbook 2015; China Statistics Press: Beijing, China, 2015.
Figure 1. Theoretical framework for researching the relationship among SO2 emissions, GDP, fossil fuel energy consumption, energy consumption intensity, and economic structure.
Figure 1. Theoretical framework for researching the relationship among SO2 emissions, GDP, fossil fuel energy consumption, energy consumption intensity, and economic structure.
Sustainability 10 00657 g001
Figure 2. The spatial distribution for SO2 discharge in China in the year of 2010 (A) and 2015 (B). Note: NX for Ningxia, SX for Shaanxi, SHX for Shanxi, HB for Hebei, SD for Shandong, JS for Jiangsu, HEN for Henan, AH for Anhui, HUB for Hubei, ZJ for Zhejiang, CQ for Chongqing, GZ for Guizhou, HUN for Hunan, JX for Jiangxi, FJ for Fujian, GX for Guangxi, GD for Guangdong, TW for Taiwan, and HN for Hainan.
Figure 2. The spatial distribution for SO2 discharge in China in the year of 2010 (A) and 2015 (B). Note: NX for Ningxia, SX for Shaanxi, SHX for Shanxi, HB for Hebei, SD for Shandong, JS for Jiangsu, HEN for Henan, AH for Anhui, HUB for Hubei, ZJ for Zhejiang, CQ for Chongqing, GZ for Guizhou, HUN for Hunan, JX for Jiangxi, FJ for Fujian, GX for Guangxi, GD for Guangdong, TW for Taiwan, and HN for Hainan.
Sustainability 10 00657 g002
Figure 3. The development trend of GDP for these five provinces (Unit: billion yuan).
Figure 3. The development trend of GDP for these five provinces (Unit: billion yuan).
Sustainability 10 00657 g003
Figure 4. The development tendency of fossil fuel energy consumption for these five provinces (Unit: million tce).
Figure 4. The development tendency of fossil fuel energy consumption for these five provinces (Unit: million tce).
Sustainability 10 00657 g004
Figure 5. The development trend of energy consuming intensity for these five provinces (Unit: tce/million yuan).
Figure 5. The development trend of energy consuming intensity for these five provinces (Unit: tce/million yuan).
Sustainability 10 00657 g005
Figure 6. The development tendency of economic structure for these five provinces (Unit: %).
Figure 6. The development tendency of economic structure for these five provinces (Unit: %).
Sustainability 10 00657 g006
Table 1. The representation forms of all variables in panel data model.
Table 1. The representation forms of all variables in panel data model.
VariablesDefinitionUnit
SO2The amount of sulfur dioxide emissionsThousand tons
GDPThe development of economyBillion yuan
FECThe amount of fossil fuel energy consumptionMillion tce
ECIEnergy consumption intensity for the whole societyTce/million yuan
ESThe proportion of secondary industry added value accounted for GDP%
Table 2. Results for cross-sectional dependence test.
Table 2. Results for cross-sectional dependence test.
Cross-Sectional Dependence TestPesaran’s Test of Cross-Sectional Independencep-Value
Pesaran cross-sectional dependence test3.86870.0147 b
Notes: The probability values less than the confidence level demonstrates the null hypothesis is rejected. b demonstrates 5% significance level.
Table 3. Results for panel unit root test.
Table 3. Results for panel unit root test.
FormVariablesL.L&CIPSADF-FisherPP-FisherConclusions
LevellnSO20.864820.402400.780351.50630Non-stationary
(0.8064)(0.6563)(0.6769)(0.4709)
lnGDP0.634422.024640.030460.00020Non-stationary
(0.7371)(0.9785)(0.9849)(0.9999)
lnFEC−1.05854−0.294162.037331.33161Non-stationary
(0.1449)(0.3843)(0.3611)(0.5139)
lnECI1.459602.511240.014160.00313Non-stationary
(0.9278)(0.9940)(0.9929)(0.9984)
lnES3.478831.859370.060780.19242Non-stationary
(0.9997)(0.9685)(0.9701)(0.9083)
First DifferencesΔlnSO2−2.01509−1.528745.665586.09979Stationary
(0.0219) b(0.0432) b(0.0488) b(0.0474) b
ΔlnGDP−3.88122−1.988577.2169015.1666Stationary
(0.0001) a(0.0234) b(0.0271) b(0.0005) a
ΔlnFEC−2.87180−2.465668.4611714.4327Stationary
(0.0020) a(0.0014) a(0.0252) b(0.0109) b
ΔlnECI−3.87144−2.526068.8446513.1489Stationary
(0.0001) a(0.0058) a(0.0120) b(0.0014) a
ΔlnES−2.24266−2.016957.591608.75067Stationary
(0.0125) b(0.0219) b(0.0225) b(0.0126) b
Notes: probability values are listed in the brackets. The probability values less than the confidence level demonstrates the variable is stationary. a illustrates 1% significance level; b demonstrates 5% significance level.
Table 4. Results for panel co-integration test.
Table 4. Results for panel co-integration test.
Test MethodTest StatisticsValueProbability
PedroniPanel v-Statistic−1.8454220.0479 b
Panel rho-Statistic−1.0253750.0410 b
Panel PP-Statistic−7.5171280.0000 a
Panel ADF-Statistic−1.9901370.0211 b
Group rho-Statistic−1.9370890.0126 b
Group PP-Statistic−8.5446970.0000 a
Group ADF-Statistic−1.5365370.0295 b
Notes: a points to 1% level of confidence; b refers to 5% confidence level.
Table 5. Results for panel data model effect test.
Table 5. Results for panel data model effect test.
LR Test Results
StatisticProb.
Cross-section F14.3823670.0001 a
Hausman Test Results
Chi-Square StatisticProb.
Cross-Section Random38.4531270.0000 a
VariableFixedRandomVar(Diff.)Prob.
GDP1.4720341.1420880.0238670.0000 a
FEC1.5039790.5053600.0274170.0015 a
ECI−1.3097880.0349780.0031470.0027 a
ES0.9763410.7364150.0041370.0059 a
Notes: a refers to 1% significance level.
Table 6. Results for F-test.
Table 6. Results for F-test.
S1S2S3F1F2
0.0210350.0578920.0903734.9279927.416710
Table 7. The estimated results for panel data model.
Table 7. The estimated results for panel data model.
VariablesHebeiHenanInner MongoliaShandongShanxi
lnGDP0.63140.65250.61980.73780.5928
(0.0368) b(0.0075) a(0.0332) b(0.0023) a(0.0022) a
lnFEC0.88940.80390.87670.93120.8905
(0.0042) a(0.0250) b(0.0451) b(0.0089) a(0.0004) a
lnECI0.39070.38850.41030.43310.3098
(0.0423) b(0.0413) b(0.0399) b(0.0385) b(0.0275) b
lnES0.44880.44250.45230.48280.3364
(0.0469) b(0.0421) b(0.0327) b(0.0069) a(0.0016) a
C1.06912.50341.00534.07753.6094
R20.9981F-statistic129.0515Prob0.0000
Notes: probability values are listed in the brackets. The values of the probability smaller than the specific confidence level demonstrates the coefficients are significant. a refers to 1% confidence level; b points to 5% confidence level.
Table 8. Results for Granger causality examination.
Table 8. Results for Granger causality examination.
Null HypothesisHebeiHenanInner MongoliaShandongShanxi
LnGDP does not Granger cause lnSO213.90915.676713.84304.15534.9914
(0.0036) a(0.0343) b(0.0037) a(0.0463) b(0.0450) b
LnSO2 does not Granger cause LnGDP2.240720.19021.18680.12430.0970
(0.1770)(0.8309) (0.1037)(0.8851)(0.1987)
LnFEC does not Granger cause lnSO211.65687.06175.069712.80586.3492
(0.0376) b(0.0209) b(0.0464) b(0.0274) b(0.0268) b
LnSO2 does not Granger cause LnFEC5.041496.40675.08364.26325.8372
(0.0496) b(0.0306) b(0.0207) b(0.0159) b(0.0221) b
LnECI does not Granger cause lnSO28.22395.64326.38695.76094.9374
(0.0426) b(0.0347) b(0.0264) b(0.0402) b(0.0214) b
LnSO2 does not Granger cause LnECI2.934882.06072.42650.44821.0361
(0.1187)(0.1978)(0.1583)(0.6559)(0.4035)
LnES does not Granger cause lnSO211.00659.547726.62375.56074.8368
(0.0129) b(0.0113) b(0.0005) a(0.0256) b(0.0252) b
LnSO2 does not Granger cause LnES1.22591.58021.31660.13970.1608
(0.3496)(0.2714)(0.1028)(0.8720)(0.8545)
Notes: probability values are listed in the brackets. a points to 1% confidence level; b refers to 5% confidence level.

Share and Cite

MDPI and ACS Style

Zhao, H.; Guo, S.; Zhao, H. Impacts of GDP, Fossil Fuel Energy Consumption, Energy Consumption Intensity, and Economic Structure on SO2 Emissions: A Multi-Variate Panel Data Model Analysis on Selected Chinese Provinces. Sustainability 2018, 10, 657. https://doi.org/10.3390/su10030657

AMA Style

Zhao H, Guo S, Zhao H. Impacts of GDP, Fossil Fuel Energy Consumption, Energy Consumption Intensity, and Economic Structure on SO2 Emissions: A Multi-Variate Panel Data Model Analysis on Selected Chinese Provinces. Sustainability. 2018; 10(3):657. https://doi.org/10.3390/su10030657

Chicago/Turabian Style

Zhao, Haoran, Sen Guo, and Huiru Zhao. 2018. "Impacts of GDP, Fossil Fuel Energy Consumption, Energy Consumption Intensity, and Economic Structure on SO2 Emissions: A Multi-Variate Panel Data Model Analysis on Selected Chinese Provinces" Sustainability 10, no. 3: 657. https://doi.org/10.3390/su10030657

APA Style

Zhao, H., Guo, S., & Zhao, H. (2018). Impacts of GDP, Fossil Fuel Energy Consumption, Energy Consumption Intensity, and Economic Structure on SO2 Emissions: A Multi-Variate Panel Data Model Analysis on Selected Chinese Provinces. Sustainability, 10(3), 657. https://doi.org/10.3390/su10030657

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop