Spatiotemporal Characteristics, Determinants and Scenario Analysis of CO2 Emissions in China Using Provincial Panel Data

This paper empirically investigated the spatiotemporal variations, influencing factors and future emission trends of China’s CO2 emissions based on a provincial panel data set. A series of panel econometric models were used taking the period 1995–2011 into consideration. The results indicated that CO2 emissions in China increased over time, and were characterized by noticeable regional discrepancies; in addition, CO2 emissions also exhibited properties of spatial dependence and convergence. Factors such as population scale, economic level and urbanization level exerted a positive influence on CO2 emissions. Conversely, energy intensity was identified as having a negative influence on CO2 emissions. In addition, the significance of the relationship between CO2 emissions and the four variables varied across the provinces based on their scale of economic development. Scenario simulations further showed that the scenario of middle economic growth, middle population increase, low urbanization growth, and high technology improvement (here referred to as Scenario BTU), constitutes the best development model for China to realize the future sustainable development. Based on these empirical findings, we also provide a number of policy recommendations with respect to the future mitigation of CO2 emissions.


Introduction
Global climate change-specifically, global warming-constitutes one of the most important issues to face human beings in the 21 st century [1]. That greenhouse gases, most notably carbon dioxide (CO 2 ) emissions from the combustion of fossil fuel, are the main influencing factor of global warming which is now a matter of global consensus and, since CO 2 emissions are closely related to socioeconomic development, climate change has gone from being an issue for pure The third category is constituted by studies that have relied on panel data analysis. Typically, this method has been used to examine the influencing factors of CO 2 emissions based on panel data. It is recognized that panel data models have many metrics over conventional cross-sectional and time series data models [34]. Recently, using the method of panel data analysis, Dinda and Coondoo [35] and Narayan and Narayan [36] analysed the relationship between economic growth and CO 2 emissions, and Wang et al. [37], Ozkan and Ozkan [38], Zhu et al. [39] and Ozcan [40] have also attempted comprehensive analyses of the impact factors of CO 2 emissions, finding population and economic level to be the main factors affecting those emissions.
The fourth category of studies is made up of research using a range of other determinant analysis models. The representative analysis methods in this group include: the input-output model [9,41], the Laspeyres method [42], the AWD method [43], the GFI method [44], the ARDL model [45][46], the LR-MEF (long-run marginal emissions factor) method [47] and the Decoupling index formulation [48]. These methods have provided new insights into the influence mechanisms in relation to CO 2 emissions, and have broadened our methodological perspectives.
These existing studies have enriched our understanding of the main human impact factors in CO 2 emissions. However, most of the literature has been inclined to analyse influence mechanisms at a cross-country level [34], at a single-city level or at a national level [18,20], and only a handful of studies have addressed the cross-province level of China. In addition, despite the fact that it is recognized that panel data sets have many metrics over conventional cross-sectional or time series data sets [33][34]49], most models used in simulating China's CO 2 emissions have been based on country or city level time series or cross-sectional data [33], with only a few studies being based on panel data models [2,33]. More importantly, existing studies have mainly focused on estimating the influencing factors of CO 2 emissions, seldom proceeding to analyse the development models that might promote a country's sustainable development.
Based on the above review, this study first calculated CO 2 emissions for China's 30 provinces over the period 1995-2011, and then analysing the spatiotemporal characteristics. We proceed to examine the key influencing factors quantitatively using panel data models. Based on the coefficients calculated in this study, and taking China as an example, this paper sets up 10 different scenarios in order to forecast future CO 2 emissions. Under a CO 2 emissions constraint force, we attempt to identify the optimal development model for China.
The remainder of this study is organised as follows. Section 2 focuses on methods and data, presenting the models and the data utilized. Section 3 displays the results of the paper and sets out a discussion of the scenarios through which the optimal development model was selected. Finally, Section 4 summaries the main conclusions and policy implications.

Methodology and Data
Estimating CO 2

emissions
In this study, CO 2 emissions are described in terms of CO 2 equivalent emissions, a measure which refers to the carbon content of greenhouse gases with the same global warming potential. Currently, in China it is very difficult to acquire provincial data for CO 2 emissions directly. However, because CO 2 emissions are mainly released through fossil fuel combustion and cement production, these two measures prove useful in approximating CO 2 emissions through energy-related statistical data. Based on methods developed in previous studies [33], this paper calculates CO 2 emissions on the basis of these two aspects (that is, the burning of fossil fuel and the cement production process) for China's 30 provinces from 1995 to 2011. Using the CO 2 emissions coefficients published by the 2006 the Intergovernmental Panel on Climate Change (IPCC) [50] and the National Coordination Committee Office on Climate Change and the Energy Research Institute under the National Development and Reform Commission [51], the CO 2 emissions from fossil fuel can be calculated through the following formula: where CE denotes the total CO 2 emissions from energy consumption; i represents the types of fossil fuel; PEC i denotes the consumption quantity of fossil fuel i; and F i is the CO 2 emissions coefficient of fossil fuels i (Table 1). CO 2 emissions from cement production can be calculated using: where CC denotes the CO 2 emissions from the production of cement; Q represents the quantity of cement production; and F is the CO 2 emissions coefficient of the cement production process (Table 1). Panel data analysis. Typically, panel data has two dimensions-time and cross-sections. However, compared to the cross-sectional or time series data models, panel data models have several advantages. Firstly, they can expand the amount of information and increase the degree of freedom of estimation and test statistics. Secondly, they can improve the effectiveness of the estimation. Thirdly, effects can be more efficiently recognised and measured than through other models. Fourth, panel data models help to provide the reliability of dynamic analysis. Fifth, they help to detect gradual change in a data structure and institution.
To explore the constraint effect of CO 2 emissions on the economic development process, a CO 2 emissions model was developed. Based on existing studies, four variables-namely, population, per capita GDP, energy intensity, and the urbanization level-were utilised as the determinants of CO 2 emissions. The CO 2 emissions model is established as follows: where ln is a natural logarithm; EM denotes CO 2 emissions; P represents the resident population; A denotes per capita GDP; T represents technological level (energy intensity); U is the urbanization level; β 1 , β 2 , β 3 and β 3 are slope coefficients; t denotes time; and i denotes the cross-section (provinces). Panel unit root tests. It is well known that panel unit root tests have higher ability than unit root tests that are based on univariate time series or cross section data. This paper will introduce three types of panel unit root test-namely, the Levin-Lin-Chu (LLC) test [52], the Breitung test [53] and the Im-Pesaran-Shin (IPS) test [54]. These tests are commonly utilised because of their high power in stationary tests among variables. The LLC and Breitung tests are widely used in common root tests. The IPS test is widely used in individual root tests.
Panel cointegration tests. Compared to normal time series cointegration, panel cointegration tests are widely used among researchers because of their high power [34]. As such, a Pedroni cointegration test, first developed by Pedroni [55], was introduced into the study to Working under a series of assumptions, Pedroni [55] develops two types of residual-based tests. In the first type, four tests are distributed as being standard normal asymptotically and are based on average test statistics for no cointegration in the time series across the cross-section for the within-group [56]. The tests that make up this group are: panel v-statistics, panel ρ-statistics, panel t-Statistic (non-parametric) and panel t-Statistic (parametric). Within the second type, three tests-group r-statistics, group t-statistics (non-parametric) and the group statistics (parametric)-are also distributed as being standard normal asymptotically but are based on limits of piecewise numerator and denominator terms for the limiting distributions for the between-group [56]. Three of the four tests in the first type of test are based on the nonparametric correction test proposed by Phillips and Perron [57]; the remaining test follows the ADF parametric test. Two of the three second-type tests are based on non-parametric tests, and the third test is based on an ADF test.
Panel data model. If variables are cointegrated, the Fully Modified OLS (FMOLS) estimator can be utilised to estimate each equation to find out the directions and magnitudes of the influences of the explanatory variables (in this study, population scale, economic level, energy intensity and urbanization level) on CO 2 emissions. FMOLS, proposed by Pedroni [58], is a single equation estimator for cointegrated relationships. FMOLS uses errors to calculate the cumulative test volume. This approach has been widely used to deal with both temporal and cross sectional dimensions data, as a result of its ability to properly deal with the long-run correlation problem between the cointegrating equation and stochastic regressor innovations [29].

Scenario simulations
Previous studies have found population, economic level, urbanization and energy intensity to have a significant impact on CO 2 emissions. Using scientific analysis and a range of different scenarios, relevant existing literature has thus also attempted to identify which scenario might be most conducive to the mitigation of CO 2 emissions. The identification of such a scenario will have great significance and meaning for the formulation of future economic development policy and the mitigation of CO 2 emissions. Taking China as an example, in this study we attempted to find the optimal development model for China, from 10 different scenarios which we developed. We assumed that the 10 different scenarios share various evolution paths, so that we might focus on analysing differences in the various development patterns.
Scenario description. (1) Economic growth scenario: Table 2 reviews the 3 economic growth scenarios. The "business-as-usual" scenario (B) assumes that the development policy on economic growth stays at the same level as in the Chinese government's 12th Five-Year Plan (2011-2015) within the forecasting horizon. According to the development experience of developed countries, after rapid economic growth, subsequent economic growth is likely to flatten. Thus, the economic growth rate of the "high economic growth" scenario and the "low economic growth" scenario increase and decrease by 1% respectively based on the economic growth rate in the business-as-usual scenario.
(2) Population scenario: According to the U.S. Energy Information Administration's (EIA) [59] forecast and the Chinese government's 12th Five-Year Plan (2011-2015), the population of China will be controlled to 1.39 billion by 2015 and 1.42 billion by 2020 respectively. Given the above, the average population growth rate will be 0.5% from 2016 to 2020. According to this rate, we therefore can predict that the population of China in 2025 will be 1.45 billion. In the "high population" scenario, the population growth rate increases by 0.1%. Table 3 reviews the predicted population of China in the following 10 years.
(3) Urbanization scenario: China has made great success in urbanization, with the urbanization level increasing from 17.9% in 1978 to 51.27% in 2011, at an average annual growth rate of 1.01% [1]. According to the National New-type Urbanization Plan for China released in 2014, the country's urbanization level will rise by 1% annually, reaching 60% by 2020. Taking the Chinese government's 12th Five-Year Plan (2011-2015) into consideration, we assumed that the urbanization level in 2015, 2020, and 2025 would be 55%, 60%, and 65% respectively. The urbanization level in the "low urbanization" scenario will decrease by 2% correspondingly. Table 4 reviews predicted levels of urbanization of China for the following 10 years.
(4) Technological progress scenario: In the content of the 11th Five-Year Plan (2005-2010), China set its goal of reducing energy intensity by 20%. Through the 12th Five-Year Plan (2011-2015), the government plans to make a decrease of 16% in energy intensity [1]. Specifically, energy intensity is to decrease at an average annual rate of 3%. Thus we referred to the targets set out in these plans as constituting the basis for our "middle technology" scenario and assumed an additional 0.5% improvement under the "high technology" scenario. Table 5 reviews China's predicted technological progress in the following 10 years.
Based on above analysis, together with the actual development of China in recent years, we set up 10 different scenarios ( Table 6). As setting up the 10 scenarios, we consider not only the permutations and combinations of the variables, but also the actual situation of China's development. Since the goal of this study is to find out the optimal development model for china's socioeconomic sustainable development and a reduction in CO 2 emissions. The requires of the best model are not only the lowest CO 2 emissions, but maintaining relatively high economic growth. Thus, the four variables (economic growth, population growth, urbanization development, and technological progress) combined 10 development scenarios (Scenario B, Scenario BTU, Scenario BH, Scenario L, Scenario LT, Scenario LH, Scenario H, Scenario HP, Scenario HT, and Scenario AH). The detailed content is listed in Table 6.
Methods for forecasting future CO 2 emissions. According to the results of the panel data analysis, the following regression model will be utilized to forecast CO 2 emissions for each province under 10 different scenarios for the year 2015 to 2025: where EM denotes the total CO 2 emissions both form the burning of fossil fuel and the cement production process. The following exponential transformation will be further implemented since we aim to forecast EM it : Once each province's CO 2 emissions are estimated, we then calculate the emissions of China using: Data Based on the availability of data, this study used annual time series data for 30 provinces in China, in the period 1995-2011. Taiwan and Tibet were excluded on the basis of missing data for most of the years. Provincial data on resident populations, per capita GDP, cement, energy intensity and urbanization levels were derived from the China Statistical Yearbook and provincial statistical yearbooks. To eliminate the price effect, per capita GDP data were adjusted by considering the official price index in the year 2000, which was also used to calculate energy intensity in tons per 10 4 Yuan. All the data on fossil fuel were derived from the China Energy Statistical Yearbook. The total energy consumption and each fossil fuel were all converted into standard coal measures (units of 10 4 tons). Annual data on CO 2 emissions for each province were estimated using (Eq 1) and (Eq 2). Table 7

Results and Discussion
Spatiotemporal characteristics of provincial CO 2 emissions Once the emissions of each province were estimated from 1995 to 2011, it is necessary to analyse the spatiotemporal characteristics in emissions. The former province produces more CO 2 emissions owing to its huge consumption of fossil energy compared to other regions. Notably, CO 2 emissions form east regions were initially larger and increased more rapidly than their neighbours in central and west regions during the period studied. From the viewpoint of geography, CO 2 emissions in China exerted noticeable regional discrepancies in 1995-2011 (Fig 3). Provinces with high CO 2 emissions were found to be mainly concentrated in eastern coastal China; in fact, CO 2 emissions decrease gradually from the eastern coastal region to the central region, and then to the western region. We then calculated global Moran's I as an index to detect spatial dependence. Fig 4 plots the trajectories of the spatial dependence of CO 2 emissions in provincial China from 1995 to 2011. As indicated in Fig 4, Moran's I fluctuates around 0.2, all are significant at 95% confidence level over the period studied. This finding provided evidence of spatial dependence taking place in provincial China. In addition, the fluctuations of Moran's I reflect the temporal dynamics of the spatial dependence of CO 2 emissions, which may weaken or strengthen the agglomeration of CO 2 emissions over time. From the viewpoint of scatter distribution (Fig 5), HH and LL clusters constitute the main types of spatial dependence seen over the studied period.  The negative slope of the linear regression equation implies a convergence trend in CO 2 emissions in Chinese provinces.

Influencing factors of provincial CO 2 emissions
Results of panel unit root tests. The analysis of casual relationships between variables has been the hot topic in modern panel data analysis. However, not all variables have cointegrating relationships. Thus, prior to cointegration estimation, it is necessary to test the stationarity of variables; this is known as a unit root test. In this study, we employed three types of panel unit root test: the Levin-Lin-Chu (LLC) test, the Breitung test and the Im-Pesaran-Shin (IPS) test. The LLC and Breitung methods are both suitable for common root tests, and the IPS method is suitable for individual root tests. Table 8 reviews the results of the panel unit root test for both types. As indicated by Table 8, the results clearly show that all the variables cannot reject the null hypothesis at level, indicating that they were not stationary at level. However, all the variables were found to be stationary at the first difference, rejecting the null hypothesis. As such, we were able to carry out a Pedroni cointegration test in order to estimate whether a cointegrating relationship existed between the variables. If the variables were found to be cointegrated, the Fully Modified OLS estimator would be used to examine the relationship quantitatively.
Panel cointegration results. Since all the variables used in this study were stationary at their first difference, we proceeded to test the dependent and independent variables respectively for Pedroni cointegration in the data, using the heterogeneous panel cointegration test first proposed by Pedroni [55]. Table 9 presents the panel cointegration test results. To examine the cointegrating relationship more accurately, we employed three trend assumptions from the Pedroni cointegration test, namely: no deterministic trend, deterministic intercept and trend, and no deterministic intercept or trend. From the results detailed in Table 9, we found that under the assumption of no deterministic trend, six statistics strongly rejected the null hypothesis of no-cointegration significantly at the 1% level; while under the assumption of deterministic intercept and trend, seven statistics were found to reject the null hypothesis of no cointegration at less than the 10% level. In addition, under the assumption of no deterministic intercept or trend, six statistics also strongly rejected the null hypothesis of no cointegration significantly at less than the 5% level. Thus we obtained clear evidence of the existence of a cointegrating relationship between the variables. Results also indicate that population scale, economic level, energy intensity and urbanization level had a long-run relationship with CO 2 emissions in China's 30 provinces during the study period. The next step was to estimate this relationship using Fully Modified OLS regression.
Parameter estimations of the panel model. Since the variables used in this study were cointegrated, we proceeded to estimate the relation. A Fully Modified OLS test was utilised based on provincial panel data. A panel data model with a Newey-West fixed effect was adopted; the Kernel option of the model was based on the Bartlett law.  Characteristics, Determinants and Simulations of CO 2 Emissions individual variable coefficients have important but various influence on CO 2 emissions. Specifically, population scale, economic level and urbanization level have positive effects on CO 2 emissions, while energy intensity produces negative, inhibitory effects. First, as expected, the increase of population directly leads to an increasing demand of energy resources and, in turn, produce a large amount of CO 2 emissions. Second, since the Chinese economic reform in 1978, China has witnessed rapid economic growth. However, this rapid growth has been achieved by consuming a large amount of energy resources, leading to a great deal of CO 2 emissions. Third, millions of people lived in rural areas moved to cities and towns across China every year. This migration creates increasing pressure on the country's urban sustainable development. Urbanization has a promoting effect on economic development, but it can also contribute to the increase of energy consumption. Finally, reducing energy intensity (per capita GDP energy consumption) is  conducive to reduce CO 2 emissions. That is because improving energy efficiency mean producing the same amount of GDP with mitigated energy use, which will indirectly migitate CO 2 emissions [18]. In addition, the significance of the relationship between CO 2 emissions and the four variables also varied across the provinces based on their scale of social and economic development. Similar results were found by Al-mulali [34] in the Middle East and Li et al. [56] at regional and national levels in China. From the perspective of geographic analysis (Fig 7),  Characteristics, Determinants and Simulations of CO 2 Emissions provinces with high coefficients of population and energy intensity were found to be mainly concentrated in northern and central China respectively. However, areas with relative high coefficients of the other two variables (per capita GDP and urbanization) were found to have decentralized distributions. Since this study used a panel data model, we will focus on the panel effect and panel analysis results in this paper. The panel FMOLS test results show that all the four variables-population, per capita GDP, energy intensity and urbanization-all maintained a positive long-run relationship with CO 2 emissions-that is: a 1% increase in total population will increase CO 2 emissions by 1.359500% and a 1% increase in per capita GDP will increase CO 2 emissions by 0.373704%, and, in addition, a 1% increase in energy intensity will increase CO 2 emission by 0.277178%. Finally, a 1% increase in urbanization will increase CO 2 emission by 0.591743%. The results of the study indicate that population, per capita GDP, energy intensity and urbanization are the main human factors that influence CO 2 emissions in the long run. These results are supported by studies undertaken by Wang et al. [18] in Beijing city, Wang et al. [20] Guangdong province, Al-mulali [34] in Middle East, Soytas et al. [60] in the United States, Hamit-Haggar [61] in Canada, Niu et al. [62] in 8 Asian-Pacific countries, Pao and Tsai [63] in BRIC countries, and Kum et al. [64] in G-7 countries.  The findings of this study indicate that in a scenario where the country's economic growth rate decreases by 1% and the process of urbanization concurrently slows (Scenario L), China's CO 2 emissions will decrease by 0.39, 0.65, and 0.97 billion tons by 2015, 2020, and 2025 respectively, in comparison with Scenario B. If the economic growth rate decreases, however, and low-carbon technologies have been promoted (Scenario LT), China's CO 2 emissions will decline by 0.41, 0.77, and 1.20 billion tons by 2015, 2020, 2025 respectively, compared to Scenario B. At the same time, of all the scenarios explored in this study, we find China to have the lowest CO 2 emissions under Scenario LT. If the country's population continues to increase on the basis of Scenario LT, CO 2 emissions will increase accordingly. If the economy has a middle level of growth, urbanization is relatively low, low-carbon technologies rapidly improve, and population growth occurs at a middle level (Scenario BTU), China's CO 2 emissions will The above analysis indicates that China would have the lowest CO 2 emissions under Scenario LT. However, Scenario LT requires that the economy maintain a relatively low level of growth, that population growth maintain a middle level, that urbanization increase at a low level, and that low-carbon technologies be promoted widely. The rapid development model currently pursued by China does not conform to these conditions. As the world's largest CO 2 emitter, China faces the dual challenge of mitigating CO 2 emissions while concurrently maintaining economic growth. Therefore, among the 10 scenarios, only the Scenario BTU conforms to the aim of maintaining a middle level of economic growth and effectively reducing CO 2 emissions.

Forecasting CO 2 emissions
It is essential to compare the results of our estimations of CO 2 emissions with other recent studies. Table 11 reviews the results of estimations made by Chen et al. [65], IEA [66], Cai et al. [67], EIA [68], ERI [69], Du et al. [33] and this study. The estimated emissions detailed in this study are higher than that of IEA [66], Cai et al. [67], and EIA [68], but are close to that of Du et al. [33]. It is notable that two sources were considered in estimating CO 2 emissions in this study: the combustion of fossil fuel and the cement production process. It is noted that the studies identified above only consider fossil fuel related emissions. Du et al. [33], however, accounted for emissions from both the sources, but set scenarios and selected impact factors which varied slightly from those used in the present study. Despite these differences, we note that our forecasting result is quite close to the result obtained by Du et al. [33].

Conclusions and Policy Implications
China's rapid economic growth has brought about increases in the scale of carbon dioxide (CO 2 ) emissions produced by the country. A growing question lies in the optimal development  and taking China as an example, this paper created 10 different scenarios to forecast future CO 2 emissions. Under a CO 2 emissions constraint force, we attempted to find the optimal development model for China. There were significant differences of CO 2 emissions among China's provinces. CO 2 emissions of all provinces have increased over the studied period and the emissions varied across the provinces. In addition, evidence of spatial dependence was identified in provincial CO 2 emissions, and the fluctuations of Moran's I may have important role in the agglomeration of CO 2 emissions over time. A graphical analysis of the initial level and average annual growth rate of CO 2 emissions suggested a convergence existing in provincial China.  Parameter estimations of the panel data model revealed that population scale, economic level and urbanization level positively influenced CO 2 emissions, while energy intensity exerted negative, inhibitory effects. In addition, the significance of the relationship between CO 2 emissions and the four variables varied across the provinces based on their scale of social and economic development. From the perspective of panel data analysis, and focusing on panel effect, a quantitative relationship between CO 2 emissions and the four independent variables was identified: a 1% increase in population, per capita GDP, energy intensity and urbanization was found to result in an increase in CO 2 emissions by 1.359500%, 0.373704%, 0.277178%, and 0.591743% respectively.
To achieve the goal of curbing CO 2 emissions while maintaining economic growth, this paper reviewed 10 different developing scenarios. The study forecasted future CO 2 emissions based on the relationship obtained from the panel data analysis. Under a CO 2 emissions constraint force, scenario simulations further showed that Scenario BTU-that is, middle economic growth, middle population increase, low urbanization growth, and high technology improvement-would provide the best development model for China to realise socioeconomic sustainable development and a reduction in CO 2 emissions.
The above findings thus contribute to the literature and suggest meaningful theoretical and policy implications [70][71][72]. First, from the evidence provided by the results of this study, population scale and urbanization are important factors that increase CO 2 emissions. Since the Chinese economic reform in 1978, China has entered a period of rapid urbanization. Results suggest that China should continue to carry out the family planning policy, which has aimed to slow the robust growth of the country's vast population and maintain healthy levels of population urbanization-according to official data, this policy has helped Asia's biggest economy reduce its population by millions. In addition, China should also continue to promote stable and relative low population urbanization and pay attention to optimizing population structure and quality. More importantly, it is essential for China to make greater efforts to improve lowcarbon awareness, advocating for low-carbon consumption and healthy lifestyles and promoting sustainable consumption modes to households. Second, the results of this study clearly indicate that economic level positively influences CO 2 emissions. Over the past thirty years, China has made great success in economic development, with an annual growth rate of 9.9%. Rapid economic growth can improve living standards, but it can also cause high levels of energy consumption, leading to high CO 2 emissions. Rapid economic growth, however, retains a high strategic importance for the Chinese government. The government therefore aims to reduce CO 2 emissions, while still fostering economic growth. The most feasible method to reduce CO 2 emissions for the Chinese government is not, as such, to sacrifice future economic growth; rather, slowing economic growth relatively constitutes the most effective measure. Third, China pledged to reduce its energy intensity by 20% on average across all provinces by 2010. In the future, China should continue to target cuts in energy intensity. First, emissions reduction indicators, such as energy and emission intensity, should be set relative to physical output rather in terms of economic growth [4]. In addition, it is also important to take various measures such as producing more renewable technologies and low-carbon technologies, boosting recycling and renewable energies, and implementing incentive policies, amongst other measures, in order to strengthen energy conservation and emission mitigation [73].
From the view of methodology, this paper underscores the promising aspects of utilizing panel data models such as panel unit root models, panel cointegration tests, and the Fully Modified OLS in understanding the main determinants of CO 2 emissions. The results of the panel data models are capable of better revealing the factors hidden in the CO 2 emissions in Chinese provinces over the period studied. Our empirical analysis of Chinese provinces also demonstrates the appropriateness of setting up different scenario simulations for analyzing the optimal development model for China by addressing the influencing variables' dynamic evolution changes.
Supporting Information S1