key: cord-1045089-o84qu6jp authors: Zhao, Chenkai; Sun, Ying; Zhong, Yaping; Xu, Senhao; Liang, Yue; Liu, Shu; He, Xiaodong; Zhu, Jinghai; Shibamoto, Takayuki; He, Miao title: Spatio-temporal analysis of urban air pollutants throughout China during 2014–2019 date: 2021-05-14 journal: Air Qual Atmos Health DOI: 10.1007/s11869-021-01043-5 sha: 330788f270c5be06c1f4442375fe2df564b62efa doc_id: 1045089 cord_uid: o84qu6jp Air pollution control has become the top priority of China’s “green development” concept since 2013. The Chinese government has enacted a range of policies and statutes to control contaminant emissions and improve air quality. On the basis of the national air quality ground observation database, the spatial and temporal distribution of air quality index value (AQI), fine particulate matter (PM(2.5)), coarse particles (PM(10)), sulfur dioxide (SO(2)), nitrogen dioxide (NO(2)), carbon monoxide (CO), and ozone (O(3)) were explored in 336 cities throughout China from 2014 to 2019. AQI and most pollutants (except O(3)) decreased in concentrations from 2014 to 2019. In 2019, all cities except Henan reached the level 2 of the ambient air quality index, and six cities had a lower ambient air quality index and reached the level 1. Spatially, higher pollutant concentrations were concentrated in large city clusters, whereas the areas with high O(3) concentration were found across the country. Furthermore, central heating was shown to have a negative impact on air quality. The observed AQI value, PM(2.5), PM(10), SO(2), NO(2), and CO concentrations were highest in north and northwest China and Henan province in central China. The correlations among pollutants suggest that the main sources of pollutants are fossil fuel combustion, industrial production, and motor vehicle emissions. The influence of meteorological factors on air quality, long-distance transportation, and the transformations of pollutants should be explored in future research. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s11869-021-01043-5. It has been widely demonstrated that particulate matter (PM) and gaseous pollutants in the atmosphere can enter the human body through the respiratory system, digestive system, or other channels and subsequently affect health. Coarse particles (PM 10 ) and fine particulate matter (PM 2.5 ) increase the risk of respiratory diseases, such as chronic obstructive pulmonary diseases (COPD) , asthma (Alotaibi et al. 2019) , bronchitis, and lung cancer. Gaseous pollutants, such as sulfur dioxide (SO 2 ), could increase the risk of cardiovascular disease (Yap et al. 2019) . Moreover, ozone (O 3 ) is associated with adverse cardiovascular effects (Xia et al. 2018) and could increase daily mortality rates (Yin et al. 2017) . Rapid industrialization, urbanization, and rural modernization in recent years have caused severe air pollution in China. Many cities are facing serious environmental problems. In the first quarter of 2013, China suffered from persistent and severe air pollution (Huang et al. 2014) . Therefore, the Chinese government adopted a series of measures to improve air quality. In September 2013, the Air Pollution Prevention and Control Action Plan (APPCAP) was disseminated. The goal of this plan was to improve air quality in heavily polluted environments (Huang et al. 2018) . In January 2015, the environmental protection law, regarded as an important milestone in China's air pollution prevention and control, was officially implemented (Song et al. 2017) . In 2016, China joined the Paris Agreement on climate change and promised to reduce hazardous emissions from coal-fired power plants by 50% by 2020, establish automatic real-time air quality monitoring systems, develop new energy technology, and promote renewable energy policies (Tambo et al. 2016) . In December 2016, China issued the thirteenth 5-year eco-environmental protection plan, which promised that SO 2 and carbon monoxide (CO) concentrations in prefecture-level cities will reach the standard by 2021, whereas the PM and nitrogen dioxide (NO 2 ) concentration would be significantly decreased (Ministry of Ecology and Environment of the People's Republic of China 2012). In 2018, China enacted the environmental protection tax law to raise awareness of environmental protection within Chinese society (The National People's Congress of the People's Republic of China 2018) . Previous studies on the temporal and spatial distribution of pollutants in China have several limitations. One issue is that several studies are limited to a short experimental period. For example, Li ) examined air pollution data for 2015-2016 and applied path analysis to reveal that temperature, wind speed, and precipitation were major factors affecting air pollutant concentrations. Another issue is that several studies only focus on a single factor. For example, Zeng (Zeng et al. 2019 ) used observations from national air quality stations to analysis the spatial-temporal variation of pollutants and found that the SO 2 , NO x and PM concentrations decreased, but the O 3 concentration increased. Moreover, China has been shown to be in transition from SO 2 -dominated to NO x -and O 3 -dominated air pollution. Analyses with larger study area have also been reported. Fan (Fan et al. 2020) conducted a study on the spatial and temporal variations of the six criteria pollutants (PM 2.5 , PM 10 , SO 2 , NO 2 , CO, and O 3 ) in 300 cities during 2014-2018 and found that PM, SO 2 , NO 2 , and CO concentrations decreased significantly. Meanwhile, the North China Plain and Central-Western Xinjiang were the areas identified to be the most seriously polluted. China is a vast country with various disparities in areas such as economic development, climate type, landform, and resource endowment. These phenomena have resulted in large differences in pollution levels among the provinces and across the seasons Lu et al. 2019; Yao et al. 2019 ). Furthermore, extremely low winter temperatures in the northern region have led to a strong demand for central heating, which could be an important source of air pollution in winter. To explore the spatial and temporal variations characteristic of air pollution in China in the present study, air pollutant concentration data for 336 prefecture-level cities in China were examined for the 2014-2019 period. The purposes of our investigation were to (1) observe the trends in variation of the air quality index (AQI) value and PM and gaseous pollutant concentrations, (2) elucidate the discrepancies in the spatial distributions of AQI values and air pollutant concentrations, (3) compare the effects of different heating methods on air quality, and (4) analyze the major sources of air pollution in China. We intend to provide some comprehensive information on air pollution and contribute to the prevention and control of air pollution in China. To evaluate the overall air quality status in China, we analyzed ambient monitoring data of PM 2.5 , PM 10 , SO 2 , NO 2 , CO, and O 3 in 336 prefecture-level cities from January 1 st , 2014, to December 31 st , covering 31 provinces, autonomous regions, and municipalities in mainland China (except Hong Kong, Macao, and Taiwan). To present the regional variation of air quality, the 31 provinces were divided into 7 regions based on their natural geographical characteristics, i.e., Central, East, North, Northeast, Northwest, South, and Southwest. The cities, regions, and their location were illustrated in Fig. 1 . According to the heating requirement, the government divides China into a central heating area and a non-central heating area by the Qinling Mountains-Huaihe River Line (Fig. S1 ). The heating season usually lasts from late autumn to early spring in central heating area. To estimate the air pollution caused by heating emissions, 14 cities were selected as the sample cities, with 7 located in the central heating area and the rest located in the non-central heating area. The sample cities are selected based on the following criterions: (a) The cities are well known in each region; (b) the cities can represent the meteorology distribution characteristic of each region; (c) the cities are the reprehensive of different levels of average pollutant concentrations in China; and (d) each geographical area has 2 sample cities. The names of sample cities were listed as follow: Beijing and Shijiazhuang in the North region, Shenyang and Changchun in the Northeast region, Shanghai and Nanjing in the East region, Zhengzhou and Wuhan in the Central region, Guangzhou and Haikou in the South region, Chengdu and Chongqing in the Southwest region, and Xi'an and Urumqi in the Northwest region. Details about 14 sample cities were shown in Table S1 . Daily AQI value and the real-time hourly concentrations of PM 2.5 , PM 10 , SO 2 , NO 2 , CO, and 8-h averaged peak O 3 (8-h O 3 ) were downloaded from the website of Ministry of Ecology and Environment of the People's Republic of China [http://113.108.142.147:20035] . Since the new technical regulation for ambient air quality index (HJ633-2012) (Ministry of Ecology and Environment of People's Republic of China 2012) has been executed, the air quality data of six criteria pollutants at individual monitoring site for major cities have been published through the website. To date, the monitoring sites have covered Mainland China, while the data of six criteria pollutants in 161 cities were lacked from January 1 st , 2014, to December 31 st , 2015 because many monitoring sites had not been established before 2015. According to the technical regulation (HJ633-2012), AQI is divided into six levels. The degree of air quality was graded on the following scale: 1, excellent; 2, good; 3, mild air pollution; 4, moderate air pollution; 5, heavy air pollution; and 6, severe air pollution. "Excellent" was defined as AQI from 0 to 50; "good" was defined as AQI from 51 to 100; "mild air pollution" was defined as AQI from 101 to 150; "moderate air pollution" was defined as AQI from 151 to 200; "heavy air pollution" was defined as AQI from 201 to 300; and "severe air pollution" was defined as AQI more than 300. We first calculated the annual average value or concentrations of AQI and six criteria pollutants (PM 2.5 , PM 10 , SO 2 , NO 2 , CO, and O 3 ) for the whole country, the seven regions, and the 31 provinces. Analysis of variance (ANOVA) was used to compare the differences in annual average concentrations on the countrywide and seven regional scales. Tukey HSD was used for multiple comparison. Based on the daily AQI values and six criteria pollutants concentrations of each city, a series of time-series boxplots were drawn to describe their temporal distribution. In order to better present the data distribution and probability density, a violin chart was superimposed on each month's box charts. A fitting line, which reflects the trend of these values over time, was added to each plot. The fitting method was generalized Subsequently, the annual average concentrations of six criteria from the 31 provinces were used to draw a series of statistical maps to describe the spatial distribution over time. The annual scale was unified to better present the variation trends in air quality. To investigate the impact of heating emissions on air quality, the Levene's test for the homogeneity of variance was used first in 14 selected sample cities, after which a multivariate analysis of variance (MANOVA) was used to analyze the influence of heating method (with or without central heating) and heating period (heating period and non-heating period) on air quality. Finally, a two-sample t-test was performed in 14 cities. The results were presented in a bar-plot with error bar. The last part of our study was to analyze the source of pollutants. First, a Shapiro-Wilk normality test was conducted for particulate matters and gaseous pollutants in each city. For those data that did not accord with the normal distribution, the Spearman rank correlation was used to calculate the correlation coefficient among variables. General linear correlation was performed for those data that accord with a normal distribution. Finally, the correlation matrix was drawn. The statistical tests were two-sided, and p values < 0.05 were considered statistically significant; in the Levene's homogeneity test, p value < 0.10 was considered statistically significant. All analyses and pictures were performed or drawn using R software (Version 3.6.1) (R Core Team 2020). On the basis of individual ANOVA analysis, the annual average variations of six criteria pollutants changed significantly (p < 0.05 prior to adjustment for multiple comparisons) over the 6-year period (Table 1) . Compared with 2014, the AQI value and the PM 2.5 , PM 10 , SO 2 , NO 2 , and CO concentrations decreased in 2019 by 28.7%, 38.0%, 34.1%, 67.9%, 24.9%, and 34.4%, respectively. The O 3 concentration increased significantly from 2014 to 2018 and decreased slightly in 2019. Overall, the O 3 concentration increased 12.4% between 2014 and 2019. The annual average variation of AQI values and the six criteria pollutants in the seven geographic regions are presented in Fig. 2 , and more detailed information were summarized in Table S2 . The variation patterns were similar among the seven regions; however, some regions showed some differences; for example, the PM concentration in northeast and central China increased significantly in 2018-2019 (p < 0.05), and the O 3 concentration in south China declined significantly from 2014 to 2015 and then increased from 2015 to 2019. The monthly average variation trends of the AQI value and six pollutants in the 14 sample cities are shown in Fig. 3 . Those of the remaining cities are shown in Figs. S2-S35. The AQI value and pollutants other than O 3 demonstrated significant seasonal variations with high values in winter and low values in summer. However, the seasonal variation of O 3 and NO 2 in Haikou was contrary to those in other cities; that is, NO 2 concentration was higher in summer and O 3 concentration was higher in winter. Fig. 4) , which marks the striking difference in the distribution of the country's population. There was a notable reduction of AQI value and PM 2.5 , PM 10 , SO 2 , NO 2 , and CO concentration over the study period in large city clusters (the Beijing-Tianjin-Hebei region (BTH); Yangtze River Delta (YRD); Pearl River Delta (PRD); Sichuan-Chongqing (CY); and Liaoning, Jilin, Shanxi, Shaanxi, and Gansu provinces), but the pollutant concentrations were still the highest compared with other areas in the same period. It is worth noting that the PM 10 concentration in Xinjiang Uygur increased significantly (p < 0.05) from 2015 to 2018. In Shandong, Henan, Jiangsu, Anhui, Shanghai, Qinghai, and Tibet, the O 3 concentration was higher than that of other provinces in the same period. Summaries of the annual average value of AQI and six criteria pollutants in the 31 provinces according to the Technical Regulation on Ambient Air Quality Index (HJ633-2012) were shown in Table S3 Table 2 presents the MANOVA results of two main effects: the central heating and heating period. The main effect of central heating showed that the AQI value and the PM 2.5 , PM 10 , SO 2 , NO 2 , and CO concentrations were significantly higher in areas with central heating than in areas without central heating (p < 0.05), while no difference in O 3 concentration was found between the two areas. The main effect of heating period demonstrated that the AQI value and the concentrations of the six criteria pollutants in the heating period were significantly higher than those in the non-heating period (p < 0.05). Furthermore, the heating period and central heating of the AQI and the six pollutants had significant interaction effects (p < 0.05). Figure 5 presents the results of the two-sample t-test, which compares the AQI values and concentrations of the six criteria pollutants between the heating and non-heating periods in 14 sample cities from 2014 to 2019. The difference between the heating period and the non-heating period in the central heating area was greater than that in the area without central heating. For example, the SO 2 concentration during the heating period in Changchun was 4.1 times that during the non-heating period, whereas the increase in Wuhan was only 0.7 times. Correlations between the five air pollutants Figure 6 shows the correlations between the five pollutants (PM 2.5 , PM 10 , SO 2 , NO 2 , and CO) in 14 sample cities. The O 3 present in the lower atmosphere is mainly generated by photochemical reactions. Therefore, relationships between O 3 and other pollutants were not independent, so the O 3 was excluded from the correlation calculations. In all 14 cities, PM 2.5 was highly correlated (R 2 > 0.5) with PM 10 , suggesting that both PMs came from the same or similar sources. Concentrations of both PM 10 and PM 2.5 were highly correlated with SO 2 , NO 2 , and CO in most cities, with R 2 of 0.54-0.72, 0.14-0.80, and 0.45-0.87, respectively. The correlation coefficients of SO 2 -NO 2 , SO 2 -CO, and NO 2 -CO were within the ranges 0.24-0.76, 0.35-0.82, and 0.21-0.82, respectively, with highly (R 2 > 0.5), moderately (0.25 ≤ R 2 ≤ 0.5), or weakly (0 ≤ R 2 ≤ 0.25) correlated. The correlation between SO 2 and CO in cities located in south China, such as Haikou, was weak or moderate, suggesting that coal combustion and industrial manufacture had a relative low contribution to air pollution in this region. Temporal variation of air pollutants is related to a series of measures taken by the Chinese government, such as shrinking the growth of harmful gas emissions, expanding production of clean energy, and continuing development of industrial structure policies and reforms (Tambo et al. 2016) . Recent research ) has indicated that the coal cap policy (a series of documents and policies released by China's State Council to support the reduction of coal consumption) resulted in significant emission reductions of PM 2.5 , PM 10 , SO 2 , NO 2 , and CO, which is also consistent with the results of the present study. For comparison, only four states in the USA failed to reach the PM 2.5 concentration standard, which was 12 μg/m 3 in 2015 (Bennett et al. 2019) . In a study covering nine countries and 16 provinces in Europe, the concentration of PM 2.5 ranged from 8.1 to 43 μg/m 3 (He et al. 2018) . Therefore, compared with developed countries, China still has higher concentrations of PM and gaseous pollutants, and thus the government needs to continue to pay attention to air pollution in the future. C h a n g c h u n S h e n y a n g B e i j i n g S h i j i a z h u a n g Z h e n g z h o u X i ' a n U r u m q i N a n j i n g S h a n g h a i W u h a n C h e n g d u C h o n g q i n g G u a n g C h a n g c h u n S h e n y a n g B e i j i n g S h i j i a z h u a n g Z h e n g z h o u X i ' a n U r u m q i N a n j i n g S h a n g h a i W u h a n C h e n g d u C h o n g q i n g G u a n g C h a n g c h u n S h e n y a n g B e i j i n g S h i j i a z h u a n g Z h e n g z h o u X i ' a n U r u m q i N a n j i n g S h a n g h a i W u h a n C h e n g d u C h o n g q i n g G u a n g C h a n g c h u n S h e n y a n g B e i j i n g S h i j i a z h u a n g Z h e n g z h o u X i ' a n U r u m q i N a n j i n g S h a n g h a i W u h a n C h e n g d u C h o n g q i n g G u a n g C h a n g c h u n S h e n y a n g B e i j i n g S h i j i a z h u a n g Z h e n g z h o u X i ' a n U r u m q i N a n j i n g S h a n g h a i W u h a n C h e n g d u C h o n g q i n g G u a n g with central heating without central heating C h a n g c h u n S h e n y a n g B e i j i n g S h i j i a z h u a n g Z h e n g z h o u X i ' a n U r u m q i N a n j i n g S h a n g h a i W u h a n C h e n g d u C h o n g q i n g G u a n g Changes in monthly averages of AQI, PM 2.5 , PM 10 , CO, SO 2 , and NO 2 conformed to a U-shaped pattern, with the highest values in late autumn and winter (November to February) and the lowest values in summer (June to August), which is consistent with the results of previous studies Ji et al. 2019) . High concentrations of pollutants in cold seasons can be explained as follows: (a) Central heating located in the north areas of the Qinling Mountains-Huaihe River Line leads to the increase in coal consumption. On the basis of the Design Code of Heating Ventilation and Air Conditioning in Civil Buildings (GB 50736-2012) (Ministry of Housing and Urban-Rural Development of the People's Republic of China 2012), central heating is only available in 13 provinces located north of the Qinling Mountains-Huaihe River Line. Residents living in Jiangsu, Anhui, and other provinces without central heating but with lower winter temperatures generally use air conditioners or individual heating equipment. Residents living in the south provinces, such as Guangdong and Hainan, do not need any heating equipment. The heating energy demands up to 8-12 kg of coal equivalent each winter when in central heating area. Only approximately 5 kg of coal equivalent is required when individual heating equipment was used (Shi et al. 2018) . Furthermore, the actual energy-saving effect in central heating is only 60.04% of the theoretical energysaving value (Lin and Lin 2019); (b) bad weather conditions, such as temperature inversions and low-powered air convection in winter, may hinder the diffusion and dilution of atmospheric pollutants (Zhan et al. 2017 ); (c) the higher atmospheric mixed layer and increased precipitation in the summer facilitated pollution dilution and deposition of atmospheric pollutants Yao et al. 2019) . However, the washout effect was affected by the precipitation threshold (Guo et al. 2016) . And it is difficult to reach the precipitation threshold due to the decrease of precipitation in winter (Fan et al. 2020 ); (d) the reduction in combustion efficiency caused by the cold start of vehicle engines in winter will also lead to increased emissions (Kim and Kim 2020) . Previous studies have found massive spatial differences in air quality and pollutant concentrations in China (Wang et al. 2017a; Ye et al. 2018) , and similar differences were also found in this study. This difference may be related to the transformation of economic development models in these regions (Lin and Zhu 2020) . Furthermore, it may also be related to the effectiveness of the implementation of environmental conservation and emission reduction policies in these regions (Song et al. 2020) . The BTH and northwest regions have higher concentrations of PM 2.5 and PM 10 . The BTH region is one of the most important economic centers in China. For the key industry in the BTH region, there were more than 6000 industrial boilers and a large number of metal basic industrial and nonmetal product sectors in 2017 . Moreover, there were over 5.5 million vehicles in Beijing Wang et al. 2017c) . High PM concentrations in northwest China were mainly from dust and aerosols in drylands . Sandstorms were also reported to have a major impact on the concentration of PM 10 in the northwest (Zhang and Cao 2015; Guan et al. 2018) . Furthermore, Miller-Schulze (Miller-Schulze et al. 2015) showed that sandstorms have no obvious effect on PM 2.5 . Therefore, PM 10 pollution was considered more serious than PM 2.5 in northwest China. Conversely, in spring and summer, the dust was transmitted to southwest through the Hexi Corridor by strong wind . As a result of the government's afforestation intervention measures, the conversion of farmland to forest, and reduction of desertification, air quality in northwest China has improved significantly (Du et al. 2019) . The SO 2 , NO 2 , and CO concentrations were higher in the north, northeast, and east regions and Xinjiang province. Moreover, the concentration of pollutants in some interior provinces was consistently higher than those of many coastal provinces. This was partly due to the transfer of emissionsintensive products, such as raw materials, from interior provinces to coastal areas. Furthermore, industries that produce pollutants tend to move to areas with less stringent pollution regulations and relatively low labor costs Sun et al. 2019) . Even when producing the same products, interior provinces usually have higher emission intensities than those of coastal provinces. For example, the coal-fired electric power industry generates 58.4 g SO 2 per US$ of electricity in Inner Mongolia, which is nearly twice that in Shanghai (Wang et al. 2017a) . Early research suggests that the production of PM and SO 2 from coal combustion in rural areas is higher than that in urban industrial areas (Zhi et al. 2017) . Therefore, the government should give more attention to underdeveloped areas in the next stage of their attempts to improve air quality nationally. Additionally, pollutants can be transported from the north and central regions to the east region through long-range meteorological transportation (Sun et al. 2019; Wang et al. 2019) . Ozone is a secondary pollutant generated by the photochemical reaction of nitrogen oxides (NOx) and volatile organic compounds (VOCs). Vehicle exhaust emissions and burning of biomass have been recognized as important sources of O 3 (Fang et al. 2019) . Areas with high NOx and VOC emissions, such as Shanghai and Guangzhou, or areas with high ultraviolet radiation, such as Tibet and Gansu, usually have higher O 3 concentrations (Wang et al. 2017b) . In a study on ozone sources and concentration variation in China, Liu proposed that large reductions in PM and SO 2 emissions also lead to increases in urban O 3 due to the complex effects of aerosols on radiation and chemical reactions, while meteorological factors also increase O 3 concentrations in the southwest (Liu and Wang 2020a; Liu and Wang 2020b) . Notably, some studies have found that simply reducing the concentration of NOx will lead to the increase of O 3 concentration. During the COVID-19 blockade, due to the shutdown of transportation industry, the reduction of NOx leads to a lower O 3 titration by NO, which will lead to the increase of O 3 concentration (Sicard et al. 2020; Siciliano et al. 2020) . It should be noted that the O 3 concentration in Hainan province is higher in the cold season, which is contrary to the seasonal characteristics of O 3 concentration in other provinces. The extremely high precipitation in summer in Hainan province not only weakens the solar radiation and thus reduce the generation of O 3 , but also precipitation scavenges and washes out the environmental O 3 concentration . Meanwhile, the NO 2 that originates from vehicle emissions in the PRD in winter forms NO 2 pollution events in Haikou city, and the O 3 concentration in Haikou city is highly correlated with the NO 2 concentration . Therefore, we suspect that the high O 3 concentration in Hainan province in winter may also be related to the meteorological transport of NO 2 , but this conjecture requires further study. The correlation coefficient of the different PM sizes is an important indicator for analyzing the source of pollutants. PM 2.5 and PM 10 are highly correlated in most areas, suggesting that these two particles come from the same or similar sources. In China, the main PM sources are combustion (such as coal-fired power plants or heaters, vehicle exhaust, biomass incineration) and dust Huang et al. 2020; Yang et al. 2020) . Therefore, the Chinese government should further promote clean energy policies aimed at reducing PM emissions from coal combustion and industrial production. The SO 2 source can be classified into natural and anthropogenic sources with the latter estimated to account for more than 70% of SO 2 global emissions (Sun et al. 2019) . Fossil fuel combustion is the main anthropogenic source of SO 2 , and thus the strong correlation with PM suggests that the main source of these pollutants was the combustion of fossil fuels. The correlation between NO 2 and PM was relatively high in the present study because of the effects of vehicle exhaust emissions, fossil fuel combustion, and agriculture production (Krotkov et al. 2016; Guo et al. 2019) . The results are consistent with the findings of the previous research ). The correlation between SO 2 and CO or NO 2 was mainly affected by coal-fired fossil fuels. The strong correlation indicates that combustion has a great influence on air quality, which is also consistent with Lu's report (Lu et al. 2019 ). The correlation between PM (both PM 2.5 and PM 10 ) and other gaseous pollutants in Shenyang and Urumqi was different from that in other cities, which was due to the extreme polluted weather record in these two cities. Taking Shenyang as an example, the concentration of PM 2.5 was recorded as 848 μg/ m 3 on November 8, 2015; however, other gaseous pollutants levels were low, so the unique correlation characteristics were generated. This study had several limitations. First, the construction of environmental monitoring stations was not fully completed in 2014, especially in the provinces of northwestern and southwestern provinces. Therefore, the representativeness of the study in 2014 in these areas may lower than those in the follow-up study. Second, this study focused on air quality in a single city or region, and the influence of meteorological, long-distance transportation, and the transformation of pollutants have not been explored. Therefore, to further promote the improvement of China's environmental conditions in the future, it is necessary to explore the influence of pollutant transmission in the atmospheric environment and to understand the influence of physical and chemical factors on air pollutant concentrations. Owing to strict control measures taken by the government, the air quality in all cities has improved significantly, but the O 3 concentrations have increased gradually. Furthermore, the spatial distribution of pollutants was variable, and the air quality of north, central, northeast, and east regions were worse than those of other regions. The central heating used in areas located in the north of the Qinling Mountains-Huaihe River Line has had a negative impact on air quality. Furthermore, the main sources of air pollutants in China were fossil fuel combustion, industrial production, and motor vehicle emissions. The air pollution in China is caused by multiple pollutants and thus shows great divergence among different regions. As a result, region-based management measures must be established. Furthermore, the influence of meteorological factors on air quality, long-distance transportation, and the transformation of pollutants should be further explored in future research. Data availability statements The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request. Competing interests The authors declare no competing interests. Traffic related air pollution and the burden of childhood asthma in the contiguous United States in 2000 and 2010 Particulate matter air pollution and national and county life expectancy loss in the USA: a spatiotemporal analysis Emission, transport, and radiative effects of mineral dust from the Taklimakan and Gobi deserts: comparison of measurements and model results Core Team R (2020) R: A language and environment for statistical computing. R Foundation for Statistical Computing A framework for investigating the air quality variation characteristics based on the monitoring data: case study for Beijing during 2013-Air Qual Atmos Health Removal of PM2. 5 and secondary inorganic aerosols in the North China Plain by dry deposition A comprehensive analysis of the spatiotemporal variation of urban air pollution in China during Rapid increase in ozonedepleting chloroform emissions from China Spatiotemporal variation of no 2 and sub-regional transport during winter pollution events in Haikou Spatial-temporal variations and mineral dust fractions in particulate matter mass concentrations in an urban area of northwestern China The washout effects of rainfall on atmospheric particulate pollution in two Chinese cities Air quality improvement and health benefit of PM 2.5 reduction from the coal cap policy in the Beijing-Tianjin-Hebei (BTH) region Spatial and temporal variations of air quality and six air pollutants in China during 2015-2017 Fine particles at a background site in Central France: chemical compositions, seasonal variations and pollution events High secondary aerosol contribution to particulate pollution during haze events in China Health impact of China's Air Pollution Prevention and Control Action Plan: an analysis of national air quality monitoring and mortality data Association between chronic obstructive pulmonary disease and PM2. 5 in Taiwanese nonsmokers Source apportionment and health risk assessment of air pollution particles in eastern district of Chengdu Spatial distribution differences in PM2.5 concentration between heating and non-heating seasons in Beijing Physical and chemical mechanisms of the daily-to-seasonal variation of PM10 in Korea Aura OMI observations of regional SO2 and NO2 pollution changes from Air pollution characteristics in China during 2015-2016: spatiotemporal variations and key meteorological factors Satellite-based estimation of full-coverage ozone (O3) concentration and health effect assessment across Hainan Island The actual heating energy conservation in China: evidence and policy implications Policy effect of the clean air action on green development in Chinese cities Worsening urban ozone pollution in China from 2013 to 2017 -Part 1: the complex and varying roles of meteorology Worsening urban ozone pollution in China from 2013 to 2017 -Part 2: the effects of emission changes and implications for multi-pollutant control Differences in concentration and source apportionment of PM2.5 between 2006 and 2015 over the PRD region in southern China Seasonal contribution of mineral dust and other major components to particulate matter at two remote sites in Central Asia Technical regulation on ambient air quality index Ministry of Ecology and Environment of the People's Republic of China (2012) Ambient air quality standards Ministry of Housing and Urban-Rural Development of the People's Republic of China (2012) Design code for heating ventilation and air conditioning of civil buildings Statistical communique of the People's Republic of China A high-resolution air pollutants emission inventory in 2013 for the Beijing-Tianjin-Hebei region Effects of climate and potential policy changes on heating degree days in current heating areas of China Amplified ozone pollution in cities during the COVID-19 lockdown Increased ozone levels during the COVID-19 lockdown: analysis for the city of Rio de Janeiro Air pollution in China: status and spatiotemporal variations Does the expansion of the joint prevention and control area improve the air quality?-evidence from China's Jing-Jin-Ji region and surrounding areas Spatio-temporal variations and factors of a provincial PM 2.5 pollution in eastern China during 2013-2017 by geostatistics Tackling air pollution and extreme climate changes in China: implementing the Paris climate change agreement The National People's Congress of the People's Republic of China Trade-driven relocation of air pollution and health impacts in China Ozone pollution in China: a review of concentrations, meteorological influences, chemical precursors, and effects Inter-regional and sectoral linkage analysis of air pollution in Beijing-Tianjin-Hebei (Jing-Jin-Ji) urban agglomeration of China Spatial-temporal evolution of PM2. 5 concentration and its socioeconomic influence factors in Chinese cities in 2014-2017 Effects of personal short-term exposure to ambient ozone on blood pressure and vascular endothelial function: a mechanistic study based on DNA methylation and metabolomics Comparison of ground-based PM2.5 and PM10 concentrations in China, India, and the US Effectively controlling hazardous airborne elements: insights from continuous hourly observations during the seasons with the most unfavorable meteorological conditions after the implementation of the APPCAP Properties of particulate matter and gaseous pollutants in Shandong, China: daily fluctuation, influencing factors, and spatiotemporal distribution Particulate air pollution on cardiovascular mortality in the tropics: impact on the elderly Spatial-temporal patterns of PM2.5 concentrations for 338 Chinese cities Ambient ozone pollution and daily mortality: a nationwide study in 272 Chinese cities Stabilization of atmospheric nitrogen deposition in China over the past decade Air pollution reduction in China: recent success but great challenge for the future Spatiotemporal variations and driving factors of air pollution in China Fine particulate matter (PM 2.5) in China at a city level Transboundary health impacts of transported global air pollution and international trade Village energy survey reveals missing rural raw coal in northern China: significance in science and policy Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations