Characteristics of Boundary Layer Structure during a Persistent Haze Event in the Central Liaoning City Cluster, Northeast China
  J. Meteor. Res.  2018, Vol. 32 Issue (2): 302-312   PDF    
http://dx.doi.org/10.1007/s13351-018-7053-6
The Chinese Meteorological Society
0

Article Information

Li, X. L., Y. F. Wang, L. D. Shen, et al., 2018.
Characteristics of Boundary Layer Structure during a Persistent Haze Event in the Central Liaoning City Cluster, Northeast China. 2018.
J. Meteor. Res., 32(2): 302-312
http://dx.doi.org/10.1007/s13351-018-7053-6

Article History

Received April 27, 2017
in final form November 3, 2017
Characteristics of Boundary Layer Structure during a Persistent Haze Event in the Central Liaoning City Cluster, Northeast China
Xiaolan LI1, Yangfeng WANG1, Lidu SHEN2, Hongsheng ZHANG3, Hujia ZHAO1, Yunhai ZHANG1, Yanjun MA1     
1. Institute of Atmospheric Environment, China Meteorological Administration, Shenyang 110166;
2. Shenyang Regional Climatic Center, Shenyang 110166;
3. Laboratory for Climate and Ocean–Atmosphere Studies, Department of Atmospheric and Oceanic Sciences, School of Physics, Peking University, Beijing 100871
ABSTRACT: The characteristics of boundary layer structure during a persistent regional haze event over the central Liaoning city cluster of Northeast China from 16 to 21 December 2016 were investigated based on the measurements of particulate matter (PM) concentration and the meteorological data within the atmospheric boundary layer (ABL). During the observational period, the maximum hourly mean PM2.5 and PM10 concentrations in Shenyang, Anshan, Fushun, and Benxi ranged from 276 to 355 μg m–3 and from 378 to 442 μg m–3, respectively, and the lowest hourly mean atmospheric visibility (VIS) in different cities ranged from 0.14 to 0.64 km. The central Liaoning city cluster was located in the front of a slowly moving high pressure and was mainly controlled by southerly winds. Wind speed (WS) within the ABL (< 2 km) decreased significantly and WS at 10-m height mostly remained below 2 m s–1 during the hazy episodes, which was favorable for the accumulation of air pollutants. A potential temperature inversion layer existed throughout the entire ABL during the earlier hazy episode from 0500 Local Time (LT) 18 December to 1100 LT 19 December, and then a potential temperature inversion layer developed with the bottom gradually decreased from 900 m to 300 m. Such a stable atmospheric stratification further weakened pollutant dispersion. The atmospheric boundary layer height (ABLH) estimated based on potential temperature profiles was mostly lower than 400 m and varied oppositely with PM2.5 in Shenyang. In summary, weak winds due to calm synoptic conditions, strong thermal inversion layer, and shallow atmospheric boundary layer contributed to the formation and development of this haze event. The backward trajectory analysis revealed the sources of air masses and explained the different characteristics of the haze episodes in the four cities.
Key words: haze event     thermal inversion layer     atmospheric boundary layer     Northeast China    
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, expansion of anthropogenic activities and urbanization (Shao et al., 2006; Fu et al., 2014; Wang et al., 2014a). In recent years, heavy and persistent haze events have been observed in different regions of China (Ji et al., 2012; Liu et al., 2015; Peng et al., 2016), resulting in atmospheric visibility reduction, increased air pollution, and risks to public health (Chan and Yao, 2008; Chen et al., 2012; Che et al., 2014; Zhang et al., 2014; Shen et al., 2015; Tang et al., 2017).

Haze formation is closely related to meteorological conditions and high aerosol mass loadings (Wang et al., 2014b). In the northern part of China, heavy air pollution often occurs in the heating season or during periods mainly due to large pollution emissions associated with agricultural residue burning and unfavorable meteorological conditions for pollutant dispersion (Fu et al., 2014; Wang et al., 2014b; Zheng et al., 2015; Chen et al., 2017). Since the heating season begins earlier and lasts longer in Northeast China as compared to most regions in China and is combined with burning of agricultural residues, haze occurs frequently during winter in this region. For example, a heavy and persistent haze event with a maximum PM2.5 concentration exceeding 1000 μg m–3 occurred during 1–8 November 2015 in Shenyang and profoundly shocked both the public and the government. Several studies have reported that haze pollution in Northeast China is serious (Cui et al., 2015), with an increasing trend in particulate matter (PM) observed at some stations in this region from 2006 to 2014 (Wang et al., 2015c).

Many efforts have been made to identify the meteorological causes of haze events involving the effects of surface meteorological conditions, synoptic patterns, climatic background, and local atmospheric circulation induced by the topography (Chen et al., 2008; Liu et al., 2009; Zhu et al., 2012; Hu et al., 2014, 2016). Zhu et al. (2012) demonstrated that the weakened monsoon circulation of recent decades has made a significant contribution to trapping pollutants over eastern China. The modeling results reported by Miao et al. (2015, 2017) revealed that local atmospheric circulation plays an important role in surface air pollution over the Beijing–Tianjin– Hebei region. In addition, the structure of the atmospheric boundary layer (ABL) is also crucial to the formation and evolution of haze (Wang et al., 2006; Yang et al., 2010; Lu et al., 2011; Tang et al., 2016, Ye et al., 2016). Wu et al. (2009) and Sun et al. (2013) found that the vertical dynamic and thermal variations within the ABL influence pollutant concentrations. Quan et al. (2013) compared the ABL evolution during clear and hazy conditions in Tianjin and considered the possibility of a positive feedback cycle between atmospheric aerosols and lower ABL height, which would induce heavy surface-level pollution in cities. Liu et al. (2015) analyzed ABL characteristics during a typical haze process in January 2014 over the Pearl River Delta region and found that weak wind, low ABL height, and a thermal inversion layer under stable atmospheric conditions were major reasons for haze formation. Hu et al. (2014) investigated the impact of the Loess Plateau on the ABL structure and air quality in the North China Plain (NCP). Based on field measurements, similar studies have also been conducted in other megacities or regions, such as the Beijing– Tianjin–Hebei region (Liao et al., 2014; Wang et al., 2014a), Xi’an in western China (Wang et al., 2016), and eastern China (Li et al., 2015; Peng et al., 2016), but such studies are still very rare in Northeast China (Zhang et al., 2010a, b; Liu et al., 2011). These previous studies have contributed to our understanding of the effects of ABL evolution on the changes in surface air pollution and provide a reference for studies of the interaction of ABL meteorology and aerosols in haze modeling and prediction (Wang et al., 2015a, b).

In this study, the characteristics of ABL structure during a persistent, severe regional haze event over the central Liaoning city cluster during 16–21 December 2016 were analyzed, by using the vertical profiles of ABL meteorological variables and surface PM concentration, combined with horizontal wind fields retrieved from the ECMWF reanalysis data and back trajectory analyses from the Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) model of NOAA (Draxler and Rolph, 2003; Wang et al., 2010). In Section 2, a brief introduction to the study area and details of the data and methods are described. The results and discussion, given in Section 3, mainly address temporal variation in PM concentration and surface meteorological variables, and the horizontal and vertical structure of the boundary layer. Finally, the conclusions of the study are presented in Section 4.

2 Data and methods 2.1 Description of the study area

Liaoning Province is an industrial region that has made a substantial contribution to the economic development of Northeast China (Che et al., 2015). Our study focused on the central Liaoning city cluster, which comprises Shenyang, Anshan, Fushun, and Benxi (Fig. 1a). These four cities are located in a heavily populated industrial region of Northeast China (38.43°–43.26°N, 118.53°–125.46°E), a major source area for aerosol pollution originating from a variety of human activities, industrial, and transportation sources (Ma et al., 2012; Zhao et al., 2013a; Che et al., 2015). Particularly in winter, emissions from coal burning for domestic heating are a well-established source of pollution (Zhao et al., 2013b; Che et al., 2015).

Figure 1 (a) Geographical locations of four cities in the central Liaoning city cluster of Northeast China, and (b) positions of 11 environmental monitoring stations (yellow circle) and the boundary layer experiment site (red triangle) in Shenyang.
2.2 PM concentrations and surface meteorological parameters

Large-scale regional haze pollution has been characterized by high concentrations of PM2.5 in recent years (Wang et al., 2014a). In this study, the hourly mass concentrations of PM10 and PM2.5 in Shenyang, Anshan, Fushun, and Benxi during a typical haze event from 16 to 21 December 2016 were obtained from the Liaoning real-time air quality monitoring system. There are eleven monitoring stations located in Shenyang, seven in Anshan, and six each in Fushun and Benxi. The ambient real-time monitoring instruments, Models 5030i and 5014i were used to measure the concentrations of PM2.5 and PM10 at each station, and then the averaged PM data from different stations in each city were used in this study.

The simultaneous hourly mean surface meteorological parameters, i.e., wind speed (WS), wind direction (WD), air temperature (Ta), relative humidity (RH), and atmospheric visibility (VIS), from the national weather stations in the four cities were used to analyze the surface meteorological conditions.

2.3 Boundary layer structure

The ECMWF reanalysis data with a spatial resolution of 0.125° × 0.125° were used to retrieve the horizontal flow fields at 10 m above ground level (AGL) over Liaoning from 16 to 21 December 2016. The ECMWF data were obtained four times a day at 0000, 0600, 1200, and 1800 UTC, and can be freely downloaded from http://apps.ecmwf.int/datasets/data/interim-full-daily/levtype=sfc/.

A boundary layer observational campaign was conducted from 1100 Local Time (LT) 17 to 1400 LT 23 December 2016 in Baitapu Town (41.6841°N,123.4160°E) about 8 km from the southern edge of the main urban zone of Shenyang (Fig. 1b). The Model Chuangzhi Tan Kong-1 (CZTK-1) sounding system, developed by the Institute of Atmospheric Physics of the Chinese Academy of Sciences, was used to measure the vertical distributions of WS, WD, Ta, and RH, with detective resolution of 0.1 m s–1, 0.1°, 0.1°C, and 0.1%, respectively. Sounding balloons were released at an open balcony of a low building (height < 3 m) eight times per day at 0200, 0500, 0800, 1100, 1400, 1700, 2000, and 2300 LT, and a total of 50 groups of profile data were obtained. The ascending velocity of the sounding balloons was about 150 m min –1, and sounding data were recorded at an interval of 1 s. The detection duration was about 20–40 min, and the detection height usually reached up to 1500–4500 m. It should be noted that no precipitation occurred from 17 to 21 December during the haze event but snow weather lasted during 22–23 December after the haze event. Following the methods used in Liu et al. (2015), vector winds at different observational heights were first calculated based on the vertical profiles of WS and WD, and then interpolated at each 50 m AGL. The profiles of Ta and RH were averaged every 10 m after data quality control, and then potential temperature (θ) was calculated.

The atmospheric boundary layer height (ABLH) can be estimated by using the vertical potential temperature profiles from sounding data based on the theory of boundary layer evolution, with the ABLH estimated at the altitude where the potential temperature first exceeds the minimum potential temperature within the boundary layer by 1.5 K (Nielsen-Gammon et al., 2008; Hu et al., 2014).

2.4 Back trajectory analyses

The HYSPLIT model from the NOAA was used to analyze the transport path of atmospheric pollutants during the haze event. Two-day back trajectory analyses were conducted for the four cities and the information is listed in Table 1.

Table 1 Information of back trajectory simulation
City Position Target altitude Starting time
Shenyang 41.77°N, 123.50°E 200, 500, and 1000 m AGL 1500 LT 16 December 2016
Anshan 41.08°N, 123.00°E 200, 500, and 1000 m AGL 0300 LT 17 December 2016
Fushun 41.88°N, 123.95°E 200, 500, and 1000 m AGL 0600 LT 17 December 2016
Benxi 41.32°N, 123.78°E 200, 500, and 1000 m AGL 1200 LT 17 December 2016
3 Results and discussion 3.1 Overview of the haze event in the central Liaoning city cluster 3.1.1 Identification of hazy episode

A haze event is characterized by an apparent decrease of VIS less than 10 km and ambient RH smaller than 80% lasting for several hours (Fu et al., 2008; Leng et al., 2016). Under the conditions of 80% < RH < 90%, the event is usually regarded as a complex of haze–fog co-occurring or transition ( Leng et al., 2014; 2016), while under the conditions of RH ≥ 90%, it is referred to as a fog event.

Figure 2 depicts the temporal variations of hourly mean PM mass concentrations and some surface meteorological parameters in Shenyang, Anshan, Fushun, and Benxi from 16 to 21 December 2016. As shown in Figs. 2ad, VIS remained below 10 km during 17–21 December in all cities except for Anshan where VIS declined to below 10 km after the midday of 18 December. RH mostly ranged from 60% to 90% in Shenyang and Fushun and from 40% to 90% in Anshan and Benxi (Figs. 2il). Since PM2.5 concentration was quite high during the whole period, this event is classified as a haze event in the present study. During the observational period, air quality in most cities reached Level 5 (daily air quality index (AQI) > 200) according to the ambient air quality standards in China, which means the central Liaoning city cluster suffered from a persistent and severe haze pollution. Here, AQI is a numerical rating and can be calculated by using air pollutant concentrations over a specified averaging period.

3.1.2 Particulate mass concentration and atmospheric visibility

Overall, PM10 and PM2.5 concentrations varied similarly during the haze event in all cities. The maximum hourly PM2.5 in different cities ranged from 276 to 355 μg m–3 and the maximum hourly PM10 in different cities ranged from 378 to 442 μg m–3. The lowest hourly VIS reached 0.14, 0.64, 0.18, and 0.63 km in Shenyang, Anshan, Fushun, and Benxi, respectively (Figs. 2ad). Taken Shenyang as an example, the PM2.5/PM10 ratio averaged during the haze event reached 75.3% (not shown), which is consistent with results from previous literatures that low visibility is influenced remarkably by fine particles in Northeast China (Ma et al., 2011; Zhao et al., 2013a, b). During the haze event, PM2.5 concentration was larger in Anshan (158.3 ± 22.2 μg m–3) and Shenyang (157.9 ± 31.0 μg m–3), and relatively smaller in Fushun (121.3 ± 48.9 μg m–3) and Benxi (135.6 ± 60.8 μg m–3). By comparing the variations of PM2.5 and VIS in Fushun with those in Anshan, it is interesting to find that Anshan had larger PM2.5 concentration but still higher visibility during 17–18 December, which was probably due to a lower RH (40%–60%) condition in Anshan. Some previous studies have indicated that an increase in RH within a certain range can cause significant VIS degradation (Stock et al., 2011; Liu et al., 2012).

Figure 2 Temporal variations of hourly mean (a–d) PM2.5 and PM10 mass concentrations from environmental monitoring stations and atmospheric visibility, (e–h) wind speed and wind direction at 10-m height, and (i–l) relative humidity and air temperature at 2-m height from meteorologicalstations in Shenyang, Anshan, Fushun, and Benxi from 16 to 21 December 2016. The dashed lines in (a–d) are the visibility at 10 km.
3.1.3 Surface meteorological conditions

Wind speed at 10 m height was lower than 2 m s–1 for most of the hazy period but increased to higher than 4 m s–1 on 21 December in all cities except for Benxi (Figs. 2eh). Wind direction in Shenyang changed frequently but was basically southeast or southwest during the haze event (Fig. 2e). In other three cities, WD also turned into southeast or southwest when PM concentrations increased rapidly on 19 December (Figs. 2fh). Ta mostly ranging from –10 to 5°C was correlated negatively with ambient RH (Figs. 2il).

3.1.4 Synoptic patterns

The 925-hPa geopotential height (unit: m) fields at 0800 LT 18–21 December are shown in Fig. 3. The central Liaoning city cluster was located in the front of a weak and slowly moving low pressure system on 18–20 December (Figs. 3ac). The weak pressure gradients represented by very sparse isobaric lines resulted in weak large-scale air circulation and low wind speed in the near-surface layer in this region. Such stable synoptic conditions favored the accumulation of air pollutants and inhibited their horizontal dilution. With the cold high pressure moving towards the southeast, the surface isobaric lines over this region became much denser on 21 December (Fig. 3d), and strong winds were conductive to the dispersion of air pollutants.

Figure 3 The 925-hPa geopotential height (unit: m) fields at 0800 LT (a) 18, (b) 19, (c) 20, and (d) 21 December 2016 retrieved from the ECMWF reanalysis data. The “SY” denotes Shenyang.
3.2 Effects of boundary layer structure on the haze event 3.2.1 Characteristics of horizontal wind fields

Vector wind field reflects the regional flow spatial characteristics; a larger vector wind indicates better dispersion conditions for air pollutants, and vice versa (Wu et al., 2007). Figure 4 shows the average daily 10-m wind fields over Liaoning during 16–21 December 2016. Liaoning was dominated by the southwesterly and westerly flows on 16 December and since the following day the central Liaoning city cluster was mainly controlled by the southeasterly flow that was favorable for the transportation of water vapor from the Bohai Sea to this region and resulted in high ambient RH condition. The daily mean wind speed over the central part of Liaoning decreased day after day from 16 to 20 December, which was favorable for the accumulation of air pollutants in this region. With the arrival of cold and dry air masses brought by the strong northeasterly flows on 21 December, the heavy haze event ultimately ended.

Figure 4 Spatial variations of the daily mean wind field retrieved from the ECMWF reanalysis data over Liaoning Province from 16 to 21 December 2016. SY, FS, AS, and BX represent Shenyang, Fushun, Anshan, and Benxi, respectively.
3.2.2 Characteristics of vertical wind distribution in Shenyang

The vertical distribution of vector wind from 0200 LT 18 December to 0200 LT 22 December 2016, measured from the boundary layer experiment in Shenyang, is shown in Fig. 5. It can be seen that during 0200–1100 LT 18 December wind speed at all heights below 2 km decreased rapidly with time, and wind direction rotated with increasing height in a clockwise direction. Meanwhile, PM concentrations at surface began to increase sharply (Fig. 2a). Weak southerly and southwesterly flows then prevailed until 1100 LT 19 December, corresponding with PM2.5 increasing to approximately 300 μg m–3 (Fig. 2a). Wind speed was still very small at low altitudes (< 600 m), whereas winds fluctuated frequently at high altitudes, contributing to the fluctuation of PM concentrations.

Figure 5 Time–height cross-section of winds from the sounding data in Shenyang from 0200 LT 18 December to 0200 LT 22 December 2016.
3.2.3 Characteristics of vertical temperature distributions in Shenyang

The vertical distribution of air temperature based on sounding data from 0200 LT 18 December to 0200 LT 22 December 2016 in Shenyang is shown in Fig. 6. With the increasing PM concentrations during the haze event, Ta in the boundary layer decreased obviously, and the depth of the warm advection declined until 1400 LT 20 December. These stable and upper-layer warm atmospheric conditions were very favorable for the maintenance of the haze pollution event.

Figure 6 Time–height cross-section of air temperature from the sounding data in Shenyang form 0200 LT 18 December to 0200 LT 22 December 2016.

To clearly outline the atmospheric stability in ABL, Fig. 7 shows potential temperature profiles in Shenyang at 0500, 1100, 1700, and 2300 LT from 18 to 22 December 2016. From 0500 LT 18 December to 0500 LT 19 December, a potential temperature inversion layer existed from the surface to the height of 1500 m, and the inversion intensity gradually increased with time. During this period, the surface PM2.5 concentration increased from 123 to 209 μg m–3. After 0500 LT 19 December, θ varied very small at the low altitudes of ABL (< 800 m), but the enhancement of the potential temperature inversion was apparent at high altitudes of ABL. The potential temperature inversion intensity gradually increased from 8.9 K km–1 at 1700 LT 19 December to 15.8 K km–1 at 2300 LT 20 December. Additionally, the bottom of inversion height gradually decreased from approximately 900 m at 1700 LT 19 December to approximately 300 m at 2300 LT 20 December, corresponding to a very high PM concentrations on the following day. After that, the bottom of inversion height increased to 900 m again at 1100 LT 22 December, which strengthened the vertical mixing of air pollutants and favored with the dilution of air pollutants at surface. It should be noted that the warm advection near the ground at 1700 and 2300 LT 19, 21, and 22 December was mainly due to the heat island effect over the urban region.

Figure 7 Potential temperature profiles from the sounding data in Shenyang at 0200, 0800, 1400, and 2000 LT from 18 to 22 December 2016.
3.2.4 Evolution of the ABLH in Shenyang

The ABLH plays an important role in the evolution of haze episodes. Decreasing the height of ABL can normally hold the pollutants within the shallow surface layer, suppress the vertical atmospheric dilution, and finally lead to regional environment shrouded by pollution (Kim et al., 2007; Leng et al., 2016). The ABLH estimated by using sounding data was plotted in Fig. 8. The regular diurnal variation of the ABLH, i.e., highest in the middle of the day and lowest at nighttime, was not so typical during heavy pollution period. From 0200 LT 18 December to 1100 LT 19 December, the ABLH was mostly below 400 m, and PM2.5 and PM10 concentrations during this period reached 265 and 370 μg m–3. When the ABLH increased rapidly to 900 m at 1400 LT 19 December and remained larger than about 600 m until 0800 LT 20 December, PM2.5 and PM10 concentrations dropped to 142 and 180 μg m–3, respectively. As ABLH remained lower than 400 m from midday on 20 December to midday 21 December, the air pollution in the near-surface layer enhanced again. On the whole, the ABLH measured by ball sounding was typically below 400 m during heavy pollution episodes. Such shallow ABL led to high PM concentration in the near-surface layer.

Figure 8 Evolution of the atmospheric boundary layer height (ABLH) from the sounding data in Shenyang.
3.3 Impacts of air mass pathways

Different aerosol transportation paths were observed in the four cities during the haze event according to the 2-day back trajectory analyses at 200, 500, and 1000 m AGL (Fig. 9). The ending time of these back trajectory simulation is listed in Table 1, corresponding with the first peak when PM2.5 > 200 μg m –3. The air masses at 200 m AGL were generated from the Bohai Sea and passed through North Korea and southeastern Liaoning before arriving at Shenyang, Fushun, and Benxi. These air masses carried water vapor and led to high RH conditions. However, the air masses at all three altitudes in Anshan originated from western China (Fig. 9b), resulting in lower RH in Anshan than in the other three cities (Fig. 2j). The transportation paths of air masses at 1000 m AGL in Fushun and Benxi were similar to that of the air mass that arrived in Anshan. The air masses at the middle altitude of 500 m AGL in Shenyang originated from the North China Plain, where anthropogenic activities are intense and severe haze pollution is common. The transport pathways may demonstrate that in addition to the high local emissions of aerosols, the trans-boundary transport of aerosols from the North China Plain could play a role in exacerbating the air pollution.

Figure 9 The HYSPLIT 48-h backward trajectory at (a) Shenyang starting from 1500 LT 16 December 2016, (b) Anshan starting from 0300 LT 17 December 2016, (c) Fushun starting from 0600 LT 17 December 2016, and (d) Benxi starting from 1200 LT 17 December 2016.
4 Conclusions

The characteristics of the ABL structure during a persistent, severe regional haze event over the central Liaoning city cluster during 16–21 December 2016 were analyzed based largely on measurements from an ABL observational experiment and PM concentrations at the surface. The mean PM2.5 concentration during the haze event was larger in Anshan (158.3 ± 22.2 μg m–3) and Shenyang (157.9 ± 31.0 μg m–3) and smaller in Fushun (121.3 ± 48.9 μg m–3) and Benxi (135.6 ± 60.8 μg m–3). The distribution of daily horizontal wind fields at 10 m that were retrieved from the ECMWF reanalysis data during 16–21 December were strongly related to allocation of the surface pressure field. A southerly flow in the near-surface layer prevailed over the most regions of Liaoning. Wind speed decreased rapidly after 16 December, and the superposition of a weak wind field with weak horizontal wind resulted in air stagnation, ultimately causing deterioration in air quality during the haze event. A potential temperature inversion layer existed within the entire ABL at the earlier hazy episode, such stable atmospheric conditions were very favorable for the formation and maintenance of the haze pollution event. After, the inversion layer enhanced apparently, with the bottom varying between 300 and 900 m, corresponding to the fluctuation of PM concentration at surface. The ABLH measured by ball sounding was typically less than 400 m during a heavy pollution episode and exhibited a negative correlation with PM2.5. Such a low ABLH would concentrate pollutants in the near-surface layer.

The back trajectory analyses indicated that the air masses at the lower altitude of 200 m AGL in all cities originated from North Korea and passed over the Bohai Sea, except for the air masses arriving in Anshan, which caused lower humidity there compared with that in the other cities. The air masses at the higher altitudes of 500 and 1000 m AGL in Shenyang were generated from the North China Plain, where anthropogenic activities are intense and haze pollution is severe. In conclusion, the unfavorable meteorological conditions within ABL played as the external reason of this haze formation, besides the anthropogenic emissions of pollutants as the basic cause.

References
Chan, C. K., and X. H. Yao, 2008: Air pollution in mega cities in China. Atmos. Environ., 42, 1–42. DOI:10.1016/j.atmosenv.2007.09.003
Che, H., X. Xia, J. Zhu, et al., 2014: Column aerosol optical properties and aerosol radiative forcing during a serious haze–fog month over North China Plain in 2013 based on ground-based sunphotometer measurements. Atmos. Chem. Phys., 14, 2125–2138. DOI:10.5194/acp-14-2125-2014
Che, H. Z., H. J. Zhao, X. G. Xia, et al., 2015: Fine mode aerosol optical properties related to cloud and fog processing over a cluster of cities in Northeast China. Aerosol Air Qual. Res., 15, 2065–2081. DOI:10.4209/aaqr.2014.12.0325
Chen, J., C. S. Zhao, N. Ma, et al., 2012: A parameterization of low visibilities for hazy days in the North China Plain. Atmos. Chem. Phys., 12, 4935–4950. DOI:10.5194/acp-12-4935-2012
Chen, W. W., D. Q. Tong, M. Dan, et al., 2017: Typical atmospheric haze during crop harvest season in northeastern China: A case in the Changchun region. J. Environ. Sci., 54, 101–113. DOI:10.1016/j.jes.2016.03.031
Chen, Z. H., S. Y. Cheng, J. B. Li, et al., 2008: Relationship between atmospheric pollution processes and synoptic pressure patterns in northern China. Atmos. Environ., 42, 6078–6087. DOI:10.1016/j.atmosenv.2008.03.043
Cui, Y., C. Y. Zhao, T. Wang, et al., 2015: Temporal and spatial characteristics of haze days with different intensities and the climatic factors during 1961–2013 in Liaoning Province. Acta Scientiae Circumstantiae, 35, 1629–1637. DOI:10.13671/j.hjkxxb.2014.0917
Draxler, R. R., and G. D. Rolph, 2003: Hybrid Single-Particle Lagrangian Integrated Trajectory Model access via NOAA ARL READY. NOAA Air Resources Laboratory, Silver Spring, MD. http://www.ready.noaa.gov/HYSPLIT.php.
Fu, G. Q., W. Y. Xu, R. F. Yang, et al., 2014: The distribution and trends of fog and haze in the North China Plain over the past 30 years. Atmos. Chem. Phys., 14, 11949–11958. DOI:10.5194/acp-14-11949-2014
Fu, Q. Y., G. S. Zhang, J. Wang, et al., 2008: Mechanism of formation of the heaviest pollution episode ever recorded in the Yangtze River Delta, China. Atmos. Environ., 42, 2023–2036. DOI:10.1016/j.atmosenv.2007.12.002
Hu, X. M., Z. Q. Ma, W. L. Lin, et al., 2014: Impact of the Loess Plateau on the atmospheric boundary layer structure and air quality in the North China Plain: A case study. Sci. Total Environ., 499, 228–237. DOI:10.1016/j.scitotenv.2014.08.053
Hu, X. M., X. L. Li, M. Xue, et al., 2016: The formation of barrier winds east of the Loess Plateau and their effects on dispersion conditions in the North China Plains. Bound.-Layer Meteor., 161, 145–163. DOI:10.1007/s10546-016-0159-4
Ji, D. S., Y. S. Wang, L. L. Wang, et al., 2012: Analysis of heavy pollution episodes in selected cities of northern China. Atmos. Environ., 50, 338–348. DOI:10.1016/j.atmosenv.2011.11.053
Kim, S. W., S. C. Yoon, J. G. Won, et al., 2007: Ground-based remote sensing measurements of aerosol and ozone in an urban area: A case study of mixing height evolution and its effect on ground-level ozone concentrations. Atmos. Environ., 41, 7069–7081. DOI:10.1016/j.atmosenv.2007.04.063
Leng, C., Q. Zhang, D. Zhang, et al., 2014: Variations of cloud condensation nuclei (CCN) and aerosol activity during fog–haze episode: A case study from Shanghai. Atmos. Chem. Phys., 14, 12499–12512. DOI:10.5194/acp-14-12499-2014
Leng, C. P., J. Y. Duan, C. Xu, et al., 2016: Insights into a historic severe haze event in Shanghai: Synoptic situation, boundary layer and pollutants. Atmos. Chem. Phys., 16, 9221–9234. DOI:10.5194/acp-16-9221-2016
Li, M. N., S. J. Niu, S. T. Zhang, et al., 2015: Comparative study of turbulent characteristics between the fog day and haze day in Nanjing. Acta Meteor. Sinica, 73, 593–608. DOI:10.11676/qxxb2015.032
Liao, X. N., X. L. Zhang, Y. C. Wang, et al., 2014: Comparative analysis on meteorological condition for persistent haze cases in summer and winter in Beijing. Environ. Sci., 35, 2031–2044. DOI:10.13227/j.hjkx.2014.06.001
Liu, J., S. J. Fan, D. Wu, et al., 2015: Boundary layer characteristics of typical haze process in the Pearl River Delta region. China Environ. Sci., 35, 1664–1674. DOI:10.3969/j.issn.1000-6923.2015.06.008
Liu, N. W., Y. J. Ma, X. M. Liu, et al., 2011: Observations on haze and fog in Shenyang area. Acta Scientiae Circumstantiae, 31, 1064–1069.
Liu, S. H., Z. X. Liu, J. Li, et al., 2009: Numerical simulation for the coupling effect of local atmospheric circulations over the area of Beijing, Tianjin, and Hebei Province. Sci. China Ser. D Earth Sci., 52, 382–392. DOI:10.1007/s11430-009-0030-2
Liu, X. G., Y. H. Zhang, Y. F. Cheng, et al., 2012: Aerosol hygroscopicity and its impact on atmospheric visibility and radiative forcing in Guangzhou during the 2006 PRIDE-PRD campaign. Atmos. Environ., 60, 59–67. DOI:10.1016/j.atmosenv.2012.06.016
Lu, C. S., S. J. Niu, P. Yue, et al., 2011: Observational research on boundary layer structure during high incidence period of winter fog in Nanjing. Trans. Atmos. Sci., 34, 58–65.
Ma, Y. J., N. W. Liu, Y. F. Wang, et al., 2011: Thedistribution characteristics of fine particles and their impact on air quality in Shenyang and surrounding areas. Acta Scientiae Circumstantiae, 31, 1168–1174.
Ma, Y. J., N. W. Liu, Y. Hong, et al., 2012: The impacts on various particle sizes and the air quality caused by a dust weather process in spring 2011 in Liaoning. Acta Scientiae Circumstantiae, 32, 1160–1167.
Miao, Y. C., S. H. Liu, Y. J. Zheng, et al., 2015: Numerical study of the effects of local atmospheric circulations on a pollution event over Beijing–Tianjin–Hebei, China. J. Environ. Sci., 30, 9–20. DOI:10.1016/j.jes.2014.08.025
Miao, Y. C., J. P. Guo., S. H. Liu, et al., 2017: Classification of summertime synoptic patterns in Beijing and their associations with boundary layer structure affecting aerosol pollution. Atmos. Chem. Phys., 17, 3097–3110. DOI:10.5194/acp-17-3097-2017
Nielsen-Gammon, J. W., C. L. Powell, M. J. Mahoney, et al., 2008: Multisensor estimation of mixing heights over a coastal city. J. Appl. Meteor. Climatol., 47, 27–43. DOI:10.1175/2007JAMC1503.1
Peng, H. Q., D. Y. Liu, B. Zhou, et al., 2016: Boundary-layer characteristics of persistent regional haze events and heavy haze days in eastern China. Adv. Meteor., 6950154. DOI:10.1155/2016/6950154
Quan, J. N., Y. Gao, Q. Zhang, et al., 2013: Evolution of planetary boundary layer under different weather conditions, and its impact on aerosol concentrations. Particuology, 11, 34–40. DOI:10.1016/j.partic.2012.04.005
Shao, M., X. Y. Tang, Y. H. Zhang, et al., 2006: City clusters in China: Air and surface water pollution. Frontiers in Ecology and the Environment, 4, 353–361. DOI:10.1890/1540-9295(2006)004[0353:CCICAA]2.0.CO;2
Shen, X. J., J. Y. Sun, X. Y. Zhang, et al, 2015: Characterization of submicron aerosols and effect on visibility during a severe haze–fog episode in Yangtze River Delta, China. Atmos. Environ., 120, 307–316. DOI:10.1016/j.atmosenv.2015.09.011
Stock, M., Y. F. Cheng, W. Birmili, et al., 2011: Hygroscopic properties of atmospheric aerosol particles over the eastern Mediterranean: Implications for regional direct radiative forcing under clean and polluted conditions. Atmos. Chem. Phys., 11, 4251–4271. DOI:10.5194/acp-11-4251-2011
Sun, Y., T. Song, G. Q. Tang, et al., 2013: PM2.5 and boundary-layer structure during summer haze in Beijing . Atmos. Environ., 74, 413–421. DOI:10.1016/j.atmosenv.2013.03.011
Tang, G. Q., J. Q. Zhang, X. W. Zhu, et al., 2016: Mixing layer height and its implications for air pollution over Beijing, China. Atmos. Chem. Phys., 16, 2459–2475. DOI:10.5194/acp-16-2459-2016
Tang, G. Q., P. S. Zhao, Y. H. Wang, et al., 2017: Mortality and air pollution in Beijing: The long-term relationship. Atmos. Environ., 150, 238–243. DOI:10.1016/j.atmosenv.2016.11.045
Wang, H., M. Xue, X. Y. Zhang, et al., 2015a: Mesoscale modeling study of the interactions between aerosols and PBL meteorology during a haze episode in Jing–Jin–Ji (China) and its nearby surrounding region. Part 1: Aerosol distributions and meteorological features. Atmos. Chem. Phys., 15, 3257–3275. DOI:10.5194/acp-15-3257-2015
Wang, H., G. Y. Shi, X. Y. Zhang, et al., 2015b: Mesoscale modelling study of the interactions between aerosols and PBL meteorology during a haze episode in China Jing–Jin–Ji and its near surrounding region. Part 2: Aerosols’ radiative feedback effects. Atmos. Chem. Phys., 15, 3277–3287. DOI:10.5194/acp-15-3277-2015
Wang, K., H. S. Zhang, Q. Wang, et al., 2006: Study on the structure and evolvement of atmospheric boundary layer of frontal fogs in spring and winter at southern suburb of Beijing. Acta Scientiarum Naturalium Universitatis Pekinensis, 42, 55–60. DOI:10.3321/j.issn:0479-8023.2006.01.010
Wang, L. L., N. Zhang, Z. R. Liu, et al., 2014a: The influence of climate factors, meteorological conditions, and boundary-layerstructure on severe haze pollution in the Beijing–Tianjin–Hebei region during January 2013. Adv. Meteor., 2014, 685971. DOI:10.1155/2014/685971
Wang, P., H. Z. Che, X. C. Zhang, et al., 2010: Aerosol optical properties of regional background atmosphere in Northeast China. Atmos. Environ., 44, 4404–4412. DOI:10.1016/j.atmosenv.2010.07.043
Wang, S., T. T. Liao, L. L. Wang, et al., 2016: Process analysis of characteristics of the boundary layer during a heavy haze pollution episode in an inland megacity, China. J. Environ. Sci., 40, 138–144. DOI:10.1016/j.jes.2015.12.008
Wang, Y. Q., X. Y. Zhang, J. Y. Sun, et al., 2015c: Spatial and temporal variations of the concentrations of PM10, PM2.5 and PM1 in China . Atmos. Chem. Phys., 15, 13585–13598. DOI:10.5194/acp-15-13585-2015
Wang, Y. S., L. Yao, L. L. Wang, et al., 2014b: Mechanism for the formation of the January 2013 heavy haze pollution episode over central and eastern China. Sci. China Earth Sci., 57, 14–25. DOI:10.1007/s11430-013-4773-4
Wu, D., X. Y. Bi, X. J. Deng, et al., 2007: Effect of atmospheric haze on the deterioration of visibility over the Pearl River Delta. Acta. Meteor. Sinica, 21, 215–223.
Wu, Z. L., A. X. Liu, C. C. Zhang, et al., 2009: Vertical distribution feature of PM2.5 and effect of boundary layer in Tianjin . Urban Environment & Urban Ecology, 22, 24–29.
Yang, J., L. Wang, R. Y. Liu, et al., 2010: The boundary layer structure and the evolution mechanisms of a deep dense fog event. Acta Meteor. Sinica, 68, 998–1006. DOI:10.11676/qxxb2010.094
Ye, X. X., Y. Song, X. H. Cai, et al., 2016: Study on the synoptic flow patterns and boundary layer process of the severe haze events over the North China Plain in January 2013. Atmos. Environ., 124, 129–145. DOI:10.1016/j.atmosenv.2015.06.011
Zhang, R. H., Q. Li, and R. N. Zhang, 2014: Meteorological conditions for the persistent severe fog and haze event over eastern China in January 2013. Sci. China Earth Sci., 57, 26–35. DOI:10.1007/s11430-013-4774-3
Zhang, Y. H., Y. J. Ma, N. W. Liu, et al., 2010a: Study on environmental characteristics of haze in Shenyang. Ecology and Environmental Sciences, 19, 2636–2641. DOI:10.3969/j.issn.1674-5906.2010.11.020
Zhang, Y. H., H. B. Yang, X. D. Zou, 2010b: Test research on the characteristics of pollution boundary layer in the river outlet area of Grand Liao River. Meteor. Environ. Res., 1, 57–61.
Zhao, H. J., H. Z. Che, X. Y. Zhang, et al., 2013a: Aerosol optical properties over urban and industrial region of Northeast China by using ground-based sun-photometer measurement. Atmos. Environ., 75, 270–278. DOI:10.1016/j.atmosenv.2013.04.048
Zhao, H. J., H. Z. Che, X.Y. Zhang, et al., 2013b: Characteristics of visibility and particulate matter (PM) in an urban area of Northeast China. Atmos. Pollut. Res., 4, 427–434. DOI:10.5094/APR.2013.049
Zheng, G. J., F. K. Duan, H. Su, et al, 2015: Exploring the severe winter haze in Beijing: The impact of synoptic weather, regional transport and heterogeneous reactions. Atmos. Chem. Phys., 15, 2969–2983. DOI:10.5194/acp-15-2969-2015
Zhu, J. L., H. Liao, and J. P. Li, 2012: Increases in aerosol concentrations over eastern China due to the decadal-scale weakening of the East Asian summer monsoon. Geophys. Res. Lett., 39, L09809. DOI:10.1029/2012GL051428