Parallel trend test
An essential prerequisite for the validity of the DID model is that the control and treatment groups should satisfy the parallel trend to ensure unbiased estimation. Since the LCCP policy was implemented in three batches, rather than one single implementation, the grouping of a city (control or treatment) could change. Therefore, this study applied the event study method, instead of the plotting trend method, to detect parallel trends more precisely (Moser and Voena,2012). The time window was set to 17 years, covering the 7 years that preceded the implementation of the LCCP and the 9 years that followed it. Based on the study of Yu and Zhang (2021), the following equation was used for the event study method:
$${y_{it}}={\alpha _0}+\sum\limits_{{ - 7}}^{9} {{\alpha _k} \times {D_{i,{t_0}+k}}} +\beta {X_{it}}+{\eta _i}+{\gamma _t}+{\varepsilon _{it}}$$
2
Where Di, t0+k represents a series of dummy variables associated with the years of implementation of the LCCP policy, t0 represents the first year of the LCCP policy implementation, and k denotes the kth year of the start of the LCCP policy. The other variables have the same meaning as in Eq. (1). Parameter αk reflects dynamic effects of the LCCP policy on CO2 emission reduction. α-7-α‐1 test the parallel trend assumption, i.e., if the hypothesis αk = 0 cannot be rejected, implying that there is no difference in CI and CP between the treatment and control groups before the implementation of the LCCP policy.
Figure 2 illustrates the test coefficients of the 95% confidence interval of Di, t0+k. The horizontal axis represents the year before and after implementation of the LCCP policy, the vertical axis indicates the difference in change of the two dependent variables. According to the coefficients in the pre-treatment period, we can deduce that the CI and CP of treatment and control groups would follow a similar trend without the LCCP policy. Hence, the parallel trend assumption could not be rejected. Figure 2 also shows that after the implementation of the LCCP policy, the carbon emission reduction effect started to be significant in the third year, supporting Hypothesis 1.
Baseline results
Table 2 demonstrates the estimation of the benchmark regression results. Columns (1) and (2) show the average impact of the LCCP policy on CI. Columns (3) and (4) show the average impact of the LCCP policy on CP. The fixed effects of city and year were controlled in columns (1)-(4). According to the results of columns (1) and (3), the LCCP policy not only significantly decreased CI but also significantly reduced CP, with regression coefficients of -0.719 and − 0.199, respectively. That is, the implementation of the LCCP policy simultaneously achieved carbon emission reduction effect from both production and consumption sides. For robustness, columns (2) and (4) contained the control variables. The regression results still showed a significance level of 1%, which further validates Hypothesis 1.
Table 2
Benchmark regression results
|
CI
|
CI
|
CP
|
CP
|
treat
|
-0.719***
|
-0.579***
|
-0.199***
|
-0.181***
|
|
(-6.212)
|
(-5.119)
|
(-9.529)
|
(-8.693)
|
constant
|
3.465***
|
46.375***
|
0.825***
|
-1.036
|
|
(36.431)
|
(11.607)
|
(48.119)
|
(-1.412)
|
lnpgdp
|
|
-6.220***
|
|
0.011
|
|
|
(-17.614)
|
|
(0.163)
|
lnpd
|
|
1.773***
|
|
0.204***
|
|
|
(4.223)
|
|
(2.639)
|
secind
|
|
5.848***
|
|
1.119***
|
|
|
(10.783)
|
|
(11.229)
|
fdip
|
|
-2.171
|
|
-0.694**
|
|
|
(-1.302)
|
|
(-2.264)
|
eximp
|
|
0.096
|
|
0.079***
|
|
|
(0.799)
|
|
(3.560)
|
scieddp
|
|
0.254
|
|
0.152
|
|
|
(0.369)
|
|
(1.208)
|
lnwrcolstu
|
|
0.270***
|
|
0.031***
|
|
|
(8.111)
|
|
(5.145)
|
Control variables
|
No
|
Yes
|
No
|
Yes
|
Year fixed effects
|
Yes
|
Yes
|
Yes
|
Yes
|
City fixed effects
|
Yes
|
Yes
|
Yes
|
Yes
|
Observations
|
4845
|
4845
|
4845
|
4845
|
R-squared
|
0.065
|
0.140
|
0.595
|
0.613
|
Note: ⁎⁎⁎ and ⁎⁎ represent significance levels of 1% and 5%, respectively. The values in parentheses are obtained by robust t-statistic. |
Robustness tests
A series of auxiliary tests were performed to ensure the robustness of the results, including mitigating the bias of non-random selection, excluding the interference from other environmental policies, placebo test, PSM-DID, and the hysteresis effect analysis of the LCCP policy.
The ideal sample case for the DID model is that pilot and non-pilot cities are randomly selected. However, the list of low-carbon pilot cities is not random, and is closely related to their corresponding geographical location, economic development and industrial structure. To control for estimation bias from these factors, the cross term of benchmark factors and temporal linear trends were added to Eq. (1), including whether it was a “two-control zone,” a provincial capital city, and a northern city requiring heating. Columns (1) and (2) of Table 3 show the regression results after introducing these variables. Although the magnitude of the coefficients differs slightly from Table 2, the direction and significance of the coefficients remained consistent with the benchmark model.
Assessment of the CO2 emission reduction effect of the LCCP policy is inevitably affected by other environmental policies, especially those implemented during the same period, leading to possible overestimation or underestimation. To address this problem, we collected and summarized environmental policies since 2010, which included the new energy vehicle subsidy pilot policy (nenervehicle), the atmospheric emission limit pilot policy (atomemi), and the carbon emission trading rights pilot policy (cemitra). The above environmental policies were included in the regression model with the cross term of time linear trend. The regression results after adding the other environmental policies’ dummy variables are displayed in columns (3) and (4) of Table 3. The coefficients of treat were similar to those of the benchmark regression. It should be noted that other environmental policies were not statistically significant, indicating that the above environmental policies did not bias the estimated results.
According to the results of the parallel trend test, the expected CO2 emission reduction effect in the current implementation year of the LCCP policy was not achieved. This led us to consider the possible temporal path-dependent characteristics of CO2 emissions. Following Chen et al. (2021), all explanatory variables were made to lag by one period to eliminate the possibility of reverse causality of the dependent variable on the independent variable. As shown in columns (5) and (6) of Table 3, the results remained robust after considering the hysteresis effect of the LCCP policy.
Table 3
The results of the robustness test
|
CI
|
CP
|
CI
|
CP
|
CI
|
CP
|
treat
|
-0.563***
|
-0.146***
|
-0.571***
|
-0.132***
|
-0.607***
|
-0.192***
|
|
(-4.796)
|
(-6.780)
|
(-4.872)
|
(-6.184)
|
(-5.024)
|
(-8.755)
|
nenervehicle×trend
|
|
|
0.040
|
-0.008
|
|
|
|
|
|
(0.877)
|
(-0.999)
|
|
|
atomemi×trend
|
|
|
0.055
|
0.005
|
|
|
|
|
|
(0.991)
|
(0.539)
|
|
|
cemitra×trend
|
|
|
-0.007
|
-0.026
|
|
|
|
|
|
(-0.045)
|
(-0.950)
|
|
|
Control variables
|
Yes
|
Yes
|
Yes
|
Yes
|
Yes
|
Yes
|
Year fixed effects
|
Yes
|
Yes
|
Yes
|
Yes
|
Yes
|
Yes
|
City fixed effects
|
Yes
|
Yes
|
Yes
|
Yes
|
Yes
|
Yes
|
Observations
|
4845
|
4845
|
4845
|
4845
|
4845
|
4845
|
R-squared
|
0.143
|
0.618
|
0.142
|
0.621
|
0.142
|
0.559
|
Note: ⁎⁎⁎ represents significance levels of 1%. The values in parentheses are obtained by robust t-statistic. |
A series of crucial observable city characteristics were added to the benchmark model, including economic growth, population density, industrial structure, openness, human capital, and science and education. However, the effect of unobservable characteristics was not controlled by the model. To solve this problem, a placebo test was used to assess whether any omitted variables would affect the obtained results. The estimated coefficients were expressed as:
$${\hat {\alpha }_1}{\text{=}}{\alpha _1}{\text{+}}\gamma \times \frac{{\operatorname{cov} (trea{t_{it}},{\varepsilon _{it}}|x)}}{{\operatorname{var} (trea{t_{it}}|x)}}$$
3
Where x comprises of all control variables and fixed effects, and γ represents the effect of unobservable factors on the explanatory variables. When γ = 0, the unobserved factors do not affect the estimation results. Based on this, we randomly generated the list of low-carbon pilot cities (by computer).
Figure 3 plots the distribution of \({\hat {\alpha }_1}^{{random}}\) (500 replications). The \({\hat {\alpha }_1}^{{random}}\) distribution was in the vicinity of zero and obeyed a normal distribution, as expected from the placebo test, and again, demonstrated the robustness of the results of this study. Meanwhile, to alleviate the impact of sample selection bias and systematic differences, the PSM-DID method was employed to estimate the robustness of the results. The results of PSM-DID further supported the previous benchmark regression (Appendix 1. and Appendix 2.). Based on the above test, it could be deduced that the observed CO2 emission reduction of the 285 cities were derived from the implementation of the LCCP policy.
Endogeneity test using the IV method
The instrumental variable approach is applied to overcome endogeneity as much as possible. Specifically, the instrumental variable is chosen to satisfy the two conditions of being correlated with the endogenous variables and uncorrelated with the random disturbance terms (Shi and Li, 2020). The air circulation coefficient (VC) and relief degree of land surface (RDLS) were selected as the instrumental variables. First, the VC and RDLS are determined by meteorological and geographical conditions, satisfying the exogeneity hypothesis. Second, cities with smaller VC typically are considered to typically adopt stricter environmental regulations, while cities with lower RDLS tend to be more densely populated and economically active. Such cities have a higher probability of being selected as low-carbon pilot cities, which is consistent with the hypothesis of correlation of instrumental variables. The VC coefficient is the product of the wind speed and atmospheric boundary layer height, and is the natural logarithm of the annual average coefficients. Table 4 reports the two-stage least squares (2SLS) regression results. Panel A reflects the first-stage regression results of VC and RDLS on the LCCP policy, and indicates that VC and RDLS were significantly correlated with the LCCP policy. In addition, we observed that the F-statistics were greater than 10, thereby rejecting the hypothesis of “weak instrument variable”. Panel B also demonstrates that implementation of the LCCP policy significantly decreased CI and CP.
Table 4 2SLS results of instrumental variables
|
treat
|
CI
|
CP
|
treat
|
CP
|
CP
|
Panel A: First-stage
|
|
|
lnVC
|
0.132***
|
|
|
|
|
|
|
(758.111)
|
|
|
|
|
|
RDLS
|
|
|
|
0.616***
|
|
|
|
|
|
|
(59.107)
|
|
|
Panel B: Second-stage
|
|
|
treat
|
|
-0.549***
|
-0.183***
|
|
-0.847***
|
-0.141***
|
|
|
(-4.835)
|
(-8.785)
|
|
(-4.934)
|
(-4.465)
|
Control variables
|
Yes
|
Yes
|
Yes
|
Yes
|
Yes
|
Yes
|
Year fixed effects
|
Yes
|
Yes
|
Yes
|
Yes
|
Yes
|
Yes
|
City fixed effects
|
Yes
|
Yes
|
Yes
|
Yes
|
Yes
|
Yes
|
R-squared
|
0.994
|
0.140
|
0.613
|
0.546
|
0.139
|
0.613
|
F-statistics
|
10.72
|
|
|
10.41
|
|
|
Note: ⁎⁎⁎ represents significance levels of 1%. The values in parentheses are obtained by robust t-statistic.
Mechanism identification
The previous staggered DID estimation results and a series of robustness tests confirmed that the LCCP policy significantly reduced CI and CP, but how could this effect be achieved? Identification of its underlying mechanism is required. Based on the implementation paths of the LCCP policy, we investigated the mechanism of CO2 emission reduction from the LCCP policy using industrial emission reduction, technological innovation, and energy usage data. Following the practice of existing literature (Li et al.,2018), we assumed that the underlying mechanisms were from variables regressed by policy variables, and constructed the following regression model:
$${X_{it}}={\alpha _0}+{\alpha _1}trea{t_{it}}+{\lambda _i}+{\tau _t}+{\sigma _{it}}$$
4
where Xit represents the matrix vector of the mechanism variables (lnindSO2, lnindfumes, wrlcpat, scireseap, energyint, renewstr). λ is the year fixed effect. τ is the city fixed effect. σ is the stochastic disturbance term. The other variables have the same meaning as in Eq. (1).
Table 5
The results of the mechanism test
|
lnindSO2
|
lnindfumes
|
wrlcpat
|
scireseap
|
energyint
|
renewstr
|
treat
|
-0.114***
|
-0.125***
|
1.054***
|
0.005***
|
-0.146***
|
0.032***
|
|
(-2.751)
|
(-3.180)
|
(11.293)
|
(3.909)
|
(-5.989)
|
(4.960)
|
Control variables
|
Yes
|
Yes
|
Yes
|
Yes
|
Yes
|
Yes
|
Year fixed effects
|
Yes
|
Yes
|
Yes
|
Yes
|
Yes
|
Yes
|
City fixed effects
|
Yes
|
Yes
|
Yes
|
Yes
|
Yes
|
Yes
|
Observations
|
4845
|
4845
|
4845
|
4845
|
4845
|
4845
|
R-squared
|
0.493
|
0.480
|
0.269
|
0.790
|
0.269
|
0.709
|
Note: ⁎⁎⁎ represents significance levels of 1%. The values in parentheses are obtained by robust t-statistic. |
The estimated results are shown in Table 5. Columns (1) and (2) show the results of the LCCP policy affecting CO2 emission reduction through industrial emission reduction. It can be observed that the LCCP policy significantly decreased industrial SO2 and fumes emissions. Columns (3) and (4) list the technological innovation results, and show that the LCCP policy significantly increased the number of low-carbon patent applications and researchers. It is worth noting that the enhancement of low carbon patent was most effective in reducing CO2 emissions. These results provide new empirical evidence for the weak Porter hypothesis. Columns (5) and (6) show the impact of the LCCP policy on energy usage, and demonstrate that the LCCP policy could indeed achieve CO2 reduction effects by improving energy intensity and using renewable energy. These findings concord with those of Hong et al. (2021), and support Hypothesis 2.
Heterogeneity analysis
Heterogenous analysis of policy tool
A text quantification method was applied to build appropriate policy tool variables to further discuss the heterogenous effect of different LCCP policy tools on CO2 emission reduction. The LCCP policy tools are usually divided into the following three types: command-mandatory, market-economic, and voluntary policy tool (Wang et al., 2015). Following Chen et al. (2018), the frequency of words related to policy tools in provincial government work reports was selected as a proxy variable for LCCP policy tools on the city scale. The specific construction steps implemented were as follows: First, 30 provincial (excluding Tibet, Hong Kong, Macau and Taiwan) government work reports from 2003 to 2019 were manually collected from official government websites. Second, the texts of these government work reports were processed for word separation. Specifically, terms related to command-mandatory LCCP policy tools included elimination, control, restriction, prohibition, compulsory, standard, emission reduction, governance, permit. Key words related to market-economy LCCP policy tools were set as tax, fee, subsidy, compensation, penalty, financing, investment, credit, market, emission trade, renewable, clean, low carbon. Terms related to voluntary LCCP policy tools included pilot, park, industrial park, nature reserve, town, green, ecology, environmental protection, public transportation, energy usage. It is worth noting that selection of provincial variables to measure policy tools at city levels alleviated endogeneity, but reduced the city’s variability of the LCCP policy. Based on Chen and Chen (2018), we multiplied the proportion of cities’ secondary industry by the natural logarithm of the frequency of policy tools in provincial government work reports, and finally obtained three policy tool variables on the city scale.
Table 6
Heterogeneity results of the policy tools
|
CI
|
CP
|
CI
|
CP
|
CI
|
CP
|
control
|
-0.980*
|
-0.387***
|
|
|
|
|
|
(-1.816)
|
(-3.891)
|
|
|
|
|
market
|
|
|
-0.304
|
0.085
|
|
|
|
|
|
(-0.412)
|
(0.626)
|
|
|
voluntary
|
|
|
|
|
-0.911***
|
-0.387***
|
|
|
|
|
|
(-2.642)
|
(-6.094)
|
Control variables
|
Yes
|
Yes
|
Yes
|
Yes
|
Yes
|
Yes
|
Year fixed effects
|
Yes
|
Yes
|
Yes
|
Yes
|
Yes
|
Yes
|
City fixed effects
|
Yes
|
Yes
|
Yes
|
Yes
|
Yes
|
Yes
|
Observations
|
4845
|
4845
|
4845
|
4845
|
4845
|
4845
|
R-squared
|
0.136
|
0.608
|
0.137
|
0.607
|
0.137
|
0.610
|
Note: ⁎⁎⁎ and ⁎ represent significance levels of 1% and 10%, respectively. The values in parentheses are obtained by robust t-statistic. |
Table 6 presents the estimated results of different LCCP policy tools on CO2 emission reduction. The results of command-mandatory LCCP policy tools in columns (1) and (2) show significantly reduced CI and CP. The estimated coefficients absolute value of market-economic LCCP policy tools in columns (3) and (4) were smaller than command-mandatory LCCP policy tools, but the results were insignificant. This might be related to the effectiveness of market-economy LCCP policy tools, which are influenced by factors such as institutional design and degree of marketization (Xu and Cui, 2020). Columns (5) and (6) illustrate voluntary LCCP policy tools’ negative CO2 emission reduction effect at the 1% significance level. Compared to market-economic LCCP policy tools, command-mandatory and voluntary LCCP policy tools were more likely to reduce CO2 emission under the LCCP policy. This partly validates Hypothesis 3.
Heterogenous analysis of city development level
Considering that the response to LCCP policy could be heterogeneous among different city development levels, we referred to the latest “Ranking of Commercial Attractiveness of Chinese Cities” released by the China New First-tier Cities Research Institute, and graded the cities based on the following five dimensions: commercial resource concentration, city hub, city activity, lifestyle diversity and future plasticity, reflecting the comprehensive city development level and CO2 emission reduction potential. The 285 cities were classified and consolidated (merging first-tier and quasi-first-tier cities), rendering a list of new first-tier to fifth-tier cities. The number of cities in each tier was 19, 30, 70, 81 and 85, respectively. This study introduced the level of city development to the benchmark model, as shown in Eq. (5):
$${y_{it}}={\alpha _0}+{\alpha _1}trea{t_{it}} \times citydevelo{p_j}+\beta {X_{it}}+{\eta _i}+{\gamma _t}+{\varepsilon _{it}}$$
5
Where citydevelopj represents the level of city development, the other variables have the same meaning as the benchmark model.
Table 7
Heterogeneity results of city development level
|
CI
|
CP
|
CI
|
CP
|
CI
|
CP
|
CI
|
CP
|
CI
|
CP
|
first-tier×treat
|
-0.319
|
-0.355***
|
|
|
|
|
|
|
|
|
|
(-1.452)
|
(-8.808)
|
|
|
|
|
|
|
|
|
second-tier×treat
|
|
|
-0.853***
|
-0.261***
|
|
|
|
|
|
|
|
|
|
(-3.789)
|
(-6.285)
|
|
|
|
|
|
|
third-tier×treat
|
|
|
|
|
-0.384*
|
0.003
|
|
|
|
|
|
|
|
|
|
(-1.723)
|
(0.082)
|
|
|
|
|
fourth-tier×treat
|
|
|
|
|
|
|
0.024
|
-0.079
|
|
|
|
|
|
|
|
|
|
(0.093)
|
(-1.639)
|
|
|
fifth-tier×treat
|
|
|
|
|
|
|
|
|
-0.921***
|
-0.030
|
|
|
|
|
|
|
|
|
|
(-3.622)
|
(-0.640)
|
Control variables
|
Yes
|
Yes
|
Yes
|
Yes
|
Yes
|
Yes
|
Yes
|
Yes
|
Yes
|
Yes
|
Year fixed effects
|
Yes
|
Yes
|
Yes
|
Yes
|
Yes
|
Yes
|
Yes
|
Yes
|
Yes
|
Yes
|
City fixed effects
|
Yes
|
Yes
|
Yes
|
Yes
|
Yes
|
Yes
|
Yes
|
Yes
|
Yes
|
Yes
|
Observations
|
4845
|
4845
|
4845
|
4845
|
4845
|
4845
|
4845
|
4845
|
4845
|
4845
|
R-squared
|
0.135
|
0.297
|
0.138
|
0.300
|
0.135
|
0.298
|
0.135
|
0.297
|
0.137
|
0.302
|
Note: ⁎⁎⁎ and ⁎ represent significance levels of 1% and 10%, respectively. The values in parentheses are obtained by robust t-statistic. |
Table 7 lists the heterogeneous results of cities with different development levels on CO2 emission reduction. Columns (1) and (2) indicate that first-tier cities’ LCCP policy significantly reduced CP, with an insignificant effect on CI, possibly related to the crowding effect in cities whereby higher levels of city development outweighed the agglomeration effect and exacerbated CO2 emissions. As seen in columns (3) and (4), second-tier cities’ LCCP policy had a significantly negative effect on CI and CP. The results in columns (5) and (6), and columns (9) and (10) indicate that the LCCP policy of third-tier and fifth-tier cities was only effective on CI. Columns (7) and (8) indicate that fourth-tier cities’ LCCP did not significantly affect CO2 emission reduction. Thus, it could be seen that study demonstrates that the LCCP policy had the most significant CO2 emission reduction effect on second-tier cities.
SDID analysis
Spatial correlation analysis
The basic assumption in classical DID is the stable unit treatment value assumption (SUTVA) (Rubin, 2014), i.e., individuals from treatment group do not affect those of the control group. Such could ignore the impact of CO2 emission reduction on neighboring regions, leading to estimation errors (Anselin and Arribas-Bel, 2013). This assumption can be broken down if the spatial correlation is taken into account, or SUTVA no longer holds when there is a correlation between different spatial units, namely, the spatial spillover effect (Kolak and Anselin, 2019). Considering that spatial regional units do not exist in isolation, this means that the LCCP policy is inevitably affected by neighboring regions. Therefore, it is necessary to incorporate spatial dependence into DID in this study.
The premise of using the SDID model was to test for spatial autocorrelation of CI and CP via the global Moran’s index. The inverse of the squared geographic distance weight matrix (W1), economic distance weight matrix (W2) and economic geography nested weight matrix (W3) were constructed to estimate the spatial effect of the LCCP policy. Table 8 shows that the estimated results of the global Moran’s index were all significantly positive at the 1% level, indicating that the spatial autocorrelation on CI and CP were significant and spatial analysis could not be ignored.
Table 8
Global Moran’s index of CI and CP under three weighting matrices
|
CI
|
CP
|
|
W1
|
W2
|
W3
|
W1
|
W2
|
W3
|
2003
|
0.355***
|
0.520***
|
0.273***
|
0.347***
|
0.634***
|
0.286***
|
2004
|
0.398***
|
0.599***
|
0.350***
|
0.374***
|
0.662***
|
0.318***
|
2005
|
0.370***
|
0.470***
|
0.319***
|
0.362***
|
0.633***
|
0.304***
|
2006
|
0.321***
|
0.465***
|
0.283***
|
0.344***
|
0.628***
|
0.286***
|
2007
|
0.293***
|
0.346***
|
0.220***
|
0.339***
|
0.609***
|
0.279***
|
2008
|
0.343***
|
0.535***
|
0.281***
|
0.351***
|
0.630***
|
0.291***
|
2009
|
0.365***
|
0.610***
|
0.295***
|
0.356***
|
0.608***
|
0.293***
|
2010
|
0.382***
|
0.614***
|
0.320***
|
0.366***
|
0.660***
|
0.312***
|
2011
|
0.395***
|
0.663***
|
0.337***
|
0.383***
|
0.673***
|
0.326***
|
2012
|
0.395***
|
0.752***
|
0.330***
|
0.373***
|
0.661***
|
0.318***
|
2013
|
0.373***
|
0.838***
|
0.280***
|
0.386***
|
0.678***
|
0.330***
|
2014
|
0.395***
|
0.855***
|
0.291***
|
0.385***
|
0.685***
|
0.329***
|
2015
|
0.399***
|
0.849***
|
0.288***
|
0.391***
|
0.684***
|
0.330***
|
2016
|
0.388***
|
0.854***
|
0.283***
|
0.392***
|
0.676***
|
0.333***
|
2017
|
0.402***
|
0.812***
|
0.290***
|
0.412***
|
0.695***
|
0.352***
|
2018
|
0.397***
|
0.795***
|
0.291***
|
0.415***
|
0.700***
|
0.362***
|
2019
|
0.412***
|
0.818***
|
0.302***
|
0.399***
|
0.707***
|
0.350***
|
Note: ⁎⁎⁎ represents significance levels of 1%. |
Spatial spillover effect analysis
After testing for the existence of spatial autocorrelation, following Chagas et al. (2016), the SDID model was constructed to capture the spillover effects of the LCCP policy on CI and CP, which was essentially adding treat to the spatial econometrics model. Before estimating model coefficients, we compared two competing models, i.e., the Spatial Lag Model and the Spatial Error Model. The LM test was significant at the 10% or 1% level, and the null hypothesis of no spatial lag term and spatial autoregressive term was rejected, indicating that the influence of spatial relationship could not be ignored in the model. In addition, the LR test and Hausman test results show that the Spatial Durbin Model was suitable for the time and space dual fixed effects of this study.
Table 9 shows the results of the SDID model. The significance and direction of treat were as expected. However, the regression coefficients of the LCCP policy did not directly reflect the degree of impact on CI and CP, and the partial differential method was applied to decompose the spatial effects into direct and spatial spillover effects (Le Sage and Pace, 2009). Columns (1)-(3) show that the spatial spillover effects of the LCCP policy on CI under W1 was significantly negative and was 3.398 times greater than the direct effect. Further, the spatial spillover effect was not significant under W2 and W3, illustrating that the spatial spillover effect of the LCCP policy on CI was generated based on geographic distance. Columns (4)-(6) indicate that the spillover effects of the LCCP policy on CI was significant under the three weight matrices. Thus, Hypothesis 4 is verified.
Table 9
Spatial effect decomposition of SDID
|
CI
|
CP
|
|
W1
|
W2
|
W3
|
W1
|
W2
|
W3
|
treat
|
-0.567***
|
-0.479***
|
-0.533***
|
-0.150***
|
-0.155***
|
-0.149***
|
|
(-5.620)
|
(-4.890)
|
(-5.240)
|
(-7.950)
|
(-8.120)
|
(-7.910)
|
LR_Direct
|
-0.615***
|
-0.492***
|
-0.572***
|
-0.166***
|
-0.164***
|
-0.159***
|
|
(-5.590)
|
(-4.540)
|
(-5.060)
|
(-8.050)
|
(-7.960)
|
(-7.670)
|
LR_Indirect
|
-1.683***
|
-0.120
|
-3.049
|
-0.564***
|
-0.083***
|
-0.679**
|
|
(0.875)
|
(-0.790)
|
(-1.580)
|
(-3.1400)
|
(-3.020)
|
(-2.200)
|
LR_Total
|
-2.298***
|
-0.612***
|
-3.621*
|
-0.730***
|
-0.247***
|
-0.838***
|
|
(0.920)
|
(-2.750)
|
(-1.830)
|
(-3.880)
|
(-6.110)
|
(-2.640)
|
LM error
|
2139.226***
|
1337.178***
|
1904.700***
|
3291.655***
|
1279.378***
|
3283.448***
|
Robust LM error
|
3.254*
|
2.303
|
33.363***
|
1570.985***
|
392.111***
|
1717.106***
|
LM lag
|
2536.263***
|
1410.723***
|
2109.699***
|
1724.61***
|
888.853***
|
1574.816***
|
Robust LM lag
|
400.290***
|
75.848***
|
238.361***
|
3.940**
|
1.585
|
8.474***
|
LR error
|
209.170***
|
226.11***
|
346.350***
|
109.560***
|
96.080***
|
224.840***
|
LR lag
|
291.230***
|
251.97***
|
143.670***
|
97.800***
|
111.450***
|
173.680***
|
Hausman
|
53.650***
|
53.650***
|
53.650***
|
100.940***
|
100.940***
|
100.940***
|
ρ
|
0.666***
|
0.350***
|
0.838***
|
0.697***
|
0.313***
|
0.818***
|
Control variables
|
Yes
|
Yes
|
Yes
|
Yes
|
Yes
|
Yes
|
Year fixed effects
|
Yes
|
Yes
|
Yes
|
Yes
|
Yes
|
Yes
|
City fixed effects
|
Yes
|
Yes
|
Yes
|
Yes
|
Yes
|
Yes
|
R-squared
|
0.105
|
0.115
|
0.048
|
0.598
|
0.596
|
0.567
|
Note: ⁎⁎⁎, ⁎⁎, and ⁎ represent significance levels of 1%, 5%, and 10%, respectively. The values in parentheses are obtained by robust t-statistic. |