Subscriber access provided by UNIV OF CAMBRIDGE
Characterization of Natural and Affected Environments
Size-related physical properties of black carbon in the lower atmosphere over Beijing and Europe Shuo Ding, Dantong Liu, Delong Zhao, Kang Hu, Ping Tian, Wei Zhou, Mengyu Huang, Yan Yang, Fei Wang, Jiujiang Sheng, Quan Liu, Shaofei Kong, Pengyi Cui, Yuandong Huang, Hui He, Hugh Coe, and Deping Ding Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.9b03722 • Publication Date (Web): 19 Aug 2019 Downloaded from pubs.acs.org on August 19, 2019
Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.
is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.
Page 1 of 28
Environmental Science & Technology
1
Size-related physical properties of black carbon in the
2
lower atmosphere over Beijing and Europe
3 4
Shuo Ding1,Dantong Liu1,*, Delong Zhao2,3,4,6*, Kang Hu1, Ping Tian2, Wei Zhou2,
5
Mengyu Huang2, Yan Yang2, Fei Wang2, Jiujiang Sheng2, Quan Liu2, Shaofei
6
Kong5,Pengyi Cui7, Yuandong Huang7, Hui He2, Hugh Coe8, Deping Ding2,3,4
7 8
1Department
of Atmospheric Sciences, School of Earth Sciences, Zhejiang University, Hangzhou, China
9
2Beijing
10 11
3Beijing
Key Laboratory of Cloud, Precipitation and Atmospheric Water Resources, Beijing, China
12 13
4Field
experiment base of cloud and precipitation research in North China, China Meteorological Administration, Beijing, China
14 15
5Department
of Atmospheric Sciences, School of Environmental Studies, China University of Geosciences (Wuhan), Wuhan, China
16
6Nanjing
17 18
7School
of Environment and Architecture, University of Shanghai for Science and Technology, Shanghai, China
19 20
8Centre
Weather Modification Office, Beijing, China
University, Nanjing, China
for Atmospheric Sciences, School of Earth and Environmental Sciences, University of Manchester, Manchester M13 9PL, UK.
21 22 23
Corresponding to:
[email protected],
[email protected] ACS Paragon Plus Environment
Environmental Science & Technology
24
Abstract
25
The size-resolved properties of atmospheric black carbon (BC) importantly determine
26
its absorption capacity and cloud condensation nuclei (CCN) ability. This study reports
27
comprehensive vertical profiles of BC size-related properties over the Beijing area (BJ)
28
and
Continental
29
Europe
(CE).
30
mass loadings over
31
CE were in the range
32
of clean background
33
over BJ. For both
34
planetary
35
layer (PBL) and lower
36
free troposphere, BC mass median core diameter over BJ during the cold season was
37
0.21±0.02μm, larger than warm season over BJ and CE (0.18±0.01μm), which may
38
reflect seasonal differences in emissions. The BC coatings were positively correlated
39
with pollution level, with background BC having smaller coated count median diameter
40
(0.19 ± 0.01μm). The modelled absorption enhancement (Eabs) due to coatings was
41
1.23±0.14 for background but in the PBL following a linear expression
42
(Eabs=0.13*MassBC,surface+1.26).The CCN ability of BC was significantly enhanced in
43
the polluted PBL, due to both enlarged size and increased hygroscopicity. In polluted
44
BJ at predicted supersaturations~0.08% half of the BC number could be activated,
45
whereas in cleaner environment needs ~0.14%. The results here suggest the highly
46
coated and absorbing BC can be efficiently incorporated into clouds and exert important
47
indirect radiative impacts over polluted East Asia region.
BC
boundary TOC art
48 49
ACS Paragon Plus Environment
Page 2 of 28
Page 3 of 28
Environmental Science & Technology
50
Introduction
51
Black carbon (BC) dominates the absorption of shortwave solar radiation in
52
atmospheric particulate matter. This alters the atmospheric thermodynamics by its
53
heating effects, especially in geographical hotspots that experience large emissions
54
from anthropogenic or biomass burning sources1. Fresh BC particles are hydrophobic
55
but can transform to be hydrophilic after aging thereby becoming cloud condensation
56
nuclei (CCN)2-4. Previous studies reported that the BC layer present within clouds could
57
cause important semi-direct effects
58
incorporated into cloud droplets8.
59
5-7
or lead to evaporation of cloud water if it was
Both the size of refractory BC (rBC) core and the size including non-absorbing
60
coatings
61
impacts of BC in the atmosphere. The combination of both terms depends on the
62
complex morphology and mixing state of BC. These have previously been
63
mathematically treated as core-shell 11, partly encapsulated 12or externally mixed 13, all
64
of which produce different results and each may provide reasonable approximations in
65
different situations. The weighting between rBC core and coated size could vary at
66
different stages of ageing process. The BC core is the essential part to determine its
67
absorbing properties, and it could be sometimes used to determine the BC source
68
signature 9, 14, 15. The amount of coating material associated with a BC particle and how
69
the coating is associated determine how efficiently the incident light interacts with the
70
absorbing core 16, 17. The addition of coatings can alter the initial hygroscopicity of BC18
71
and also increase the particle size. The evolution of coated BC size could thus increase
8-10,
which are the uncoated and coated size, determine the properties and
ACS Paragon Plus Environment
Environmental Science & Technology
Page 4 of 28
72
its ability of acting as cloud condensation nuclei, thereby determining its atmospheric
73
lifetime 19.
74
Although the mass concentration of BC has been widely measured globally20, there
75
is still lack of measurements of the BC size distribution in both horizontal and vertical
76
directions. The uncertainties in modelling BC mass loadings could largely result from
77
the unresolved size distributions of BC 21. To incorporate the mixing state information
78
of BC at different sizes into models could improve the realistic treatments of BC in
79
terms of radiative effects 22 and CCN-activity 23. In addition, it is important to study the
80
vertical distribution of BC properties because the location of BC in the atmospheric
81
column and the vertical gradient of BC properties could importantly determine its
82
impacts on the thermodynamic balance of the planetary boundary layer (PBL)
83
Depending on the location of BC relative to cloud layers, the BC could either promote
84
or inhibit cloud development 5.
85
The aerosol properties over polluted megacity have been raised much attention, where
86
their vertical profiles such as over the North China Plain were measured using a range
87
of platforms
88
properties over Beijing and compares these with those above continental Europe
89
throughout the boundary layer and lower free troposphere. The results are used to
90
examine the size distribution and size-related optical and hygroscopic properties of BC
91
in detail.
25, 26.
24.
This study characterizes the vertical distributions of BC physical
ACS Paragon Plus Environment
Page 5 of 28
Environmental Science & Technology
92
Materials and methods
93
Aircraft platform
94
The flights conducted over Beijing region used a King-Air350 aircraft with an
95
average speed about 250 km/hr27, 28. The flights were conducted in 2018 during both
96
warm (April to September) and cold (November to December) seasons. The sampling
97
inlet system is a passive Isokinetic Aerosol Sampling Inlet (BMI, Brechtel
98
Manufacturing Inc), which can deliver 150 lpm of sample flow at 100 m s-1 air speed
99
and transmit particles with collection efficiency over 95% over a particle size range of
100
0.01-6 μm. A silica drier was used in the sampling line before the instruments. This,
101
combined with the constant cabin temperature (~25℃), guarantees the sampled aerosol
102
was dry. Most flights were performed to avoid clouds but data collected when the
103
aircraft was in cloud has been eliminated based on the measurements of cloud liquid
104
water content. Most flights were conducted up to 3 km altitude. The flight tracks are
105
shown in Fig. 1(a), and the aircraft was typically operated over the Shahe rural area to
106
the northeast of central Beijing away which was ~20 km. The flights over continental
107
Europe took place between April and September 2008 during the EUCAARI-
108
LONGREX (European integrated Project on Aerosol Cloud Climate and Air Quality
109
Interactions - Long Range Experiment) campaign
110
campaign, a high pressure system continuously impacted western Europe and hence
111
regionally-influenced urban pollutants were sampled.
29.During
112
ACS Paragon Plus Environment
the EUCAARI aircraft
Environmental Science & Technology
113
Meteorology
114
Meteorological parameters including ambient pressure, temperature, relative
115
humidity and wind speed/direction were in-situ characterized by an AIMMS-20
116
(Aventech Research Inc), which was calibrated on an annual basis. The planetary
117
boundary layer height (PBLH) is determined by a combination of factors (illustrated in
118
Fig. S1): a low vertical gradient of potential temperature indicating the effective
119
convective mixing and a temperate inversion on top of the PBL.
120
Backward trajectory analysis was performed by HYSPLIT 4.0 model30 in ensemble
121
mode using 1◦×1◦, 3-hourly GDAS1 reanalysis products. A 24h backward time is
122
analyzed since the focus is on recent potential source influences.
123 124
Measurements on black carbon physical properties
125
The mass of the refractory core (rBC) of an individual BC particle is measured by a
126
single particle soot photometer (SP2, DMT Inc.)31, 32. The SP2 was calibrated using
127
Aquadag (Acheson Inc., USA) and a factor of 0.75 was then applied to correct the
128
calibration so that it is representative of ambient BC 33. The mass equivalent diameter
129
of the rBC core (Dc) is obtained by assuming sphericity and a rBC density of 1.8 g cm-3
130
34. The BC mass outside of the detectable range is obtained by extrapolating a lognormal
131
distribution on the measured mass size distribution of the BC cores (Fig. 3(b)). The
132
vertical profiles in this study are classified as low, medium and heavy pollution
133
(denoted as LP, MP, HP respectively) according to the occurrence frequency of surface
134
rBC mass loading, i.e. 0-0.4 μg m-3as LP, 0.4-1 μg m-3 as MP and >1 μg m-3 as HP.
ACS Paragon Plus Environment
Page 6 of 28
Page 7 of 28
Environmental Science & Technology
135
The coated BC size (Dp) is obtained by applying a Mie-look up table to match the
136
modelled scattering cross section to the measured values9, 35. For a given time window,
137
the mass (or count) median diameter (MMD or CMD) from an uncoated or coated BC
138
size distribution is obtained as above and below which the mass (or number) population
139
is equal.
140
Fig. 2 shows an example of the distribution of the diameters of the coated BC
141
particles as a function of uncoated BC core size of single particles for the
142
flight20180426 in the PBL. Dp=Dc means there is no coating associated with BC, and
143
the coating increases as Dp becomes larger than Dc. The bulk coating/rBC volume ratio
144
(Vcoating/VrBC) is calculated as the total volume of coated BC divided by the total volume
145
of uncoated rBC cores in a given time window, expressed by
146
𝑉𝑐𝑜𝑎𝑡𝑖𝑛𝑔 𝑉𝑟𝐵𝐶
∑𝑖𝐷3𝑝,𝑖
= ∑ 𝐷3 ―1(1), 𝑖 𝑐,𝑖
147
where Dp,i and Dc,iare the coated and uncoated diameters of ith BC particle respectively.
148
In addition, the mass ratio of coating and rBC (MR) could also be obtained for each
149
particle. Note that there is a retrieval rate of successfully SP2-measured Dp and this
150
uncertainty is estimated by assuming two extreme scenarios for the particles missing
151
Dp information (bottom panel of Fig. 2), as detailed in the Supplement. This uncertainty
152
analysis is performed throughout the following calculations of optical and hygroscopic
153
properties of BC.
154
If the volume ratio of coating and rBC is known, the hygroscopicity parameter of a
155
given BC-containing particle (κBCc) could be obtained by assuming a κcoating and κrBC=0,
156
based on volume-weighted Zdanovskii–Stokes–Robinson (ZSR) rule36, expressed as,
ACS Paragon Plus Environment
Environmental Science & Technology
157
𝜅𝐵𝐶𝑐 = 𝑉𝑐𝑜𝑎𝑡𝑖𝑛𝑔 × 𝜅𝑐𝑜𝑎𝑡𝑖𝑛𝑔(2),
158
The direct measurement on coating composition was not available in this study hence
159
κcoating is not directly determined, thus we applied a value of κcoating=0.3 measured from
160
a previous study37 and a reduced value of 0.2 for a sensitivity test. We only tested the
161
sensitivity for κcoating3 is considered to have a core-shell structure leading to an absorption
166
enhancement caused by a lensing effect. The calculated MAC of each particle is
167
weighted by the total rBC mass to obtain the MAC in bulk, expressed as:
168
𝑀𝐴𝐶 =
169
where 𝑀𝐴𝐶𝑖 and 𝑚𝑟𝐵𝐶,𝑖 are the MAC and rBC mass for each particle respectively.
∑𝑖𝑀𝐴𝐶𝑖 × 𝑚𝑟𝐵𝐶,𝑖 ∑𝑖𝑚𝑟𝐵𝐶,𝑖
(3) ,
170 171
Results and discussion
172
rBC mass loading
173
Fig. 1(c) shows the vertical profiles of rBC mass loading for all continental Europe
174
and Beijing flights under different pollution levels. The rBC mass loading over
175
continental Europe was generally of the same order of magnitude as that in light
176
pollution levels over Beijing, with mean rBC mass loadings in the PBL and FT