1. Introduction
China has, all at once, achieved tremendous economic growth but has also witnessed severe environmental pollution under the heavy industrialization mode characterized by high consumption of resources and energy, which has attracted a great deal of focus at home and abroad. For the recent decade, however, China has adopted a series of creative environmental protection measures to control the environmental pollution caused by the heavy industrialization mode, such as the “Action plan for prevention and control of air pollution” released in 2013 and the “Three-year plan on defending the blue sky” announced in 2018. Meanwhile, China has stepped into a new energy-saving and material-saving recycling industrial system, which not only improves the utilization efficiency of resources but also reduces the waste emissions in the way of economic output. Benefiting from the effective energy policy measures and the active transformation of the industrial system, China’s environmental pollution has been controlled to some degree. According to the “China Air Quality Improvement Report (2013–2018)” released in 2019, the atmospheric pollutant concentrations were significantly reduced in 2018, which brings an improvement of the overall environmental air quality. As for the 74 cities which implemented the first batch of “Environmental Air Quality Standards”, the average concentration of the sulfur dioxide (SO2) emissions decreased by 68%, and the average concentration of PM2.5 decreased by 42% between 2013 and 2018. In addition, the total emissions of nitrogen oxides and SO2 have fallen by 28% and 26% since 2013, respectively.
Despite the initial achievement, SO
2 emissions in China still surpass the sum of the members of the Organization for Economic Cooperation and Development (OECD) and the U.S. [
1]. According to the statistics in China Statistical Yearbook (2019), by the year 2018, in China’s total energy consumption, the share of fossil energy (including coal, crude oil, and natural gas) has remained above 85%, while the share of renewable energy (including hydropower, nuclear power, and wind power) has been below 15%. More specifically, in the year 2018, the consumption of coal represents 59% of the total level of energy consumption in China, which leads to the fact that the SO
2 emissions per capita have continuously increased. As the third-largest SO
2 emitter in the world, China is facing mounting domestic and external pressure to tackle the increasingly serious SO
2 pollution.
High SO
2 emission is a tremendous threat to the environment change. It is one of the most important on-going anxieties for both emerging countries and developed countries. As a typical traditional contaminant, SO
2 brings many adverse effects to the human body, such as breathing difficulty, pulmonary edema, eye irritation, asthma attacks, cardiopulmonary diseases [
2,
3,
4]. Meanwhile, SO
2 has considerable negative impacts on the ecological environment. For example, a high concentration of SO
2 will change the potential of hydrogen (pH) value of plants, which will lead to agricultural production reduction and forest death. Wei et al. [
5] found that in 899 Chinese counties, the estimated cost of agricultural losses induced by SO
2 air pollution reached USD 1.43 billion. According to the data from the Ministry of Ecology and Environment of China, in 2018, the proportion of acid rain (average precipitation pH less than 5.6) was 18.9%, 0.1 percentage point higher than the previous year; the proportion of cities with heavier acid rain (average precipitation pH below 5.0) was 4.9%, 1.8 percentage points lower than the previous year. Moreover, the type of acid rain was generally the sulfuric acid type, which means that SO
2 is still the main cause of acid rain.
Finding effective ways to control SO2 emissions levels has become one of the most invoking issues for environmental policymakers. China’s effective regulation of SO2 emissions not only improves its own environmental conditions but also contributes to the sustainable development of the entire world. In order to achieve this purpose and shape the effective energy policy to reduce SO2 emissions, policymakers should have a clear knowledge of whether the SO2 emission series is convergent or not. That is because convergence implies that the impact from the SO2 emissions per capita reduction policy is temporary and the SO2 emissions series would revert to a trend path in the long run. In view of this, the aim of the present paper is to investigate the time series property of the SO2 emissions per capita in 74 cities of China.
Numerous studies have been conducted to examine whether harmful gas is characterized as a random walk or mean-reverting process. The first study on pollutant emissions in the field of environmental economics is List [
6], who ascertained the convergence of SO
2 and nitrogen oxides emissions of ten Environmental Protection Agency regions of the U.S. over the period of 1929–1994. After that, the convergence of pollutant emissions between countries has been investigated, essentially for CO
2 and SO
2 emissions. Briefly speaking, the extant studies can be divided into two strands by surveying the convergence through various unit root tests and modeling the distribution for air pollutant emissions.
The first strand of the literature which focuses on mean reversion of air pollutants mostly implements the traditional unit root tests, such as the ADF test [
7], the DF–GLS test [
8], the PP test [
9], the KPSS test [
10], and the MZa test [
11]. For example, Strazicich and List [
12] utilize panel unit root tests and cross-section regressions to draw a conclusion that CO
2 emissions converge for 21 OECD nations between 1960 and 1997. Aldy [
13] utilized the univariate unit root ADF test to survey the CO
2 emissions in 88 selected countries. The results showed that the unit root null hypothesis of divergence in relative emissions can only be rejected in 13 countries, of which three are OECD countries. Using over a century of data across 28 developed and developing countries, Westerlund and Basher [
14] obtained support in favor of CO
2 emission convergence at the international level through the panel unit root test. Lee and Chang [
15] presented that the CO
2 emissions in 14 out of 21 OECD countries exhibited divergence using the DF test. Recent work also relies on some nonlinear unit root tests to examine the time series property of hazardous gas emissions. For instance, Payne et al. [
16] explored the stochastic convergence of SO
2 emissions per capita among U.S. states by conducting the residual augmented least squares–Lagrange multiplier (RALS–LM) unit root test with structural breaks. Yavuz et al. [
17] tested the CO
2 per capita emissions using the TAR panel unit root test and found that the United Kingdom is the transition country whose CO
2 per capita emissions determines the switch from one regime to the other. Li et al. [
18] utilized the panel KSS unit root test with a Fourier function, and the results indicated that CO
2 emissions only converge in 12 out of the 50 U.S. states in our analysis.
The other strand of the literature contributes to model the distribution of the data series to investigate the convergence. Wu et al. [
19] employed a continuous dynamic distribution approach and panel data of 286 cities at the prefecture and above-prefecture level; the results showed that per capita CO
2 emissions tend to converge during the sample period of 2002 to 2011. Herrerias [
20] utilized the distribution dynamics approach to analyze the convergence of CO
2 emissions per capita among the EU-25 countries between 1920 and 2007 and found that the convergence patterns differ strongly before and after World War II. Burnett [
21] utilized the two-stage procedure to examine the convergence of the states in the U.S. and obtained the conclusion that 26 states converged to a unique steady-state equilibrium. Yang et al. [
22] examined the SO
2 geographical distributions in 113 main cities of China and found that the cities located in the north of the country are heavily polluted, while cities with low pollution levels mainly agglomerate in the south. Zhou et al. [
23] studied the nexus of SO
2 emissions and economic development by employing the spatial panel model and suggested an inversely N-shaped environmental Kuznets curve. Yu et al. [
24] investigated the carbon emissions intensity convergence of 24 industrial sectors in China between 1995 and 2015, and based on an environmental performance index method and the convergence model, the results indicated that find the carbon intensities of all sectors converged to different steady levels. Ulucak and Apergis [
25] employed the club clustering approach to test for the convergence of ecological footprint by employing the annual data for the case of the European Union countries; the empirical results documented the presence of certain convergence clubs.
There is an increasing consensus that the conventional unit root tests do not consider the existence and the number of structural breaks in the data series, such as the ADF test and the traditional quantile unit root test proposed by Koenker and Xiao [
26]. The neglect of structural breaks may cause the efficiency of detecting the mean reversion of the data, and thus the ensuing results may not be convincing. In other words, ignoring the structural breaks sway the analysis toward accepting the null hypothesis of a unit root [
27]. The structural break tests were first introduced to the literature by Perron [
27] and tested by several authors [
28,
29,
30]. Lee and Chang [
31] further found that the CO
2 emissions per capita of 21 OECD countries were stationary using the panel unit root test with multiple breaks. Nonetheless, it is an enormous obstacle to accurately detect the locations of the estimated structural breaks and the number of breaks in the series. Taking the number and the specific dates for structural breaks into consideration, Becker et al. [
32] proposed a stationary test with a Fourier function. Christopoulos and Leon-Ledesma [
33] developed tests for unit roots that account jointly for structural breaks and nonlinear adjustment. Lee et al. [
34] and Meng and Lee [
35] developed a two-step LM and a three-step RALS-LM Fourier unit root test, respectively; they both utilize the Fourier function to control for a small number of smooth breaks of an unknown functional form and nonlinearity. Considering the merits and demerits, Bahmani-Oskooee et al. [
36] proposed a newly Fourier quantile unit root test to examine the integrational properties. This approach can solve the inaccurate inference generated by structural breaks and is able to test the unit root hypothesis in each quantile of data distribution and capture the type of asymmetric dynamics by allowing different speeds of adjustment at various quantiles of data distribution.
As a preliminary practice, we first utilize five conventional univariate unit root tests and the quantile unit test to investigate the stationarity of the SO
2 emissions per capita of 74 cities in China. Then we utilize the newly Fourier quantile unit test proposed by Bahmani-Oskooee et al. [
36] to re-investigate the convergence of the SO
2 emissions per capita from the perspective of both particular quantiles and overall conditions. It is able to capture the asymmetric dynamics by allowing different speeds of adjustment at various quantiles of the SO
2 emissions per capita distribution, regardless of whether the SO
2 emissions per capita of a city are above or below its steady-state. Therefore the economic implications would be suggested not only relying on whole quantiles, but also at each quantile. In addition, different from the quantile unit root test, the Fourier quantile unit root test could estimate the optimal frequency, which makes it able to deal with smooth breaks in time series without identifying beforehand the breaking numbers, breaking dates, or breaking forms.
The main contribution of this paper lies in two aspects. First, we conduct both the conventional and the newly proposed Fourier unit root tests to examine the SO2 emissions per capita convergence for 74 Chinese cities. As such, a robust conclusion regarding the convergence of SO2 emissions of different cities can be obtained through the comparison between the results of different unit root tests. Moreover, the 74 cities will be divided into two groups according to the time-series properties of the SO2 emissions. In doing so, differentiated policy measures could be proposed for different groups to combat the SO2 pollution. Second, we concentrate on the mean-reverting properties of SO2 emissions both at selected quantiles and at overall conditions. As such, shedding new light on previous literature that treats the data-generating process of the SO2 emissions as being linear, the potential asymmetric behavior of the SO2 emissions can be clearly detected and the possible smooth breaks in the data series can be fully accounted for. Empirical results indicate that the SO2 emissions per capita converge in 72 out of the 74 Chinese cities in our analysis during the period of December 2014 to June 2019.
3. Results and Discussion
For comparative purposes, we firstly implemented standard unit root tests, including the ADF test, the DF–GLS test, the PP test, the KPSS test, and the MZa test to assess the convergence of SO
2 emissions per capita of 74 cities in China.
Table 2 reports the corresponding results. As shown in this table, the DF–GLS and the MZa tests indicate that the SO
2 emissions per capita in 2 out of 74 cities are convergent at a 10% significant level. The KPSS test shows that the SO
2 emissions in 3 out of 74 cities are convergent. The ADF test supports that the SO
2 emissions per capita in 24 out of 74 cities are convergent at the usual significant level. These cities are Tianjin, Langfang, Cangzhou, Hengshui, Xingtai, Shijiazhuang, Baoding, Lianyungang, Zhenjiang, Quzhou, Jiaxing, Hangzhou, Shaoxing, Guangzhou, Chongqing, Xi’an, Shenyang, Changchun, Harbin, Yinchuan, Ji’nan, Zhengzhou, Dalian, and Qingdao. On the contrary, the PP test indicates that only the SO
2 emissions per capita in Zhenjiang are not converging, and all of the other cities are converging. Overall, the empirical results from the five conventional unit root tests are somewhat contradictory. The possible reason may lie in that these tests tend to fail to refuse the null hypothesis of the unit root when the data series exhibit structural breaks and/or non-normal distribution. Moreover, the traditional unit root tests can only provide the convergence over the whole sample; thus the mean-reverting properties cannot be revealed at particular quantiles of the series.
It is well known that univariate unit root tests may be inefficient when applied to finite samples. Next, this paper utilizes the quantile unit root test [
26] to revisit the convergence of SO
2 emissions per capita in concerned cities. This approach has the following two advantages. On the one hand, the quantile unit root test is more suitable to test the unit root hypothesis for the non-Gaussian series. On the other hand, the quantile unit root test provides unit root behavior both over the whole quantiles and at each selected quantile. As such, the asymmetric persistency could be uncovered through the quantiles.
To conduct the quantile unit root test beforehand, some parameters should be declared. Specifically, as mentioned earlier, the quantiles are determined by the range of , aiming to reveal the unit root behavior at different quantiles. In doing so, we can investigate whether the unit root behaviors of SO2 emissions per capita are asymmetric at different quantiles. In addition, the QKS statistics are employed to check the stochastic convergence over the whole quantiles. However, given that no standard distribution of and QKS statistics are available, we get help from bootstrap techniques with 5000 replications to generate critical values.
Table 3 and
Table 4 report the results from the quantile unit root tests.
Table 3 reports the results of the unit root hypothesis at particular quantiles. It is clear that from quantile 0.1 to 0.9,
is not rejected for the cities of Pearl River Delta; however, from the national level, the Pearl River Delta is undoubtedly far ahead in tackling SO
2 emissions and environmental governance. It is worth noting that the SO
2 emissions per capita are stationary at all quantiles for only one city, Lanzhou. On the contrary, the SO
2 emissions are divergent at all quantiles for 14 cities, that is, Shanghai, Lianyungang, Suzhou, Yancheng, Suqian, Zhenjiang, Guangzhou, Shenzhen, Zhaoqing, Kunming, Chengdu, Hefei, Changsha, and Xiamen. Second, as indicated by the QKS statistics from
Table 4, nine cities of the Pearl River Delta are all diverging. Nevertheless, although the quantile unit root test has significant superiority over the traditional unit root tests, it may have low efficiency of detecting the mean reversion of the data when structural breaks exist. In other words, if structural breaks exist in our concerned series, then the above results may not be persuasive.
To gain more insight, we also display the actual time paths of SO
2 emissions per capita and the fitted time paths by using the Fourier function in
Figure 1,
Figure 2,
Figure 3,
Figure 4,
Figure 5 and
Figure 6. It is worth noting that the actual nature of break(s) is generally unknown, and there is no specific guide as to where and how many breaks to use in the process of producing fitted SO
2 emissions per capita series. Despite this, we still clearly observe that the fitted series closely followed the actual series for a vast majority of the cities we considered. This implies that the quantile unit root test with a Fourier function has high power in detecting the potential structural breaks and assessing the mean-reverting properties and asymmetric behavior of the SO
2 emissions per capita series. Moreover, the distributions of the SO
2 emissions per capita in concerned 74 cities differ significantly from the normal distribution and the presence of smooth breaks. As a consequence, this paper shifts to test the stochastic properties of the SO
2 emissions per capita using the Fourier quantile unit root test. Similarly, the critical values are generated through the bootstrap technique with 5000 replications.
We first estimate the coefficients of the intercept
and autoregressive coefficient
, with the corresponding results reported in
Table 5. Note that
and
values are key indicators in determining the permanent/temporary effects of negative and positive shocks,
refers to the size of the shocks on each quantile, while
plays a crucial role in deciding on the mean reversion of SO
2 emissions in each quantile.
To compare the degree of persistence among quantiles, we plot the aggressive coefficient
in
Figure 7. In this figure, the solid line represents the values of
and dashed line represents a 95% confidence interval. Overall, we observe three types of patterns. First, the values of
for SO
2 emissions per capita series of Haikou, Shenzhen, Huizhou, Fuzhou, Ningbo, Qingdao, Jiaxing, Jinhua, and Ji’nan display concave or the straight-line upward patterns. The results indicate that the positive shocks to high quantiles of the SO
2 emissions per capita of these cities have a more persistent effect than negative shocks to low quantiles, indicating that the positive shocks to the long-run path of the urban–rural income gap will be unbound. The second type pattern can be approximated by the concave or the straight-line download pattern for Suqian and Kunming, which means positive shocks to the SO
2 emissions per capita have temporary effects while negative shocks have permanent effects for these two cities. In the third type pattern, the values of
for other cities display upward patterns in low quantiles and downward patterns in high quantiles. The results indicate that the high positive and negative shocks to SO
2 emissions per capita of these cities have a transitory effect.
Table 6 reports
, which are used to measure the degree of persistence in each quantile. We can clearly observe that at quantiles between 0.2 and 0.7, the statistic
is converging for 62 cities. However, at the extremely low quantile of 0.1,
cannot be rejected for Tianjin, Zhangjiakou, Chengde, Tangshan, Xingtai, Shanghai, Nanjing, Lianyungang, Changzhou, Yancheng, Suqian, Yangzhou, Zhuhai, Dongguan, Zhongshan, Zhaoqing, Chongqing, Kunming, Xi’an, Hefei, and Jinan, which means that the SO
2 emissions per capita in these 21 cities are divergent when the emissions are high. At the extremely high quantile of 0.9, the SO
2 emissions per capita are stationary for 33 cities. These results clearly show that the impact of external changes on SO
2 emissions per capita is non-linear and asymmetric and that there is more evidence to support the fact that it is more efficient when using the Fourier function in the model.
Table 7 reports the results of the unit root hypothesis over the whole quantiles. The optimal frequency and its
F statistics are reported at the end of the two columns. The results of Fourier QKS test statistics indicate that the null of unit root over the range of quantiles [0.1, 0.9] is rejected for 13 cities of Beijing–Tianjin–Hebei Urban Agglomeration, 24 cities of Yangtze River Delta, 9 cities of Pearl River Delta, the Chongqing municipality, 22 provincial capital cities, and 3 cities under separate state planning. Only the SO
2 emissions per capita of Nanjing and Hefei were not stationary. The optimum frequency varies from 0.1 to 4.8, and the computed
F statistics are in all cases greater than the critical values even at 1% significance level, indicating the validity of the choice of the optimal. Hence, the mean-reverting function with the nonlinear component is accepted in favor of the one without the nonlinear component.
Several important policy implications emerge from our study. First, the convergence of the SO2 emissions per capita for those 72 out of 74 cities does not mean that the environmental improvement goals have been achieved; it only suggests that the past environmental protection work has been effected. More to the point, stationarity indicates that the impact from the SO2 emissions per capita reduction policy is temporary and the SO2 emissions per capita series would revert to a trend path in the long run. Especially for 13 cities in Beijing–Tianjin–Hebei Urban Agglomeration, though the SO2 emissions per capita are converging, the SO2 emission levels are still larger than the levels of the Yangtze River and Pearl River Deltas. Therefore, the policymakers should find more creative and effective ways to curb the SO2 emissions, given that shifting the SO2 emissions from the current level to another would eventually return to its equilibrium level. The environmental situation is not optimistic. A large number of high-pollution industries, such as cement, steel, refining, and petrochemicals, are gathered and the industrial energy consumption is still dominated by coal. Therefore, it is still necessary to actively implement pollution control measures and adjust the structure of industrial sectors.
Second, Nanjing and Hefei are the only two cities where the SO2 emissions per capita are divergent. The environmental protection must be further strengthened because the environmental protection can leave permanent effects. The electricity, heat production and supply industries produce the largest amount of SO2 in these two cities. For a long time, there has been a phenomenon of uncoordinated development of light and heavy industries. The local government should appropriately adjust the industrial structure, increase the proportion of the tertiary industry, and reduce the proportion of the secondary industry, especially the proportion of low-end manufacturing, high-consumption, and high-pollution industries, so that its energy efficiency can be significantly improved and to reduce the SO2 emissions per capita.
Last, our findings suggest that the negative/positive shocks to the SO2 emissions per capita can produce three different effects. Negative shocks to the SO2 emissions per capita in Haikou, Shenzhen, Huizhou, Fuzhou, Ningbo, Qingdao, Jiaxing, Jinhua, and Ji’nan have transitory effects, while positive shocks have permanent effects. On the contrary, positive shocks to the SO2 emissions per capita in Suqian and Kunming have temporary effects, while negative shocks have permanent effects. For the other cities, neither positive nor negative shocks have permanent effects since the SO2 emissions per capita for these cities show stationary behavior. In light of the evidence, local governments should formulate environmental protection policies that are consistent with local practices based on the permanent or transient characteristics of the shocks to the SO2 emissions. This section may be divided by subheadings. It should provide a concise and precise description of the experimental results, their interpretation as well as the experimental conclusions that can be drawn.