Heating rate of light absorbing aerosols: time-resolved measurements

Feb 23, 2018 - Heating rate of light absorbing aerosols: time-resolved measurements, the role of clouds and source-identification. Luca Ferrero, Grisa...
0 downloads 14 Views 20MB Size
Subscriber access provided by UNIV OF DURHAM

Article

Heating rate of light absorbing aerosols: time-resolved measurements, the role of clouds and source-identification Luca Ferrero, Grisa Mocnik, Sergio Cogliati, Asta Gregoric, Roberto Colombo, and Ezio Bolzacchini Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.7b04320 • Publication Date (Web): 23 Feb 2018 Downloaded from http://pubs.acs.org on February 25, 2018

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.

Environmental Science & Technology 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 26

Environmental Science & Technology

1

Heating rate of light absorbing aerosols: time-

2

resolved measurements, the role of clouds and

3

source-identification

4

Luca Ferrero1,2*, Griša Močnik3,4, Sergio Cogliati2,5, Asta Gregorič3,6, Roberto Colombo2,5, and Ezio

5

Bolzacchini1,2

6

1

POLARIS Research Centre, Department of Earth and Environmental Sciences, University of Milano-

7

Bicocca, 20126, Milano, Italy 2

8

GEMMA Centre, Department of Earth and Environmental Sciences, University of Milano-Bicocca,

9

20126, Milano, Italy 3

10 4

11

Aerosol d.o.o., Kamniška 41, SI-1000 Ljubljana, Slovenia

Department of Condensed Matter Physics, Jozef Stefan Institute, Jamova 39

12 13

SI-1000 Ljubljana, Slovenia 5

Remote Sensing of Environmental Dynamics Lab., Department of Earth and Environmental Sciences,

14 15

University of Milano-Bicocca, 20126, Milano, Italy 6

Center for Atmospheric Research, University of Nova Gorica, SI-5000 Nova Gorica, Slovenia

16 17

Luca Ferrero: [email protected]

18

Luca Ferrero, POLARIS Research Centre, GEMMA Centre, Department of Earth and Environmental

19

Sciences, University of Milano-Bicocca, 20126, Milano, Italy.

20

Tel. ++39 02/64482814; Fax. ++39 02/64482839

21 22

ABSTRACT

23

Light Absorbing Aerosols (LAA) absorb sunlight and heat the atmosphere.

24

This work presents a novel methodology to experimentally quantify the heating rate (HR) induced by

25

LAA into an atmospheric layer. Multi-wavelength aerosol absorption measurements were coupled with ACS Paragon Plus Environment

1

Environmental Science & Technology

Page 2 of 26

26

spectral measurements of the direct, diffuse and surface reflected radiation to obtain highly-time

27

resolved measurements of HR apportioned in the context of LAA species (Black Carbon, BC; Brown

28

Carbon, BrC; dust), sources (fossil fuel, FF; biomass burning, BB), and as a function of cloudiness.

29

One year of continuous and time-resolved measurements (5 min) of HR were performed in the Po

30

Valley. We experimentally determined: 1) the seasonal behavior of HR (winter: 1.83±0.02 K day-1;

31

summer: 1.04±0.01 K day-1); 2) the daily cycle of HR (asymmetric, with higher values in the morning

32

than in the afternoon); 3) the HR in different sky conditions (from 1.75±0.03 K day-1 in clear sky to

33

0.43±0.01 K day-1 in complete overcast); 4) the apportionment to different sources: HRFF (0.74±0.01 K

34

day-1) and HRBB (0.46±0.01 K day-1); 4) the HR of BrC (HRBrC: 0.15±0.01 K day-1, 12.5±0.6% of the

35

total) and that of BC (HRBC: 1.05±0.02 K day-1; 87.5±0.6% of the total).

36 37

KEYWORDS: Light Absorbing Aerosol, Radiation, Black Carbon, Brown Carbon, Sources, Heating

38

Rate, Field spectroscopy, Po Valley.

39 40

1.

41

Light absorbing aerosol (LAA), such as black carbon (BC), brown carbon (BrC) or dust, directly

42

interacts with the solar radiation (direct radiative effect; DRE) absorbing it and exerting a positive

43

atmospheric forcing. This process releases energy (heat) to the surrounding air with a specific heating

44

rate (HR).[1-9]

45

Among the various LAA species, BC, generated from incomplete combustion of both fossil fuels (FF,

46

i.e. diesel, coal) and biomass (BB)[2,10,11], represents the second most important anthropogenic climate-

47

forcing agent after CO2.[2,12] Bond et al.[2] report a mean DRE at the top of atmosphere (TOA) in the

48

range of 0.2-1.0 W m−2 (average: +0.71 W m-2) even if regionally should be larger (up ~10 W m-2 over

49

Asia; 0.0-4.0 W m-2 over Europe).[7,13] The LAA do not only influence the DRE at TOA, but also its

50

vertical distribution. In this respect, BC heats the atmosphere (average +2.6 W m-2, 10-20 W m-2 within

51

regional plumes and instantaneously up to +50-80 W m-2) and dims the surface (mean value -1.7 W m-2;

52

instantaneously down to -40 – -70 W m-2).[6,7,12,14] Thus, BC HR varies as a function of height (range:

53

0.3 - 8 K day-1)[6,11,15,16]; over Asia, the LAA HR was found larger than that of CO2.[16]

54

Considering the other LAA species, the BrC (the absorbing component of organic matter) is receiving a

55

growing attention as it contributes ~10-30%[17,18] to the total LAA absorption on a global scale. It causes

56

an atmospheric DRE of 2.6 W m-2 over the Indo-Gangetic Plain that can rise up to 3.6 W m-2 during

57

high BB days.[18] Finally, also dust aerosol can exert a positive atmospheric DRE; instantaneously it was

58

estimated up to ~80-100 W m-2,[3,19] and HRs of ~2-3 K day-1 were found in elevated dust layers over

59

Korea.[20]

INTRODUCTION

ACS Paragon Plus Environment

2

Page 3 of 26

Environmental Science & Technology

60

Unlike CO2, which is distributed quite homogeneously in the atmosphere, the LAA species are

61

characterized by a significant spatial variation, due to their relatively short lifetime (i.e., weeks).[21,22]

62

Thus, on a regional scale, the radiative effects of LAA can be larger than that of greenhouse gases[14,16]

63

triggering different local/regional feedbacks: alterations of the atmospheric thermal structure, of the

64

atmospheric stability, of the cloud distribution (cloud “burn-off”, and cloud cover) and of synoptic

65

winds such as the monsoon.[2,3,6,10,12,14,15,23-27] Furthermore, the retreat of important glaciers (i.e.

66

Himalayan) can be also influenced by the LAA induced HR at elevated altitudes.[14,28] Due to the

67

importance of the aforementioned feedbacks on climate, an experimental determination of the HR

68

induced by LAA is called for.

69

Actually, most of the forcing values and atmospheric feedbacks reported in literature rely on the use of

70

several theoretical models.[2,10] Model outputs are continuous in space (horizontally and vertically) and

71

time, allowing the description of the impact of LAA on climate at the global scale. However,

72

uncertainties in the modelling chain (i.e. LAA mass concentrations, optical properties and mixing

73

state)[13,24] result in large discrepancies among model outputs. This amplifies the uncertainty of the LAA

74

influence on the global and regional climate.[1] As a conclusion, Andreae and Ramanathan[1] showed the

75

need of better observational constraints at all scales (local/regional/global), while Chung et al.[17] called

76

for model independent and observationally based estimates of the LAA direct radiative forcing. One

77

step further in this direction is given by hybrid approaches, based on experimental measurements of

78

LAA properties (mass/number concentrations) coupled with optical/radiative transfer calculations.[6,15,29-

79

31]

80

be affected by the same uncertainties of the modelling approach. Moreover, radiative transfer

81

calculations are usually carried out under clear-sky approximation, which represents a strong limitation

82

of the technique (it is useful only if the experimental campaign is conducted during such sky

83

conditions).[6,16] Hybrid approaches based on vertical profile data describe both the atmospheric

84

behavior and the related feedbacks, but they are usually limited to short/intense campaigns.[6,15,29,32-37]

85

On the contrary, hybrid approaches based on ground observations only are usually continuous in time,

86

but they are limited in the vertical direction.[30,31] A complete experimental determination of the

87

radiative power absorbed into the atmosphere (and the related HR) requires the measurement of the

88

vertical divergence of the net radiative flux with altitude.[9,16] With this approach, the obtained HR are

89

usually limited to short/intense campaigns and to specific locations; moreover they are comprehensive

90

of the whole atmospheric composition (greenhouse gases and LAA) thus again theoretical models need

91

be used (coupled with concentrations measured in situ) to unravel the contribution of each single

92

atmospheric component.[11,16]

They have the advantage that are based on experimental measurements but, at the same time they can

ACS Paragon Plus Environment

3

Environmental Science & Technology

Page 4 of 26

93

Thus, a method to provide continuous, LAA specific, experimental determination of atmospheric

94

forcing and HR is required. Theoretically, the method should be simple and accessible enough to be

95

used in reference stations belonging to widely distributed monitoring networks (i.e. the European

96

Monitoring and Evaluation Programme, EMEP network) and should provide an estimation of the

97

contribution of each aerosol source to the atmospheric HR.

98

This paper fills this gap, describing a new methodology to experimentally determine both the radiative

99

power density absorbed into a near-surface atmospheric layer, and the consequent HR induced by LAA.

100

The main goal of this method is to provide continuous experimental observations of LAA HR by means

101

of a simple and cost-effective technique (see supporting information) applicable in established aerosol

102

monitoring networks.

103

Other very important goals of the new method are: 1) the apportionment of the HR as a function of the

104

LAA species (i.e. BC, BrC, dust) and the LAA sources (i.e. FF, BB); 2) the experimental determination

105

of HR in any sky condition; 3) the apportionment of different radiation components (i.e., downward

106

direct/diffuse fluxes and the reflected) on the HR, giving the possibility to determine the influence of

107

clouds on the HR.

108 109

2.

110

This section describes both the theoretical considerations (section 2.1) and the experimental design

111

(section 2.2) of the novel method for the experimental determination of the LAA induced HR.

METHODOLOGY

112 113

2.1

114

The instantaneous aerosol direct radiative effect (DRE; W m-2) at a z altitude can be quantified as the

115

change in the net radiative flux between the atmospheric conditions with and without the aerosols (Qaer,z

116

and Qw/o-aer,z, respectively) across the surface at z altitude:[6,9]

117

 =  , −  / ,

118

Considering an atmospheric layer of thickness ∆z, the difference between the DRE at the top and the

119

bottom of this atmospheric layer represents the instantaneous radiative power density absorbed by the

120

aerosol (∆DRE; W m-2):[6-9]

121

 =  − 

122

From ∆DRE, the instantaneous heating rate (HR; K day-1) of the same atmospheric layer can be

123

computed as follows:[6-9]

124

 =

THEORETICAL CONSIDERATIONS

 

(1)

(2)

 

= −



(3)



ACS Paragon Plus Environment

4

Page 5 of 26

Environmental Science & Technology

125

where ∂T/∂t represents the instantaneous HR, g is the gravitational acceleration constant, Cp (1005 J kg-1

126

K-1) is the isobaric specific heat of dry air, ∆P is the pressure difference between the top and the bottom

127

of the considered layer.

128

A more useful definition of the HR is based on the thickness of the atmospheric layer (∆z), and can be

129

obtained introducing the hydrostatic equation (dp = -ρgdz) in eq. 3:[9]

130

 =

131

where ρ represents the air density (kg m-3). The last term of eq. 4 (∆DRE/∆z) represents the radiative

132

power absorbed by the aerosol for unit volume of the atmosphere (W m-3) and was previously defined in

133

Ferrero et al.6 as the absorptive direct radiative effect (ADRE) of LAA; the ADRE is the vertical

134

derivative of the DRE (dDRE/dz; the implications of the derivative are detailed and discussed in section

135

4). Hence, the HR becomes:[6]

136

 =

137

Thus, any method able to determine the ADRE at high time resolution will produce continuous time

138

series of HR. This possibility is given by the simple observation that ADRE is expressed per unit volume

139

of air, coherent with the concentrations of LAA routinely measured in monitoring networks. Below the

140

rigorous methodology used for the experimental determination of the ADRE is reported.

141

Let us consider a near-surface atmospheric layer (but can be applied to any atmospheric layer; see

142

section 4 for further details) of thickness ∆z on which a n-th kind (direct or diffuse or reflected) of

143

monochromatic radiation ray Fn(λ,θ) of wavelength λ strikes with a zenith angle θ (Figure S1a, supporting

144

information). The amount of radiation absorbed by the aerosol within the present layer is as follows:[9,38]

145

#(%,&) = )#(%,&) (1 − +% ),1 − - .//0 1

146

where (1- ωλ)(1-e-τλ/µ) represents the fraction of light absorbed within the layer[9,38] and is function of: ωλ

147

which is the single scattering albedo of the aerosol within the atmospheric layer, τλ the aerosol optical

148

depth and µ, the cosine of θ (cosθ). The ωλ and τλ terms can be computed from the aerosol extinction,

149

scattering and absorption coefficients (bext(λ), bsca(λ) and babs(λ)):

150

+% =

151

7% = 8 9:(%) ;