Analysis of factors influencing air quality during COVID-19
Interannual differences of air quality
The average AQI value and concentration of the six pollutants observed in different periods are shown in Table S1. From Jan. 24 to Jun. 30, 2020, the average AQI value in Tangshan was 93.20, which corresponded to the second air quality level (good), and the average concentrations of PM2.5, PM10, SO2, NO2, CO, and O3-8h were 50.03 µg/m3, 93.04 µg/m3, 19.86 µg/m3, 41.36 µg/m3, 1.25 mg/m3, and 110.17 µg/m3, respectively, which were all below the Ambient Air Quality Standard Grade II (PM2.5 75 µg/m3, PM10 150 µg/m3, SO2 150 µg/m3, NO2 80 µg/m3, CO 4 mg/m3, O3-8h 160 µg/m3).
Compared with the data from Jan. 24 to Jun. 30, 2017–2019, the daily average concentration of all pollutants in 2020 decreased to varying degrees between 6.13% and 42.78%, with SO2 and O3-8h showing the largest (42.78%) and smallest (6.13%) reductions, respectively. There were significant differences in the concentrations of PM2.5 at the 95% confidence interval between 2017–2019 and 2020, and significant differences in the concentrations of PM10, SO2, NO2, and CO at the 99% confidence interval. No significant decreased O3 in Tangshan was mainly attributed to the enhanced atmospheric oxidation capacity (Zhu et al., 2021). During the same period in 2021, the concentrations of PM2.5, PM10, and NO2 increased slightly compared with those observed in 2020, potentially meaning that the Level I and II response had a greater impact on air quality than III response. While the SO2 and CO values in 2021 were lower than those in 2020, which were mainly related to the continuous implementation of the “coal to gas” and “coal to electricity” policies (Liu et al. 2020).
The 2019–2021 period was selected as the study scope to analyze the temporal trend of six atmospheric pollutant concentrations (Fig. 3). PM2.5, PM10, SO2, NO2, and CO concentrations showed similar variation patterns, suggesting that they might have common sources and influencing factors. O3 is a secondary pollutant affected by the combined action of nitrogen oxides (NOx) and volatile organic compounds (VOCs) (Li et al. 2019; Roy et al. 2021), and due to the complexity of its mechanism of generation, it shows a different pattern of variation compared to the other pollutants. Overall, the range of variation in the concentration of each pollutant in 2020 was smaller than that in 2019 and 2021, and the number of clean air days was significantly higher. The proportion of days in which PM2.5 concentration exceeded the standard (75 µg/m3) in 2019, 2020, and 2021 was 19.2%, 13.29%, and 25.95%, respectively, indicating that the health risks associated with the exposure to particulate matter were clearly lower in 2020 (Shen et al. 2017).
Analysis of significant difference based on DID model
This section explored the long-term effect of COVID-19 control measures on the air quality in 2020. The DID model used the data from 23 days (January 1) before the Level I response and 159 days (June 30) after the response. The descriptive statistical results of the main variables were reported in Table S2A. Table S3A reported the estimation results of model. The coefficient \(\alpha\) was found statistically significantly negative at the 1% level, which indicated that air quality improved significantly during COVID-19 in Tangshan. The coefficient \(\delta\) significantly negative at the 1% level, which confirmed the previous hypothesis, i.e., Tangshan has a trend of improving air quality year by year.
However, the parallel trend test was conducted, and the result opposed the use of the DID model (Figure S1A). In 2017–2019 and 2021, the coefficients were not significantly different from 2020, which shown that the seasonal influences played an important role in improving air quality. In conclusion, the results of the DID model cannot demonstrate that the control measures have a significant impact on air quality after the beginning of the Level III response for a period of time.
Impact of meteorological factors
Meteorological factors have a strong influence on the variation of air pollution levels (Dantas et al. 2020; Shenfeld 2011). To explore the impact of these factors on air quality in 2020, the Spearman correlation analysis was conducted among AQI, six pollutants, and six meteorological parameters (i.e., temperature, humidity, wind speed, atmospheric pressure, visibility, and precipitation) using SPSS software. The data from Jan. 24 to Jun. 30, 2017–2021 was selected for the analysis because the use of longer time series could reduce the impact of uncertainties. The results were shown in Fig. 4: temperature was strongly positively correlated with O3 (rs = 0.842, P < 0.01), and negatively correlated with PM2.5, PM10, NO2, and CO; humidity was positively correlated with pollutants other than O3, and moderately positively correlated with PM2.5 and CO (rs = 0.441, P < 0.01; rs = 0.372, P < 0.01); wind speed was negatively correlated with pollutants other than O3, as wind could transfer, diffuse, and dilute pollution; air pressure was strongly negatively correlated with O3 (rs = − 0.682, P < 0.01), and was positively correlated with PM2.5, NO2, SO2, and CO; visibility was negatively correlated with PM2.5, PM10, NO2, and CO, and was more significantly correlated with particulate matter; precipitation was negatively correlated with pollutants other than CO, but the overall correlation was low.
From Jan. 24 to Jun. 30, 2020, the average temperature in Tangshan was 10.56°C, humidity was 54.43%, wind speed was 1.92m/s, atmospheric pressure was 1015 hPa, visibility was 16.09 km, and precipitation was 0.73 mm (Table S4). These were mostly unfavorable meteorological conditions compared to those recorded in 2017–2019 and 2021, indicating that meteorological factors in 2020 had a negative effect on the improvement of air quality. However, overall, the air quality in 2020 was better than that in 2017–2019 and 2021, showing that anthropogenic factors might played a positive role in improving the air quality and masked the impact of meteorological factors.
Impact of COVID-19 control measures
Based on the Air Pollution Source Emission Inventory of Tangshan (2017), it is clear that industrial processes, dust, fossil fuel combustion, and vehicles are the main sources of conventional air pollutants (Figure S2). In particular, the emission contribution rates of industrial processes to PM2.5, PM10, SO2, and CO were 53.55%, 38.03%, 78.16%, and 90.68% in 2017, respectively. The measures adopted to control the epidemic had an impact on the emission intensity of the above-mentioned pollution sources. Therefore, the available economic indicators closely related to them were selected as parameters reflecting the impact of the variation of anthropogenic emissions on air quality during the 2019–2021 COVID-19 period (Table 1).
Table 1
Main economic indicators for the city of Tangshan in the first and second quarters of the 2019–2021 period
Economic Indicator | Unit | 2019 | 2020 | 2021 |
1st | 2ed | 1st | 2ed | 1st | 2ed |
GDP of the Secondary Industry | Billion Yuan | 82.21 | 188.79 | 73.69 | 171.49 | 91.91 | 208.74 |
Total Profit of Manufacturing | Billion Yuan | 3.79 | 16.46 | 1.57 | 8.35 | / | / |
Product Output | | | | | | | |
Steel | Million Tons | 34.12 | 40.41 | 34.74 | 44.46 | 38.74 | 37.53 |
Coal | Million Tons | 5.98 | 5.73 | 5.60 | 5.47 | 5.07 | 4.80 |
Cement | Million Tons | 3.92 | 9.05 | 3.42 | 10.88 | 5.64 | 10.90 |
Power Generation | Billion kWh | 16.6 | 15.6 | 16.8 | 16.9 | 19.7 | 16.4 |
Electricity Consumption | | | | | | | |
Industry | Billion kWh | 16.72 | 18.08 | 15.18 | 18.18 | 17.26 | 18.29 |
Construction | Billion kWh | 0.19 | 0.13 | 0.16 | 0.14 | 0.23 | 0.17 |
Urban and Rural Life | Billion kWh | 1.36 | 1.08 | 1.52 | 1.14 | 1.75 | 1.21 |
Operating Income of Transportation, Warehousing and Postal Industry | Billion Yuan | 8.52 | 13.22 | 7.77 | 15.18 | 8.51 | 14.15 |
The GDP of the secondary industry in the first and second quarters of 2020 was significantly lower than that in 2019 and 2021, and the total profit of manufacturing in 2020 was about 50% of that of 2019, which indicates that COVID-19 had a certain impact on industrial production in 2020. However, the product output related to people's livelihood, such as steel, coal, cement, and power generation, was not significantly different from that reported in 2019–a year not affected by COVID-19–which demonstrates that the drop in air pollutant concentrations was likely due to the closure of industries not related to livelihood (Hu et al. 2021). The secondary industry includes industry and construction. The electricity consumption of industry dropped significantly in the first quarter of 2020 and that of the construction, which is a major source of particulate matter, showed the same trend. As work and production resumed in late February, the levels of power usage gradually returned to normal. The increase in electricity consumption by urban and rural residents and the apparent decline in the level of transportation in the first quarter of 2020 may be due to the lockdowns and travel controls imposed in this year, which led to a reduction of NO2, SO2, and CO concentrations. In summary, COVID-19 control measures had a significant impact on air pollutant emissions, which in turn affected pollutant concentrations (Rocha et al. 2022).
Analysis of factors influencing air quality during the Level I response
Interannual differences of air quality
Existing research shows that the COVID-19 outbreak has improved China’s air quality in the short term (Wang and Su 2020), and the impact of related measures on pollutant concentrations was gradually weakened as the prevention and control levels were lowered (Huang et al. 2021b). In this section, we attempted to conduct an in-depth comparative analysis of AQI and pollutant concentrations during the Level I response to examine the impact of COVID-19 control measures on air quality. The average concentrations of PM2.5, PM10, SO2, NO2, CO, and O3-8h from Jan. 24 to Apr. 30, 2020, were 56.22 µg/m3, 96.93 µg/m3, 19.38 µg/ m3, 41.73 µg/m3, 1.32 mg/m3, and 87.80 µg/m3 (Table S1), respectively. Compared with the same period in 2017–2019, the concentrations decreased by 21.26%, 24.58%, 46.02%, 26.08%, 23.70%, and 0.54%, respectively; of these, the variations of PM10, NO2, and SO2 were significantly different (P < 0.05). In the same period in 2021, the concentrations of PM2.5, PM10, and NO2 increased by 19.57%, 38.86% and 19.39%, respectively, and the differences were significant (P < 0.05).
Analysis of significant difference based on DID model
The DID model was used to study the impact of measures on air quality during Level I response. The model used the data from 23 days before and after the Level I response. The descriptive statistical results of the main variables and the estimation results of model were reported in Table S2B and Table S3B, respectively. Similar to the study results of long-term, the coefficient \(\alpha\) was found statistically significantly negative at the 1% level, which indicates that air quality improved significantly during the Level I response.
The parallel trends test was conducted, and Figure S1B presented the estimated coefficients and their 95% confidence intervals. In 2017–2019, the coefficient \(\alpha\) was significantly positive, indicating that air quality after January 24 was significantly worse than before. In 2020, the coefficient \(\alpha\) changed from significantly positive to significantly negative, which means that 23 days after the Level I response, the air quality suffered an impact. Obviously, this result further supported the use of the DID model.
Additionally, the results of the DID model illustrated that the air quality improved was mainly attribute to anthropogenic factors during the Level I response, while natural factors might have played a counterproductive role. Current environmental protection measures include basic environmental protection measures and indirect COVID-19 control measures, and good air pollution control strategies can provide crucial impacts in reducing air pollution (Othman and Latif 2021).
Impact of COVID-19 control measures based on MLR model
Assuming that basic environmental protection measures remained unchanged from 2017 to 2021, regression models were established to identify the relationship between meteorological parameters and AQI, and to quantify the impact of COVID-19 control measures on air quality. After the comparative analysis of the MLR model and the principal component regression (PCR) model, the former, which included the meteorological parameters as the independent variable and logAQI as the dependent variable, was finally adopted for subsequent simulations (Text S1). In studies simulating air quality using meteorological data, the difference between the simulated and observed values can be considered as the influence of COVID-19 control measures. Data from Feb. 1 to Apr. 30, 2017–2019 were used to build MLR models simulating air quality from Feb. 1 to Apr. 30, 2019–2021.
The model simulation results for February, March, and April 2019–2021 are shown in Fig. 5. In 2019 (the year without COVID-19), the simulated monthly AQI average values for the three months were 88.52 µg/m3, 100.23 µg/m3, and 99.66 µg/m3, respectively, and compared to them, the observed values increased by 12.54% for February, and decreased by 8.83% and 11.73% for March and April, respectively. In 2020, the observed monthly AQI average values in February, March, and April were 96.79 µg/m3, 73.19 µg/m3, and 80.50 µg/m3, respectively, while the simulated values were 136.45 µg/m3, 106.74 µg/m3, and 100.68 µg/m3, respectively. Therefore, the reduction in AQI caused by COVID-19 control measures in February, March, and April during the Level I response in 2020 was estimated at approximately 29.07%, 31.43%, and 20.04%, respectively. While the reduction in February, March, and April during the Level III response in 2021 was estimated at approximately 20.68%, 18.37% and 22.04%, respectively.
The above results showed that COVID-19 control measures played a role in improving air quality. The impact of anthropogenic factors on the AQI in February and March was significantly higher in 2020 than in 2021, which might be attributable to the difference in the Level I and Level III response. This result also shown that, to a certain extent, air pollution could be reduced by strengthening control efforts. However, the AQI in April was slightly lower in 2020 than in 2021, which may be due to the fact that the difference in COVID-19 control measures during the month of April in these two years was not significant enough to influence the air quality (Li et al. 2021a). However, it cannot be ignored that basic environmental protection measures have been continuously upgraded, such as prevention and regulations implemented by environmental protection departments, and the adjustment of the industrial structure and energy sources. Therefore, the effect of COVID-19 control measures on the AQI reduction should be smaller than that estimated by the simulations for each month.
The AQI can only reflect changes in the concentration of primary pollutants. From February to April 2020, the number of days in which PM2.5 or PM10 was the primary pollutant accounted for 70%, indicating that the AQI decrease was largely related to the decrease in particulate matter (PM) concentrations. In other words, restrictive measures implemented to combat the COVID-19 pandemic had a positive impact on reducing PM concentrations. The changes in the PM concentrations were mainly attributable to reduction in anthropogenic emissions as a result of COVID-19 control measures, although seasonal influences might also have contributed in part (Hu et al. 2015).
Analysis of factors influencing air quality during the Spring Festival
Interannual differences of air quality
Due to the large population movement and fireworks being set off during the Spring Festival period, short-term air pollution episodes are usually inevitable (Pang et al. 2021). The survey results showed that the level of anthropogenic activity during this period in 2020 dropped significantly; the average daily indices of immigrant population and emigrant population decreased by 0.629 and 0.982, respectively, compared with the values observed during the Spring Festival of 2019 (https://qianxi.baidu.com/); and the level of congestion on highways decreased significantly compared with the same period in previous years (https://jiaotong.baidu.com/).
However, during the Spring Festival in 2020, the average concentrations of PM2.5, PM10, SO2, NO2, CO, and O3-8h were 122.43 µg/m3, 158.29 µg/m3, 36.57 µg/m3, 50.29 µg/m3, 2.37 mg/m3, and 69.14 µg/m3, which were significantly higher than those recorded in 2019 and 2021. From the point of view of numerical analysis, the phenomenon that the T-test was not significant with the large difference between the means exists, which was considered to be affected by outliers (Table S1).
The hourly changes in pollutant concentrations and some meteorological elements during the Spring Festival between 2019 and 2021 may provide answers to the above questions (Fig. 6). The results showed that two heavy pollution episodes occurred on Jan. 22–23 and Jan. 25–27, 2020, in Tangshan. As anthropogenic activities were restricted during the Spring Festival of 2020, it is excluded that population movements and firework set off had an impact. Therefore, the heavy pollution episodes might have been caused by relatively unfavorable meteorological conditions (Sulaymon et al. 2021). Compared with the Spring Festival period of 2019 and 2021, the high average humidity and low average wind speed in 2020 (Table S4) promoted the formation and growth of aerosols and the formation of haze (Le et al. 2020; Wang et al. 2021a).
Impact of zone transfers
The possible sources of pollutants causing heavy pollution episodes may be local anthropogenic emissions and transportation from surrounding cities. This section simulated the backward trajectory of air masses from 0:00 on Jan. 24 to 23:00 on Jan. 30, 2020, and identified possible transport routes for pollutants through cluster analysis. Through the total spatial variance analysis, all backward trajectories could be aggregated into six categories (Fig. 7). The trajectory-3 from Qinhuangdao-Tangshan and the trajectory-5 from Beijing-Tianjin corresponded to air masses with a short-distance transport and accounted for 25.60% and 26.79% of the total trajectories, respectively. These air masses move slowly, which is not conducive to the diffusion of pollutants, and they are most likely to carry the pollution to the receiving point. In contrast, the trajectory-1 from the northeastern region and the trajectory-6 from central Inner Mongolia corresponded to air masses with a long-distance transport and accounted for 16.67% and 5.36% of the total trajectories, respectively. Due to their high speed and low possibility of occurrence, this type of air masses may have little effect on pollutant concentration in Tangshan.
Based on the analysis of the transportation path, PSCF analysis was used to determine the main potential areas representing pollution sources (Fig. 8). The results showed that local emissions and short-distance transportation were the main sources of air pollution in Tangshan. The main potential sources of PM2.5 and PM10 were located in the northern parts of Beijing and Tianjin, and this is specifically related to the increase in particulate matter concentration caused by the frequent occurrence of sandstorms in Beijing in winter and spring (Meo et al. 2021). The main potential sources of SO2 and CO were located in Tangshan and Qinhuangdao, which were the leading cities in terms of growth rates of the industrial output value in Hebei Province in 2020, and the main sources of the two above-mentioned pollutants were the industries producing steel, energy, glass, and building materials. As an important indicator to evaluate the contribution of mobile sources, the potential sources of NO2 were mainly distributed between Beijing and Tangshan, suggesting that high-intensity anthropogenic activities may occur in these areas (Aneja et al. 2001). O3 may be greatly affected by strong sources, mainly through long-distance transportation, and the potential source of this pollutant was located in northeastern China.