Next Article in Journal
Variability of Soil Water Heat and Energy Transfer Under Different Cover Conditions in a Seasonally Frozen Soil Area
Previous Article in Journal
Characterization of Excavated Waste of Different Ages in View of Multiple Resource Recovery in Landfill Mining
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Convergence of Sulphur Dioxide (SO2) Emissions Per Capita in China

Department of Finance, School of Economics, Qingdao University, Qingdao 266071, China
*
Author to whom correspondence should be addressed.
Sustainability 2020, 12(5), 1781; https://doi.org/10.3390/su12051781
Submission received: 6 February 2020 / Revised: 23 February 2020 / Accepted: 23 February 2020 / Published: 27 February 2020
(This article belongs to the Section Energy Sustainability)

Abstract

:
As the third-largest SO2 emitter in the world, China is facing mounting domestic and external pressure to tackle the increasingly serious SO2 pollution. Figuring out the convergence and persistence of sulfur dioxide (SO2) emissions matters much for environmental policymakers in China. This study mainly utilizes the Fourier quantile unit root test to survey the convergence of the SO2 emissions per capita in 74 cities of China during the period of December 2014 to June 2019, by conducting five traditional unit root tests and a quantile root unit test as a comparative analysis. The empirical results indicate that the SO2 emissions per capita in 72 out of 74 cities in China are convergent in the sample period. The results also suggest that the unit root behavior of the SO2 emissions per capita in these cities is asymmetrically persistent at different quantiles. For the cities with the convergent SO2 emissions, the government should consider the asymmetric mean-reverting pattern of SO2 emissions when implementing environmental protection policies at different stages. For Hefei and Nanjing, the local governments need to enact stricter environmental protection policies to control the emission of sulfur dioxide.

Graphical Abstract

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, SO2 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 SO2 emissions per capita have continuously increased. As the third-largest SO2 emitter in the world, China is facing mounting domestic and external pressure to tackle the increasingly serious SO2 pollution.
High SO2 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, SO2 brings many adverse effects to the human body, such as breathing difficulty, pulmonary edema, eye irritation, asthma attacks, cardiopulmonary diseases [2,3,4]. Meanwhile, SO2 has considerable negative impacts on the ecological environment. For example, a high concentration of SO2 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 SO2 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 SO2 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 SO2 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 CO2 and SO2 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 CO2 emissions converge for 21 OECD nations between 1960 and 1997. Aldy [13] utilized the univariate unit root ADF test to survey the CO2 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 CO2 emission convergence at the international level through the panel unit root test. Lee and Chang [15] presented that the CO2 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 SO2 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 CO2 per capita emissions using the TAR panel unit root test and found that the United Kingdom is the transition country whose CO2 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 CO2 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 CO2 emissions tend to converge during the sample period of 2002 to 2011. Herrerias [20] utilized the distribution dynamics approach to analyze the convergence of CO2 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 SO2 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 SO2 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 CO2 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 SO2 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 SO2 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 SO2 emissions per capita distribution, regardless of whether the SO2 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.

2. Data and Method

2.1. Data Source

The analysis uses the monthly data of SO2 emissions per capita for 74 cities of China, which implements the first batch of “Environmental Air Quality Standards” from December 2014 to June 2019. According to the "Environmental Air Quality Standards", we further divided the 74 cities into 6 categories: Beijing–Tianjin–Hebei Urban Agglomeration, Yangtze River Delta, Pearl River Delta, municipalities, provincial capital cities, and under separate state planning. All of the data used are retrieved from the Chinese Environmental Monitoring Station.
Table 1 reports summary statistics of the SO2 emissions per capita in concerned 74 cities. As we can see, the maximum for the SO2 emissions per capita belongs to Yinchuan, where coal is the main energy resource. The minimum of the SO2 emissions per capita belongs to Beijing, Hainan, Suzhou, Yancheng, Jiaxing, Zhuhai, Zhongshan, and Changchun. We can also observe that most of the cities experience wide volatility of SO2 emissions. For instance, the maximum per capita SO2 emissions of Shenyang is 200, but the minimum is 10. In the penultimate column, we report the Jarque–Bera test statistic [37] and its significance. Clearly, all cities except Zhenjiang, Zhoushan, Shenzhen, Zhaoqing, Lhasa, Kunming, and Xiamen exhibit a clear sign of non-normal distribution, and it is strong evidence that favors the use of the Fourier quantile regression unit root test of Bahmani-Oskooee et al. [36].

2.2. Methodology

Koenker and Xiao [26] first propose the unit root test based on quantile regression. However, this method does not fully consider the impacts of structural breaks. Given that, Bahmani-Oskooee et al. [36] develop a quantile based unit root test with smooth breaks, which could approximate the unknown breaks in the series.
Firstly, assume y t , a data series of the SO2 emissions level which could be determined by a time-varying deterministic component d ( t ) and a stationary error term with variance δ 2 and zero mean, as follows:
y t = d ( t ) + ε t
In order to obtain a global approximation from the smooth transition and equip deterministic components with unknown breaks, the term d ( t ) could be expressed, as follows:
d ( t ) = c + a t + α k sin ( 2 π k t T ) + β k cos ( 2 π k t T )
The reason to conduct the Equation (2) in the model is based on the fact that a Fourier expression is capable of approximating absolutely integrable functions to any desired degree of accuracy. k represents the frequency of the Fourier function, α k and β k measures the amplitude and displacement of the frequency component, respectively, and π = 3.1416.
Then, the Equation (1) can be written as follows:
y t = c + a t + α k sin ( 2 π k t T ) + β k cos ( 2 π k t T ) + ε t
A desired feature of Equation (2) is that the standard linear specification emerges as a special case by setting α k = β k = 0 . It also follows that at least one frequency component must be present if there is a structural break. Here, to reject the null hypothesis α k = β k = 0 , the series must have a nonlinear component. Becker et al. [32] created a more powerful test to detect structural breaks under an unknown form. We set the maximum of K = 5 when we determine an optimal k. For any K = k, we estimate Equation (3) by employing the ordinary least squares (OLS) method and save the sum of squared residuals (SSR). Frequency is setting as the optimum frequency at the minimum of SSR. With the above assumption and respect to the deterministic components, we can test the following null hypothesis:
H 0 :   ξ t = ν t , ν t = ν t 1 + u t
To test the null hypothesis, Bahmani-Oskooee et al. [36] compute the OLS residuals as
e t = y t c ^ a ^ t α ^ k sin ( 2 π k t T ) β ^ k cos ( 2 π k t T )
Next, Bahmani-Oskooee et al. [36] used the quantile unit root test proposed by Koenker and Xiao [26] to investigate the stationarity of e t . The test is an extension of the ADF unit root test and has much more power than a standard ADF test when a given shock exhibits heavy-tailed behavior. Another advantage of the test is that it allows for different adjustment mechanism towards the long-run equilibrium at different quantiles. The ADF regression model on e t can be presented as follows:
e t = α 0 e t 1 + i = 1 p α i Δ e t i + ε t
where P is the lag order. α 0 is used to reflect the persistence degree. As usual, α 0 = 1 means that e t contains a unit root with persistency, and | α 0 | < 1, e t is required for the mean-reverting properties of SO2 emissions per capita and for ruling out explosive behavior. Equation (6) could be rewritten based on quantile regression, as follows:
Q e t ( τ | e t 1 , , e t p ) = Q ε ( τ ) + θ ( τ ) e t 1 + i = 1 p φ i Δ e t i
where Q ε ( τ ) is τ th conditional quantile of ν t and its estimated value captures the magnitude of the shock in each quantile, θ ( τ ) measures the speed of mean reversion within each quantile; here, the quantiles are set to be τ i ( 0.1 , 0.2 , 0.3 , 0.4 , 0.5 , 0.6 , 0.7 , 0.8 , 0.9 ) . Optimum lags are selected by the AIC information criteria. To obtain the coefficient θ ( τ ) and i = 1 p φ i , we minimize the following equation:
min t = 1 n ( τ I t ( e t < Q ε ( τ ) + θ ( τ ) e t 1 + i = 1 p φ i Δ e t i ) ) | e t Q ε ( τ ) θ ( τ ) e t 1 i = 1 p φ i Δ e t i |
where, I t ( ) = 1 if e t < Q ε ( τ ) + θ ( τ ) e t 1 + i = 1 p φ i Δ e t i , otherwise I t ( ) = 0 . Koenker and Xiao [26] further proposed a t-ratio statistic with the null non-stationary hypothesis α ( τ ) = 1 against different alternative hypotheses, α ( τ ) < 1 , α ( τ ) > 1 and α ( τ ) 1 , to check the unit root hypothesis at specific quantiles, which can be expressed as
t n ( τ i ) = f ( F ^ 1 ( τ i ) ) τ i ( 1 τ i ) ( Y 1 P ( 1 , Δ e t 1 , , Δ e t p ) Y 1 ) 1 2 ( θ ^ ( τ ) 1 )
where f ( ) is probability functions of e t , and F ( ) is the cumulative density function of series e t , Y 1 is the vector of lagged dependent variables ( e t 1 ) , and P X is the projection matrix onto the space orthogonal to X = ( 1 , Δ e t 1 , , Δ e t p ) . f ( F ^ 1 ( τ i ) ) is a consistent estimator of f ( F 1 ( τ i ) ) indicated by Koenker and Xiao [26], which can be expressed as
f ( F 1 ( τ i ) ) = ( τ i τ i 1 ) G ( ω ( τ i ) ω ( τ i 1 ) )
where ω ( τ i ) = ( c ( τ i ) , θ ( τ i ) , φ 1 ( τ i ) , , φ P ( τ i ) ) and τ i [ λ _ , λ ¯ ] . In this paper, we set λ _ = 0.1 and λ ¯ = 0.9 with interval 0.1. Obviously, we test the unit root hypothesis at different quantiles in comparison with the traditional ADF test, which only emphasizes on the conditional central tendency. To assess the unit root behavior over a range of quantiles, Koenker and Xiao [26] recommend the following Kolmogorov–Smirnov (QKS) test which could be presented as
Q K S = s u p τ i [ λ _ , λ ¯ ] | t n ( τ i ) |
We select the maximum of t n ( τ i ) to build the QKS–Fourier statistics over the quantiles τ i ( 0.1 , 0.2 , 0.3 , 0.4 , 0.5 , 0.6 , 0.7 , 0.8 , 0.9 ) . Though the limiting distributions of both t n ( τ i ) and the QKS tests are nonstandard, Koenker and Xiao [26] suggest using the re-sampling procedure to generate the critical values. Hence, to derive critical values for the above-mentioned test, we implement the re-sampling procedures of Koenker and Xiao [26] as follows:
We run following k-order autoregression by ordinary least square:
Δ e t = θ ( τ ) Δ e t 1 + t
We save fitted values Δ e ^ t = θ ( τ ) Δ e ^ t 1 and residuals ^ t , and then create bootstrap residuals ( t b ) with replacement from the centered residuals ^ t = ^ t 1 n 1 t = 1 + 1 n ^ t .
We calculate the bootstrap sample of observations as follows:
e t b = e t 1 b + Δ e t b
with { Δ e t b = k = 1 l θ ( τ ) Δ e t k b + t b ; Δ e t b = Δ e j   for   j = 1 , 2 , , l ; Δ e t b = e 1
In this paper, 5000 bootstrap iterations were used to accurately calculate the critical values, and the bootstrap variance estimator is inferior to more classical estimators. In order to carry out the empirical analysis, we need critical values for our tests. These are not available and must be constructed via the Monet Carlo simulation.

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 SO2 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 SO2 emissions per capita in 2 out of 74 cities are convergent at a 10% significant level. The KPSS test shows that the SO2 emissions in 3 out of 74 cities are convergent. The ADF test supports that the SO2 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 SO2 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 SO2 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 τ i ( 0.1 , 0.2 , 0.3 , 0.4 , 0.5 , 0.6 , 0.7 , 0.8 , 0.9 ) , 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 t _ n ( τ _ i ) 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, t n ( τ i ) 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 SO2 emissions and environmental governance. It is worth noting that the SO2 emissions per capita are stationary at all quantiles for only one city, Lanzhou. On the contrary, the SO2 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 SO2 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 SO2 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 SO2 emissions per capita series. Moreover, the distributions of the SO2 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 SO2 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 α 0 ( τ ) and autoregressive coefficient ρ 1 ( τ ) , with the corresponding results reported in Table 5. Note that α 0 ( τ ) and ρ 1 ( τ ) values are key indicators in determining the permanent/temporary effects of negative and positive shocks, α 0 ( τ ) refers to the size of the shocks on each quantile, while ρ 1 ( τ ) plays a crucial role in deciding on the mean reversion of SO2 emissions in each quantile.
To compare the degree of persistence among quantiles, we plot the aggressive coefficient ρ 1 ( τ ) in Figure 7. In this figure, the solid line represents the values of ρ 1 ( τ ) and dashed line represents a 95% confidence interval. Overall, we observe three types of patterns. First, the values of ρ 1 ( τ ) for SO2 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 SO2 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 SO2 emissions per capita have temporary effects while negative shocks have permanent effects for these two cities. In the third type pattern, the values of ρ 1 ( τ ) 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 SO2 emissions per capita of these cities have a transitory effect.
Table 6 reports t n ( τ i ) , 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 t n ( τ i ) is converging for 62 cities. However, at the extremely low quantile of 0.1, t n ( τ i ) 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 SO2 emissions per capita in these 21 cities are divergent when the emissions are high. At the extremely high quantile of 0.9, the SO2 emissions per capita are stationary for 33 cities. These results clearly show that the impact of external changes on SO2 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 SO2 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.

4. Conclusions

In this study, we mainly utilized the Fourier quantile unit root test to survey the convergence of the SO2 emissions per capita in 74 cities of China from December 2014 to June 2019, by conducting five traditional unit root tests and a quantile root unit test as a comparative analysis. The empirical results indicate that the SO2 emissions per capita in 72 out of 74 cities in China are converged in the sample period. Only the SO2 emissions per capita in Nanjing and Hefei are divergent. The results also suggest that unit root behavior of the SO2 emissions per capita in these cities is asymmetrically persistent at different quantiles.
Our findings have great implications for the environmental improvement of all cities in China. For the cities with the converging SO2 emissions, the government should consider the asymmetric mean-reverting pattern of SO2 emissions when implementing environmental protection policies at different stages. Moreover, it is worth noting that the environmental protection policies in these cities will not always make a consistent impact on SO2 emissions moving from lower to upper quantiles. For Hefei and Nanjing, the local governments need to enact stricter environmental protection policies to control the emission of sulfur dioxide. In addition, given that the old industrialization route is not desirable, it is thus necessary to implement an energy-saving and material-saving recycling industrial system, which reduces waste discharge per unit of economic output, improves resource utilization efficiency, and achieves low resource consumption and less environmental pollution. Only when the SO2 emissions per capita from cities reach convergence can the SO2 emissions per capita from the country reach convergence.

Author Contributions

Conceptualization, D.-K.S.; methodology, D.-K.S.; formal analysis, B.Z.; data curation, Y.-C.Z.; writing—original draft preparation, D.-K.S. and B.Z.; writing—review and editing, Y.-C.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Social Science Foundation of China (Grant number: 18CJY056).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Liu, C.; Hong, T.; Li, H.; Wang, L. From club convergence of per capita industrial pollutant emissions to industrial transfer effects: An empirical study across 285 cities in China. Energy Policy 2018, 121, 300–313. [Google Scholar] [CrossRef]
  2. Lin, C.; Li, C.; Yang, G.; Mao, I. Association between maternal exposure to elevated ambient sulfur dioxide during pregnancy and term low birth weight. Environ. Res. 2004, 96, 41–50. [Google Scholar] [CrossRef] [PubMed]
  3. Khan, R.R.; Siddiqui, M. Review on effects of particulates; sulfur dioxide and nitrogen dioxide on human health. Int. Res. J. Environ. Sci. 2014, 3, 70–73. [Google Scholar]
  4. Goudarzi, G.; Geravandi, S.; Idani, E.; Hosseini, S.A.; Baneshi, M.M.; Yari, A.R.; Vosoughi, M.; Dobaradaran, S.; Shirail, S.; Marzooni, M.B. An evaluation of hospital admission respiratory disease attributed to sulfur dioxide ambient concentration in Ahvaz from 2011 through 2013. Environ. Sci. Pollut. Res. 2016, 23, 22001–22007. [Google Scholar] [CrossRef] [Green Version]
  5. Wei, J.; Guo, X.; Marinova, D.; Fan, J. Industrial SO2 pollution and agricultural losses in China: Evidence from heavy air polluters. J. Clean. Prod. 2014, 64, 404–413. [Google Scholar] [CrossRef] [Green Version]
  6. List, J.A. Have air pollutant emissions converged amongst U.S. regions? Evidence from unit-root tests. South. Econ. J. 1999, 66, 144–155. [Google Scholar] [CrossRef]
  7. Dickey, D.; Fuller, W. Distribution of the estimates for autoregressive time series with unit root. J. Am. Stat. Assoc. 1979, 74, 427–431. [Google Scholar] [CrossRef]
  8. Mackinnon, J.G. Numerical distribution functions for unit root and cointegration tests. J. Appl. Econom. 1996, 11, 601–618. [Google Scholar] [CrossRef] [Green Version]
  9. Phillips, P.C.B.; Perron, P. Testing for a unit root in time series regression. Biometrika 1988, 75, 335–346. [Google Scholar] [CrossRef]
  10. Kwiatkowski, D.; Phillips, P.; Schmidt, P.; Shin, Y. Testing the null hypothesis of stationarity against the alternative of a unit root: How sure are we that economic time series have a unit root? J. Econ. 1992, 54, 159–178. [Google Scholar] [CrossRef]
  11. Ng, S.; Perron, P. Lag length selection and the construction of unit root tests with good size and power. Econometrica 2001, 69, 1519–1554. [Google Scholar] [CrossRef] [Green Version]
  12. Strazicich, M.C.; List, J.A. Are CO2 emission levels converging among industrial countries? Environ. Resour. Econ. 2003, 24, 263–271. [Google Scholar] [CrossRef]
  13. Aldy, J.E. Per capita carbon dioxide emissions: Convergence or divergence? Environ. Resour. Econ. 2006, 33, 533–555. [Google Scholar] [CrossRef] [Green Version]
  14. Westerlund, J.; Basher, S.A. Testing for convergence in carbon dioxide emissions using a century of panel data. Environ. Resour. Econ. 2008, 40, 109–120. [Google Scholar] [CrossRef]
  15. Lee, C.C.; Chang, C.P. New evidence on the convergence of per capita carbon dioxide emissions from panel seemingly unrelated regressions augmented Dickey-Fuller tests. Energy 2008, 33, 1468–1475. [Google Scholar] [CrossRef]
  16. Payne, J.E.; Miller, S.; Lee, J.; Cho, M.H. Convergence of per capita sulphur dioxide emissions across US states. Appl. Econ. 2014, 46, 1202–1211. [Google Scholar] [CrossRef]
  17. Yavuz, C.N.; Yilanci, V. Convergence in per capita carbon dioxide emissions among g7 countries: A TAR panel unit root approach. Environ. Resour. Econ. 2013, 54, 283–291. [Google Scholar] [CrossRef]
  18. Li, X.L.; Tang, D.P.; Chang, T. CO2 emissions converge in the 50 US states—Sequential panel selection method. Econ. Model. 2014, 40, 320–333. [Google Scholar] [CrossRef]
  19. Wu, J.; Wu, Y.; Guo, X.; Cheong, T.S. Convergence of carbon dioxide emissions in Chinese cities: A continuous dynamic distribution approach. Energy Policy 2016, 91, 207–219. [Google Scholar] [CrossRef]
  20. Herrerias, M.J. CO2 weighted convergence across the EU-25 countries (1920–2007). Appl. Energy 2012, 92, 9–16. [Google Scholar] [CrossRef]
  21. Burnett, J.W. Club convergence and clustering of the US energy-related CO2 emissions. Resour. Energy Econ. 2016, 46, 62–84. [Google Scholar] [CrossRef] [Green Version]
  22. Yang, X.; Wang, S.; Zhang, W.; Zhan, D.; Li, J. The impact of anthropogenic emissions and meteorological conditions on the spatial variation of ambient SO2 concentrations: A panel study of 113 Chinese cities. Sci. Total Environ. 2017, 584, 318–328. [Google Scholar] [CrossRef]
  23. Zhou, Z.; Ye, X.; Ge, X. The impacts of technical progress on sulfur dioxide kuznets curve in China: A spatial panel data approach. Sustainability 2017, 9, 674. [Google Scholar] [CrossRef] [Green Version]
  24. Yu, S.; Hu, X.; Fan, J.L.; Cheng, J. Convergence of carbon emissions intensity across Chinese industrial sectors. J. Clean. Prod. 2018, 194, 179–192. [Google Scholar] [CrossRef]
  25. Ulucak, R.; Apergis, N. Does convergence really matter for the environment? An application based on club convergence and on the ecological footprint concept for the EU countries. Environ. Sci. Policy 2018, 80, 21–27. [Google Scholar] [CrossRef] [Green Version]
  26. Koenker, R.; Xiao, Z. Unit root quantile auto-regression inference. J. Am. Stat. Assoc. 2004, 99, 775–787. [Google Scholar] [CrossRef] [Green Version]
  27. Perron, P. The great crash, the oil price shock, and the unit root hypothesis. Econometrica 1989, 57, 1361–1401. [Google Scholar] [CrossRef]
  28. Zivot, E.; Andrews, D.W.K. Further evidence on the great crash, the oil-price shock, and the unit-root hypothesis. J. Bus. Econ. Stat. 1992, 10, 251–270. [Google Scholar] [CrossRef] [Green Version]
  29. Lee, J.; Strazicich, M.C. Minimum Lagrange multiplier unit root test with two structural breaks. Rev. Econ. Stat. 2003, 85, 1082–1089. [Google Scholar] [CrossRef]
  30. Lee, J.; Strazicich, M. Minimum LM Unit Root Test with One Structural Break. Econ. Bull. 2013, 33, 2483–2492. [Google Scholar]
  31. Lee, C.C.; Chang, C.P. Stochastic convergence of per capita carbon dioxide emissions and multiple structural breaks in OECD countries. Econ. Model. 2009, 26, 1375–1381. [Google Scholar] [CrossRef]
  32. Becker, R.; Enders, W.; Lee, J. A stationarity test in the presence of an unknown number of smooth breaks. J. Time Ser. Anal. 2006, 27, 381–409. [Google Scholar] [CrossRef]
  33. Christopoulos, D.K.; León-Ledesma, M.A. Smooth breaks and non-linear mean reversion: Post-Bretton-Woods real exchange rates. J. Int. Money Financ. 2010, 29, 1076–1093. [Google Scholar] [CrossRef] [Green Version]
  34. Lee, J.; Strazicich, M.C.; Meng, M. Two-step LM unit root tests with trend-breaks. J. Stat. Econom. Methods 2012, 1, 81–107. [Google Scholar]
  35. Meng, M.; Lee, J. More Powerful LM Unit Root Tests with Trend-Breaks in the Presence of Non-Normal Errors; University of Alabama: Tuscaloosa, AL, USA, 2012. [Google Scholar]
  36. Bahmani-Oskooee, M.; Chang, T.; Elmi, Z.; Ranjbar, O. Re-testing Prebisch-Singer hypothesis: New evidence using Fourier quantile unit root test. Appl. Econ. 2017, 1–14. [Google Scholar] [CrossRef]
  37. Jarque, C.M.; Bera, A.K. Efficient tests for normality, homoscedasticity and serial independence of regression residuals. Econ. Lett. 1980, 6, 255–259. [Google Scholar] [CrossRef]
Figure 1. The Plots of the SO2 emissions per capita and fitted nonlinearities of Beijing-Tianjin-Hebei Urban Agglomeration. (am) respectively represent the plots of the SO2 emissions per capita and fitted nonlinearities of Beijing, Tianjin, Zhangjiakou, Chengde, Qinhuangdao, Langfang, Tangshan, Cangzhou, Handan, Hengshui, Shijiazhuang, Baoding, and Xingtai.
Figure 1. The Plots of the SO2 emissions per capita and fitted nonlinearities of Beijing-Tianjin-Hebei Urban Agglomeration. (am) respectively represent the plots of the SO2 emissions per capita and fitted nonlinearities of Beijing, Tianjin, Zhangjiakou, Chengde, Qinhuangdao, Langfang, Tangshan, Cangzhou, Handan, Hengshui, Shijiazhuang, Baoding, and Xingtai.
Sustainability 12 01781 g001aSustainability 12 01781 g001b
Figure 2. The plots of the SO2 emissions per capita and fitted nonlinearities of Yangtze River Delta. (ay) respectively represent the plots of the SO2 emissions per capita and fitted nonlinearities of Shanghai, Nanjing, Wuxi, Nantong, Lianyungang, Suzhou, Huaian, Changzhou, Tàizhou, Yancheng, Suqian, Zhenjiang, Yangzhou, Xuzhou, Zhoushan, Lishui, Ningbo, Wenzhou, Quzhou, Jiaxing, Jinhua, Huzhou, Hangzhou, Shaoxing, and Táizhou.
Figure 2. The plots of the SO2 emissions per capita and fitted nonlinearities of Yangtze River Delta. (ay) respectively represent the plots of the SO2 emissions per capita and fitted nonlinearities of Shanghai, Nanjing, Wuxi, Nantong, Lianyungang, Suzhou, Huaian, Changzhou, Tàizhou, Yancheng, Suqian, Zhenjiang, Yangzhou, Xuzhou, Zhoushan, Lishui, Ningbo, Wenzhou, Quzhou, Jiaxing, Jinhua, Huzhou, Hangzhou, Shaoxing, and Táizhou.
Sustainability 12 01781 g002aSustainability 12 01781 g002bSustainability 12 01781 g002cSustainability 12 01781 g002d
Figure 3. The plots of the SO2 emissions per capita and fitted nonlinearities of Municipality-Chongqing.
Figure 3. The plots of the SO2 emissions per capita and fitted nonlinearities of Municipality-Chongqing.
Sustainability 12 01781 g003
Figure 4. The plots of the SO2 emissions per capita and fitted nonlinearities of Pearl River Delta. (ai) respectively represent the plots of the SO2 emissions per capita and fitted nonlinearities of Guangzhou, Zhuhai, Foshan, Dongguan, Zhongshan, Huizhou, Zhaoqing, Jiangmen, and Shenzhen.
Figure 4. The plots of the SO2 emissions per capita and fitted nonlinearities of Pearl River Delta. (ai) respectively represent the plots of the SO2 emissions per capita and fitted nonlinearities of Guangzhou, Zhuhai, Foshan, Dongguan, Zhongshan, Huizhou, Zhaoqing, Jiangmen, and Shenzhen.
Sustainability 12 01781 g004aSustainability 12 01781 g004b
Figure 5. The plots of the SO2 emissions per capita and fitted nonlinearities of provincial capital cities. (aw) respectively represent the plots of the SO2 emissions per capita and fitted nonlinearities of Lhasa, Haikou, Kunming, Fuzhou, Chengdu, Xi’an, Guiyang, Urumqi, Nanning, Xining, Shenyang, Changchun, Harbin, Wuhan, Yinchuan, Lanzhou, Nanchang, Hefei, Taiyuan, Ji’nan, Hohhot Changsha, and Zhengzhou.
Figure 5. The plots of the SO2 emissions per capita and fitted nonlinearities of provincial capital cities. (aw) respectively represent the plots of the SO2 emissions per capita and fitted nonlinearities of Lhasa, Haikou, Kunming, Fuzhou, Chengdu, Xi’an, Guiyang, Urumqi, Nanning, Xining, Shenyang, Changchun, Harbin, Wuhan, Yinchuan, Lanzhou, Nanchang, Hefei, Taiyuan, Ji’nan, Hohhot Changsha, and Zhengzhou.
Sustainability 12 01781 g005aSustainability 12 01781 g005bSustainability 12 01781 g005c
Figure 6. The plots of the SO2 emissions per capita and fitted nonlinearities of cities under separate state planning. (ac) respectively represent the plots of the SO2 emissions per capita and fitted nonlinearities of Xiamen, Dalian, and Qingdao.
Figure 6. The plots of the SO2 emissions per capita and fitted nonlinearities of cities under separate state planning. (ac) respectively represent the plots of the SO2 emissions per capita and fitted nonlinearities of Xiamen, Dalian, and Qingdao.
Sustainability 12 01781 g006
Figure 7. Autogressive coefficient ρ 1 ( τ ) . (avvv) respectively represent the autogressive coefficient of Lhasa, Haikou, Shenzhen, Huizhou, Zhuhai, Kunming, Fuzhou, Zhangjiakou, Xiamen, Zhoushan, Jiangmen, Guiyang, Dongguan, Chengde, Zhongshan, Táizhou, Lishui, Guangzhou, Nanning, Ningbo, Dalian, Wenzhou, Zhaoqing, Foshan, Shanghai, Xining, Quzhou, Qingdao, Jiaxing, Changchun, Jinhua, Nanjing, Wuxi, Shaoxing, Nantong, Chongqing, Lanzhou, Nanchang, Lianyungang, Hohhot, Suzhou, Huzhou, Changsha, Shenyang, Hangzhou, Qinhuangdao, Huaian, Beijing, Changzhou, Harbin, Wuhan, Yinchuan, Yancheng, Tàizhou, Zhenjiang, Chengdu, Suqian, Ji’nan, Tianjin, Yangzhou, Hefei, Taiyuan, Langfang, Xuzhou, Cangzhou, Tangshan, Handan, Hengshui, Zhengzhou, Xingtai, Xi’an, Shijiazhuang, Baoding, and Urumqi.
Figure 7. Autogressive coefficient ρ 1 ( τ ) . (avvv) respectively represent the autogressive coefficient of Lhasa, Haikou, Shenzhen, Huizhou, Zhuhai, Kunming, Fuzhou, Zhangjiakou, Xiamen, Zhoushan, Jiangmen, Guiyang, Dongguan, Chengde, Zhongshan, Táizhou, Lishui, Guangzhou, Nanning, Ningbo, Dalian, Wenzhou, Zhaoqing, Foshan, Shanghai, Xining, Quzhou, Qingdao, Jiaxing, Changchun, Jinhua, Nanjing, Wuxi, Shaoxing, Nantong, Chongqing, Lanzhou, Nanchang, Lianyungang, Hohhot, Suzhou, Huzhou, Changsha, Shenyang, Hangzhou, Qinhuangdao, Huaian, Beijing, Changzhou, Harbin, Wuhan, Yinchuan, Yancheng, Tàizhou, Zhenjiang, Chengdu, Suqian, Ji’nan, Tianjin, Yangzhou, Hefei, Taiyuan, Langfang, Xuzhou, Cangzhou, Tangshan, Handan, Hengshui, Zhengzhou, Xingtai, Xi’an, Shijiazhuang, Baoding, and Urumqi.
Sustainability 12 01781 g007aSustainability 12 01781 g007bSustainability 12 01781 g007cSustainability 12 01781 g007d
Table 1. Descriptive statistics of the SO2 emissions per capita.
Table 1. Descriptive statistics of the SO2 emissions per capita.
CityMeanMaxMinSkewnessKurtosisJarque-Bera
Beijing–Tianjin–Hebei Urban Agglomeration
Beijing9.0733531.7865.80247.222 ***
Tianjin19.3647752.2838.301112.182 ***
Zhangjiakou19.6736271.6705.18036.445 ***
Chengde17.1095061.6606.01146.019 ***
Qinhuangdao27.3828091.5635.25834.094 ***
Langfang16.1275441.6665.05335.113 ***
Cangzhou31.69193131.6645.87244.300 ***
Tangshan40.12785161.2184.74920.596 ***
Handan35.45513191.6565.51939.687 ***
Hengshui24.3827971.7165.32739.394 ***
Xingtai42.127146131.6675.40138.667 ***
Shijiazhuang34.618127101.7576.00849.030 ***
Baoding35.673166102.3438.539120.657 ***
Yangtze River Delta
Shanghai12.7273551.6126.40150.343 ***
Nanjing15.3823380.9403.96210.223 ***
Wuxi16.5274761.4345.12129.168 ***
Nantong22.0005171.0424.74916.955 ***
Lianyungang20.3645351.1383.97314.050 ***
Suzhou14.3274231.2245.33826.266 ***
Huaian15.2554161.2405.00823.336 ***
Changzhou20.0734881.1533.79913.640 ***
Tàizhou15.9644551.1404.35916.144 ***
Yancheng13.5644231.5446.60751.684 ***
Suqian16.0004641.2154.55719.095 ***
Zhenjiang17.5823960.3492.3232.169
Yangzhou18.6734481.0333.84811.429 ***
Xuzhou26.4006971.0674.04212.926 ***
Zhoushan7.6731440.4182.7451.749
Lishui8.8182653.06514.176372.366 ***
Ningbo11.5453352.0379.038121.567 ***
Wenzhou11.6912961.5075.52035.366 ***
Quzhou13.7454441.8608.600103.589 ***
Jiaxing13.1095132.3579.794156.709 ***
Jinhua14.2185252.58111.418223.446 ***
Huzhou14.8733751.3025.23226.959 ***
Hangzhou11.8553441.9458.803111.852 ***
Shaoxing15.5456052.1248.696115.719 ***
Táizhou15.9644551.1404.35916.144 ***
Pearl River Delta
Guangzhou10.8912150.5054.2605.978 **
Shenzhen7.618124−0.1943.8131.860
Zhuhai7.6001730.8273.3136.491 **
Foshan13.1643561.7338.14688.208 ***
Jiangmen11.7823752.0558.730113.957 ***
Dongguan11.4182271.0854.31114.725 ***
Zhongshan9.9642430.9905.01918.330 ***
Huizhou8.6551661.3746.06038.768 ***
Zhaoqing14.8732670.4132.3592.5066
Municipality
Chongqing12.1272761.2384.57119.710 ***
Provincial capital cities
Lhasa7.7821340.0792.9580.061
Haikou5.582940.8573.8208.266 **
Kunming15.218228−0.0632.8300.103
Fuzhou6.3451040.7403.3705.334 *
Chengdu11.4732350.6863.7745.682 *
Xi’an18.2735651.3804.47522.437 ***
Guiyang13.3644451.5264.82729.006 ***
Urumqi13.2364462.1396.98878.370 ***
Nanning11.3642561.3125.53230.476 ***
Xining25.23662111.3384.94825.096 ***
Shenyang42.691200102.0286.56066.742 ***
Changchun25.20010431.4444.30022.972 ***
Harbin28.76415762.3699.199139.483 ***
Wuhan11.8003541.6956.57055.546 ***
Yinchuan45.63620191.6165.01633.245 ***
Lanzhou20.2185370.9432.7898.254 **
Nanchang14.5453440.9984.99518.256 ***
Hefei12.0002940.9613.3808.797 **
Taiyuan52.21818581.3884.00519.979 ***
Jinan30.83610181.5955.15733.976 ***
Hohhot26.3099681.5034.89728.964 ***
Changsha13.5453451.5186.64951.647 ***
Zhengzhou23.5097741.4695.23031.181 ***
Cities under separate state planning
Xiamen9.527154−0.0432.7000.223
Dalian20.6556861.5234.07923.918 ***
Qingdao17.8736041.5505.11732.299 ***
J–B denotes the Jarque–Bera Test for Normality. ***, * *, and * indicate significance at the 1%, 5%, and 10% levels, respectively.
Table 2. Results for univariate unit root test.
Table 2. Results for univariate unit root test.
CityADF TestDF–GLS TestPP TestKPSS TestMZa Test
Beijing–Tianjin–Hebei Urban Agglomeration
Beijing−5.107 ***1.156−3.124 **0.603 **0.776
Tianjin−2.2251.395−3.921 ***0.603 **0.803
Zhangjiakou−2.684 *2.938−3.193 **0.702 **0.960
Chengde−5.307 ***−1.329−2.992 **0.479 **−3.306
Qinhuangdao−2.656 *2.008−3.355 **0.550 **1.038
Langfang−1.3341.660−3.407 **0.611 **0.961
Cangzhou−0.3772.544−3.782 ***0.578 **1.016
Tangshan−3.033 **−0.757−3.069 **0.742 ***−1.293
Handan−4.272 ***−1.123−4.278 ***0.606 **−2.217
Hengshui−1.4332.245−3.339 **0.715 **0.965
Xingtai−0.2592.394−2.780 *0.671 **0.956
Shijiazhuang−0.6622.339−3.493 **0.575 **0.575
Baoding−1.9462.853−4.921 ***0.642 **0.777
Yangtze River Delta
Shanghai−4.217 ***−0.537−4.037 ***0.811 ***−0.355
Nanjing−3.330 **−0.899−3.330 **0.861 ***−1.587
Wuxi−3.860 ***−0.630−3.771 ***0.851 ***−0.576
Nantong−3.340 **−0.392−4.020 ***0.969 ***−0.568
Lianyungang−0.6211.510−3.363 **0.617 **1.239
Suzhou−3.863 ***−0.491−3.789 ***0.921 ***−0.303
Huaian−3.869 ***−0.690−3.829 ***0.907 ***−0.797
Changzhou−2.818 *−0.470−3.340 **0.898 ***−0.790
Tàizhou−3.554 **−0.312−3.554 **0.928 ***−0.319
Yancheng−4.402 ***−0.455−4.197 ***0.897 ***−0.171
Suqian−3.633 ***−0.630−3.621 ***0.732 **−0.674
Zhenjiang−1.870−0.158−2.5530.984 ***−0.169
Yangzhou−3.454 **−0.950−3.454 **0.896 ***−1.721
Xuzhou−3.002 **−0.469−3.012 **0.903 ***−0.319
Zhoushan−2.853 *−0.910−3.617 ***0.503 **−2.643
Lishui−5.909 ***−1.027−6.457 ***0.718 **−1.724
Ningbo−5.171 ***−0.667−4.949 ***0.845 ***−0.637
Wenzhou−4.612 ***−1.466−4.592 ***0.924 ***−3.925
Quzhou−2.0440.434−5.206 ***0.947 ***0.388
Jiaxing−2.271−0.667−5.064 ***0.777 ***−0.592
Jinhua−5.414 ***0.576−5.566 ***0.804 ***0.575
Huzhou−4.495 ***−1.236−4.538 ***0.629 **−3.018
Hangzhou−2.092−0.594−4.918 ***0.802 ***−0.490
Shaoxing−1.363−0.703−5.091 ***0.890 ***−0.655
Táizhou−4.250 ***−0.736−4.300 ***0.759 ***−1.764
Pearl River Delta
Guangzhou−1.6600.044−3.651 ***0.675 **−0.024
Shenzhen−2.921 **−1.142−2.835 *0.708 **−3.214
Zhuhai−4.040 ***−1.259−3.752 ***0.438 *−3.024
Foshan−5.303 ***−1.021−5.272 ***0.965 ***−1.805
Jiangmen−5.066 ***−0.773−5.066 ***0.771 ***−0.912
Dongguan−3.901 ***−1.459−3.934 ***0.472 **−4.051
Zhongshan−4.146 ***−0.986−4.162 ***0.631 **−1.912
Huizhou−4.687 ***−1.136−4.692 ***0.446 *−2.418
Zhaoqing−3.201 **−0.558−2.963 **0.902 ***−1.499
Municipality
Chongqing−0.7663.774−2.928 **0.836 ***1.415
Provincial capital cities
Lhasa−3.949 ***−1.792 *−3.940 ***0.725 ***−5.274 *
Haikou−5.305 ***−1.426−5.254 ***0.179−4.341
Kunming−3.728 ***−0.326−3.728 ***0.875 ***−0.391
Fuzhou−3.865 ***−1.807 *−3.811 ***0.168−6.351 *
Chengdu−3.504 **2.012−3.528 **0.914 ***1.726
Xi’an−0.1643.034−3.439 **0.486 **1.073
Guiyang−5.014 ***0.711−3.673 ***0.354 *0.706
Urumqi−5.174 ***−0.955−4.317 ***0.444 **−2.516
Nanning−4.653 ***−1.193−4.682 ***0.560 **−2.728
Xining−4.331 ***−1.361−4.221 ***0.732 **−3.370
Shenyang−1.5073.305−3.336 **0.489 **0.800
Changchun0.2684.082−2.960 **0.412 *0.916
Harbin−2.3401.988−4.129 ***0.360 *0.716
Wuhan−3.620 ***−0.799−3.558 ***0.681 **−1.040
Yinchuan0.5752.847−3.341 **0.489 **0.837
Lanzhou−5.535 ***0.364−3.257 **0.0790.905
Nanchang−4.351 ***0.748−4.343 ***0.971 ***0.743
Hefei−2.733 *−0.629−2.790 *0.833 ***−0.775
Taiyuan−4.617 ***−1.414−3.260 **0.489 **−3.783
Jinan−1.3862.821−3.372 **0.749 ***0.982
Hohhot−7.009 ***1.535−3.001 **0.357 *0.900
Changsha−3.731 ***−0.641−3.799 ***0.913 ***−0.670
Zhengzhou−1.9641.076−3.293 **0.687 **1.007
Cities under separate state planning
Xiamen−3.365 **0.146−3.358 **0.468 **0.661
Dalian−0.5981.113−2.997 **0.563 ***0.787
Qingdao−1.2021.456−3.538 **0.810 ***0.984
***, **, and * indicate significance at the 1%, 5%, and 10% levels.
Table 3. Results for quantile unit root at particular quantiles.
Table 3. Results for quantile unit root at particular quantiles.
Quantiles0.10.20.30.40.50.60.70.80.9
Beijing–Tianjin–Hebei Urban Agglomeration
Beijing−10.52 ***−5.475 ***−4.639 ***−4.235 ***−2.24 **−1.776−1.1940.150.953
Tianjin−6.84 ***−6.192 ***−5.882 ***−5.245 ***−4.496 ***−2.855 ***−2.29 **−0.2870.597
Zhangjiakou−3.112 *−2.451 *−3.526 ***−3.357 ***−2.345 **−2.129−1.322−0.465−0.577
Chengde−3.479 *−3.411 ***−4.269 ***−4.635 ***−3.589 ***−3.526 ***−2.32−1.957−1.405
Qinhuangdao−5.57 **−6.345 ***−5.41 ***−3.182 ***−2.798 ***−2.083*−1.266−0.657−0.526
Langfang−10.609 ***−6.047 ***−4.684 ***−5.247 ***−4.13 ***−3.373 ***−2.565 **−1.147−0.046
Cangzhou−4.782 **−3.376 ***−2.964 **−1.651−1.884−1.736−1.807−1.427−1.03
Tangshan−2.519−2.443 *−3.034 **−2.439 *−2.576 **−2.802 **−2.27 *−2.604 **−0.134
Handan−7.307 ***−6.266 ***−5.361 ***−4.84 ***−3.744 ***−3.139 **−1.423−0.1540
Hengshui−10.509 ***−6.028 ***−3.332 ***−3.445 ***−3.714 ***−3.03 ***−1.908 *−1.486−0.103
Xingtai−5.389 ***−2.55 *−1.904−3.142 ***−3.86 ***−3.69 ***−3.834 ***−1.635 *−1.265
Shijiazhuang−5.175 **−3.774 ***−2.43 **−2.528 *−2.147 *−0.822−0.135−0.1−1.345
Baoding−8.481 ***−3.831 ***−3.355 ***−4.714 ***−3.869 ***−2.718 **−1.711−0.729−0.95
Yangtze River Delta
Shanghai−2.274−1.158−1.232−1.503−1.753−1.656−1.862−1.289−1.173
Nanjing−1.466−1.991−2.325 **−2.071−1.896−0.902−0.722−0.813−1.592
Wuxi−5.215 *−4.998 ***−4.127 ***−2.116−3.057 ***−2.123 *0.040.0761.14
Nantong0.017−1.423−1.779−2.361 *−2.431 *−1.382−0.432−0.0630.981
Lianyungang−3.703−2.057−1.664−1.163−1.588−1.241−1.074−0.5360.352
Suzhou−2.019−1.341−0.95−1.14−1.606−0.795−0.407−0.536−0.578
Huaian−1.367−2.524 **−2.641 ***−2.64 **−2.498 **−2.145 *−1.314−0.5720.696
Changzhou−5.183 *−3.006 **−2.646 **−2.271 *−0.962−1.093−0.5150.8242.896
Tàizhou−4.208 ***−1.799−1.139−0.794−0.854−0.687−0.724−0.9610
Yancheng−1.499−1.008−0.941−0.939−1.059−0.672−0.7460.7332.43
Suqian−2.221−0.589−0.761−1.1280.027−0.344−0.836−0.677−0.272
Zhenjiang−1.092−1.101−0.685−0.393−1.042−0.862−0.793−0.3060
Yangzhou−3.867 *−3.353 *−4.563 ***−1.508−0.999−0.1861.0372.0221.525
Xuzhou−6.415 **−1.922−1.139−1.1110.449−0.1151.0242.2273.719
Zhoushan−4.635 **−2.731 **−2.826 ***−1.979−0.669−0.6020.1920.281−0.677
Lishui−3.547 **−3.252 ***−2.666 **−2.163 *−1.486−1.57−1.878−1.347−1.305
Ningbo−1.251−1.925−1.475−2.115−3.368 ***−2.573 **−2.098−0.7260.526
Wenzhou−2.931 *−3.05 **−2.586 **−3.953 ***−3.554 ***−2.203 *−1.563−0.295−1.094
Quzhou−3.528−2.49 **−1.103−0.2190.259−0.388−0.745−0.804−0.776
Jiaxing−2.592−3.14 **−2.291 **−2.549 **−2.163 *−0.606−0.236−0.5070.171
Jinhua−7.144 ***−4.733 ***−4.079 ***−2.023−1.422−1.511−0.067−0.6110.11
Huzhou−1.795−2.894 ***−2.117 *−2.335 *−2.611 *−2.12 *−2.205 *−1.573−1.832
Hangzhou−2.474−2.182−2.085−0.969−0.476−0.188−0.926−0.6160.289
Shaoxing−7.672 ***−6.291 ***−3.063 ***−2.786 **−2.405 **−1.83−1.115−1.1990.515
Táizhou−2.889 *−2.786 *−2.63*−3.461 ***−2.719 **−1.347−1.536−1.2720.729
Pearl River Delta
Guangzhou−2.015−1.565−1.388−0.8920.081−0.0810.1230.82−0.42
Shenzhen−0.948−0.2170.5390.025−0.164−0.1110.9271.0550.483
Zhuhai−1.551−1.919−1.732−2.636 **−2.239 *−2.161 *−2.006−2.67 *−2.801
Foshan−3.191 *−2.11 **−1.811−1.757−2.121 **−2.294 *−0.936−0.2140.405
Jiangmen−5.17 *−2.563−1.581−1.233−0.753−0.775−0.0850.692.02
Dongguan−3.915 *−1.811−0.1−1.048−0.98−2.054 *−1.713−1.299−0.572
Zhongshan−2.926 *−1.149−1.879−1.794−1.149−0.878−1.148−0.811−3.235 **
Huizhou−1.53−1.444−1.321−1.892−2.561 *−2.294 *−1.705−1.651−2.497
Zhaoqing−1.981−1.256−0.995−0.653−0.9550.081−0.53−1.062−0.87
Municipality
Chongqing−3.672−3.549 ***−1.732−0.67−0.402−0.41−0.096−0.2861.339
Provincial capital cities
Lhasa−0.93300−1.68−1.608−1.947−2.236−2.885 ***−3.531 **
Haikou−2.962 *−2.255 *−2.58 **−3.314 **−3.409 ***−2.713 ***−1.4−1.538−1.545
Kunming1.8981.071−0.375−0.331−0.364−1.093−0.418−0.4380.572
Fuzhou−2.071−1.3−1.802−1.3−1.355−1.322−1.313−1.496−4.72 *
Chengdu−2.245−1.164−0.582−0.929−0.257−0.2620.7891.4180.216
Xi’an−9.242 ***−4.219 ***−4.575 ***−3.717 ***−2.762 **−2.635 **−3.056 **−1.323−2.221
Guiyang−8.762 ***−4.739 ***−4.5 ***−3.671 ***−2.782 **−1.778−1.973−0.548−1.056
Urumqi−14.573 ***−12.138 ***−9.883 ***−7.421 ***−5.396 ***−3.289 ***−1.162−1.0961.076
Nanning−5.557 **−3.668 ***−4.065 ***−3.968 ***−3.905 ***−3.021 ***−1.248−1.406−0.7
Xining−3.17 **−3.879 ***−4.413 ***−4.489 ***−4.1 ***−3.192 ***−1.89−1.84−1.692
Shenyang−9.234 ***−6.22 ***−4.56 ***−4.521 ***−4.247 ***−2.642 **−1.99 *−0.902−1.059
Changchun−4.409 **−4.076 **−4.351 ***−3.019 ***−3.241 ***−3.017 ***−1.07−0.895−1.757
Harbin−11.648 ***−4.805 ***−4.886 ***−4.012 ***−4.083 ***−2.644 **−2.369 **−1.667−0.542
Wuhan−3.652 **−2.796 ***−2.901 ***−1.69−1.665−1.467−1.166−1.674−2.894 **
Yinchuan−8.893 **−8.574 ***−6.84 ***−4.584 ***−2.668 **−1.168−0.977−0.001−2.456
Lanzhou−10.972 ***−6.794 ***−4.568 ***−3.926 ***−2.56 **−2.344 **−1.885 *−2.517 **−6.588 **
Nanchang−1.075−1.695−2.842 **−2.508 **−1.633−0.926−0.726−0.9922.525
Hefei−0.86−0.832−0.51−0.058−0.0320.335−0.611−0.695−1.148
Taiyuan−9.072 ***−7.119 ***−5.191 ***−5.652 ***−4.674 ***−2.959 ***−2.214 **−0.86−0.525
Jinan−10.373 ***−6.973 ***−4.522 ***−2.879 **−1.829 *−0.4570.0660.9972.837
Hohhot−17.668 ***−7.861 ***−6.927 ***−3.603 ***−2.949 **−2.473 **−3.135 ***−4.034 **−1.13
Changsha−3.025−2.052−1.610−1.647−0.346−0.323−0.489−0.522−1.301
Zhengzhou−4.585 *−1.908−2.177−2.096 *−2.188 *−0.691−0.755−0.493−0.021
Cities under separate state planning
Xiamen−1.832−1.051−0.775−0.565−0.427−0.2410.23600.288
Dalian−5.167 **−5.263 ***−6.797 ***−4.222 ***−3.998 ***−2.330 *−1.835 *−1.044−0.457
Qingdao−5.175 ***−3.374 ***−3.967 ***−3.024 **−3.776 ***−3.241 ***−1.132−0.984−0.051
***, **, and * indicate significance at the 1%, 5%, and 10% significant levels. The critical value for QKS statistics generated by bootstrap techniques with 5000 iterations.
Table 4. Results for overall quantile unit root test.
Table 4. Results for overall quantile unit root test.
CityQKS statisticCV 10%CV 5%CV 1%
Beijing–Tianjin–Hebei Urban Agglomeration
Beijing10.520 ***4.4675.3828.125
Tianjin6.840 **5.0216.38510.958
Zhangjiakou3.5264.9176.06510.408
Chengde4.6354.9536.0989.079
Qinhuangdao6.345 ***4.8025.8598.781
Langfang10.609 ***3.2993.8215.359
Cangzhou6.415 *5.6266.86011.124
Tangshan3.0343.8754.5296.218
Handan7.307 ***3.1723.6704.841
Hengshui10.509 ***3.3253.8485.046
Xingtai5.389 ***3.4103.8915.310
Shijiazhuang5.1755.5156.71510.049
Baoding8.481 ***3.8384.4616.227
Yangtze River Delta
Shanghai2.2744.5225.6079.027
Nanjing2.3253.3423.8715.236
Wuxi5.2155.5427.05411.447
Nantong2.4313.5223.9475.254
Lianyungang3.7035.4236.64910.699
Suzhou2.0195.3256.52710.989
Huaian2.6413.4033.8895.177
Changzhou5.183 *5.1726.47710.237
Tàizhou4.208 **3.4863.9655.240
Yancheng2.4305.5856.86910.477
Suqian2.2215.9427.33511.768
Zhenjiang1.1013.4723.9615.371
Yangzhou4.5635.3986.81710.275
Xuzhou6.415 *5.6266.86011.124
Zhoushan4.635 **3.7514.3406.433
Lishui3.5473.4713.9315.079
Ningbo3.3683.4994.0205.740
Wenzhou3.953 *3.9374.7466.786
Quzhou3.5285.4506.76410.200
Jiaxing3.1405.0136.3379.946
Jinhua7.144 ***3.4833.9155.459
Huzhou2.8944.2094.9856.907
Hangzhou2.4745.2686.5379.889
Shaoxing7.672 **4.5755.5878.526
Táizhou3.4614.2365.1117.533
Pearl River Delta
Guangzhou2.0153.8474.4746.300
Shenzhen1.0555.4756.84810.967
Zhuhai2.8014.2465.0577.525
Foshan3.1913.4383.9455.179
Jiangmen5.1705.8137.34212.113
Dongguan3.9155.2696.62210.757
Zhongshan3.2353.5294.1075.526
Huizhou2.5615.4026.73110.846
Zhaoqing1.9814.4485.4227.924
Municipality
Chongqing3.6725.1686.52310.111
Provincial capital cities
Lhasa3.5313.7884.4176.187
Haikou3.4093.7934.4386.207
Kunming1.8985.5927.01610.946
Fuzhou4.7204.7785.9049.215
Chengdu2.2455.9067.34011.347
Xi’an9.242 ***4.6105.7768.667
Guiyang8.762 **4.3705.4109.004
Urumqi14.573 ***3.4954.1856.248
Nanning5.557 **4.1654.9427.466
Xining4.489 **3.2593.6254.569
Shenyang9.234 ***4.1165.0108.102
Changchun4.4094.6225.6358.748
Harbin11.648 ***4.3435.2318.410
Wuhan3.652 *3.4853.9905.448
Yinchuan8.893 **4.7936.0729.976
Lanzhou10.972 ***4.8336.26910.536
Nanchang2.8423.7844.4126.042
Hefei1.1485.0636.3989.815
Taiyuan9.072 ***3.3683.9205.471
Jinan10.373 ***4.8786.1729.586
Hohhot17.668 ***4.4855.4248.622
Changsha3.0254.9426.0168.984
Zhengzhou4.5855.4356.88810.591
Cities under separate state planning
Xiamen1.8323.4944.0205.257
Dalian6.797 **4.2525.3418.122
Qingdao5.175 **3.4253.9685.524
***, **, and * indicate significance at the 1%, 5%, and 10% significant levels. The critical value for QKS statistics generated by bootstrap techniques with 5000 iterations.
Table 5. Intercept—α0(τ) and autoregressive coefficient—ρ1(τ).
Table 5. Intercept—α0(τ) and autoregressive coefficient—ρ1(τ).
Cities τ 0.10.20.30.40.50.60.70.80.9
Beijing–Tianjin–Hebei Urban Agglomeration
Beijingα0(τ)−3.895 **−2.173 ***−1.468 ***−0.957 **−0.3560.1030.758 **1.416 ***1.993 ***
ρ1(τ)0.491 ***0.599 ***0.597 ***0.639 ***0.691 ***0.729 ***0.647 ***0.599 ***0.580 ***
Tianjinα0(τ)−6.349 ***−4.166 ***−2.445 ***−1.888 **0.3530.911 *1.653 ***2.407 ***4.265 ***
ρ1(τ)0.701 ***0.654 ***0.662 ***0.699 ***0.639 ***0.650 ***0.609 ***0.646 ***0.514 ***
Zhangjiakouα0(τ)−5.389 ***−3.558 ***−2.387 ***−1.573 **−0.3480.1501.1772.776 ***5.197 ***
ρ1(τ)0.615 **0.662 ***0.756 ***0.725 ***0.797 ***0.828 **0.829 *0.801 **0.563 ***
Chengdeα0(τ)−4.725 ***−3.233 ***−2.275 ***−1.519 ***−1.178 **0.2211.194 **1.911 ***3.375 ***
ρ1(τ)0.741 **0.723 ***0.609 ***0.593 ***0.589 ***0.550 ***0.526 ***0.581 ***0.526 ***
Qinhuangdaoα0(τ)−6.836 ***−6.011 ***−4.069 ***−2.999 ***−1.390*0.1181.871 **2.635 ***5.437 ***
ρ1(τ)0.714 ***0.724 ***0.657 ***0.631 ***0.594 ***0.609 ***0.630 ***0.609 ***0.676 **
Langfangα0(τ)−6.214 ***−3.348 ***−1.805 **−1.337 **−0.722−0.2741.265 *3.046 ***3.746 ***
ρ1(τ)0.644 ***0.638 ***0.707 ***0.730 ***0.701 ***0.694 ***0.689 ***0.755 **0.805
Cangzhouα0(τ)−7.181 ***−5.402 ***−4.145 ***−2.558 **−0.6060.9791.607 **2.973 ***6.390 ***
ρ1(τ)0.708 **0.675 ***0.647 ***0.597 ***0.629 ***0.580 ***0.572 ***0.523 ***0.650 ***
Tangshanα0(τ)−11.385 ***−6.315 ***−3.173 **−1.006−0.0151.1212.256 **3.919 ***7.754 ***
ρ1(τ)0.612 **0.526 ***0.595 ***0.619 ***0.61 ***0.638 ***0.629 ***0.639 ***0.439 ***
Handanα0(τ)−20.679 ***−11.365 ***−7.262 ***−2.777 *−1.1770.5984.748 ***6.194 ***10.815 **
ρ1(τ)0.522 ***0.486 ***0.627 ***0.618 ***0.597 ***0.574 ***0.627 ***0.632 ***0.651 **
Hengshuiα0(τ)−7.285 ***−6.004 ***−4.627 ***−2.676 ***−1.255 *0.0301.917 *3.495 ***6.990 ***
ρ1(τ)0.635 ***0.651 ***0.575 ***0.577 ***0.638 ***0.692 ***0.787 **0.727 ***0.617 ***
Xingtaiα0(τ)−14.876 ***−11.235 ***−6.421 ***−3.568 **−1.4482.824 **3.420 **4.994 ***10.098 ***
ρ1(τ)0.675 **0.614 ***0.704 ***0.655 ***0.614 ***0.517 ***0.524 ***0.529 ***0.591 ***
Shijiazhuangα0(τ)−14.283 ***−8.549 ***−3.520 **−2.569 **−1.5321.3263.337 **5.803 ***11.716 ***
ρ1(τ)0.497 ***0.577 ***0.641 ***0.611 ***0.656 ***0.665 ***0.686 ***0.725 ***0.802 **
Baodingα0(τ)−12.584 ***−8.933 ***−5.448 ***−2.003−0.1301.2093.247 **5.072 **12.236 ***
ρ1(τ)0.669 ***0.697 ***0.738 ***0.808 ***0.797 ***0.752 ***0.734 ***0.714 ***0.647 ***
Yangtze River Delta
Shanghaiα0(τ)−2.261 ***−1.853 ***−1.415 ***−0.983 ***−0.524 *0.0580.5311.418 ***2.234 ***
ρ1(τ)0.631 ***0.677 ***0.614 ***0.592 ***0.692 ***0.724 ***0.633 ***0.481 ***0.565 ***
Nanjingα0(τ)−3.514 ***−2.388 ***−1.526 ***−1.146 ***−0.4860.3121.084 **1.478 **3.424 ***
ρ1(τ)0.647 **0.512 ***0.595 ***0.624 ***0.727 **0.8250.656 **0.604 **0.155 ***
Wuxiα0(τ)−4.595 ***−2.872 ***−1.603 ***−1.017 **−0.4080.0980.814 *1.826 ***3.125 ***
ρ1(τ)0.352 ***0.613 ***0.547 ***0.549 ***0.624 ***0.640 ***0.514 ***0.635 ***0.697*
Nantongα0(τ)−6.060 ***−4.289 ***−2.303 ***−0.718−0.0640.4911.391 ***2.378 ***4.946 ***
ρ1(τ)−0.377 ***−0.005 ***0.468 **0.228 ***0.236 ***0.277 ***0.238 ***0.25 ***0.071 ***
Lianyungangα0(τ)−7.285 ***−6.004 ***−4.627 ***−2.676 ***−1.255 *0.0301.917 *3.495 ***6.99 ***
ρ1(τ)0.755 *0.667 ***0.729 ***0.769 ***0.756 ***0.760 ***0.741 ***0.742 ***0.662 ***
Suzhouα0(τ)−3.589 ***−2.442 ***−1.316 ***−0.662 *−0.1400.1060.3761.920 ***2.907 ***
ρ1(τ)0.384 ***0.451 ***0.556 ***0.562 ***0.579 ***0.562 ***0.546 ***0.708 **0.677 *
Huaianα0(τ)−3.018 ***−2.272 ***−1.766 ***−1.480 ***−0.665−0.3541.420 ***2.379 ***3.537 ***
ρ1(τ)0.558 ***0.567 ***0.620 ***0.593 ***0.585 ***0.651 **0.449 ***0.524 ***0.366 ***
Changzhouα0(τ)−4.688 ***−2.891 ***−1.556 ***−1.009 ***−0.5150.1831.017 *2.719 ***4.598 ***
ρ1(τ)0.631*0.505 ***0.495 ***0.475 ***0.507 ***0.480 ***0.472 ***0.324 ***0.388 ***
Tàizhouα0(τ)−4.021 ***−3.096 ***−1.847 ***−1.163 **−0.4240.0941.388 **2.582 ***3.910 ***
ρ1(τ)0.245 ***0.379 ***0.681 ***0.554 ***0.620 ***0.633 ***0.598 ***0.419 ***0.360 ***
Yanchengα0(τ)−2.672 ***−2.325 ***−1.992 ***−1.569 ***−0.202−0.0870.6872.492 ***3.480 ***
ρ1(τ)0.534 ***0.529 ***0.567 ***0.566 ***0.682 ***0.663 ***0.533 ***0.537 ***0.312 ***
Suqianα0(τ)−4.864 ***−3.396 ***−2.068 ***−0.9990.3161.306 ***1.427 ***2.520 ***4.080 ***
ρ1(τ)0.8440.749 **0.593 ***0.644 ***0.617 ***0.629 ***0.621 ***0.494 ***0.735 **
Zhenjiangα0(τ)−4.262 ***−2.242 ***−1.561 ***−0.866 *0.0270.4161.416 ***1.729 ***2.882 ***
ρ1(τ)0.243 ***0.556 **0.526 ***0.425 ***0.34 ***0.303 ***0.448 ***0.404 ***0.242 ***
Yangzhouα0(τ)−4.166 ***−3.631 ***−1.321 *−0.760−0.0800.2450.6712.223 ***4.721 ***
ρ1(τ)0.518 **0.451 ***0.356 ***0.296 ***0.370 ***0.327 ***0.265 ***0.403 ***0.520 *
Xuzhouα0(τ)−6.634 ***−3.091 ***−2.352 ***−2.143 ***−0.7280.3571.476 **2.485 **6.965 ***
ρ1(τ)0.578 ***0.628 ***0.671 ***0.685 ***0.668 ***0.522 ***0.522 ***0.588 ***0.424 ***
Zhoushanα0(τ)−1.796 ***−1.012 ***−0.739 ***−0.476 ***−0.2130.0730.2840.624 **1.491 ***
ρ1(τ)0.197 ***0.355 ***0.361 ***0.511 ***0.563 ***0.644 ***0.615 ***0.561 ***0.423 ***
Lishuiα0(τ)−2.948 ***−1.551 ***−0.936 ***−0.535 **−0.305 **−0.0700.4941.129 ***1.973 ***
ρ1(τ)0.065 ***0.375 ***0.487 ***0.510 ***0.493 ***0.417 ***0.458 ***0.542 ***0.403 ***
Ningboα0(τ)−2.599 ***−1.508 ***−1.299 ***−0.850 ***−0.3240.1890.636 **1.030 ***2.261 ***
ρ1(τ)0.315 ***0.573 ***0.538 ***0.598 ***0.578 ***0.449 ***0.472 ***0.505 ***0.806
Wenzhouα0(τ)−4.207 ***−1.840 ***−1.056 **−0.755 *−0.1480.1830.867 ***1.557 ***3.457 ***
ρ1(τ)−0.180 ***0.175 ***0.374 ***0.430 ***0.355 ***0.314 ***0.197 ***0.05 ***−0.262 ***
Quzhouα0(τ)−4.441 ***−3.018 ***−2.391 ***−1.338 ***−0.993 **0.6191.638 ***2.631 ***3.276 ***
ρ1(τ)−0.287 ***0.077 ***−0.027 ***0.420 ***0.409 ***0.069 ***0.212 ***0.303 ***0.280 ***
Jiaxingα0(τ)−4.853 ***−3.111 ***−1.468 **−0.756 *−0.1980.3551.175 ***1.621 ***3.022 ***
ρ1(τ)0.259 ***0.585 ***0.572 ***0.602 ***0.592 ***0.596 ***0.583 ***0.637 ***0.809
Jinhuaα0(τ)−5.475 ***−2.969 ***−2.266 ***−0.928 **−0.4940.1341.088 *2.420 ***4.250 ***
ρ1(τ)−0.15 ***0.095 ***0.303 ***0.432 ***0.414 ***0.394 ***0.573 ***0.578 ***0.947
Huzhouα0(τ)−4.768 ***−3.213 ***−2.02 ***−1.215 ***−0.5580.1641.200 **2.002 ***3.615 ***
ρ1(τ)0.499 **0.533 ***0.538 ***0.428 ***0.327 ***0.363 ***0.477 ***0.41 ***0.223 ***
Hangzhouα0(τ)−2.584 ***−1.997 ***−0.833 **−0.597 **−0.2370.1730.876 ***1.485 ***2.256 ***
ρ1(τ)0.336 ***0.457 ***0.573 ***0.586 ***0.607 ***0.652 ***0.693 ***0.566 ***0.694 *
Shaoxingα0(τ)−4.377 ***−3.228 ***−1.99 ***−1.32 **−0.3860.4471.128 ***1.447 ***3.420 ***
ρ1(τ)0.355 ***0.401 ***0.512 ***0.576 ***0.52 ***0.277 ***0.345 ***0.333 ***0.346 ***
Táizhouα0(τ)−1.971 ***−1.520 ***−0.923 ***−0.652 ***−0.2610.0930.421 *1.135 ***1.792 ***
ρ1(τ)0.278 ***0.443 ***0.691 ***0.656 ***0.599 ***0.544 ***0.516 ***0.594 ***0.616 **
Pearl River Delta
Guangzhouα0(τ)−2.567 ***−1.541 ***−1.069 ***−0.649 ***−0.1630.0580.451 *1.213 ***2.022 ***
ρ1(τ)0.167 ***0.109 ***0.178 ***0.095 ***0.132 ***0.072 ***0.111 ***0.212 ***0.367 ***
Shenzhenα0(τ)−0.951 ***−0.799 ***−0.613 ***−0.439 ***−0.343 **0.0970.312 *0.564 ***1.095 ***
ρ1(τ)0.474 ***0.425 ***0.383 ***0.479 ***0.470 ***0.576 ***0.574 ***0.501 ***1.013
Zhuhaiα0(τ)−1.802 ***−1.384 ***−0.967 ***−0.391 *0.0640.408 **0.523 **1.271 ***1.758 ***
ρ1(τ)0.521 ***0.481 ***0.488 ***0.517 ***0.561 ***0.648 ***0.659 ***0.672 **0.451 ***
Foshanα0(τ)−3.109 ***−2.610 ***−1.565 ***−0.505−0.272−0.0090.942 **1.237 ***2.418 ***
ρ1(τ)−0.129 ***0.016 ***0.049 ***0.279 ***0.359 ***0.326 ***0.419 ***0.427 ***0.563 *
Jiangmenα0(τ)−3.535 ***−2.809 ***−2.096 ***−1.145 **−0.747 *0.5250.954 **1.594 **3.832 ***
ρ1(τ)0.473 ***0.492 ***0.514 ***0.484 ***0.436 ***0.411 ***0.458 ***0.532 ***0.73 *
Dongguanα0(τ)−2.202 ***−1.819 ***−1.338 ***−0.728 **−0.1750.1460.682 **1.207 ***1.739 ***
ρ1(τ)0.420 ***0.409 ***0.493 ***0.590 ***0.501 ***0.453 ***0.307 ***0.244 ***0.363 ***
Zhongshanα0(τ)−2.599 ***−1.953 ***−1.261 ***−0.716 **−0.2890.1660.910 **1.309 ***2.372 ***
ρ1(τ)0.587 **0.593 **0.651 **0.616 ***0.621 ***0.551 ***0.654 **0.654 **0.433 **
Huizhouα0(τ)−1.593 ***−1.020 ***−0.701 ***−0.405 **−0.0340.227 **0.333 ***0.477 ***1.412 ***
ρ1(τ)0.353 ***0.637 ***0.546 ***0.552 ***0.538 ***0.521 ***0.527 ***0.502 ***0.857
Zhaoqingα0(τ)−3.788 ***−1.768 ***−1.518 ***−0.709 *0.2720.514 *0.923 ***1.671 ***3.435 ***
ρ1(τ)0.489 ***0.306 ***0.263 ***−0.132 ***0.033 ***−0.005 ***−0.064 ***−0.179 ***0.339 ***
Municipality
Chongqingα0(τ)−2.920 ***−1.251 ***−0.954 ***−0.607 ***−0.3910.1070.497 *1.062 **2.354 ***
ρ1(τ)0.342 ***0.593 ***0.636 ***0.606 ***0.657 ***0.509 ***0.489 ***0.545 ***0.511 **
Provincial capital cities
Lhasaα0(τ)−1.367 ***−1.073 ***−0.861 ***−0.406 **−0.276 *0.1600.626 **1.340 ***1.722 ***
ρ1(τ)0.699 **0.632 **0.684 **0.435 ***0.385 ***0.313 ***0.158 ***0.295 ***0.286 ***
Haikouα0(τ)−1.220 ***−0.957 ***−0.774 ***−0.450−0.1980.1880.419 **0.626 ***1.133 ***
ρ1(τ)−0.017 ***−0.048 ***−0.027 ***0.168 ***0.289 ***0.505 **0.417 ***0.452 ***0.977
Kunmingα0(τ)−2.620 ***−2.259 ***−1.502 ***−0.4560.0020.4741.548 ***1.772 ***2.606 ***
ρ1(τ)0.5790.505 **0.275 ***−0.050 ***−0.079 ***0.039 ***0.016 ***0.054 ***−0.119 ***
Fuzhouα0(τ)−1.244 ***−0.960 ***−0.634 ***−0.198−0.0830.0940.349 ***0.679 ***1.151 ***
ρ1(τ)0.274 ***0.401 ***0.354 ***0.515 ***0.659 ***0.592 ***0.520 ***0.575 ***0.982
Chengduα0(τ)−2.367 ***−1.587 ***−1.180 ***−0.786 **−0.1130.386 *0.638 ***1.074 **2.254 ***
ρ1(τ)0.359 ***0.446 ***0.400 ***0.459 ***0.429 ***0.419 ***0.385 ***0.326 ***0.183 ***
Xi’anα0(τ)−5.653 ***−4.220 ***−1.762 **−1.169 ***−0.679 *0.1511.040 *2.106 ***3.409 ***
ρ1(τ)0.599 **0.612 ***0.512 ***0.542 ***0.523 ***0.470 ***0.388 ***0.302 ***0.355 ***
Guiyangα0(τ)−5.629 ***−3.194 ***−2.275 ***−1.467 **−0.2230.4191.392 ***2.139 ***2.955 ***
ρ1(τ)0.407 ***0.566 ***0.538 ***0.510 ***0.470 ***0.432 ***0.448 ***0.381 ***0.348 ***
Urumqiα0(τ)−4.522 ***−3.654 ***−1.325 **−1.101 **−0.3880.4620.914 ***1.648 ***2.256
ρ1(τ)0.474 ***0.417 ***0.550 ***0.533 ***0.621 ***0.552 ***0.578 ***0.744 ***0.744
Nanningα0(τ)−2.283 ***−1.711 ***−1.41 ***−0.986 ***−0.446 *−0.0050.4741.1852.615 ***
ρ1(τ)0.488 ***0.450 ***0.495 ***0.509 ***0.440 ***0.413 ***0.366 ***0.296 ***0.157 ***
Xiningα0(τ)−7.846 ***−5.151 ***−3.787 ***−2.762 ***−0.9480.9312.837 ***4.156 ***6.816 ***
ρ1(τ)0.541 ***0.433 ***0.396 ***0.479 ***0.542 ***0.613 ***0.625 ***0.538 ***0.385 ***
Shenyangα0(τ)−14.438 ***−10.199 ***−7.654 ***−6.133 ***−4.659 ***−1.7741.8426.919 ***14.989 ***
ρ1(τ)0.678 ***0.641 ***0.619 ***0.603 ***0.609 ***0.579 ***0.560 ***0.506 ***0.479 ***
Changchunα0(τ)−10.151 ***−7.356 ***−4.499 ***−2.756 ***−2.278 **−0.0012.339*4.384 ***8.613 ***
ρ1(τ)0.669 ***0.659 ***0.593 ***0.592 ***0.577 ***0.461 ***0.484 ***0.480 ***0.486 ***
Harbinα0(τ)−11.739 ***−9.598 ***−4.847 ***−3.252 **−0.850.2443.349 ***4.392 ***5.679 ***
ρ1(τ)0.470 ***0.507 ***0.538 ***0.589 ***0.510 ***0.496 ***0.465 ***0.415 ***0.329 ***
Wuhanα0(τ)−3.777 ***−2.607 ***−1.693 ***−1.141 **−0.1750.5790.982 **2.327 ***2.983 ***
ρ1(τ)0.430 ***0.676 ***0.583 ***0.495 ***0.635 ***0.640 ***0.638 ***0.764 *0.813 *
Yinchuanα0(τ)−22.699 ***−14.753 ***−10.862 ***−8.093 ***−3.529−1.0627.217 **10.43 ***18.557 ***
ρ1(τ)0.446 ***0.662 ***0.652 ***0.626 ***0.603 ***0.528 ***0.579 ***0.594 ***0.53 ***
Lanzhouα0(τ)−5.699 ***−3.987 ***−2.576 ***−1.838 ***−1.227 **−0.0161.702 **2.886 ***5.694 ***
ρ1(τ)0.250 ***0.167 ***0.061 ***0.127 ***0.130 ***0.305 ***0.303 ***0.359 ***0.178 ***
Nanchangα0(τ)−4.104 ***−2.712 ***−1.334 ***−0.65*−0.2060.2710.5701.822 ***3.115 ***
ρ1(τ)−0.127 ***−0.101 ***−0.086 ***−0.115 ***−0.13 ***−0.102 ***−0.077 ***−0.243 ***−0.569 ***
Hefeiα0(τ)−3.387 ***−1.632 ***−0.888 **−0.426−0.0660.4571.229 ***1.573 ***2.593 ***
ρ1(τ)0.544 **0.613 ***0.644 ***0.67 ***0.696 ***0.718 ***0.813 *0.827 *0.763 **
Taiyuanα0(τ)−19.09 ***−16.517 ***−13.061 ***−7.071 ***−5.128 **0.3472.5528.722 **24.697 ***
ρ1(τ)0.580 ***0.566 ***0.617 ***0.618 ***0.638 ***0.570 ***0.615 ***0.562 ***0.111 ***
Ji’nanα0(τ)−7.89 ***−4.715 ***−3.252 ***−1.155 *−0.3660.5490.986 *1.2928.675 ***
ρ1(τ)0.783 **0.652 ***0.623 ***0.681 ***0.658 ***0.691 ***0.678 ***0.675 ***0.895
Hohhotα0(τ)−7.422 ***−5.408 ***−4.193 ***−3.018 ***−2.014 **−0.0252.221 **3.253 ***5.526 ***
ρ1(τ)0.554 ***0.516 ***0.465 ***0.508 ***0.538 ***0.535 ***0.61 ***0.577 ***0.473 ***
Changshaα0(τ)−2.982 ***−1.882 ***−1.310 ***−0.860 **−0.280.2280.897 **1.632 ***2.561 ***
ρ1(τ)0.404 ***0.318 ***0.265 ***0.362 ***0.317 ***0.249 ***0.495 ***0.434 ***0.360 ***
Zhengzhouα0(τ)−6.349 ***−4.166 ***−2.445 ***−1.888 **0.3530.911 *1.653 ***2.407 ***4.265 ***
ρ1(τ)0.695 ***0.685 ***0.683 ***0.856 **0.801 **0.855*0.760 ***0.729 ***0.752 **
Cities under separate state planning
Xiamenα0(τ)−6.349 ***−4.166 ***−2.445 ***−1.888 **0.3530.911*1.653 ***2.407 ***4.265 ***
ρ1(τ)0.148 ***0.217 ***0.139 ***0.197 ***0.337 ***0.154 ***0.243 ***0.166 ***−0.146 ***
Dalianα0(τ)−7.389 ***−4.762 ***−3.17 ***−2.273 ***−1.375 **−0.7021.4074.213 ***5.294 **
ρ1(τ)0.633 ***0.621 ***0.714 ***0.698 ***0.752 ***0.722 ***0.810 **0.693 ***0.670 *
Qingdaoα0(τ)−4.180 ***−2.814 ***−2.442 ***−1.638 ***−0.649 *−0.2390.4422.310 ***3.529 **
ρ1(τ)0.625 ***0.645 ***0.644 ***0.647 ***0.649 ***0.686 ***0.752 ***0.839*0.923
***, **, and * indicate significance at the 1%, 5%, and 10% levels.
Table 6. Results for Fourier quantile unit root at particular quantiles.
Table 6. Results for Fourier quantile unit root at particular quantiles.
Quantiles0.10.20.30.40.50.60.70.80.9
Beijing–Tianjin–Hebei Urban Agglomeration
Beijing−4.431 ***−3.271 ***−3.737 ***−3.477 ***−3.494 ***−2.993 **−3.804 ***−2.275 *−2.192 *
Tianjin−1.792−2.176*−3.486 ***−4.336 ***−4.995 ***−6.152 ***−6.389 ***−2.676 **−3.072 **
Zhangjiakou−1.455−1.751−2.334 **−3.408 ***−2.264 **−1.657−1.554−1.166−1.815
Chengde−2.104−2.284 *−3.933 ***−3.476 ***−3.858 ***−4.152 ***−4.027 ***−3.220 ***−3.189 **
Qinhuangdao−3.920 ***−2.846 **−3.838 ***−3.665 ***−3.707 ***−3.424 ***−2.369 *−2.065−1.660
Langfang−3.039 **−2.926 **−3.041 **−3.258 ***−3.405 ***−3.639 ***−3.253 ***−1.656−1.343
Cangzhou−2.384 *−2.306 *−3.000 **−4.256 ***−3.755 ***−4.407 ***−3.947 ***−4.273 ***−2.579 *
Tangshan−1.963−2.405 *−2.143 *−2.209 *−3.126 ***−3.237 ***−2.711 **−1.957−2.401 *
Handan−2.615 **−3.395 ***−2.633 **−3.772 ***−4.095 ***−4.523 ***−3.969 ***−1.637−1.121
Hengshui−5.771 ***−5.049 ***−5.887 ***−4.646 ***−4.101 ***−3.560 ***−2.051 *−1.504−0.967
Xingtai−1.657−3.367 ***−2.715 **−3.198 ***−4.359 ***−5.634 ***−5.998 ***−6.918 ***−1.025
Shijiazhuang−2.594 *−2.357 **−2.958 **−3.489 ***−3.414 ***−3.707 ***−2.735 **−2.552 **−1.798
Baoding−4.118 ***−2.980 **−3.093 **−2.277 *−2.717 **−2.956 **−2.563 **−2.217 *−2.490 *
Yangtze River Delta
Shanghai−1.588−3.965 ***−4.744 ***−4.276 ***−2.743 **−2.164 *−3.032 **−3.135 ***−1.012
Nanjing−1.974−2.877 **−3.037 **−2.934 **−2.100 *−0.932−1.864 *−1.660−2.972 *
Wuxi−3.113 **−2.655 **−3.104 ***−4.017 ***−3.148 ***−2.857 **−3.066 ***−1.095−0.798
Nantong−5.365 ***−2.858 **−2.065 *−2.906 **−3.299 ***−4.075 ***−3.517 ***−2.597 **−3.406 **
Lianyungang−1.149−2.553 **−2.318 **−3.214 ***−3.126 **−2.559 **−2.681 **−2.517 **−0.971
Suzhou−2.771 **−4.425 ***−3.179 ***−3.855 ***−3.261 ***−3.925 ***−3.157 ***−1.611−0.753
Huaian−1.963*−2.950 **−3.445 ***−3.023 **−2.874 **−2.355*−3.171 ***−2.337 *−1.993
Changzhou−1.025−1.999 *−3.275 ***−4.541 ***−4.067 ***−3.817 ***−3.531 ***−2.461 **−1.565
Tàizhou−2.000 *−4.374 ***−2.504 **−3.346 ***−2.514 **−2.468 **−2.467 *−3.473 ***−2.671 **
Yancheng−1.584−4.341 ***−4.339 ***−3.880 ***−2.110−2.052−2.824 **−2.519 **−3.229 **
Suqian−0.815−1.684−2.934 **−2.824 **−3.120 **−2.914 **−3.679 ***−5.163 ***−1.267
Zhenjiang−3.759 ***−2.375*−2.490 **−2.910 **−4.897 ***−4.397 ***−4.157 ***−2.663 **−2.858 **
Yangzhou−1.890−2.366 *−3.512 ***−4.960 ***−3.608 ***−3.986 ***−4.269 ***−2.519 **−1.030
Xuzhou−2.108 *−2.414 **−2.660 **−2.473 **−3.735 ***−5.057 ***−3.128 ***−1.731−2.120
Zhoushan−3.741 ***−4.082 ***−4.145 ***−4.317 ***−4.024 ***−2.339 **−2.245 *−1.810−1.553
Lishui−3.595 **−2.676 **−2.633 **−3.100 ***−3.990 ***−4.136 ***−3.545 ***−2.350 **−2.996 **
Ningbo−3.736 **−2.980 **−2.968 **−2.778 **−3.871 ***−3.965 ***−3.569 ***−2.639 **−0.600
Wenzhou−3.109 **−2.816 **−2.877 **−3.226 ***−4.668 ***−4.047 ***−4.364 ***−3.640 ***−4.837 ***
Quzhou−6.526 ***−4.819 ***−5.251 ***−2.836 **−2.880 **−4.890 ***−3.953 ***−3.393 ***−3.729 ***
Jiaxing−2.343 *−1.994 *−3.003 ***−3.625 ***−3.658 ***−3.793 ***−4.052 ***−3.245 ***−0.396
Jinhua−5.078 ***−4.466 ***−4.264 ***−3.421 ***−4.130 ***−4.673 ***−2.516 **−2.201 *−0.271
Huzhou−2.374 *−2.002 *−2.299 *−3.767 ***−4.522 ***−4.303 ***−3.126 ***−1.984−1.850
Hangzhou−2.874 **−3.534 ***−3.045 **−2.861 **−2.741 **−2.782 **−2.308 *−3.116 **−1.498
Shaoxing−2.747 **−2.561 **−2.134 **−3.079 **−3.858 ***−4.830 ***−4.303 ***−2.959 **−1.367
Táizhou−3.735 ***−3.460 ***−2.353 **−2.617 **−3.202 ***−2.964 **−3.126 **−1.606−1.550
Pearl River Delta
Guangzhou−3.012 **−3.991 ***−4.614 ***−5.417 ***−5.302 ***−6.002 ***−4.105 ***−3.747 ***−1.347
Shenzhen−2.153 *−4.308 ***−7.644 ***−4.326 ***−4.198 ***−2.425 **−2.448 *−2.233 *0.031
Zhuhai−1.551−2.860 **−3.493 ***−3.533 ***−3.529 ***−2.295 *−2.443 *−2.512 **−4.548 ***
Foshan−4.444 ***−4.717 ***−5.051 ***−4.960 ***−4.104 ***−4.972 ***−3.125 ***−1.530−0.987
Jiangmen−3.066 **−3.658 ***−3.357 ***−3.963 ***−4.340 ***−3.716 ***−3.093 ***−2.177 *−1.088
Dongguan−1.788−2.838 **−2.968 **−2.966 **−3.170 ***−3.414 ***−3.953 ***−2.975 **−2.438*
Zhongshan−2.165−1.994−2.015−2.298 *−2.241 *−2.678 **−1.802−1.648−3.293 **
Huizhou−3.648 ***−2.158 *−2.849 **−3.160 ***−4.163 ***−3.686 ***−3.491 ***−3.225 ***−0.379
Zhaoqing−1.276−2.482 **−3.125 ***−4.606 ***−6.966 ***−4.288 ***−4.505 ***−3.274 **−3.464 **
Municipality
Chongqing−2.029−1.866 *−1.912 *−3.249 ***−3.300 ***−3.290 ***−3.261 ***−1.977−1.487
Provincial capital cities
Lhasa−1.646 *−2.143 **−1.985 **−3.441 ***−2.940 **−3.140 ***−3.740 ***−3.868 ***5.554 ***
Haikou−8.688 ***−6.417 ***−5.434 ***−3.652 ***−3.344 ***−2.061−2.331*−1.692−0.057
Kunming−0.943−1.769−3.106 **−4.492 ***−4.126 ***−3.907 ***−4.379 ***−3.282 ***−4.063 ***
Fuzhou−5.844 ***−3.291 ***−4.203 ***−3.038 **−2.383 **−2.482 **−2.200−1.659−0.059
Chengdu−3.663 ***−3.027 **−3.287 ***−3.365 ***−4.284 ***−3.181 ***−2.939 **−3.136 **−2.692 **
Xi’an−1.853−2.434 **−4.252 ***−4.296 ***−4.063 ***−6.378 ***−4.492 ***−5.009 ***−1.820
Guiyang−3.091 **−3.001 **−3.064 ***−3.649 ***−4.248 ***−4.476 ***−4.717 ***−3.549 ***−2.650*
Urumqi−5.919 ***−4.259 ***−3.797 ***−3.787 ***−3.237 ***−5.535 ***−4.766 ***−0.943−0.601
Nanning−2.634 **−3.875 ***−4.503 ***−4.073 ***−4.140 ***−3.673 ***−3.105 **−3.200 ***−2.547 **
Xining−2.878 **−4.882 ***−4.406 ***−3.397 ***−3.187 ***−2.407 *−1.887−2.747 **−2.175 *
Shenyang−3.119 **−3.920 ***−5.963 ***−7.012 ***−5.409 ***−5.491 ***−4.107 ***−4.181 ***−2.139 ***
Changchun−4.222 ***−3.361 ***−4.842 ***−4.664 ***−4.372 ***−5.844 ***−4.496 ***−3.142 **−3.766 **
Harbin−7.021 ***−6.353 ***−5.377 ***−4.720 ***−5.555 ***−7.051 ***−8.833 ***−11.338 ***−5.662 ***
Wuhan−3.593 **−2.346 **−2.881 ***−3.556 ***−2.372 *−2.495−2.294−1.585−1.578
Yinchuan−5.568 ***−2.886 **−2.875 **−3.002 **−3.638 ***−3.667 ***−3.143 **−2.388 **−2.606*
Lanzhou−5.398 ***−5.548 ***−6.695 ***−6.973 ***−5.670 ***−4.467 ***−3.204 ***−2.148 *−2.225 *
Nanchang−3.521 **−3.788 ***−4.615 ***−6.417 ***−5.654 ***−6.634 ***−4.305 ***−4.119 ***−4.239 ***
Hefei−1.632−1.905−2.559 **−2.273 *−2.948 **−2.405 *−1.679−1.775−1.592
Taiyuan−4.166 ***−3.798 ***−3.220 ***−4.176 ***−3.324 ***−3.477 ***−2.141 *−1.837−1.937
Jinan−1.435−2.470 **−4.412 ***−4.649 ***−6.768 ***−3.981 ***−3.028 **−2.074 *−0.396
Hohhot−5.727 ***−6.320 ***−6.228 ***−5.221 ***−5.281 ***−4.474 ***−3.432 ***−4.263 ***−2.286 **
Changsha−2.546 *−3.002 **−5.320 ***−4.225 ***−4.666 ***−4.595 ***−3.294 ***−3.250 ***−3.203 **
Zhengzhou−3.442 **−3.725 ***−3.505 ***−1.465−2.012 *−1.427−2.846 **−3.219 ***−1.640
Cities under separate state planning
Xiamen−3.660 **−3.487 ***−4.534 ***−4.052 ***−3.154 ***−4.103 ***−3.161 ***−2.633 **−3.782 **
Dalian−3.723 **−3.632 ***−3.379 ***−3.240 ***−2.766 **−3.140 **−2.062 *−1.605−0.965
Qingdao−3.765 **−6.113 ***−5.171 ***−5.241 ***−5.403 ***−3.669 ***−2.444 **−1.016−0.142
***, **, and * indicate significance at the 1%, 5%, and 10% levels.
Table 7. Results for overall the Fourier quantile unit root test.
Table 7. Results for overall the Fourier quantile unit root test.
CityFQKS statisticCV 10%CV 5%CV 1%Optimal FrequencyOptimal F-statistic
Beijing–Tianjin–Hebei Urban Agglomeration
Beijing4.431 **3.2573.6914.8424.721.229
Tianjin6.389 ***3.2473.6794.8050.118.423
Zhangjiakou3.408 *3.0763.5674.6700.128.983
Chengde4.152 **3.3323.7224.9884.824.667
Qinhuangdao3.920 **3.3023.6814.7714.833.107
Langfang3.639 **3.2713.6384.6834.722.617
Cangzhou4.407 **3.3113.7234.8384.728.438
Tangshan3.237 *3.2233.6594.7040.122.792
Handan4.523 ***2.9843.4004.3054.719.901
Hengshui5.887 ***3.1603.6824.9590.121.210
Xingtai6.918 ***3.2643.6854.9034.817.478
Shijiazhuang3.707 **3.1883.6024.6024.724.438
Baoding4.118 **3.2633.6694.7850.119.500
Yangtze River Delta
Shanghai4.744 ***3.2453.6394.7030.129.304
Nanjing3.0373.1863.5904.6980.135.179
Wuxi4.017 **3.2273.6804.9130.153.096
Nantong5.365 ***3.2223.6224.8260.160.307
Lianyungang3.214 *3.0903.4634.3174.728.714
Suzhou4.425 **3.2103.6784.7440.147.641
Huaian3.445 *3.2153.6314.7300.145.100
Changzhou4.541 ***3.0873.4754.3200.159.147
Tàizhou4.374 ***3.1213.4844.3080.171.957
Yancheng4.341 **3.2893.6634.7880.140.810
Suqian5.163 ***3.3513.7674.9710.122.540
Zhenjiang4.897 ***3.1343.4834.2910.5122.378
Yangzhou4.960 ***3.2373.6624.6920.148.814
Xuzhou5.057 ***3.2363.6434.7380.147.864
Zhoushan4.317 **3.1363.5184.3830.117.966
Lishui4.136 **3.1033.4894.8500.114.357
Ningbo3.965 **3.2533.6864.9290.128.160
Wenzhou4.837 ***3.2293.6524.7930.120.234
Quzhou6.526 ***3.1873.6504.7980.140.906
Jiaxing4.052 **3.0983.6114.9330.129.824
Jinhua5.078 ***3.1633.5834.8420.135.089
Huzhou4.522 **3.2423.6754.8464.711.500
Hangzhou3.534 *3.2593.7374.9240.127.361
Shaoxing4.830 **3.1743.6464.9840.158.567
Táizhou3.735 **3.1313.4954.3954.818.381
Pearl River Delta
Guangzhou6.002 ***3.2103.6334.6690.119.932
Shenzhen7.644 ***3.2283.5744.4490.124.866
Zhuhai4.548 ***3.1143.4164.1784.740.026
Foshan5.051 ***3.2943.7494.8700.124.451
Jiangmen4.340 **3.3053.7074.8840.119.891
Dongguan3.953 **3.2403.7034.8224.720.705
Zhongshan3.293 *3.2313.6214.6994.825.200
Huizhou4.163 **3.1463.4614.3224.710.425
Zhaoqing6.966 ***3.2443.6574.7550.139.582
Municipality
Chongqing3.300 *3.1983.6324.8690.132.420
Provincial capital cities
Lhasa5.554 ***3.0663.4134.2750.123.913
Haikou8.688 ***3.2243.6734.6094.68.313
Kunming3.593 **3.1833.5774.7090.135.461
Fuzhou5.844 ***3.1133.5144.3881.59.918
Chengdu4.284 **3.1623.4914.3120.142.587
Xi’an6.378 ***3.2733.6994.8864.840.747
Guiyang4.717 **3.2593.6204.8554.838.941
Urumqi5.919 ***3.1833.6374.8234.819.204
Nanning4.503 ***3.1403.5174.3984.825.52
Xining4.882 ***3.2573.6704.7264.713.739
Shenyang7.012 ***3.2523.6955.0564.822.382
Changchun5.844 ***3.2603.6624.7624.742.081
Harbin11.338 ***3.2323.7055.1264.727.913
Wuhan3.593 **3.1833.5774.7090.135.461
Yinchuan5.568 ***3.2233.6124.8564.829.116
Lanzhou6.973 ***3.2643.6644.7344.770.174
Nanchang6.634 ***3.2823.6834.6300.143.119
Hefei2.9483.2553.6594.8500.137.914
Taiyuan4.176 **3.2223.7034.7934.732.412
Jinan6.768 ***3.3483.7674.8300.130.034
Hohhot6.320 ***3.2673.7254.9874.844.381
Changsha5.320 ***3.2603.6594.7820.152.570
Zhengzhou3.725 **3.2383.6724.6154.821.613
Cities under separate state planning
Xiamen4.534 **3.2663.6424.7370.116.879
Dalian3.723 **3.2643.7224.8564.825.841
Qingdao6.113 ***3.1363.5994.8070.131.960
***, **, and * indicate significance at the 1%, 5%, and 10% significant levels. The critical value for QKS statistics generated by bootstrap techniques with 5000 iterations.

Share and Cite

MDPI and ACS Style

Zhang, Y.-C.; Si, D.-K.; Zhao, B. The Convergence of Sulphur Dioxide (SO2) Emissions Per Capita in China. Sustainability 2020, 12, 1781. https://doi.org/10.3390/su12051781

AMA Style

Zhang Y-C, Si D-K, Zhao B. The Convergence of Sulphur Dioxide (SO2) Emissions Per Capita in China. Sustainability. 2020; 12(5):1781. https://doi.org/10.3390/su12051781

Chicago/Turabian Style

Zhang, Yu-Chen, Deng-Kui Si, and Bing Zhao. 2020. "The Convergence of Sulphur Dioxide (SO2) Emissions Per Capita in China" Sustainability 12, no. 5: 1781. https://doi.org/10.3390/su12051781

APA Style

Zhang, Y. -C., Si, D. -K., & Zhao, B. (2020). The Convergence of Sulphur Dioxide (SO2) Emissions Per Capita in China. Sustainability, 12(5), 1781. https://doi.org/10.3390/su12051781

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

Article Metrics

Back to TopTop