Dome effect of black carbon and its key influencing factors : A one-dimensional modelling study

Black carbon (BC) has been identified to play a critical role in aerosol-planet boundary layer (PBL) interaction and further deterioration of near-surface air pollution in megacities, which has been named as its “dome effect”. However, the impacts of key factors that influence this effect, such as the vertical distribution and aging processes of BC, and also the underlying land surface, have not been quantitatively explored yet. Here, based on available in-situ measurements of meteorology and 15 atmospheric aerosols together with the meteorology-chemistry online coupled model, WRF-Chem, we conduct a set of parallel simulations to quantify the roles of these factors in influencing the BC’s dome effect and surface haze pollution, and discuss the main implications of the results to air pollution mitigation in China. We found that the impact of BC on PBL is very sensitive to the altitude of aerosol layer. The upper level BC, especially those near the capping inversion, is more essential in suppressing the PBL height and weakening the turbulent mixing. The dome effect of BC tends to be significantly intensified 20 as BC aerosol mixed with scattering aerosols during winter haze events, resulting in a decrease of PBL height by more than 15%. In addition, the dome effect is more substantial (up to 15%) in rural areas than that in the urban areas with the same BC loading, indicating an unexpected regional impact of such kind of effect to air quality in countryside. This study indicates that China’s regional air pollution would greatly benefit from BC emission reductions, especially those from the elevated sources from the chimneys and also the domestic combustions in rural areas, through weakening the aerosol-boundary layer 25 interactions that triggered by BC.


Introduction
Air pollution, particularly haze pollution, has been one of the key environmental challenges for China, especially in the more developed regions like northern and eastern China (Wang et al., 2017).Haze pollution in these regions is generally characterized by extremely low visibility and dramatically rising surface aerosol concentration (Cai et al., 2017;Zhao et al., 2011).For example, in January 2013, a longlasting episode of severe haze occurred in central and eastern China with a maximum PM 2.5 (particles with dynamic diameter less than 2.5 µm) mass concentration in Beijing reaching up to 680 µg m −3 (Wang et al., 2014a, c, d).In addition to deteriorating air quality in megacities like Beijing, large-scale regional haze pollution also covers rural and suburban areas (Xu et al., 2011;Chen and Wang, 2015).During serious haze pollution in that month, thick haze engulfed 1.4 million km 2 of land area, affecting up to 800 million people in 17 provinces.Such severe and aggravating regional haze pollution has triggered extensive public alarm due to PM 2.5 -associated adverse health effects (e.g., cardiovascular and respiratory diseases) (Kim et al., 2015;Mauderly and Chow, 2008;Gao et al., 2015a;Pope et al., 2002).Consequently, air pollution mitigation has been one of the top priority for China's central and local governments.However, even though a series of control measures, such as the Action Plan on Prevention and Control of Air Pollution, have been carried out in order to reduce emissions and mitigate fine particle pollution (Zhang et al., 2016), the frequencies of severe pollution events have continued to increase in recent years and the intensity of pollution episodes has not shown a sig-Z.Wang et al.: Dome effect of black carbon and its key influencing factors nificant decrease yet (Ding and Liu, 2013;Niu et al., 2010;Cai et al., 2017).
During hazy days, concentration of BC, one of the most important aerosol components from both environmental and climate perspectives, could exceed 20 µg m −3 in China's megacities (Yang et al., 2007;Sun et al., 2014), higher than those in other regions across the world.Such a high level of BC concentration in China is primarily attributed to intensive residential combustions and a coal-dominant energy structure (Qin and Xie, 2011;Zhang et al., 2009b).With a high light-absorbing efficiency, BC can exert a substantial impact on climate change on a regional or even global scale (Bond et al., 2013;Menon et al., 2002;Ramanathan and Carmichael, 2008).At the same time, BC has been proven to be inextricably linked to short-term changes in public health such as cardiovascular mortality and cardiopulmonary hospital admissions (Eklund et al., 2014;Janssen et al., 2011).A recent study by Ding et al. (2016) has revealed the vital role of BC in enhancing near-surface haze pollution by the combined effects of heating by the lightabsorbing BC aerosols in the upper PBL and the reduction of surface heat flux, which substantially suppresses the development of the PBL and consequently causes extreme haze pollution episodes in China's megacities.Such an effect has been referred to as the "dome effect" of BC (Ding et al., 2016).A similar effect has also been found for BC over the Indian Ocean (Wilcox et al., 2016) and dust aerosols in northern and eastern China (Liu et al., 2016;Yang et al., 2016).
In terms of the aerosol-PBL interaction induced by BC, the vertical distribution of the aerosol layer is expected to be an important influencing factor since the upper-level BC can alter the air temperature stratification much more efficiently (Ding et al., 2016).Located in one of the main monsoon regions, East China is subjected to large-scale monsoon circulations, where frequent cyclones, fronts and convections tend to lift near-surface air pollutants to the middle and even upper troposphere (Ding et al., 2009;Zhang et al., 2009a).In addition, elevated sources like power plant plumes and biomass burning smoke and subsequent long-range transport of air masses frequently lead to vertical heterogeneity of BC profile (Ding et al., 2013;Huang et al., 2016;Yang et al., 2015;Guinot et al., 2006;J. M. Chen et al., 2017).Thus, it is of great importance to understand the sensitivity of the dome effect to vertical distribution of BC in China.On the other hand, during haze events, it is highly possible that the coexistence of concentrated air pollutants and complex physicochemical interactions among them to lead to a dramatic increase in secondary aerosols (Huang et al., 2014a, b).Hence, freshly emitted BC has usually been observed to undergo notable aging and hygroscopic growth and become almost internally mixed with scattering secondary aerosols like sulfate during hazy days (Bond et al., 2013;Cui et al., 2016;Huang et al., 2013), thereby remarkably enhancing its light-absorbing properties (Peng et al., 2016;B. Chen et al., 2017;Cappa et al., 2012;Yang et al., 2012;Shen et al., 2017).
Subsequently, the dome effect of BC might be significantly strengthened by mixing with scattering aerosols.However, such an influence of BC on the dome effect has not been quantitatively examined yet.As mentioned above, haze pollution has generally occurred on a regional scale in East and North China.Previous studies have highlighted the importance of aerosol-PBL interaction in cities (Petaja et al., 2016;Ding et al., 2016;Wang et al., 2014a;Gao et al., 2015b;Cai et al., 2017;Li et al., 2017b).It is noteworthy that more than half of the population lives in rural areas, with a majority in the plain areas in North and Central China, which has also been exposed to fine particulate pollution.The rural areas in East China are usually covered by cropland with different land-surface properties to urban areas.The difference of the dome effect over the two regions with distinct land cover categories remains to be further explored.In consequence, this study aims at identifying and quantitatively assessing the dependence of the dome effect on several key factors, including vertical distribution, aging of BC and the different underlying land surfaces (i.e., urban and rural areas).
Numerical simulation with meteorology-chemistry online coupled models has served as a practicable and effective way to characterize the aerosol radiative effect and its impact on PBL evolution (Grell et al., 2005;Baklanov et al., 2014;Yu et al., 2002).To disentangle the impacts of various factors on aerosol-PBL interaction, one-dimensional meteorology-chemistry online coupled model is applied in this study for the purpose of excluding influences from, for example, synoptic processes and regional transport.Additionally, one-dimensional modeling with high vertical resolution enables better representation of PBL evolution and also allows flexible initial and boundary conditions.Therefore, in the present work, a single-column version of WRF-Chem (Weather Research and Forecasting model coupled with Chemistry) driven by available observations is employed to investigate the aerosol-PBL interaction and its influencing parameters.The paper is organized as follows.Section 2 discusses major aspects of model simulations, including the WRF-Chem model and its configurations, parameterizations of aerosol properties, initial and boundary meteorological and chemical conditions, and the design of numerical experiments.The dome effect due to BC and its dependence on different conditions are presented and discussed in Sect.3. Specifically, discussions in Sect. 3 include the effects of vertical distribution and aging process of BC, as well as different underlying land surfaces.Main results and the possible implications in future policy of air pollution mitigation are discussed in Sect. 4.
2 Data and method

Model configuration
The simulations were conducted with the WRF-Chem version 3.6.1 (Grell et al., 2005) single-column model (SCM).Except for advection, the physical and chemical processes are exactly the same as in the three-dimensional version, which is a fully coupled online meteorology-chemistry model including emission and deposition of pollutants, gaseous and aqueous chemical transformation, aerosol chemistry, and dynamics.The WRF-Chem SCM runs on a 3 × 3 grid with periodic lateral boundary conditions in both zonal and meridional directions.We used a spatial resolution of 4 km and a vertically stretched sigma coordinate with the model top set at a constant pressure, corresponding to about 6000 km.One hundred vertical levels were placed equidistantly with a height of 60 m from the ground surface to model top to better resolve the vertical structure of the atmosphere and ensure identical BC mass loadings.
The parameterization schemes were selected following the work by Ding et al. (2016), in which the model configuration showed good performance on boundary layer meteorology.The RRTMG short-and longwave radiation schemes (Iacono et al., 2008) were used to couple with aerosol scheme in order to reproduce aerosol-radiation interactions.The YSU nonlocal-K boundary layer scheme (Hong et al., 2006) and the Noah land surface scheme (Tewari et al., 2016) were applied for boundary layer evolution and land-atmosphere interactions, respectively.For representation of cloud and precipitation processes the Lin microphysics scheme (Lin, 1983) together with the Grell-Deveny cumulus parameterization (Grell and Devenyi, 2002) were employed.The CBMZ (Carbon Bond Mechanism version Z) photochemical mechanism combined with MOSAIC (Model for Simulating Aerosol Interactions and Chemistry) aerosol model (Fast et al., 2006;Zaveri and Peters, 1999) was applied to represent atmospheric chemistry.Aerosols were assumed to be spherical particles.The size distribution was divided into four discrete size bins defined by their lower and upper dry particle diameters (0.039-0.156, 0.156-0.625,0.625-2.5 and 2.5-10.0µm).Aerosol water content is calculated using the Zdanovskii-Stokes-Robinson method according to relative humidity (Zaveri et al., 2008).In each bin, aerosols are presumed to be spherical and internally mixed.The volumeaveraged mixing rule is adopted to calculate the bulk refractive index by considering both aerosol water content and chemical compositions, including BC, sulfate, nitrate, and ammonium.Mie theory is then used to calculate the extinction coefficient, single-scattering albedo (SSA) along with asymmetry factor for each bin.Model domains and configuration selections are summarized in Table 1.
The WRF-Chem SCM was initialized with monthly averaged radio-sounding measurements at 12:00 UTC taken in Beijing (116.28 • E, 39.93 • N) to represent the typical atmo- The soil temperature profile was taken from the monthly averaged WRF regional simulation results for December 2013.
The initial condition of atmospheric pollutants, which is mainly concerned with BC will be given in details in the following section.Vertical mixing of aerosol is switched off to ensure that the concentration and altitude of aerosol layer do not vary with PBL evolution.Each numerical experiment was conducted for a time period of 72 h with the first 48 h as model spin-up time.

Design and analysis of numerical experiments
We conducted a set of multidimensional experiments.Each individual experiment contains hundreds of simulations and can be divided into two groups.These two experiment groups shared exactly the same model settings and configurations except that one was with aerosol radiation interaction (ARI) and the other without (noARI).Each group consists of several parallel experiments with various initial conditions of airborne pollutants or surface parameters.To investigate the impact of heterogeneity in vertical distribution of BC on the dome effect, BC plumes with a width of 300 m were settled at different altitudes, ranging from the ground surface to about 2000 m.BC concentration was assumed to be 0-30 µg m −3 according to existing field measurements across China (Zhang et al., 2013;Sun et al., 2004;Zhao et al., 2013a).BC concentrations in the plume were presumed to be a Gaussian distribution vertically, with a maximum value ranging from 0 to 30 µg m −3 in the central axis, the altitude of which also ranges from 150 to 2250 m with an interval of 300 m.Furthermore, to better understand the difference of the dome effect over different underlying ground surfaces (i.e., urban and rural areas) and its roles in regional air pollution, 18.9 × 10 5 25.0 × 10 5 we conducted similar parallel numerical experiments over both urban and rural land surfaces.The distinct differences in important surface parameters between them are listed in Table 2.It is well known that BC emission is mainly related to fossil fuel combustion in China; hence, it is usually coemitted with other gaseous pollutants like SO 2 and NO x and then mixed with their oxidation products, i.e., sulfate and nitrate, during transportation and aging processes (Wang et al., 2014b;Huang et al., 2013;Cheng et al., 2006).To quantify the impact of mixing with these scattering aerosols, we also performed another two parallel simulations with different SNA (sulfate, nitrate and ammonium) aerosol concentration levels in the BC plume.The mass ratio of BC to SNA was derived from year-round measurement in Beijing (Zhang et al., 2013).
3 Results and discussion 3.1 One-dimensional modeling of the dome effect of BC The dome effect of black carbon was first revealed by Ding et al. (2016) based on integrating online coupled regional simulation and corresponding observations during the winter haze event in December 2013, when severe PM 2.5 pollution covered East China with hourly concentrations up to ∼ 900 µg m −3 and visibility less than 100 m (Zheng et al., 2015).In this study, it was proven that BC plays an important role in aggravating haze pollution via aerosol-PBL interaction.To clearly demonstrate BC-induced aerosol-PBL interaction, we selected one typical episode during 23-24 December 2013 in Beijing and initialized the WRF-Chem SCM with an averaged BC profile taken from three-dimensional simulation results for the same case by Ding et al. (2016).
The regional WRF-Chem simulation of winter haze in December 2013 in Ding et al. (2016) had been verified to perform well in capturing the temporal variations of BC.The inputted BC vertical profile in Fig. 1 indicated that, during the abovementioned case, BC concentration reached up to approximately 30 µg m −3 near the ground surface and decreased rapidly with altitude to a relatively constant value of less than 5 µg m −3 above 800 m.However, although the BC profile featured its maximum concentration near the surface, the heating efficiency of BC due to shortwave radiation absorption peaked around 600-800 m, indicating that BC in the upper PBL is more efficient in terms of absorbing shortwave radiation and heating surrounding air masses.Such a heating profile is partly caused by the fact that incident solar radiation is attenuated by aerosols, trace gases and cloud when transferring within the atmosphere.Therefore, BC at higher altitudes tended to be subjected to higher incident radiation flux and absorb more solar energy.At the same time, lower air density caused the upper air to be more readily heated compared with that near the surface.Due to light absorption caused by BC together with additional extinction caused by scattering aerosols, solar radiation reaching the surface was diminished to a relatively large extent, resulting in less sensible heat flux and thus lower temperature in the near-surface atmosphere.
Overall, upper-level warming and surface cooling substantially modified the temperature stratification.As illustrated in Fig. 2, the afternoon upper-level heating and morning surface cooling could be clearly identified.The strongest upperair warming had remained at over 1.0 • C between 1000 and 1200 m since midday (12:00 LT) because the incident shortwave radiation was most intense and the heating effect had already been accumulated before midday.The strongest cooling effect appeared in the morning (09:00-11:00 LT) and decreased very sharply after noon.The reason is that, with in- creasingly intense turbulence during morning boundary evolution, part of the surface cooling was gradually compensated for by the entrainment of more warmed air into the turbulent PBL from above when the PBL was developed in the late morning.The warming and cooling of different atmospheric levels remarkably altered the stratification, thereby weakening convective motions.Stable stratification combined with decreased sensible heat flux at the ground surface greatly suppressed vertical turbulence in the boundary layer (Wilcox et al., 2016), contributing to a delay of PBL development and an earlier drop as well as a substantial decrease in PBL height, which hindered the air pollutants from being further dispersed vertically.These modifications in temperature stratification were relatively extensive in East China during winter according to the simulation conducted in Ding et al. (2016).For the whole month of December in 2013, the occurrence probability of enhanced atmospheric stability (PBL height decreasing over 10 %) due to the dome effect could reach up to 66, 73, 69 and 66 % in Zhengzhou, Shijiazhuang, Shanghai and Nanjing, respectively (Ding et al., 2016).

Impacts of altitude and concentration of BC aerosol layer
As previously mentioned, specific synoptic conditions, chimney plumes and regional transport of air pollution could result in non-uniformly distributed pollutant profiles (Xu et al., 2014;Guinot et al., 2006;Ding et al., 2009).Vertically inhomogeneous distribution of BC aerosol has been frequently observed by in situ aircraft and tethered balloon measurements (Zhao et al., 2015;Li et al., 2015).Given that BC was measured to rise around 800-1000 m during daytime in a haze episode (Li et al., 2015) and the upper-PBL BC has higher light-absorbing efficiency where solar heating usually maximizes, a BC plume with a maximum concentration of 10 µg m −3 at the central axis is taken as a typical example to illustrate the perturbations on detailed physical processes related to PBL evolution due to a vertically non-uniform BC profile.As displayed in Fig. 3, with the existence of absorptive BC, incident solar radiation at the surface is diminished by about 5.9 W m −2 at 12:00 LT, leading to a decline of 0.2 • C in surface temperature and a decline of 3.6 W m −2 in surface sensible heat flux, which would otherwise heat the lower atmosphere, promote turbulent motions and enhance the PBL development.At the same time, absorbed radiation energy is converted to thermal energy and leads to substantial heating near the top of the PBL, as reflected by a rise in air temperature around 1.0 • C. Accordingly, modified temperature profile results in a more stable stratification with its exchange coefficient falling by 15 % (Fig. 3b), i.e., less turbulence mixing.Similarly, a substantial decrease in PBL height is attributed to less surface heat flux accompanied by stable stratification, which is showed in Fig. 3c.
To shed more light on the impacts of various BC vertical distributions on the dome effect, we conducted hundreds of parallel numerical experiments by changing the altitude and concentration magnitude of BC from the simulation shown in Fig. 3 to a larger extent in order to figure out how the PBL evolution responds to varying BC profiles and mass loading.Specifically, we manually increase BC concentrations from 0 to 30 µg m −3 with an increment of 2 µg m −3 , all of which are placed at altitudes from 150 to 2250 m with an increment of 300 m.Thus, the multidimensional experiment consists of 105 parallel numerical simulations.All the individual simulations (marked by black dots in Fig. 4) and resultant perturbations on PBL height and the turbulence coefficient under different altitude and concentration loadings of BC are shown in Fig. 4. As several studies have pointed out, BC near the ground surface warms the Earth-atmosphere system and favors the development of the PBL by trapping more solar radiation that should be reflected by land and then promoting convective motions (Huang et al., 2015;Barbaro et al., 2014).The increased air temperature near the surface weakens the capping inversion, acting to cancel the effect of reduced buoyancy flux at the ground surface and raising the top of the PBL (Yu et al., 2002).However, this enhancement in PBL development is not that noticeable in magnitude while compared with PBL suppression due to BC at higher altitude.Consistent with several existing observational and numerical studies, absorbing smoke aloft is capable of remarkably changing the energy balance between the surface and the atmosphere in a way that stabilizes the boundary layer and suppresses convection (Koren et al., 2004;Ackerman et al., 2000).As demonstrated in Fig. 4, the PBL top could be decreased by about 15 % and the turbulent exchange coefficient dropped by over 20 % due to high-altitude BC plume.The substantial suppression effect induced by the BC plume maximizes around the top of the PBL and can be attributed to two main reasons.Firstly, incident solar radiation flux at the top of the PBL is usually most intense.After being absorbed by BC, it can effectively heat the surrounding air and change the strength of the capping inversion.Secondly, the turbulent exchange coefficient falls to a relatively small value when close to the top of the PBL (Fig. 3b), indicating that vertical turbulence exchange for heat at this height is rather weak.Therefore, the warming layer is prone to being kept at that level rather than diffusing vertically, further strengthening the capping inversion and consequently lowering the PBL height.
Another interesting fact is that, while BC layer is located above the original top of the PBL (i.e., above the 1300 m in Fig. 4), the boundary layer also becomes shallower but relatively less notable compared with upper-PBL BC layer (i.e., 600-1000 m).In this case, BC also blocks part of the incoming solar radiation and diminishes surface fluxes proportionally; nevertheless, it no longer alters the temperature stratification of the PBL below.It is well known that changes in the height of the PBL are determined by the surface buoyancy flux and the capping inversion, and both of them are affected by BC plumes (Ding et al., 2017(Ding et al., , 2016)).Since the impact of column loading of BC to surface buoyancy flux will not change substantially in the lower troposphere, here the different impact of BC on PBL height around the capping inversion (i.e., 600-1200 m) with the above altitudes indicates that the upper-air warming due to BC layer plays a more dominant role in depressing the PBL height and hence enhancing the near-surface air pollution.

Amplified dome effect by mixing with scattering aerosols
During wintertime heavily polluted episodes, when both primary and secondary pollutants increase dramatically, BC will be coated by scattering aerosols through condensation of low-volatility gases and coagulation with secondary aerosols and thus becomes hydrophilic and more internally mixed (Shiraiwa et al., 2007;Chen et al., 2016).Theoretically, its absorption properties can be significantly enhanced by internal mixing with other compounds because the coatings act as a lens and enlarge its mass absorption cross section effectively (Bond et al., 2013;B. Chen et al., 2017).Here, although organic matter also grows rapidly during hazy days, we do not include it in this work due to the uncertainties in its optical properties as well as less notable hygroscopicity.We analyzed this absorption amplification effect based on consistent in situ measurements during the development of severe haze pollution in Beijing in January 2013 (Sun et al., 2014).Compared with the clean period, the secondary aerosol increased significantly during haze episodes.Take sulfate aerosol for instance, its concentrations during hazy days were approximately 55 times higher than those under clean condition.Based on the ratio of BC to SNA during different periods and assuming BC concentration is 5 µg m −3 (Zhang et al., 2013), absorption at the simulated wavelength (i.e., 300, 400, 600 and 1000 nm) was amplified by factors of 1.8, 1.7, 1.6 and 1.4, respectively (Fig. 5).These absorption coefficient amplification factors are comparable with in situ observations during clean period and haze event in the North China Plain (B.Chen et al., 2017).To investigate the relative importance of amplified absorption for each band, the solar spectral irradiance at sea level is also given for reference.The absorptive magnification for radiation at the wavelength of 400-600 nm tends to have a dominant effect since the solar irradiance peaks in this band.Accordingly, shortwave heating rate also increased from about 0.15 to over 0.22 K h −1 , which considerably accelerates the warming effect induced by BC (Fig. 5).
In addition to directly increasing the absorption cross section, SNA could undergo deliquescence and lead to an increase in aerosol diameter to a large extent, especially under humid conditions, indirectly enhancing the light-absorbing Figure 5. Aerosol absorption coefficient at different wavelengths (300, 400, 600, and 1000 nm) and shortwave heating rate at 12:00 LT for runs during a clean period and haze episode when the ratio of BC and SNA is 1 : 3 and 1 : 8, respectively.The yellow shaded area schematically shows solar spectral irradiance at sea level.
capacity (Tsai and Kuo, 2005;Zheng et al., 2015;Liu et al., 2011).To comprehensively understand this enhancement effect induced by SNA, we conducted tens of sensitivity experiments under different levels of SNA concentration and relative humidity (RH).As presented in Fig. 6, aerosol absorption extinction coefficient increases with SNA concentration.This dependence on SNA is much notable under lower RH (50 %) before aerosol starts to deliquesce.By contrast, when the air becomes more humid, aerosol water uptake becomes increasingly more important.According to year-round observational data (Zhang et al., 2012), SNA concentrations and humidity condition vary a lot across China.For typical cities in China, the annual-averaged SNA level ranges from 40-60 µg m −3 in Nanjing (118.95 • E, 32.12 • N) and Shanghai (121.45 • E, 31.22 • N) to almost 140 µg m −3 in Beijing (116.30• E, 39.99 • N), but the annual mean BC concentrations are approximately 5 µg m −3 (Ye et al., 2003;Zhang et al., 2013).In Nanjing, with a lower SNA level but higher RH, SNA-induced hygroscopic growth may play the dominant role in the enhancement of aerosol light-absorbing efficiency.Instead, relatively higher SNA concentration in Beijing makes the main contribution according to Fig. 6.Specifically, the absorption could be elevated by 63 % when RH rises from 50 to 90 % in Nanjing, while for Beijing the corresponding enhancement is only 18 %.Such disparities indicate that aerosol absorption and further impact on boundary layer evolution have a closer link to humidity and will be more important in coastal region at lower latitudes in China.
Enhanced light absorption and heating efficiency by mixing with SNA certainly perturb aerosol-PBL interactions (Fig. 7).Overall, the PBL responses tend to be magnified in terms of both increased PBL top by near-surface BC and suppressed convections by upper-level BC.For typical wintertime meteorological conditions in Beijing, aging processes could lead to a decline of 5 W m −2 in surface heat flux and an additional 15 % decrease in the PBL top at a maximum when compared with those caused by freshly emitted BC.What is more, when BC concentration is fixed at 5 µg m −3 and SNA concentration exceeds 80 µg m −3 , this absorption enhancement leads to the phenomenon that even BC at lower altitudes may increase the stratification and exert a negative effect on PBL development, as shown in Fig. 7.However, this abrupt change is not observed in simulations with only SNA aerosols included.In other words, there probably exists a critical point of the BC/SNA ratio for its associated dome effect, from which the impact of lower-level BC on PBL suppression would be oppositely attributed to a large decrease in surface sensible heat flux.

A comparison of the dome effect over urban and rural areas
In the past, BC-induced aerosol-PBL interaction has mostly been investigated for cities.However, urban areas only take up approximately 5 % of the area in China, while cropland is the dominant surface type and also features dense population and high PM 2.5 exposure (Liu, 2005;Xu et al., 2011).Frequent regional-scale haze pollution and cropland-dominant land cover in China make it necessary to understand this interaction over rural areas (Yang et al., 2015;Zhao et al., 2013b).The differences in surface parameters between cities and rural area are listed in Table 2. Comparatively, rural sur- 3 : NH + 4 = 3 : 2 : 1, which is based on year-round observations in Beijing (Zhang et al., 2012).The x and y axes are the same as Fig. 4. faces, including cropland and pasture, are characterized by larger surface albedo, higher soil moisture and greater heat capacity.
The radiation energy balance without any influence from aerosol already shows great disparities for the two kinds of land surface.Under the same level of incident solar radiation intensity, net shortwave radiation at the ground surface is mainly determined by land surface albedo.In comparison with urban areas, less net downward solar radiation on rural surfaces is ascribed to larger albedo reflecting more solar radiation, accompanied by higher upwards shortwave radiation flux at the top of the atmosphere, as shown in Fig. 8. Correspondingly, the surface temperature in rural areas at noon is almost 1.5 • C colder than that in urban areas, indicating a strong urban heat island effect (Oke, 1982).As for the atmospheric heating, there is little difference over both surface types when excluding the radiative perturbation of aerosols.The resultant height of the PBL is expected to be 250 m lower in rural areas.
If taking the aerosol radiative effect into account when BC is 10 µg m −3 around an altitude of 1000 m, over both urban and rural areas, net downward shortwave radiation flux, surface temperature and upwelling shortwave radiation show a notable decrease.At the same time, atmospheric heating rate in the upper air increases in response.The most distinctive difference between these two kinds of land surface is that upwelling shortwave radiation flux at the top of the atmosphere over rural surfaces shows a reduction twice high as over urban surfaces, indicating that a larger portion of so-  lar energy is blocked in the atmosphere, which corresponds to a much stronger atmospheric heating rate due to radiation absorption.At the same time, it suggests a slightly larger decrease in surface temperature over rural surfaces, followed by a lower sensible heat transfer into the atmosphere.The greater surface cooling together with more intense upper warming jointly enhances stable stratification.Therefore, the decrease in the PBL top in rural areas is double that in cities, as shown in Fig. 8.
A comprehensive analysis of responses of the PBL to various magnitudes and altitudes of BC aerosol over urban and rural land surfaces is illustrated in Fig. 9.As mentioned be-fore, since the PBL is less well developed in rural area, the altitude of the BC plume with maximum suppression effect on the PBL was slightly lower than that over cities.Additionally, BC-induced heating is more intensive over rural surfaces due to the fact that larger albedo results in more shortwave radiation reflected by the ground surface, part of which is absorbed by BC aerosol when transferring upwards.Consequently, the promotion and suppression effect on the PBL attributed to BC aerosol over rural surfaces is about 15 % larger than that over urban surfaces.That is, BC would exert a more intense dome effect over rural surfaces.It should be noted that residential combustion of raw coal and biofuel Z. Wang et al.: Dome effect of black carbon and its key influencing factors in a small domestic stove in rural area is responsible for the majority of BC emission in China (Zhi et al., 2008;Li et al., 2017a).A more stable PBL and more significant dome effect over rural surfaces would further favor the formation and accumulation of regional air pollution.

Summary and implications
The dome effect of BC plays a crucial role in air pollution deterioration, which is expected to be highly dependent on many factors like vertical distribution, mass loading, and aging processes of BC as well as the underlying land surface.By integrating available in situ observations of meteorology and atmospheric aerosols together with simulations using meteorology-chemistry online coupled model WRF-Chem SCM, we conduct a set of multidimensional experiments, each of which contains hundreds of parallel simulations, to quantify the impacts of these factors in enhancing the surface haze pollution.We found that the dome effect of BC is extremely sensitive to the altitude of aerosol layer.In more detail, near-surface BC tends to promote the PBL development, while BC aloft is more essential in PBL suppression, especially near the capping inversion.At the same time, this work indicates that the dome effect of BC can be significantly intensified when BC gets internally mixed with scattering aerosol during winter haze events, which could further decrease PBL height by up to 15 %.In terms of different underlying surface, the dome effect is more substantial in rural areas.Under the same conditions of BC, the PBL top decrease in rural areas could be 15 % greater than the corresponding value in urban areas.
Haze pollution in China is becoming increasingly frequent and thus has received extensive attention from the public and government.In spite of a series of emission control measures, the air quality has yet to be improved substantially, especially for the extreme haze events in winter.Our study highlights the importance of BC in worsening air pollution and further provides clues for both long-term policy making on air quality improvement and short-term pollution emergency plan.Our study has shown that the long-range transport of absorption aerosols like BC could cause a larger impact on regional air quality.In fact, the elevated emission sources, i.e., point sources with chimneys like coal-fired power plants that are typically higher than 200 m, could be easily transported and deteriorate dispersion conditions by suppressing the PBL on larger scales.Therefore, an efficient way to mitigate nearsurface air pollution is to preferentially reduce BC emission from these sources before and during regional severe haze pollution.In the longer run, clean energy substitution and extension of BC emission reduction technologies should be encouraged.In addition, it is noteworthy that residential combustion is responsible for the majority of pollution sources in the countryside/rural areas and is the most important contributor to BC emissions in China, which could be lifted to the upper PBL through convective motions.Considering the significant dome effect of BC over rural surfaces, the residential combustion sources would further favor the formation and deterioration of regional air pollution.As a result, it can be implied that abatement technology for domestic stoves in rural areas could be another cost-effective way to reduce regional air pollution.Of course, as one of the key short-lived climate forcers, BC has been found to be one of the most important components contributing to global warming.These efforts devoted to reducing BC emission in both urban and rural areas in China definitely result in a substantial reduction in the national carbon emission and co-benefit the mitigation of global warming.
In addition, our study here highlights the significance of investigating the vertical structure of atmospheric processes in detail in the lower atmosphere, which is very important for improving the understanding of the interaction of atmospheric physics and chemistry.However, in China, most of the existing and current efforts of field measurements, including both routinely operated monitoring networks and research based field measurement stations, have been mainly focused on the ground surface.There is very limited information on the vertical profiles of atmospheric aerosols and physical parameters related to key processes in the PBL.For numerical models, there is also a lack of vertical datasets for model evaluation and parameterization improvement.To gain a more comprehensive understanding of the causes and evolution of air pollution, more vertical measurements -including in situ measurements based on aircraft and tagged balloon platforms and remote sensing measurements using ground-based instruments and satellites -are urgently needed in the near future.
Data availability.Radiosonde observations used in this work can be acquired from http://weather.uwyo.edu/upperair/sounding.html.Model outputs are stored in the server of School of Atmospheric Sciences at Nanjing University and are available upon request from the corresponding author.

Figure 1 .
Figure 1.Vertical profiles of BC and shortwave heating rate included by unit BC mass in the afternoon (12:00-16:00 LT) during a heavily polluted episode on 23-24 December 2013 in Beijing.The black dashed line denotes the average PBL height during the same daytime period.The BC profile was extracted from regional WRF-Chem modeling by Ding et al. (2016).

Figure 2 .
Figure 2. Diurnal variations of the air temperature change caused by aerosols during haze episodes in Beijing and of PBL height for simulation scenarios with (solid line) and without (dashed line) ARI.

Figure 3 .
Figure 3. (a) Vertical distribution of BC aerosols (contour) and shortwave heating rate at 14:00 LT for runs with (solid line) and without (dashed line) ARI.(b) Vertical potential temperature profile (red) and exchange coefficient profile (blue) for runs with (solid line) and without (dashed line) ARI.The black dashed line represents the averaged PBL height during 12:00-16:00 LT.(c) Diurnal variations of PBL height for runs with (solid line) and without (dashed line) ARI and shortwave radiation attenuation at the surface induced by BC absorption.

Figure 4 .
Figure 4. PBL height (contour map) and turbulence exchange coefficient (grey dashed isoline) variations in percentage as a function of the altitude and mass concentration of BC aerosol layers.Note: each dot on the figure represents an experiment with different BC input.The x axis gives concentrations in the center of the inputted BC plume, while the y axis gives the altitude of the residing plume.The yellow dashed line represents the original PBL height without any BC plume.

Figure 6 .
Figure 6.Variation of absorption extinction coefficient at 400 nm as a function of different mass concentrations of scattering aerosols mixing with BC.The BC concentration is fixed at 5 µg m −3 under various moisture conditions of relative humidity, i.e., 50, 70 and 90 %.Three main cities (BJ, SH and NJ) are marked in the figure with their annual mean PM 2.5 concentrations.

Figure 7 .
Figure 7. Enhanced decreasing/increasing of PBL height (contour map) and enhanced reduction of sensible heat flux at surface (dashed isoline) caused by amplified absorption of BC internally mixed with scattering aerosols.BC concentration is fixed at 5 µg m −3 , and the ratio of scattering aerosol components is SO 2− 4 : NO − 3 : NH + 4 = 3 : 2 : 1, which is based on year-round observations in Beijing (Zhang et al., 2012).The x and y axes are the same as Fig. 4.

Figure 8 .
Figure 8. Net downward shortwave radiation at surface (SWnet_Surf), surface temperature (Temp_Surf), upwelling shortwave radiation at model top (SWupward_Top) and shortwave heating rate, PBL height for urban and rural underlying surfaces without aerosol effect (a-e), and their corresponding perturbations due to aerosol-boundary layer interaction (f-j).

Figure 9 .
Figure 9. Different PBL height variations due to specific BC distribution over (a) urban surface and (b) rural land surface.The x and y axes are the same as Fig. 4.

Table 1 .
WRF-Chem SCM model settings and configuration options.

Table 2 .
Distinctive differences of surface parameters between urban and rural land cover.