Next Article in Journal
Automatic Detection of Electric Field VLF Electromagnetic Wave Abnormal Disturbance on Zhangheng-1 Satellite
Next Article in Special Issue
The Hot Topics, Frontiers and Trends about Research on the Relationship between Air Pollution and Public Health—Visual Analysis Based on Knowledge Map
Previous Article in Journal / Special Issue
Simulation of Summer Rainfall in Thailand by IAP-AGCM4.1
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Investigating Spatial Heterogeneity of the Environmental Kuznets Curve for Haze Pollution in China

1
Institute of Tropical Agriculture & Food Security (ITAFoS), School of Business and Economics (SPE), Universiti Putra Malaysia, Serdang 43400, Malaysia
2
Department of Public Administration, Ewha Womans University, Seoul 03760, Korea
3
Faculty of Economics and Management, Universiti Putra Malaysia, Serdang 43400, Malaysia
*
Author to whom correspondence should be addressed.
Submission received: 18 April 2022 / Revised: 9 May 2022 / Accepted: 10 May 2022 / Published: 14 May 2022
(This article belongs to the Special Issue Advanced Climate Simulation and Observation)

Abstract

:
This study investigates the environmental Kuznets curve (EKC) for haze in 31 cities and provinces across China using the spatial data for a period of 15 years, from 2000 to 2014. We utilized the geographically weighted regression (GWR) model to consider the spatial non-stationary characteristics of the air quality in a vast territory. This approach allowed us to verify the region-specific characteristics, while the global model estimated the average relationship across the entire nation. Although the EKC for haze was statistically significant in the global models, the results only confirmed the existence of an EKC between the overall air quality and economic performance. Thus, it was difficult to determine the regional differences in an EKC. The results of the GWR model found the spatial variability of each variable and showed significant spatial heterogeneity in the EKC across regions. Although six regions—Beijing, Gansu, Heilongjiang, Jiangxi, Jilin, Liaoning, Shanghai, Tianjin, Xinjiang, and Zhejiang—showed inverted U-shaped EKCs, these were only statistically significant in three big cities—Beijing, Tianjin, and Shanghai. The results demonstrated no EKCs in the other 25 provinces and cities. These results provide strong empirical evidence that there is significant spatial heterogeneity in the EKC of China. Thus, a more regionally specialized air pollution control policy is required to create an effective policy for balanced economic growth in China.

1. Introduction

China has emerged as one of the largest and most robust economies in the world. The ruling Chinese Communist Party made the decision to open up the country for trade and investment, as well as adopting a limited number of capitalist strategies. This is a natural progression toward a better model for regional trade and local growth [1]; however, regional discrepancies in sustainable development in China are among the barriers to future success [2].
The rapid industrialization in China has caused side effects, with the leading issue being environmental challenges. Currently, the most prominent environmental issue is air pollution. The haze phenomenon has significantly and negatively affected economic performance and contributes to environmental pollution in China. Haze is defined as the introduction of particulate matter (PM), biological molecules, and other contaminants into the atmosphere. The presence of excess suspended matter in the atmosphere can disrupt the daily life of humans. These contaminants usually negatively affect the health of humans and other living organisms. Furthermore, they damage crops, natural environments, and built environments [3,4]. It can also affect the economic performance of a region by causing declining employee performance, haze-related illnesses, slower crop growth, low visibility, and other issues [4].
In recent years, large cities in China have frequently experienced severe fog and haze episodes that disrupt visibility. There have been rapid changes in the patterns of haze and fog occurrence in these areas over the past several decades. For example, in research about the occurrence of haze and fog in the North China Plain, the occurrence of haze and fog was quite low from the 1950s to the 1980s but reached a peak from 1981 to 1998 [4,5]. Since 2000, haze pollution has shown a rapid increase; one of the reasons for this is that China’s policy of opening up stimulates the evolution of industries. The haze problem became more noticeable in the mid-2000s, when environmentalists observed that the haze was not associated with the changing of the seasons. Instead, these particles remained in the air perpetually. In 2013, the Chinese Academy of Sciences documented the occurrence of a severe episode of haze, which affected northern and eastern China.
One of the most affected places was Beijing, which is the administrative, economic, and cultural center of China, meaning that the negative effects there far outweighed those in other regions. Beijing’s high population density means that the people affected outnumber those affected anywhere else in the country. In Beijing, air pollution index (API) values of 900 have been recorded, making the environment toxic to its local residents (an API of 0–50 is considered good, 50–100 is moderate, 101–200 is unhealthy, 201–300 is very unhealthy, and anything above 300 is classified as hazardous). In terms of economic losses, Beijing lost approximately 254 million USD, which was about 0.008% of the total gross domestic product (GDP) of Beijing in 2013 [6]. Previous studies have examined the effects of haze on China’s economy and most found that haze reduces economic performance as well as human resources [7,8,9,10].
Previous studies have tried to investigate the existence of an environmental Kuznets curve (EKC) in China to show the relationship between economic development and environmental performance and understand the current pollution status [11,12,13,14,15,16,17,18,19]. Scholars have analyzed the EKC for haze using data regarding the concentrations of PM, carbon dioxide (CO2), sulfur dioxide (SO2), and nitrous oxide in the atmosphere [20,21,22,23]. (PM is classified into different sizes, the most commonly measured of which are PM10 and PM2.5. There has been a significant debate between researchers as to which is the better measurement system, but measuring PM2.5 appears to be the preferred method in China as it has been found to be the dominant polluting agent). To investigate the heterogeneity of the EKC across regions in China, some studies have considered the spatial differences in the EKC and found significant spatial heterogeneity [16,20,21,22,23,24,25]. However, most studies focused on a certain pollutant, such as SO2 or PM; thus, it was difficult to discover the EKC between overall air quality and economic performance.
Scholars have also attempted to identify the determinants of haze [7,8,9,10,26,27]. The energy structure is one of the crucial factors affecting haze. Scholars have suggested that the heavy usage of coal, renewable energy, and energy consumption have an influence on haze; however, the empirical evidence is still controversial [27,28,29]. Other socio-economic factors, such as income, human capital, trade openness, transportation, and population density, have also been considered as determinants of haze [16,17,18,19]. Moreover, if we widen the perspective, we should not overlook the impact of meteorology and stable atmospheric conditions on atmospheric pollution. For example, the amount of rainfall and temperature have been regarded as factors that influence the haze problem [25].
This study seeks to determine which economic factors influence haze and which regions have an EKC and turning point. Considering the spatial heterogeneity of economic growth and the haze phenomenon, this study applies the geographically weighted regression (GWR) model to investigate the EKC between haze and economic determinants for 31 Chinese cities. The overall air-quality data API and factors including the gross regional product (GRP), trade openness, and population density are used for analysis. After determining that the EKC has an inverted U, we verify and calculate the turning point, i.e., the decline in pollutants in line with increased economic performance [30]. The results from the GWR approach, which is suitable for considering the spatial non-stationary characteristics in vast territories, can provide a better understanding of the spatial heterogeneity of the haze problem in provinces in China with visual evidence. Moreover, the estimation of the turning points of the EKCs using the GWR estimations can be used to identify trends in sustainable development after the analysis periods. Through these processes, this study determines whether any stronger measures should be undertaken by the Chinese government to curb the haze issue, as it may affect the livelihood of the nation.

2. Materials and Methods

In this study, we utilized an EKC model to study the haze in China using a database of 31 cities and provinces across China from the period 2000 to 2014. To verify whether there was a significant difference in the regional economic capacity or innovation ability during the same period, we first applied global regression models to grasp the overall EKC of the entire country. We adopted the EKC model, as modified by Kang et al., to study the relationships among CO2 emissions in China [31]. Interestingly, this study adopted the API, which is an overall measure of pollution that encompasses particles of all sizes, to determine the level of haze for each region, unlike previous studies, which used data regarding each pollutant. The API has been recommended as a sustainable indicator in air-quality assessment since it can overcome difficulties in interpretation of research results and problems of subjective classification of air quality [32].
In addition to the EKC analysis of haze in the form of API readings, we included the possible determinants of the environmental elements of haze—per capita GRP, trade openness, and population density––in the analysis. To verify the inverted U shape, per capita GRP and its square were included in the model. Trade openness represents the degree of economic openness, which reflects economic competitiveness and the flow of goods and services. Previous research insisted that the opening up of the economy by the Chinese government affected the degree of pollution [33]. This study also considers population density as the indirect measure of active economic activities [25]. Data regarding the GRP, trade openness, and population density were taken from the China Compendium of Statistics for the period 2000–2014 and the China Statistical Yearbook [34]. The API levels for all 31 cities in China were obtained from China’s Environmental Agency.
The EKC model used in this study is given as follows:
ln ( Y i t ) = α i t + β 1 ln ( PGRP i t ) + β 2 ln ( PGRP i t ) 2 + β 3 ln ( TR i t ) + β 4 ln ( PD i t ) + ε i t
where Y i t is the haze level measured in API per capita, PGRP stands for per capita GRP, TR is the trade openness expressed as ratio of the gross export and import value to the GRP, and PD represents the population density in each of the 31 cities in China. i and t represent the region- and time-specific data.
The following models can be utilized for analysis to consider autocorrelation within panels and cross-sectional correlation and heteroskedasticity across panels: the pooled ordinary least squares (OLS) model, the fixed- and random-effects models, and the generalized least squares (GLS) model. This is because there can be subtle or no changes in a geographical area over time, which can cause homogeneity in the parameters of each spatial unit. An F-test, a Breush–Pagan Lagrange multiplier (LM) test, and a Hausman diagnostic test were conducted to determine which model provided the best fit for the data. Pesaran’s test was also conducted to determine cross-sectional dependence.
Next, considering the spatial heterogeneity of each city and province, we performed a local linear regression using the GWR model. The GWR model has been used in earlier studies, such as in the Hedonic Pricing Model. The GWR model is also known as the local model and shows regression coefficients that vary across space, showing the spatial variations and relationships between the environmental determinants and their environmental factors via a local estimation done. The results allow us to understand how the regional determinants at a provincial level affect EKCs [2,21].
According to Wheeler [35], the GWR model is based on the linear regression model and assumes that the regression coefficient is a function of the observation point location. The GWR model focuses on the local geography, making it locally weighted. Thus, its performance from the perspective of spatial econometrics can be illustrated by using spatial heterogeneity. For the GWR model, we followed Fotheringham et al.’s approach [36] to consider regions expressed differently from the global models. This model expands the global model as follows:
Y i = β 0 ( u i , v i ) +   β k ( u i , v i ) X i k + ε i
where Y i and X i k are the dependent and independent variables, respectively. ( u i , v i ) represents the location i of the observation and ε i is the error term. The estimation for the GWR model is given as follows:
β = [ X T W ( u i , v i ) X ] 1 X T W ( u i , v i ) Y .
where W ( u i , v i ) is the square matrix of the weight assigned to a location ( u i , v i ) . X and Y are the geographically weighted matrices of the values of the independent and dependent variables. The matrix of the geographical weights, W ( u i , v i ) , is as follows:
W ( u i , v i ) = [ W 1 ( u 1 , v 1 ) 0 0 0 0 0 0 W n ( u 1 , v 1 ) ]
The GWR model estimates the coefficients for each region, assuming that a closer observation has a greater effect on the parameters than a more distant observation. Based on the equations above, the specific model for the EKC can be expressed as follows:
Y i t = α i t + β 1 ( u i , v i ) ln ( PGRP i t ) + β 2 ( u i , v i ) ln ( PGRP i t ) 2 + β 3 ( u i , v i ) ln ( TR i t ) + β 4 ( u i , v i ) ln ( PD i t ) + ε i t
where Yit is the per capita volume of the API in province or city i for the year t. ( u i , v i ) denotes the location i’s longitude and latitude coordinates. PGRP and PGRPSQ denote the provincial gross regional product per capita and its squared term. TR is the trade openness expressed as the ratio of gross export and import value to the GRP. PD is population density (in people/m2). This study utilized STATA for the global models, GWR 4 software for the spatial models, and QGIS for the mapping process. Table 1 presents the descriptive statistics.
The determinant with the largest standard deviation is the PGRP, which indicates how varied the per capita income is for different regions in China. This is also seen in the trade openness variable, where some regions have almost no trade whatsoever and other regions trade almost all of their products. There is also a surprisingly significant difference in the population density across the regions, which can be attributed to the fact that China has a large land area and some places have a more inhospitable topography, which makes it difficult to farm and earn a living. Most of China’s population is concentrated along coastal cities, with lower population density further inland.

3. Results

3.1. Subsection Results of Global Models

In this study, first, we tested the global models using the EKC formulas to understand the determinants of haze. Table 2 shows the results of the global models and the verifications of their fitness. The results demonstrate that the F-test value rejects the null hypothesis with a value of 210.19 with high statistical significance, and the Breush–Pagan LM test value also rejects the null hypothesis with a value of 2727.87 with a p-value of 0.000. Thus, both the fixed- and the random-effects model showed a better fit than the simple OLS model. However, the results of the Pesaran’s test for the cross-sectional dependence show cross-sectional dependence in the models. In this regard, this study applied the GLS model to consider the autocorrelation within panels, cross-sectional correlation, and heteroskedasticity across panels.
The results from the GLS model confirm the existence of an EKC in overall cities and provinces across China for the period 2000–2014. The PGRP variable shows a positive coefficient of 0.122, while the PGRP square’s coefficient is negative, and both variables are statistically significant. These results indicate that significant EKCs have inverted U-shapes. Moreover, the positive value of the PGDP coefficient indicates that production activities contribute to haze concentrations; this finding supports those of previous studies, which suggested that production activities are the main contributors to haze concentrations [25,37]. This result shows the necessity of balanced environmental regulation, even if China’s internal policy still prioritizes economic growth.
The trade openness has a positive coefficient at the 1% confidence level. The empirical evidence reveals that increasing trade openness results in an increase in API across China. This result is in line with previous research, which found a negative impact of the opening-up policy in China on haze pollution [33]. The population density variable has a negative coefficient at the 1% confidence level. The population density result shows that increases in population density cause a drop in API. This was an unexpected result, since population density is one of the indirect measures of active economic activity, which contributes to the increase in industrial emissions, including the haze phenomenon [25].
Even though the GLM model found the existence of the EKC for haze pollution, the global model only confirmed the existence of an EKC between overall air quality and economic performance in China. Thus, it is difficult to grasp the regional differences in the EKC [2,21,25]. Considering the possibility of significant spatial heterogeneity in the economic status of regions in China, it is necessary to investigate the heterogeneity of the EKC for each region.

3.2. Results of Geographical Weighted Regression Models

To verify the necessity of the spatial approach, we tested the spatial variability. Table 3 presents the results of the test for the spatial variability of the GWR coefficients. The differential criterion values are negative and highly significant, meaning that the variance of the regression values across the different provinces in China are also extremely high. These results demonstrate that there is significant spatial variability in terms of the model’s selection criteria.
The spatial variability test results provide strong evidence that the EKCs are not always constant but, instead, vary among provinces or cities in China. The main contributors to the PM2.5 in haze are coal and fossil-fuel burning, especially from vehicles and energy-intensive industries, and they are thus likely to be linked to the different levels of economic activity in a certain city [38]. However, the previous global OLS models could not capture these spatial differences, since the global model estimates the average relationships among all the provinces or cities in China to be in line with the findings of previous studies [21,25].
Table 4 presents the estimation results of the GWR model. The GWR parameter estimates for the independent variables show the distribution of the coefficients across regions that vary widely over space. The Akaike information criterion (AIC) value of the GWR model is also lower, at −1461.025, than the OLS value, of −324.234. Thus, the GWR model is superior to the global model.
The GWR model results confirm the spatial variability and indicate that there is significant spatial heterogeneity in the EKCs. In this regard, we investigated the existence of the EKC between the economic performance and the haze using the local coefficients from the GWR model for each of the 31 regions. One of the advantages of the GWR model is that the estimated results are spatially displayed, based on the resolution of the data used in the study. Since, in this this study, we conducted province- or city-level local regression, it is possible to display all the results using a provincial or municipal map of China generated by QGIS. Figure 1 presents the illustrated spatial variations in the EKC for haze in China.
Figure 1 shows the significant spatial heterogeneity of the EKC for haze in China. We colored the regions based on the significance and shape of their EKC (i.e., the symbols β1 and β2 in Equation (5)), which was calculated by the GWR estimation results. We verified the EKCs based on the sign and significance of the coefficients of the PGRP and the squared term of the PGRP variables. There were four cases based on the statistical significance and sign of the coefficients. Specifically, if the coefficient of β1 had a plus sign and that of β2 had a minus sign, the EKC had an inverted U-shape. Moreover, if β1 and β2 were statistically significant, it was possible to confirm the existence of EKCs in certain regions.
Six regions—Beijing, Gansu, Heilongjiang, Jiangxi, Jilin, Liaoning, Shanghai, Tianjin, Xinjiang, and Zhejiang—show inverted U-shaped EKCs. However, only three big cities (dark green colored regions in Figure 1)—Beijing, Tianjin, and Shanghai—show a statistically significant inverted U-shape between air quality and economic performance. According to the results, there are no EKCs in the other 25 provinces and cities, and only three regions—Anhui, Jiangsu, and Ningxia—are statistically significant. These findings support the necessity of balanced environmental regulation, even if China’s internal policy still prioritizes economic growth. When considering these spatial analysis results, it is important to note that the initial control policies are currently working in certain big cities, such as, Beijing, Shanghai, and Tianjin, albeit to a lesser degree of effectiveness. However, it is necessary to perform a long-term analysis of the EKCs for haze, especially in statistically non-significant regions.

3.3. Estimation of the Turning Points of the EKCs for Haze

Based on the empirical results from the GWR analyses, it was possible to calculate the turning points of the EKCs for the haze in China. We used local coefficients of the three regions that showed a significant EKC for haze. Table 5 presents the ratio of the actual GRP per population in 2017 to the calculated turning point based on the analytical results regarding the EKC for each province or city.
From these estimates, it is possible to simulate the trend in sustainable development after the analysis periods. Specifically, a ratio over 1.00 indicates that the GRP per population reached a calculated turning point. The pollutants start to decrease after this point, along with economic growth; therefore, it is possible to understand the condition of pollutants on the EKC of each province or city. Since the GWR results indicate significant spatial heterogeneity, only three big cities—Beijing, Tianjin, and Shanghai—passed their turning points. Their ratios were 1.334, 1.306, and 1.114, respectively. In these developed cities, significant political efforts and investment were made to improve the environmental quality. Consequently, these cities are more likely to achieve sustainable development than other regions.
In light of these findings, it is still difficult to define the relationship between economic growth and haze for other provinces or cities. It is possible that there are different shapes, not only for EKCs, but also for other relationships, such as linear relationships. Thus, a long-term analysis is needed for these provinces or cities, in order to capture the graph shapes of the relationship between economic growth and environmental pollutants.

4. Conclusions and Policy Recommendation

This study examined the EKC in China, focusing on both average relationships and spatial heterogeneity. The study applied the global models and the GWR model using the data regarding the haze in about 31 provinces or cities in China from 2000 to 2014. Using the GWR model, we explored the non-stationary spatial characteristics, since this model can be fitted to each data point and weighs all the observations as functions of distance from the regression point. As a result, it was possible to obtain local coefficients that varied spatially [20,21,22,23,24,25]. Based on the results of the GWR model, this study verified the existence of the EKC among regions and calculated the turning point by applying the real 2017 GRP per population data.
This study provided convincing evidence that the GWR model is valid for explaining the spatial heterogeneity of the EKC in China more effectively than the global models. The spatial variability test also confirmed that all the variables used in this study have spatial variability. The GWR model could determine the existence of an EKC for the haze and its spatial variations in China. Based on the empirical results from the GWR analysis, this study visualized the local coefficients and estimated the turning points of the EKCs. The results provided empirical evidence that only the three most developed regions in China—Beijing, Tianjin, and Shanghai—have statistically significant EKCs. Furthermore, the estimated turning points of the EKCs show that Beijing, Tianjin, and Shanghai have already passed their turning points. However, it was difficult to find an EKC for the remaining 28 provinces and cities. Moreover, the results indicated that there was no EKC for the haze in the three regions—Anhui, Jiangsu, and Ningxia.
The empirical results from the research period show that although three developed cities—Beijing, Tianjin, and Shanghai—showed a positive relationship between economic and environmental performance, while other regions were still suffering from the vicious cycle of air pollution. To achieve harmony between t environmental and economic development, it is necessary to find an effective way to achieve balanced economic growth by considering the spatial heterogeneity of the levels of development and pollution. The Chinese authorities are trying to find a way to balance productivity with pressing environmental issues through their policies and subsequent economic plans, in order to ensure the sustainability of China’s economy. For example, there are environmental regulations regarding air pollution in China, such as the Two-Control Zone policy, which was established in 1998 to define acid rain control zones in the southern regions and SO2 control zones in the north-eastern regions in China. Moreover, the Chinese government announced regulations on environmental protection, such as the Taxation Law on Environmental Protection and the Law on the Prevention and Treatment of Water Pollution to consider environmental degradation [39].
The empirical results of the current study show that the policies and laws regarding the prevention and control of atmospheric pollution need to be improved. Stronger macroeconomic regulations and control measures to promote the transformation of enterprise development can be effective. This study found significant spatial heterogeneity, confirming the EKC relationship between air quality and economic growth only in three big cities. In other words, most regions in China are still far from achieving sustainable development. From these findings, region-specific policies should be established to consider regional differences. For example, it would be effective for policymakers to tighten environmental regulations on air pollution, especially in statistically significant no EKC regions, such as Anhui, Jiangsu, and Ningxia. Additionally, monitoring and evaluation systems, including the management of haze indexes, are needed to ensure long-term effectiveness. Enhancing financial support and updating haze indexes in a timely manner can be also effective methods to improve air quality.
In addition, the Chinese government is trying to shift the development focus to balanced economic growth by converting the country’s energy structure [40]. As part of these efforts, renewable energy laws and policies were established by the National Energy Administration to improve China’s energy consumption structure. The Renewable Energy Law is an example of these efforts [41]. Since one of the main sources of haze is coal consumption, reforming the old pattern of long-term economic transformation from a “high-carbon” to a “low-carbon” economy will be effective [28]. To change the original energy consumption structure, active step-by-step actions must be taken to explore new low-pollution energy sources to reduce the consumption of fossil fuels [42].
Although this study provides empirical evidence for the spatial heterogeneity of the EKC for the haze in China, there is still room for improvement. This study applied the haze levels measured in API per capita in one representative station of each region. Estimation using API microdata in specific monitoring stations can be more effective for grasping the spatial heterogeneity of the EKC in China. Moreover, the factors that affect haze were not comprehensively evaluated in this study. For example, energy structure, industrial structure, capital intensity, and environmental management capacities can be considered. These factors were not evaluated mostly due to the existence of multicollinearity and data availability. Thus, it may be necessary to investigate more socio-economic factors and specify them both theoretically and statistically to design effective policies to achieve balanced economic growth in China. Furthermore, further studies could investigate the relationship between environmental quality and economic performance with more advanced models and data, such as through the use of the Distance between Indices of Simulation and Observation, which is a new comprehensive index presenting the overall performances of different models [43,44]. A comprehensive approach to determining regional trends and the current state of pollution will help to inform future policy on sustainable development in China.

Author Contributions

Conceptualization, A.S.A.-R., Y.K. and L.Y.; methodology, A.S.A.-R., Y.K. and L.Y.; software, Y.K. and L.Y.; validation, Y.K.; formal analysis, A.S.A.-R., Y.K. and L.Y.; investigation, A.S.A.-R., Y.K. and L.Y.; resources, A.S.A.-R. and L.Y.; data curation, A.S.A.-R., Y.K. and L.Y.; writing—original draft preparation, A.S.A.-R., Y.K. and L.Y.; writing—review and editing, Y.K.; visualization, Y.K.; supervision, A.S.A.-R.; project administration, A.S.A.-R. and Y.K. All authors have read and agreed to the published version of the manuscript.

Funding

No funding was received for conducting this study.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data in this study are available upon reasonable request.

Conflicts of Interest

The authors declare that they have no conflict of interest.

References

  1. Ma, L.J. Urban transformation in China, 1949–2000: A review and research agenda. Environ. Plan A 2002, 34, 1545–1569. [Google Scholar] [CrossRef]
  2. Kim, Y. Spatial analysis of factors affecting location decisions of Korean companies in China. Reg. Sci. Policy Pract. 2019, 11, 383–401. [Google Scholar] [CrossRef]
  3. World Health Organization (WHO). Ambient Air Pollution: A Global Assessment of Exposure and Burden of Disease; World Health Organization: Geneva, Switzerland, 2016; Available online: https://apps.who.int/iris/handle/10665/250141 (accessed on 15 April 2022).
  4. Zhuang, X.; Wang, Y.; He, H.; Liu, J.; Wang, X.; Zhu, T.; Ge, M.; Zhou, J.; Tang, G.; Ma, J. Haze insights and mitigation in China: An overview. J. Environ. Sci. 2014, 26, 2–12. [Google Scholar] [CrossRef]
  5. Quan, J.; Zhang, Q.; He, H.; Liu, J.; Huang, M.; Jin, H. Analysis of the formation of fog and haze in North China Plain (NCP). Atmos. Chem. Phys. 2011, 11, 8205–8214. [Google Scholar] [CrossRef] [Green Version]
  6. Gao, M.; Guttikunda, S.K.; Carmichael, G.R.; Wang, Y.; Liu, Z.; Stanier, C.O.; Saide, P.E.; Yu, M. Health impacts and economic losses assessment of the 2013 severe haze event in Beijing area. Sci. Total Environ. 2015, 511, 553–561. [Google Scholar] [CrossRef]
  7. Jalil, A.; Mahmud, S.F. Environment Kuznets curve for CO2 emissions: A cointegration analysis for China. Energy Policy 2009, 37, 5167–5172. [Google Scholar] [CrossRef] [Green Version]
  8. Hao, Y.; Niu, X.; Wang, J. Impacts of haze pollution on China’s tourism industry: A system of economic loss analysis. J. Environ. Manag. 2020, 295, 113051. [Google Scholar] [CrossRef]
  9. Long, X.; Naminse, E.Y.; Du, J.; Zhuang, J. Nonrenewable energy, renewable energy, carbon dioxide emissions and economic growth in China from 1952 to 2012. Renew. Sustain. Energy Rev. 2015, 52, 680–688. [Google Scholar] [CrossRef]
  10. Wang, Y.; Han, R.; Kubota, J. Is there an environmental Kuznets curve for SO2 emissions? A semi-parametric panel data analysis for China. Renew. Sustain. Energy Rev. 2016, 54, 1182–1188. [Google Scholar] [CrossRef]
  11. Shen, J. A simultaneous estimation of environmental Kuznets curve: Evidence from China. China Econ. Rev. 2006, 17, 383–394. [Google Scholar] [CrossRef]
  12. Tao, S.; Zheng, T.; Lianjun, T. An empirical test of the environmental Kuznets curve in China: A panel cointegration approach. China Econ. Rev. 2008, 19, 381–392. [Google Scholar]
  13. Shu, L.; Fantang, Z.; Huaiyang, F.; Zhencheng, X. An empirical test of the environmental Kuznets curve in Guangdong Province, China. APCBEE Procedia 2012, 1, 204–209. [Google Scholar] [CrossRef] [Green Version]
  14. Wang, Q.; Yuan, X.; Lai, Y.; Ma, C.; Ren, W. Research on interactive coupling mechanism and regularity between urbanization and atmospheric environment: A case study in Shandong Province, China. Stoch. Environ. Res. Risk Assess. 2012, 26, 887–898. [Google Scholar] [CrossRef]
  15. Song, M.-L.; Zhang, W.; Wang, S.-H. Inflection point of environmental Kuznets curve in Mainland China. Energy Policy 2013, 57, 14–20. [Google Scholar] [CrossRef]
  16. Du, G.; Liu, S.; Lei, N.; Huang, Y. A test of environmental Kuznets curve for haze pollution in China: Evidence from the penal data of 27 capital cities. J. Clean. Prod. 2018, 205, 821–827. [Google Scholar] [CrossRef]
  17. Xu, T. Investigating environmental Kuznets curve in China—Aggregation bias and policy implications. Energy Policy 2018, 114, 315–322. [Google Scholar] [CrossRef]
  18. Xie, Q.; Xu, X.; Liu, X. Is there an EKC between economic growth and smog pollution in China? New evidence from semiparametric spatial autoregressive models. J. Clean. Prod. 2019, 220, 873–883. [Google Scholar] [CrossRef]
  19. Pata, U.K.; Caglar, A.E. Investigating the EKC hypothesis with renewable energy consumption, human capital, globalization and trade openness for China: Evidence from augmented ARDL approach with a structural break. Energy 2021, 216, 119220. [Google Scholar] [CrossRef]
  20. Hao, Y.; Liu, Y.-M. The influential factors of urban PM2.5 concentrations in China: A spatial econometric analysis. J. Clean. Prod. 2016, 112, 1443–1453. [Google Scholar] [CrossRef]
  21. Kim, Y.; Tanaka, K.; Ge, C. Estimating the provincial environmental Kuznets curve in China: A geographically weighted regression approach. Stoch. Environ. Res. Risk Assess. 2018, 32, 2147–2163. [Google Scholar] [CrossRef]
  22. Ma, Y.-R.; Ji, Q.; Fan, Y. Spatial linkage analysis of the impact of regional economic activities on PM2.5 pollution in China. J. Clean. Prod. 2016, 139, 1157–1167. [Google Scholar] [CrossRef]
  23. Zeraibi, A.; Ahmed, Z.; Shehzad, K.; Murshed, M.; Nathaniel, S.P.; Mahmood, H. Revisiting the EKC hypothesis by assessing the complementarities between fiscal, monetary, and environmental development policies in China. Environ. Sci. Pollut. Res. 2022, 29, 23545–23560. [Google Scholar] [CrossRef] [PubMed]
  24. Chuai, X.; Huang, X.; Wang, W.; Wen, J.; Chen, Q.; Peng, J. Spatial econometric analysis of carbon emissions from energy consumption in China. J. Geogr. Sci. 2012, 22, 630–642. [Google Scholar] [CrossRef]
  25. Kim, Y.; Tanaka, K.; Zhang, X. A spatial analysis of the causal factors influencing China’s air pollution. Asian J. Atmos. Environ. 2017, 11, 194–201. [Google Scholar] [CrossRef]
  26. Brajer, V.; Mead, R.W.; Xiao, F. Searching for an environmental Kuznets curve in China’s air pollution. China Econ. Rev. 2011, 22, 383–397. [Google Scholar] [CrossRef]
  27. Fang, M.; Chan, C.K.; Yao, X. Managing air quality in a rapidly developing nation: China. Atmos. Environ. 2009, 43, 79–86. [Google Scholar] [CrossRef]
  28. Beşe, E.; Friday, H.S.; Özden, C. Coal consumption environmental Kuznets curve (EKC) in China and Australia: Evidence from ARDL model. J. Appl. Bus. Econ. 2020, 22, 25–36. [Google Scholar]
  29. Koshta, N.; Bashir, H.A.; Samad, T.A. Foreign trade, financial development, agriculture, energy consumption and CO2 emission: Testing EKC among emerging economies. Indian Growth Dev. Rev. 2020, 14, 50–80. [Google Scholar] [CrossRef]
  30. Stern, D.I. The rise and fall of the environmental Kuznets curve. World Dev. 2004, 32, 1419–1439. [Google Scholar] [CrossRef]
  31. Kang, Y.Q.; Zhao, T.; Yang, Y.Y. Environmental Kuznets curve for CO2 emissions in China: A spatial panel data approach. Ecol. Indic. 2016, 63, 231–239. [Google Scholar] [CrossRef]
  32. Payus, C.M.; Syazni, M.N.; Sentian, J. Extended air pollution index (API) as tool of sustainable indicator in the air quality assessment: El-Nino events with climate change driven. Heliyon 2022, 8, e09157. [Google Scholar] [CrossRef]
  33. Dong, F.; Zhang, S.; Long, R.; Zhang, X.; Sun, Z. Determinants of haze pollution: An analysis from the perspective of spatiotemporal heterogeneity. J. Clean. Prod. 2019, 222, 768–783. [Google Scholar] [CrossRef]
  34. National Bureau of Statistics of China. Available online: http://www.stats.gov.cn (accessed on 15 April 2022).
  35. Wheeler, D.C. Geographically weighted regression. In Handbook of Regional Science; Fischer, M.M., Nijkamp, P., Eds.; Springer: Berlin/Heidelberg, Switzerland, 2014; pp. 1435–1459. [Google Scholar]
  36. Fotheringham, A.S.; Brunsdon, C.; Charlton, M. Geographically Weighted Regression: The Analysis of Spatially Varying Relationships; John Wiley & Sons: Chichester, UK, 2002. [Google Scholar]
  37. Shi, T.; Zhang, W.; Zhou, Q.; Wang, K. Industrial structure, urban governance and haze pollution: Spatiotemporal evidence from China. Sci. Total Environ. 2020, 742, 139228. [Google Scholar] [CrossRef] [PubMed]
  38. Huang, R.-J.; Zhang, Y.; Bozzetti, C.; Ho, K.-F.; Cao, J.-J.; Han, Y.; Daellenbach, K.R.; Slowik, J.G.; Platt, S.M.; Canonaco, F.; et al. High secondary aerosol contribution to particulate pollution during haze events in China. Nature 2014, 514, 218–222. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Sun, W.; Yang, Q.; Ni, Q.; Kim, Y. The impact of environmental regulation on employment: An empirical study of China’s Two Control Zone policy. Environ. Sci. Pollut. Res. 2019, 26, 26472–26487. [Google Scholar] [CrossRef]
  40. Zhang, K.; Dearing, J.A.; Tong, S.L.; Hughes, T.P. China’s degraded environment enters a new normal. Trends Ecol. Evol. 2016, 31, 175–177. [Google Scholar] [CrossRef]
  41. Liu, J. China’s renewable energy law and policy: A critical review. Renew. Sustain. Energy Rev. 2019, 99, 212–219. [Google Scholar] [CrossRef]
  42. Bekun, F.V. Mitigating emissions in India: Accounting for the role of real income, renewable energy consumption and investment in energy. Int. J. Energy Econ. Policy. 2022, 12, 188–192. [Google Scholar] [CrossRef]
  43. Zhou, Q.; Chen, D.; Hu, Z.; Chen, X. Decompositions of Taylor diagram and DISO performance criteria. Int. J. Climatol. 2021, 41, 5726–5732. [Google Scholar] [CrossRef]
  44. Hu, Z.; Chen, X.; Zhou, Q.; Chen, D.; Li, J. DISO: A rethink of Taylor diagram. Int. J. Climatol. 2019, 39, 2825–2832. [Google Scholar] [CrossRef]
Figure 1. EKC for the haze in China.
Figure 1. EKC for the haze in China.
Atmosphere 13 00806 g001
Table 1. Descriptive Statistics of Variables.
Table 1. Descriptive Statistics of Variables.
VariableUnitObs.MeanMedianMaximumMinimumStd. Dev.
APIIndex465294.518307.000365.000168.00048.783
PGRP100 million CNY46525,076.28019,580.230103,671.3002742.07019,679.650
TR%46513.9293.600100.2200.00122.425
PDPerson/m2465414.479263.0004349.0002.000621.088
Note: The data observations for each region for 15 years are consistent with the calculation of 15 years multiplied by 31 regions, which equals 465.
Table 2. Estimation Results of Global Models.
Table 2. Estimation Results of Global Models.
VariablePool OLSFEREGLS
Intercept5.518 ***
(1.167)
6.356 ***
(0.539)
5.938 ***
(0.376)
5.675 ***
(0.165)
ln(PGRPit)0.156
(0.236)
−0.008
(0.078)
0.032
(0.069)
0.122 ***
(0.034)
ln(PGRPit)2−0.013
(0.007094)
−0.003
(0.004)
−0.005
(0.004)
−0.011 ***
(0.002)
ln(TRit)0.014 *
(0.007)
0.001
(0.004)
0.001
(0.004)
0.013 ***
(0.001)
ln(PDit)−0.012
(0.009)
−0.059
(0.046)
−0.013
(0.019)
−0.011 ***
(0.001)
R20.1680.0900.150
F-test/Breush–Pagan LM test 210.19 ***2727.87 ***
Pesaran’s test 4.041 ***4.948 ***
Note: ***, and * indicate significance at the 1% and 10% levels, respectively.
Table 3. Spatial Variability Test.
Table 3. Spatial Variability Test.
VariableFDiff of Criterion
Intercept68,192.880 ***−3960.470
ln(PGRPit)250.132 ***−1049.020
ln(PGRPit)21473.370 ***−1915.418
ln(TRit)54.240 ***−379.452
ln(PDit)1516.960 ***−1548.354
Note: Positive value of diff-criterion indicates that there is no spatial variability in terms of model selection criteria. *** indicates significance at the 1% level.
Table 4. Estimation Results of GWR Model.
Table 4. Estimation Results of GWR Model.
VariableMinLQMedUQMax
Intercept0.1964.2407.3489.66716.167
ln(PGRPit)−1.585−0.726−0.2060.4431.120
ln(PGRPit)2−0.053−0.0260.0080.0320.079
ln(TRit)0.031−1.099−0.0740.0130.622
ln(PDit)−1.099−0.3250.0160.0470.147
N465
Adjusted R20.939
AIC−1461.025
Table 5. Calculated Turning Points of EKC; unit: 100 million CNY.
Table 5. Calculated Turning Points of EKC; unit: 100 million CNY.
Region2017 PGRPLn (2017 PGRP)Turning PointRatio
Beijing128,99411.7682.9441.334
Tianjin118,94411.6862.7351.306
Shanghai126,63411.7491.2051.114
Note: The turning point is the calculated inflection point of EKCs from the GWR results. The ratio is calculated by equation: 1 + (PCGRP2017 − Turning point)/Turning point. The ratio 1.000 means that the PGRP in 2017 reached the turning point. The turning point is calculated only for the regions with significant EKCs—Beijing, Tianjin, and Shanghai—that show a statistically significant inverted U-shape.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Abdul-Rahim, A.S.; Kim, Y.; Yue, L. Investigating Spatial Heterogeneity of the Environmental Kuznets Curve for Haze Pollution in China. Atmosphere 2022, 13, 806. https://0-doi-org.brum.beds.ac.uk/10.3390/atmos13050806

AMA Style

Abdul-Rahim AS, Kim Y, Yue L. Investigating Spatial Heterogeneity of the Environmental Kuznets Curve for Haze Pollution in China. Atmosphere. 2022; 13(5):806. https://0-doi-org.brum.beds.ac.uk/10.3390/atmos13050806

Chicago/Turabian Style

Abdul-Rahim, Abdul Samad, Yoomi Kim, and Long Yue. 2022. "Investigating Spatial Heterogeneity of the Environmental Kuznets Curve for Haze Pollution in China" Atmosphere 13, no. 5: 806. https://0-doi-org.brum.beds.ac.uk/10.3390/atmos13050806

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