1. Introduction
Biomass burning aerosol mainly originates from wildfires and constitutes an important natural source of the atmospheric aerosol that influences the Earthâs radiation budget and affects atmospheric circulation and regional climate [
1,
2]. It is primarily composed by light-absorbing components, as brown and black carbon, representing the main source of organic aerosol in the atmosphere [
3]. Biomass burning aerosol plays an important role in radiation budget, because it influences the scattering or absorption of the solar radiation and acts as cloud condensation nuclei or ice nuclei. Moreover, its radiative effect can be either direct or indirect [
4]. The chemical and physical properties of biomass burning aerosol can vary during the transport in the atmosphere, as a consequence of both ageing effects and mixing with aerosol originated by other sources, and significantly influence the impact of wildfires on climate changes and air quality [
5]. The vertical transport of smoke particles through the tropopause remains hitherto unclear, and hence still under investigation; recent studies have reported observation or modeling analyses of smoke particles layers produced by late wildfire episodes in the lower stratosphere [
6,
7,
8,
9,
10]. In this altitude range, the particles are confined for a long time and can reach areas located many kilometers far away from their source region. During the transport, the particles influence the atmospheric circulation and interact with clouds, eventually reducing the precipitation rate [
11]. Moreover, major fire events can even lead to stratospheric ozone depletion in Arctic and Antarctic regions [
12].
Lidars have been extensively used to study long range smoke transport, being able to retrieve the vertical profile of their optical and microphysical properties. Lidar measurements are usually analyzed in synergy with data provided by other instruments, like sun photometers [
13,
14,
15], spaceborne instruments [
7,
16] or nephelometers [
17].
In the Summer of 2017, huge wildfires in the British Columbia region (Canada) led to the injection of an unusually high concentration of biomass burning aerosol in the atmosphere, which was observed even over Europe in the following weeks. An uncommon load of smoke was released on 12
th August 2017 in the troposphere and lower stratosphere over the western area, and during the successive weeks it was transported by air flow all over the northern hemisphere, reaching also the tropics [
18]. Its presence was detected by many European observation stations during the following weeks. In Leipzig on 22
nd August 2017, three lidars analyzed the optical and microphysical properties of two atmospheric smoke layers: one in the troposphere and one in the lower stratosphere; the former showed smaller and less depolarizing particles than the latter one [
19]. In Kosetice, Czech Republic, a smoke layer raising from 12 km to 16 km altitude was detected by a lidar from 20
th to 23
rd August 2017 [
6]. Ref. [
7] reported on the observation of a smoke layer from the spaceborne CALIOP lidar and the ground-based lidars located at the Observatoire de Haute-Provence, in southern France, whereas Ref. [
9] characterized the smoke plume with simultaneous lidar and sun-photometer observations in Lille and Palaiseau. Smoke layers were observed up to 20 km altitude and persisted in the atmosphere until the first months of 2018 [
18,
20].
Between 29
th December 2019 and 4
th January 2020, another noteworthy outbreak of pyrocumulonimbus, due to the Australian bushfire season, was observed in the Southern Hemisphere; it was dubbed as the Black summer due to the unprecedent intensity of the wildfires, which caused vast land damages [
21]. The smoke particles arrived up to the lower stratosphere and the transportation event was observed by lidars located in Punta Arenas (Chile) and RĂo Grande (Argentina) [
22,
23].
Here we report on aged biomass burning aerosol, coming from North America during the Summer of 2017, observed over the observation station of Naples, Italy. The aerosol optical and microphysical properties were characterized by exploiting both a multi-wavelength elastic-Raman lidar and a sun-sky-lunar photometer.
3. Data Analysis
During August 2017 intense fires developed in the British Columbia region (Canada), leading to the injection of an unusually high concentration of biomass burning aerosol in the atmosphere. Large amounts of fire smoke were lifted to the tropopause, even reaching the lower stratosphere. The lifting of these particles to the high troposphere was very fast due to the high radiative power of the wildfire event, preventing mixing with other aerosol or atmospheric constituents and leading to the peculiar compositional and shape properties observed for this event [
19]. In the following weeks, higher altitudes layers of this aerosol were transported by winds from the source site throughout Europe, where they were observed in many locations and with different methodologies.
A first signature of a stratospheric aerosol layer was observed over the Naples station, on 31
st August 2017. This aerosol layer is clearly addressed in the map of the range corrected lidar signals (RCS) registered at 532 nm reported in
Figure 1(a) and extends over an altitude range going from 18.2 to 19.4 km. In the same figure a clear signature of cirrus clouds is present at about 10-11 km of altitude. Unfortunately, it was not possible to characterize this aerosol layer at 355 nm and 1064 nm, due to the lower signal to noise ratio achieved for these two wavelengths at such a high altitude. Aerosol depolarization ratio at 532 nm was (16.7 ± 5.5)%; such a value, which is rather unusual for inter-continental transported fire smoke particles, is in agreement with the observations described in other works reporting measurements carried out in the stratosphere over Europe [
6,
7,
8,
9].
The aerosol wildfire transport episode continued in the successive period, and a very interesting situation was registered on 4
th September 2017, as displayed in
Figure 1(b). The map of
Figure 1(b) evidences the presence of various aerosol layers at different altitudes in the atmosphere over the station. Besides the persistence of a smoke aerosol layer at about 19 km, already addressed on 31
st August 2017 in panel (a) of
Figure 1, it is possible to recognize other three weak layers: the first extends from 11 to 13 km, the second is in the lower troposphere at an altitude between 2 and 4.5 km and the third below 2 km. A thorough characterization of these smoke aerosol layers can be carried out for the data reported in the map of
Figure 1(b), since Raman measurements are also available in addition to the elastic ones. Hereafter, we will focus on this last observation with the aim of characterizing the features of these various aerosol layers.
Figure 2 illustrates the back-trajectories of the air masses that arrived over the measurement area on 4
th September 2017 at the different altitude levels of interest for the three layers identified above. The HYSPLIT model clearly identifies Canada as the common origin of the observed layers; the symbols in the three panels of
Figure 2 show FIRMS maps highlighting the locations of wildfires over the interested area, which was severely damaged by their intense activity during the August of 2017.
The fire locations reported in the three maps of
Figure 2 correspond to those observed in the days at which the back-trajectories, eventually ending over the observation station on 4
th September 2017 at 18:00 UTC, passed over Canada. This analysis confirms that all the three layers observed by MALIA came from Canada and, therefore, are expected to present an important biomass burning component.
With the aim of characterizing the aerosol properties in each layer, we analyzed the lidar profiles measured at three wavelengths in terms of aerosol backscattering (ÎČ), aerosol extinction (α) and depolarization ratio (ÎŽ). 30 minutes time integrated lidar profiles (from 18:10 to 18:40 UTC) are reported in
Figure 3.
In the left panel of
Figure 3, only the profile at 532 nm reaches an altitude up to 19 km due to the higher signal to noise ratio at these heights for such a wavelength, whereas the maximum level for the signals at 355 nm and 1064 nm is around 12-13 km. From the backscatter profiles in the left panel, the layers already identified in the false color map of the RCS signal in
Figure 1(b) are clearly discernible. The backscatter and extinction profiles (left and right panels of
Figure 3) allow estimating the average values of the Lidar Ratio (LR) for the layers located in the troposphere (range R
1 and R
2) and low stratosphere (range R
3), respectively, for both 355 nm and 532 nm. For the tropospheric layers we estimate LR
R1,355=(32 ± 10) and LR
R1,532=(82 ± 17), LR
R2,355=(20 ± 17) and LR
R2,532=(32 ± 9), whereas for the stratospheric one the values are LR
R3,355=(62 ± 16) and LR
R3,532=(82 ± 47). The measured LR are in the range of values reported in the literature for biomass burning aerosols [
19]; the lower LR values in the range R
2 could be indicative of less absorbing particles in this layer [
8]. The average aerosol depolarization ratios at 532 nm (central panel of
Figure 3) for the layers located at R
1, R
2 and R
3 are ÎŽ
R1=(3.9 ± 0.3) %, Ύ
R2=(2.5 ± 0.3) % and Ύ
R3=(6.6 ± 1.8) %, respectively; such values are suggestive of almost spherical smoke particles. In the range R
4 (17.0-19.0 km) a higher value of ÎŽ
R4=(29.8 ± 7.6) % was measured, which is in agreement with those reported in previous works that investigated the summer 2017 inter-continental transport of smoke particles from Canada to Europe [
6,
8,
9], addressing the higher value mainly to different transport ways and relative humidity conditions that lead to different aging processes.
Mean values of ÎČ(z) and α(z) measured at 355 nm and 532 nm in the observed layers were used to retrieve the backscatter-related Ă
ngström (BAE) and extinction-related Ă
ngstrom (EAE) exponents, as well as the average values of the linear depolarization ratio (ÎŽ) and lidar ratio (LR). The estimated values of these lidar-derived aerosol parameters are summarized in
Table 1. The values of BAE range from (1.1 ± 0.7) to (4.4 ± 4.9) and EAE varies from -(1.0 ± 0.5) to (2.9 ± 4.0), in agreement with those reported in earlier works characterizing the properties of the biomass burning aerosol layers observed in Europe. Ref. [
51], analyzed the transport of biomass burning aerosol over Thessaloniki from regions with intense fire activity, reporting BAE values in the range 0.5-2.4. Moreover, in a recent study of Ref. [
52] related to long term lidar observations of biomass burning aerosol over Warsaw, the EAE values range from -0.8 to 2.8, depending on the air masses origin, the lowest values corresponding to aged smoke particles. In the same study, the BAE values range from 1 to 2.7. In the present case, the retrieved values of BAE and EAE increase with altitude, although affected by a larger uncertainty for higher altitudes. This may indicate that the smoke particles in the ranges R
1 and R
2 have larger dimensions compared to those in the R
3 range. Conversely, Ref. [
19] in their study on western Canadian wildfire smoke observed in Leibniz measured higher values in the troposphere than in the stratosphere, which is suggestive of larger particles in the stratosphere probably related to different aging processes.
Mean values of ÎČ(z) and α(z) measured at different wavelengths were also used to retrieve vertically resolved volume particle size distributions đđ(đ)/đđđ(đ) (expressed in a.u.) [
53]. The size distribution averaged over all the measured atmospheric column obtained from the lidar data is reported in
Figure 4(a). The blue solid line refers to the mean value, whereas the two grey dotted curves define the error thresholds, corresponding to the minimum and maximum retrieved values. The lidar-derived size distribution, obtained from the profiles measured at about 18:00 UTC and discussed above, shows a bimodal shape with two peaks, whose mode radius values lie at about 0.20 ”m and 3.46 ”m. The lidar-derived size distribution results in fairly good agreement with the columnar particle size distribution provided by the AERONET sun-photometer data reported in
Figure 4(b). In particular, for the AERONET size distribution measured at the time closest to the lidar observations, i.e. 14:42 UTC, the peak mode radius values sit at about 0.15 ”m and 3.86 ”m, respectively. The differences between the mode radii of the two distributions are likely due to the different measurement time intervals.
The overall consistency of the columnar size distributions derived from lidar and sun-photometer data demonstrates the reliability of the lidar inversion algorithm, thus encouraging gaining further insights on the features of the three aerosol layers recognized in the map of
Figure 1(b). In this respect, it is worth recalling that the size distribution derived from sun-photometer data is not vertically resolved; therefore, the possibility of gaining space resolved size distribution offered by lidar is a unique and complementary source of information to characterize complex aerosol multi-layered conditions, as that occurred on 4
th September 2017.
The lidar derived size distributions for the ranges R
1, R
2 and R
3 are reported in
Figure 5. All the three size distributions are bimodal, similarly to the columnar one discussed above (
Figure 4). However, some interesting differences can be distinguished. For the higher layers at ranges R
2 and R
3, a larger concentration of fine particles is observed, whose peak mode radius is about 0.2 ”m in both cases; the second, less populated fraction of particles is characterized by a peak mode radius of 5.2 ”m and 3.5 ”m, respectively, for the ranges R
2 and R
3. A different situation is, instead, observed for the layer located at the range R
1, i.e. below 2 km of altitudes. For such a layer, the size distribution evidences a predominance of particles in the coarse mode fraction, with a peak radius of about 5.5 ”m. The comparison between the three size distributions in
Figure 5 show that the relative fraction of the fine mode particles progressively increases with the altitude. The larger fraction of coarse particles observed in the lower troposphere can be rationalized as mainly due to deposition phenomena towards the ground.
It is worth to notice that the predominance of the fine mode particles in the atmosphere revealed by the lidar data from 18:10 to 18:40 UTC agrees with the columnar aerosol properties derived by the sun-photometer. In fact, as shown in
Figure 6, AERONET data evidence a progressive increase of the fine mode aerosol fraction in the atmospheric column since 14:00 UTC and an even larger dominance of the small particles (Fine Mode Fraction > 0.8) after 16:00 UTC.
Further information on the aerosol features can be gained by resorting to the columnar parameters provided by the sun-photometer. The columnar AOD
440 and α
440/870 measured by diurnal and nocturnal sun-photometer observations shows AOD
440 and α
440/870 values larger than 0.1 and 1.3, respectively, namely in the afternoon (not shown). This observation, in turn, suggest a higher content of fine particles in the atmospheric column, in agreement with the analyses reported above. In addition, as asserted by Ref. [
54], the observed value of α
440/870 can be considered as representative of pure brown carbon or black carbon cores coated in non-absorbing matter.
The dominance of fine particles in the atmospheric column is further supported by the spectral dependence of the asymmetry factor reported in
Figure 7a, which shows a decreasing trend with the wavelength in the afternoon that can be ascribed to a more efficient scattering at lower wavelengths typical of the fine aerosol. The observed values of the asymmetry factor below 0.7 at 440nm agree with those reported in the literature for biomass burning aerosol; for example, Ref. [
8] registered a mean value of 0.7 at 440nm for aged smoke plumes observed in the Iberian Peninsula on 7
th and 8
th September 2017. Panels (b) and (c) of
Figure 7 also show SSA and AOD absorption evidencing large values of the SSA (> 0.97) and low AOD absorption (< 0.0025) at 440 nm. These values are coherent with the results obtained in previous studies related to biomass burning aerosol, that reported SSA values in the range 0.91-0.99. For example, Ref. [
55] estimated an average SSA at 440 nm ranging from 0.95 to 0.97 for smoke aerosol, with slightly larger values for more aged plumes. Moreover, studies on the biomass burning aerosol observed in Europe and originated from North America [
8,
56,
57] attributed these values mainly to the transformation of the biomass burning aerosol and its mixing with non-absorbing aerosol species during their transport. Ref. [
8] asserted that the low absorption (SSA>0.98) of the smoke particles observed over the Iberian Peninsula on 7
th and 8
th September 2017 might result from a smoldering combustion at the source, that gives rise to aerosol whose absorption properties are lower than those generated during a flaming phase.
The high SSA values discussed above are associated to small values of the imaginary part of the refractive index, as shown in
Figure 8, which since 13:32 UTC lowers to less than 0.001, suggesting a low black carbon fraction in the smoke particles probably due a predominance of smoldering combustion [
58]. The spectral dependence of the imaginary part of the refractive index show a somewhat decreasing trend with the wavelength, possibly due to the absorption by organic carbon [
59]. The real part of the refractive index in the afternoon in
Figure 8 levels at a value of about 1.60 independently of the wavelength, whereas in the morning hours varies between 1.4 and 1.6. These values are consistent with those reported in the literature; for example, Ref. [
60] observed values in the range 1.54-1.69 at 475 nm, with an average of (1.635 ± 0.056) for the low black carbon fires, Ref. [
61] measured values in the range 1.47-1.64 for biomass burning organic aerosol, showing that its variations could depend on the combustion efficiency. Moreover, values of the real part of the refractive index varying from 1.47 to 1.58 were also reported by Ref. [
48] for biomass burning aerosol.
Finally, we would like to anticipate some results on the refractive index estimated by using an improved version of our lidar data inversion method [
27] that will be the subject of a future publication [
62]. These data are useful to further characterize the different aerosol layers recognized by the lidar measurement, since they provide slightly different values of the refractive index for the three aerosol layers located at the R
1, R
2 and R
3 ranges. At a wavelength of 355 nm, the retrieved values of the refractive index are (1.4 ± 0.03i) for R
1, (1.5 ± 0.008i) for R
2, and (1.4 ± 0.01i) for R
3. The higher real refractive index in the layer at R
2 could be due to the ageing processes occurred during the air masses transport; the lower imaginary part of the refractive index suggests a lower black carbon fraction in R
2, in agreement with lower LR values measured by lidar that suggested particles with a lower absorption at those altitudes [
57]. In addition, the columnar value of the refractive index obtained by using our inversion method resulted to be (1.5 ± 0.03i), addressing a good consistency with the columnar data provided by the sun-photometer and also agreeing with the columnar values reported by other authors [
63]. These observations eventually highlight that long-range transported aged smoke particles are less absorptive.
4. Conclusions
Biomass burning aerosols arriving over Naples, Italy, from North America in 2017 were characterized in terms of their optical and microphysical properties with both a multi-wavelength elastic-Raman lidar of the ACTRIS research infrastructure and an AERONET sun-sky-lunar photometer.
Lidar profiles highlighted the presence of four aerosol layers. Two layers were located in the lower troposphere, namely below 2 km (R1) and in the range 2-4.5 km (R2), and the other two weak layers extended from 11 to 13 km (R3) and from 17 to 19 km (R4), respectively. Mean LR values at 355nm and 532nm in the first three layers were obtained and resulted to be equal to (32 ± 10) and (82 ± 17) for the R1 tropospheric layer, (20 ± 17) and (32 ± 9) for the R2 tropospheric layer and (62 ± 16) and (82 ± 47) for the lower stratospheric layer R3, respectively. The aerosol depolarization ratio showed an increasing dependence on the altitude with averaged values of (3.9 ± 0.3)%, (2.5 ± 0.3)%, (5.5 ± 0.5)% and (29.8 ± 7.6)% for the different layers located at progressively larger heights, which is indicative of almost spherical smoke particles at the lower altitudes and aspheric particles at the highest range probably because of aging processes.
Volume particle size distributions were derived from the lidar data for the first three aerosol layers and highlighted that the particles observed in the troposphere have a slightly larger dimension than those present at higher altitude, mainly due to deposition phenomena towards the ground. Moreover, the layers located at ranges R2 and R3 related to smoke particles evidence a larger content of fine mode aerosol with peak mode radius values at 0.2 ”m and 0.1 ”m, respectively. Also, AERONET columnar properties revealed a dominance of fine particles in the atmosphere with low absorption properties (absorption AOD < 0.0025 and SSA > 0.97) and addressed aerosol properties in good agreement with those evidenced by the lidar observations.
Finally, the refractive index for the three lower aerosol layers recognized by the lidar measurements was estimated by resorting to an inversion method capable of retrieving the microphysical particles properties. Our results, although preliminary, show a higher real part and a lower imaginary part of the refractive index in the R2 later, suggesting a lower black carbon fraction and, in agreement with lower LR values measured by Lidar, less absorbing particles in this altitude range. Moreover, the space-resolved results on the refractive index are consistent with AERONET sun-photometer retrieved values, thus evidencing the reliability and capability of lidar characterization of atmospheric aerosol in a very interesting condition of multiple aerosol layers originating from Canadian fires overpassing the observation station.