Abstract

The air-pollution episodes in China in January 2013 were the most hazardous in the Beijing-Tianjin-Hebei (BTH) region. PM2.5, AOD, and long-term visibility data, along with various climate and meteorological factors and the boundary-layer structure, were used to investigate the cause of the heavy-haze pollution events in January 2013. The result suggests that unfavorable diffusion conditions (weak surface winds and high humidity) and high primary-pollutant emissions have induced heavy-haze pollution in the BTH region over the past two decades. A sudden stratospheric warming (SSW), weak East Asian winter monsoon, a weak Siberian High, weak meridional circulation, southerly wind anomalies in the lower troposphere, and abnormally weak surface winds and high humidity were responsible for the severe haze pollution events, rather than an abrupt increase in emissions. Heavy/severe haze pollution is associated with orographic wind convergence zones along the Taihang and Yanshan Mountains, slight winds (1.7∼2.1 m/s), and high humidity (70%∼90%), which limits the diffusion of pollutants and facilitates the hygroscopic growth of aerosols. Recirculation and regional transport, along with the poorest diffusion conditions and favorable conditions for hygroscopic growth of aerosols and secondary transformation under the high emission, led to explosive growth and the record high hourly average concentration of PM2.5 in Beijing.

1. Introduction

Haze pollution in China has increased over the past three decades, particularly in city clusters, as a result of the rapidly developing economy, expanding anthropogenic activities and urbanization [1, 2]. The Beijing-Tianjin-Hebei (BTH) city cluster located in the North China Plain, which includes two megacities (Beijing and Tianjin) and Hebei province, has one of the largest populations (Figure 1(a)) and the third largest gross domestic product in China [35]. In recent years, large-scale regional haze pollution characterized by high concentrations of PM2.5 (particulate matter with aerodynamic diameters ≤ 2.5 μm) has frequently occurred [68], and the levels of PM2.5 continue to frequently exceed the national ambient air quality standards of China [9]; as a result, serious environmental, climatic, and health problems are rampant [10, 11].

A relatively high level of PM2.5 is usually associated with the close proximity of high precursor emissions [12]. In addition, weather/climate features play a significant role in air quality through a complex combination of processes and influences, such as emissions, transport, chemical transformations, and wet or dry removal processes. Severe pollution episodes in the urban environment mainly result from unfavorable air dispersion related to weather conditions [13, 14]. Researchers have investigated that the effects of climate change on haze/PM2.5 pollution and revealed that a significant contribution came from the weakening monsoon circulation in past decades, trapping more pollutants over Eastern China [15, 16]. On the synoptic scale, high pressure and successive low pressure systems are associated with the formation of air pollution, and pressure systems with strong gradients lead to the decrease of the air-pollution index (API) over North China [17]. Moreover, the BTH region is surrounded by mountains in the north and west directions and by Bohai Sea in the east direction (Figure 1(b)). Unique local topographic factor, together with stagnant weather with weak wind and temperature inversion, and regional transported contribution by south wind can induce heavy air pollution [1820]. Due to the combustion of coal for heating and the unfavorable dispersion attributable to geographical and meteorological conditions, haze pollution in winter over the BTH region is more serious [21, 22]. The governments, especially the Beijing municipal government, have produced a series of emission reduction measures to control air pollution [4, 23], such as changes of industrial structure and improved energy efficiency, the use of clean energy and preferred use of clean coal and gases, vehicle limitation measures, and implementation of advanced environmental standards. However, the pollution levels are still very high all the year round. Hazardous air-pollution episodes occurred in January 2013 and affected the BTH region of China the most (Figure 1(a)); thus, this event has attracted extensive attention [6, 19, 24].

Scientists have investigated the effect of regional transport and local accumulation during stagnant weather and the secondary formation of particulate on severe regional PM2.5 pollution events in January 2013 [6, 19, 25, 26]. Notably, the emission source did not increase rapidly in January 2013. Zheng et al. elaborated that [27] the changes of monthly averaged emissions in January over the Beijing-Tianjin-Hebei region were increased in rates, that is, 2.1%, 1.5%, and 2.5% for primary PM2.5, SO2, and NOx, respectively, compared with emissions in January 2012, but they were still not significant compared to the changes in pollutant concentrations and the weather/climate elements was principally responsible for these heavy pollution episodes. Zhang et al. [28] pointed out that the weakened surface winds, the anomalous southerly winds in the middle and lower troposphere, and the anomalous inversion were associated with the strong haze/fog event in Eastern China by diagnosing both its atmospheric background fields and daily evolution in January 2013. However, few studies have explored the characteristics and formation of heavy PM2.5 pollution in regard to climate factors, meteorological conditions, and boundary-layer structure over the BTH region; these factors will be important for implementing effective control measures of air pollution on the regional scale.

In this study, PM2.5, AOD, and long-term visibility data, along with various climate and meteorological factors and the boundary-layer structure, were used to investigate the cause of serious haze pollution events in January 2013 and to determine the relevant meteorological features during the most hazardous PM2.5 pollution episode (10–14 January, 2013).

2. Measurement Sites and Data Descriptions

2.1. Sites and PM2.5 Measurements

The locations of the 22 meteorological observation sites and 7 air quality sites selected for this study over the BTH region are displayed in Figure 1(b). Detailed information on the air quality sites is given in Table 1. Fourteen meteorological sites are located in the plains area, and 8 sites are located in the mountainous regions. Seven air quality sites (including 4 urban sites: Beijing Tower (BJT), Shuangqing Lu (SQL), Tianjin (TJ), and Shijiazhuang (SJZ); 1 suburban site: Xianghe (XH); 1 rural site: Yucheng (YC); and 1 regional background site: Xinglong (XL)) are local in the southern plain, where the air pollution is most severe. Air quality monitoring was conducted via an observation network established by the Institute of Atmospheric Physics (IAP), Chinese Academy of Sciences (CAS). The online measurement of the mass concentration of PM2.5 or PM10 was obtained by a Tapered Element Oscillating Microbalance (Model 1400A, R&P), and data were collected every five minutes. The instrument’s measuring principle coincides with that of the US Environmental Protection Agency (EPA) (http://www.epa.gov/ttnamti1/inorg.html). The hourly average concentrations were calculated with the 5 min data; and daily average concentrations were acquired based on 1 h data. Monthly mean deep-blue AOD at 550 nm wavelength in January 2013 over China derived from Aqua MODIS 1° × 1° level-3 daily products was shown in Figure 1(a), which revealed that the BTJ region was the most pollution area in China.

The air-pollution index (AQI) is a quantitative measure describing air-pollution levels in China based on the conversion of air-pollution data, mainly PM10, PM2.5, SO2, NO2, CO, and O3 concentrations, into a single value ranging from 0 to 500. The API is divided into six ranks representing different air quality levels with respect to their impacts on human health. The first and second air quality levels (I: AQI = 0–50 and II: AQI = 51–100) represent good air quality and no human-health risk. The third and fourth levels (III: AQI = 101–150 and IV: AQI = 151–200) indicate light pollution affecting human health to some degree. The fifth and sixth levels (V: AQI = 201–300 and VI: AQI > 300) denote heavy pollution and serious effects on human health. The AQI of the city is defined by the maximum subindex among all the six pollutants. In January 2013, the major pollutant was PM2.5; thus, air quality levels were divided by air quality index (AQI) according to PM2.5 concentrations. The corresponding concentration at the breakpoint in different air quality levels is 35, 75, 115, 150, and 250 μg/m3, respectively.

2.2. Meteorological Data

The meteorological data (daily average and monthly average) at the plains sites and mountainous sites, the visibility and relative humidity (RH) at 14:00 LST at six plains sites since 1980 (Figures 1(b) and 2), and synoptic maps were obtained from the National Meteorological Information Center and China Meteorological Administration. The visibility and RH at 14:00 LST from 1960 to 1979 were from NOAA’s National Climatic Data Center (NCDC); data was unavailable during 1965–1972. In this work, haze events were defined as days with the visibility (at 14:00 LST) <10 km and the RH < 90%, excluding days with fog, precipitation, dust storms, smoke, snow storms, and so forth. This definition is more appropriate for quantifying haze days compared with the definition using daily mean meteorological data in the BTH area [2, 21, 29]. The hourly meteorological data, wind direction, wind speed, and RH are measured with automatic weather stations (AWSs) in BJ, TJ, SJZ, XH, YC, and XL. The wind parameters are observed at a 10 m height, and the moisture information is collected at 1.5 m height. In addition, 370 AWSs with wind data and humidity information within 36°–42.5°N and 113.5°–120°E are used to obtain the gridded hourly surface-wind field and relative humidity using the Grid Analysis and Display System (GrADS) software. Gridded surface RH and wind data from the National Centers for Environmental Prediction (NCEP)/National Center for Atmospheric Research (NCAR) final global forecast system (FNL) reanalysis dataset (1° horizontal resolution) were also used to reveal the characteristic wind field and moisture distribution that have led to haze events in the BTH region. The wind directions in this paper are referred to as the directions the winds are coming from. In the Beijing urban area and Tianjin suburban area, the wind profiles with a temporal resolution of 6 min were observed with a boundary-layer wind-profile Lidar (Figure 1(b)). The radiosonde data at the Beijing meteorological station (54511) were used to analyze the vertical structure of the temperature and humidity.

Zonal circulation index for Asia (IZ), meridional circulation index for Asia (IM), East Asian winter monsoon (EAWM) intensity index, Siberian High (SH) intensity index, 1000 hPa wind anomaly, and integrated moisture flux anomaly were provided by National Climate Center (http://cmdp.ncc.cma.gov.cn/cn/index.htm). The anomalies are calculated from their departures from 1981–2010 mean values. NCEP GDAS-CPC zonal temperature anomaly downloaded from http://www.cpc.ncep.noaa.gov/products/stratosphere/strat-trop/ was used to analyze a sudden stratospheric warming (SSW).

2.3. Satellite Data

Moderate resolution imaging spectroradiometer (MODIS) onboard the Terra and Aqua satellites provides a high spatial resolution and near daily global coverage aerosol optical depth (AOD) data [30, 31]. The MODIS AOD products (MOD04/MYD04) over land comprise two datasets: dark-target AOD data and deep-blue AOD data (http://modis-atmos.gsfc.nasa.gov/MOD04_L2/index.html). The former is applicable to dark surfaces, and the latter is applicable to bright desert surfaces [3235]; these datasets have been released as MODIS aerosol products, Collection 051, with spatial resolutions of 10 km × 10 km or 1° × 1° (http://modis-atmos.gsfc.nasa.gov/MOD04_L2/atbd.html). Because less vegetation cover exists in winter over the BTH region, the deep-blue AOD data at 550 nm from MYD04_L3 C051 were used. Aqua MODIS 1° × 1° level-3 daily products were averaged over the region 35°–43° N and 112°–120° E in January from 2004 to 2013.

3. Results and Discussion

3.1. Distributions of PM2.5 and AOD in the BTH Region

Variations in the hourly average and daily average concentrations of PM2.5 at different sites and the regional average AOD in January 2013 over the BTH region are described in Figure 2. The most polluted site was SJZ, with a monthly average value of ~214 μg/m3; and the second most polluted site was YC, with a value of ~203 μg/m3. BJT, SQL, TJ, and XH were moderately polluted, with values of ~132–143 μg/m3. The least polluted site was XL, with a value of ~49 μg/m3. The YC site had a very high PM2.5 level, although it was a rural site, which indicated the influence of regional transport. Overall, the variation of PM2.5 concentrations at the different sites showed regional characteristics and periodic changes. Except XL (a regional background site), the PM2.5 concentrations far exceeded the Grade II National Ambient Air Quality Standard (75 μg/m3, GB3095-2012). The frequencies of the regional air quality at levels I and II, III and IV, and V and VI (the regional air quality level was based on the level occurring at a minimum of 4 sites) were 3 days, 13 days, and 15 days, respectively, and the corresponding regional mean PM2.5 concentrations were μg/m3, μg/m3, and μg/m3, respectively. The regional average AOD was also very high, with a range from 0.35 to 1.5. The regional background site (XL) had 6 days that exceeded the Grade II standard. Two consecutive and regionally severe haze episodes were significant (10–14 and 26–30) but with different characteristics, that is, “explosive growth” and “sustained growth” of PM2.5, respectively. However, the highest regional mean concentration (PM2.5 of ~250 μg/m3 and an AOD of ~1.3) was observed during 10–14; meanwhile, the highest hourly concentration (881 μg/m3) was found at the SQL site in Beijing.

3.2. Effects of Climatic Factors on Haze Pollution

Several datasets were used to analyze the effects of climatic factors on haze pollution over the BTH region (Figure 3). Interannual variations were calculated for precipitation, wind speed, temperature, relative humidity, and sunshine hours’ anomalies in the plains and mountainous sites. Temperature, depression of the dew point, and wind speed at 850 hPa level anomalies at BJ and XT were analyzed. IZ and IM for Asia were also displayed. Additionally, the occurrences of haze at five urban sites and a suburban site, the mean visibility for all days and haze days, and the mean AOD over the BTH region were computed for January. The precipitation presented an irregular fluctuation. The change in the averaged surface wind speed and sunshine hours over all the stations, particularly over the plains stations, showed a rapid decrease from 1960 to around 1990, whereas the wind speed remained constant in the plains sites and slightly increased in the mountainous sites; the sunshine hours exhibited a sustained decline in the plains sites and a steady oscillation in the mountain sites from 1990 to 2013; this trend is consistent with that found in many areas in China [36, 37]. Therefore, increasingly poor dispersion conditions superimposed with high aerosol loadings contributed to regional dimming. The temperature and RH increased over time as a result of industrial activities and urbanization. High moisture facilitates the hygroscopic growth of aerosols, which increases aerosol mass concentrations and extinction; these conditions lead to low-visibility events. However, the trends in the temperature, moisture, and wind speed at the 850 hPa level were not significant. IZ and IM for Asia changed dramatically over different years, and the amplitude of the variation increased. A weak negative correlation existed between the mean visibility and IZ, but a significantly positive correlation existed between the mean visibility and IM; thus good dispersion conditions resulted in high visibility when IM was high. The meridional circulation was strong, which facilitates ventilation due to the frequent exchange of southern and northern air masses.

On average, low-visibility and haze days at the urban sites were less frequent before the early 1970s; then the frequency increased and remained a constant high value before 1998 due to the rapid economic development starting in the late 1970s. Afterward, a sudden decrease in the frequency occurred, while the amplitude of the variation changed drastically. At the suburban site, an increase in haze days occurred from the early 2000s (when this region became a special economic zone) to the present, indicating that the timing of economic development plays an important role in haze pollution. Surprisingly, the interannual variation in the mean visibility during the haze days showed a significant decrease, which indicated that the primary type of haze changed from slight haze (5–10 km) to severe or serious haze (<5 km) since the beginning of the 1990s. The results suggest that unfavorable diffusion conditions (weak surface winds and high surface humidity) combined with high primary-pollutant emissions [2] have induced very heavy-haze pollution in the BTH region over the past two decades.

Under the influences of abnormally weak wind, high humidity near the ground, high upper-air temperature, low surface temperatures, weak meridional circulations, haze pollution levels, and AOD loadings were very strong in the years 2006 and 2013, whereas the air quality in 2011 was rather good due to the opposite meteorological conditions. In addition, a sudden stratospheric warming (SSW) (Figure 4) kept the stratosphere anomalously warm in 2006 and 2013 compared with the zonal temperature anomaly distribution in 2011. This phenomenon can contribute to the reversed direction of the westerly winds of the polar vortex; these conditions weaken the polar vortex. In addition, the increasing temperature in high-latitude regions may influence the circulation and contribute to a weaker pressure gradient force. As shown in Figure 5, during the course of its rotation and collapse, a weakened polar vortex in northern Asia was associated with a weak EAWM intensity index and SH intensity index; the monsoon occasionally transported weak cold air to the eastern regions of China. When cold air transport was weak, Eastern China was under the influence of southerly wind anomalies in the lower troposphere (Figure 6(a)) and abnormal water vapor divergence (less precipitation) (Figure 6(b)), which resulted in reduced wind velocity and increased air pollution.

3.3. Influence of Meteorological Conditions on PM2.5 Concentrations in January 2013

The concentration level of PM2.5 is dependent on wind and relative humidity; wind is related to dispersion and the transport of air pollution, while RH is associated with the hygroscopicity and scattering of particles. The distributions of the average surface RH and wind field for different periods in January 2013 are depicted in Figure 7. Wind speeds in the plains area of the BTH region are weaker over North China, and the lowest values are found in the cities near the Yanshan and Taihang Mountains because of the orography and surface roughness in large cities (densely distributed high buildings). A weak southerly wind prevailed over the piedmont plain in the BTH region in January 2013. Moreover, an orographic wind convergence zone was set up along the plain-mountain transition area along the Taihang Mountains, where a concentrated contamination zone was present. Therefore, the highest concentration in SJZ was partly attributed to the convergence zone. In addition, the wind convergence zones were located in the north along the Yanshan Mountains during the two heavy-haze pollution episodes. These two orographic wind convergence zones resulted in the pollution accumulation in the piedmont plain and restrained the diffusion of pollutions; as a result, severe regional haze pollution developed. In addition, haze events are strongly influenced by the ambient RH. The RH values over the plains region of the BTH area were 65%~80% in January, 70%~90% on 10–14 January, and 80%~95% on 26–30 January. Thus, an increasing number of haze events were possibly caused by the hygroscopic growth of aerosols. Because the decreased wind speed and weakened southerly winds resulted in more stable atmospheric conditions and weaker dispersion, more effort should be paid to control emissions and prevent haze events.

Meteorological factors in January 2013 for different air quality levels over the BTH region are summarized in Table 2. Compared with good air quality levels (I and II), the characteristics of the surface meteorological variables during moderate and severe haze pollution episodes in the BTH region were as follows: weaker pressure, higher temperature, particularly in the plains region, higher relative humidity, weak wind (less than 1.7~2.1 m/s), fewer sunshine hours, and lower visibility. Wind roses, the relationships between the hourly average PM2.5 concentration and wind speed, and wind direction in January at different sites are presented in Figure 8. The lowest wind speed was found in SJZ, which is consistent with the result in Figure 7. The wind speeds associated with southeast to southwest winds were very low at most of the sites, except at YC. The northwest winds were strong in BJ and XL but weak at the remaining sites. Southwest wind, northeast wind, and calm winds were dominant. Higher PM2.5 values were associated with weak winds (less than 3 m/s), except at XL and YC sites. At XL site, higher PM2.5 values were associated with strong southwest, south, and north winds; and the PM2.5 concentration at the YC site increased with the elevated wind speed, except for the northeast and southeast winds, which indicated the influence of the regional transport as the surrounding regions have high emissions.

3.4. Influence of Weather Systems, Horizontal Wind Field, and Structures of the Atmospheric Boundary Layer (ABL) on Consecutive and Regionally Severe Haze Episode

The heaviest regional haze pollution episode (January 9 to 14) is described in Figure 9. In contrast to other winter regional haze events reported by Zhao et al. [21]; the most significant characteristic of this episode was the “explosive growth” of PM2.5 and PM10 [6] on 12 January in Beijing. Synoptic maps and wind fields (Figure 10), vertical profiles of wind (Figure 11), and temperature/humidity (Figure 12) were used to explore the formation of this heavy pollution from a meteorological perspective. On 10 January, the BTH region was located in the preceding part of low pressure and experienced southerly winds; the weak mixing and reduced dispersion under these conditions led to high levels of aerosol concentrations over the entire region. A north wind arrived north of Beijing due to moving low pressure at 14:00 LST on 11 January, and the particles in Beijing were quickly dispersed. However, the north wind only affected the air pollution in Beijing, and the south winds were still prevalent in the remainder of the region. Shortly afterward, a south wind occupied the entire region, and the PM2.5 concentration in Beijing immediately rebounded due to the recirculation and transport of high-concentration pollutants from the south. Then, because of the combined effect of fast-moving weak low pressure in the north and steadily building high pressure in the south (Figure 10(c)), the winds in Beijing suddenly veered to the north and dispersed the pollutants at 02:00 LST on the 12th; the winds were soon southerly again at 14:00 LST. The pollutants were initially transported from Beijing but later returned. At the same time, the high-concentration pollutants in the south were transported into Beijing. Because the north wind only arrived in Beijing, the PM2.5 concentrations in the remaining areas increased or stayed constant under the prevailing southerly winds. This recirculation and repeated accumulation and transport from the south over a short period were responsible for the explosive growth of the PM2.5 and PM10. A few hours later, east winds prevailed in the region, and the PM2.5 concentrations at most of the sites began to decrease. However, Beijing was the convergence point of the weak southwest wind, southeast wind, and north wind (Figure 10(j)), mainly because of the urban heating effect under the weak pressure field. In addition, because the southeasterly winds transported warm humid air from the Bohai Sea, high humidity facilitated the hygroscopic growth of aerosols; as a result, severe low-visibility events ensued.

The vertical wind profile (Figure 11) revealed very weak horizontal and vertical wind speeds from 14:00 LST on the 12th to 02:00 LST on the 13th in the Beijing urban region compared with those in the Tianjin suburban region. Moreover, two layers above and below 500 m were distinguishable; therefore, the mixed layer heights were approximately less than 500 m. Figure 12 shows the thickness of the inversion layer (600–800 m), the temperature range (5–7°C), and the high humidity forcing below 500 m on 12 January. Thus, recirculation and regional transport, along with the poorest diffusion conditions and favorable secondary transformation conditions under high emissions and the hygroscopic growth of aerosols, led to the explosive growth and the record high hourly average concentration of PM2.5 in Beijing.

4. Conclusions

Haze pollution in winter over the Beijing-Tianjin-Hebei (BTH) region has become worse in recent decades. The most serious events occurred in January 2013, with a monthly regional average PM2.5 value of ~145 μg/m3 and a visibility of ~8.3 km, which were most likely due to unfavorable meteorology, rather than an abrupt increase in emissions. The interannual variation in the mean visibility during haze days showed a significant decrease, which indicated that the primary type of haze has changed from slight haze (5–10 km) to severe or serious haze (<5 km) since the beginning of the 1990s due to unfavorable diffusion conditions (weak surface winds and high humidity) combined with high primary-pollutant emissions. The characteristics of the surface meteorological variables during moderate and severe haze pollution episodes in the BTH region are as follows: weaker pressure, higher temperature, particularly in the plains region, higher relative humidity, weak winds (1.7~2.1 m/s), fewer sunshine hours, and lower visibility. Orographic wind convergence zones resulted in the pollution accumulation in the piedmont plain and restrained the diffusion of pollutions; as a result, severe regional haze pollution developed. Recirculation and regional transport, along with the poorest diffusion conditions and favorable secondary transformation conditions under high emissions and the hygroscopic growth of aerosols, led to the explosive growth and the highest hourly average concentration of PM2.5 in Beijing. Considering that decreasing wind speeds and weakened southerly winds resulted in more stable atmospheric conditions and weaker dispersion abilities, an effort should be made to control emissions and prevent haze events.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This work was supported by National Natural Science Foundation of China (nos. 41230642 and 41175107) and the CAS Strategic Priority Research Program Grant (no. XDB05020203). The authors thank Hu Bo, Guiqian Tang, Yue Wang, Yang Yang, Na Chao, and Yangchun Yu for working for aerosol measures. They acknowledge National Meteorological Information Center, NOAA’s National Climatic Data Center (NCDC), National Climate Center, and MODIS team for the data used in their work.