Enhanced light absorption by mixed source black and brown carbon particles in UK winter.

Black carbon (BC) and light-absorbing organic carbon (brown carbon, BrC) play key roles in warming the atmosphere, but the magnitude of their effects remains highly uncertain. Theoretical modelling and laboratory experiments demonstrate that coatings on BC can enhance BC's light absorption, therefore many climate models simply assume enhanced BC absorption by a factor of ∼1.5. However, recent field observations show negligible absorption enhancement, implying models may overestimate BC's warming. Here we report direct evidence of substantial field-measured BC absorption enhancement, with the magnitude strongly depending on BC coating amount. Increases in BC coating result from a combination of changing sources and photochemical aging processes. When the influence of BrC is accounted for, observationally constrained model calculations of the BC absorption enhancement can be reconciled with the observations. We conclude that the influence of coatings on BC absorption should be treated as a source and regionally specific parameter in climate models.

Quantitative prediction of the direct warming effect of absorbing carbonaceous particles is challenging but crucial to assess their climate forcing1234. The major atmospheric absorbing carbonaceous components are black carbon (BC) and brown carbon (BrC). BC strongly absorbs across the solar spectrum, whereas BrC preferentially absorbs at short wavelengths25. Uncertainties in the absorptive properties of BC and BrC contribute substantially to the total uncertainty estimates of their direct warming effects, in addition to other uncertainties associated with emission estimates, atmospheric transport simulations, and removal rates. Bond et al.6 suggest that the mass absorption cross-section (MAC) of uncoated BC has an uncertainty of 16% at 550 nm. Additional uncertainties in absorption can be introduced as the optical features of BC are continuously modified during atmospheric aging processes. These aging processes mix BC with non-BC species, including BrC. Such mixing, and more specifically formation of coatings on BC particles, can occur over a few hours after BC emission and continue for days7. The coatings can enhance the absorption above that of uncoated BC particles through the so-called ‘lensing effect.' The lensing effect has been confirmed by theoretical calculations891011 and laboratory measurements1213141516 but is not in line with field observations, discussed further below. The inconsistency suggests that using a lensing effect calculated from Mie theory may introduce ∼50% uncertainty in BC's absorption, which can translate to ∼50% uncertainty of modelled BC's direct warming effect. An additional uncertainty of ∼20% may be introduced when including the absorption of BrC4, which is ignored by most current climate models. Accordingly, it is important to constrain the dynamic absorptive properties of BC and BrC to improve model predictions.

The influence of BC coatings on BC's absorption has been studied theoretically and experimentally in the laboratory and in the field. Mie theory calculations that assume spherical BC-containing particles with a core–shell configuration, such as those employed by many current climate models, indicate that enhancement factors of up to 3 are plausible891011. Laboratory studies demonstrate that such absorption enhancement can occur1213141516. Therefore, current climate models typically either calculate an enhancement factor based on simplified mixing-state assumptions1718 or assume a constant enhancement value of ∼1.5 for calculating BC absorption1920. In contrast, field observations of BC absorption enhancement in urban and coastal California (USA) and Japan142122 yield small values (on average no more than 1.1 at 405, 532, and 781 nm), even when substantial coatings are present1421. If true globally, this result questions BC's role as the second most important anthropogenic climate warmer after carbon dioxide13. The difference may be due to a preponderance of ambient BC being only partially encapsulated or attached to the host particle surface23. However, larger absorption enhancements of 1.4 (at 532 nm) have been observed in air masses strongly impacted by biomass burning24. Because of the paucity of comprehensive ambient measurements, the mechanisms responsible for the model-observation discrepancies remain elusive, underscoring the need for field studies in other representative environments, in particular those with mixed BC and BrC sources.

We perform an in-depth field study of light absorption enhancement factors (Eabs) for BC particles emitted from fossil fuel and residential burning sources in wintertime UK. The measurements are part of the 2012 Clean Air for London (ClearfLo) project25 at Detling, a rural site 45 km southeast of central London. We demonstrate substantial light absorption enhancement of BC from our measurements. Through detailed chemical and microphysical analyses, we find that the magnitude of the observed Eabs strongly depends on the amount of coatings on BC, the BC sources, and the extent of particle aging. We also show that the Eabs at 405 nm is affected by low-volatility BrC. In addition to the observations, we systematically examine the theoretically calculated Eabs and determine the refractive indices for BC and BrC that result in agreement between calculated and observed Eabs. Furthermore, single-particle morphological analysis provides mechanistic insight of Eabs in comparison to previous studies. In light of these analyses, we conclude that the Eabs is source and regionally dependent. Therefore, parameters used for Mie modelling should be calibrated with observations to account for this effect.

Results

Direct observation of E abs

Observed Eabs values were calculated using two independent and complementary methods. In one, Eabs was taken as the ratio between the MAC observed for ambient particles and for particles heated at 250 °C. BC is typically considered refractory at these temperatures and the heating is therefore used to remove the part of the coating material that is volatile. Particle heating was achieved by passing particles through a heated thermodenuder (TD). MAC is defined as the ratio of the absorption coefficient (babs) to refractory BC (rBC) concentration ([rBC]): MAC=babs/[rBC]. babs at 405 and 781 nm were quantified using photoacoustic spectroscopy26 and [rBC] was measured using a single-particle soot photometer (SP2)2728. Normalizing babs by [rBC] accounted for particle losses in the TD and temporal variations of babs during each ambient and TD measurement cycle. The observed Eabs calculated using this method is termed (λ denotes wavelength). The campaign average and values were 1.3 and 1.4, respectively, with their time series and frequency distributions shown in Supplementary Fig. 1. The second method was based on the absolute MAC values for ambient particles and the estimated enhancement is termed . was calculated as , where MACobs is the absolute MAC for ambient particles (shown in Supplementary Fig. 2.), and MACref is the literature reference value for uncoated BC (7.5±1.2 m2 g−1 at 550 nm)6. Using the inverse wavelength dependence relationship for uncoated BC6, we derived MACref values of 10.2±1.6 and 5.3±0.8 m2 g−1 at 405 and 781 nm, respectively.

Since BrC absorbs at 405 nm but not at 781 nm, the observed Eabs at 405 nm (using both methods) includes the lensing-driven enhancement and the enhancement induced by semi-volatile BrC, whereas the observed Eabs at 781 nm represents the lensing-driven enhancement only. A glossary of the abbreviations used in this study is listed in Supplementary Table 1.

Effects of internal mixing on E abs

We measured the composition of total non-refractory organic mass (OM) in submicron particles using a high-resolution time-of-flight aerosol mass spectrometer (HR-ToF-AMS)2930 and the non-refractory OM internally mixed with BC with a soot particle aerosol mass spectrometer (SP-AMS)31. The SP-AMS was operated in laser-only mode, which allows for detection of only rBC-containing particles (Methods). The ratio of non-rBC mass to rBC mass in rBC-containing particles is termed RBC. Assuming a core–shell particle configuration, RBC is a measure of relative coating thickness.

The observed Eabs at Detling increased substantially with increasing RBC (Fig. 1a,b) unlike the previous studies in California1421, even though the maximum RBC values observed here were three times smaller. In our study, there was good agreement between and , where BrC is thought to absorb negligibly. This result indicates that a lensing-driven absorption enhancement was, at times, substantial in Detling. It also suggests that residual non-rBC materials did not have a notable impact on absorption by the thermodenuded particles. In contrast, was notably higher than . Furthermore, was larger than at a given RBC. These features resulted from the direct contribution of BrC on the absorption at 405 nm.

Figure 1

Variation of Eabs and particle composition as a function of RBC.

(a,b) show measured and calculated absorption enhancement values versus RBC at 405 and 781 nm, respectively. kBrC of 0.004 and 0.032 were used for modelling absorption of ambient and thermodenuded organic mass at 405 nm, respectively. (c) Particle average oxidation state versus RBC. The points are coloured by the organic mass fraction of the oxygenated organic aerosol factor, which was derived from factor analysis of the BC-associated OM. (d) Mass fraction of the non-refractory components internally mixed with BC (shaded areas) and BC core median volume-weighted diameter (open circles) as a function of RBC. The colours represent nitrate (blue), ammonium (orange), sulfate (red), chloride (purple), oxygenated organic aerosol factor (pink), solid fuel organic aerosol factor (brown) and hydrocarbon-like organic aerosol factor (grey). The error bars in (a–d) represent s.d. of the values for each RBC interval.

As BrC contributes to the absorption at 405 nm, it is necessary to separate the absorption by BrC and evaluate the effect of BrC on Eabs at 405 nm. babs at 405 nm was apportioned into absorption by residual BC (BC and low-volatility coatings on BC that remain on the particles after thermodenuding at 250 °C), low-volatility BrC that did not fully evaporate in the TD (BrCLV), semi-volatile BrC that evaporated in the TD (BrCSV), and the increment of BC absorption due to the lensing effect, with their respective fractions of 53, 24, 7 and 16% (Methods). The apportioned values were similar to the apportionment results of Lack et al.24 at 404 nm inside a wildfire plume (Supplementary Table 2), even though we sampled a mixture of regional fossil fuel and residential burning sources with much lower absorption and OM concentrations. Given that the absorption of BrCLV was 3.4 times the absorption of BrCSV and that the mass of BrCLV was only 9.3% of the mass of BrCSV (measured by the HR-ToF-AMS assuming all OM is BrC), we estimated that the MAC of the BrCLV was ∼36 times the MAC of BrCSV. The much higher MAC of BrCLV was likely because BrCLV contained larger molecules with more conjugated carbon bonds (compared with BrCSV) that have lower volatility. This result is consistent with a recent study that shows BrC absorption from biomass burning is mostly associated with compounds of extremely low volatility32. Since the total non-refractory OM was substantially greater than the non-refractory OM that was associated with rBC, it is likely that much of the BrC was externally mixed with rBC.

The observed Eabs at 405 nm using the TD method was affected by BrCLV. The increase of BrCLV with RBC (Supplementary Fig. 3a) suggests a larger contribution of BrCLV to the absorption of thermodenuded particles as the rBC coating thickness increased, resulting in lower at high RBC. As a result, leveled off for RBC values greater than ∼3 (Fig. 1a). In contrast, continuously increased with RBC (Fig. 1b) as the absorption of BrC was negligible at 781 nm.

Dependence of E abs on sources and aging

To elucidate the mechanisms causing the enhanced BC absorption and formation of BrCLV, the BC-associated OM (measured by the SP-AMS) and total OM (measured by the HR-ToF-AMS) were apportioned into linearly independent components using positive matrix factorization3334. The analysis resulted in three factors for both BC-associated OM and total OM, namely oxygenated organic aerosol, solid fuel organic aerosol, and hydrocarbon-like organic aerosol. Based on the mass spectra, correlations with source tracers, and correlations with air mass sources from back trajectories, the oxygenated organic aerosol components were associated with long-range transport and aging of particles, the solid fuel organic aerosol components were associated with solid fuel burning for residential heating including biomass and coal, and the hydrocarbon-like organic aerosol components were associated with traffic emissions. The organic carbon in the total oxygenated organic aerosol component (oxygenated organic carbon) was further separated into fossil and non-fossil fractions using a combination of radiocarbon (14C), factor analysis, and BC source analyses (Methods). The apportionment showed that total oxygenated organic carbon (and likely BC-associated oxygenated organic carbon) was largely (73–90%) derived from the non-fossil sources. In addition, we find that the absorption by BrCLV correlated with the total oxygenated organic aerosol component (Supplementary Fig. 3b), indicating that BrCLV may be secondary, as previous laboratory studies proposed3536. In contrast, such correlation was not found for BrCSV.

The variation of Eabs with RBC was associated with changing sources as well as photochemical processing (Fig. 1c,d). The overall composition of rBC-containing particles varied with RBC (Fig. 1d). The relative abundance of secondary components (inorganics and oxygenated organic aerosol component) increased with RBC, while the relative abundance of primary components (hydrocarbon-like organic aerosol and solid fuel organic aerosol components) decreased. At lower RBC, the BC-associated OM was dominated by traffic sources and the Eabs value (∼1.1) was similar to the value observed at much higher RBC by Cappa et al.14, in which fossil fuel combustion is the dominant source for BC. Since the total oxygenated organic carbon was largely from non-fossil sources, the increase in RBC was correlated with changing sources that could produce BC particles in distinct size regimes. Given that BC particles from biomass burning are typically larger than those from fossil fuel combustion37, the correlation of BC core median volume-weighted diameter (measured by the SP2) with RBC (Fig. 1d) was consistent with the oxygenated organic aerosol components being increasingly produced by biomass combustion related sources as the oxygenated organic aerosol concentration increased. The increase of BC diameter with RBC may also indicate that coagulation during transport increased the BC size. However, since coagulation rates scale with the square of particle number (which decreases away from sources), and only coagulation between two BC-containing particles can lead to an increase in the BC core size, this hypothesis seems unlikely.

The increase in RBC was correlated with an increasing average carbon oxidation state of BC-containing particles (Fig. 1c), which was calculated from the SP-AMS measured atomic O:C and H:C ratios (oxidation state=2 × O:C—H:C)38. The atomic O:C and H:C ratios were calculated using the parameterization in Canagaratna et al.39 Particle oxidation state increases with atmospheric oxidation and was used as a proxy for particle photochemical age here (we do not have the necessary measurements to calculate photochemical age). The correlation of RBC with oxidation state suggests that photochemical processes over long timescales were likely important for producing both organic and inorganic coatings on BC.

Comparison of calculated and observed E abs

An important consideration is the extent to which the observations are consistent with Mie theory calculations of BC absorption, as core–shell Mie theory is used in many current climate models. The Eabs values were calculated using observationally constrained core–shell Mie theory (Methods) by either assuming that no (Eabs-calc-MAC) or some (Eabs-calc-TD) residual material remained both on the BC particles and as externally mixed particles after heating21. The Eabs-calc-MAC and Eabs-calc-TD were conceptually comparable to Eabs-obs-MAC and Eabs-obs-TD, respectively, assuming core–shell configuration for both ambient and thermodenuded BC-containing particles. An additional set of calculations were performed under the Rayleigh–Debye–Gans (RDG) approximation: The absorption by both uncoated and coated BC was determined by the properties of the individual spherules that made up a single BC particle40. The total absorption by a single particle is approximated as the sum of the absorption from individual spherules comprising a single particle, as calculated from spherical-particle Mie theory. As the size of the spherule was not explicitly known, and further might vary with the BC source, calculations were performed assuming spherules of 40 and 70 nm (diameter) to assess the sensitivity of the results to the assumed spherule size. To account for BrC contributions at 405 nm, calculations were performed using a range of imaginary refractive indices for the OM to find solutions that match the observed values. The base case assumed that OM internally mixed with BC had the same refractive index as that externally mixed from BC, and that the refractive index of the ambient OM was the same as the residual OM after heating. Alternative cases were also considered, and results of the calculations are shown in Fig. 2.

Figure 2

Comparison of observed and modelled Eabs,MAC and Eabs,TD at 405 nm.

Mie and RDG calculations (assuming BC spherule sizes of 40 and 70 nm) were performed using a range of kBrC values that are shown in the legend. kBrC indicates the cases in which a single kBrC value was assumed for OM; kBrC,ext and kBrC,int represent the kBrC values for OM that was externally and internally mixed with BC, respectively; kBrC,amb and kBrC,TD represent the kBrC values for ambient and thermodenuded OM, respectively. The error bars in (a–f) represent s.d. of the measured absorption enhancement values for each RBC interval.

Recalling that Eabs-obs-TD∼Eabs-obs-MAC at 781 nm, the calculated Eabs generally reproduced the observed trend with RBC (Fig. 1b). Comparison with the observed values at 781 nm showed that: (1) Eabs-calc-MAC overestimated observations by 14% when RBC was smaller than 2.9; (2) Eabs-calc-TD overestimated observations by 14% for RBC<1.4 and underestimated observations by 17% for RBC>3.8; and (3) Eabs-calc-TD was less dependent on RBC compared to Eabs-obs-TD. However, the differences between calculated and observed absorption enhancement were comparable to the measurement uncertainty (15%) of Eabs-obs-TD. The difference between modelled Eabs-calc-MAC and Eabs-calc-TD at 781 nm resulted from (theoretical) contributions of residual coatings to the absorption by thermodenuded particles. As the observations suggested negligible contributions of residual coatings, this indicates that the Eabs-calc-MAC overestimated the absorption enhancement due to the lensing effect for thermodenuded particles at a given RBC compared with the observations.

The Eabs-calc-MAC at 405 nm were lower than the observations when OM was assumed to be non-absorbing for both the Mie and RDG cases (Fig. 2), suggesting a role for BrC absorption. Direct absorption by the internally mixed BrC was very small compared with externally mixed BrC, which was primarily a reflection of the much larger mass concentration of externally mixed OM compared with BC-associated OM. This was assessed by comparing two cases, one with a bulk, campaign average imaginary refractive index (kBrC) of 0.008 and the other with kBrC=0 for the internally mixed OM while leaving kBrC=0.008 for the externally mixed OM. The difference in observed and calculated Eabs-MAC can be reconciled if the ambient OM was assumed to be slightly absorbing, with a bulk, campaign average kBrC=0.004 for the Mie case, kBrC∼0.008 for RDG(70 nm) case, and kBrC∼0.012 for the RDG(40 nm) case at 405 nm. The larger kBrC required for RDG (40 nm) and RDG (70 nm) compared with Mie resulted from the larger initial difference between the Eabs-MAC values for RDG (40 nm) and RDG (70 nm) when kBrC=0.

The Eabs-TD at 405 nm was well-simulated when the OM was assumed non-absorbing, which is inconsistent with the above finding that the OM needed to be absorbing at 405 nm to match the Eabs-MAC measurements. Alternatively, reasonable quantitative agreement between the observed and calculated Eabs-TD at 405 nm was obtained when the small amount of residual, low-volatility OM after thermodenuding was assumed to be more absorbing than the ambient OM. Specifically, good agreement (on average within 5%) was obtained for the Mie calculation when kBrC of thermodenuded OM (kBrC,TD) equals 0.032 at 405 nm, while kBrC=0.004. The MAC calculated for the BrC scales approximately linearly with kBrC, and thus these results indicate the residual OM was approximately 8 times as absorbing per mass as the ambient OM. This finding indicates that the OM components that evaporated were comparably much less absorbing, with an estimated imaginary refractive index of 0.001 at 405 nm assuming volume mixing and that the average residual OM was 8.5% of the total (measured by the HR-ToF-AMS). This result gives a ratio between the kBrC,TD and the kBrC for the evaporated components of ∼32, consistent with the results from the apportionment of above. For the RDG calculation, good observation-model agreement can be achieved when the residual OM was six times as absorbing as the ambient OM, that is, kBrC,TD=0.048 for RDG(70 nm) and kBrC,TD=0.072 for RDG(40 nm). The overall results indicate that there was a non-negligible lensing-driven absorption enhancement of BC in Detling, and BrC contributed substantially to the observed particle absorption, and thus as well to the observed Eabs at 405 nm.

Mechanistic insight of E abs from particle morphology

To develop a mechanistic understanding of the notable dependence of Eabs on RBC at Detling, single-particle images from scanning electron microscopy were analysed. Individual BC particles were classified visually based on their coating and morphologies as ‘embedded', ‘partly coated', ‘thinly coated', and ‘partially encapsulated and/or surface attached' (Fig. 3a), similarly to China et al.4142 The number fraction of each type was determined (Supplementary Table 3) and the averages are shown in Fig. 3b. (The electron microscopy analysis may exhibit some composition-dependent biases due to possible losses of semi-volatile materials; Methods.)

Figure 3

Morphology and statistics of BC-containing particles.

(a) Representative electron microscopy images of BC-containing particles collected at the Detling site for embedded (top left), partly coated (top right), thinly coated (bottom left) and partially encapsulated and/or surface attached (bottom right) BC particle types. The size of each panel is 1 μm by 1 μm. (b) Average number fraction of particle types with the colours indicating embedded (dark green), partly coated (green), thinly coated (black) and partially encapsulated and/or surface attached (purple) BC particle types for the high- and low-Eabs-obs-TD samples at Detling during ClearfLo and at Sacramento during CARES. The ‘D' and ‘A' in the parenthesis indicate thermodenuded and ambient samples, respectively. (c) Box plot of for the high- and low-Eabs-obs-TD single particle sampling periods in Detling. (d) Box plot of area equivalent diameter of the embedded BC particles for the high-Eabs-obs-TD samples (#1 and #2 in Supplementary Table 3) and low-Eabs-obs-TD samples (#3 and #4 in Supplementary Table 3) for the Detling measurements. In the box plots (c,d), the bounds of each box represent quartiles, and the whiskers extend to the most extreme data points that are within 1.5 interquartile range of the box. The horizontal bars inside the boxes represent the median values.

Two samples were selected for analysis when the Eabs-obs-TD at 781 nm was high and another two samples were selected during low-Eabs-obs-TD periods. The values for the high- and low-Eabs-obs-TD samples were 1.67±0.42 (s.d.) and 1.23±0.07 (s.d.), respectively (Fig. 3c). The embedded, partly coated, thinly coated and partially encapsulated and/or surface attached particle types had number fractions of 17%, 60%, 12% and 11% (out of 2187 particles) for the high-Eabs-obs-TD and 6%, 66%, 12% and 16% (out of 2017 particles) for the low-Eabs-obs-TD samples. Therefore, the largest difference between the high- and low-Eabs-obs-TD samples was that the high-Eabs-obs-TD samples contained ∼3 times more embedded BC particles. The morphological configuration of embedded BC particles is closest to the core–shell structure commonly used in models and it is theoretically expected to show the largest absorption enhancement. In addition, the embedded particles for the high-Eabs-obs-TD samples had significantly larger area equivalent diameter than those for the low-Eabs-obs-TD samples (Fig. 3d) with the difference ranging between 58 and 159 nm at the 95% confidence level, suggesting that the BC particles associated with high-Eabs-obs-TD underwent substantial coating. Moreover, the thermodenuded particles during the high-Eabs-obs-TD period (Supplementary Table 3) showed a substantially larger fraction (36%) of thinly coated BC particles and a smaller fraction (3%) of embedded BC particles compared with the simultaneously collected ambient sample (Fig. 3b), consistent with the presence of semi-volatile components in the coatings on BC particles.

A similar analysis for particles collected at the urban T0 site (14 km northeast of the Sacramento downtown, CA, USA) during the Carbonaceous Aerosols and Radiative Effects Study (CARES) campaign is also shown in Supplementary Table 3 and the fractions of each type are shown in Fig. 3b. The embedded fraction was substantially higher at Detling (12%) than in Sacramento during the CARES campaign (1%), where only small Eabs is observed even at relatively large photochemical age (up to ∼15 h)14. In addition, the convexity42, which is associated with the compactness of the BC-containing particles, was higher for Detling (0.82±0.005) than CARES (0.76±0.005), suggesting different particle morphology at the two sites. This result indicates that mixing state may play an important role in determining Eabs. However, comparable electron microscopy measurements are not available for the CalNex campaign14, where small Eabs values are observed even though similar SP-AMS measurements show substantially larger RBC compared to the Detling measurements. Therefore, the hypothesis that the morphological details of the internal mixing-state significantly influence the absorption enhancement needs additional evidence with further field studies covering a range of different environments.

Discussion

Our field results demonstrate clearly that coatings can substantially enhance light absorption by BC emitted from mixed fossil fuel and residential solid fuel combustion sources. Photochemical aging can play an important role in BC's absorption enhancement by coating BC particles with condensable materials. Comparison with results from other field campaigns indicates that the influence of coatings on BC absorption may be source and regionally specific. We hypothesize that mixed combustion sources and atmospheric transformations can produce favourable conditions of enhanced BC absorption by coatings, including BrC. Our analysis of absorption enhancement as a function of morphology, mixing state, coating thickness and the amount of low-volatility and semi-volatile BrC at wintertime Detling supports our hypothesis. We highlight that the effects of BrC, in particular the low-volatility BrC, on BC's absorption enhancement at short wavelengths should be considered in future studies. Our findings should be evaluated by future measurements of atmospherically processed BC particles from mixed combustion sources at other locations. We recommend that models treat the variability of BC's absorption to account for the dynamic aging processes and regional variations in sources by calibrating Mie models to measurements. Continued assessment of the absorption enhancement by BC coatings (including BrC) through comparing observations and observationally constrained calculations is necessary to develop robust numerical models, in order to accurately quantify the past, present and future contributions of absorbing particles to climate change.

Methods

Thermodenuder

A thermodenuder (TD, Aerodyne Research, Inc., Billerica, MA, USA) was installed upstream of all the instruments used in this study to evaporate semi-volatile particulate components. The flowrate through the TD was 2.3 l.p.m., resulting in a residence time of 5.3 s, which was shorter than the 8.5 s residence time in the study by Cappa et al.14

Ambient particles and heated particles were alternately sampled during 5–15 February 2012 with TD temperatures of 120 and 250 °C. Each ambient-TD cycle included 10-min ambient measurements and 10-min TD measurements. Unless stated otherwise, the Eabs values derived from 250 °C thermodenuding were reported.

Photoacoustic soot spectrometer

A three-wavelength photoacoustic soot spectrometer (PASS-3, Droplet Measurement Technologies, Inc., Boulder, CO, USA) was used to continuously measure the absorption and scattering coefficients of fine particles at 405, 532 and 781 nm. In the PASS-3, the particle laden sample is illuminated with a modulated laser. The energy absorbed from the incident light is thermally transferred to the air surrounding the particle and the subsequent air expansion produces a sound wave26. The sound wave is measured by a microphone that was calibrated at 532 nm at Droplet Measurement Technologies, Inc. using NO2 (ref. 43). The microphone calibration is independent of wavelength and the lasers at 405, 532 and 781 nm were calibrated using a laser power meter. The calibration was checked using size selected fullerene soot particles (Alpha Aesar, Ward Hill, Massachusetts) following the procedures described by Flowers et al.44 An NO2 scrubber was installed upstream of the PASS-3 to remove the interference of NO2 absorption at 405 nm. Background measurements of filtered air were made every 2.5 min to allow for correction of changes in gas composition. During the Detling measurements, low 532 nm laser power and drifts caused large variations of the signal in this channel. Therefore, was used as a qualitative measurement only and is not included here. Allan variance plots were made before the measurements to evaluate the precision and stability of the performance of the PASS-3 (refs 45, 46). The precision at 2.5 min integration time was 0.5 M m−1 for and 0.3 M m−1 for . Zero tests using high-efficiency particulate arrestance filters suggested that the gas-phase absorption was within the precision of the PASS-3 measurements. The upper limit of the gas-phase contribution to absorption was estimated to be 4 and 6% at 405 and 781 nm, respectively. The relative humidity inside the PASS-3 instrument was smaller than 30%, suggesting that the particles were dry by the time they were measured. The low relative humidity inside the PASS-3 instrument was due to warm up of the cold air brought into the warmer container and the relatively high temperature (20–30 °C) in the PASS-3 measurement cell induced by the electronics.

To verify the PASS-3 measurements, the absorption measured by the PASS-3 was compared with the absorption calculated from a 7-wavelength Aethalometer (MAGEE Scientific, model AE31) that operated simultaneously at the same site47. The PASS-3 measured absorption correlated well with the Aethalometer-derived absorption with Pearson's r≥0.94 at 405 and 781 nm. Major axis line-fitting suggests that the PASS-3 and Aethalometer measurements agreed well, with slopes close to 1 (0.934±0.004 at 405 and 0.918±0.004 at 781 nm) and intercepts slightly negative (−2.724±0.085 and −1.847±0.037 M m−1).

SP2

The refractory BC (rBC) mass concentration was quantified using an SP2 (Droplet Measurement Technologies, Inc., Boulder, CO, USA). Measurement principles of the SP2 have been described in detail in previous publications2728. In brief, the SP2 selectively measures the rBC-containing particles using a diode-pumped 1064, nm Nd:YAG laser on a single-particle basis. The rBC-containing particles absorb the energy emitted by the laser and are heated to their vaporization temperature and incandesce. The peak intensity of the incandescence signal is proportional to the mass of rBC, allowing for the quantification of rBC mass concentration. The SP2 measures single rBC-containing particles with an rBC mass range of ∼3–300 fg (ref. 27). The size of the rBC was calculated assuming a density of 1.8 g cm−3 and assuming spherical geometry with no internal voids48. Size-selected fullerene soot (Alfa Aesar, Ward Hill, Massachusetts; Stock# 40971, Lot# L18U002) was used for the SP2 mass calibration49. Contributions of rBC particles outside the detection window—which can vary with rBC source—were quantitatively estimated by fitting time-dependent rBC size distributions assuming a bimodal lognormal shape. The good correlation between absorption and [rBC] (Supplementary Fig. 4) suggests that the lognormal correction was sufficient to account for the rBC mass outside the SP2 detection limit.

HR-ToF-AMS

An Aerodyne HR-ToF-AMS (Aerodyne, Research Inc., Billerica, MA, USA) was deployed to measure the non-refractory organics, sulfate, nitrate, ammonium and chloride in submicron particles with 2-min time resolution. The HR-ToF-AMS configuration and principle of operation have been described in detail previously2930. Briefly, particles that are focused and accelerated by an aerodynamic lens impact a heated vaporizer (600 °C). The non-refractory components in the particles are flash vapourized and ionized, and the ionized fragments are measured by a high-resolution time-of-flight mass spectrometer.

SP-AMS

The SP-AMS (Aerodyne Research, Inc., Billerica, MA, USA) incorporates a 1064, nm Nd:YAG laser vaporizer (the same laser as used in the SP2 instrument) into the standard Aerodyne HR-ToF-AMS to detect rBC via laser-induced heating and vaporization31. During 5–15 February 2012, the tungsten thermal vaporizer was removed from the SP-AMS so that only rBC and the rBC-associated non-refractory particulate matter (including organics, sulfate, nitrate, ammonium and chloride) were detected. By comparing the ambient and thermodenuded OM that was associated with BC (normalized to [rBC] to account for particle losses in the heated TD), we find that on average 24% of OM that was internally mixed with rBC remained associated with the rBC particles after heating to 250 °C. This result is consistent with a previous study50 that shows large contribution (∼11%) of organic components to particulate residuals after particle heating at 300 °C.

Electron microscopy analysis

Single particles were collected on Nuclepore filters (100 nm pores) using an aspiration sampling technique. The sample collection time was 3–6 h depending on the particle concentration. For each collection time period, two filter samples, one representing ambient and the other representing heated (at 120 °C through a TD) particles, were collected simultaneously. This TD was different from the TD used for the PASS-3, HR-ToF-AMS and SP-AMS, and therefore it may remove a different amount of material. After collection, the samples were stored in the dark before analysis. The particle images were acquired using a field emission scanning electron microscope (Hitachi S-4700). We note that only particles larger than 100 nm were analysed to reduce uncertainties in particle classification. In addition, the number of embedded BC-containing particles was likely underestimated because (1) electron microscopy images only the surface features of the particles, thereby the BC-containing particles that were entirely engulfed by other materials and completely lost their original BC morphology could not be identified and (2) the electron microscopy analysis was conducted under vacuum, so evaporation of semi-volatile coatings on BC particles could not be avoided. Therefore, the number fraction for each particle type derived from the electron microscopy analysis was semi-quantitative.

Calculation of E abs-obs-TD

The TD method is described by the following equation

where and represent the absorption coefficients of ambient and thermodenuded particles, respectively; [rBC] is the mass concentration of rBC measured by the SP2; t represents the TD measurement, and t−1 and t+1 represent the times corresponding to ambient measurements before and after the TD measurement, respectively. Over 95% of the Eabs-obs-TD values were >1, with the values smaller than 1 resulting from temporal atmospheric variability during the sequential ambient and TD measurements.

Calculation of E abs-obs-MAC

In the calculation of Eabs-obs-MAC, the measurements were binned by the RBC intervals, and the MAC of BC was determined as the slope of the line-fitting of as a function of [rBC] for each measurement section (Supplementary Fig. 4). We examined the linear regression, major axis (also called orthogonal regression), and standardized major axis (also called reduced major axis) least square line-fitting methods. The linear regression method only accounts for the errors in the Y direction, while the other two methods take into account errors in both X and Y directions. The standardized major axis is the major axis calculated on standardized data. The standardized major axis is typically used when X and Y data are not measured on comparable scales51, which is appropriate for the babs versus [rBC] line-fitting in this study. In addition, comparison of the three fitting methods results showed that the slopes derived from standardized major axis were in between those derived from the linear regression and major axis fitting methods (Supplementary Fig. 2). An additional method was to take the average point-to-point babs/[rBC] as the MAC value for each RBC interval, defined as the ‘ratio method.' The ratio method did not give the same results as those derived from the slope methods. We suspect that this difference was due to a systematic offset in the PASS-3 measurements that affected the ratio method but not the slope methods. For these reasons, we used the slopes derived from the standardized major axis fitting as the MAC values in this study. Eabs-obs-MAC was calculated as the ratio of the calculated MAC to the reference MAC of uncoated BC.

Apportionment of the absorption at 405 nm

The measured was the sum of the absorption of BC, BrC (both internally and externally mixed with BC), and the enhanced absorption due to the lensing effect of coatings on BC. Upon heating at 250 °C, the semi-volatile OM and semi-volatile inorganic species evaporated, leaving only BC and low-volatility non-BC materials. Unless stated otherwise, the temperature threshold for separating semi-volatile and low-volatility components was defined as 250 °C in this study. can be apportioned using and (at 250 °C) based on the following assumptions: (1) BrC absorbs at 405 nm but not at 781 nm (ref. 24); (2) the absorption of the residual BC particles (after heating) follows the inverse wavelength dependence; and (3) to a first-order approximation the lensing-driven enhancement for BC particles is wavelength independent, that is, for the same BC-containing particle52. Note that the measured ambient included the lensing-driven enhancement and the enhancement induced by evaporation of BrCSV, and further was affected by BrCLV as discussed in the main text. Therefore is different from used in the equations below: excludes the effects of BrC (BrCLV and BrCSV) and represents the lensing-driven enhancement only. In comparison with , the measured was due to the lensing-driven enhancement only, that is, . The apportionment is described by the following equations with the variables described in Supplementary Table 1:

From equations (2, 3, 4), and were solved; was then derived using equations (5) and (6).

Similar apportionment of was performed using the for particles that were thermodenuded at 120 °C. In this case, the BrCLV was referred to as the fraction of BrC that remained in the particles after thermodenuding at 120 °C. The apportioned at 120 °C and 250 °C correlated in time with a Pearson's r of 0.8, and the latter was 19% lower on average, which was likely because of the absorption of BrC that evaporated between 120 and 250 °C. Since charring of OM is unlikely to take place at 120 °C, the correlation of apportioned at 120 and 250 °C suggests that at 250 °C was a real signal rather than an artefact resulting from pyrolyzed OM. This result is consistent with previous studies that confirm that the interference of pyrolyzed OM on the absorption measurements is <1.7% at 300 °C (ref. 53).

Non-fossil fractions of oxygenated organic carbon

Radiocarbon (14C) analysis was used to distinguish the fossil and non-fossil carbon fractions of the oxygenated organic aerosols retrieved from the HR-ToF-AMS analysis. Non-fossil carbon fraction (fNF) was calculated as the ratio of 14C/12C for the sample to 14C/12C in 1950 after correcting for field blanks and nuclear weapon tests in the 1950s and 1960s (ref. 54). The fNF for total carbon (fNF,TC) in Detling has been reported by Crilley et al.55, in which the total carbon is the sum of organic carbon and elemental carbon measured with a Sunset thermal-optical carbon analyser. fNF of BC (fNF,BC) in Detling was assumed to be the wood burning fraction of BC, apportioned using the Aethalometer measurements and a source apportionment model56 as reported in Crilley et al. The average fNF,TC and fNF,BC at Detling during ClearfLo are 64±7% and 30±13% (ref. 55), respectively.

Here, we expanded the Crilley et al. analysis to analyse the fNF for oxygenated organic aerosols using the method described in Zotter et al.54 The equation we used for the fNF calculation is as follows

where HOC, SFOC and OOC are the organic carbon concentrations for hydrocarbon-like organic aerosol, solid fuel organic aerosol, and oxygenated organic aerosol factors, respectively. The organic carbon concentration of the factor was converted from the corresponding factor OM concentration, using the OM to organic carbon ratios derived from the factor spectra. rBC is the ambient refractory BC concentration measured by the SP2. TC is the total carbon concentration, equivalent to the sum of HOC, SFOC, OOC, and BC concentrations. fNF,TC and fNF,BC values are the same as reported in Crilley et al.55 We used the fNF,TC and fNF,BC values for 5–12 February 2012 that overlapped the Eabs measurements. fNF,HOC and fNF,SFOC are the non-fossil carbon fractions for HOC and SOFC, respectively. Since fNF,HOC and fNF,SFOC were not explicitly known, we performed sensitivity analysis using a pseudo Monte Carlo method. In this method, fNF,HOC was varied between 0 and 0.3, representing its substantial fossil emissions and accounting for the potential contribution of non-fossil cooking and biofuel combustion (for example, biodiesel) emissions that were not separated from HOC; fNF,SFOC was varied between 0.9 and 1.0 to consider the dominance of wood burning and account for the contribution of residential (fossil) coal burning to SFOC.

The apportionment result (Supplementary Fig. 5) showed that the fNF,OOC ranged from 0.75 to 0.93, suggesting that the majority of oxygenated organic aerosols was from non-fossil sources in Detling during ClearfLo.

Mie and RDG model calculations

The Mie calculations were performed using the same Igor code (Wavemetrics Inc) as in Cappa et al.14, which was adapted from the Fortran code of Bohren and Huffman57. The input size distributions for the calculation were constrained by observations using the same set of measurements (from the same types of instruments) and methods as detailed in Cappa et al., that is, the core BC size distribution was constrained by the SP2 measurements, the coated BC-containing particle size distribution was determined by RBC and the core BC size distribution assuming that the coated components were evenly distributed around BC particles, and the non-BC particle size distribution was derived by subtracting the SP2-measured BC size distribution from the total particle size distribution measured by a scanning mobility particle sizer.

In the RDG calculation, the number of spherules per particle was calculated by the assumed spherule size and the SP2-measured spherical equivalent volumes of the overall BC core. It was assumed that the amount of coating on each spherule was proportionally the same as for the total particle, that is, that RBC was the same for the spherule as for the total particle. This assumption is equivalent to having the ratio between the coated diameter and the core diameter conserved between the standard Mie (total particle) formulation and the RDG formulation. The absorption enhancements for the RDG cases were calculated analogously to the Mie case.

The refractive index of BC was assumed to be n=1.88+0.8i (ref. 14). Since the imaginary refractive index of BrC (kBrC) is highly uncertain and can vary greatly, we performed calculations to estimate the kBrC that resulted in good agreement between predicted and measured Eabs-MAC and Eabs-TD at 405 nm. It was assumed that the refractive index was a bulk average for OM.

Additional information

How to cite this article: Liu, S. et al. Enhanced light absorption by mixed source black and brown carbon particles in UK winter. Nat. Commun. 6:8435 doi: 10.1038/ncomms9435 (2015).