Size-related physical properties of black carbon in the lower

5 days ago - In polluted BJ at predicted supersaturations ~0.08% half of the BC number could be activated, whereas in cleaner environment needs ~0.14%...
0 downloads 0 Views 1MB Size
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