Abstract

As an important field for human activities, cities play a critical role in PM2.5 reductions. Among the determinants for PM2.5 concentration, technological progress is considered to exhibit significant inhibitory effects. Although most extant research has focused on energy technologies or total factor productivity, due to limitations in data and methods, few scholars have focused on emission reduction technological changes at a city-level scale. Therefore, based on the combination of k-means clustering and the log-mean Divisia index method, this study estimates and explores the impact of PM2.5 emission reduction technology (PME) on the temporal changes and spatial differences of 262 Chinese cities’ PM2.5 concentration during 2003–2017. The findings show the following: (1) although the results based on econometric methods indicate that emission reduction technological changes decreased China’s city-level PM2.5 concentration, there were turning points in the yearly impacts, indicating that the improvements to emission reduction efficiency were not stable; (2) compared with PME, energy intensity played a more stable role in PM2.5 emissions reductions, implying that the improvement of energy efficiency was still very important in controlling PM2.5 concentrations; (3) based on the classified groups after clustering, most cities’ PME contributed to negative differences, but the PME of a small number of cities was very weak to largely lower the average level of their group; and (4) distributions of the spatial decomposition of the three classified groups were stable in the period of 2003–2017, implying that the catch-up and transcendence effects of PME within the group were limited. Thus, policymakers should focus on the impact of different policies on PME differences between cities.

1. Introduction

In recent years, PM2.5 pollution has seriously concerned countries around the world due to the threats it represents to human survival and to the sustainable development of society [13]. In light of previous studies, PM2.5 emissions can not only lead to high haze reduction costs but also indirectly cause health costs through medical expenses and work-time loss. In particular, the emissions control cost in Shanghai would be about 1.01%–2.26% of the corresponding gross domestic product (DDP) values in 2030 [4]. And the direct external cost of residential health issues caused by PM2.5 pollution in Beijing and Changsha was as high as 0.3%–2.69% of the GDP during 2012–2017 [5, 6].

As one of the world’s largest developing countries, China is undergoing a process of accelerating urbanization and industrialization, facing increasingly severe PM2.5 pollution [7, 8]. Thus, the government has focused on formulating a series of policy regulations, such as the Air Pollution Prevention and Control Law, switching from coal consumption to natural gas consumption, and setting emission reduction targets for PM2.5 concentration in the 12th Five-Year Plan (2011–2015) [9, 10]. At the same time, many actions have also been taken to monitor PM2.5. For example, China has established thousands of PM2.5 monitoring observatories, which have been reporting hourly and daily air quality since 2012. Additionally, China’s grassroots government has intensified supervision of straw burning in rural areas to curb haze pollution.

Concurrently, an increasing number of scholars have begun to pay attention to the field of PM2.5 pollution and have provided some corresponding policy implications with respect to curbing it. Given that PM2.5 concentration has mainly been driven by human activities, such as vehicle emissions, power generation, and industrial production, many studies have paid specific attention to the impacts of socioeconomic drivers on PM2.5 pollution [1114]. For example, Ma et al. [15] analyzed the relationships among GDP per capita, the price of refined oil, vehicle population, energy intensity, and PM2.5 concentration in 152 Chinese cities based on spatial linkage and found that the impacts of economic, social, and energy efficiency on different cities varied significantly. Similarly, Zhang et al. [16] used the log-mean Divisia index (LMDI) method to break down changes in the PM2.5 concentration in these 152 cities. Their results showed that the decline of PM2.5 concentration was mainly driven by energy intensity and emission intensity, while GDP per capita and population were responsible for the increasing PM2.5 concentration. Among the driving factors of PM2.5 concentration, many socioeconomic factors have been found to contribute to the increase in PM2.5 emissions and concentration. These factors, such as increasing urbanization, economic growth, and population growth, are difficult to reduce in order to decrease PM2.5 pollution. Therefore, several studies have focused on factors that directly contribute to reducing PM2.5 pollution.

Among many economic factors, technological progress has been considered to play an important inhibitory role [17, 18]. In particular, Chen et al. [18] studied the effects of foreign direction investment, export learning effects, research and development, and import technology on the reduction of PM2.5 concentrations in 48 Chinese cities during 2000–2015. Likewise, Li et al. [19] used spatial econometric models to study the potential relationship between environmental total factor productivity and 283 Chinese cities during 2000–2013. Similarly, Xu et al. [20] analyzed the impacts of total factor productivity on 281 prefecture-level cities in China from 2007 to 2017 based on a spatial dynamic panel model, concluding that technological progress has significant positive effects on alleviating PM2.5 concentration. In total, these studies have only considered energy technologies or total factor productivity, ignoring the effects of PM2.5 emission reduction technology (PME). There is no doubt that improvements to energy efficiency can reduce energy input and decrease PM2.5 emissions; however, there is a need to pay more attention to PME changes, since these determine the PM2.5 pollution produced under the same unit of energy consumption and economic output and directly influence the level of cities’ sustainable development.

In the light of existing literature, although scholars have studied the impacts of drivers on PM2.5 pollution in China’s prefecture-level cities from a socioeconomic perspective, the current literature on PM2.5 concentrations has several shortcomings. (1) Some scholars have analyzed the impacts of technological progress on PM2.5 pollution, but few studies have examined the impacts of PME changes on China’s city-level PM2.5 pollution; (2) many studies have concluded that there is significant spatial and temporal heterogeneity in the determinants of PM2.5 pollution based on spatial econometrics, but have not outlined detailed temporal changes; and (3) although some scholars have used the LMDI method to analyze determinants of temporal PM2.5 pollution changes, they have not considered spatial differences, especially within classified groups.

Hence, this study explores the impacts of PME changes on PM2.5 concentrations of 262 Chinese cities during 2003–2017. First, we applied three spatial econometric models (SDM, SAR, and SEM) to identify the spatial autocorrelation and potential relationship between PME changes and PM2.5 concentrations. Next, we decomposed the temporal and spatial changes in the PM2.5 concentrations of the 262 cities based on the combination of LMDI and k-means clustering methods. Our study makes the following contributions and findings: (1) we estimate the PM2.5 PME of 262 cities and analyze its impacts on China’s city-level PM2.5 concentration; (2) the results show that there have been turning points in the yearly impacts, indicating that improvements to the emission reduction efficiency have been limited and unstable; and (3) based on the combination of k-means clustering and spatial LMDI methods, we find that the gaps among the spatial decomposition in our classified three groups were stable in the period 2003–2017, implying that the ranking of the emission reduction technology development level of each prefecture-level city in each group was relatively stable.

2. Materials and Methods

2.1. Spatial Impacts of PME on PM2.5 Concentrations

In light of the modified production-theoretical decomposition analysis method proposed by Wang et al. [21], we first estimated the PME during 2001–2017. We selected three input factors, one desirable output, and one undesirable output: fixed capital stocks, employed population, energy consumption, real GDP, and PM2.5 concentrations. The calculation of PME was based on the Shephard distance functions:where represents fixed capital stocks, represents labor force, represents energy consumption, represents PM2.5 concentration, represents gross domestic output, and represents PM2.5 reduction efficiency.

Thus, the PME can be estimated using the following equation:where denotes emission reduction technological changes.

Based on previous studies [22], we further used spatial econometric methods to explore the potential spatial impacts of PME on PM2.5 concentration. The spatial models of Durbin (SDM), autoregression (SAR), and error (SEM) were selected. The normal form of spatial econometric model is given aswhere represents PM2.5 concentration, represents explanatory variables, is the spatial distance weight matrix, and is the error term. We used the reciprocal of the geographical distance of the city center to construct the spatial econometric weight.

2.2. The Temporal and Spatial LMDI Methods

Although spatial regression methods have the ability to reveal potential relationship between the expected values of PME and PM2.5 concentrations, they ignore temporal and spatial changes. Thus, we used the index decomposition analysis (IDA) method to analyze the yearly and spatial impacts of PME on PM2.5 concentration. In light of previous studies (e.g., [11, 16]), the IDA identity of PM2.5 concentrations can be described as follows:where denotes the ith city’s PM2.5 concentration in period , represents the ith city’s energy consumption in period , is the ith city’s real gross domestic productivity in period , and denotes the ith city’s population in period .

Based on a proposition by Zhou and Ang [23], we can obtain the following equation:

To simplify equation (5), it can be written as follows:where denotes the ith city’s potential emission intensity of PM2.5 concentration without technological change’s impacts in period ; represents the ith city’s energy intensity in period , which reflects energy efficiency; represents the ith city’s GDP per capita in period ; and represents the ith city’s PM2.5 emissions efficiency in period , whose changes reflect the changes in PME.

Next, considering that there may be strong spatial heterogeneity in the effects of PME, we adopted a spatial clustering decomposition analysis method in this study. As described by Cheng et al. [24], this method entails a combination of spatial index decomposition analysis and clustering methods. In particular, we adopted the LMDI method to study the impacts of the five driving forces on PM2.5 concentrations because it is robust, easily understood, and has clear economic meaning. Appendix A presents the formula of the temporal and spatial LMDI method, which can indicate the temporal and spatial impacts of the five driving forces on PM2.5 concentration.

2.3. The Cluster Method

Regarding traditional spatial LMDI method, it is necessary to select a reference for the spatial index decomposition analysis. Normally, the average of total samples’ may be used as a reference. However, it will ignore heterogeneity between and within groups. Thus, the clustering method was adopted in this study to combine with the traditional spatial LMDI method, better capturing differences between and within groups.

With regard to the clustering method, hierarchical and nonhierarchical cluster methods are the two main approaches used. The algorithms of hierarchical cluster analysis are related to the construction of the tree structure of clusters, while that of nonhierarchical cluster analysis aims to classify objects into a predetermined number of disjointed clusters. Thus, we use a nonhierarchical cluster to make classifications among the observed cities. In particular, we selected k-means clustering because the algorithm has exhibited high efficiency and empirical success [25, 26] and has been widely adopted in many fields of study, such as environmental analysis, image classification, and so on [24, 27]. However, the operations of k-means clustering analysis need to assume the number of center points in advance, leading to some uncertainty. To overcome this disadvantage, we first used hierarchical cluster analysis based on the shortest Euclidean distance to identify the dendrogram and classifications [28] so as to provide a more intuitive and reasonable reference for the number of center points selected for k-means methods.

Additionally, with regard to the cluster indicators, the cities’ PM2.5 concentrations and GDP in 2017 were considered. Given that the PM2.5 emission intensity (i.e., ratio of PM2.5 concentration to GDP), like carbon intensity, has often been used to represent the level of PME [16], PM2.5 concentrations and GDP would perform well in revealing the clustering and classifications of PME. Thus, we selected cities’ PM2.5 concentrations and GDP in 2017 for this purpose.

The resulting dendrogram based on hierarchical cluster analysis is presented in Figure 1. From the dendrogram, we can clearly observe the clustering process of each sample. In total, 262 cities can be divided into 28 categories based on the shortest Euclidean distance. From the right to left, the 18th category can be regarded as a single big category, and the 17th category can also be regarded as a separate category. In addition, the 23rd category and the remaining categories can be combined into one big category. Therefore, we speculated that the cities may be classified into three categories. Thus, we preset the number of center points as three. Next, based on the main idea of the k-means clustering algorithm, we continuously optimize the selection of the center points according to the Euclidean distance between each point and the center point, until the center points are stable.

2.4. Data

Using the equations above, city-level data were derived for PM2.5 concentrations, energy consumption, fixed capital stock, population, and real GDP. Given the availability of PM2.5 concentration and fixed capital stock data, this paper focused on 262 cities, and the research period spanned from 2003 to 2017.

The data for PM2.5 concentrations were obtained from the real-time monitoring of the Ministry of Ecology and Environment monitoring site (https://106.37.208.233:20035). We estimated the data for city-level energy consumption based on city-level CO2 emissions, as proposed by Chen et al. [29], due to the significant relationship between energy use and CO2 emissions. The economic output was taken from the China Statistical Yearbook (2003–2017). To avoid the influence of price, we used the price in 2001 as the constant price and obtained the real GDP, which is consistent with the approach used by Chen et al. [30]. Additionally, we calculated the fixed capital stock based on the perpetual inventory method, which is consistent with Chen et al. [31].

3. Results and Discussion

3.1. The Results of Spatial Regressions

Based on the spatial econometric method outlined in Section 2.1, we used fixed effects, spatial Durbin (SDM), autoregression (SAR), and error (SEM) models to analyze the relationship between PM2.5 concentrations and PME. In line with Hao et al. [32], we selected GDP and population as the control variables. Table 1 shows the results. Among them, P denotes population and Y denotes GDP.

Column (1) shows that there was a significant negative relationship between PM2.5 concentrations and PME, indicating that PME helpfully reduced PM2.5 concentrations. Columns (2)–(4) imply that there should be spatial heterogeneity. At the same time, the SDM and SAR models individually report the direct, indirect, and total effects of PME progress, indicating that the impacts of PME on different regions’ PM2.5 concentrations may vary by regions. In total, the regressions results indicate that China’s PME has generally reduced the cities’ PM2.5 concentration.

3.2. Drivers of the Temporal Changes in PM2.5 Concentrations

Based on the temporal LMDI method outlined in Section 2.2, we categorized the city-level changes in PM2.5 concentrations into five main driving forces that pertain to the influences of energy use. Given that there are no official data about city-level energy use, some scholars have used electronic power consumption to replace total energy use. However, such replacement may cause significant errors, since fluctuations in electronic power use are inconsistent with total energy combustion. Therefore, we calculated the energy use of the 262 cities during 2003–2017 based on data on city-level CO2 emissions, as proposed by Chen et al. [29].

Considering that the calculation of city-level CO2 emissions followed a top-down approach (i.e., provinces and cities), the corresponding city-level energy use can be estimated according to city-level CO2 emissions and the relationship between provincial CO2 emissions and energy consumption. To capture the provincial differences in the ratio of CO2 emission to energy use, we used the varied coefficient model to conduct a regression between provincial CO2 emissions and energy consumption during 1997–2017. At the same time, to avoid the negative values of city-level energy use, we used a no-constant model. Table 2 in Appendix A presents the results.

Evidently, the coefficient of determination was 0.986, implying that the estimated city-level energy consumption was highly accurate. Based on the city-level energy use, we can obtain the impacts of driving forces on the PM2.5 concentrations of the 262 cities.

To capture more temporal fluctuations of PME impacts, we divided the period of 2003–2016 into four parts: 2003–2006, 2006–2009, 2009–2012, and 2012–2017. These results are presented in Figure 2. They show that not all cities’ PME continuously decreased the PM2.5 concentrations during 2003–2017, and detailed information thereon cannot be found based on regressions. At the same time, there were many turning points in the impacts of PME changes on PM2.5 concentrations, which again cannot be captured by regressions.

Given that serious PM2.5 emissions always occurred in large cities, we selected PME impacts in Beijing, Shanghai, Guangzhou, and Shenzhen for analysis. In particular, Beijing’s PME impacts changed by about −0.59, −14.84, and −4.50 ug/m3 in the periods of 2003–2006, 2006–2009, and 2013–2017, respectively; however, it increased by about 4.49 ug/m3 during 2009–2013. PME changes reduced Shanghai’s PM2.5 concentrations from 2003 to 2006 (−0.43 ug/m3), while they increased by about 0.30, 9.30, and 8.18 ug/m3 during 2006–2009, 2009–2013, and 2013–2017, respectively. Guangzhou’s PME changes had a negative impact during 2009–2013 (−6.29 ug/m3), while it stimulated PM2.5 concentrations in the other periods (+4.34, +0.56, and +7.81 ug/m3). Shenzhen decreased its PM2.5 concentrations in 2003–2006 and 2009–2013 (−0.82 and −5.66 ug/m3), but increased them in the periods 2006–2009 and 2013–2017 (+2.68 and + 2.91 ug/m3, resp.). Evidently, Beijing fared better than the other three cities in developing PME, especially after 2013; this is consistent with Yang et al. [14] and Zíková et al.’s [33] findings. The reasons for this may derive from the issuance and implementation of more stringent policies related to smog control in the Beijing–Tianjin–Hebei region in 2012, such as “Beijing–Tianjin–Hebei regional environmental protection takes the lead in breaking through the cooperation framework agreement,” switching from coal consumption to natural gas consumption, and implementation of the Atmospheric Pollution Prevention and Control Action Plan [11, 34].

Figure 3 presents the impacts of the other four driving forces on PM2.5 concentrations. Potential emission intensity in 210 cities curbed PM2.5 concentrations, reflecting the optimization of the energy consumption structure. In particular, Jinchang (−42.93 ug/m3), Hengshui (−40.87 ug/m3), Anshan (−38.13 ug/m3), Daqing (−37.39 ug/m3), and Baoding (−33.40 ug/m3) were the top five cities whose potential emission intensity decreased PM2.5 concentrations during 2003–2017. Energy intensity had reduction effects on PM2.5 concentrations in all observed cities during 2003–2017, which is consistent with Chen et al. [11], Zhang et al. [16], and Li et al. [35]. The results indicate that the improvement in energy efficiency may be responsible for energy conservation, leading to the decline in PM2.5 concentrations. GDP per capita and population in almost all cities contributed to increasing PM2.5 concentrations during 2003–2017. Among them, Zhengzhou’s GDP per capita and Shengzhen’s population individually contributed to the highest increases, of 139.21 and 29.85 ug/m3, respectively.

3.3. PME for the Spatial Changes in PM2.5 Concentrations among Cluster Cities

Based on the above-described clustering methods, we divided the 262 cities into three groups. The first group was concentrated on the cities with the largest emissions intensity (0.026 ug/m3million yuan, averagely); the second group had the middle emissions intensity (0.006 ug/m3million yuan, averagely); and the third group’s average emissions intensity were the lowest (0.002 ug/m3million yuan). The detailed classification is presented in Appendix A, Table 3.

Then, we studied the impacts of PME changes on within-group differences of PM2.5 concentrations. Figure 4 reports the impacts of PME on PM2.5 concentrations between each city and the classified group’s average in 2003 and 2017.

Regarding the cities in the first group, 141 cities’ PME contributed to negative differences within the first group’s PM2.5 concentrations and reduced the within-group differences when the average PM2.5 concentrations in the first group in 2003 were set as the reference, indicating that these cities’ PME did better than their average level. And there were 75 cities’ PME stimulated within-group differences, implying that they failed to reach their groups’ average level. Taking six cities in 2003 as examples, Daqing (−176.31 ug/m3), Taizhou (−161.50 ug/m3), and Xiamen (−147.90 ug/m3) were the top three cities with negative differences. Meanwhile, Jiayuguan (+270.48 ug/m3), Jinchang (+152.6 ug/m3), and Tongchuan (+151.10 ug/m3) were the top three cities that contributed to positive within-group differences. Compared with 2003, there were several changes in the spatial decomposition results in 2017: 147 cities’ PME lowered the between-emissions differences within the first group’s PM2.5 concentrations, but PME changes in 69 cities caused a positive within-group gap. For example, Ankang city’s PME reduced the within-group gap in 2013 (−54.28 ug/m3) but failed in 2017 (+9.82 ug/m3).

In the second group, the PME in 20 cities led to reductions in the within-group differences in 2003, and 16 cities’ emission reductions caused positive within-group differences. Particularly, PME in the cities of Quanzhou (−51.10 ug/m3), Hangzhou (−50.30 ug/m3), and Fuzhou (−47.35 ug/m3) contributed to a decline in the within-group difference, while Hefei (+54.07 ug/m3), Taizhou (+52.57 ug/m3), and Jining (+51.09 ug/m3) saw an increase of the within-group gap. Moreover, there were 18 cities in which PME helped to reduce the within-group differences in 2017, implying that the spatial differences did not change much in these cities.

With regard to the third group, PME in the first-tier cities (Shanghai, −53.98 ug/m3; Guangzhou, −35.47 ug/m3; Shenzhen, −30.03 ug/m3; and Beijing, −14.17 ug/m3) contributed to reductions in the within-group differences in 2003. However, the PME of Tianjin, Suzhou, and Chongqing increased to 55.82, 24.16, and 19.30 ug/m3 of the within-group differences, respectively. In 2017, the four first-tier cities’ and Chongqing city’s PME changes were responsible for the negative within-group differences (Shanghai, −18.89 ug/m3; Guangzhou, −16.15 ug/m3; Shenzhen, −27.79 ug/m3; Beijing, −24.41 ug/m3; and Chongqing, −10.41 ug/m3), and the remaining cities still contributed to positive within-group gaps (Tianjin, +18.89 ug/m3; Suzhou, +35.19 ug/m3).

In summary, since most cities’ PME contributed to negative impacts on the within-group differences, indicating that most cities’ PME was better than their group’s average level. However, it also should be noted that the PME of a small number of cities, such as Jiayuguan city and Tongchuan city, were lower than their groups’ average level and had a very weak impact on curbing PM2.5 concentrations, which is consistent with the results proposed by Zhang et al. [16]. At the same time, the distributions of the spatial decomposition of the three groups were stable in the period of 2002–2017, implying that the catch-up effect within the group was not obvious. Thus, the governments should focus more on regional heterogeneity when encouraging the promotion of PM2.5 emission reduction technologies.

4. Conclusions

This study explored the impacts of PME changes on the PM2.5 concentrations of 262 Chinese cities during 2002–2017. First, we applied three spatial econometric models (SDM, SAR, and SEM), identifying the spatial autocorrelation and potential relationship between PME changes and PM2.5 concentrations. Furthermore, we decomposed the temporal and spatial changes in the 262 cities’ PM2.5 concentrations based on a combination of LMDI and k-means clustering methods. In terms of the temporal changes, we found that most cities’ PME impacts fluctuated during 2002–2017, and there were many turning points, indicating that the reduction effects in many cities were not stable. However, energy intensity continuously played a significant role in decreasing PM2.5 concentrations due to energy conservation caused by the improvement in energy efficiency. With regard to the spatial decomposition based on clusters, we found that most cities’ PME contributed to negative differences, implying that the PME in a small number of cities, such as Jiayuguan City and Jinchang City, was not strong and lowered the average level of their group. At the same time, the distributions of the spatial decomposition of the three groups were stable in the period 2002–2017, implying that a catch-up effect within the group was not obvious.

Based on the analysis, we provide several policy implications for curbing PM2.5 concentrations from a socioeconomic perspective. First, although the results, based on econometric methods, show that PME progress had significant curbing impacts on China’s city-level PM2.5 concentrations, the temporal and regional heterogeneity cannot be ignored. Given that the PME failed to continuously decrease China’s city-level PM2.5 concentration, China’s cities should pay more attention to increased fiscal expenditures in technological innovations, such as composite carbon slurry clean combustion technology, supercritical flue gas desulfurization (FGD) systems, and ultra-supercritical FGD systems [21, 35].

Second, in addition to increasing PME, it is important to promote the application of energy-saving technologies and the optimization of energy use structures. In light of the effects of energy intensity on PM2.5 concentrations, the promotion of energy efficiency will reduce the energy input, thus reducing PM2.5 emissions and concentrations under the same conditions of PME. In addition, since PM2.5 pollution is mainly driven by fossil fuel combustion [36], speeding up the use of clean energy would also significantly reduce PM2.5 concentration.

Third, considering that the spatial impacts of PME were relatively stable over time, regional advanced emission reduction technological innovations should be promoted around and selected as demonstrations so that the catching-up effects within group will be highlighted.

Appendix

The temporal impacts of drivers on PM2.5 concentration can be estimated based on equations (A2.1)–(A2.5):where denotes the impacts of potential emission intensity on the ith city’s PM2.5 concentration from period to ; denotes the impacts of energy intensity on the ith city’s PM2.5 concentration from period to ; denotes the impacts of GDP per capita on the ith city’s PM2.5 concentration from period to ; denotes the impacts of population on the ith city’s PM2.5 concentration from period to ; and denotes the impacts of PM2.5 emission efficiency on the ith city’s PM2.5 concentration from period to .

Similarly, based on the definite reference, the spatial LMDI method can also be used for spatial decomposition; the corresponding formulas are presented as follows:where denotes the gap between the impacts of potential emission intensity on the ith city’s PM2.5 concentration and the reference; denotes the gap between the impacts of energy intensity on the ith city’s PM2.5 concentration and the reference; denotes the difference between effects of GDP per capita on the ith city’s PM2.5 concentration and the reference; denotes the gap between the effects of population on the ith city’s PM2.5 concentration; and the reference denotes the gap between the impacts of PM2.5 emission efficiency on the ith city’s PM2.5 concentration and the reference.

Data Availability

The data can be obtained upon request from the corresponding author.

Additional Points

PM2.5 emission reduction technology (PME) was estimated for 262 cities. The 262 cities were classified into three groups based on the k-means cluster method. The curbing impacts of PME on PM2.5 were fluctuated and unstable. Energy intensity played a more stable role in PM2.5 emissions reductions. The catch-up and transcendence effects of PME within the group were limited.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work was supported by the National Key Natural Science Foundation of China (Grant no. 71934001), the National Natural Science Foundation of China (Grant nos. 71471001, 41771568, 71533004, and 71503001), the National Key Research and Development Program of China (Grant no. 2016YFA0602500), Sichuan Province Social Science High-Level Research Team Building Program, and the Program for Major Projects in Philosophy and Social Science Research under China’s Ministry of Education (Grant no. 14JZD031).