Study of Chemical and Optical Properties of Biomass Burning Aerosols during Long-Range Transport Events toward the Arctic in Summer 2017

Biomass burning related aerosol episodes are becoming a serious threat to the radiative balance of the Arctic region. Since early July 2017 intense wildfires were recorded between August and September in Canada and Greenland, covering an area up to 4674 km2 in size. This paper describes the impact of these biomass burning (BB) events measured over Svalbard, using an ensemble of ground-based, columnar, and vertically-resolved techniques. BB influenced the aerosol chemistry via nitrates and oxalates, which exhibited an increase in their concentrations in all of size fractions, indicating the BB origin of particles. The absorption coefficient data (530 nm) at ground reached values up to 0.6 Mm–1, highlighting the impact of these BB events when compared to average Arctic background values, which do not exceed 0.05 Mm–1. The absorption behavior is fundamental as implies a subsequent atmospheric heating. At the same time, the AERONET Aerosol Optical Depth (AOD) data showed high values at stations located close to or in Canada (AOD over 2.0). Similarly, increased values of AODs were then observed in Svalbard, e.g., in Hornsund (daily average AODs exceeded 0.14 and reached hourly values up to 0.5). Elevated values of AODs were then registered in Sodankylä and Andenes (daily average AODs exceeding 0.150) a few days after the Svalbard observation of the event highlighting the BB columnar magnitude, which is crucial for the radiative impact. All the reported data suggest to rank the summer 2017 plume of aerosols as one of the biggest atmosphere related environmental problems over Svalbard region in last 10 years.


Introduction
The ongoing climate change is leading to an increase in surface temperatures all over the globe [1] and as a consequence numbers of wild forest fires increase [2]. Biomass burning (BB) aerosols contribute to climate change causing a further increase in temperatures that affects the frequency of fires by triggering a feedback mechanism. According to Gillet et al. [3], the wildfire area burned in Canada, for example, is continuously increasing due to the human-induced climate change.
The chemical composition and the optical properties of BB aerosols depend not only on the type of forest burned, but also on the environmental conditions in which combustion takes place [4,5]. Generally, organics and black carbon are the two main BB aerosol components followed by sulphate, ammonium, nitrate and chloride [6][7][8][9].
Biomass burning aerosols are responsible for many environmental issues [1]. They have serious influence on human health and due to their optical properties and consequent impact on radiation budget, BB aerosols affect the planet climate [7][8][9][10][11][12][13][14][15][16]. These aerosol particles may both absorb and scatter solar radiation, thus play a key role both in modifying atmospheric visibility and photochemistry in the boundary layer [17,18] and they can act as cloud condensation nuclei [19][20][21]. Despite this climatic impact, estimations of the quantitative effects of BB aerosols on radiation balance are difficult to make and rather of a complex matter [22]. BB particles may cause short-term cooling, but their gases cause warming, the effect, which over the extended time, is of greater magnitude than the cooling effect. The subsequent IPCC reports provide information about the level of radiative forcing of BB aerosols, and they range from −0.2 W/m 2 in the Second and Third Assessment Reports, through +0.03 W/m 2 in the Fourth Assessment Report, to 0.0 (−0.2 to +0.2) W/m 2 in the Fifth Assessment Report [1]. However, recently, Ferrero et al. [23] experimentally showed that the BB aerosol particles are 1.63 ± 0.03 times more efficient than pure black carbon, when heating of the atmosphere is considered. This proves that an increasing number of BB events is an emerging issue, which we will have to deal with in the near future.
In the most extreme cases, wildfires can produce aerosols and trace gases in large quantities. These can be transported and distributed across the globe very fast (days/weeks) and at long distances [24]. In some cases, aerosols produced from large wildfire activities can, via clouds, reach stratosphere. The impact of such particles on the radiative budget as well as chemistry of stratosphere has been discussed by Fromm et al. [25], however, it still needs further studies. BB aerosols are constituted by a wide range of carbonaceous substances, which are abundant in carbon oxidation features, which in turn, have impact on optical properties (e.g., brown carbon) [26,27]. Black carbon (BC) is a main constituent of BB events, which is predominantly associated with absorption of solar radiation and is typically addressed as the most important optical absorbing fraction of aerosols. BC is typically a component of air pollutants, which are produced as a result of several events, e.g., boreal forest fires or shipping activities [28]. BB aerosols are usually internally mixed with other aerosol chemical components, such as ammonium nitrate and ammonium sulphate [29][30][31][32].
Due to its pristine environment, its remoteness and the consequent low anthropogenic influence, the Arctic atmosphere has been regarded as extremely clean, however, this view is being changed [33]. Airborne pollutants can be transported from lower latitudes into the Arctic region or directly emitted into the Arctic area from local sources.
Aerosol properties (chemical and optical) seasonally vary in the Arctic due to changes both in sources (outside or inside the Arctic region) and in meteorological conditions that may allow or inhibit the transport/local formation from source regions [34,35]. Stohl et al. [36] described the late winter/springtime period as affected by the Arctic haze: an inflow of air pollution dominated by accumulation mode aerosols that can be transported into the Arctic both at low-level, or with an uplift outside the Arctic. Conversely, during summer, the high continental pressure disappears and the transport from medium latitudes becomes less relevant.
Atmosphere 2020, 11, 84 3 of 27 Thus, during summer, any emission related to forest fires/biomass burning [37,38] or to local anthropogenic sources (e.g., increasing of shipping emissions/gas flaring) [35,39,40] becomes very important as far as the climate change is concerned.
So far, the Arctic haze has been regarded as the most important aerosol phenomenon in the Arctic. However, recent information on air pollution events in the Arctic, show that the BB episodes are possibly becoming a very important source of atmospheric pollution for the Arctic region, including Svalbard ( [41][42][43][44][45][46][47]).
High levels of BB aerosols in summertime can be very problematic as far as the climate change is concerned, since interactions with solar radiation are at their maximum and BB aerosol particles are very effective in heating the atmosphere [23]. This is related to the fact that interactions with solar radiation are then at their maximum [23] contributing locally to an intensification of the Arctic amplification and globally to an average increase in temperatures that can trigger again higher frequency of forest fires [3]. BB aerosols contain a large amount of light-absorbing material, such as black carbon, and other organic matter, such as vanillic acid, isovanillic acid, homovanillic acid, syringic acid, syringaldehyde, ferulic acid, and p-coumaric acid [48] which, substantially, make up brown carbon particles (BrC).
In general, not only the BC, BrC, and BB aerosols interact with solar radiation directly in the atmosphere but settling on snow and ice reduce albedo and thus facilitate melting. Recent measurements have shown the noticeable increase of CO, ozone, and aerosol particles advected from wildfires in Eastern Europe into the European Arctic [49]. Warneke et al. [50] showed contribution of wildfires in central Asia and Russia to an increased pollution in Alaska. Additionally, in North America, a number of strong wildfire episodes, has been increasing due to anthropogenic related warming [42,44].
Moreover, from an optical point of view, during last ten years (2004, 2005, 2006, 2008, 2011, and 2015), AODs in the Arctic reached high values, ranging from 0.3 up to 1. Such high values resulted both from significant springtime Arctic haze events and, what is important for the aim of this study, transport of BB related aerosols during summertime BB events [37,38]. According to several reports from experimental measurements and models, the year 2017 was particularly important for atmospheric pollution events in the Arctic, with BB particles reaching different altitudes in the atmosphere, up to 20 km a.s.l. [51][52][53][54].
In this paper we present the results of comprehensive studies of aerosol properties during a long-range transport of biomass burning from Canada and Greenland to Svalbard region in summer 2017. The research team included researchers who carried out measurements at different locations on Spitsbergen and outside of the Svalbard Archipelago taking advantage of the Oceania research vessel (r/v Oceania) which enabled to perform measurements at sea. An ensemble of passive and active remote sensing methods, modeling of aerosol optical properties, simulations of variability of AOD, have been described in this work in order to properly evaluate impact of that Biomass Burning event on Svalbard.

Period and Area of Investigation
The main event of BB aerosol transport under investigation in this study was observed between 28 July and 26 August 2017 in the Svalbard region. Later, a less pronounced event was also observed in the first half of September.
Measurements of several parameters related to the aerosol properties were carried out in two stations of the Spitsbergen Island: Ny-Ålesund and Hornsund.
Ny-Ålesund (78 • 55 25" N, 11 • 55 21" E) is a research town in the Oscar II Land on the island of Spitsbergen. The research stations are placed on the southern shore of Kongsfjorden [55]. Aerosol in situ measurements were made in the Gruvebadet laboratory (78.918 • N, 11.895 • E; 61 m a.s.l.), which is placed 800 m south-west of Ny-Ålesund. The Polish Polar Station in Hornsund (77 • 0 0" N, 15 • 33 0" E) is located in Isbjørnhamna Bay, on the northern shore of the Hornsund Fjord. It is placed 10 m a.s.l., at the mouth of the fjord.
Additional measurements (optical properties of aerosols) were also planned and only partially made due to bad weather conditions (mostly clouded, foggy, and high swell) in the area west and south-west of Spitsbergen from board of r/v Oceania. Air masses, which move in the direction of Svalbard may be disturbed by local orographic conditions, related to the location of the Spitsbergen stations. Thus, the aerosol particle distribution is strongly related to the specific station. Additionally, water areas, which surround Svalbard archipelago have their own specific features in terms of aerosol production.

Instruments and Models
During the entire August 2017 event an ensemble of aerosol instruments were employed for measurements (Table S1). For the purpose of this paper we divided them into four groups: ground-based, columnar integrated, vertically resolved, and models.

Ground-Based Data
Aerosol properties were obtained at a ground from a large number of in-situ measurements, which were made at the Gruvebadet station in Ny-Ålesund. Bulk and size-segregated aerosol samples were collected. PM 10 bulk aerosol was collected through two parallel TECORA SkyPost samplers (EN 12341; PM 10 sampling head, flow 2.3 m 3 h -1 ; PTFE and Quarz fiber filters, Ø = 47 mm) in actual conditions. Size-segregated aerosol (<1, 1-2.5, 2.5-10, and >10 µm) was sampled by means of a DEKATI 4-stage impactor; polycarbonate membranes (25 mm diameter and 0.1 µm nominal porosity) were used as substrate for the first three stages (1-2.5, 2.5-10, and >10 µm), while a Teflon filter collected the sub-micrometric particles (<1 µm) in the last stage). The sampling flow rate was 1.7 m 3 h -1 and for all the samplings they were kept constant thanks to a continuous pressure and temperature monitoring.
The first TECORA SkyPost collected PM 10 sampled for 24 h on Teflon filters (Pall R2PJ047). The PM10 collection efficiency is reported and described in Ferrero et al. [8]. Here, we briefly resume the concept and we refer to the above quoted publication. PM 10 , which is usually described as the mass concentration of PM having aerodynamic diameters below 10 µm, is more precisely defined as the mass concentration of the PM collected with a~50% efficiency at 10 µm of aerodynamic diameter following the European/EPA regulations. As the collection efficiency represents the probability that a particle would be collected onto the filter, PM 10 may contain particles larger than the 'cut off' aerodynamic diameter of 10 µm as reported in Ferrero et al. [8]. They were conditioned for 48 h (at 25 • C and 50% relative humidity) before and after the sampling. After conditioning, filters were weighted by a 5-digit microbalance (Sartorius ME235P) with a reproducibility error lower than 5% (experimentally evaluated) [40]. Those 24 h Teflon PM 10 samples were analyzed to determine the water-soluble ionic fraction. Half of each filter was extracted for 20 min, by ultrasonic bath, in 10 mL of ultrapure water (Milli-Q, 18.3 MΩ cm resistivity). Inorganic and organic anions together with inorganic cations were measured by a triple Dionex ion-chromatography system, equipped with electrochemical-suppressed conductivity detectors. A Dionex AS4A-4 mm analytical column with a 1.8 mM Na 2 CO 3 /1.7 mM NaHCO 3 eluent, was used for anions (Cl - 10 samples were analyzed to determine the organic carbon (OC) and elemental carbon (EC) by means of a thermo-optical transmission (TOT) method following the NIOSH protocol as reported by Ferrero et al. [40]. From the OC, the organic matter (OM) was determined multiplying the OC by 2.1 [58], typical for remote sites with a large fraction of secondary aerosols. All the filter treatments were carried out under a class-100 laminar-flow hood, in order to prevent any contamination risks.
In addition to the chemical measurements, aerosol optical properties were also determined using a particle soot absorption photometer (PSAP), that allowed the measurement of absorption coefficient (b a ) at three wavelengths (440, 530, and 660 nm) [59], and a M903 Nephelometer, measuring the scattering coefficient (b s ) at 530 nm, both from Radiance Research. The manufacturer declares a detection limit for the two instruments equal to 0.2 and 1 Mm -1 , respectively, but both of them have demonstrated higher performances if specific precautions are taken [60,61]. The PSAP measurements have been corrected for multiple scattering and non-purely absorbing aerosol effects following Haywood and Osborne [62] and the three values at different wavelengths were used to calculate the Ångström absorption exponent (AAE). The combination of absorption and scattering properties allows the evaluation of the single scattering albedo (SSA) at 530 nm. Both these two quantities are intensive properties of the particles. The AAE is the slope of the absorption coefficient vs. wavelength on a log-log scale and is related both to size, age, and carbonaceous content of the particles [63]. SSA is evaluated as scattering over extinction (scattering plus absorption) and is related to the overall fraction of absorbing elements in the aerosols [64]. Being the two derived quantities, their accuracy strongly depends on the absolute b a and b s values. It is worth noting here that the measured aerosol properties should be considered to pertain to dry conditions, due to the exsiccation they experience entering the instruments. If the environmental original relative humidity is above about 60%, this leads to a decrease in the scattering coefficient and in turn affects also SSA.

Columnar Integrated Data
The AOD and related parameters were measured using sun photometers, which provide column integrated aerosol multi-wavelength extinction of the solar radiation through the atmosphere using the Langley method. These values are measured within the VIS and near-IR parts of the solar spectrum [65]. From sun photometer measurements several important aerosol optical parameters can be obtained, including: aerosol optical depth (AOD), Ångström exponent (AE), and the asymmetry parameter (g). AODs obtained from Hornsund and Ny-Ålesund were acquired with Cimel sun photometers: CE-318 and CE-318-T, respectively operated within the AERONET network. These sun and sky radiometers are equipped with spectral interference filters centered at the following wavelengths: 340, 380, 440, 500, 670, 870, 1020, and 1640 nm [66][67][68]. These instruments work in real time with a time-resolution of 15 min in free-cloud conditions [67]. The AOD measurement uncertainty at level 2.0 for the Cimel sun photometer is ±0.01-0.02 [68].
While sun photometers placed in Hornsund and Ny-Ålesund allowed to acquire continuous time-series of measurements at a fixed location, data obtained on board of r/v Oceania were supposed to complement the stationary measurements with data from different locations around Svalbard. Microtops II sun photometers were used due to their high accuracy and easy use (hand-held devices) for measurements of aerosol optical depth, irradiance (direct) at different wavelengths as well as water vapor column. The instruments provided AOD measurements at 5 channels: 440, 500, 675, 870, and 936 nm. AOD measurements obtained using Microtops II were taken as routines of 5 single measurements, which lasted around 1 min (~10 s for a single measurement run), that followed the AERONET procedures and the measurement uncertainty is less than ±0.02 [69]. In the described study 1.5 and 2.0 level data were considered and for analyses we used the lowest value out of these routine set of 5 measurements, connected to least error.
Finally, all of these parameters have been supported by the additional data, such as optical thickness of clouds, satellite AOD, cloud fraction from MODIS (moderate resolution imaging spectroradiometer), which widely used for analyzing the aerosol distributions from regional to global scales [70]. The MODIS sensors take~5 min to collect a~2300 km scene during the overpass [8,71] that, occurs two timed a day (~10.30 and~13.30) above Svalbard region. MODIS data are gridded for aerosol particles to a 10 × 10 km resolution, while for clouds, to a 5 × 5 km resolution. We have used level 2.0 data [72][73][74][75].

Vertically Resolved Data
The Koldeway Aerosol Raman Lidar (KARL) was operating within the Atmosphere Observatory of the German Base in Ny-Ålesund. The lidar is managed by the Alfred Wegener Institute for Polar and Marine Research in Potsdam. KARL is a "3β+2α+2δ+2wv" Raman-lidar, which measures aerosols at 3 wavelengths (355, 32, and 1064 nm); backscatter coefficient (β) and extinction coefficient (α) are measured at 2 wavelengths (355 and 532 nm); the depolarization (δ) in another 2 colors (407 and 660 nm); and the water vapor (WV) [75]. Due to the high-level (~600 m) overlap of the Raman lidar also data from a CHM-15 Ceilometer (from Lufft), operating at 1064 nm, was also considered in this work. The Ceilometer has a resolution of 15 m, from 30 m to 15 km altitude (with vertical resolution of 15 m) [76]. Even though ceilometers are most often applied in cloud measurements, they have been used to study vertical distribution of aerosol optical properties and they can also be applied (using time integration techniques) to estimate total sky coverage and they can show the presence of aerosol plumes.
In addition to KARL and ceilometer measurements (continuous in time but at a fixed location), a set of additional data was acquired from the CALIPSO products-the Lidar Images, Version 4.10 [77,78]. These data include information about 532 nm total (parallel + perpendicular) attenuated backscatter, 532 nm perpendicular attenuated backscatter, depolarization ratio and aerosol subtype chosen based on the orbit tracks. The CALIPSO orbit tracks overpass the investigated region (Svalbard) from 2 to 3 times a day, while surrounding region-over 6 times a day, thus providing a complementary, spatially resolved information.

Preparation of Modelled Datasets
Climate models, such as NAAPS [79] or CAMS [80,81] were also used to simulate aerosol features in order to fulfill the spatial description of the BB events.
The Navy Aerosol Analysis and Prediction System (NAAPS) is a system used for the prediction of the distribution of aerosols in the troposphere. The modified version of the model is fed with global meteorological fields from the Navy Operational Global Atmospheric Prediction System (NOGAPS) [82] analyses and forecasts on a 1 • × 1 • degree grid. Results are provided at 6 h intervals and 24 vertical levels reaching 100 mb (12 levels for initial runs then 18 levels beginning 1998.06.24.12 and 24 levels beginning 2002.09.17.00.) Current features of the NAAPS model include [80]: the use of operational dynamics (e.g., 'REAL' weather), a forecast of 120 h, a near-real-time operation, global coverage, dust and smoke simulations.
Full information on meteorological data (wind speed and wind direction, precipitation, air pressure, cloudiness, including dominating cloud types, relative humidity, and horizontal visibility) were obtained from the World Meteorological Organization (WMO) stations in Hornsund (01003), Ny-Ålesund (01007) and Longyearbyen (01008). We separated wind directions into 16 classes (with a 22.5 • step).
The FLEXPART ("FLEXible PARTicle dispersion model") model has been applied to describe air mass paths (with a 1 • × 1 • grid) at vertical levels not exceeding 3000 m a.s.l. This was made in polar stereographic projection.
The Copernicus Atmosphere Monitoring Service (CAMS) model, which was used to simulate aerosol features has originated from the EU Global Earth-system Monitoring using Space and in-situ data (in short: GEMS) as well a series of Monitoring Atmospheric Composition and Climate (MACC) programs of the European Centre for Medium-Range Weather Forecast (ECMWF). The CAMS provides reanalysis of the AODs (550 nm) and atmospheric aerosol composition for 5 types of tropospheric aerosol particles: sea salt, sulphates, dust, black carbon, and organic matter.

Results
A detailed description of the obtained results begins here with a description of the biomass burning event (Section 3.1), which is followed by a description of the chemical tracers and related optical properties measured at ground in Ny-Ålesund (Section 3.2). Then the BB impact on aerosol optical properties in the atmospheric column is addressed (Section 3.3).

Biomass Burning Events
At the beginning of July 2017, the Canadian wildfire season started reaching just half of the intensity registered during the last 10 years of activity (Canadian Wildland Fire Information System) [83]. However, later in July those fires spread and 24 fires were recorded and the area burned increased rapidly. By mid-August, over 1000 separate fires had been recorded and the burned area grew to almost 9000 km 2 . In addition, no precipitation was observed (Canadian Wildland Fire Information System) [83]. On 21 August 2017, a large number of independent fires merged into a single large wildfire covering an area of 4674 km 2 in size. By late August 2017, the fire activity begun to recede, however, its effects (smoke and dust) were still observed in September 2017 ( Figure 1). Additionally, a wildfire in Greenland was spotted by satellites on 31 July 2017 and it was burning some 40 miles from the Greenland ice sheet ( Figure 2) [84,85].
Atmosphere 2019, 10, x FOR PEER REVIEW 7 of 27 [83]. However, later in July those fires spread and 24 fires were recorded and the area burned increased rapidly. By mid-August, over 1000 separate fires had been recorded and the burned area grew to almost 9000 km 2 . In addition, no precipitation was observed (Canadian Wildland Fire Information System) [83]. On 21 August 2017, a large number of independent fires merged into a single large wildfire covering an area of 4674 km 2 in size. By late August 2017, the fire activity begun to recede, however, its effects (smoke and dust) were still observed in September 2017 ( Figure 1). Additionally, a wildfire in Greenland was spotted by satellites on 31 July 2017 and it was burning some 40 miles from the Greenland ice sheet ( Figure 2) [84,85]. These BB events generated huge quantities of aerosols, whose atmospheric transport was observed between 28 July and 26 August 2017 in the Arctic (Sections 3.2 and 3.3). According to the NASA database, the Visible Infrared Imaging Radiometer Suite (VIIRS) on the Suomi NPP satellite captured particularly heavy smoke originating in northern Canada as of August 15 and 16, 2017, which was transported to the Arctic over Greenland ( Figure 2).
In order to confirm the source regions of the aerosol particles transported over Svalbard during the BB events, air mass trajectories were calculated at altitudes of 500, 1000, and 1500 m a.s.l., using the FLEXPART model ( Figure S1).
The trajectories have been calculated for a period between 18 and 26 August 2017. The results are presented in Figure S1 and they demonstrate that in August 2017, Greenland and Canada wildfires were the sources of aerosols whose evolution and properties are discussed hereinafter.

Chemical Speciation and Optical Properties at Ground-Level (Svalbard)
During the entire August 2017, an ensemble of aerosol in-situ chemical and optical measurements were carried out at the Gruvebadet station in Ny-Å lesund. These measurements are representative of the portion of the atmospheric aerosols that can reach the surface as demonstrated by vertical aerosol profiles measured at Gruvebadet station using a tethered balloon system since 2011 [8,37,[86][87][88]. Atmospheric PM10 concentrations and the related chemical composition, measured at Gruvebadet are presented in Figures 3 and 4. Figure 4 shows the temporal profile of PM10 concentrations measured at Gruvebadet with daily resolution [89] in the period from 28 July to 24 September 2017. While, a number of peaks can be observed in this period, none, especially that of 20-24 August, appears to be influenced by BB aerosol. In fact, the largest atmospheric load was measured between 5-9 September and 13-17 September and peaks of similar intensity can be also seen in previous and following days making this parameter not discriminating for possible BB events. However, these results suggest to investigate the aerosol chemical composition and the vertical aerosol structure to understand both the presence of chemical tracers of BB and the location of maximum concentrations of the transported plumes (see Section 3.3).
Focusing on the chemical composition, additional information appears. In this respect, Figure 4 reports the temporal behavior of three chemical species that can be used as markers of BB (nitrate and oxalate) and of transport from anthropic areas (non-sea salt sulphate, nssSO4 2-). They were analyzed both in bulk PM10 samples and size-segregated aerosols (<1, 10-2.5, and 2.5-1 µm; Section 2.2). Both nitrate and oxalate exhibit a peak in the 21-24 August period, in all the size fractions as well as in bulk PM10. This represents the second largest peak in August-September, with 6-9 September being characterized by the highest concentration of both markers in each class, except for the sub These BB events generated huge quantities of aerosols, whose atmospheric transport was observed between 28 July and 26 August 2017 in the Arctic (Sections 3.2 and 3.3). According to the NASA database, the Visible Infrared Imaging Radiometer Suite (VIIRS) on the Suomi NPP satellite captured particularly heavy smoke originating in northern Canada as of August 15 and 16, 2017, which was transported to the Arctic over Greenland (Figure 2).
In order to confirm the source regions of the aerosol particles transported over Svalbard during the BB events, air mass trajectories were calculated at altitudes of 500, 1000, and 1500 m a.s.l., using the FLEXPART model ( Figure S1).
The trajectories have been calculated for a period between 18 and 26 August 2017. The results are presented in Figure S1 and they demonstrate that in August 2017, Greenland and Canada wildfires were the sources of aerosols whose evolution and properties are discussed hereinafter.

Chemical Speciation and Optical Properties at Ground-Level (Svalbard)
During the entire August 2017, an ensemble of aerosol in-situ chemical and optical measurements were carried out at the Gruvebadet station in Ny-Ålesund. These measurements are representative of the portion of the atmospheric aerosols that can reach the surface as demonstrated by vertical aerosol profiles measured at Gruvebadet station using a tethered balloon system since 2011 [8,37,[86][87][88]. Atmospheric PM10 concentrations and the related chemical composition, measured at Gruvebadet are presented in Figures 3 and 4. Figure 4 shows the temporal profile of PM 10 concentrations measured at Gruvebadet with daily resolution [89] in the period from 28 July to 24 September 2017. While, a number of peaks can be observed in this period, none, especially that of 20-24 August, appears to be influenced by BB aerosol. In fact, the largest atmospheric load was measured between 5-9 September and 13-17 September and peaks of similar intensity can be also seen in previous and following days making this parameter not discriminating for possible BB events. However, these results suggest to investigate the aerosol Atmosphere 2019, 10, x FOR PEER REVIEW 9 of 27 micrometric one, where the two peaks (20-24 August and 5-9 September) are comparable in intensity. Since the 21-24 August peak does not appear in the nssSO4 2-profile (while the 6-9 September increase is clearly visible), it is likely that the two concentration enhancements are related to two different transport events in terms of source region. Indeed, nssSO4 2-usually marks "haze" from lower latitude areas characterized by a higher anthropic impact while nitrates and oxalates are more closely linked to combustion processes and related chemical reactions during long-range transport [40,56,90].   Atmosphere 2019, 10, x FOR PEER REVIEW 9 of 27 micrometric one, where the two peaks (20-24 August and 5-9 September) are comparable in intensity. Since the 21-24 August peak does not appear in the nssSO4 2-profile (while the 6-9 September increase is clearly visible), it is likely that the two concentration enhancements are related to two different transport events in terms of source region. Indeed, nssSO4 2-usually marks "haze" from lower latitude areas characterized by a higher anthropic impact while nitrates and oxalates are more closely linked to combustion processes and related chemical reactions during long-range transport [40,56,90].   Focusing on the chemical composition, additional information appears. In this respect, Figure 4 reports the temporal behavior of three chemical species that can be used as markers of BB (nitrate and oxalate) and of transport from anthropic areas (non-sea salt sulphate, nssSO 4 2-). They were analyzed both in bulk PM 10 samples and size-segregated aerosols (<1, 10-2.5, and 2.5-1 µm; Section 2.2). Both nitrate and oxalate exhibit a peak in the 21-24 August period, in all the size fractions as well as in bulk PM 10 . This represents the second largest peak in August-September, with 6-9 September being characterized by the highest concentration of both markers in each class, except for the sub micrometric one, where the two peaks (20-24 August and 5-9 September) are comparable in intensity.
Since the 21-24 August peak does not appear in the nssSO 4 2profile (while the 6-9 September increase is clearly visible), it is likely that the two concentration enhancements are related to two different transport events in terms of source region. Indeed, nssSO 4 2usually marks "haze" from lower latitude areas characterized by a higher anthropic impact while nitrates and oxalates are more closely linked to combustion processes and related chemical reactions during long-range transport [40,56,90]. Ratios between average concentration of each determined chemical species and PM 10 load in bulk and size-segregated aerosol during the event are presented in Figure 5 below.
In order to highlight possible differences in chemical composition of aerosols collected during the event (20-23 August 2017) and in the previous and following days, here we show the ratio between corresponding averages for each of the measured species. Figure 5 shows that, among all of the measured chemical markers, the ones more related to biomass burning processes show the highest ratios in the period 20-23 August, when compared with the rest of August. In particular, both nitrate and oxalate show the largest enhancements, with ratios above 2 for the submicrometric fraction and bulk PM 10 and ratios above 1.5 for 1-2.5 and 2.5-10 µm for oxalate. Nitrate alone exhibit marked increases in all size classes. Interestingly, other components that can be associated with combustion processes as an intermediate product, acetates, and formate, show substantial enhancements in selected size classes. In particular formate shows an increase larger than 100% in the submicron size range and acetate show very high enhancements in two size classes, <1 µm (increase of about 350%) and 2.5-10 µm (increase of about 500%).
Among the analyzed ionic markers, K+ is often used as a proxy of biomass burning processes, and it actually shows an increase during the event, with respect to the rest of the month, in all the size classes, except for the submicrometric one. Except for this size fraction in fact, K+ concentration is more than 50% higher than the average of the rest of the month. Nevertheless, it is likely that in this case such an enhancement is due to the arrival of marine air masses, also delivering the products of the biomass burning event. Indeed, together with the markers of biomass burning, also typically primary marine components (Na + , Cl -, Mg 2+ , K + ) show notable increases, with ratios above 1.5 and in a few cases also above 2 (Na + and Clin submicrometric fraction, Na + and Mg 2+ in bulk PM 10 ). This support the hypothesis that air masses supplying the products of biomass burning probably travelled at low altitudes over the ocean, for a relatively long time. It is likely that air masses delivering the products of biomass burning were enriched in primary marine aerosols, probably travelling at low altitudes over the ocean, for a relatively long time as also demonstrated by back trajectories reported in Figure 4.
As mentioned above, also PM 10 load exhibits an increase during the event, with a ratio above 1.5. By this general look at the chemical composition "inside" and "outside" the event, we can conclude that aerosol chemistry studies in Ny-Ålesund can provide information on biomass burning events through concentration peaks of specific markers, especially nitrate and oxalate. Nevertheless, the chemical signature for this event cannot be completely defined at this stage, showing also a role of acetate and formate, in selected size classes. These results suggest that the biomass burning aerosols are quite inhomogeneous in terms of chemical composition. Ratios between average concentration of each determined chemical species and PM10 load in bulk and size-segregated aerosol during the event are presented in Figure 5 below. In order to highlight possible differences in chemical composition of aerosols collected during the event (20-23 August 2017) and in the previous and following days, here we show the ratio between corresponding averages for each of the measured species. All of the aforementioned data highlighted the impact of BB at ground level. Simultaneously, aerosol optical properties (i.e., scattering and absorption coefficients) were measured at ground in the same site. Their temporal behavior for the entire for the period between 27 July and 24 September 2017 is reported in Figure 6. It is particularly evident that the scattering and absorption coefficients were quite stable till the beginning of September (around 2 ± 1 Mm -1 for b s and lower than 0.05 Mm -1 for b a ). Only narrow peaks of b a were detected during August, with the exception of the 10-12 August period, when it reached values above 0.4 Mm -1 for a relatively long period (almost one day). This corresponds well with the mid-August period when over 1000 separate fires had been recorded in Canada and the burned area grew to almost 9000 km 2 (Section 3.1). Moreover, a second significant absorption peak (again up to 0.6 Mm −1 ) was recorded in the period between 21 and 23 August, which is in agreement with the Canadian fire development when a large number of independent fires merged into a single larger wildfire covering an area of 4674 km 2 in size (Section 3.1). During neither of these episodes, the increasing of the aerosol absorption was concurrent with an increase in the scattering, indicating an increase of the fraction of radiation-absorbing particles, like BB ones. This is confirmed also by the SSA values, dropping down to values well below 0.9, from the average background value above 0.95. After the beginning of September, very high values of b s (up to 15 Mm -1 on 14-15 September) and b a (up to 0.4 Mm -1 for more than one day on the 5-6 September) were measured. Despite these high values, only the first one (5-6 September) should be partly associated with a BB event, showing a SSA value slightly below the usual one, while during the second period (14-15 September) SSA remained stable above 0.95. On the overall considered period, AAE didn't give meaningful information on the relative contribution of different types carbonaceous particles, presenting very scattered values not correlated with the BB intrusions. The ground-based data (chemical and optical) have been impacted by BB aerosols even though Ny-Å lesund is often affected by local micrometeorology [91]. It has been previously demonstrated that high altitude long range transport events may be driven to the ground trough atmospheric boundary layer dynamics [40]. However, as reported by Warneke et al. [92], layers of biomass burning aerosols are often present above the top of the inversion layer. Due to the magnitude of the Canadian and Greenland wildfires, the ground-based values suggested to also investigate the presence of aerosol layers aloft in the atmosphere, and the results are discussed below.

Optical Impact in the Arctic: Vertical and Columnar Measurements
In order to highlight the impact of BB events in the Arctic, summer background aerosol conditions are recalled first as a base-line level for comparison. Svalbard summers are usually characterized by long periods of high cloudiness and fairly common precipitation episodes [93].
Usual summer values of AOD should not exceed 0.1 at 500 nm [47]. Such situation was observed at the beginning of summer and until early August 2017 in the investigated region. Particularly, average The ground-based data (chemical and optical) have been impacted by BB aerosols even though Ny-Ålesund is often affected by local micrometeorology [91]. It has been previously demonstrated that high altitude long range transport events may be driven to the ground trough atmospheric boundary layer dynamics [40]. However, as reported by Warneke et al. [92], layers of biomass burning aerosols are often present above the top of the inversion layer. Due to the magnitude of the Canadian and Greenland wildfires, the ground-based values suggested to also investigate the presence of aerosol layers aloft in the atmosphere, and the results are discussed below.

Optical Impact in the Arctic: Vertical and Columnar Measurements
In order to highlight the impact of BB events in the Arctic, summer background aerosol conditions are recalled first as a base-line level for comparison. Svalbard summers are usually characterized by long periods of high cloudiness and fairly common precipitation episodes [93]. Usual summer values of AOD should not exceed 0.1 at 500 nm [47]. Such situation was observed at the beginning of summer and until early August 2017 in the investigated region. Particularly, average daily AOD values (500 nm) measured onboard r/v Oceania between June and early August 2017 varied from 0.032 to 0.067. Unfortunately, due to long cloudiness periods over the ocean, only very few data could be acquired from sun photometer measurements from board of r/v Oceania in the study period. Lidar based measurements, carried out using the AWI KARL lidar (Ny-Ålesund-June 2017) also showed a small presence of aerosol particles in the air above Svalbard ( Figure S2) confirming the representativeness of the investigated Arctic summer prior to the BB events (details in Supplementary Material).
Moving from the background conditions prior the Canadian/Greenland fires towards August (when the BB events were measured at ground, Section 3.2) both the columnar values of AOD and the vertical structure of aerosols appeared seriously impacted.
Looking at the distribution (both temporal and spatial) of AOD at 500 nm based on AERONET data it is clear what magnitude the biomass burning had in summer 2017 in the Arctic region. Table 1 shows AOD data for the investigated period, available from AERONET. All the reported values exceeded 0.15 and reached values up to 2.1 near the source region at Churchill underlying the spatial impact of such wildfires. The AERONET based AODs for the period of early August up to the end of the month show that the AODs at stations located close to or in Canada present very high daily average values and some are clearly a few days ahead of the event observation in the Svalbard region. Increased values of AOD in Ny-Ålesund at all wavelengths started on 18/19 August 2017. Elevated values of AODs are then registered in Sodankylä and Andenes a few days after the Svalbard observation of the event. An interesting case is attributed to 19 August over r/v Oceania (not in the table), who was west of Spitsbergen. Very low values of AOD were registered using a Microtops II sun photometer during the entire day and the daily average AOD at 500 nm amounted to 0.075. The smoke plume was just outside the range of Oceania measurements while at the same time the AOD values at 500 nm in other locations showed significant increase starting 19 August.
Focusing on Ny-Ålesund and Hornsund AOD values, during transport events, two outbreaks occurred in August over the low aerosol background typical of the Arctic area (<0.05 at 550 nm). The first pronounced peak occurred between 20 and 22 August. The AOD values (at 500 nm) recorded in Hornsund exhibit a sharp increase on 21 August (Figure 7a), reaching the maximum close to 0.4. Unfortunately, cloudy conditions occurred during the BB event in Ny-Ålesund and thus only few measurements were made at that period (Figure 7b). However, they reached values up to 0.5 suggesting that the BB plume affected the whole archipelago. Later, in the middle of September (between 14 and 17), there is an increase in AOD levels at both sites. In this case, data from Hornsund are cloud contaminated until 17 September (Figure 7c), when the AOD (at 500 nm) is beyond 0.2. In Ny-Ålesund (see Figure 7d), there is a value increase on 14 September but clouds also appeared again. These later September peaks relate to the last BB events in Canada.   The analysis of the Ångström Exponent together with AOD suggested the presence of very highly concentrated, tiny, particles in the atmosphere. Figure 8    In addition, looking at the columnar size distribution of aerosol particles measured over Hornsund during the August event (Figure 9a), aerosol volume increased in the fine fraction along the day while the coarse mode was unaffected, which is in agreement with a BB signature, especially because the volume mode radius shifted from 0.15 to 0.26 µm along with aerosol size increasing due to both aging and condensation of organic material. Moreover, single scattering albedo (SSA; Figure 9b) decreased, especially in the infrared region (from 0.75 to 0.7) very sensible to an increase of the absorbing black carbon fraction generated by the event. The absolute values of SSA at different wavelengths were very low for the Arctic (below 0.85) stressing the absorbing nature of aerosols loaded along the atmospheric column and in agreement with ground-based observations reported in Figure 6. Finally, the asymmetry parameter ( Figure 9c) increased along the day describing a situation in which the forward scattering was slightly enhanced by the volume mode radius increasing (Figure 9a).
Due to the weather related relatively small number of sun photometer data, we have run both NAAPS and CAMS models to obtain results to describe spatial distribution of the BB plume over the Arctic. Figure 10 reports simulations from the NAAPS aerosol transport model for selected days in August 2017. The results confirm that the sources of aerosols detected over Svalbard were located in Canada and Greenland but, most important, the BB smoke plume (blue color) was spreading from over Canada and/or Greenland towards Svalbard region on each of the investigated days. Air mass transport from Canada and Greenland towards Svalbard was observed during the entire BB event in July and August 2017 at all altitudes a.s.l. (Figure 10). Moreover, the CAMS model provided optical property results for the entire August and for all locations. The results clearly show a plume of aerosols generated around Canada and spreading towards the Arctic region, including Svalbard on 18 August 2017 ( Figure 11).
Values of aerosol optical depth presented in Figure 12 and calculated for Ny-Ålesund and Hornsund between 16 and 27 August 2017 (550 nm) show relatively low values between 17 and 19 August (below 0.1) at both locations and then a dramatic increase between 19/20 and 26 August 2017 (up to 1). There is a distinct shift in time, which is obvious due to the location of both stations, Ny-Ålesund is further north, thus an aerosol event started there earlier than in the southern station of Hornsund (Figure 12). The AOD values reached their peak at c. 0.9 (550 nm) on 22/23 August (Figure 12). Comparison with chemical components from the CAMS model shows that during events parameters connected to the BB episodes, especially BC and OM were showing increasingly higher values simultaneously with increasing AOD, while non-anthropogenic component-SS gets lower.
The dust component associated with small particles confirms changes in chemistry of atmosphere during the BB event. Domination of OM component (up to 65%) is simultaneous with decreasing participation of SS between 22 and 27 August (1%-3%) at both locations (Figures 13 and 14).
Hornsund during the August event (Figure 9a), aerosol volume increased in the fine fraction along the day while the coarse mode was unaffected, which is in agreement with a BB signature, especially because the volume mode radius shifted from 0.15 to 0.26 µ m along with aerosol size increasing due to both aging and condensation of organic material. Moreover, single scattering albedo (SSA; Figure  9b) decreased, especially in the infrared region (from 0.75 to 0.7) very sensible to an increase of the absorbing black carbon fraction generated by the event. The absolute values of SSA at different wavelengths were very low for the Arctic (below 0.85) stressing the absorbing nature of aerosols loaded along the atmospheric column and in agreement with ground-based observations reported in Figure 6. Finally, the asymmetry parameter (Figure 9c) increased along the day describing a situation in which the forward scattering was slightly enhanced by the volume mode radius increasing ( Figure  9a). Figure 9. AERONET based results of particle size distribution (a), single scattering albedo (b) and asymmetry parameter (c) obtained from Cimel sun photometer during the aerosol event over Hornsund.
Due to the weather related relatively small number of sun photometer data, we have run both NAAPS and CAMS models to obtain results to describe spatial distribution of the BB plume over the Arctic. Figure 10 reports simulations from the NAAPS aerosol transport model for selected days in August 2017. The results confirm that the sources of aerosols detected over Svalbard were located in Canada and Greenland but, most important, the BB smoke plume (blue color) was spreading from over Canada and/or Greenland towards Svalbard region on each of the investigated days. Air mass transport from Canada and Greenland towards Svalbard was observed during the entire BB event in July and August 2017 at all altitudes a.s.l. (Figure 10). Moreover, the CAMS model provided optical property results for the entire August and for all locations. The results clearly show a plume of aerosols generated around Canada and spreading towards the Arctic region, including Svalbard on 18 August 2017 ( Figure 11).
Values of aerosol optical depth presented in Figure 12 and calculated for Ny-Å lesund and Hornsund between 16    The CAMS model provides AOD values every three hours. As AERONET instantaneous data cannot match with that time resolution, the modelled values have been temporally interpolated with the two closest model points. A simple linear evolution is assumed in the 3 h step. In this way, all the experimental points recorded by AERONET can be used in the assessment of the model. This improves the statistics and the significance of the results obtained. To assess model calculations, the statistical methods explained by Willmot [95] are considered.     The European Arctic presents very low AOD values which can be a major source of discrepancies between model and experimental data. Conditions of AOD (550 nm) <0.10 are well characterized by the two techniques. When the aerosol load increases in experimental values model data increase too. In most of the cases, the slope of the increase is close to the 1:1 line. The largest discrepancy occurs in Hornsund when AERONET data notably increase until 0.3 meanwhile modelled values are still below 0.1.  The CAMS model provides AOD values every three hours. As AERONET instantaneous data cannot match with that time resolution, the modelled values have been temporally interpolated with the two closest model points. A simple linear evolution is assumed in the 3 h step. In this way, all the experimental points recorded by AERONET can be used in the assessment of the model. This improves the statistics and the significance of the results obtained. To assess model calculations, the statistical methods explained by Willmot [95] are considered.
The CAMS modeling AOD results with distinction to 5 aerosol chemical components are presented in Figures 13 and 14.
Results presented in Figure 14 clearly show an increasing role of organic matter in the composition of aerosols between 16 and 27 August in Hornsund. On 16 and 17 August sea salt was still a dominating component (33% to 42%) and starting on 18 August it was organic matter that took over in the aerosol composition and on 18 August its percentage was 37 and it reached its peak on 21 August with 65% which related to a very high AOD (500 nm) of 0.62. Also, the role of sulphate and black carbon was increasing starting 18 August and it peaked around 21 August. This indicates the presence of BB aerosols over the region. The same pattern is clear in Figure 14, where the same data are presented for Ny-Ålesund.  Figure 15. The value for AERONET at that wavelength is calculated using as input the measured AOD value at 500 nm and the Å ngström coefficient. The 1:1 line is also given in the plot. In general, the CAMS model compares quite well to the observational data if the sparse meteorological information in the Arctic is considered. Only for a few cases for Hornsund, the model underestimates the air pollution. The mean bias error (MBE) and mean absolute bias error (MABE) present values below 0.03 and 0.05, respectively. The root mean square error (RMSE) also shows low values below 0.06 at both sites. The index of agreement (d) between both variables is very high, >0.95. There is a good agreement between measurements and model using 1154 and 182 points in the evaluations for Hornsund and Ny-Å lesund, respectively. Further evaluations with a larger number of points in a longer time should improve the statistics of this comparison.
Moving from the columnar properties to the vertical aerosol structure, during the BB event, ceilometer (Lufft CHM-15k) measurements were made between 19 and 21 August 2017. Results are presented in Figure 16. The plot shows multiple aerosol layers at altitudes up to 8 km a.s.l. These layers are mostly pronounced on 20/21 August, which is in agreement with AOD values measured by sun photometers in Hornsund.  Figure 15. The value for AERONET at that wavelength is calculated using as input the measured AOD value at 500 nm and the Ångström coefficient. The 1:1 line is also given in the plot. In general, the CAMS model compares quite well to the observational data if the sparse meteorological information in the Arctic is considered. Only for a few cases for Hornsund, the model underestimates the air pollution. The mean bias error (MBE) and mean absolute bias error (MABE) present values below 0.03 and 0.05, respectively. The root mean square error (RMSE) also shows low values below 0.06 at both sites. The index of agreement (d) between both variables is very high, >0.95. There is a good agreement between measurements and model using 1154 and 182 points in the evaluations for Hornsund and Ny-Ålesund, respectively. Further evaluations with a larger number of points in a longer time should improve the statistics of this comparison. This result is very important when considering that the same type of aerosols can produce different climatic effects (warming or cooling) and local feedbacks (snow/ice albedo, clouds) depending on its vertical location [95][96][97][98][99][100][101]. For example, even counterintuitive, if BC is located in the free troposphere it may reduce the surface air temperature and promote the increase in the sea-ice fraction. However, it may potentially warm the Arctic if it is located immediately close to the ground, [98,100,101] as described in Section 3.2.

Discussion and Conclusions
In this paper, the authors present the results of aerosol studies over the European part of the Arctic, with special emphasis on the Svalbard region. The results refer to the summer 2017 (July-September) biomass burning episode which took place in Canada and Greenland. Special attention is given to the middle August, when the episode had the greatest strength, however, we provide some hints that the September part of the episode could have been also very significant.
The peak of the BB event, in case of the Svalbard region, took place around 18-21 August 2017. The analyses of air mass back trajectories show that in this period air masses were blown to Svalbard from over Canada and Greenland. This is also supported by the analyses made with the NAAPS aerosol transport model.
Measurements made with an ensemble of various instruments reveled that this aerosol episode was one of the worst cases in last decade [45,46,[51][52][53][54]. Daily averages of level 2.0 aerosol optical Moving from the columnar properties to the vertical aerosol structure, during the BB event, ceilometer (Lufft CHM-15k) measurements were made between 19 and 21 August 2017. Results are presented in Figure 16. The plot shows multiple aerosol layers at altitudes up to 8 km a.s.l. These layers are mostly pronounced on 20/21 August, which is in agreement with AOD values measured by sun photometers in Hornsund. This result is very important when considering that the same type of aerosols can produce different climatic effects (warming or cooling) and local feedbacks (snow/ice albedo, clouds) depending on its vertical location [95][96][97][98][99][100][101]. For example, even counterintuitive, if BC is located in the free troposphere it may reduce the surface air temperature and promote the increase in the sea-ice fraction. However, it may potentially warm the Arctic if it is located immediately close to the ground, [98,100,101] as described in Section 3.2.

Discussion and Conclusions
In this paper, the authors present the results of aerosol studies over the European part of the Arctic, with special emphasis on the Svalbard region. The results refer to the summer 2017 (July-September) biomass burning episode which took place in Canada and Greenland. Special attention is given to the middle August, when the episode had the greatest strength, however, we provide some hints that the September part of the episode could have been also very significant.
The peak of the BB event, in case of the Svalbard region, took place around 18-21 August 2017. The analyses of air mass back trajectories show that in this period air masses were blown to Svalbard from over Canada and Greenland. This is also supported by the analyses made with the NAAPS aerosol transport model.
Measurements made with an ensemble of various instruments reveled that this aerosol episode This result is very important when considering that the same type of aerosols can produce different climatic effects (warming or cooling) and local feedbacks (snow/ice albedo, clouds) depending on its vertical location [95][96][97][98][99][100][101]. For example, even counterintuitive, if BC is located in the free troposphere it may reduce the surface air temperature and promote the increase in the sea-ice fraction. However, it may potentially warm the Arctic if it is located immediately close to the ground, [98,100,101] as described in Section 3.2.

Discussion and Conclusions
In this paper, the authors present the results of aerosol studies over the European part of the Arctic, with special emphasis on the Svalbard region. The results refer to the summer 2017 (July-September) biomass burning episode which took place in Canada and Greenland. Special attention is given to the middle August, when the episode had the greatest strength, however, we provide some hints that the September part of the episode could have been also very significant.
The peak of the BB event, in case of the Svalbard region, took place around 18-21 August 2017. The analyses of air mass back trajectories show that in this period air masses were blown to Svalbard from over Canada and Greenland. This is also supported by the analyses made with the NAAPS aerosol transport model.
Measurements made with an ensemble of various instruments reveled that this aerosol episode was one of the worst cases in last decade [45,46,[51][52][53][54]. Daily averages of level 2.0 aerosol optical depth (500 nm) were as high as over 2.0 (AERONET Churchill station), with values high above 0.2 over Greenland (AERONET stations Narsarsuaq and Kangerlussuaq) while in Spitsbergen they reached values of 0.147 (AERONET Hornsund station). However, momentary values (500 nm) were much higher and could be exceeded e.g., 0.3 in Kangerlussuaq, 0.4 in Ny-Ålesund, or 0.3 in Hornsund. The Ångström Exponent values indicate presence of small particles, presumably of BB origin. This statement is also supported by the absorption, scattering, SSA, and asymmetry parameter analyses for the Spitsbergen studies. These results are also presented for a September period, however, due to not complete analyses the authors leave the potential discussion for a future paper.
The BB aerosol plumes were observed over Hornsund at altitudes up to 8 km a.s.l., with their peak on 20/21 August, which is in agreement with AOD values measured by sunphotometers.
Aerosol analyses of aerosols made at the Gruvebadet laboratory in Ny-Ålesund show that the scattering values are quite homogeneous during August 2017, however, absorption data show several significant peaks throughout August 2017, including the period between 19 and 23 August, associated with absorbing aerosols such as those from biomass burning.
Temporal profiles of PM 10 load at Gruvebadet with daily resolution in the period 28 July-24 September 2017 show a number of peaks. Even though an increase in concentration appears in 20-24 August, peaks of similar intensity can be seen also in previous and following days making this parameter not discriminating for possible biomass burning events. However, temporal behavior of two main chemical species that can be used as markers of biomass burning (nitrate and oxalate) show a peak between 21 and 24 August in all the size fractions and in bulk PM 10 .
Based on the analyses we realize that it is non-trivial to bring ground-based in situ and remote sensing data together, e.g., the pollution pathways in the boundary layer and in the free troposphere may differ. Therefore, we need dedicated campaigns to understand this. A ship, such as e.g., r/v Oceania in situ and with at least Microtops II measurements can be of great help since the boundary layer over the ocean is less disrupted than that in the fjords due to mountain orography, as it is e.g., in Ny-Ålesund.
This biomass burning event was among two worst aerosol events recently reported for Svalbard region (previous in 2015) and the climatic results need further investigations. However, the authors, who were also involved in the 2015 study and previous events, see that the summer biomass burning related aerosol episodes over Svalbard may be becoming a serious threat to the radiative balance of the region, hence may have a serious impact on the strength of the so-called Arctic amplification effect [47].