In this Millennium, black carbon (BC) in the Arctic has gained huge public interest  . 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  .
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     . 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  . Consequently, cloud and precipitation formation as well as atmospheric removal processes might change with changing BC concentrations  . Changes in cloudiness, cloud-life time, coverage and precipitation again impact the energy cycle via changes in the radiation budget (e.g.   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 (PM2.5) is health adverse (e.g.     ). PM2.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 PM2.5 see   ; see  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 (NO2) and sulfur dioxide (SO2), as well as atmospheric water. These species may deposit on hygroscopic carbonaceous particles; also other aerosol components can join BC during particle growth  . 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  . Transport of BC emitted in mid latitudes contributes a substantial fraction to Arctic BC     .
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  typically have better air quality than cities in continental microthermal climates  .
Given the huge variability of BC-containing aerosol, and the fact that the chemical composition of PM may affect light-absorbing (soot/BC/EC) and/or may cause light-scattering (organics, inorganics) many different indirect optical measurement techniques and thermal optical analysis methods were developed. Some of them permit continuous, or semi-continuous (online) observations, while others work offline. Those instruments have been successfully applied in mid- and low-latitudes.
As of today, no overall Pan-Arctic comprehensive assessment of BC from emissions to removal exists. Emission inventories and measurements of concentration baselines, however, are pre-requisites to assess impacts from anticipated increases in Arctic shipping on human and ecosystem health as well as subsistence lifestyle. The goal of this paper is to 1) assess Arctic BC emissions; 2) review, evaluate current state-of-the-art BC measurement methods for Arctic application and consolidate knowledge on past; 3) Arctic BC concentrations at breathing level, 4) BC removal in the Arctic; and 5) identify gaps in knowledge and steps for further investigations. Here, we consider the region north of 59.9N as “Arctic”.
2. Black Carbon Emissions
Often countries report emission data at the state, province, borough or country level. Generally, estimates of BC emissions base on PM2.5 emissions from combustion processes multiplied by the fraction of BC in the PSM2.5 for the particular emission source  . All eight Arctic countries (Canada, Denmark, Finland, Iceland, Norway, Russia, Sweden, United States) have also territory outside the Arctic. Consequently, assessment of Arctic BC emissions is difficult. In the following, we discuss emission estimates for the Arctic countries at-large, and restrict us to the region north of 59.9 N when possible.
2.1. Arctic Emission Sources
Major Arctic BC-emission sources are shipping      , boreal wildfires and biomass burning  . In the new Millennium, Arctic ship traffic and hence BC emissions from shipping have seen notable increases due to the long record of, on average, decreasing sea-ice extend, increased last-chance tourism, and shipping of supply for offshore oil/gas extraction activities     . Further emission sources for BC are space heating  , solid fuel based power generation  , traffic  , on and off road diesel vehicles    , and flaring   . Small contributions stem from exploratory drilling, fire places, cook stoves, municipal waste-burning, as well as meat and fish smoking by indigenous people. In the European Arctic, sauna stoves are strong contributors to [BC]. In boreal regions, residential wood combustion is a major source of BC in winter  .
2.2. Emission Inventories
Generally, BC-emission measurements barely exist and lack a common standard. 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   for reviews on US and European BC source apportionment studies, and  for BC emissions in Russia.
Many emissions inventories only report PM10 and PM2.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    . For differences between bottom-up and top-down emissions inventories see    .
Figure 1 shows exemplarily gridded Arctic BC emission totals for 2008 from different inventories      . The Emissions Database for Global Atmospheric Research version 4.1 (EDGARv4.1) is based on worldwide consistent emission and activity factors, while EDGAR-HTAP (Hemispheric Transport of Air Pollution) inserted nationally reported emissions where available.
Discrepancies among emissions inventories occur from 1) different national classifications of BC sources, data collection methods (bottom-up vs. top-down); 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 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 1. Comparison of 2008 annual total anthropogenic BC emissions except from flares and fires as reported by (a) The EDGAR-HTAP  , and (b) Technology based  emission inventories; 2008 annual total BC emissions from (c) Flares  , and (d) Fires   . Annual emission totals of other years look similar (therefore not shown). One metric ton equals 0.001 Gg.
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
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 (Figure 1(a)). Note that similar spatial discrepancies exist for other emissions sectors and emission inventories using reported emissions data (therefore not shown).
example above, the total off-road emissions are equally distributed over the region while the actual real world emissions are released in an open-pit mine. Note that inversions in open-pit mines are a major air quality issue. In case of lumping areal sources like residential/commercial heating with point sources like generation of electrical power, local maxima are smoothed out. Lumping of emission sources also leads to incorrect emission factors and temporal allocations for n − 1 of the n lumped emission sources.
Besides fuel-type, ship emissions depend on speeds, engine load and size. These factors cause uncertainty in estimated BC emissions from Arctic shipping   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  .
Russia’s 2010 annual anthropogenic BC emissions were estimated as 223.7 Gg with 81.0, 56.0, 45.3, 29.3, and 12.1 Gg from gas ﬂaring, residential emissions, transportation, industries, and power plants, respectively  . The Evaluating the CLimate and Air Quality ImPacts of ShortlivEd Pollutants version 4.0 (ECLIPSE) reports the 2010 Russian BC emissions from residential heating, transportation, power, industry, and agriculture waste as 22.4, 52.0, 20.6, 5.64 and 24.6 Gg, respectively. These inventories differ mainly in their estimates of contributions from gas ﬂaring, residential, and industrial emissions. Russia’s 2014 anthropogenic BC emissions were estimated at 688 Gg with uncertainty between 401 and 1453 Gg; the mean annual total BC emitted by the 2002-2015 wildfires (forests, grasslands, other natural land) and 2003-2009 agricultural burning amounted 569 Gg (345 - 790 Gg) and 9 Gg (2 - 22 Gg)  . Here the values in brackets give the range of uncertainty. Russian space heating uses raw coal, light brown coal, briquettes, kerosene, residual fuel oil, liquefied petroleum gas, natural and refinery gas, coke-oven coke, firewood, peat, and charcoal  . The 2014 BC emissions from this sector and power generation were estimated between 11.9 and 26.7 Gg with a mean of 16.4 Gg  . BC emissions from the domestic sector were between 15 and 468 Gg centered at 26.7 Gg. The 2010 and 2014 BC emissions from other stationary sources were estimated at 238.1 Gg and 192.2 Gg, respectively. BC emissions from industry were 5.2 Gg with huge uncertainty (1.0 - 51.3 Gg)  . Estimates of total BC emissions in Russia due to transportation (diesel and gasoline vehicles, railroad, ships) range from 16.3 Gg to 41.2 Gg with a mean of 29.8 Gg  . For a comprehensive discussion of Russian BC emissions uncertainties see  .
According to Swedish government reports  , their total annual BC emissions were 4.81 Gg, 5.1 Gg, 4.2 Gg and 3.67 Gg in 2000, 2005, 2012, and 2013, respectively. Industrial processes, road traffic, off-road vehicles, energy sectors, and residential biomass burning contributed about 0.210 Gg, 0.714 Gg, 0.756 Gg, 0.798 Gg, and 1.134 Gg in 2012.
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  .
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  .
In Canada, BC emissions are calculated using PM2.5 emissions and BC/PM2.5 ratios. In 2015, they were assessed as 38 Gg  .
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   . Hydrocarbons can cause soot; impurities cause incomplete combustion. Flare gases with higher energy content than methane produce more BC. Russian flares are larger than Alaska flares. In recent years, flaring declined in Russia. In 2014, BC emissions from flaring in Russia were between 10 and 54 Gg centered at 32 Gg  .
Appreciable differences in total wildfire BC emissions exist among years depending on the severity of the wildfire season, type of land-cover and area burned  . Boreal fires, for instance, can emit BC from the belowground organic layer (e.g. peat, moss, lichen), litter layer, and the forest; black spruce BC emissions differ from white spruce, larch, birch or alder. The Global Fire Emissions Database version 3s (GFED3s), for instance, reported the 2000 biogenic BC emissions north of 60N in Alaska as 50.414 Gg; GFED4s reports BC emissions from North American, Asian and Finnish boreal forest fires as 0.5 Gg  .
3. Measuring Near-Surface Black Carbon Concentrations
As of today, no standardized method to measure [BC] exists    . [BC] are operationally defined as the light absorbing and thermal refractory fraction of carbonaceous aerosol, respectively  . 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  .
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  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  .
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)  . 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  . Unfortunately, dust concentrations are often high in the North American Arctic due to unpaved roads and non-managed rivers  .
Uncertainty in [EBC] is about 36% for aethalometer in the Arctic, 28% (Arctic haze) to 80% (summer) for thermal desorption techniques, 28% - 40% for SPSs, and 10% - 20% for PSAP  .
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 use optical absorption for in-situ (online) EBC detection. They 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%  .
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.
The sample attenuation is
where Is and Ir 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 Dt are the filter spot size area (m2), sample flow rate (m3・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.
Despite aethalometer measure , the output is [EBC]. Equivalent black carbon encompasses BC plus other light-absorbing material. To convert to [EBC] a wavelength-dependent mass-specific absorption coefficient or short specific absorption coefficient (m2・g−1) is used. It is defined as absorption per mass concentration. The aethalometer mass-specific absorption coefficient relates optical attenuation through the filter with equivalent black carbon mass concentration,
Typically, the specific absorption coefficient is based on calibrations and theoretical calculations performed during instrument development. The recommended conversion for the Magee aethalometer, for instance, is = 14625/λ. This value corresponds to 39.5, 31.1, 28.1, 24.8, 22.2, 16.6 and 15.4 m2・g−1 EBC at λ = 370, 470, 520, 590, 660, 880, and 950 nm, respectively.
Figure 3 shows data collected at College, Alaska (64.8594N, 147.8508W) between September 11 and October 30, 2017 at 370 and 880 nm. The aethalometer pulled air from a parking lot behind the International Arctic Research Center. [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 range from 5 to 20 m2・g−1, but are constant for a given atmospheric environment  . During the 1989 to 2007 [EBC] measurements at Alert, Canada (82.45083N, 62.34167W), for instance, the manufacturer-recommended conversion factor of = 19 m2・g−1 was applied  . Multiplication of aethalometer data with the applied conversion factor and dividing by a conversion factor X permits using X instead. See   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, accounts for this shadowing effect at low to moderate filter loads. At high filter loads, optical measurements of light absorption require additional corrections  .
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.
High values of the filter-loading parameter in the infrared indicate EBC of local sources, while low filter-loading parameters indicate aged BC coated by ammonium sulfate and/or secondary organics during long-range transport. This means the filter-loading effect can provide a proxy for particle coating, and differentiation between local/fresh and transported/aged particles  .
Various empirical correction methods exist to calculate the absorption coefficient, from the attenuation coefficient,  . A simple approach derived from observations at three Finnish sites is  .
where , and are the corrected and non-corrected , and k is the shadowing correction aka aethalometer-correction factor. Negative k-values mean overestimation of by . At the three Finnish sites, k varied between 0 and 0.015 with a mean of 0.0018  .
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  . In the European Arctic, annual shadowing correction is around 20%  .
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  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  . 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.
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). 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].
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.
3.2. 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  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 (Is) and clean (non-sampled) area (Ir) 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.
Table 1. September to October statistics for selected [BC] mean 5-min values derived from aethalometer measurements in the Arctic. R2, 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.
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  . Non-absorbing aerosols and suspended particles can yield overestimates of absorption at 550 nm of about 20% - 30% that can be corrected for  .
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  . 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 m2・g−1 and 6.5 m2・g−1, respectively  .
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  . The MAAP seems to be less prone to interpret light scattering ( ) as light absorption than a PSAP or aethalometer.
3.3. 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  . 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 (O2) mix to ≥ 550˚C. IR absorption or flame ionization serves to detect carbon dioxide or methane formed from the released gases  .
Overestimates of [EC] may occur when some OC undergoes pyrolysis or charring, which convert OC to EC at high temperatures  . 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 O2 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  .
Some thermal methods monitor the optical reflectance of the sample filter  to correct for charring and pyrolysis effects. Thermal optical reflectance (TOR) or thermal optical transmittance (TOT) methods provide apparent elemental carbon (ECa). 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  . 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  . Uncertainty results from the assumptions of no scattering losses and homogeneous optical filter properties.
3.4. 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 NO2 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%  . A recent study demonstrated that the sensitivity of PAS is too coarse to capture typical [BC] of the high Arctic.
3.5. 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  . 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; largest discrepancy occurred for large contributions of brown carbon to total carbon  . 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  . 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  , 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]  . This means at the typical Arctic [BC], distinguishing between signals, noise and inter-instrument uncertainty may become challenging.
Aethalometer, thermal desorption and SPS measurements at Alert between March 2011 and December 2013 showed that typically, [EBC] was 2.7 and [EC] was 3.1 times higher than [rBC]  . At Zeppelin station, Norway, annual mean [BC] from aethalometer observations in 1991, and 2001-2007  , and an optical method on W41-filters between 1991 and 2004  agreed well during overlapping times  .
[BC] determined by COSMOS, thermal-optical measurements, and SPS laser-induced incandescence agreed within 10%  . 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  .
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  . 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  .
Figure 5 compares measurements from two calibrated, co-located aethalometer performed at College, Alaska. For the time shown, at 370 nm, [BC] mean, standard deviation, skewness, and kurtosis were 8.4 ± 8.8 ng・m−3, 7.5, and 116.8 for aethalometer 545, and 7.7 ± 7.5 ng・m−3, 7.2, and 98.1 for aethalometer 361, 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.
Figure 5. Comparison of [EBC] as obtained by two co-located aethalometer in College. Data shown are for 1-8-2018 1230 LT to 1-16-2018. Temporal evolution as obtained by aethalometer 361 (purple) and 545 (green) at (a) 370 nm, (b) 880 nm, scatter plot of [EBC] obtained from aethalometer 545 vs. 361 for (c) 370 nm, and (d) 880 nm.
Observations at 370 and 880 nm correlated 83% (R = 0.923) for aehalometer 545, and 87% (R = 0.939) for aehalometer 361. During the time period, 5-min mean and maxima of [BC]370nm were 196.1 ng・m−3 and 148.5 ng・m−3 (24% difference) and those of [BC]880nm were 208.9 ng・m−3 and 220.7 ng・m−3 (5% difference) for instrument 545 and 361, respectively.
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  . 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).
3.6. Normalization of Data from Various Sites
When applied to different sites, aethalometer-correction algorithms often fail to remove the ATN-dependence from shadowing effects  . Values of 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  . 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  . 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  .
Temporary co-location of filter-based absorption photometers (e.g. PSAP, MAAP, CLAP) can provide a site-dependent normalization factor
That 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.
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 Cf 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 m2・g−1 and no loading corrections were used  . 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  . 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  . PSAP light absorption coefficients exceeded those of COSMOS by 22% and 43%, at Barrow (PM1) and Ny-Ålesund (PM10), respectively  . COSMOS σATN served to derive [BC] from the PSAP σATN measurements made since 1998  .
4. Arctic Black Carbon Observations
4.1. Black Carbon Vertical Profiles
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  : 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  .
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  . 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.
4.2. Black Carbon in the Near-Surface Atmosphere
Arctic BC measurements are scarce. At Kevo, atmospheric [BC] were derived from weekly filter samples using optical and thermal-optical methods from 1964 to 2010. This dataset is the longest Arctic [BC] record. At the research stations Alert, Barrow, Ny-Ålesund, Summit, Nord (Greenland), Pallas, and Tiksi, data exist from long-term monitoring of [BC] with aethalometer. At some places, data exist from individual research projects or field campaigns. At five background sites in Finland, [BC] observations from MAAPs and aethalometer exist; annual mean [BC] was highest at Virolahti (385 - 460 ng・m−3), followed by Hyytiälä (250 - 370 ng∙m−3), Utö (230 - 270 ng∙m−3), Puijo (225 - 230 ng∙m−3), and Pallastunturi (60 - 70 ng∙m−3)  . Typically, PM2.5 contained between 5 and 10% BC with highest percentage at Virolahti located close to the Russian border.
Figure 6 summarizes mean [BC] in the sense of a climatology following   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]  . 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  . 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 6. Composite of mean near-surface [BC] derived from monitoring, field campaigns or projects with various devices as reported in the literature        -  and EBAS. Data are means over periods of different lengths and/or years. Data of the RV Akademik Mstislav Keldysh cruises were digitized from Figure 2 in  .
burning contributed about 62% and 36% to total [BC], respectively; biomass burning occurred most often in spring  .
In the summers and winters from 1990 to 1992, [BC] were about 11 and 93 ng∙m−3 at Zeppelin (475 m ASL); in summer/fall and winter/spring 1979 to 1990, [BC] was 5 and 66 ng∙m−3 at Gruvebadet (0 m ASL), respectively. It is unclear whether the different elevation, time, changes in circulation, and/or emissions caused the quite different [BC]; generally, measurements at Zeppelin represent the free troposphere  . At Zeppelin, the 1998-2007 annual mean and median [BC] based on thermo-optical calibration were 39 and 27 ng∙m−3, respectively.
Calculation of mean annual courses of [BC] revealed stronger day-to-day variation and an order of magnitude or so higher concentrations in winter than summer (Figure 7(a)). [BC] data indicated strong relationships to upwind local emissions, meteorological conditions, and climatological regimes. In the temperate maritime climate along the Norwegian Coast, for instance, frontal passages occur frequently and remove particles by scavenging; in sub-arctic climates, frequent inversions yield accumulation of particles  . Multi-year annual mean diurnal courses differ strongly from monthly or seasonal mean diurnal courses (Figure 7(b), Figure 7(c)). Distinct differences exist in weekend and weekday hourly means and diurnal courses (Figure 7(d), Figure 7(e), Table 1).
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(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 seasonal, daily, and diurnal variability, tropical and mid-latitude long-term regimes govern Arctic [BC]. At Alert, for instance, observed [BC] was 40% higher during positive than negative phases of NAO  .
Various studies indicated changing [BC]. In Finland, for instance, annual mean [BC] decreased by 9.5 ng∙m−3 per decade between 2001 and 2007; this decrease is similar to the detection limit  . According to a Mann-Kendall trend test, Kevo annual mean [BC] decreased significantly by 78% with about 1.8% yr−1 from about 300 to 82 ng・m−3 between 1970 and 2010. This decrease showed peaks around 1976-1977, 1985-1987, and 1999 that well correlated with observed nickel in air concentrations. Obviously, emissions from extensive ore smelting on the Kola Peninsula contributed notably to PM during these years  . At Alert, [BC] decreased 61% with 3.4% yr−1 from 90 ng・m−3 to 35 ng・m−3 between 1989 and 2007; in the 1990s, net decreases of [BC] were about 54% and 27% at Alert and Barrow, respectively  . Various authors attribute the decreasing [BC] to the collapse of the economy in the former USSR   .
At Alert and Barrow, [BC] began to increase around 2000-2001  . At Alert, relative contributions from Eurasia to observed [BC] have declined from greater than 90% to about 75%; concurrently, those from North America have increased from less than 10% to nearly 25%  .
4.3. Black Carbon in Snow
Nucleation-scavenging is the major process for removal of BC-containing aerosol from the atmosphere  . Thus, cloud and precipitation formation and ultimately precipitation affect BC’s lifetime. Examining aqueous deposition for [BC] is uncommon.
Sampling snow to determine its BC mixing ratios started in the 1980s; it became more frequent since the 1990s. Figure 8 summarizes BC in snow mixing ratios found in the literature. In the European Arctic, 484 surface samples and 24 column samples (covering the accumulation season) were collected from 2007 to 2009. Mixing ratios were up to 88 ng of black carbon per gram of snow (ng・g−1) in Scandinavia; mixing ratios decreased with increasing latitude: 11 - 14 ng・g−-1 in Svalbard (74N - 81N, 10E - 34E), 7 - 42 ng・g−1 in Fram Strait, and 9 ng・g−1 in Barrow  . In Svalbard, the 2007 median of 81 samples was 4.1 ng・g−1 with values from 0 to 80.8 ng・g−1. The 2010 snow samples collected over the Canada Basin and Arctic Ocean north of 65N showed mean and median of 13.6 ng・g−1 and 10.7 ng・g−1 with mixing ratios ranging from 1.4 to 164.6 ng・g−1  . Snow samples on sea-ice collected in spring 2008-2013 between Greenland, Ellesmere Island, and the North Pole (82N - 89N, 0W - 100W) had medians of 4 ± 3 ng・g−1  . Observed [BC] in snow on sea-ice decreases from Arctic coastal regions to the center of the Arctic Ocean  (cf. also Figure 8).
According to 36 snow samples from sites across the Arctic (Alaska, Canada, Greenland, Russia, Arctic Ocean near the North Pole) from 2007 to 2009, often more than 75% of the BC in snow stemmed from biomass or biofuel combustion  .
Figure 8. Composite of BC in snow mixing ratios found in    .
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.
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.
5. 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.
We thank C.F. Cahill, G. Kramm, and the anonymous reviewers for fruitful discussion and helpful comments. The Alaska Department for Environmental Conservation supported this study with two of their aethalometer. Other data were downloaded from EBAS  , ARM data site, EDGAR site, and GFEDv4.1  ,     . Research computational support came from the University of Alaska Fairbanks’ Geophysical Institute. NEFCO and CATG provided financial support.
 Bond, T.C., Doherty, S.J., Fahey, D.W., Forster, P.M., Berntsen, T., DeAngelo, B.J., et al. (2013) Bounding the Role of Black Carbon in the Climate System: A Scientific Assessment. Journal of Geophysical Research: Atmospheres, 118, 5380-5552.
 Haywood, J. and Boucher, O. (2000) Estimates of the Direct and Indirect Radiative Forcing Due to Tropospheric Aerosols: A Review. Review of Geophysics, 38, 513-543.
 Dou, T.-F. and Xiao, C.-D. (2016) An Overview of Black Carbon Deposition and Its Radiative Forcing over the Arctic. Advances in Climate Change Research, 7, 115-122.
 Raju, M.P., Safai, P.D., Sonbawne, S.M. and Naidu, C.V. (2015) Black Carbon Radiative Forcing over the Indian Arctic Station, Himadri during the Arctic Summer of 2012. Atmospheric Research, 157, 29-36. https://doi.org/10.1016/j.atmosres.2015.01.013
 Moore, D., Copes, R., Fisk, R., Joy, R., Chan, K. and Brauer, M. (2006) Population Health Effects of Air Quality Changes Due to Forest Fires in British Columbia in 2003: Estimates from Physician-Visit Billing Data. Canadian Journal of Public Health, 97, 105-108.
 Liu, J.C., Pereira, G., Uhl, S.A., Bravo, M.A. and Bell, M.L. (2015) A Systematic Review of the Physical Health Impacts from Non-Occupational Exposure to Wildfire Smoke. Environmental Research, 136, 120-132.
 Asikainen, A., Pärjälä, E., Jantunen, M., Tuomisto, J.T. and Sabel, E.C. (2017) Effects of Local Greenhouse Gas Abatement Strategies on Air Pollutant Emissions and on Health in Kuopio, Finland. Climate, 5, 43. https://doi.org/10.3390/cli5020043
 Quarato, M., De Maria, L., Franca Gatti, M., Caputi, A., Mansi, F., Lorusso, P., et al. (2017) Air Pollution and Public Health: A PRISMA-Compliant Systematic Review. Atmosphere, 8, 183. https://doi.org/10.3390/atmos8100183
 Liu, H.-Y., Dunea, D., Iordache, S. and Pohoata, A. (2018) A Review of Airborne Particulate Matter Effects on Young Children’s Respiratory Symptoms and Diseases. Atmosphere, 9, 150. https://doi.org/10.3390/atmos9040150
 Janssen, N.A.H., Gerlofs-Nijland, M.E., Lanki, T., Salonen, R.O., Cassee, F., Hoek, G., et al. (2012) Health Effects of Black Carbon. In: Bohr, R., Ed., WHO.
 Sinkemani, R., Sinkemani, A., Li, X. and Chen, R. (2018) Risk of Cardiovascular Disease Associated with the Exposure of Particulate Matter (PM2.5): Review. Journal of Environmental Protection, 9, 607-618. https://doi.org/10.4236/jep.2018.96038
 Segersson, D., Eneroth, K., Gidhagen, L., Johansson, C., Omstedt, G., Nylén, A.E., et al. (2017) Health Impact of PM10, PM2.5 and Black Carbon Exposure Due to Different Source Sectors in Stockholm, Gothenburg and Umea, Sweden. International Journal of Environmental Research and Public Health, 14, 742. https://doi.org/10.3390/ijerph14070742
 Sharma, S., Andrews, E., Barrie, L.A., Ogren, J.A. and Lavoué, D. (2006) Variations and Sources of the Equivalent Black Carbon in the High Arctic Revealed by Long-Term Observations at Alert and Barrow: 1989-2003. Journal of Geophysical Research: Atmospheres, 111, D14208. https://doi.org/10.1029/2005JD006581
 Hirdman, D., Burkhart, J.F., Sodemann, H., Eckhardt, S., Jefferson, A., Quinn, P.K., et al. (2010) Long-Term Trends of Black Carbon and Sulphate Aerosol in the Arctic: Changes in Atmospheric Transport and Source Region Emissions. Atmospheric Chemistry and Physics, 10, 9351-9368.
 Matsui, H., Kondo, Y., Moteki, N., Takegawa, N., Sahu, L.K., Zhao, Y., et al. (2011) Seasonal Variation of the Transport of Black Carbon Aerosol from the Asian Continent to the Arctic during the Arctas Aircraft Campaign. Journal of Geophysical Research: Atmospheres, 116, D05202.
 Law, K.S., Stohl, A., Quinn, P.K., Brock, C.A., Burkhart, J.F., Paris, J.-D., et al. (2014) Arctic Air Pollution: New Insights from Polarcat-Ipy. Bulletin of the American Meteorological Society, 95, 1874-1895. https://doi.org/10.1175/BAMS-D-13-00017.1
 Kottek, M., Grieser, J., Beck, C., Rudolf, B. and Rubel, F. (2006) World Map of the Köppen-Geiger Climate Classification Updated. Meteorologische Zeitschrift, 15, 259-263.
 Mölders, N., Porter, S.E., Cahill, C.F. and Grell, G.A. (2010) Influence of Ship Emissions on Air Quality and Input of Contaminants in Southern Alaska National Parks and Wilderness Areas During the 2006 Tourist Season. Atmospheric Environment, 44, 1400-1413.
 Corbett, J.J., Lack, D.A., Winebrake, J.J., Harder, S., Silberman, A.J. and Gold, M. (2010) Arctic Shipping Emissions Inventories and Future Scenarios. Atmospheric Chemistry and Physics, 10, 10271-10311. https://www.atmos-chem-phys.net/10/9689/2010/
 Zhan, J., Gao, Y., Li, W., Chen, L., Lin, H. and Lin, Q. (2014) Effects of Ship Emissions on Summertime Aerosols at Ny-Ålesund in the Arctic. Atmospheric Pollution Research, 5, 500-510. https://doi.org/10.5094/APR.2014.059
 Winther, M., Christensen, J.H., Plejdrup, M.S., Ravn, E.S., Eriksson, ó.F. and Kristensen, H.O. (2014) Emission Inventories for Ships in the Arctic Based on Satellite Sampled AIS Data. Atmospheric Environment, 91, 1-14.
 McCarty, J.L., Krylov, A., Prishchepov, A.V., Banach, D.M., Tyukavina, A., Potapov, P., et al. (2016) Agricultural Fires in European Russia, Belarus, and Lithuania and Their Impact on Air Quality, 2002-2012. In: Garik, G. and Volker, R., Eds., Land-Cover and Land-Use Changes in Eastern Europe after the Collapse of the Soviet Union in 1991, Springer, Berlin, 193-221.
 Berkman, P.A., Vylegzhanin, A.N. and Young, O.R. (2016) Governing the Bering Strait Region: Current Status, Emerging Issues and Future Options. Ocean Development & International Law, 47, 186-217. https://doi.org/10.1080/00908320.2016.1159091
 Walsh, J.E., Fetterer, F., Scott Stewart, J. and Chapman, W.L. (2017) A Database for Depicting Arctic Sea Ice Variations Back to 1850. Geographical Review, 107, 89-107.
 Tran, H.N.Q. and Mölders, N. (2012) Numerical Investigations on the Contribution of Point Source Emissions to the PM2.5 Concentrations in Fairbanks, Alaska. Air Pollution Research, 3, 199-210. https://doi.org/10.5094/APR.2012.022
 Roiger, A., Thomas, J.L., Schlager, H., Law, K.S., Kim, J., Schäfler, A., et al. (2014) Quantifying Emerging Local Anthropogenic Emissions in the Arctic Region: The Access Aircraft Campaign Experiment. Bulletin of the American Meteorological Society, 96, 441-460.
 Pirjola, L., Paasonen, P., Pfeiffer, D., Hussein, T., Hämeri, K., Koskentalo, T., et al. (2006) Dispersion of Particles and Trace Gases nearby a City Highway: Mobile Laboratory Measurements in Finland. Atmospheric Environment, 40, 867-879.
 Robinson, M.A., Olson, M.R., Liu, Z.G. and Schauer, J.J. (2015) The Effects of Emission Control Strategies on Light Absorbing Carbon Emissions from a Modern Heavy Duty Diesel Engine. Journal of the Air & Waste Management Association, 65, 759-766.
 Kholod, N., Evans, M., Gusev, E., Yu, S., Malyshev, V., Tretyakova, S., et al. (2016) A Methodology for Calculating Transport Emissions in Cities with Limited Traffic Data: Case Study of Diesel Particulates and Black Carbon Emissions in Murmansk. Science of the Total Environment, 547, 305-313.
 Evans, M., Kholod, N., Kuklinski, T., Denysenko, A., Smith, S.J., Staniszewski, A., et al. (2017) Black Carbon Emissions in Russia: A Critical Review. Atmospheric Environment, 163, 9-21. https://doi.org/10.1016/j.atmosenv.2017.05.026
 Calvello, M., Esposito, F., Lorusso, M. and Pavese, G. (2017) A Two-Year Database of BC Measurements at the Biggest European Crude Oil Pre-Treatment Plant: A Comparison with Organic Gaseous Compounds and PM10 Loading. Atmospheric Research, 164-165, 156-166.
 Quincey, P., Butterfield, D., Green, D., Coyle, M. and Cape, J.N. (2009) An Evaluation of Measurement Methods for Organic, Elemental and Black Carbon in Ambient Air Monitoring Sites. Atmospheric Environment, 43, 5085-5091.
 Briggs, N.L. and Long, C.M. (2016) Critical Review of Black Carbon and Elemental Carbon Source Apportionment in Europe and the United States. Atmospheric Environment, 144, 409-427. https://doi.org/10.1016/j.atmosenv.2016.09.002
 Wang, R., Tao, S., Shen, H., Huang, Y., Chen, H., Balkanski, Y., et al. (2014) Trend in Global Black Carbon Emissions from 1960 to 2007. Environmental Science & Technology, 48, 6780-6787. https://doi.org/10.1021/es5021422
 Kurokawa, J., Ohara, T., Morikawa, T., Hanayama, S., Janssens-Maenhout, G., Fukui, T., et al. (2013) Emissions of Air Pollutants and Greenhouse Gases over Asian Regions during 2000-2008: Regional Emission Inventory in Asia (REAS) Version 2. Atmospheric Chemistry and Physics, 13, 11019-11058. https://doi.org/10.5194/acp-13-11019-2013
 Van der Werf, G.R., Randerson, J.T., Giglio, L., van Leeuwen, T.T., Chen, Y., Rogers, B.M., et al. (2017) Global Fire Emissions Estimates during 1997-2016. Earth System Science Data, 9, 697-720. https://doi.org/10.5194/essd-9-697-2017
 Wiedinmyer, C., Akagi, S.K., Yokelson, R.J., Emmons, L.K., Al-Saadi, J.A., Orlando, J.J., et al. (2011) The Fire Inventory from NCAR (FINN): A High Resolution Global Model to Estimate the Emissions from Open Burning. Geoscientific Model Development, 4, 625-641.
 Griffis, T.J., Lee, X., Baker, J.M., Russelle, M.P., Zhang, X., Venterea, R., et al. (2013) Reconciling the Differences between Top-Down and Bottom-Up Estimates of Nitrous Oxide Emissions for the U.S. Corn Belt. Global Biogeochemical Cycles, 27, 746-754.
 Amnuaylojaroen, T., Barth, M.C., Emmons, L.K., Carmichael, G.R., Kreasuwun, J., Prasitwattanaseree, S., et al. (2014) Effect of Different Emission Inventories on Modeled Ozone and Carbon Monoxide in Southeast Asia. Atmospheric Chemistry and Physics, 14, 12983-13012. https://doi.org/10.5194/acp-14-12983-2014
 Van der Gon, H.D., Beevers, S., D’Allura, A., Finardi, S., Honoré, C., Kuenen, J., et al. (2012) Discrepancies between Top-Down and Bottom-Up Emission Inventories of Megacities: The Causes and Relevance for Modeling Concentrations and Exposure. Springer, Dordrecht, 199-204.
 Timmermans, R.M.A., Denier van der Gon, H.A.C., Kuenen, J.J.P., Segers, A.J., Honoré, C., Perrussel, O., et al. (2013) Quantification of the Urban Air Pollution Increment and Its Dependency on the Use of Down-Scaled and Bottom-Up City Emission Inventories. Urban Climate, 6, 44-62. https://doi.org/10.1016/j.uclim.2013.10.004
 Lack, D.A. and Corbett, J.J. (2012) Black Carbon from Ships: A Review of the Effects of Ship Speed, Fuel Quality and Exhaust Gas Scrubbing. Atmospheric Chemistry and Physics, 12, 3985-4000. https://doi.org/10.5194/acp-12-3985-2012
 Mölders, N., Gende, S. and Pirhalla, M.A. (2013) Assessment of Cruise-Ship Activity Influences on Emissions, Air Quality, and Visibility in Glacier Bay National Park. Atmospheric Pollution Research, 4, 435-445. https://doi.org/10.5094/APR.2013.050
 Hitzenberger, R., Petzold, A., Bauer, H., Ctyroky, P., Pouresmaeil, P., Laskus, L., et al. (2006) Intercomparison of Thermal and Optical Measurement Methods for Elemental Carbon and Black Carbon at an Urban Location. Environmental Science & Technology, 40, 6377-6383.
 Petzold, A., Ogren, J.A., Fiebig, M., Laj, P., Li, S.M., Baltensperger, U., et al. (2013) Recommendations for Reporting “Black Carbon” Measurements. Atmospheric Chemistry and Physics, 13, 8365-8379.
 Backman, J., Schmeisser, L., Virkkula, A., Ogren, J.A., Asmi, E., Starkweather, S., et al. (2017) On Aethalometer Measurement Uncertainties and an Instrument Correction Factor for the Arctic. Atmospheric Measurements and Technology, 10, 5039-5062. https://doi.org/10.5194/amt-10-5039-2017
 Winiger, P., Andersson, A., Eckhardt, S., Stohl, A. and Gustafsson, &OUML;. (2016) The Sources of Atmospheric Black Carbon at a European Gateway to the Arctic. Nature Communications, 7, Article No. 12776. https://doi.org/10.1038/ncomms12776
 Petzold, A. and Schönlinner, M. (2004) Multi-Angle Absorption Photometry—A New Method for the Measurement of Aerosol Light Absorption and Atmospheric Black Carbon. Journal of Aerosol Science, 35, 421-441.
 Liousse, C., Cachier, H. and Jennings, S.G. (1993) Optical and Thermal Measurements of Black Carbon Aerosol Content in Different Environments: Variation of the Specific Attenuation Cross-Section, Sigma (Σ). Atmospheric Environment. Part A. General Topics, 27, 1203-1211.
 Dutkiewicz, V.A., de Julio, A.M., Ahmed, T., Liang, J., Hopke, P.K, Skeire, et al. (2014) Forty-Seven Years of Weekly Atmospheric Black Carbon Measurements in the Finnish Arctic: Decrease in Black Carbon with Declining Emissions. Journal of Geophysical Research: Atmospheres, 119, 7667-7683. https://doi.org/10.1002/2014JD021790
 Moosmüller, H., Chakrabarty, R.K. and Arnott, W.P. (2009) Aerosol Light Absorption and Its Measurement: A Review. Journal of Quantitative Spectroscopy and Radiative Transfer, 110, 844-878. https://doi.org/10.1016/j.jqsrt.2009.02.035
 Drinovec, L., Gregorič, A., Zotter, P., Wolf, R., Bruns, E.A., Prévôt, A.S.H., et al. (2017) The Filter-Loading Effect by Ambient Aerosols in Filter Absorption Photometers Depends on the Coating of the Sampled Particles. Atmospheric Measurement Techniques, 10, 1043-1059.
 Virkkula, A., Mäkelä, T., Hillamo, R., Yli-Tuomi, T., Hirsikko, A., Hämeri, K., et al. (2007) A Simple Procedure for Correcting Loading Effects of Aethalometer Data. Journal of the Air & Waste Management Association, 57, 1214-1222.
 Hagler, G.S.W., Yelverton, T.L.B., Vedantham, R., Hansen, A.D.A. and Turner, J.R. (2011) Post-Processing Method to Reduce Noise while Preserving High Time Resolution in Aethalometer Real-Time Black Carbon Data. Aerosol and Air Quality Research, 11, 539-546.
 Ogren, J.A., Wendell, J., Andrews, E. and Sheridan, P.J. (2017) Continuous Light Absorption Photometer for Long-Term Studies. Atmospheric Measurement Techniques, 10, 4805-4818. https://doi.org/10.5194/amt-10-4805-2017
 Bond, T.C., Anderson, T.L. and Campbell, D. (1999) Calibration and Intercomparison of Filter-Based Measurements of Visible Light Absorption by Aerosols. Aerosol Science and Technology, 30, 582-600. https://doi.org/10.1080/027868299304435
 Kondo, Y., Sahu, L., Moteki, N., Khan, F., Takegawa, N., Liu, X., et al. (2011) Consistency and Traceability of Black Carbon Measurements Made by Laser-Induced Incandescence, Thermal-Optical Transmittance, and Filter-Based Photo-Absorption Techniques. Aerosol Science and Technology, 45, 295-312.
 Müller, T., Henzing, J.S., de Leeuw, G., Wiedensohler, A., Alastuey, A., Angelov, H., et al. (2011) Characterization and Intercomparison of Aerosol Absorption Photometers: Result of Two Intercomparison Workshops. Atmospheric Measurement Techniques, 4, 245-268.
 Cavalli, F., Viana, M., Yttri, K.E., Genberg, J. and Putaud, J.P. (2010) Toward a Standardised Thermal-Optical Protocol for Measuring Atmospheric Organic and Elemental Carbon: The EuSAAR Protocol. Atmospheric Measurement Techniques, 3, 79-89. https://doi.org/10.5194/amt-3-79-2010
 Chow, J.C., Watson, J.G., Chen, L.W.A., Arnott, W.P., Moosmüller, H. and Fung, K. (2004) Equivalence of Elemental Carbon by Thermal/Optical Reflectance and Transmittance with Different Temperature Protocols. Environmental Science & Technology, 38, 4414-4422. https://doi.org/10.1021/es034936u
 Hitzenberger, R., Dusek, U. and Berner, A. (1996) Black Carbon Measurements Using an Integrating Sphere. Journal of Geophysical Research: Atmospheres, 101, 19601-19606. https://doi.org/10.1029/95JD02412
 Lack, D.A., Lovejoy, E.R., Baynard, T., Pettersson, A. and Ravishankara, A.R. (2006) Aerosol Absorption Measurement Using Photoacoustic Spectroscopy: Sensitivity, Calibration, and Uncertainty Developments. Aerosol Science and Technology, 40, 697-708.
 Reisinger, P., Wonaschütz, A., Hitzenberger, R., Petzold, A., Bauer, H., Jankowski, N., et al. (2008) Intercomparison of Measurement Techniques for Black or Elemental Carbon under Urban Background Conditions in Wintertime: Influence of Biomass Combustion. Environmental Science & Technology, 42, 884-889.
 Mölders, N., Tran, H.N.Q., Quinn, P., Sassen, K., Shaw, G.E. and Kramm, G. (2011) Assessment of WRF/Chem to Capture Sub-Arctic Boundary Layer Characteristics During Low Solar Irradiation Using Radiosonde, Sodar, and Station Data. Atmospheric Pollution Research, 2, 283-299.
 Husain, L., Dutkiewicz, V.A. and Maenhaut, W. (2011) Variation in Aerosol Black Carbon in Ny-Ålesund, Spitsbergen, Norway, from 1991 to 2004. The Arctic as a Messenger for Global Processes—Climate Change and Pollution, Copenhagen.
 Kramm, G. and Meixner, F.X. (2000) On the Dispersion of Trace Species in the Atmospheric Boundary Layer: A Re-Formulation of the Governing Equations for the Turbulent Flow of the Compressible Atmosphere. Tellus, 52A, 500-522.
 Collaud Coen, M., Weingartner, E., Apituley, A., Ceburnis, D., Fierz-Schmidhauser, R., Flentje, H., et al. (2010) Minimizing Light Absorption Measurement Artifacts of the Aethalometer: Evaluation of Five Correction Algorithms. Atmospheric Measurement Techniques, 3, 457-474. https://doi.org/10.5194/amt-3-457-2010
 Sinha, P.R., Kondo, Y., Koike, M., Ogren, J.A., Jefferson, A., Barrett, T.E., et al. (2017) Evaluation of Ground-Based Black Carbon Measurements by Filter-Based Photometers at Two Arctic Sites. Journal of Geophysical Research: Atmospheres, 122, 3544-3572.
 Ferrero, L., Cappelletti, D., Busetto, M., Mazzola, M., Lupi, A., Lanconelli, C., et al. (2016) Vertical Profiles of Aerosol and Black Carbon in the Arctic: A Seasonal Phenomenology along 2 Years (2011-2012) of Field Campaigns. Atmospheric Chemistry and Physics, 16, 12601-12629.
 Schwarz, J.P., Weinzierl, B., Samset, B.H., Dollner, M., Heimerl, K., Markovic, M.Z., et al. (2017) Aircraft Measurements of Black Carbon Vertical Profiles Show Upper Tropospheric Variability and Stability. Geophysical Research Letters, 44, 1132-1140.
 Hyvärinen, A.P., Kolmonen, P., Kerminen, V.M., Virkkula, A., Leskinen, A., Komppula, M., et al. (2011) Aerosol Black Carbon at Five Background Measurement Sites over Finland, a Gateway to the Arctic. Atmospheric Environment, 45, 4042-4050.
 Viidanoja, J., Sillanpää, M., Laakia, J., Kerminen, V.-M., Hillamo, R., Aarnio, P., et al. (2002) Organic and Black Carbon in PM2.5 and PM10: 1 Year of Data from an Urban Site in Helsinki, Finland. Atmospheric Environment, 36, 3183-3193.
 Pakkanen, T.A., Mäkelä, T., Hillamo, R.E., Virtanen, A., Rönkkö, T. and Hämerti, K. (2006) Monitoring of Black Carbon and Size-Seggregated Particle Number Concentrations at 9-m and 65-m Distances from a Major Road in Helsinki. Boreal Environment Research, 11, 295-309.
 Enroth, J., Saarikoski, S., Niemi, J., Kousa, A., Ježek, I., Močnik, G., et al. (2016) Chemical and Physical Characterization of Traffic Particles in Four Different Highway Environments in the Helsinki Metropolitan Area. Atmospheric Chemistry Physics, 16, 5497-5512.
 Shevchenko, V., Starodymova, D., Vinogradova, A., Lisitzin, A., Makarov, V., Popova, S., et al. (2015) Elemental and Organic Carbon in Atmospheric Aerosols over the Northwestern Coast of Kandalaksha Bay of the White Sea. Doklady Earth Sciences, 461, 242-246.
 Quinn, P.K., Shaw, G., Andrews, E., Dutton, E.G., Ruoho-Airola, T. and Gong, S.L. (2007) Arctic Haze: Current Trends and Knowledge Gaps. Tellus, 59B, 99-114.
 Gong, S.L., Zhao, T.L., Sharma, S., Toom-Sauntry, D., Lavoué, D., Zhang, X.B., et al. (2010) Identification of Trends and Interannual Variability of Sulfate and Black Carbon in the Canadian High Arctic: 1981-2007. Journal of Geophysical Research: Atmospheres, 115, D07305. https://doi.org/10.1029/2009JD012943
 Ching, J., West, M. and Riemer, N. (2018) Quantifying Impacts of Aerosol Mixing State on Nucleation-Scavenging of Black Carbon Aerosol Particles. Atmosphere, 9, 17. https://doi.org/10.3390/atmos9010017
 Doherty, S.J., Dang, C., Hegg, D.A., Zhang, R. and Warren, S.G. (2014) Black Carbon and Other Light-Absorbing Particles in Snow of Central North America. Journal of Geophysical Research: Atmospheres, 119, 12,807-12,831.
 Dou, T., Xiao, C., Shindell, D.T., Liu, J., Eleftheriadis, K., Ming, J., et al. (2012) The Distribution of Snow Black Carbon Observed in the Arctic and Compared to the GISS-PUCCINI Model. Atmospheric Chemistry and Physics, 12, 7995-8007.
 Forsström, S., Isaksson, E., Skeie, R.B., Ström, J., Pedersen, C.A., Hudson, S.R., et al. (2013) Elemental Carbon Measurements in European Arctic Snow Packs. Journal of Geophysical Research: Atmospheres, 118, 13,614-13,627.
 Hegg, D.A., Clarke, A.D., Doherty, S.J. and Ström, J. (2011) Measurements of Black Carbon Aerosol Washout Ratio on Svalbard. Tellus B, 63, 891-900.
 Tørseth, K., Aas, W., Breivik, K., Fjæraa, A.M., Fiebig, M., Hjellbrekke, A.G., et al. (2012) Introduction to the European Monitoring and Evaluation Programme (EMEP) and Observed Atmospheric Composition Change during 1972-2009. Atmospheric Chemistry and Physics, 12, 5447-5481. https://doi.org/10.5194/acp-12-5447-2012
 Randerson, J.T., van der Werf, G.R., Giglio, L., Collatz, G.J. and Kasibhatla, P.S. (2018) Global Fire Emissions Database, Version 4.1 (GFEDv4). ORNL DAAC, Oak Ridge.
 Doherty, S.J., Warren, S.G., Grenfell, T.C., Clarke, A.D. and Brandt, R.E. (2010) Light-Absorbing Impurities in Arctic Snow. Atmospheric Chemistry and Physics, 10, 11647-11680.