Review of Black Carbon in the Arctic — Origin , Measurement Methods , and Observations

Current knowledge about black carbon (BC) emission estimates, state-of-the-art measurement methods, near-surface BC concentrations ([BC]), and mixing ratios in snow is consolidated for the Arctic. Since no direct method exists to measure [BC], results from modern indirect methods differ among devices. Pan-Arctic wide [BC] and changes are hard to access; monitoring often ends once national ambient air quality standards are met. Few remote sites have long records. Past measurements showed distinct differences among the various Arctic climate regions. Past and own observations in communities permit qualitative discussion of the diurnal course, response to weather, season, or different emission situations like weekdays and weekends at a given site and/or among sites. Comparison of data from collocated aethalometer indicated more similar accuracy than found in midand low-latitudes despite of much lower ambient temperatures and [BC]. Snow samples give an incomplete glimpse at the removal and input into ecosystems.


Introduction
In this Millennium, black carbon (BC) in the Arctic has gained huge public interest [1].In general, BC results from incomplete combustion of carbonaceous fuels.Different fuel types produce different amounts of BC.Natural gas and diesel fuels are on the low and high amount side, respectively.Bottom-up inventory methods estimated the 2000 total global BC emissions to be 7500 Gg•yr −1 with an uncertainty of 2000 -29,000 Gg•yr −1 [2].
The major concern of the public is the impacts of BC on the local radiation budget.Numerous studies discussed its direct and indirect impacts on the radiation budget in depth [2] [3] [4] [5].In the interest of brevity, we renounce repeating them here.
BC is also long known to affect the water cycle; while BC alone is insoluble in water, aged BC aerosols can act as cloud and ice nuclei [6].Consequently, cloud and precipitation formation as well as atmospheric removal processes might change with changing BC concentrations [7].Changes in cloudiness, cloud-life time, coverage and precipitation again impact the energy cycle via changes in the radiation budget (e.g.[6] [7] and work cited therein).
An aspect of BC, which the public is often unaware of, is its role in the formation of particulate matter (PM), which at sizes less than 2.5 µm in diameter (PM 2.5 ) is health adverse (e.g.[8] [9] [10] [11] [12]).PM 2.5 can enter the lungs, pass thru human skin, and may cause cardiovascular problems, respiratory diseases, and even cancer.For a comprehensive review on the health adverse effects of short and long-term exposure to PM 2.5 see [13] [14]; see [15] for a BC impact study on health with focus on the Sub-Arctic.
PM is a composite of various species which can include BC, elemental carbon (EC) and organic carbon (OC), among other things.BC and EC themselves have negligible impacts on human health.However, upon emission, BC particles undergo complex processes with co-emitted OC, toxic semi-volatile organics, inorganic salts from both nitrogen dioxide (NO 2 ) and sulfur dioxide (SO 2 ), as well as atmospheric water.These species may deposit on hygroscopic carbonaceous particles; also other aerosol components can join BC during particle growth [6].This so-called aging process means that a particle rarely is pure BC, and BC containing aerosols may carry toxic and/or health adverse constituents.
Ambient air samples of BC-aerosol contain up to more than 90% other constituents.Thus, BC mass concentration (in ng•m −3 ) cannot be measured directly by weighing the BC in an air sample; BC mass must be determined indirectly by measuring optical, and/or thermal properties of BC.For instance, BC has high sp2-bonded carbon content and Raman spectroscopic responses similar to graphite.The physical relationship between the measured property and the amount of BC then permits deriving the BC concentration, hereafter written as [BC].Typical atmospheric residence times of BC-containing PM are 4 to 7 days depending on weather conditions; during transport, atmospheric species from emission sources along the way may attach to BC particles [6].Transport of BC emitted in mid latitudes contributes a substantial fraction to Arctic BC [16] [17] [18] [19].
Recent research showed that air quality in Arctic cities differs not only due to their emissions, size, and the advection of pollutants from downwind regions, but also due to their local climate.Cities located in temperate mesothermal climate according to the Köppen-Geiger classification [20] typically have better air quality than cities in continental microthermal climates [21].
Given the huge variability of BC-containing aerosol, and the fact that the Sweden, for instance, measures in stacks, while Denmark, Norway, and Finland apply emission factors based on dilution tunnels.This means available emissions measurements are hard to compare.See [40] [41] for reviews on US and European BC source apportionment studies, and [37] for BC emissions in Russia.
Many emissions inventories only report PM 10 and PM 2.5 without specification of constituents.Out of the inventories reporting carbon emissions, only few distinguish between EC, OC and BC.Even then, emission estimates differ among inventories due to the assumed partitioning of PM, assumed emission factors, and whether these data are activity based or from reported data.Emissions inventories differ due to the method of creating them (e.g.bottom-up, top-down), their spatial and temporal resolution, and the base-years used [46] [47] [48].For differences between bottom-up and top-down emissions inventories see [47] [49] [50].
Figure 1  2) differences in accounting for sources as point-, area-, or line-sources; and/or 3) lumping of emission sectors.Some inventories lump, for instance, residential and commercial space heating, while others lump residential space heating and power generation as one emissions sector.
Even within an emissions inventory spatial discrepancies of emissions from the same sector may occur due to data reporting requirements prescribed by regional and/or national agencies, and how they allocate (point-, area-, or line-sources) and/or lump data.Figure 2   Accounting for emissions over the region instead of at their actual area of occurrence leads to too high emissions everywhere, but too low emissions where they actually occur.Concentration ratios of competing reactants affect various chemical reactions.Thus, reaction products will differ notably if, as in the Besides fuel-type, ship emissions depend on speeds, engine load and size.These factors cause uncertainty in estimated BC emissions from Arctic shipping [51] [52] especially when ships cruise in an area with partial sea-ice occurrence and/or ice-breaking mode.The use of automatic identification system marine traffic data (position, speed, ship characteristics) could improve accuracy of BC ship emissions [52 The same source lists the Norwegian total annual BC emissions at 5.1 Gg in 2011.Wood stoves are the primary device for space heating with 1.581 Gg.National shipping, diesel vehicles, and the transport sector contributed about 0.867, 2.04, and 3.06 Gg.In the Norwegian Barents region, BC emissions were about 0.408 Gg [53].
In the Finish Barents region, BC emissions were about 0.77 Gg.The total annual mean Finish BC emissions were 6.46 Gg.Besides boilers and stoves, masonry heaters and sauna stoves are common in Finland [53].
In Canada, BC emissions are calculated using PM 2.5 emissions and BC/PM 2.5 ratios.In 2015, they were assessed as 38 Gg [54].
Many emission inventories lack BC emissions from flaring.These emissions depend on the gas flared (e.g.ethane, propane), its fraction of heavy hydrocarbons, its purity, flare size, and velocity across the flare tip [37]
[BC] are operationally defined as the light absorbing and thermal refractory fraction of carbonaceous aerosol, respectively [57].Health and atmospheric scientists applied aethalometers, particle soot absorption photometers (PSAP), single particle soot photometers (SPS), Continuous Light Absorption Photometers (CLAP), Integrating Spheres (IS), Hybrid Integrating Plate Systems (HIPS), Multi-Angle Absorption Photometers (MAAP), continuous soot-monitoring systems (COSMOS), thermal desorption techniques, and acoustic detection.Not all of these devices have been applied in the Arctic.Data from these instruments are hard to compare.Results from thermal desorption techniques depend on the applied protocol.All filter-based instruments provide time-integrated values for the chosen sampling interval, which may differ among sites.Unit-to-unit variability can reach up to 30% for PSAPs and aethalometer, and is lowest (<5%) for MAAP; PSAP and MAAP have lower instrument noise than aethalometer.Sources of systematic errors differ among instrumental designs; inlets may cutoff particle size at different diameters [58].
Light scattering by accumulated particles affects all filter-based absorption photometers (e.g.aethalometers, PSAP, CLAP, SPS, HIPS, MAAP).These effects are (incorrectly) termed apparent absorption.Their impact on light transmission causes overestimated absorption.Often empirical factors serve to correct for these effects reducing uncertainty to 20% -30% for the PSAP, IS, HIPS, and aethalometers, and 12% for the MAAP [59] in labs and/or outside the Arctic.Uncertainty also exists due to electronic noise, instrument variability, and calibration.Optical interaction between particles and filter matrix may yield positive bias for uncorrected filter-based measurements; otherwise, bias can be positive or negative [2].
Brown carbon and dust influence absorption at the same wavelengths as BC.Aging of BC alters the absorption properties of BC-containing aerosol.Increasing coating of BC enhances particle size and causes uncertainty in the mass attenuation cross section (MAC) [16].The bias in optical measurements of [BC], equivalent black carbon mass concentrations ([EBC]) and [EC] caused by coating material depends on the wavelength used; typically 630 nm is least affected except for high dust concentrations [2].Unfortunately, dust concentrations are often high in the North American Arctic due to unpaved roads and non-managed rivers [21].
In the following, we elucidate the principles of 1) filter-based absorption instruments exemplarily by the aethalometer, briefly describe; 2) thermal desorption techniques; 3) acoustic detection; 4) summarize results from instrument comparisons; and 5) evaluate the suitability of these methods for Arctic applications.

Aethalometer
Aethalometer use optical absorption for in-situ (online) EBC detection.They The sample attenuation is where I s and I r are the light transmitted through the sample and blank reference, respectively.Their ratio is the filter transmittance, τ .The temporal change in attenuation is proportional to the attenuation coefficient Here, A, Q, and ∆t are the filter spot size area (m 2 ), sample flow rate (m 3 •s −1 ), and time between light intensity measurements (collection time).The attenuation coefficient depends on the cross-section of light interactions with the collected particles.
The filter spot must advance to a clean filter spot after 2 to 5 min of sampling to avoid that transmission drops too much.This time span is the inner data processing cycle.The internal software provides a temporal mean at the end of this cycle.When collection time is set to values above an aethalometer's upper time limit of the inner cycle, sampling occurs on various spots; data are averaged at the end of the set collection time.This so-called boxcar averaging increases the detection limit with one over the square root of the sampling time at the cost of temporal resolution.Best practice is to sample at the upper limit, and perform post-processing offline as needed for improved signal-to-noise ratio.[EBC] maximum, mean, standard deviation, skewness and kurtosis were 65.9 ng•m −3 , 6.0 ± 6.2 ng•m −3 , 2.1, and 7.0 at 370 nm, and 137.3 ng•m −3 , 12.3 ± 12.5 ng•m −3 , 2.2, and 8.1, at 880 nm.Data correlation between the two channels was 88% (0.936).The much higher kurtosis than 3 indicates that the size distribution was non-Gaussian and outliers occurred more frequently for high concentrations than in a normal distribution.The spread in the scatter plot of measurements at 370 nm and 880 nm (Figure 3) indicates the presence of cofounding aerosol at increasing concentrations.Besides wavelength, the specific absorption coefficient depends, among other things, on season, filter interaction, environment, and filter loading effects related to particle aging.Values of ap σ range from 5 to 20 m 2 •g −1 , but are constant for a given atmospheric environment [60].During the 1989 to 2007 [EBC] measurements at Alert, Canada (82.45083N, 62.34167W), for instance, the manufacturer-recommended conversion factor of ap σ = 19 m 2 •g −1 was applied [61].Multiplication of aethalometer data with the applied conversion factor and dividing by a conversion factor X permits using X instead.See [62] [63] for comprehensive reviews of absorption coefficient measurements.

Quality Assurance/Quality Control
Even though all scientists apply quality assurance/quality control (QA/QC) and flag data following their national agency's requirements, no Pan-Arctic standardized method exists for flagging BC data.The same is true for corrections related from scattering.Particles in the upper layers of a sample scatter the light beam.
The reduced beam diminishes absorption in deeper layers, and [EBC] is underestimated.The manufacturer recommended attenuation coefficient, ap σ accounts for this shadowing effect at low to moderate filter loads.At high filter loads, optical measurements of light absorption require additional corrections [61].means the filter-loading effect can provide a proxy for particle coating, and differentiation between local/fresh and transported/aged particles [64].
Various empirical correction methods exist to calculate the absorption coefficient, ap σ from the attenuation coefficient, ATN σ [65].A simple approach derived from observations at three Finnish sites is [65].
Unfortunately, k-values vary seasonally because of higher EBC to total aerosol ratios in winter than summer; in summer, k depends on wavelength.Correction factors also differ with filter type.At Kevo, Finland (69.75N, 27.033E), for instance, k = 0.025 and k = 0.0021 were found for Whatman-42 and glass-fiber filters, respectively [61].In the European Arctic, annual shadowing correction is around 20% [61].
In laboratory studies, aethalometers have shown notable drift on timescales from hours to days.Unfortunately, magnitudes of drift differ among instruments.To determine instrument drift in the field, an absolute filter could serve to capture pressure changes in the sampling line relative to changes in ambient pressure, and relative humidity on a pristine filter.This method might be applicable also to loaded filters of various loadings.Drift testing for 24 h or more on a pristine filter minimizes influences from semi-volatile species; the absolute filter may itself impact the sample.Uncertainty due to in-field drift found for five instruments in the Arctic was about 0.01 -0.1 Mm −1 .
Measurements of [BC] are challenging.BC varies highly at high altitude as well as in coastal and Polar Regions.During clean periods, concentrations can be below the detection limit.Post-processing methods serve to increase the signal-to-noise ratio.They lower the detection limit at the cost of temporal resolution.Most often boxcar averaging [66] is applied.The averaging time varies among instruments, and sometimes with application.In aethalometer applications, for instance, the typical averaging time is 5 min.Some data-processing methods are instrument specific.To keep temporal resolution as high as possible, some scientists use the aethalometer equation and a temporal attenuation change threshold; post-processing only occurs when the threshold is exceeded [67].This method known as Optimized Noise-reduction Averaging (ONA) expands (reduces) averaging windows at low (high) filter loads.In contrast to the boxcar averaging, the temporal resolution at high aerosol concentrations is barely affected; temporal resolution is only low at small concentrations.ONA permits displaying notable trends at minimum distortion.
Applying an attenuation-threshold criterion, aethalometer were able to detect attenuation > 2.1 -6.7 Mm −1 at six Arctic sites in 2012 to 2014.Correlations of 5-min means at 370 and 880 nm were 88%, 89% and 88% over the entire period, on weekdays, and weekends, respectively (Table 1).
On week days, [EBC] were higher during the day and lower during the night than for these times on weekends (Figure 4).Overall, weekdays saw lower 5-min means than weekends.Differences between [EBC] 370nm and [EBC] 880nm were greater on weekdays, while variance was largest on weekends.These differences can be explained by the different traffic patterns and heating behaviors on weekdays and weekends.On weekends and on weekdays in the evening, households with furnace and woodstove prefer using the less expensive wood instead of heating fuel.The higher than 3 kurtosis (Table 1) indicated that the size distribution approached zero frequency slower than a Gaussian normal distribution meaning more outliers.The distribution was skewed towards few, but large BC particles.Skewness on weekdays and weekends barely differed indicating similar co-founders for the College site.
Similar was found, for instance, for September-October data collected at Tiksi (Russia 71.633N, 128.867E) in 2010 (Table 1).Despite Tiksi's population was only 5063 at that time, its [EBC] was about six times higher than at College.At Tiksi, ship emissions in port and along the shore of the Buor-Khaya Gulf of the Laptev Sea contribute to the [BC].

Other Photometers (PSAP, CLAP, COSMOS, MAAP)
Like the aethalometer, PSAP and CLAP instruments measure transmission of light through a light-diffusing filter while particles accumulate on the filter [68] at various wavelengths (467 to 660 nm).The optical designs are quite similar.
Silicon photodiodes measure the intensity of diffuse light transmitted through the sample spots (I s ) and clean (non-sampled) area (I r ) of the filter.
Both PSAP and CLAP typically use E70-2075W filters.Solenoid valves advance to the next sample spot once filter transmittance reaches 0.7.The PSAP uses a single sample spot on a 10-mm filter.The CLAP has two reference clean filter spots and eight sample filter spots on a 47-mm filter.This eight-spot design permits the CLAP to run at ideal conditions (τ < 0.7) eight times as long as the single-spot PSAP making the CLAP more suitable for remote Arctic sites.Typical flow rates and averaging times are 0.5 -1 liter•min −1 and 3 s −1 , respectively.The response depends on PM size and cross-sensitivity to particle scattering, which can be measured simultaneously by a nephelometer.
Redistributed liquid-like OC can cause bias as it affects light scattering and absorption.Bias in PSAP absorption coefficients can reach up to 50% -80% at high (15 -20) organic aerosol to BC ratios [2].Non-absorbing aerosols and suspended particles can yield overestimates of absorption at 550 nm of about 20% -30% that can be corrected for [69].
In contrast to PSAP and/or CLAP, COSMOS samples ambient air thru a heated inlet.The heat removes most of the volatile aerosol components prior to particle accumulation on the filter [70].Charring of low volatility organic species may cause bias in COSMO measurements.
The MAAP is the only real-time absorption photometer that corrects for artifacts from scattering by its design.It namely measures both the radiation transmitted through and scattered back from particles on a glass-fiber tape, and irradiation at various detection angles.Doing so, determines radiation fields in forward and backward direction and permits correction for enhanced absorption by filter loading, back-and multiple-scattering by PM and the filter matrix.A two-stream radiative transfer model calculates the absorption coefficient σ ATN .
MAAPs typically sample with a flow rate of 16.7 liter•min −1 .
Mass-specific absorption coefficient σ ap can be determined by dividing the MAAP absorption coefficient by the [EC] obtained from co-located thermal methods; at 550 nm and 670 nm, for instance, the recommended specific absorption coefficients are 8 m 2 •g −1 and 6.5 m 2 •g −1 , respectively [55].
The manufacturer specified MAAP minimum detection limit is <0.1 µg•m −3 for 2-min sampling times.Organic carbon can modify particle morphology, which causes bias.
Absolute uncertainty of PSAP and MAAP measurements are typically 0.02 and 0.06 Mm −1 for 5-min averages [71].The MAAP seems to be less prone to interpret light scattering ( sp σ ) as light absorption than a PSAP or aethalometer.

Thermal Techniques
The SPS uses laser-induced incandescence, i.e. visible thermal radiation to detect the refractory BC (rBC) mass of individual particles.An IR intra-cavity laser heats individual particles to vaporization temperature (about 400˚C).The SPS permits continuous, real-time observation of rBC over a wide range of concentrations.Calibration to ambient rBC material is the major cause of bias.Limitation of the size range may add uncertainty.In remote areas, the SPS captures the size range, i.e. most of the rBC mass and about 50% of the rBC number [2].The SPS may fail to detect particles less than 160 nm in diameter.
Several thermal desorption techniques exist.All of them are offline.They measure total OC, as they fail to distinguish between organic material (e.g.pollen, spores, yeasts) and primary emitted combustion-related OC aerosol.The most common thermal methods for BC-OC separation use that BC is non-volatile while OC is volatile.First, the sample filter is heated in helium gas to volatilize OC at temperatures ≤ 550˚C.Some protocols ask for a cooling time thereafter.Next the sample is heated in a 98% helium, 2% oxygen (O 2 ) mix to ≥ 550˚C.IR absorption or flame ionization serves to detect carbon dioxide or methane formed from the released gases [2].
Overestimates of [EC] may occur when some OC undergoes pyrolysis or charring, which convert OC to EC at high temperatures [2].Further uncertainty exists from 1) interpreting OC as EC, or vice versa, 2) sample compounds retarding OC volatilization, or 3) sample compounds facilitating EC release.Obviously, analytical results are an operational definition depending on the measurement protocol.Any data comparison would require using the same protocol at all sites.Unfortunately, used protocols differ by temperature ramping, correcting for OC charring during pyrolysis, and sample treatment.In the VDI 2465/1 protocol, for instance, filter samples reside for 24 h in a 1:1 mixture of toluene and 2-propanol to extract OC.After drying, filters are heated for 1 min at 200˚C and 7 min at 500˚C.Carbon evolving during this process counts as non-extractable OC.Combustion at 650˚C in an O 2 atmosphere provides the [EC].Minimum detectable carbon mass loading is 9 µg C. The scientific community has started a discussion on a standardized thermal-optical protocol for measuring OC and EC [72].
Some thermal methods monitor the optical reflectance of the sample filter [2] to correct for charring and pyrolysis effects.Thermal optical reflectance (TOR) or thermal optical transmittance (TOT) methods provide apparent elemental carbon (EC a ).They collect PM on quartz-fiber filters.Thermal-optical analysis monitors OC charring by the change in a laser signal either reflected from or transmitted through a filter punch.These methods define the carbon evolving after the detected optical signal attains the value it had prior to heating as EC.All other carbon count as OC [73].Uncertainty in thermal measurements may results from inaccurate correction for charred OC and catalytic oxidation of BC when the sample contains metals or metal oxides.
The integrating sphere method measures attenuation.A calibration curve established with commercial BC (Elftex 124, Cabot Corporation) permits conversion of the decrease in signal to [BC].The analysis requires homogeneous samples either from a rotating impactor, or suspensions in liquid.In the former case, the sample filter is placed at the entrance port between the light source and sphere.The sphere behaves as a diffuse detector.The change in transmission between the loaded and reference filters is interpreted as caused by absorption.
In the latter case, the loaded filter is dissolved in a suitable solvent (e.g.chloroform for polycarbonate filters).The refractive indices of the solvent and most aerosols are similar.Thus, non-BC aerosols cause no enhanced light absorption even when not fully dissolved.After placement of the dissolved sample in the center of the sphere, light from the source and the light uniformly distributed within the sphere irradiate the sample.Any absorbing substance in the sample decreases the signal.In this kind of spheres, the detector is at the bottom port and a second diffuse detector is at the entrance port of the light beam.This setup also permits analysis of samples on filters [74].Uncertainty results from the assumptions of no scattering losses and homogeneous optical filter properties.

Acoustic Detection Method
Photoacoustic spectrometers (PAS) draw particles into an acoustic cavity at typical flow rates of 1 liter•min −1 .Here a power-modulated laser irradiates them.
When particles absorb laser light, they transfer heat to the surrounding gas.Sensitive microphones detect the sound wave caused by heating and cooling.The intensity of the wave is interpreted as PM light absorption by calibration with NO 2 absorption.Typical averaging times are 3 to 4 s.Uncertainty results from gas-phase absorbers interfering with the BC detection.Overall uncertainty is about 5% [75].A recent study demonstrated that the sensitivity of PAS is too coarse to capture typical [BC] of the high Arctic.

Instrument Comparison and Cross-Calibration
Due to their design the various instruments may provide different concentrations.Recommendations for calibration and inter-comparison of filter-based visible light absorption instruments can be found in [69].Most instrument inter-comparisons took place in mid-latitude urban environments.A comparison of thermal and thermal-optical methods (TOT) using Sunset instruments, the IS method, and MAAP in Vienna, for instance, revealed 44% and 17% lower [EC] with the TOT than the mean of all observed [EC] including the TOT data; larg- est discrepancy occurred for large contributions of brown carbon to total carbon [76].Mean [BC] from filter and impactor samples collected for 3 weeks analyzed with a thermal method, a thermal-optical method, the VDI method, a filter transmission method and IS method agreed with the mean [BC] from co-located aethalometer and MAAP within their standard deviations [55].Co-located (sidewalk level in Vienna in June 1998) [BC] measurements revealed that the aethalometer provided only about 50% of the [BC] obtained with the IS method [55], i.e. a factor of 2 difference.
In urban environments, [BC] can be several orders of magnitude higher than in the Arctic.At low concentrations, a small difference among measurements can already mean a factor of 2 differences, while at high concentrations, a factor two difference requires a large difference in [BC] [77].This means at the typical Arctic [BC], distinguishing between signals, noise and inter-instrument uncertainty may become challenging.
[BC] determined by COSMOS, thermal-optical measurements, and SPS laser-induced incandescence agreed within 10% [70].CLAP and PSAP attenuation coefficients co-measured at 17 sites agreed within 8%, i.e. measurement uncertainty; CLAP and PSAP noise levels are similar within a factor of 2 [68].
Corrected TOR data provide the same OC/EC splits for most temperature protocols.For identical temperature protocol, simultaneous thermal/optical transmittance (TOT) corrections provide 30% lower [EC] than TOR; for protocols with high temperatures and short residence times, [EC] can be 70% -80% lower than TOR.This behavior occurs for samples dominated by anthropogenic combustion and wildfires [73].Aethalometer derived [BC] were only 30% of those obtained by an IS method; the IS method suggested a 21% overestimation compared to a thermal method [55].respectively.This means that these moments were 10%, 17%, 5% and 20% lower for 361 than 545.Measurements had a correlation coefficient of 0.909, i.e. 83% correlation.At 880 nm, correlation was 87% (R = 0.932).Mean, standard deviation, skewness, and kurtosis were 15.0 ± 16.1 ng•m −3 , 5.5, and 41.5 for aethalometer 545, and 13.9 ± 13.8 ng•m −3 , 5.5, and 49.1 for aethalometer 361; i.e. aethalometer 545 provided an 8% and 16% higher mean and standard deviation, a 15% lower kurtosis, but same skewness than aethalometer 361.Discrepancies may be due to local emissions, heterogeneous distribution of BC in air, discrepancies in the onset of a new measurement cycle, and aforementioned differences from one instrument to the other.Heterogeneous distributions may occur due to segregation effects [79].The skewness, and kurtosis values hint at discrepancies from size distribution.The tails approach zero frequency more slowly than a Gaussian, i.e. there are more outliers.Once a large particle enters one device, the other device cannot pick it up and large differences occur (Figure 5).

Normalization of Data from Various Sites
When applied to different sites, aethalometer-correction algorithms often fail to remove the ATN-dependence from shadowing effects [80].Values of ap σ are about constant in the same environment, but differ among regions.Aging of BC aerosol namely not only depends on co-emissions, but also on meteorological conditions due to temperature and/or moisture dependent reactions, transport and removal processes [6].Unfortunately, climatological conditions in the Arctic range from temperate maritime climate along the Norwegian coast, over warm summer continental (hemiboreal), sub-arctic and boreal climates to tundra climate along the North American and Siberian coasts of the Arctic Ocean [20].Consequently, even when measurement protocols, filters, sampling times, flow rates, and instrument types are the same, [BC] data from sites in different Arctic environments are not necessarily comparable due to site-dependent corrections [80].Using data from co-located filter-based absorption photometers at Alert (Canada), Summit, Barrow (Alaska), Tiksi, Pallas (Finland), and Zeppelin (Norway) yielded a normalization factor of 3.45 for Arctic aethalometer measurements at low elevation; this factor fails for high elevation (Summit).Typically, uncertainty of C f exceeds uncertainty from measurement noise.At Alert, for instance, [BC] was determined for 1989 to 2007 using an aethalometer cross-calibrated to a two-step thermal method.Consistent with common practice, an effective attenuation coefficient of 19 m 2 •g −1 and no loading corrections were used [55].At Barrow, and Ny-Ålesund, for instance, long-term PSAP observations of light absorption coefficient σ ATN exist.However, σ ATN has been converted to BC mass concentration ([BC] = σ ATN /σ ap ) using different σ ap values [81].During 2012-2015, COSMOS measured σ ATN at these sites.COSMOS derived [BC] (using a σ ap determined in previous studies) agreed within 9% with [EC] at Barrow over the 11-mon of PSAP and COSMOS co-location [81].PSAP light absorption coefficients exceeded those of COSMOS by 22% and 43%, at Barrow (PM 1 ) and Ny-Ålesund (PM 10 ), respectively [2].COSMOS σ ATN served to derive [BC] from the PSAP σ ATN measurements made since 1998 [81].These mixing ratios give a glimpse at the input into ecosystems after snowmelt from both sedimentation and removal by precipitation related processes.Some of the snow mass sublimates over winter.Thus, the mixing ratios found in snow samples at the end of winter, most likely exceed the ratios present in snowflakes at the time of snowfall.

Black Carbon Vertical Profiles
Like atmospheric [BC], BC in snow mixing ratios or [BC] in snow decreased over time.In Alaska, Canada, and on Arctic Ocean sea-ice, for instance, [BC] in snow decreased from 15 -30 ppb to 5 -10 ppb within 20 years.

Conclusions and Recommendations
Meteorological (e.g.temperature, humidity, inversions, stability, snow, storms, haze, prevailing wind systems) and oceanic conditions (e.g.open water vs. sea-ice) influence BC emissions, [BC], aerosol formation, accumulation, aging and removal.Since these conditions occur with regularity in a region, we conclude that [BC] mitigation measures in response to the anticipated increase of Arctic shipping have to be assessed with air quality models, which need spatial and temporal emission data as input.
Our study revealed that the actual amount of BC emitted in the Arctic was among minor reasons, uncertain due to 1) lack or inconsistent legislature for reporting requirements among Arctic countries, and even within countries at state/province levels; 2) no common standards for stack emissions measurements, 3) classification of emission sources, and 4) spatial allocation of emissions.Consequently, emission inventories using reported data show artifacts.We conclude that standardized definitions of emission sectors, spatial allocation, as well as emission measurement and reporting protocols are needed when air quality models are to be deployed in the search for Arctic-wide mitigation.Enforcement of Pan-Arctic wide consistent reporting would involve (time-consuming) negotiations among the eight Arctic countries.In the meanwhile, a comparison campaign could deploy the different emissions measurement devices concurrently at the same sample of stacks to provide a basis for homogenizing existing BC emissions datasets.
Unfortunately, [BC] can only be derived by indirect methods.Consequently, existing [BC] data are inconsistent due to 1) systematic differences in measurement errors among different types of devices, 2) correction procedures, 3) different measurement protocols, 4) inlet size, 5) different QA/QC methods, and 6) methods to improve the signal-to-noise ratio.We conclude that establishing generally accepted measurement protocols, and QA/QC standards for each of the indirect methods could improve consistency of measurements performed by the same equipment.If raw data from past measurements still exist, they should be stored together with the processed data in public databases including metadata on instrument settings, changes, QA/QC, corrections applied online, as well as information on co-located instruments (if available).If possible, raw data from the past have to be re-flagged according to an agreed-upon standard.Doing so would allow scientists to create a homogenized picture of past Arctic [BC] as a baseline for health and other impact studies, as well as determination of trends.
illustrates the consequences of such agency/law related discrepancies by two examples.Figure 2(a) reveals that in some Canadian provinces, agencies assign BC emissions from off-road transport to the entire province (see purple area), while others assign them to the actual mining area.In the North American Arctic, inland and coastal ship/boat traffic is accounted for as off-road transport.The Yukon, Mackenzie, and the shipping lanes along the coast and into Hudson Bay are clearly visible; some frozen rivers serve as ice roads in winter.On the contrary, Russian and European agencies report BC release from coastal shipping as ship emissions; the shipping lanes are clearly visible in the annual total emissions shown in Figure 1(a) (mostly green lines).Figure 2(b) illustrates consequences of lumping and spatial allocation.Including space heating into the energy sector makes the emissions of power generation facilities to areal sources.Reporting at county/borough/city levels leads to spatially heterogeneous distribution of emissions (e.g.Scandinavia); reporting at the province/state level leads to more homogeneous, and seemingly low emissions over a large area (e.g.southern part of the Northwest Territories).

Figure 2 .
Figure 2. Examples of artifacts in emission inventories due to regional differences in classification, spatial allocation, and lumping of emission sources.Shown are the total 2008 BC emissions associated with (a) Transportation, and (b) The energy sector in the EDGAR-HTAP inventory (Figure1(a)).Note that similar spatial discrepancies exist for other emissions sectors and emission inventories using reported emissions data (therefore not shown).
sample ambient air with typical flow rates of ~6.7 liter•min −1 .At high flow rates (>3.5 liter•min −1 ), diffusion losses remain below 4%.Non-size specific inlets permit sampling the total aerosol-size distribution during fog periods.In the Arctic, typically, a snow-hood and stainless steel-duct inlet are used.The sampled air should reach laboratory temperature (∼21˚C) prior to entering the actual measurement device to ensure relative humidity is less than 20%[57].Particles accumulate over a given time interval (e.g. 5 min) on a quartz fiber filter.A light source illuminates the sample and a blank filter on one side.On the other side, detectors measure real-time attenuation of light transmitted through the two filters.Typical operational effective wavelengths are 370, 470, 520, 590, 660, 880, and 950 nm.Multispectral (370 -950 nm) absorption coefficients can give information about the chemical composition.

Figure 3 .
Figure 3. Black carbon concentrations derived from aethalometer measurements at College, Alaska between September 11 and October 30, 2017.(a) Temporal evolution of 5-min mean [EBC].Data for 880 nm look similar (therefore not shown), but are higher than at 370 nm.(b) Scatter plot of 5-min mean [EBC] as obtained at 370 and 880 nm.
non are the corrected and non-corrected [ ] EBC , and k is the shadowing correction aka aethalometer-correction factor.Negative k-values mean overestimation of [ ] EBC by [ ]

Figure 4
Figure 4 exemplarily shows the mean diurnal course of 5-min [EBC] averaged over the aforementioned 50-d period.The means reveal peaks around noon and at night that correspond to lunch break and the nighttime inversion that frequently builds in the shoulder season.The differences indicate coated BC.For the 50-d period, hourly means of [EBC] 370nm and [EBC] 880nm correlated 90% (R = 0.949); daily means of [EBC] 370nm and [EBC] 880nm correlated 63% (R = 0.795).

Figure 4 .
Figure 4. Examples of offline post-processing of observations at College, Alaska.(a) Timeseries of hourly means of [EBC] for a 50-d period from September 11 to October 30, 2017 shown at 880 nm, (b) 50-d mean diurnal course of [EBC] at 370 and 880 nm, (c) Mean diurnal course of [EBC] averaged over all weekdays and weekends of the 50-d period both at 370 and 880 nm, (d) Period averaged 5-min means of [EBC] at 370 nm vs. 880 nm on weekdays and weekends.
Temporary co-location of filter-based absorption photometers (e.g.PSAP, MAAP, CLAP) can provide a site-dependent normalization factor permits comparison of BC measurements by identical aethalometer at different sites.This factor gives how much greater the attenuation coefficient is when compared to the light absorption coefficient of the co-located filter-based absorption photometer that had been corrected for loading and scattering artifacts (see Section 3.2).This normalization avoids statistics to be concentration-weighted.Recall filter-based absorption photometers are accurate within 20% -30% of the true σ ap .

Few
vertical profiles of [BC] exist.A 2-years (2011-2012) campaign with 200 vertical profiles up to 1200 m above ground level from tethered balloons at Ny-Ålesund showed the following major behaviors [82]: Background conditions showed homogeneous [BC] profiles, while transport caused strong vertical gradients.Surface inversions and/or local emissions led to high [BC] close to ground.Ground-based plumes of secondary aerosol formed from local emissions of pollutants led to decreasing [BC] with height dependent on size [82].During the Arctic Climate Change Economy and Society (ACCES) campaign, vertical profiles of BC mass mixing ratios (MMR) up to 11 km height were derived from SP2s flown onboard of two aircrafts over Iceland, West Greenland both in 2012 and 2013, as well as Scandinavia up north to Svalbard in 2012; BC MMR were ~10 ng/kg and nearly independent of height as expected given weak local sources/sinks and far distant source regions [83].Data collected in the North American Arctic during the April and June-July 2008 Arctic Research of the Composition of the Troposphere from Aircraft and Satellites (ARCTAS) campaign also showed that transport efficiency, source region and season determined [BC] profiles; on average over all Arctic flights, [BC] was 107 ± 55 ng•m −3 and 9.9 ± 5.7 ng•m −3 in spring and summer, respectively.

Figure 6
Figure6summarizes mean[BC]  in the sense of a climatology following[21] [85] with samples of [BC] taken for periods of different length, at different times and locations between 1969 and 2017 as found in the literature and calculated from publically available data.At Tiksi, 2-yr observations showed distinct seasonality in [BC] (8 -302 ng•m −3 ); biomass burning contributed 19% -73% to [BC][58].At Kevo, winter, spring, summer, and fall 1964-2010 mean [BC] were 339, 199, 127, and 213 ng•m −3 , respectively.Wood burning contributions to [BC] were marginal[61].At the five Finish sites, spring and winter concentrations were highest due to long-range-transport of BC, increased domestic wood burning and reduced atmospheric boundary-layer heights.Summer [BC] was lowest due to more effective vertical mixing.At all sites, highest [BC] occurred under southerly winds.Backward trajectories indicated Central and Eastern Europe as source regions.During these episodes, non-fire anthropogenic sources and open biomass

Figure 7 .
Figure 7. Examples of temporal variability as obtained from aetholometer measurements at selected sites.Alert 1989-2012 (a) Mean annual course of hourly means of [BC], (b) Mean diurnal course of hourly means of [BC], and (c) Mean September October (SO) 1989-2012 diurnal course of hourly means of [BC].Mean annual cycles at Barrow, and Ny-Ålesund show simular behavior as Alert.Mean diurnal course of 5-min means of [BC] on weekdays and weekends both at 370 and 880 nm at (d) Tiksi during SO 2010, and (e) Circle, Alaska 1-23-2018 to 3-13-2018 (LT).

Figure 7 (
Figure 7(d), Figure 7(e) and Figure 4(c) also reveal that Tiksi, Circle and College have similar mean diurnal courses with higher concentrations around noon and at night than in the morning and evening.Hourly means of [BC] also reflect the different typical work hours in the US (9 to 5 LT) vs. Russia (7-4 LT).Besides

Table 1 .
September to October statistics for selected [BC] mean 5-min values derived from aethalometer measurements in the Arctic.R 2 , R, Max, and StDev are the correlation between between 5-min [BC] 370nm and [BC] 880nm in percent, correlation coefficient, maximum concentration and the standard deviation for the respective sites.Regular and Italic values are at 370 nm and 880 nm, respectively.Data for Tiksi are 2010, College 2017.