3.1. Benchmark regression
In this paper, distribution regression is adopted to obtain the DID results with lntq and lnpec as explained variables to analyze the effects of WEGTP on CO2 emissions. Table 2 shows seven models, whose construction is already shown above. Model 1 is derived from Formula (1) when there are only three dummy variables: policy point T, treatment group D, and policy effect, namely, the interaction item . The results in Table (2) show the net impact of WEGTP when no control variable is added. The coefficient of the interaction item is negative, but the significance of the net impact does not meet the standard of 1% to 10%. It indicates that the model without any control variable cannot explain CO2 emissions sufficiently. Therefore, control variables must be added on that basis in order to observe the performance. “Year-end resident population” is added to Model 2. The coefficient of the interaction item is negative just like Model 1, but the 10% level of significance test is failed. After the new control variable “government intervention” is added to Model 3, the 10% level of significance test is passed, and the coefficient of the interaction item is still negative. After the addition of control variables population and government intervention, Model 3 can better explain the explained variable. The coefficients of population and government intervention are both positive and significant, indicating that the increase in population has caused an increase in total CO2 emissions. Each person increased leads to 0.761 units of carbon emissions to meet the needs of industry, transportation and so on in daily life. Each unit GDP of fiscal revenue increase leads to 0.761 units of carbon emissions.
Table 2 PSM-DID estimation of carbon dioxide emissions
Variables
|
Model 1
|
Model 2
|
Model 3
|
Model 4
|
Model 5
|
Model 6
|
Model 7
|
T
|
0.629***
(8.290)
|
0.621***
(13.620)
|
0.551***
(12.160)
|
0.427***
(10.220)
|
0.261***
(5.580)
|
0.253***
(5.420)
|
0.414***
(7.920)
|
D
|
-0.001
(-0.010)
|
0.097*
(1.740)
|
0.108**
(2.020)
|
0.024
(0.490)
|
0.042
(0.900)
|
0.065
(1.400)
|
0.076*
(1.700)
|
|
-0.046
(-0.350)
|
-0.090
(-1.150)
|
-0.126*
(-1.660)
|
-0.134**
(-1.990)
|
-0.115**
(-1.780)
|
-0.120**
(-1.880)
|
-0.111**
(-1.800)
|
LNTP
|
|
0.711***
(28.430)
|
0.761***
(29.950)
|
0.736***
(32.430)
|
0.727***
(33.550)
|
0.808***
(22.690)
|
0.911***
(23.770)
|
GI
|
|
|
4.302***
(6.100)
|
6.099***
(9.430)
|
2.535***
(3.120)
|
2.857***
(3.510)
|
4.757***
(5.620)
|
IS
|
|
|
|
3.025***
(11.030)
|
1.964***
(6.430)
|
1.895***
(6.230)
|
1.529***
(5.110)
|
LNAGDP
|
|
|
|
|
0.234***
(6.740)
|
0.332***
(6.810)
|
0.745***
(8.950)
|
LNRD
|
|
|
|
|
|
-0.075***
-2.840)
|
-0.106***
(-4.070)
|
LNACDI
|
|
|
|
|
|
|
-0.747***
(-6.010)
|
F-test
|
32.780
(0.000)
|
270.730
(0.000)
|
241.380
(0.000)
|
275.180
(0.000)
|
265.550
(0.000)
|
237.020
(0.000)
|
231.190
(0.000)
|
Adjust R2
|
0.172
|
0.703
|
0.725
|
0.783
|
0.802
|
0.805
|
0.819
|
Note: The regression coefficient is t value in parentheses; F - test: Prob > F; "***", "**" and "*" represent the significance level of 1%, 5% and 10% respectively.
IS (industrial structure) is added to Model 4 on the basis of Model 3. The coefficient of the interaction item is the same as the expectation and the negative value is smaller, indicating that CO2 emissions still show a downward trend. As a result, the 5% significance level test is passed. The value of adjusted R2 (coefficient of determination) is 0.783, up by 0.05 points on the basis of Model 3. Moreover, the secondary industry ratio coefficient is 3.025, indicating a very significant positive effect on total CO2 emissions. Due to the different industrial departments, the structure and types of energy consumption are also different. The secondary industry mainly refers to the processing and manufacturing industries. Compared with the primary industry and the tertiary industry, the secondary industry has higher requirements for energy consumption and a higher carbon emission coefficient in production and life. The increased amplitude of CO2 emissions is smaller than that of production output, so CO2 emissions increase.
Average GDP reflects the affluence of the people of a country, and AGDP is added to Model 5 on the basis of Model 4. The coefficient is -0.115, passing the level of significance test. The value of adjusted R2 (coefficient of determination) is 0.802, and the AGDP coefficient becomes 0.234, indicating that AGDP has a positive effect on CO2 emissions. The reason is that when AGDP increases, energy consumption also increases with improved living standards. When the actual energy consumption becomes bigger, and the standard coal conversion coefficient and carbon emission coefficient remain unchanged, the total CO2 emissions become bigger. After the interaction item of “R&D investment” is added to Model 6 on the basis of Model 5, the coefficient becomes -0.120, which once again proves the emissions reduction effect of WEGTP. The improvement of technology level is largely determined by investment in research and development. The improvement of technology level has both positive and negative effects on CO2 emissions: on the one hand, it drives economic growth, thus increasing energy consumption and CO2 emissions; on the other hand, it improves energy utilization efficiency and optimizes energy structure, thus reducing CO2 emissions. According to the table, when the value of adjusted R2 (coefficient of determination) is adjusted to 0.805, the coefficient of R & D investment is -0.075, indicating a negative effect of R & D investment on CO2 emissions. This result is consistent with the theoretical possibility of reducing CO2 emissions. Increased R & D investment promotes technological progress. In the past 10 to 20 years, China's scientific research units, taking into account market demand, have increased research and development in ECER and achieved considerable achievements. People use high-level production tools in daily life, or industrial production, or technological improvements lead to the optimization of energy structure, thus reducing CO2 emissions.
Model 7 puts together six control variables, with per capita disposable income newly added on the basis of Model 6 The value of the interaction item becomes -0.111, and the result is significant. Throughout the above process, the regression coefficient of PSM-DID remains negative, indicating that the conclusion “WEGTP reduces CO2 emissions” is stable after adding the control variables. From the table, we can see that per capita disposable income reduces CO2 emissions, possibly because the increase of ACDI enhances purchasing power, so that people have a more solid material basis to improve consumption structure-- they become more willing to use clean energy like natural gas, thus increasing the demand for clean energy and contributing to carbon emissions reduction.
3.2 Heterogeneity test
For energy importer provinces and energy exporter provinces, the impacts of WEGTP on energy consumption and carbon emissions are different. In energy importer provinces, the supply of natural gas increases total energy consumption, but changes in the total quantity and intensity of carbon emissions are uncertain. In energy exporter provinces, natural gas takes a substitution effect, the primary energy consumption structure sees a substantial reduction in the proportion of coal, and there should be a significant decline in carbon emissions. According to the above characteristics, the paper will next conduct a heterogeneity test on the energy importer provinces and the energy exporter provinces in the treatment group.
Table 3 Estimates of energy consumption in import provinces
Variables
|
Energy consumption in imported provinces
|
Model 8
|
Model 9
|
Model 10
|
Model 11
|
|
-0.071***
(-3.150)
|
-0.077***
(-3.510)
|
-0.066***
(-2.960)
|
-0.065***
(-2.860)
|
LNTP
|
|
0.423***
(3.830)
|
0.478***
(4.250)
|
0.489***
(4.260)
|
GI
|
|
|
-1.351***
(-2.220)
|
-1.324***
(-2.170)
|
IS
|
|
|
|
0.110
(0.500)
|
LNAGDP
|
0.319***
(4.920)
|
0.459***
(6.280)
|
0.468***
(6.480)
|
0.448***
(5.380)
|
F-test
|
378.57
(0.000)
|
375.800
(0.000)
|
361.310
(0.000)
|
342.330
(0.000)
|
Adjust R2
|
0.958
|
0.960
|
0.961
|
0.961
|
Note: The regression coefficient is t value in parentheses; F - test: Prob > F; "***", "**" and "*" represent the significance level of 1%, 5% and 10% respectively.
This paper still adopts the abovementioned method of substituting control variables step by step. According to the regression results of four models on the energy consumption of importer provinces, as shown in Table 3, the interaction items of the four models are all negative at a significance level under 10%. The coefficient of interaction item of Model 8 is negative, indicating that WEGTP brings a downward trend to energy consumption in the importer provinces. However, the AGDP coefficient is positive, probably because the importer provinces are mostly located in the eastern coastal areas which were densely populated. With the increase of per capita income of residents, the narrowing gap between urban areas and rural areas and the realization of common prosperity, and the overall balanced and coordinated local development, the energy consumption levels of urban and rural residents began to rise, gradually raising the total energy consumption of local society as a whole. Year-end resident population, government intervention, and industrial structure are added to Model 9, Model 10 and Model 4 respectively. It can be seen that the interaction item result and control variables in Model 9 and Model 10 are significant, indicating that energy consumption was reduced in importer provinces after the opening of WEGTP. Although the interaction item is negative and significant after the addition of industrial structure to Model 11, the impact of industrial structure alone is not significant. The reason may be that although the proportion of the secondary industry in GDP is gradually declining, the industrial sector itself is a large energy consumer, and that WEGTP has a structural substitution effect on industrial energy consumption, however the energy conservation effect is subject to technical conditions. The manufacturing sectors of the eastern importer provinces undertake a large amount of export commodity production, and there is a trend that high energy consumption commodities are produced domestically and then exported to foreign countries.
Table 4 Heterogeneity analysis of energy export provinces
Variables
|
Emission of CO2(lnTQ)
|
Energy consumption(lnPEC)
|
|
-0.154**
(-1.820)
|
-0.111*
(-1.740)
|
LNTP
|
0.884***
(19.100)
|
0.964***
(27.550)
|
GI
|
3.068***
(2.840)
|
3.906***
(4.780)
|
IS
|
1.866***
(4.740)
|
1.677***
(5.630)
|
LNAGDP
|
0.566***
(5.720)
|
0.888***
(11.86)
|
LNACDI
|
-0.437***
(-2.570)
|
-0.755***
(-5.880)
|
F-test
|
161.340
(0.000)
|
329.400
(0.000)
|
Adjust R2
|
0.797
|
0.889
|
Note: The regression coefficient is t value in parentheses; F - test: Prob > F; "***", "**" and "*" represent the significance level of 1%, 5% and 10% respectively.
This paper discusses the heterogeneity of energy exporter provinces from the two explained variables of CO2 emission and energy consumption. It can be seen from Table 4 that when CO2 emission is used as an explained variable of energy exporter provinces, the coefficient of interaction item is -0.154. The significance test is passed, and the adjusted R2 is 0.797, indicating that CO2 emissions are reduced in energy exporter provinces under. The reason may be that China's primary energy structure is dominated by coal. After the comprehensive operation of WEGTP, coal became uneconomical in terms of price and thermal efficiency, while natural gas is superior in price, thermal efficiency and carbon emissions. Coal consumption in the market, subject to the government's ECER targets stipulated by the 11th Five-Year Plan and the 12th Five-Year Plan, was first replaced. In the same period, many northern provinces gradually replaced coal with natural gas to fuel urban central heating boilers, which also directly reduced CO2 emissions in winter. When energy consumption is used as the explained variable of energy exporter provinces, the coefficient of interaction item is still negative at the significance level, and the adjusted R2 is 0.889, indicating that the energy consumption in exporter provinces was declining, which is consistent with the CO2 emission results obtained above.
3.3 Robustness test
In order to prove the credibility of the above PSM-DID regression results, this paper provides a robustness test on the above regression results through a randomness test of policy intervention time, a test of the control group’s immunity to policy, and a test of policy implementation uniqueness listed in Table 5.
1. Randomness test of policy intervention time
First, a randomness test of policy intervention time was conducted. The settings of the treatment group and the control group kept unchanged, the time node of policy intervention was advanced to 2001. It can be seen from table 5 that the interaction item result is not significant, indicating that CO2 emissions did not decrease significantly in 2001. Similarly, the time node was advanced to 2000. Then, the interaction item result is still not significant. To sum up, the policy intervention time is random.
2. Placebo test of the policy treatment
By supposing a place that is not affected by the policy as a placebo sample, to estimate the treatment of hypothetical policy as if WEGTP accessed in there. That is, when the time node of policy intervention is kept unchanged, provinces are randomly selected from the original control group as the treatment group to observe whether the coefficient of interaction item is significant. If it is significant, it indicates that the control group is affected by policy, and the test is failed; if not, it indicates that the control group is not affected by policy, and the test is passed. It can be seen from Table 5 that the coefficient of interaction item is not significant, so the selection of the control group passed the test.
3. Test of policy implementation uniqueness
WEGTP was opened roughly at the end of 2004. In the same period, in addition to the policy of WEGTP, the Great Western Development Strategy was also underway. It is one-sided to ignore the impact of the latter policy and simply attribute CO2 emissions reduction in that period to WEGTP. Therefore, this paper eliminates Xinjiang, Gansu, Qinghai, Yunnan and other provinces involved in the Great Western Development Strategy, and re-conducts PSM-DID estimation on the processed samples for the purpose of eliminating the impact of the Great Western Development Strategy on the explained variable. From the new regression results, we find the coefficient of interaction item is -0.097 and significant, which is only a gap of 0.014 from the previous result. This shows that the Great Western Development Strategy has little impact on CO2 emissions, and WEGTP is the main policy reason affecting CO2 emissions.
3.4 Discussion on mechanisms
The mechanisms of action of WEGTP on China's energy consumption and CO2 emissions reduction are as follows:
First, the use of natural gas increased greatly. The natural gas consumption in the involved provinces structurally exceeded that in the non-involved provinces. Figure 1 shows that the average consumption of the involved provinces exceeds that of the non-involved provinces, solving the energy shortage problem in the economic development process of central and eastern regions.
Second, an interprovincial structural adjustment mechanism was generated. The natural gas substitution effect and the overuse of natural gas under price control occurred in the importer provinces. Natural gas has the advantages of high thermal efficiency, low price and low emission intensity, encouraging the enterprises in central and eastern energy importer provinces to change their industrial energy from coal to natural gas, thus bringing about a continuous rise in the total energy consumption in the importer provinces.
Third, the low emission property of natural gas is conducive to the reduction of CO2 emissions in central and eastern importer provinces, but due to the rapid economic growth in the eastern region, the effect of total energy consumption increase and the CO2 emissions reduction effect brought by energy structure changes neutralized each other, so that the reduction effect of total CO2 emissions was not obvious.
Fourth, the energy exporter provinces of WEGTP show a trend of increased ECER pressure and increasingly insignificant energy consumption. Because natural gas has a more significant substitution effect on high pollution energy, under the ECER policy constraints of local governments, reducing CO2 emission intensity brought by reducing coal consumption should be more obvious.
Table 5 Robustness test
Variables
|
Randomness of policy intervention time
|
The spillover of policy effects
|
Uniqueness of policy implementation
|
|
T1=2001
|
T2=2000
|
|
GWDS policy
|
|
-0.021
(-0.330)
|
0.007
(0.170)
|
0.043
(1.02)
|
-0.097***
(-2.940)
|
F-test
|
98.330
(0.000)
|
121.730
(0.000)
|
122.080
(0.000)
|
109.570
(0.000)
|
Control variables
|
Controled
|
Controled
|
Controled
|
Controled
|
Adjust R2
|
0.873
|
0.868
|
0.868
|
0.903
|
N
|
459
|
459
|
459
|
300
|
Note: The regression coefficient is t value in parentheses; F - test: Prob > F; "***", "**" and "*" represent the significance level of 1%, 5% and 10% respectively.