{'fraction_non_alphanumeric': 0.02967786549876102, 'fraction_numerical': 0.01630841929349392, 'mean_word_length': 5.087879319417646, 'pattern_counts': {'":': 0, '<': 2, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Spatial Econometric Analysis on the Exports-PM2.5 Emissions Nexus: Evidence from the Chinese Smog Provinces

The study presents an empirical analysis to investigate the effect of export volume and significant factors on PM2.5 emissions in the highest concentrations of air pollution provinces in China during 1998-2015. This paper proposes a spatial econometric model based on an extended Stochastic Impact by Regression on Population, Affluence, and Technology (STIRPAT) model. The results revealed evidence that spatial spillover effects of PM2.5 emissions exist among major provinces. Provinces in the northeast, south-central China, and north China regions have more potential in reducing PM2.5 emissions. Exports per capita have negative and significant effects on PM2.5 emissions, which revealed that an increase in own-province exports per capita would reduce the PM2.5 concentrations of both own province and the neighboring provinces. Also, the results confirm the existence of an inverted-U shaped Environmental Kuznets Curve (EKC) between PM2.5 emissions and economic growth. Surprisingly, the study doesn’t validate the pollution haven hypothesis (PHH). Along a similar line, the empirical evidence has shown that the private vehicle has a significant impact on PM2.5 concentrations, whereas energy intensity exerts a negative effect. Other policy suggestions put forth for policymakers on reducing PM2.5 emissions.


Introduction
Extensive evidence has suggested that China is one of the most rapid economic growth countries in the world over the last decade. Meanwhile, air pollution and industrialization with smokestacks problems have become increasingly prominent [1] [2], and recently become increasingly with a rapid increase in industrialization and urbanization [3].
The word "smog" a mixture of smoke and fog trapped beneath a temperature inversion [4], gases, and particulate matter, which produces haze and decreases visibility. Smog is mostly caused by emissions of industrial and vehicular, which responded within high concentrations of fine particles in the atmosphere [5]. Fine particulate matter (PM2.5 emissions) originates primarily from many sources such as severe smog and haze episodes, vehicle exhaust, and road dust, these pollutants decrease the atmospheric visibility by scattering and absorbing ambient light and cause serious harm to ecosystem health, and massive economic losses in many developing countries [6], and threaten the human health [7].
According to the World Health Organization (WHO) [8], "China PM2.5 median equivalent 54, currently, it becomes one of the largest emitters of global emissions, mainly evident by the infamous smog that has happened frequently in many provinces particularly in eastern and northern regions". It's evident that various air pollution problems have emerged in China associated with the expansion of trade such as CO 2 emissions, carbon monoxide (CO), and sulfur dioxide (SO 2 ). And PM2.5 emissions have become increasingly serious in recent times [9] [10]. After China joined the World Trade Organization (WTO) in 2001, it is exports continued to grow steadily [11] [12]. Moreover, a series of environmental problems have arisen to an unprecedented level affecting seriously the quality of life.
To control air pollution Chinese government has adopted a series of policies deemed to solve the problem, among those policies such as subsidies for use of new energy-efficient cars to decrease automobile exhaust emission, provided PM2.5 real-time air quality monitoring data for the public, in addition to the Act for controlling atmospheric emissions from coal combustion. The control policy measures and prevention of PM2.5 emissions have always been an important task in China due to the repetitive occurrence of pollution. Therefore, to provide a significant research basis for effectively preventing PM2.5 concentrations and reducing air pollution, it is of great theoretical and realistic significance to study the impact of exports on PM2.5 emissions.
Globally, several studies are attempting to investigate the effects of trade on air pollution which is useful in designing an effective climate change program and reducing global emissions in international trade [13]. Research-related interest has recently demonstrated that exports have been an inseparable correlate with air pollution, and the pollutants emitted in China as a result of the production of goods for export to foreign consumption [14] [15].
Many existing recent-studies provided evidence for the impact of exports on air pollution, but empirical results are mixed. For example, Xu et al. [15] revealed that exports promoted China's air pollutant emissions. Zhao et al. [16] argued that about 55% -62% of production emissions due to exports and most provinces' exports go to the developed east coast and bordering provinces in China. In addition, Lin [17] found that trade openness had a significant positive effect on increasing air pollution. Recently, Ning et al., [18] verified that the interregional trade played a significant role in spillover and feedback effects relating to emissions, and examined that the emission spillover effect was strong in the border provinces.
In contrast, Du et al. [19]  proven that the rapid growth in exports and industrialization among the provinces was the most significant factor that can explain the air pollution status [23] [24] [25]. Regardless of the type of emissions, Prell et al. [26] examined the impact of structural changes in international trade on the environment and found that a country with a higher level of integration has much more pollution. We conclude that the existing studies on the role of international trade on air pollution are inconclusive and have not reached a consensus and still exist some deficiencies.
Along a similar line, many studies have investigated the role of socioeconomic indicators on environmental pollutants by expanded in the Stochastic Impacts by Regression on Population, Affluence, and Technology (STIRPAT) identity by adding some independent variables in the basic model [3] [27], and others by used the spatial econometric approach with extending the STIRPAT model [28] [29]. But Authors do not discuss their results in terms of specification tests of the spatial econometric along the lines developed by [30]. During different historical periods, the spatial effects of exports on PM2.5 emissions in the provinces should not be ignored due to the development stages. However, the recent literature review has not given sufficient attention to these approaches to PM2.5 emissions. Cheng et al. [28] is the only group that has used the expanded STIRPAT model in a spatial approach, in which the authors identify the key factors that influence PM2.5 pollution in China. Their findings are consistent with our study. But they used data from 2001-2012 in the city-level, not provincial level, which raises the concern that their results may be partial. While some studies ignored the role of the spatial effect and correlation of PM2.5 emissions between different provinces, our study aims and attempting to fill the gaps. As a result, the study investigates and detects the spatial autocorrelation of PM2.5 emissions among the highest concentrations of air pollution provinces in China, and depending on previous research, the study provides significant contributions as follows: First, employ a novel spatial econometric technique, to explore the role of export volume and the important factors that influence PM2.5 concentrations from economic aspects. To the best of our knowledge, the kind of this research has not given sufficient attention in these provinces. Second, more considering for developed and incorporate between the expanded STIRPAT model and spatial econometrics model in PM2.5 emission in this distinguished study; will avoid the bias estimation caused due to use the traditional panel econometric models, such as fixed effects, random effects, GMM method, and co-integration method. Third, although the cities-level data provide observations within the narrowly geographical area in the specific regions. Our study used provincial-level long-term panel data (1998)(1999)(2000)(2001)(2002)(2003)(2004)(2005)(2006)(2007)(2008)(2009)(2010)(2011)(2012)(2013)(2014)(2015) to avoid the erroneous estimates of the spatial spillover effect and the neglect of the complex links between PM2.5 emissions and economic activities resulting from the lack of data for PM2.5 emission of the city-level before 2013. Furthermore, by using the spatial model, we will explore whether the PM2.5 emissions Kuznets curve (PMKC) exists among the provinces or not.
The remainder of the study is structured as follows: In Section 2 we introduce the methodology including model specification and data sources. While Section 3 provides empirical results and discussion and Section 4 concludes.

Model Specification
As a global issue, PM2.5 emissions have been studied widely from different perspectives. The STIRPAT model is one of the most recent measures that has been developed and applied to explore the influence of anthropogenic factors on air pollution indicators, such as PM2.5 emissions. According to the literature review, the original formulation of the IPAT identity presented by Ehrlich's [31], who first found evidence of the existence of the relationship between population and environment, known as I = P•F(P). Thus, I represent environmental impact, P is the population size, and F is a function that measures per capita impact. Commoner [32] played an essential role in the formulation of the IPAT identity in the algebraic formulation, to address policies related to increasing environmental degradation [33]. Consequently, the IPAT identity has become well-known in 1972, by [34] as follows: where, A for affluence, and T for technology.
Recently, several reformations of the IPAT identity had conducted by adding new factors, for instance, [35] [36] [37]. Later, the IPAT model and its all refor- mulations was considered to be too simple because it is only a mathematical identity, i.e., it does not allow testing a certain hypothesis, and it does not permit the non-propositional and non-monotonic changes in human indicators. To address these limitations, Dietz and Rosa [38] proposed that the IPAT accounting identity would be recast into a stochastic form labeled STIRPAT model, to allow random errors in the estimation of parameters and provide a testable model to estimate the contributions of anthropogenic factors on greenhouse gas (GHG) and other emissions, they assume that the STIRPAT model will treat the estimated parameters of the original I PAT = identity, which suppose the elasticity coefficients equal 1, (i.e., where I is for environmental impact, and T for technology. Furthermore, a is the constant term, b, c, and d, which are the coefficients for population, affluence, and technology, respectively and e indicate the residual term. This version allows tests of hypotheses regarding other factors than population and affluence that may modify environmental impact. The model can be expressed in logarithmic form as: Thus, the coefficients in Equation (3) can be thought of as elasticities. Therefore, the elasticity coefficient in the STIRPAT model is not always equal to 1, and this is the significant difference between STIRPAT and the IPAT model.
Several studies have built on the STIRPAT model by expanding decomposition to include more terms to capture greater complexity, such as [37] [39] expanded by adding quadratic terms, to assess the relationship between IPAT, ImPACT, and STIRPAT formulations and to discuss an inverted U-shaped EKC.
Fan et al. [40]  First, the structure of energy consumption; many scholars had to focus on energy consumption that directly results in PM2.5 emission to explain the fundamental relationship, with this perspective in consideration, China accelerated industrialization with energy consumption in the secondary industry much higher than in other sectors [28]. Furthermore, most energy consumption in China comes from coal, and it is becoming the main energy source with rich re- serves and low prices [42]. Therefore, consume a lot of coal will emit large amounts of PM2.5. Thus, the structure of energy consumption has been included in the framework of our model. Because of this, we use the ratio of province coal consumption to total energy consumption to measure the structure of energy consumption.
where PM2.5 is the fine particulate matter (PM2.5 emissions) of province i, POPD is the population density (expressed as the population divided by the area of each province), GDPP is real GDP per capita represent affluence, GDP per capita is commonly used to verify the inverted-U shaped the EKC hypothesis this study seeks to test this hypothesis and its spatial pattern at the provincial level in China, ENI is energy intensity represent technology level, we included energy intensity proxied with total energy consumption divided by total outputs to represent the technology level. EX denotes exports per capita, SEC is the structure of energy consumption, PVEH is the amount of possession of private vehicles, and ε is the standard error term.

Global Spatial Autocorrelation Model
The model begins with the global spatial autocorrelation index to explain the degree of spatial correlation among the provinces. For the robustness of our model, the study performs global Moran's I index [46] as follows: where w ij represents the spatial weight matrix element to measure the distance between provinces i and j, if the provinces have a mutual border the value of the element in the matrix equals 1 otherwise it equals 0, the term Y denotes the ob-servation of PM2.5 emissions and S 2 is the standard deviation of observation.
The values of the global Moran's I index ranges from −1 to 1. The negative values indicate that air pollution exists negative spatial autocorrelation, while the positive values imply that air pollution exists positive spatial autocorrelation, the values −1 and 1 indicate perfect dispersion and perfect correlation respectively. The zero value means that there is no spatial autocorrelation or random spatial pattern.
Typically, Moran's I index scatter reflects the existence of autocorrelation between observations and the neighbors [47]. The local Moran's I scatter divided into four quadrants. The first and third quadrant reflects the positive autocorrelation, inversely; the second and fourth quadrant denotes negative autocorrelation.

Spatial Econometric Models
Based on the extended STIRPAT model in Equation (4) we will construct the spatial econometric model, by taking into consideration the fact that PM2.5 emissions are heterogeneous and spatially correlated among provinces. Therefore, the study follows the spatial econometric models proposed by [48] which allow describing the spatial correlation base on three basic models.
First, the spatial lag model (SLM), the SLM model aims to explores whether the dependent variable in a specific province is partially determined by a spatially weighted average of neighboring dependent variables [49] [50], this means that PM2.5 emission in province i affected by the surrounding conditions of PM2.5 emission intensity in province j. The SLM is formulated as follows: where Y it denotes the dependent variable for the cross-sectional province i at time t (i.e. vector of dependent variables), X is a matrix of explanatory variables, W ij is the element of a spatial weighting matrix, is the N × N weight matrix, ρ is the scalar spatial autoregressive coefficient, which explores the effects of the spatial factors on the dependent variable, β is a column vector parameter and matching with X i , the term μ indicates individual-specific effect for each province, η represent a time-period specific effect, ε is the error term assumed to be ε it ≈ N(0, σ 2 I) and the term ΣW ij Y jt is the interaction impact of the dependent variable in province i with the dependent variables in the neighboring provinces. Second is a spatial error model (SEM), in SEM the error term in province i is affected by neighboring error term in province j according to spatial weights matrix W and the random error component, the SEM is expressed as follows: where ϕ it is the spatially autocorrelated error term, the term λ is the spatial au- where γ, like β, is column vector parameters and matching with the explanatory variables.

Estimation Methods and Relevant Tests
This study primarily focuses on the specification tests outlined by [48], as follows: Firstly, to examine which type of the model is the best to fits the data, we estimate and test the non-spatial fixed models and will apply the likelihood ratio (LR) test to determine fixed effects, and then test the spatial autocorrelation of the residual error to determine what kind of fixed effect should be included using Lagrange multiplier [51] tests and their robustness [51]. Ordinary least-squares (OLS) are usually used to estimate the traditional mixed panel data models, while LM tests are usually used to produce consistent estimates for spatial autocorrelation of the residual error [52]. Besides the traditional fitting goodness (R2) and LR test, the corrected R2 was also used to test the spatial panel model. Second: to verify which spatial model is the suitable specification, we will test either the SDM model, the SLM model, or the SEM model according to the values of the Wald test and LR test.
The null hypothesis of Wald test (H0: γ = 0) to determine if the SDM model can be simplified to the SLM model or not, against the null hypothesis of LR test (H0: γ + ρβ = 0) to examines whether if the SDM model can be simplified to the SEM model [30]. Both tests follow a chi-square distribution. If we reject both hypotheses, then the SDM model represents the best fit to the data, and if we accept the first hypothesis suggests, then the SLM model best describes the data. Similarly, accept the second hypothesis indicates that the SEM model best describes the data.

Study Area
The study uses annual data of 21 provinces (from the highest concentrations of air pollution regions in China) which have relatively large PM2.5 emissions on average and unfavorable meteorological conditions for pollution dispersion [53], due to the lack of data Taiwan province, Hong Kong and Macao special administrative regions are not included. The sample selection represents about 68% of the total Chinese provinces. We depend on the official division method to divide the 21 provinces of China into four regions as follows (Table 1). Figure 1 presents the spatial characteristics of PM2.5 emissions in our sample provinces, it is clear that the dominant pattern of PM2.5 emissions increases Open Access Library Journal

Data Acquisition
The main variable of interest in the analysis is the PM2.5 emissions, which come from the National Centers for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR), we used provincial panel data for 1998-2015 for 21 provinces in China. The mean value of PM2.5 concentration is 47.50 mg/m 3 , which exceeded the level of annual PM2.5 standard (35 mg/m 3 ). The overall maximum and minimum also high during this period and highlights the severity of the air quality problem faced by China provinces since the 21st century.

Spatial Distribution of PM2.5 Emissions in China Provinces
First, we examine the possible existence of PM2.5 emissions' spatial autocorrelation. According to Equation (5) Table 2 demonstrates the global Moran's I index, the test indicates that the values of global Moran I vary from 0.556 to 0.495, and PM2.5 emissions display positive spatial autocorrelation at a 1% significant level, which point out that smog pollution in China exists in an extremely spatial autocorrelation. Furthermore, the overall Moran I over eighteen years is 0.501, which displays also a positive spatial correlation at the same significant level.
The positive values of Moran I indicate that the areas with high PM2.5 emissions (provinces in high-high groups) tend to cluster together, whereas the low PM2.5 emissions intensities, cluster together (provinces in low-low groups) ( Table 1). These provinces covered by PM2.5 concentration higher than (35 µg/m 3 ) were heavily polluted regions (see Figure 1). For provincial data, a coincidence of attribute similarity spatial autocorrelation is almost unavoidable, and if ignored like this autocorrelation, it can lead to misleading outcomes [54]. Our results similar to [55] also demonstrated that the economically developed and eastern coastal areas in China with a higher degree of urbanization will suffer more from emissions concentration than others. Therefore, emissions spillover Furthermore, Moran scatter plots were constructed to examine the spatial clustering pattern of emissions. Figure 2 shows that most provinces are located in the first (high-high) and third (low-low) quadrants which indicates the positive spatial autocorrelation.
The result indicates that the provinces in the high-high quadrants have more potential for exporting emission to their neighboring provinces, in contrast, provinces with low-low air pollution "cluster together" with other low polluted provinces. The

Empirical Results of Econometric Models and Discussion
To determine which model of those defined in Sections 2.3 fits the data best, and due to the existence of significant spatial autocorrelation as explained above, the study utilizes non-spatial panel model analysis to examine exists of spatial correlation among spatial units using classic Lagrange multiplier [51] tests before building the spatial model for the effect of exports on air pollution. If we reject the non-spatial models, the spatial model should be used to capture the spatiality through the procedures mention in Section 2.4. Table 3 shows the results of the non-spatial models.  To examine the spatial dependence variables, LM tests were employed and its robust tests, all the four kinds of fixed effects support the spatial lag model at a 1% significance level. Similarly, the error spatial lag LM test and the robust spatial lag LM tests with pooled-effects, spatial, and time-fixed effect and time-fixed effect are significant. Thus, the results reject the hypothesis that spatial correlations do not exist, which confirm exist spatial correlation among the data.
The results of the LR test reject the null hypotheses that spatial fixed effects and time fixed effects do not exist at a 1% significance level, indicating the significance of introducing spatial dependence items [52], and also justify the extension of the model by including both the space-fixed effects and time-period fixed effects.
To investigate whether the study depends on the SLM model or SEM model instead of the SDM model through performing both the Wald test and the LR test. The results of the Wald test (33.50, with 7-degree freedom, P < 0.01) to determine if the SDM model can be simplified to the SLM model is rejected at a 1% significance level. Likewise, the results of the LR test (27.10, with 7-degree freedom, P < 0.01) to explain if the SDM model can be simplified to the SEM model also rejected at the same significance level. Such a finding means the SDM model is the convenient specification for this regression.
Furthermore, to determine the best model, we conducted a Hausman test. The estimation results in Table 4, denoted to reject random effects at a 1% significance level (104.81), indicating that the more appropriate specification model is the model with fixed effects. Table 4 shows that the spatial fixed effects model has the highest goodness-of-fit statistics value (0.6146) and log-likelihood value (439.7567) suggesting that the spatial fixed effects model outperforms the other three models. Therefore, the interpretation will limit on it is coefficients.
The first column in Table 4, shows the coefficient of the spatial lagged term of the dependent variable is positive and significant, this result indicates that PM2.5 emissions in the neighboring provinces have a positive effect on local PM2.5 concentrations. This value is consistent with the outcome of Moran's I Table 4. Estimation results of the spatial Durbin model.  Note: The symbols *, **, and *** denotes a significance level of 10%, 5% and 1% respectively. index in Figure 2 and LR tests for spatial autocorrelation in Table 5. The result indicates that for every 1% increase in the average of PM2.5 emission of neighborhood provinces, the PM2.5 emission of home provinces increases by more than 0.76% net of other explanatory variables, this coefficient implies that the provinces with similar air pollution emissions generally gather together. Thus due to the existence of spatial autocorrelation, the estimated coefficients of the independent variables in the SDM model cannot be expressed as marginal effects [48], and cannot accurately reflect the spatial spillover effect of provincial PM2.5 emissions. Therefore, the study depends on the direct, indirect, and total effects to quantify the impacts of explanatory variables and their spatial spillover effect on air pollution.
The results in Table 6 summarize the decomposition estimated direct and indirect effects based on the coefficients of the SDM model in Table 4. The results of direct effects are close to the corresponding regression results of the spatial fixed effects in the SDM model, which explains that the results are steady and effective. But the difference between them in the values of the coefficients estimated due to exit feedback effects which pass through adjacent provinces, and then back to the provinces themselves, which contains two parts; the impacts through the value of a spatially lagged dependent variable (ρ∑W*lnPM2.5) and the other part comes from the coefficients of the spatially lagged independent variables (∑W ij X itγ ).
With respect to the indirect effects, the results show that only three explanatory variables such as (population density, GDP per capita, and private vehicle) have significant effects, the result indicates that the coefficients estimated of Table 5 Note: The symbols *, **, and *** denotes a significance level of 10%, 5% and 1% respectively. these variables in one specific province have a spatial spillover effect on the corresponding factors in the neighboring provinces. Therefore, we can argue that spillover effects exist.
The direct effect of the export variable appears to be −0.0526 with a significant level of 5%, due to the spatial aggregation, while the coefficient estimated by the SDM in the fixed-effects model was −0.0422. Hence, the feedback effects amount to −0.0104 or 19.8 percent of the direct effect, indicating that the elasticity of exports per capita is underestimated. Moreover, the indirect effect of a change in the export variable in the SDM model seems to be 184.98% of the direct effect.
These results demonstrate that an increase in own-province exports per capita would reduce PM2.5 emissions in their own province and neighboring provinces, hence the relationship between exports and PM2.5 emissions doesn't validate the pollution haven hypothesis (PHH), which states that high-income per capita countries export their pollution to the lowest countries.  [57]. Zhang and Tang [58] argue that the reduction of China's emissions mainly depends on the improvement of technical efficiency, whereas export changes will have little effect on emissions reduction. Likewise, the effect of Chinese exports on embodied air pollutants has varied in various sectors [23].
The findings of the study are similar and consistent to those of [10] that has been reviewed earlier, in which the authors examine the effect of trade liberalization on air pollution in 279 Chinese cities. They find that trade is negatively correlated with air pollution. Their findings are consistent with our study. Compared with [10], our results have a greater spillover effect, since we use a larger dataset at the provincial level, but they used city-data.
The findings provide clear evidence that substantiates the relationship between exports per capita and PM2.5 concentrations across the highest air pollution provinces in China.
Therefore, air pollutants in China would be overestimated without considering exports per capita. The difference between our findings concerning the previous studies, the findings are robust in dealing with controls for spatial autocorrelation and provides strong support for the hypothesis that increases in exports per capita will reduce PM2.5 emissions, this might be explained by numerous possible reasons among them. Firstly, an improvement in the open trade at the province level will lead to improving environmental quality via a decrease of PM2.5 emission, the spatial dependence model indicate an improvement in pollution levels by higher exports per capita. Secondly, in recent, the type of export products of high emission in China gradually changes to low emission and high value-added. In addition, maybe due to China's strict regulations against pollution-intensive in-dustries. Third, after entered in the WTO, the Chinese government does not encourage the establishment of high-emission firms; however, to reach international standards it instructed the firms to reduce pollution emissions.
On the other hand, the feedback effect of the population density shows that the population aggregation is a significant factor in local air pollution and has a positive influence on PM2.5 emissions in the neighboring provinces, maybe due to the dramatic expansion of capital cities of the provinces with high population density which will emit more urban pollutants.
The findings also indicate that PM2 Moreover, the estimated coefficient of direct effects of energy intensity is negative and statistically significant at 10%, implying that the effect of a 1% growth in local energy intensity will reduce 0.0991% in local PM2.5 emissions. This implies that a province with higher energy intensities tends to have lower PM2.5 emissions. This finding is consistent with Luo et al. [3] who argued that cities in Northwest China and Inner Mongolia provinces have low PM2.5 emissions with high energy intensities. This result also implies that China's technological innovation and the associated increase in energy intensity contribute to reducing PM2.5 concentrations. Thus, as noted above, energy intensity is cleaner than exports per capita.
Furthermore, both the direct effect values and the spillover effects estimation of the structure of energy consumption are insignificantly different from zero, which indicates that a specific province's energy consumption barely affects local and other provinces' emissions. In other words, the coefficient of the structure of energy consumption indicates that a 1% reduction will lead to a 0.122% reduction in emission intensity, this means the increase in the ratio of coal consumption to total energy consumption will lead to an increase in PM2.5 concentration. However, the problem with China's power production indicates that the proportion of coal is very high while the proportion of natural gas is very low.
When the proportion of coal used in energy production for a long time remained high, this means that consume more coal, and produce large amounts of dust, minute particles, and carbon dioxide which further aggravate urban smog, a finding which is consistent with theoretical expectations and it is supported by many works [43] [50]. Additionally, replacing coal consumption with clean energy sources is an effective mechanism to decrease PM2.5 emissions.
Similarly, the private vehicle also has a positive impact on PM2.5 concentrations on the direct effect, whereas the coefficient of the spillover effect is negative and significant. This means that own-province private vehicles increases will increase its own PM2.5 emission, but decrease the emission intensity of neighboring provinces. This outcome is consistent with [60]. Thus, increases in the number of vehicles will induce to consume more vehicle fuel, which is reflected in motor vehicle exhaust. The private vehicle in China has risen dramatically, in 2015; the average growth rate of private vehicles owned in Chinese provinces was 12.42% [11].

Conclusion and Policy Implications
This study has detected the spatiotemporal variations of export per capita and the main factors that drive the scale and intensity of PM2.5 emissions in the highest concentrations of air pollution provinces in China. The study used spatial econometric methods and extended the STIRPAT model. Based on the above results and discussion, important findings and conclusions are as follows.
The provinces in the northeast, south-central China, and the North China regions have more potential to reduce PM2.5 concentrations than the other provinces, and the Moran I index displays a declining tendency in the sample period indicating that the spatial autocorrelation is strengthened. Therefore, the main way to reduce emissions is by narrowing the difference between the province's emissions. This finding suggests that it is necessary for policy-makers should address provincial differences, and the policies issued should be specific provincial, exactly in the highest emissions provinces to reduce the emissions.
The results show that an increase in own-province exports will lead to reducing the PM2.5 emissions of both own province and neighboring provinces respectively. The results also provided evidence to verify the hypothesis of inverted U-shaped PMKC curves. Although for implementation of the Act for controlling automobile exhaust emission, private vehicles are still an important source of PM2.5 concentrations, therefore, the government should be aware. Thus, in light of these findings, and to move China into the low emissions pathway, policymakers need to adopt policies which aim to reduce PM2.5 emissions.
Firstly, to maintain low levels of PM2.5 emissions, China should continue to extend opening-up policies for trade and change the comparative advantage of the trade-in favor of cleaner production, and promote interprovincial technology cooperation, including both production and emission control technologies.
Furthermore, to prohibit the provinces from becoming more emissions in the future, policymakers should implement more strict regulations, such as, imposing more advanced technological processes also will make it possible to suppress

Conflicts of Interest
The authors declare no conflict of interest.