Fang Zhang1,2, Yuan Wang3, Jianfei Peng2, Lu Chen1, Yele Sun4, Lian Duan2,5, Xinlei Ge6, Yixin Li7, Jiayun Zhao7, Chao Liu8, Xiaochun Zhang9,10, Gen Zhang9,10, Yuepeng Pan4, Yuesi Wang4, Annie L Zhang11, Yuemeng Ji12, Gehui Wang13, Min Hu14, Mario J Molina15, Renyi Zhang16,7. 1. College of Global Change and Earth System Science, Beijing Normal University, 100875 Beijing, China. 2. Department of Atmospheric Sciences, Texas A&M University, College Station, TX 77843. 3. Division of Geological and Planetary Sciences, California Institute of Technology, Pasadena, CA 91125. 4. State Key Laboratory of Atmospheric Boundary Layer Physics and Atmospheric Chemistry, Institute of Atmospheric Physics, Chinese Academy of Sciences, 100080 Beijing, China. 5. State Environmental Protection Key Lab of Environmental Risk Assessment and Control on Chemical Processes, East China University of Science and Technology, 200237 Shanghai, China. 6. Jiangsu Key Laboratory of Atmospheric Environment Monitoring and Pollution Control, Nanjing University of Information Science and Technology, 210044 Nanjing, China. 7. Department of Chemistry, Texas A&M University, College Station, TX 77843. 8. School of Atmospheric Physics, Nanjing University of Information Science & Technology, 210044 Nanjing, China. 9. State Key Laboratory of Severe Weather, Chinese Meteorological Administration, Chinese Academy of Meteorological Sciences, 100081 Beijing, China. 10. Key Laboratory of Atmospheric Chemistry, Chinese Meteorological Administration, Chinese Academy of Meteorological Sciences, 100081 Beijing, China. 11. Department of Chemistry, College of Natural Sciences, The University of Texas at Austin, Austin, TX 78712. 12. Guangzhou Key Laboratory of Environmental Catalysis and Pollution Control, School of Environmental Science and Engineering, Institute of Environmental Health and Pollution Control, Guangdong University of Technology, 510006 Guangzhou, China. 13. Key Laboratory of Geographic Information Science of the Ministry of Education, School of Geographic Sciences, East China Normal University, 200241 Shanghai, China. 14. State Key Joint Laboratory of Environmental Simulation and Pollution Control, College of Environmental Sciences and Engineering, Peking University, 100871 Beijing, China. 15. Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, CA 92093 mjmolina@ucsd.edu renyi-zhang@tamu.edu. 16. Department of Atmospheric Sciences, Texas A&M University, College Station, TX 77843; mjmolina@ucsd.edu renyi-zhang@tamu.edu.
Abstract
Although regional haze adversely affects human health and possibly counteracts global warming from increasing levels of greenhouse gases, the formation and radiative forcing of regional haze on climate remain uncertain. By combining field measurements, laboratory experiments, and model simulations, we show a remarkable role of black carbon (BC) particles in driving the formation and trend of regional haze. Our analysis of long-term measurements in China indicates declined frequency of heavy haze events along with significantly reduced SO2, but negligibly alleviated haze severity. Also, no improving trend exists for moderate haze events. Our complementary laboratory experiments demonstrate that SO2 oxidation is efficiently catalyzed on BC particles in the presence of NO2 and NH3, even at low SO2 and intermediate relative humidity levels. Inclusion of the BC reaction accounts for about 90-100% and 30-50% of the sulfate production during moderate and heavy haze events, respectively. Calculations using a radiative transfer model and accounting for the sulfate formation on BC yield an invariant radiative forcing of nearly zero W m-2 on the top of the atmosphere throughout haze development, indicating small net climatic cooling/warming but large surface cooling, atmospheric heating, and air stagnation. This BC catalytic chemistry facilitates haze development and explains the observed trends of regional haze in China. Our results imply that reduction of SO2 alone is insufficient in mitigating haze occurrence and highlight the necessity of accurate representation of the BC chemical and radiative properties in predicting the formation and assessing the impacts of regional haze.
Although regional haze adversely affects human health and possibly counteracts global warming from increasing levels of greenhouse gases, the formation and radiative forcing of regional haze on climate remain uncertain. By combining field measurements, laboratory experiments, and model simulations, we show a remarkable role of black carbon (BC) particles in driving the formation and trend of regional haze. Our analysis of long-term measurements in China indicates declined frequency of heavy haze events along with significantly reduced SO2, but negligibly alleviated haze severity. Also, no improving trend exists for moderate haze events. Our complementary laboratory experiments demonstrate that SO2 oxidation is efficiently catalyzed on BC particles in the presence of NO2 and NH3, even at low SO2 and intermediate relative humidity levels. Inclusion of the BC reaction accounts for about 90-100% and 30-50% of the sulfate production during moderate and heavy haze events, respectively. Calculations using a radiative transfer model and accounting for the sulfate formation on BC yield an invariant radiative forcing of nearly zero W m-2 on the top of the atmosphere throughout haze development, indicating small net climatic cooling/warming but large surface cooling, atmospheric heating, and air stagnation. This BC catalytic chemistry facilitates haze development and explains the observed trends of regional haze in China. Our results imply that reduction of SO2 alone is insufficient in mitigating haze occurrence and highlight the necessity of accurate representation of the BC chemical and radiative properties in predicting the formation and assessing the impacts of regional haze.
Because of rapid urbanization, industrialization, and economic growth, many developing countries have experienced frequent haze pollution (1–4), which profoundly impacts human health, weather, and climate (5–9). Noticeably, air quality in the developing world moves along a trajectory that had been previously encountered in most developed countries (1, 4). Air quality is inherently controlled by the synergetic effects of emissions, chemistry, transport, and removal of pollutants (1, 4). While atmospheric circulation is likely subjected to the influence of climate change, fine aerosols also exert direct and indirect feedbacks to climate, via the aerosol–radiation interaction (ARI) and aerosol–cloud interaction, respectively (8–11). Regional haze, which is characterized by high levels of fine aerosols with large spatial and temporal coverage, occurs frequently in Asia, such as China and India (1–4). For example, a month-long haze extreme in January 2013 covered more than 10% (1.3 million km2) of the territory and affected more than 50% (about 800 million) of the population in China (4).Improved understanding of the chemical/physical processes leading to haze formation is crucial to devising effective mitigation strategies to protect the public and reduce uncertainties in climate predictions (5–7). The frequent occurrence of severe haze in China has prompted a variety of legislative actions for air-quality improvement, which have resulted in significantly reduced emissions of SO2 and primary aerosols (1, 4). Satellite top-down estimates showed decreased SO2 emissions from over 30 Mt y−1 in 2005 to less than 10 Mt y−1 in 2016 (4, 12). However, the effectiveness of most regulatory policies has yet to be rigorously evaluated, and implementation of regulatory measures without a scientific basis may also incur undesirable consequences (1, 4). Specifically, it has been suggested that regional haze may be effectively minimized by controlling NH3 emission, but reduction of NH3 without simultaneously regulating emissions of SO2 and NO2 likely poses a threat for the occurrence of highly acidic haze events, analogous to the catastrophic 1952 London Fog (13, 14). Also, while lower levels of air pollutants benefit the public health, most climate models showed that the shield of sunlight by fine aerosols (mainly ARI associated with sulfate, nitrate, and organics) has kept the planet cooler by as much as 0.7 °C globally, suggesting that reduction in aerosol pollution likely accelerates global warming (15). In addition, the black carbon (BC) direct radiative forcing represents another central issue on climate, because the global solar absorption by BC is estimated to be as much as 0.9 W m−2, second only to that by carbon dioxide (16, 17). Currently, large uncertainties exist on BC optical properties, since its ability to absorb light varies considerably during atmospheric processes (18, 19). Under polluted environments, rapid aging of BC particles results in significantly enhanced surface cooling/atmospheric heating and consequently diminished diurnal variation of the atmospheric planetary boundary layer (PBL) and exacerbates haze development (10, 20). In this work, we evaluated the underlying mechanisms for regional haze formation by combining field measurements, laboratory experiments, and model simulations ().
Field and Laboratory Measurements
We analyzed long-term measurements of PM2.5, sulfate, BC, and gaseous concentrations of SO2, NO2, and NH3 in the Northern China Plain (NCP). In our work, air quality was categorized as clean/light haze, moderate haze, and heavy haze conditions, with the corresponding daily mean PM2.5 (particulate matter smaller than 2.5 µm) mass concentrations of <100, 100–200, and >200 µg m−3, respectively (Figs. 1 and 2). Haze occurrence in NCP exhibits an improving trend since 2013 (Fig. 1), as reflected by a decreasing number of heavy haze days and an increasing number of clean/light haze days. For example, the numbers of clean/light haze and heavy haze days are 176 and 58 in 2013 and 227 and 15 in 2018, respectively. On the other hand, there exists no obvious trend for moderate haze events with an average of 113 d from 2004 to 2018 and no improving trends for the daily average mass concentrations of PM2.5 and sulfate during moderate and heavy haze days (Fig. 1 ). The decreasing (increasing) number of heavy (clean/light) haze days correlates to the SO2 trend, which exhibits a constant decrease from 70 parts per billion (ppb) in 2013 to a few ppb in 2018 during heavy haze events (Fig. 1). The declining SO2 level is consistent with other ground-based and satellite measurements (4, 12). In addition, slightly decreasing trends in BC and NO2 during moderate and heavy haze periods are evident in recent years (Fig. 1 ), while no apparent trend exists for NH3 (). Emissions of SO2, NOx, and BC in China are mainly associated with combustion from gasoline or diesel vehicles, residential sources, power plants, and industrial facilities (4). Satellite measurements showed that NH3 emissions from agricultural sources increased substantially over China in recent years, because of increasing temperature and fertilizer use (21).
Fig. 1.
Long-term trends of air quality in Beijing. (A) Trends of the numbers of clean/light haze (blue), moderate haze (red), and heavy haze (black) days since 2004. (B–F) Measurements of PM2.5 (B), sulfate (C), BC (D), SO2 (E), and NO2 (F).
Fig. 2.
An update of air quality from Beijing: no alleviated haze severity. (A) Pictures of the Forbidden City in Beijing during clean/light haze, moderate haze, and heavy haze periods. (B) Measurements of sulfate/BC (Left axis) and PM2.5 (Right axis) mass concentrations during evolutions for clean/light haze (Left), moderate haze (Middle), and heavy haze (Right) episodes in Beijing. (C) Measurements of gaseous SO2/NO2/NH3 (Left axis) and RH (Right axis) during the similar period as in B.
Long-term trends of air quality in Beijing. (A) Trends of the numbers of clean/light haze (blue), moderate haze (red), and heavy haze (black) days since 2004. (B–F) Measurements of PM2.5 (B), sulfate (C), BC (D), SO2 (E), and NO2 (F).An update of air quality from Beijing: no alleviated haze severity. (A) Pictures of the Forbidden City in Beijing during clean/light haze, moderate haze, and heavy haze periods. (B) Measurements of sulfate/BC (Left axis) and PM2.5 (Right axis) mass concentrations during evolutions for clean/light haze (Left), moderate haze (Middle), and heavy haze (Right) episodes in Beijing. (C) Measurements of gaseous SO2/NO2/NH3 (Left axis) and RH (Right axis) during the similar period as in B.The declined frequency but not severity for heavy haze is mostly evident from recent measurements during the 2018–2019 winter. There is a total number of 15 heavy haze events during this period, with the daily mean PM2.5 and sulfate concentrations of 250 and 45 µg m−3, respectively. For comparison, there are on average 29 heavy haze events, with the daily mean PM2.5 and sulfate concentrations of 230 and 40 µg m−3, respectively, during 2004–2018. Heavy haze typically develops with a sharp rise in the sulfate concentration, which occurs concurrently with elevated BC, NO2, NH3, and RH, even at a low SO2 level (Fig. 2 ).To elucidate the chemical mechanism leading to sulfate formation, we performed laboratory experiments by exposing fresh BC particles to variable levels of SO2, NO2, NH3, and RH in a dark environmental chamber (). Monodispersed BC particles with an initial size of 100–150 nm grow rapidly upon exposure to these gaseous species (Fig. 3), as quantified from the simultaneously measured size, mass, and light extinction/scattering (22, 23). Considerable mass growth (defined by dGfm/dt, where Gfm is the mass growth factor) occurs when BC particles are exposed simultaneously to SO2, NO2, and NH3 (Fig. 3 ). The mass growth is insensitive to SO2 in the range of 7–375 ppb (Fig. 3), but depends strongly on the concentrations of NO2 and NH3 (Fig. 3 ). Also, the mass growth increases monotonically with RH in the range of 30–70% and exhibits a slight decreasing trend at high RH (>70%) (Fig. 3). The BC mass growth is also accompanied by changes in the optical properties (Fig. 3). The BC absorption (denoted by absorption enhancement, E) increases significantly on heavily coated (aged) particles. The mass absorption cross-section (MAC) increases from 7.5 to 12.8 m2 g−1 for fresh and aged BC, respectively, with a corresponding E value of 1.7 (Fig. 3). To identify the species responsible for the measured mass growth, we analyzed the particle chemical compositions using a thermal-desorption ion-drift chemical ionization mass spectrometer (TD-ID-CIMS) (24, 25). The contribution of sulfate formation to the particle growth is confirmed by significant presence of several sulfate-related ions, i.e., SO4-2 at m/z = 96 and HSO4− at m/z = 97, on the exposed particles (Fig. 3). In addition, there is negligible growth of BC particles in the absence of one of the three gaseous species (SO2, NO2, or NH3) or at low RH (<10%), suggesting that SO2, NO2, NH3, and H2O are indispensable to sulfate production on BC particles. The sulfate formation on BC is distinct from a previously proposed mechanism that involves aqueous SO2 oxidation by NO2 with NH3 neutralization at high RH (>70%) (13, 14). The BC-catalytic chemistry involves the initial production of nitrous acid (HONO) from the reaction of NO2 on BC and subsequent aqueous reactions with hydrogen sulfite/bisulfite ions to yield sulfate. Uptake of NO2 occurs through physical adsorption at BC surface sites, {S}, and HONO formation occurs via extraction of an allylic hydrogen by NO2 at reducing surface sites, {C-H}red, and to form oxidized carbon, {C}ox, and {HONO} (26),where the bracket denotes adsorbed species or surface sites. NH3 effectively stabilizes {HONO} via the equilibrium with NO2−,The conversion of sulfide to sulfate is subsequently facilitated by {HONO} (27),The rate-limiting step of this mechanism corresponds to HONO production (i.e., Eqs. and ), which explains the strong dependence on NO2, NH3, and RH, but weak dependence on SO2 for the measured sulfate formation (Fig. 3 ). Previous measurements of NO2 uptake on BC surfaces showed that coating of organic acids (glutaric acid) significantly increases the uptake coefficient and HONO yield, because organic acids form strong hydrogen-bonded complexes with polar organic molecules from the occupying {C–H}red sites (26). Also, sulfuric acid coating was shown not to alter the uptake NO2 coefficient but to reduce the amount of HONO released (26). In contrast, our work reveals an insignificant effect of sulfate formation on BC particle growth upon exposure to SO2, NO2, and NH3. Since our laboratory experiments are conducted to mimic ambient conditions, the mechanism of BC-catalyzed sulfate formation is applicable to other polluted regions worldwide.
Fig. 3.
BC-catalyzed sulfate formation. (A–D) Measured mass growth rate, dGfm/dt), as a function of SO2, NO2, NH3, and RH. (E) Absorption enhancement (Eabs) and MAC for fresh and aged BC particles. (F) Single-ion monitoring of collected particles for the sulfate ion SO4− at m/z = 96, the bisulfate ion HSO4− at m/z = 97, and the oxygen adduct of SO3O2− at m/z = 112 measured by ID-TD-CIMS. BC particles during the experiments have an initial diameter (D) of 100–150 nm.
BC-catalyzed sulfate formation. (A–D) Measured mass growth rate, dGfm/dt), as a function of SO2, NO2, NH3, and RH. (E) Absorption enhancement (Eabs) and MAC for fresh and aged BC particles. (F) Single-ion monitoring of collected particles for the sulfate ion SO4− at m/z = 96, the bisulfate ion HSO4− at m/z = 97, and the oxygen adduct of SO3O2− at m/z = 112 measured by ID-TD-CIMS. BC particles during the experiments have an initial diameter (D) of 100–150 nm.Using the laboratory-derived sulfate growth rate on BC particles () and the field-measured concentrations of gaseous species and BC (), we quantified sulfate production during pollution episodes in Beijing (Fig. 4). During moderate haze events, the measured sulfate mass concentration by an aerosol mass spectrometer (AMS) (28) increases steadily and reaches a peak mass concentration of about 15 µg m−3. Incorporation of the BC catalytic SO2 reaction reproduces the measured sulfate mass concentration (Fig. 4 and ). During heavy haze events, the mass concentration of sulfate increases rapidly to about 60–100 µg m−3 (Fig. 4 and ), while the calculated BC catalytic sulfate formation (17–48 µg m−3) accounts for less than 50% of the measured sulfate production. In our work, fresh BC mass concentrations were estimated from a previously measured ratio of fresh to aged BC particles in Beijing (29), to derive the range of sulfate production rates (Fig. 4 and ). To further gain insight into the role of BC particles in sulfate production, we determined the number and mass fractions of BC-containing aerosols using a single-particle AMS (SP-AMS) and a combined SP2 (single-particle soot photometer) and AMS (), respectively. During moderate haze events, the number and mass fractions of BC-containing particles are 76% and 78%, respectively, showing significant BC presence (Fig. 4). Using the measured mass fraction of BC-containing particles from SP2-AMS, along with the PM2.5 mass concentration and the mass fraction of BC particles, we derived a sulfate concentration of about 12 µg m−3, consistent with that measured by AMS (Fig. 4). During heavy haze events, the number and mass fraction of BC-containing particles are 64% and 44% (Fig. 4), respectively. The derived sulfate concentration of 13 µg m−3 is much smaller than that measured by AMS (Fig. 4). Hence, the BC catalytic sulfate production accounts dominantly (about 90–100%) and nonnegligibly (about 30–50%) for sulfate production during the moderate and heavy haze events, respectively. The coexistence of BC, sulfate, nitrate, and organics (Fig. 4 ) on PM2.5 also implies a mutually promoting mechanism among the various aerosol ingredients to produce high levels of PM2.5 (4, 13). In particular, formation of hygroscopic sulfate is essential to drive multiphase chemistry of other aerosol species (1).
Fig. 4.
Quantifying BC-catalyzed sulfate formation. (A and B) Calculated (green lines) and measured (circles and dashed line) sulfate mass concentrations during a moderate day and a heavy haze day, respectively. (C and D) BC mass concentration (Left axis) and sulfate formation rate (d[SO42-]/dt, Right axis) during a moderate haze day and a heavy haze day, respectively. The ranges of sulfate mass concentration (green lines in A) and sulfate formation rate (blue line in B) are derived by assuming that 20% and 40% of the measured BC mass concentration are freshly emitted (34). The peak PM2.5 and sulfate mass concentrations for the two episodes (27–30 November 2018 and 1–3 December 2018) are 154.0 and 14.2 µg m−3 and 330.0 and 55.8 µg m−3, respectively. (E and F) Number and mass fractions of BC-containing and non-BC particles in ambient PM2.5 during moderate and heavy haze events, respectively. The mass fraction of BC-containing particles are derived according to our field-measured BC mass concentrations and previous measurements of mass fraction of each species in BC-containing particles (represented by the pie charts) by a combined single-particle SP2 and AMS during typical moderate and heavy haze days in urban Beijing (37). The number concentration is measured from February 1 to March 11, 2019, during which the moderate and heavy haze events correspond to PM2.5 mass concentrations of 160.2 and 355.0 µg m−3 and BC mass concentrations of 10.4 and 17.8 µg m−3, respectively.
Quantifying BC-catalyzed sulfate formation. (A and B) Calculated (green lines) and measured (circles and dashed line) sulfate mass concentrations during a moderate day and a heavy haze day, respectively. (C and D) BC mass concentration (Left axis) and sulfate formation rate (d[SO42-]/dt, Right axis) during a moderate haze day and a heavy haze day, respectively. The ranges of sulfate mass concentration (green lines in A) and sulfate formation rate (blue line in B) are derived by assuming that 20% and 40% of the measured BC mass concentration are freshly emitted (34). The peak PM2.5 and sulfate mass concentrations for the two episodes (27–30 November 2018 and 1–3 December 2018) are 154.0 and 14.2 µg m−3 and 330.0 and 55.8 µg m−3, respectively. (E and F) Number and mass fractions of BC-containing and non-BC particles in ambient PM2.5 during moderate and heavy haze events, respectively. The mass fraction of BC-containing particles are derived according to our field-measured BC mass concentrations and previous measurements of mass fraction of each species in BC-containing particles (represented by the pie charts) by a combined single-particle SP2 and AMS during typical moderate and heavy haze days in urban Beijing (37). The number concentration is measured from February 1 to March 11, 2019, during which the moderate and heavy haze events correspond to PM2.5 mass concentrations of 160.2 and 355.0 µg m−3 and BC mass concentrations of 10.4 and 17.8 µg m−3, respectively.
Radiative Transfer Calculations
Using a radiative transfer model (30), we quantified the variations of the direct radiative forcing (ΔDRF) with and without inclusion of this BC catalytic chemistry on the top of atmosphere (TOA), in the atmosphere (ATM), and at the surface (SFC) (Fig. 5 and ). The model simulations show nearly invariant TOA ΔDRF values of 0.3, 2.1, and 3.4 W m−2 from clean/light haze, moderate haze, to heavy haze days, respectively, indicating that the positive ATM ΔDRF is nearly offset by negative SFC ΔDRF. On the other hand, the BC-catalyzed sulfate formation causes large heating in ATM and large cooling at SFC. The ΔDRF values range from 38.3, 126.8, to 175.9 W m−2 and from −38.0, −124.8, to −172.6 W m−2, respectively, from clean/light, moderate, to heavy haze days. The large increases in ATM heating and SFC cooling during the haze evolution are attributed to increased coating by sulfate (Figs. 4 and 5), which substantially increases MAC or E for BC particles (Fig. 3). The large ATM and SFC ΔDRF values also indicate significant atmospheric stabilization, which exacerbates accumulation of PM2.5 within the PBL (Fig. 5) (4, 10). Air stagnation is clearly reflected from the measured PBL height, which decreases from 1 to 2 km on clean/light haze days to below 0.5 km on heavy haze days (Fig. 5 and ). The BC-induced air stagnation imposes strong feedbacks to haze formation and evolution (1, 10). Note that the ΔDRF results in our work likely represent a lower bound, since the radiative forcing is expected to be enhanced by coating of other non-BC materials, i.e., nitrate and organics, in addition to sulfate (16). Clearly, the topic of aging and variations in the optical properties and radiative forcing associated with the BC-catalyzed sulfate formation warrants further investigation.
Fig. 5.
BC-catalyzed chemistry on regional haze. (A) Calculated (∆DRF) of aged BC particles at the TOA, yellow), in the atmosphere (ATM, red), and at the surface (SFC, gray) during clean/light haze, moderate haze, and heavy haze events. The initial fresh BC diameter is ∼150 nm. The colored boxes denote the vertical extends of the free troposphere (dark blue) and the PBL (light blue, gray, or dark gray) during haze evolution from clean/light haze, moderate haze, and heavy haze events, respectively. (B) Schematic representation of the morphology evolution from fresh to aged BC particles because of the BC-catalyzed reaction, which significantly impacts light absorption. (C) Evolution in the PBL height, absorption enhancement (Eabs), and PM2.5 mass concentration driven by the BC-catalyzed and aqueous sulfate formation (13, 14). The green, red, and black colors correspond to PBL height, Eabs, and PM2.5 mass concentration, respectively.
BC-catalyzed chemistry on regional haze. (A) Calculated (∆DRF) of aged BC particles at the TOA, yellow), in the atmosphere (ATM, red), and at the surface (SFC, gray) during clean/light haze, moderate haze, and heavy haze events. The initial fresh BC diameter is ∼150 nm. The colored boxes denote the vertical extends of the free troposphere (dark blue) and the PBL (light blue, gray, or dark gray) during haze evolution from clean/light haze, moderate haze, and heavy haze events, respectively. (B) Schematic representation of the morphology evolution from fresh to aged BC particles because of the BC-catalyzed reaction, which significantly impacts light absorption. (C) Evolution in the PBL height, absorption enhancement (Eabs), and PM2.5 mass concentration driven by the BC-catalyzed and aqueous sulfate formation (13, 14). The green, red, and black colors correspond to PBL height, Eabs, and PM2.5 mass concentration, respectively.
Discussion
Our analysis of long-term measurements of gases and aerosols in NCP shows an improved air-quality trend since 2013, as evident from increasing days of clean/light haze and decreasing days of heavy haze, along with a significantly decreased SO2 (Fig. 1). The latter is consistent with implementation of strict air-quality standards in China, particularly associated with reduction in coal burning from power plants, industrial facilities, and winter residential heating (4). However, despite considerably lower levels of SO2 and reduced (increased) frequency of heavy haze (clean/light haze), heavy haze events with high levels of daily PM2.5 and sulfate content persist in recent years, showing little alleviated severity (Fig. 2 ). In addition, moderate haze occurs annually on about 30% of the days between 2004 and 2018, with no reduced frequency or sulfate level.Haze development in China is characterized by rapid sulfate formation (Fig. 2). SO2 oxidation occurs through gas-phase oxidation or aqueous reactions, but the detailed chemical mechanisms remain controversial (1, 3, 4). The gas-phase oxidation of SO2 is slow, with a corresponding lifetime of about 1 wk at typical tropospheric levels of hydroxyl radicals (1). For example, the relatively long lifetime of SO2 was invoked as an explanation for significantly reduced SO2 emission but without a commensurate reduction in aerosol pollution in the eastern United States (31). Also, atmospheric measurements revealed rapid sulfate production during severe haze events in China (1, 3, 4, 13), which cannot be explained by current atmospheric models (4) and suggests missing sulfur oxidation mechanisms (13).Our work offers explanations for the puzzling trends of regional haze from the long-term measurements in China (Fig. 5). Haze formation in NCP involves two distinct secondary processes, i.e., new particle formation and growth (1, 2, 13). During haze evolution from clean to polluted conditions, photochemistry is considerably inhibited because of reduced UV radiation by aerosols (13, 32, 33). While new particle formation is mainly driven by photochemistry and occurs consistently under clean conditions (1, 2), rapid particle growth under polluted conditions requires efficient multiphase chemistry (1, 2, 14). Our work unravels a sulfate formation mechanism: SO2 oxidation is efficiently catalyzed by BC in the presence of NO2 and NH3, even at low SO2 levels (down to a few ppb) and an intermediate RH range (30–70%). This BC catalytic pathway accounts for about 90–100% and 30–50% of sulfate formation during moderate and heavy haze events in Beijing, respectively. Sulfate formation enhances BC aging and light absorption, leading to atmospheric stabilization (Fig. 5 ). Also, BC-containing particles likely engage in coagulation growth with nucleation-mode particles, which consist dominantly of secondary organics (1, 2, 34), to further contribute to BC aging and light absorption (16). On the other hand, our results indicate that this BC catalytic reaction yields only a modest sulfate amount. Other mechanisms, such as the aqueous SO2 oxidation by NO2 with NH3 neutralization (13, 14, 35), are necessary to account for large sulfate formation during heavy haze events. The BC-catalyzed sulfate formation represents the most critical step to increase aerosol hygroscopicity and aqueous chemistry for secondary aerosol formation, i.e., sulfate, nitrate, and organics (1, 4, 13, 14). Furthermore, air stagnation caused by BC light absorption traps moisture within PBL (4), as evident from considerably increased RH, i.e., from <30%, 30–60%, to 60–90% during the evolution from clean/light haze, moderate haze, to heavy haze conditions, respectively (Fig. 2). Elevated RH results in a positive feedback to multiphase chemistry to further enhance BC-catalyzed and aqueous SO2 oxidation (4).The current air-pollution status in China corresponds to a significantly reduced SO2 level, yet still elevated BC, NOx, and NH3 levels (Fig. 2 ). Under such a circumstance, BC-catalyzed SO2 oxidation occurs efficiently and is responsible for no improved trend for moderate haze events in China. Reduced levels of primary gases and aerosols, as reflected by significantly declined SO2 levels and slightly declined NOx and BC levels (Fig. 1), from emission control measures in China inevitably prevent achieving the optimal conditions (i.e., NOx, NH3, and BC levels as well as RH and air stagnation) for efficient aqueous chemistry to produce high levels of PM2.5 and sulfate (13, 14). This likely explains the decreasing frequency of heavy haze events and increasing frequency of clean/light haze days. However, under unfavorable meteorological conditions, the BC-catalyzed sulfate pathway facilitates the transition to aqueous chemistry to produce high levels of PM2.5 during heavy haze events, explaining the long-term measurements of no alleviated haze severity in China. Furthermore, the BC-catalyzed sulfate formation results in an invariant TOA radiative forcing of nearly 0 W m−2 throughout haze development, leading to little net climatic cooling or warming by regional haze (Fig. 5).A most recent work indicated significantly a declined annual mean PM2.5 concentration in China over the last 5 y because of anthropogenic emission abatements (36), consistent with our analysis of decreasing heavy haze days and increasing clean/light haze days. On the other hand, our results reveal no alleviated haze severity and little improved trend for moderate haze events, despite lower SO2 levels and reduced frequency of heavy haze events. We conclude that the BC-catalyzed chemistry plays a dominant role in driving the formation, trend, and radiative forcing of regional haze, highlighting the necessity of simultaneous reduction in emissions of SO2, NOx, NH3, and BC for air-quality improvement.
Methods
Field Measurements.
Field measurements were performed to measure the nonrefractory PM2.5/chemical components including sulfate (SO42-), nitrate (NO3−), ammonium (NH4+), chloride (Cl), and organics (Org) by a PM2.5-Q-ACSM during 11 November to 20 December 2018 in Beijing and Gucheng, which is a suburban site about 120 km south of Beijing. At the both sites, a seven-wavelength aethalometer (model AE33, Magee Scientific Corp.) was used to measure BC. BC Measurements using different methods typically yielded variable results, because of lack in standards (37). For example, measurements by the aethalometer could suffer from interference from non-BC coatings (16). Gaseous species of NO2 and SO2 were measured by gas analyzers (ECOTECH). In Beijing, gaseous NH3 was measured using a cavity ring-down spectrometer (Picarro G2103). The number fractions of BC-containing particles were measured during February 1 to March 11, 2019 by an SP-AMS (Hexin, Inc.) (38). The mass fraction of BC-containing particles was derived according to field measurements of chemical composition of BC-containing particles by a combined single-particle SP2 and AMS (39). The long-term sulfate mass concentration in submicrometer particles during 2011–2018 was measured at the Institute of Atmospheric Physics, Chinese Academy of Sciences in Beijing (40). The concentrations of SO2, NO2, and PM2.5 during 2009–2012 were measured on the campus of Peking University (2), and the data during 2013–2018 were measured at a site (Wanliu) close to the campus of Peking University (http://beijingair.sinaapp.com/). The NH3 measurements during winter 2016–2018 were collected using a NH3 analyzer (Ecotech 9842, Australia) on the campus of China Meteorological Administration. Long-term trends of BC and NH3 were taken from previously published measurements (41–43). The PBL height during 11–12 January of 2019 in urban Beijing was retrieved from Raman Lidar measurements (44).
Laboratory Experiments.
The laboratory experiments were conducted by exposing seed particles to SO2, NO2, and NH3 at variable RH and measuring the dry size variation and sulfate formation on the exposed particles in a 1.2-m3 Teflon reaction chamber covered with aluminum foil as described previously (13, 14). BC particles from incomplete combustion of propane were heated by a thermal denuder to 300 °C to remove coated organics and then were continuously sampled, diluted, dried, and size-classified by a differential mobility analyzer (DMA) (3081, TSI). The BC density generated from propane combustion in our laboratory experiments was 0.5–0.6 g cm−3, which is within the range of ambient measured density for fresh BC particles in Beijing (16, 18, 22, 30) and representative of fresh BC in the atmosphere. An integrated instrument including the DMA, an aerosol particle mass analyzer (APM, 3600, Kanomax), a condensation particle counter (3760A, TSI), and a thermal denuder was used for measurements of the particle size and mass. A commercial Nephelometer (TSI, 3563) and a cavity ring-down spectrometer was used to measure optical properties, as described elsewhere (45). The chemical composition of the exposed particles was measured by a TD-ID-CIMS (13, 24, 25). The particle mass growth factors (Gfm) are expressed as mp/m0, respectively, where m0 and mp are the particle mass before and after aging. The material density of particles is calculated using those for BC and sulfate (ρBC = 1.77 g cm−3, ρsul = 1.77 g cm−3). The particle absorption (Eabs) cross-sections were derived by dividing the corresponding optical coefficient by particle number concentration. The TD-ID-CIMS was used to measure sulfate formation on BC. The mp and m0 were determined from the size of BC particles measured by DMA and APM, respectively. When determining sulfate formation on BC particles, the bulk measured BC mass concentration by the Aethalometer was used by assuming size-independent mass growth on BC particles due to sulfate formation according to our laboratory results, which showed similar growth for 100- and 150-nm BC particles.
Parameterization of BC-Catalyzed Sulfate Formation from Laboratory Experiments.
Based on the laboratory results, we derived an exponential parameterization for the mass growth rate, denoted as dGfm/dt, which was used to calculate ambient sulfate production rate from measured trace gases concentrations and RH as follows:where [SO2], [NO2], and [NH3] are the concentrations of gaseous SO2, NO2, and NH3 respectively. We employed a polynomial function f(RH) to parameterize the RH dependence of the laboratory results shown in Fig. 2,We derived the power exponential constants of α, β, γ, and δ using the power exponential relationship for each variable according to dGfm/dt (dGfm/dt = f(X) = a[X]), where a is the coefficient and [X] is the concentration of the gas precursor, with n = α, β, γ, or δ. A is a prefactor constant derived from the linear fit of dGfm/dt and the product of the power exponential for each variable. An iteration method was further applied to locate the residual value, which is the average of the differences of experimentally measured dGfm/dt and the fitted dGfm/dt calculated from A, . We obtained the following relationship between the gaseous concentrations and RH (>30%):The parameterization to estimate the growth rate was evaluated by comparing with the experimental measurements, showing the growth rate with a slope and correlation coefficient of 1.03 and 0.87, respectively, between the predicted and experimental values (). Since the concentrations of 7–750, 37–750, and 100–1,000 ppb for the gaseous SO2, NO2 and NH3, respectively, in our laboratory experiments are close to the ambient levels, this parameterization is applicable broadly to atmospheric models. On combining the observed BC mass concentration and the laboratory results, the rate of sulfate formation (d [SO42-]/d t in μg m−3⋅h−1) and sulfate mass concentration are calculated byOnly the episodes with stable meteorological conditions were selected for the estimation. In this study, we assumed that 20% and 40% of the observed BC particles are freshly emitted, to derive the upper and lower limits, respectively (29).
Calculations of Radiative Forcing.
The Mie theory with the core–shell assumption for homogeneous spherical particles was used to calculate MAC (46). Because the core−shell Mie model assumes that all particles are compact and both the BC core and outer shell are homogeneous, the core of aged BC particles is represented as the mass-equivalent diameter (Dme) of fresh BC particle (44). For aged BC particles, the size distribution of the core (Dcore) is equal to the size distribution of the mass-equivalent diameter of fresh BC particle (Dme,0). With the assumption that the increases of mass-equivalent diameter (ΔDme) are identical for BC particles within the entire size range, the mass-equivalent diameter of aged BC particles is calculated by adding Dme,0 and ΔDme,0. With an insignificant organics fraction for fresh BC particles (<30%) during the experiments, we adopted a refractive index of 1.75 + 0.44i for fresh BC particles to determine E. Since the analysis of particle composition by TD-ID-CIMS indicates the presence of dominant sulfate as the coating materials in our experiments, the refractive index of coating sulfate used in this study is 1.60–0.13i (47). The ΔDRF for TOA, ATM, and SFC due to BC-catalyzed sulfate formation is calculated and evaluated by using the Santa Barbara DISORT Atmospheric Radiative Transfer model (30).
Data and Materials Availability.
All data needed to evaluate the conclusions in the paper are present in the paper and/or the . Additional data related to this paper may be requested from the authors.
Authors: Guoyao Wu; Jacob Brown; Misti L Zamora; Alyssa Miller; M Carey Satterfield; Cynthia J Meininger; Chelsie B Steinhauser; Gregory A Johnson; Robert C Burghardt; Fuller W Bazer; Yixin Li; Natalie M Johnson; Mario J Molina; Renyi Zhang Journal: Proc Natl Acad Sci U S A Date: 2019-05-28 Impact factor: 11.205
Authors: Song Guo; Min Hu; Misti L Zamora; Jianfei Peng; Dongjie Shang; Jing Zheng; Zhuofei Du; Zhijun Wu; Min Shao; Limin Zeng; Mario J Molina; Renyi Zhang Journal: Proc Natl Acad Sci U S A Date: 2014-11-24 Impact factor: 11.205
Authors: Kristal A Rychlik; Jeremiah R Secrest; Carmen Lau; Jairus Pulczinski; Misti L Zamora; Jeann Leal; Rebecca Langley; Louise G Myatt; Muppala Raju; Richard C-A Chang; Yixin Li; Michael C Golding; Aline Rodrigues-Hoffmann; Mario J Molina; Renyi Zhang; Natalie M Johnson Journal: Proc Natl Acad Sci U S A Date: 2019-02-11 Impact factor: 11.205
Authors: Lei Yao; Xiaolong Fan; Chao Yan; Theo Kurtén; Kaspar R Daellenbach; Chang Li; Yonghong Wang; Yishuo Guo; Lubna Dada; Matti P Rissanen; Jing Cai; Yee Jun Tham; Qiaozhi Zha; Shaojun Zhang; Wei Du; Miao Yu; Feixue Zheng; Ying Zhou; Jenni Kontkanen; Tommy Chan; Jiali Shen; Joni T Kujansuu; Juha Kangasluoma; Jingkun Jiang; Lin Wang; Douglas R Worsnop; Tuukka Petäjä; Veli-Matti Kerminen; Yongchun Liu; Biwu Chu; Hong He; Markku Kulmala; Federico Bianchi Journal: Environ Sci Technol Lett Date: 2020-09-25