Introduction

Haze is a weather phenomenon in which dust, smog, and other dry particles affect the atmospheric visibility (Cai et al. 2017; Chang et al. 2009). The USA and European countries began to experience haze in the 1950s, but after decades of governance, air quality has been greatly improved. Unlike the haze in London in 1952 caused by coal burning and Los Angeles smog caused by vehicle emissions in the 1950s, haze in China is the consequence of diverse and high primary emissions, and efficient secondary generations (An et al. 2019). Although the Chinese government has made many years of efforts to control haze by transferring heavy industries from core cities to neighboring cities recently, the North China Plain has a high incidence of heavy haze in winter, and air pollutions there are more serious during the Spring Festival (Xiao et al. 2015; An et al. 2019). The Spring Festival is the most important and lively traditional festival in China, marking the beginning of the Lunar New Year by setting off fireworks and firecrackers. Previous studies have shown that residential heating, industrial production, and vehicle emissions caused severe haze pollution in North China Plain in winter (Li et al. 2014; Zhang et al. 2010; Kong et al. 2015).

Meanwhile, the topography and emission characteristics in the North China Plain (Zhang et al. 2019) are conducive to the regional transport and accumulation of pollutants because its north and west sides are Yanshan and Taihang Mountains, respectively. For example, in the winter of 2013, a nationwide haze broke out in China, especially in the North China Plain. Due to the fact that North China Plain is the largest industrial urban agglomeration in China, North China Plain bear the brunt of nationwide haze in 2013. The peak mass concentration of PM2.5 (aerodynamic equivalent diameter less than or equal to 2.5 μm) in Beijing in January 2013 reached 855 μg m−3 (Jin et al. 2020; Sun et al. 2017; Yang et al. 2015; Zheng et al. 2015). In addition to the impacts from local emissions, the average contribution rates of regional transport to Beijing's PM2.5 concentrations were found to exceed 40% during the whole year of 2013 (Ge et al. 2018). Several recent studies have shown that areas outside the North China Plain could bring pollutants into cities in the North China Plain through regional transports (Pu et al. 2020; Zhu et al. 2021).

Meteorological factors are another key factor causing air pollution in North China Plain (NCP) in winter. Stable weather conditions are not conducive to the diffusion and sedimentation of pollutants. When haze frequently occurs in winter, the North China Plain is mostly stable weather, which is generally characterized by low wind speed, low mixed layer height (MLH), and high relative humidity (Beyrich 1997; Jin et al. 2020; Sun et al. 2016; Yu et al. 2018; Kanawade et al. 2020; Zhu et al. 2018). The mixed layer is regarded as the link between the near surface and the free atmosphere and a low mixing layer height is not conducive to the diffusion of pollutants (Seidel et al. 2010).

Due to the outbreak of the COVID-19 pandemic, the Spring Festival in 2020 was very special with imposing extreme measures of nationwide lockdown by the Chinese government. The national emergency responses include extension of the Spring Festival holiday, suspension of all public transportation, closes of schools and entertainment venues, prohibition of public gatherings, enhancement of social distancing, and stay-at-home orders and closes of non-essential businesses (Wang et al. 2020a, b; Wang et al. 2021). Some recent scientific studies have shown that the improvement of air quality was one of the social benefits (Dutheil et al. 2020; He et al. 2020), while other studies have shown that limiting human activities by lockdown were not the fundamental way to alleviate air pollution (Wang et al. 2020a, b; Chang et al. 2020). It is worth studying why the unexpected haze problem still existed over the North China Plain when the air pollutant emissions were significantly reduced during the COVID-19 pandemic. The reduction of emissions from traffic and industrial sources during the COVID-19 lockdown in China provided a unique nature experiment for evaluating the efficiency of air pollution control measures (Le et al. 2020). The overall goal of this study is to investigate the regional severe haze over the North China Plain during the COVID-19 lockdown period from January 24 to January 29, 2020, by conducting a detailed analysis of the characteristics of haze at three representative cities in the North China Plain and digging the origins of the haze formation. The study explored the impact of regional transport on the formation of haze by identifying the pollution source areas of the target cities, and assessed the role of meteorological factors and ground emissions in the special atmospheric pollution during the lockdown.

Experimental

Observational data

The study period was from January 24 to January 29, 2020, of which January 24 is the New Year’s Eve of 2020 and the first day when China adopted lockdown measures. Figure 1a shows the North China Plain area (Zhang et al. 2019), which spans 5 provinces (Hebei, Shandong, Henan, Anhui, and Jiangsu) and two municipality cities (Beijing, Tianjin), and is one of the three plains and the most populous plain in China. The three receptor cities selected in this study are Beijing, Hengshui and Baoding, all located in the northern part of the North China Plain. One of the reasons for this selection is that the observational data in the three cities were available. Figure S1 shows that the North China Plain (NCP) area borders Yanshan Mountain to the north, Taihang Mountain to the west, and Bohai Sea to the east. Hourly PM2.5, NO2, O3, PM10, SO2, and CO concentrations in three representative cities (there are 12, 4 and 6 monitoring stations in Beijing, Hengshui and Baoding cities, respectively) were obtained from China National Environmental Monitoring Center (CNEMC, http://www.cnemc.cn). The hourly chemical component concentrations of PM2.5 at three representative cities were obtained from the National Air Pollution Prevention and Control Joint Center (https://view.inews.qq.com/media/6474488). The chemical components include water-soluble inorganic ions (WSIIs) (Cl, NO3, SO42−, NH4+, Na+, K+, Mg2+ and Ca2+), elemental carbon (EC) and organic matter (OM). According to the hourly concentrations of PM2.5, the pollution levels were divided into light haze (75 ≤ PM2.5 < 115 μg m−3), moderate haze (115 ≤ PM2.5 < 150 μg m−3), heavy haze (150 ≤ PM2.5 < 250 μg m−3), and severe haze (PM2.5 ≥ 250 μg m−3). The observed meteorological data were derived from the National Meteorological Data Center website (http://www.nmic.cn/site/index.html). The meteorological data with a high spatial resolution of 0.5° × 0.5° required for the Hybrid single-particle Lagrangian integral trajectory (HYSPLIT) model were obtained from the NOAA website (ftp://arlftp.arlhq.noaa.gov/pub/archives/gdas0p5).

Fig. 1
figure 1

a Map of China showing the boundary of North China Plain (NCP) and locations of Beijing (BJ), Hengshui (HS), and Baoding (BD) cities. b Hourly mass concentrations of five air pollutants; c Hourly mass concentrations of PM2.5 and chemical components in Beijing (BJ), Hengshui (HS), and Baoding (BD) cities during the study period (January 24 to January 29, 2020). The mass concentrations of chemical components were stacked in the figure. The highest values of PM10, O3, NO2, SO2 and CO were 692 μg m−3 in Baoding, 125 μg m−3 in Beijing, 90 μg m−3 in Baoding, 81 μg m−3 in Baoding, and 2.7 μg m−3 in Baoding, respectively. The maximum PM2.5 concentrations in Beijing, Hengshui and Baoding were 266 μg m−3 (13:00, January 28), 286 μg m−3 (11:00, January 26), and 641 μg m−3 (3:00, January 25), respectively. The concentrations of NO3, SO42−, NH4+, organics and K+ were the main components of fine particles (PM2.5) in Beijing and Hengshui, while organics, K+, NO3, SO42−, and NH4+ were the main components of PM2.5 in Baoding

The HYSPLIT model

We used the HYSPLIT model to calculate backward trajectories and concentration weighted trajectories (CWT). The backward trajectory provides a transport path for air pollutants, which is used to determine whether the high-concentration pollutants in one location originate from the pollution transport of another location (Draxler and Hess 1998). This study calculated the 48 h backward trajectory of PM2.5 at a height of 100 m with the starting calculation time points for each day during the study period on 00:00, 03:00, 06:00, 09:00, 12:00, 15:00, 18:00 and 21:00 (UTC time). One back trajectory for one monitoring station were calculated each time. Trajectory cluster analysis based on the trajectory space similarity is used to group a large number of trajectories. In this work, we used the angular distance to do cluster analysis which is based on the trajectory space similarity to group a large number of trajectories. This is mainly because the trajectories can be applied to locate the directions and sources from which the air masses reaching the receptor site had transported (Yan et al. 2015). The CWT method assigns a weighted concentration to each grid cell, which is obtained by averaging the sample concentrations associated with the trajectory passing through the grid cell (Reizer and Orza 2018). This method calculates the CWTij value in each grid cell (i, j) as follows (Hsu et al. 2003; Yan et al. 2015):

$$CWT_{ij} = \frac{1}{{\mathop \sum \nolimits_{l = 1}^{M} \tau_{ijl} }}\mathop \sum \limits_{l = 1}^{M} C_{l} \tau_{ijl}$$

where l and M represent the index and the total number of trajectories in the grid cell, respectively. Cl is the concentration detected at the receptor position on the arrival of trajectory l. \(\tau_{ijl}\) is the time that the trajectory l spent in the grid cell (i, j). A high CWTij value indicates that the air mass with attributions of high potential pollution has a significant contribution to the receptor site (Yan et al. 2015). In order to find the specific location of the pollution sources, the common pollution source area calculation method of Chen et al (2020) was adopted to determine the common source area of PM2.5 in the three cities during the study period.

Calculation of meteorological parameters

This study used mixed layer height (MLH), 10 m wind speed (WS), relative humidity (RH) and static weather index to represent the meteorological conditions in the three recipient cities. The MLH values can be calculated using an empirical method as follows (Nozaki 1973; Beyrich 1997):

$$MLH = \frac{121}{6}\left( {6 - P} \right)\left( {T - T_{d} } \right) + \frac{{0.169P\left( {U_{z} + 0.257} \right)}}{{12fln\left( {Z/z_{0} } \right)}}$$
$$f = 2\Omega \sin \varphi$$

where P is Pasquill's stability coefficient, T is the surface temperature, \(T_{d}\) is the dew point temperature, \(U_{z}\) is the average wind speed (m/s) at the height of Z (Z = 10 m), f is the Coriolis parameter (s−1), \(z_{0}\) is the surface roughness (0.03–0.2 m for rural area, 0.8–2 m for urban area, and 1.4 m is used in this study), \({\Omega }\) is the Earth’s angular velocity (7.29 × 10–5 rad s−1), and \(\varphi\) is the latitude value. Pasquill divides atmospheric stability into six levels from A to F (Pasquill 1961), and the six levels from A to F are replaced by 1–6 numbers in this study.

By correlating the air pollutant levels with a combination of meteorological factors in a specific region or city, Zhang et al (2017) developed the stable weather index (SWI) which included the advection, vertical diffusion and humidity and other meteorological factors (Gong et al. 2021; Liu et al. 2021). A higher SWI value means a weaker diffusion of air pollutants. The detailed information about the calculation method of SWI can be found in Zhang et al (2017). In this study, the SWT values were calculated and used to assess the atmospheric stable states of the meteorological conditions.

Modelling meteorological and chemical parmeters

The offline weather research and forecasting (WRFv3.9.1)-Community Multiscale Air Quality Modeling (CMAQv5.3.2) model was used to simulate spatiotemporal variations in meteorological and chemical parameters. Model configurations are listed in Table S1. In this study, the Carbon Bond 6 (CB6r3) scheme and AERO7 module were responsible for gas and aerosol chemistry simulations, respectively. The CMAQ version 5.3.2 (released in 2020) is the latest version and its major science advances in detail can be found in Murphy et al. (2021), Appel et al. (2021) and US EPA (2020). A brief summary relevant to the present study is presented here. Major science advances in the new aerosol module (AERO7) includes significant updates to secondary organic aerosol (SOA) chemistry, updates for chlorine, bromine, iodine chemistry, and halogen chemistries, and the addition of dimethyl sulfide (DMS) chemistry in the CB6r3 chemical mechanism (Appel et al. 2021). The CB6r3 is the most recent version of CB implemented in CMAQ53 with 4 heterogeneous reactions, 38 gas-phase reactions, and 44 gas-phase reactions and 10 heterogeneous reactions for iodine species, and 5 aqueous-phase reactions for bromine species (Appel et al. 2021).

The model domain covering China and a portion of East Asia with a horizontal resolution of 12 km × 12 km and 345 × 395 grid cells is shown in Fig. S5. The WRFv3.9.1 (Skamarock et al. 2008) model was used to provide meteorological fields for chemical simulations. The configurations and components for the WRF-CMAQ model used in this study are the same as those in Yu et al. (2014) and are summarized in Table S1. The anthropogenic emissions for China were generated by the Emission Inventory of Air Benefit and Cost and Attainment Assessment System (EI-ABaCAS) developed by Tsinghua University (Zheng et al. 2018). The EI-ABaCAS includes 16 anthropogenic sectors: fertilizer application, livestock, domestic bio-fuel, combustion, fossil fuel, solvent and other use, industry combustion, open burning, power plant, cement steel, other industry process, industry solvent use, road and non-road transport. The natural sources for biogenic emissions were calculated inline using the Biogenic Emission Inventory System version 3.14 (BEISv3.14).

Results and discussion

Air pollutants and meteorological conditions during the severe haze period

Figure 1 shows the time series of the concentrations of PM10, O3, NO2, SO2, CO, PM2.5 and chemical components of PM2.5 in the three cities during the study period. The results show that the changing trends of air pollutants in the three cities were obviously different, indicating different PM2.5 formation mechanisms. The highest levels were 692 μg m−3 for PM10 in Baoding, 125 μg m−3 for O3 in Beijing, 90 μg m−3 for NO2 in Baoding, 81 μg m−3 for SO2 in Baoding, and 2.7 mg m−3 for CO in Baoding. The maximum PM2.5 concentrations in Beijing, Hengshui and Baoding were 266 μg m−3 at 13:00 January 28 in Beijing, 286 μg m−3 at 11:00, January 26 in Hengshui, and 641 μg m−3 at 3:00, January 25 in Baoding. Baoding experienced a one-day severe pollution process from January 24 to 25, during which PM10, NO2 and SO2 all reached their peak values. The K+, Mg2+, Cl and Ca2+ ions from fireworks as indicated in Fig. 1 increased sharply after January 24, and reached their peaks on January 25, being consistent with the fact that the air pollution emissions from fireworks and firecrackers flooded into the air on the New Year’s Eve. This indicates that the fireworks and firecrackers still affected the formation of haze in the North China Plain during the lockdown period. In addition to being derived from fireworks and firecrackers, K+ was also a tracer of biomass combustion (Cheng et al. 2013). The increases of K+ also show that the suburban biomass burning contributed significantly to the K+ in the three cities.

In Beijing, after experiencing two peaks and valleys from January 25 to January 26, PM2.5 rebounded after 12 noon on January 26 and then remained at high concentrations above 147 μg m−3. This shows that after January 26, when local emissions were almost unchanged, regional transport made a certain contribution to Beijing's PM2.5. The main component of PM2.5 in Beijing was water soluble inorganic ions (WSIIs) with an average proportion of 53.6% for PM2.5. The average concentration of NO3 in Beijing was 32.0 μg m−3, accounting for 12.1–29.7% of PM2.5 with an average proportion of 20.9%, while SO42− in Beijing ranged from 10 μg m−3 to 27.8 μg m−3, accounting for 12.2% of PM2.5 on average. The average relative humidity in Beijing was 60.7 ± 16.4% (Fig. S2) during the study period, which promoted the formation of a large amount of SO42− from the heterogeneous chemical reactions and liquid phase oxidations of SO2. The concentrations of NH4+ in in Beijing ranged from 8.9 to 21.9 μg m−3, with an average proportion of 10.4% for PM2.5. Since emissions from motor vehicles and industries have declined to a certain extent under the COVID-19 lockdown (Le et al. 2020), residential heating and non-stop industrial productions in Beijing emitted a certain amount of NO3, SO42− and NH4+ precursors.

In Hengshui, there were four peak values of PM2.5 concentrations in the early stage with decreasing trends in the later stage. The average percentage of WSIIs in PM2.5 in Hengshui was 71.6%. Similar to Beijing, the main components of PM2.5 in Hengshui were NO3, SO42− and NH4+, with average mass concentrations of 36.1 ± 20.4 μg m−3, 22.3 ± 13 μg m−3, and 16.9 ± 7.6 μg m−3, respectively. The proportion of WSIIs in PM2.5 changed over time. For example, on January 25, the main component of PM2.5 was SO42− and on January 26, the main component of PM2.5 turned to be NO3, which peaked at night on January 26. Considering that the local emissions were relatively unchanged, this indicates that the contributions of regional transports to PM2.5 in Hengshui were considerable, and the formation mechanism of PM2.5 in Hengshui has changed. Figure S2 shows that the changes of relative humidity in Hengshui on January 25 were consistent with those of PM2.5. It should be noted that there were diurnal differences in the concentrations of NO3 in Hengshui with the average concentrations of 27.2 μg m−3 (daytime) and 40.7 μg m−3 (night). The higher night NO3 concentrations in Hengshui might be affected by external sources and stable weather conditions.

In Baoding, unlike Beijing and Hengshui, organic matter dominated PM2.5, accounting for 25 ± 6.4% on average. The average concentrations of chemical components of PM2.5 in Baoding in the descending order were 60.3 ± 19 μg m−3 for organic matter (OM), 35.2 ± 20.1 μg m−3 for K+, 33.9 ± 20.4 μg m−3 for SO42−, 26.6 ± 7.0 μg m−3 for NO3, 25.5 ± 6.8 μg m−3 for NH4+, 16.6 ± 11.7 μg m−3 for Cl, 10.7 ± 3.3 μg m−3 for elemental carbon (EC), 8.3 ± 3.4 μg m−3 for Ca2+, 3.9 ± 1.4 μg m−3 for Mg2+, and 2.9 ± 0.5 μg m−3 for Na+. PM2.5 in Baoding peaked on January 25 with sharp drops in the next 24 h by 380 μg m−3 because of the diffusion caused by the increases of the mixed layer height (MLH, Fig. S2). The concentration of organic matter (OM) in Baoding was extremely high, indicating that organic gas precursors played a key role in the formation of PM2.5 in Baoding. Note that the OM concentrations were still high after Jan 25 in Baoding because of the effects of the regional transports. The very high concentrations of K+ and Cl in Boading’s PM2.5 indicate substantial contributions from the emissions of fireworks and firecrackers. This is consistent with the findings of Tang et al (2016), who found that the K+, Cl, SO42−, NH4+, and NO3 concentrations in Tianjin’s PM2.5 during the firework were 65.8, 42.9, 38.9, 8.7, and 1.9 μg m−3, respectively, and that the K+ and Cl concentrations during the peak period were 23.5 and 15.5 times higher than those during Chinese Lunar New Year, respectively. Since fireworks contained charcoal and organic materials used as adhesives, they were responsible for elevated OC and EC concentrations during the firework and firecracker events (Tsai et al. 2012; Kong et al. 2015). Since potassium compounds, e.g., KNO3, KClO3, KClO4, K2CrO4, and K2Cr2O7, in black powder are used as the main oxidizers, the corresponding chemical reactions, e.g., 2KClO3 = 2KCl + 3O2 and KClO4 = KCl + 2O2, can lead to high concentrations of K+ and Cl during the firework and firecracker events (Azhagurajan et al. 2011; Tang et al. 2016).

Figure S2 shows the time series of hourly meteorological elements in the three cities. The meteorological conditions in the three cities during the study period can be summarized as low MLH and wind speed (WS) but high stable weather index (SWI) and relative humidity. For example, the MLH mean values in Beijing, Hengshui, and Baoding were 448.2 ± 294.1, 497.3 ± 209.1, and 330.0 ± 232.2 m, respectively. Comparing to the annual average MLH values in Beijing and Shijiazhuang being 594 ± 183 m and 464 ± 183 m, respectively (Zhu et al. 2018), the MLH values in the three cities during the study period can be considered as relatively low. The average WS in Beijing, Hengshui, and Baoding were 1.4 ± 0.8, 0.8 ± 0.5, and 0.8 ± 0.8 m s−1, respectively. The low wind speeds in these three cities limit the diffusion of local sources. On the other hand, the maximum relative humidity in Beijing, Hengshui and Baoding were 83%, 92% and 97%, respectively. The SWI ranges in Beijing, Hengshui and Baoding were 8–19, 11–20 and 9–21, respectively. The high SWIs in the three cities indicate that the meteorological conditions were stable during the study period. In conclusion, the stable meteorological conditions in the three cities exacerbated the formation of the haze.

Regional transport impacts based on the HYSPLIT analysis

Table 1 summarizes the average concentrations of PM2.5 and its chemical components, and the percentages for each cluster and three cluster trajectories for each receptor city are shown in Fig. S3. The results of all backward trajectories in Fig. S3 indicate that most of the air masses came from Inner Mongolia, Liaoning, Hebei and Tianjin.

Table 1 Number of backward trajectories and average concentrations of PM2.5 and chemical components for different clusters in Beijing, Hengshui, and Baoding during the study period

For Beijing, Cluster 1, which accounted for 54.1% and carried 51.0% of the pollution trajectories, transported from the neighboring areas to Beijing through a small area in northern Hebei province. Cluster 1 had the highest concentration of PM2.5 (177.3 ± 61.2 μg m−3), which was the most important transport path for PM2.5 in Beijing. Cluster 2 (3) accounted for 37.7% (8.2%) and carried 39.5% (9.5%) of the pollution trajectories. According to the concentrations of chemical components, Cluster 2 carries the most K+ and Cl probably generated by firework because Cluster 2 was not a cross-sea transport as shown in Fig. S3.

For Hengshui, all cluster trajectories were cross-sea transports in the northeast direction with obvious long-distance transport characteristics as indicated in Fig. S3. Cluster 1 with the longest transport path came from western Inner Mongolia and Liaoning Province, accounting for 43.8% and carrying 44.3% of the pollution trajectories (see Table 1). Clusters 2 and 3 passed through western Liaoning, southern Tianjin, and northern Shandong. The PM2.5 concentrations (157.6 μg m−3) carried by the air mass in Cluster 2 was the highest. The Na+ concentrations in Cluster 2 were the highest, further indicating that cross-sea transports had a certain influence on the formation of PM2.5 in Hengshui.

For Baoding, except for Cluster 3, the other two Clusters all passed through the western Liaoning and Inner Mongolia. Cluster 3, which accounted for 21.1% and carried PM2.5 with an average concentration of 216.7 μg m−3, moved around Baoding, indicating a great impact of local pollution sources on PM2.5 in Baoding. The percentages of Clusters 1 and 2 were 45.4% and 33.5%, respectively, and Cluster 2 carried the highest concentration of PM2.5 (295.9 μg m−3).

To further study the impact of regional transports on the formation of haze, Fig. 2 shows the 48-h backward trajectories in three cities, superimposed by the observed PM2.5 concentrations over the North China Plain every 12 h within 48 h. As can be seen, the spatial variations of PM2.5 in the three cities were consistent with those of the monitoring stations on the transport paths. The pollution incidents in Beijing were a typical full-time severe pollution event. PM2.5 concentrations in Beijing were higher at 13:00 on January 28 than those in other four time points. From January 26 to January 28, there was a rising process of PM2.5 in the western cities of Liaoning, which were on the pollution transport path. The concentrations of PM2.5 in Beijing remained at a high level above 147 μg m−3, indicating that PM2.5 in western Liaoning was imported into Beijing along the pollution transport path. Figure 2 indicates that nearby cities located in the pollution transport path had a greater impact on the increases of Beijing’s PM2.5. Figure 2 also shows that the spatial variations of PM2.5 concentrations in Hengshui and Baoding were consistent with those at the monitoring stations on the transport path. For the pollution event from January 24 to January 26, the cities in Liaoning province located on the northeast transport path for Hengshui experienced a substantial increase of PM2.5 within 24 h from 11:00 on January 24 to 11:00 on January 25. The PM2.5 concentrations in Hengshui began to rise from 23:00 on January 25, and reached a peak on January 26. Baoding’s PM2.5 experienced a sharp rise from 15:00 on January 24. At the same time, the concentrations of PM2.5 in cities of western Liaoning and Hebei Provinces on the transport path for Baoding also experienced the same rising processes.

Fig. 2
figure 2

Maps of 48 h polluted backward trajectories overlaid with the observations of PM2.5 concentrations (μg·m−3) at different times (in the format of YYYYMMDD_HH). The spatial variations of PM2.5 in the three cities were consistent with those of the monitoring stations on the transport paths. The pollution incidents in Beijing were a typical full-time severe pollution event. As can be seen, the PM2.5 concentrations of the three recipient cities partly depend on the cities located on pollution transport paths for them. The northeast transport path through Liaoning province had the greatest impact on PM2.5 in the three recipient cities. The results show that the northeast transport path dominates the regional transport for the haze formation in NCP

In summary, the PM2.5 concentrations of the three receptor cities mainly depend on the cities located in the pollution transport paths for them. The northeast transport path, especially through Liaoning province, had the substantial impact on PM2.5 in the three recipient cities. The previous studies found that the regional transports were mainly a southerly transport path for NCP (Zhu et al. 2016; Du et al. 2020; Lv et al. 2020). However, this study shows that the northeast transport path dominates the regional transport for the haze formation in the North China Plain.

Identification of common source areas for pollution

The concentrations of PM2.5 in the receptor cities were combination effects of local emissions and regional transports. Figure 3 shows the concentration weighted trajectories (CWT) for the entire and severe pollution periods and the CWT for the light, medium and heavy pollution periods are shown in Fig. S4. Figure 3 indicates that different cities had different degrees of contributions to PM2.5 in the three cities. The regions with the CWT values higher than 250 μg m−3 (see Fig. 3) included Chaoyang and Huludao in Liaoning province, Qinhuangdao, Tangshan, Langfang and Chengde in Hebei province, and Tianjin, all of which contributed high concentrations of PM2.5 to Beijing. For Hengshui, the potential sources of pollution during the severe pollution period were Panjin, Anshan and Yingkou in Liaoning province. It was worth noting that Liaoning province had a large amount of PM2.5 input for Hengshui for all the pollution periods. For Baoding, Cangzhou, Langfang, Chaoyang, Jinzhou and Tianjin were potential source areas during the heavy pollution period. Note that some of the air masses from Liaoning were transported across the Bohai Sea, which is located in the east of the North China Plain.

Fig. 3
figure 3

source areas of PM2.5 (PM2.5 ≥ 75 μg m−3) during the entire study period (24th January–29th January, 2020). The full names of the abbreviation of cities in the figure are as follows: AS-Anshan, BT-Baotou, CD-Chengde, CF-Chifeng, CY-Chaoyang, CZ-Cangzhou, HHHT-Hohhot, HLD-Huludao, JZ-Jinzhou, LF-Langfang, PJ-Panjing, QHD-Qinhuangdao, TJ-Tianjijng, TL-Tongliao, TS-Tangshan, WLCB-Wulanchabu, YK-Yingkou, ZJK-Zhangjiakou. For Hengshui, the potential sources of pollution during the severe pollution period were Panjin, Anshan and Yingkou in Liaoning province. For Baoding, Cangzhou, Langfang, Chaoyang, Jinzhou and Tianjin were potential source areas during the heavy pollution period

Concentration weighted trajectory (CWT) maps of PM2.5 for Beijing, Hengshui, Baoding during the entire and severe haze periods (PM2.5 ≥ 250 μg m−3), and common

Figure 3 shows the common source areas of PM2.5 in the three cities during the study period. Controlling the common source area is an effective choice for air pollution control (Yu et al. 2018). As can be seen, western Liaoning, northern Hebei, and southern Tianjin were the common source areas for PM2.5 pollution events in the three recipient cities. This is consistent with the fact that Hebei province is the largest industrial cluster in China (Wang et al. 2014) and Liaoning province is a region with industry as its main industry (Shi et al. 2020). Although a temporary industrial shutdown was implemented to curb the spread of COVID-19, industrial emissions have not dropped significantly due to the fact that the main industrial sectors have not yet been completely closed (Wang et al. 2020a, b).

Regional transport impacts based on the WRF-CMAQ model simulations

To determine if the northeast transport path had the substantial impact on PM2.5 in the three receptor cities during the study period, the WRF-CMAQ model was applied to simulate this regional severe haze episode in the retrospective mode. The concentration weighted trajectories (CWT) values for PM2.5 were calculated on the basis of the air mass back trajectories and their associated PM2.5 concentrations observed at the three receptor sites to locate the source areas with the largest potential contributions to high concentrations. Figure S6 shows the spatial distributions of the CWT values created by combining all results at the three cities during the study period. Followed Yu et al. (2016), an emission reduction scenario was created on the basis of the CWT values in Fig. S6 by removing all anthropogenic emissions for grid cells with CWT ≥ 75 μg m−3, while those will remain uncontrolled for grid cells with CWT < 75 μg m−3. Figures S7 and S8 show the comparisons of the spatial distributions of anthropogenic emissions for primary PM2.5, NH3, NOX, SO2, VOCS, and POC between the base and reduction cases. Figure S9 shows the daily mean comparisons of the simulation results for the spatial distributions of PM2.5 concentrations for the base and reduction cases and their differences for each day from January 24 to 29, 2020. The maps of the PM2.5 reduction amounts between the base and reduction cases overlaid with the wind fields for different times in Fig. S10 clearly show that the prevailing northeasterly winds during study period brought air pollution continuously from the northeast areas to the North China Plain regions, causing the severe haze formation there. Comparisons of the time series of PM2.5 and its chemical components concentrations in the three receptor cities for the base and reduction cases in Fig. S11 also show that the peak values for all these pollutants in the reduction cases in the three receptor cities were effectively decreased by more than 60%. As summarized in Tables S2, S3 and S4, the PM2.5, NH4+, NO3, OM, and SO42− concentrations in the reduction case decreased at most by − 106.73 μg m−3, − 8.93 μg m−3, − 25.96 μg m−3, − 13.75 μg m−3 and − 5.51 μg m−3 in Beijing, respectively, while they were − 144.38 μg m−3, − 12.01 μg m−3, − 34.15 μg m−3, − 21.77 μg m−3 and − 7.79 μg m−3 in Hengshui, respectively, and were − 148.32 μg m−3, − 12.21 μg m−3, − 35.86 μg m−3, − 23.68 μg m−3 and − 8.56 μg m−3 in Baoding, respectively. The mean PM2.5, NH4+, NO3, OM, and SO42− concentrations in the reduction case could decrease by more than 60% in the three receptor cities relative to the base case. The model sensitive simulations results highlight the fact that regional severe haze in North China Plain during the coronavirus lockdown was caused mainly by the anthropogenic emissions through the northeast transport paths.

Conclusion

This study investigated the causes of the regional severe haze over the North China Plain during the COVID-19 lockdown period on the basis of chemical observations in the three cities (Beijing, Baoding and Hengshui) and model simulations. The concentrations of NO3, SO42−, NH4+, organic matter and K+ were found to be the main components of PM2.5 in Beijing and Hengshui, while they were OM and K+, NO3, SO42−, and NH4+ in Baoding. The results indicate that the emissions from firework, firecrackers, residential heating, and non-stop industries under the stable meteorological conditions with low wind speed, low mixing layer height, and high relative humidity caused this severe haze pollution in North China Plain. The results of backward trajectories in the three cities indicate that most of the air masses influencing the cities came from Inner Mongolia, Liaoning, Hebei and Tianjin. For Hengshui, all cluster trajectories are cross-sea transports in the northeast direction with obvious long-distance transport characteristics by passing the western Inner Mongolia and Liaoning Province. For Baoding, the two cluster trajectories also passed through the western Liaoning and Inner Mongolia. It was found that the common source areas for the haze formation in the three cities were located in western Liaoning, northern Hebei, and southern Tianjin as the common source cities. The stagnant airflow and uninterrupted emissions from non-stop industries (e.g., power plants and petrochemical facilities) can still cause regional severe haze formation with high humidity to promote aerosol heterogeneous chemistry. The results of the cluster trajectories and model simulations show that different from the previous studies that mainly found a southerly transport path of air pollution for North China Plain, the northeast transport path dominates the regional transport for this haze formation. This indicates that the transport paths can be variable for a city for different haze events in regional scale. This implies that the variabilities of transport paths should be taken into account in the regulation of air pollution in the North China Plain, especially in the case of short-term severe pollutions.