Black carbon aerosol size in snow.

The effect of anthropogenic black carbon (BC) aerosol on snow is of enduring interest due to its consequences for climate forcing. Until now, too little attention has been focused on BC's size in snow, an important parameter affecting BC light absorption in snow. Here we present first observations of this parameter, revealing that BC can be shifted to larger sizes in snow than are typically seen in the atmosphere, in part due to the processes associated with BC removal from the atmosphere. Mie theory analysis indicates a corresponding reduction in BC absorption in snow of 40%, making BC size in snow the dominant source of uncertainty in BC's absorption properties for calculations of BC's snow albedo climate forcing. The shift reduces estimated BC global mean snow forcing by 30%, and has scientific implications for our understanding of snow albedo and the processing of atmospheric BC aerosol in snowfall.

BC, a byproduct of combustion, has a unique role as a climate forcing agent. It absorbs the most solar light per unit mass of any atmospheric aerosol material, as indicated by its large mass absorption cross section (MAC). The effect of anthropogenic BC aerosol on snow albedo and lifetime is the subject of enduring interest12 due to its consequences for regional and global climate forcing3. When deposited to snow, BC produces a definitive decrease in albedo45 that is believed to have directly contributed, for example, to the reduction of Himalayan glaciers6. BC's small size and initial tendency to hydrophobia enhances its propensity to be transported to remote regions, such as the Arctic7, where its impact on snow albedo is predicted to significantly alter the radiative balance of the planet3. Human activity produces a large fraction of emitted BC, which as an atmospheric pollutant has serious health8 impacts, making BC mitigation for both climate change and public health an important subject of scientific and political consideration. Hence, the need to understand the uncertainties and mechanisms of BC's climate interactions is great.

MAC is often used to quantify the extent to which BC absorbs solar radiation because it allows convenient translation between BC mass and its absorption and, as such, is a useful parameter for modelers and experimentalists. BC MAC in the atmosphere depends on the extent of internal mixing of the BC with other materials9, a relationship that is the focus of much research10. MAC also depends on the size of the BC particles because particles tend to absorb light in amounts proportional to their mass when their diameters are small with respect to the light's wavelength, and in the limit of large size, only in amounts proportional to their surface area. Atmospheric BC size distributions have only recently become commonly111213, available, and almost uniformly show that atmospheric BC is small enough to absorb light in amounts linearly proportional to its mass.

Generally, BC's atmospheric MAC values are adopted for use in snow. For example, measurements of light absorption in snow have been translated to BC mass concentration in snow via an assumed MAC value14 based on atmospheric BC, and an atmospheric BC MAC value has been used as the basis for a global model estimate of the radiative forcing of snow-bound BC3. In the results shown here, BC's size is found to be larger in snow than in the atmosphere, in part due to the influences of the processes associated with removal of BC from the atmosphere. Mie theory analysis indicates that BC's atmospheric MAC value may be inappropriately large for use in snow.

Results

Our results were obtained with a Single Particle Soot Photometer (SP2), an instrument that detects the refractory BC mass content and mixing state of individual aerosol particles15. The SP2 quantifies BC mass, not the physical size of the BC material, so each observed mass is translated to a volume-equivalent diameter (VED) by assuming a 1.8-g/cc void-free density. Hence, VED size distributions are referred to as “mass size distributions”. The SP2 can be used in conjunction with a nebulizer for determination of BC mass concentration and mass size distribution in liquid samples (in this case melted snow), even those with significant dust contamination, as discussed in detail elsewhere16, and briefly in the Methods section.

A system that characterizes BC mass size distribution in snowmelt will only give accurate information about BC size in snow if the melting process does not change the mass size distribution. Here, we assume that neither melting of the snow, nor nebulization of the resulting melt, causes BC agglomeration that would shift the BC distribution to larger sizes. This assumption is well supported: 1) the mass concentrations we observed in snow were low (<~100 ng-BC/g-snow), so the BC in snow inhabited less than 0.1 × 10−6 of the volume of water present and agglomeration in liquid was unlikely; 2) laboratory testing showed that the BC mass size distribution was fairly stable within melted snow over timescales of ~1 d, directly indicating low enough concentration to avoid significant agglomeration in the liquid phase on this time scale, and at least 24 times as long as was taken to perform the measurements discussed here; 3) laboratory tests also showed that individual thaw/refreeze cycles of laboratory BC standards in water only changed the BC mass size distribution slightly, in a manner suggesting that even when BC was “pressed” into smaller volumes by the exclusionary effect of ice crystal formation, only a small mass fraction (< 5%) agglomerated. This result is valid for the concentrations observed in the snow melt and suggests <5% shift in MAC for a single freeze/thaw cycle of fallen snow; and 4) laboratory BC standards in water revealed no sign (< 5% of total BC mass) of agglomeration even at much higher (x 50) concentrations than typical snow samples, validating the performance of the nebulizer.

Based on the observed small agglomeration effect in thaw/refreeze cycles, we hypothesize that, with repetition, such cycles occurring in ambient snow could significantly affect the BC mass size distribution. It follows that the first laboratory melting of a snow sample is expected to provide the most accurate size distribution of the BC contained in the snow at the time of sampling, and the earliest sampling after snowfall will provide the mass size distribution closest to that in the snow when deposited on the surface. For this reason, a single average mass size distribution of snow-bound BC is presented as determined only from optimally handled samples.

Figure 1 shows a mass size distribution averaged from five local snow samples collected from three Colorado snowfalls within 60 km of Denver, CO, in both semi-rural and rural areas (black line, with whiskers showing the standard deviation of the five individual distributions around the average). These were directly sampled within minutes of their first melt in the laboratory with the nebulizer/SP2 system (for up to 40 minutes). Also shown in Figure 1 are two mass size distributions for atmospheric BC. Although there is some variability in the size distribution of BC observed in the atmosphere with SP2s, the distribution is typically close to a pure lognormal function with a mass median diameter of ~220 nm VED or smaller. The distribution with the smallest BC median size in Figure 1 (red circles and fit, labeled “urban”) is from ambient atmospheric BC observed in the Los Angeles basin during an entire flight of the NOAA WP-3D research aircraft in the CalNex 2010 project (unpublished data obtained with the NOAA SP2). In this dataset, the BC median size was small enough that a significant fraction of the total accumulation mode mass could not be detected by the SP2, as shown by the lognormal fit. Less than 1% of the mass of this mode was associated with BC cores larger than 600 nm VED based on the fit. The next larger median size (blue line with squares, labeled “remote”) is associated with the distribution of atmospheric BC observed over the remote regions of the Pacific Ocean during the 2009 HIPPO campaign17. This size distribution is very similar to that measured in fresh biomass burning plumes observed in 2006 in Texas13 (not shown), and also has less than 1% of its mass associated with BC > 600 nm VED based on log-normal fitting. An exceptionally large atmospheric distribution was measured with an SP2 in a Southern China18. There, a secondary mode of atmospheric BC at larger size was observed in the polluted boundary layer in the region of a mega city, yet still only 6% of BC mass was found in sizes larger than 600 nm. This Chinese distribution probably provides an upper bound for the contribution of non-lognormal “shoulders” of BC at masses larger than the primary peak, as seen in the urban data of Figure 1. Urban-sourced atmospheric BC measured during the Texas campaign13 was intermediate in size between the atmospheric distributions shown. BC mass size distributions measured in the air in Colorado (not shown) also typically have this intermediate size, and hence fit comfortably within the small particle range where MAC is not strongly dependent on BC size.

Figure 1

BC mass size distributions and mass absorption cross section (MAC).

The normalized BC mass size distributions from the atmosphere (remote and urban) and from the five snow samples, as a single average with whiskers representing the variation in the underlying distributions (standard deviation) are plotted against the left vertical axis. The MAC is plotted with a dashed line against the right vertical axis, and was calculated for 460 nm light and BC with index of refraction (2.26, −1.26).

BC mass size distributions are not available in the air masses that generated the snowfall examined here. The snowmelt-BC mass size distribution has a primary mode fairly similar to the remote atmospheric distribution, but deviates from a log-normal distribution with a significant shoulder of mass in refractory BC particulate larger than ~500 nm VED, with 28% of its mass in BC larger than 600 nm VED. Included in this set of snow samples was a single measurement of extremely fresh snow, collected near our laboratory the morning after the snowfall (which occurred at night). This sample revealed that larger BC particles than typical were already present in the snow. It is unlikely that this snow experienced even a single thaw/freeze cycle before its laboratory evaluation; hence, this indicates that the shift in BC size occurred in the atmosphere, and likely was due to significant agglomeration and/or size selection of BC during the formation of the snow in the atmosphere. Note that although BC is dilute in snow, its concentration in freshly nucleated ice crystals or in freshly frozen particles in the atmosphere may be much higher. Thus, BC agglomeration via freezing of BC-containing droplets in the atmosphere may be more likely than via freeze cycles of fallen snow. We are not able to distinguish growth or selection of BC that occurred in the atmosphere from agglomeration that may have occurred after temperature cycles of fallen snow, or during dry deposition of BC independent of snowfall as in19.

In addition to the local snow samples, 31 snow samples from Arctic areas sampled far from sources were provided by Stephen Warren (Dept. of Atmospheric Science, Univ. of Washington, Seattle, WA USA) and Sarah Doherty (JISAO, Univ. of Washington, Seattle, WA USA). As discussed in the Methods section, although the Arctic mass size distributions are not directly presented here due to less than optimal handling, the majority of these samples showed some evidence of size distribution shift to larger sizes, although some did not; on average they revealed 17% of the BC mass from cores larger than 600 nm with a standard deviation of 24%, substantially more than seen even in the exceptional Chinese mass size distribution18 mentioned above. Hence these remote results support our basic conclusions based on the small set of local samples and suggest that even very remote snow may contain BC particles larger than expected based on atmospheric size distributions.

We have run simple calculations of particle coagulation and size-dependent removal to assess the likelihood that these processes could explain the large BC particles observed in snow. For high atmospheric concentrations of BC aloft in Colorado we modeled 200 accumulation-mode particles (without BC content) with 100 accumulation mode BC cores (log-normally distributed around number median diameter of 100 nm VED, equivalent to ~200 nm VED mass median diameter) in each cubic centimeter of air. Over twelve hours, a time scale that serves as an upper bound for aerosol interaction during precipitation events, only negligible changes in BC size distribution were predicted due to coagulation (≪ 1 nm shift in mass median diameter).

Size dependent removal of atmospheric BC has been observed13, in which BC mass distributions (based on the published number size distributions) had a median mass diameter of ~200 nm in air that underwent wet removal, in contrast with a mass median diameter of 219 nm in air that had not experienced wet removal (with geometric standard deviations of 1.51 and 1.64, respectively). The precipitated BC mass size distribution can be calculated from these values along with the reported ~80% removal efficiency that caused the shift in these distributions. This analysis shows that the precipitated BC would be shifted only slightly to larger size with mass median diameter ~220 nm and still less than 1% of its mass in BC cores larger than 600 nm. Hence, neither size dependent removal nor coagulation are likely to completely explain the larger BC found in snow. Hence, we hypothesize that one or more complex interaction are also at play, perhaps involving some combination of the following: 1) the agglomeration in the atmosphere of water-bound BC; 2) long-duration in-cloud residence times and prolonged exposure to collision-coalescence in a small fraction of cloud parcels (as has been modeled in Arctic stratus20) that could enhance BC content within individual droplets; 3) agglomeration of snow bound BC due to freeze cycles on the ground. Note that BC collected via impaction of snowflakes is not likely to agglomerate without repeated freeze cycles due to the extreme volume dilution in the ice (as discussed above). To test this hypothesis, both rigorous modeling and extensive measurements of BC size in snow of different ages will be needed.

The difference observed in BC's size range between the snow and the ambient atmosphere can cause a reduction of BC MAC in snow because larger BC particles tend to absorb light less efficiently per unit mass than smaller particles (Figure 1, dashed line). To quantify the magnitude of this reduction, MAC was calculated for the different BC size distributions shown in Figure 1 with Mie theory21. This approach assumes a simple picture of BC particles as spherules with values for density (1.8 g cm3) and complex index of refraction (n = (2.26, −1.26)) based on optimized index and density for optical modeling determined from ambient BC in Tokyo22. These values provide reasonable agreement between Mie theory scattering and that observed for fullerene soot, a calibration material with properties very close to that of ambient Tokyo BC and ambient BC observed in Switzerland2223. The Mie calculations were performed for light at 460 nm, the wavelength at which BC is estimated to most strongly affect snow albedo24. For the urban mass size distribution in Figure 1, the log-normal fit was used as the basis to calculate the MAC.

The sensitivities of the calculated MAC values to the index of refraction and light wavelength were explored as shown in Table 1. The MAC values from the remote atmosphere and first-melt snow distributions are shown as ratios to that from the urban distribution (shown separately in the second row). Substantial decreases in MAC are calculated for BC in the size range observed in snow compared to that seen in the atmosphere, up to a 40% decrease. In contrast, the atmospheric MAC values differ by only ~10%. The exceptional Chinese atmospheric BC size distribution18 discussed above had a MAC ~25% lower than that of the urban distribution. The differences between the typical atmospheric and first-melt snow MAC values are much larger than the differences that occur with differing choice of index and wavelength. Larger uncertainties associated with the Mie estimate of the MAC changes arise from the fact that the actual morphology of the BC particles in snow is not known. For example, if the BC were in highly fractal shapes, as it tends to be upon emission from fossil fuel sources, then it is likely that there would be no shift in MAC between snow and atmospheric BC (as in the Raleigh/Debye/Gans approximation for atmospheric BC25). This extreme is unlikely, as BC tends to collapse from a fractal to a denser morphology with increasing age and processing2627. Nevertheless, the Mie theory analysis presented here likely provides an upper estimate on the impact of BC size on snow-bound BC MAC for BC without optically significant coatings.

Table 1

MAC values and ratios calculated with Mie theorya

Index of refraction	(2.26, −1.26)	(2.26, −1.26)	(2, −1)	(2.26, −1.26)
Wavelength, nm	550	500	550	460
Urban MAC, m2 g−1	7.1	7.7	6.2	8.35
Remote/urban atmosphere ratio	0.97	0.92	0.94	0.88
First-melt snow/urban ratio	0.67	0.63	0.67	0.58

aThe ratios of refractory BC MAC calculated for remote or first-melt snow size distributions to that calculated from the urban size distribution (CalNex). The calculation was repeated with different values for absorbing wavelength and the index of refraction of bare BC spheres in Mie theory. The right-most column represents our best estimate.

Discussion

Increased BC size and decreased BC MAC in snow have implications for both modeling and measurement communities. BC MAC is often used as a fundamental model input, and has a strong effect on model estimates of BC's radiative forcing. BC's global mean snow albedo forcing, as calculated by the Snow, Ice, and Aerosol Radiative (SNICR) model coupled with a GCM3, is linearly dependent on BC MAC over a ~±20% range, with a constant of proportionality of 0.75. This indicates that an overestimate of BC MAC of 40% would result in a ~30% overestimate in the forcing. Determinations of BC mass loading in snow that depend on absorption measurements (e.g. 11) will underestimate the loadings if an atmospherically relevant MAC value is used.

This sensitivity highlights the need for new work to harmonize measurements with model treatments of BC's size transformation in snowfall, and to improve knowledge of both BC size and MAC, and their variability, in snow around the world, as well as the processes that influence them.

Methods

To use the SP2 to quantify particulate concentration and size distribution in liquid, the efficiency with which a nebulizer aerosolizes the various sized particulates in the liquid must be determined. We used a new technique16 to characterize this efficiency based on polystyrene latex sphere (PSL) concentration standards The PSLs were aerosolized from known concentrations in water, and the efficiency with which the nebulizer aerosolized them into a given volume of air was observed to depend on particle diameter, with larger diameter particles less well aerosolized. We assumed that losses scaled with particle stopping distance, not some other particle size/shape parameter. Thus, particle stopping distance, a measure of the inertia of particles that is often used in calculations of particle losses in tubing/sampling/etc., was used to connect the results from the PSLs to BC particles. Stopping distance S is defined as: where ρ is the density of the particle, η is the air dynamic viscosity, C is the Cunningham slip correction factor, and D is the particle diameter. BC stopping distance was calculated from the SP2-measured BC masses taken in conjunction with published BC density estimates2829, and the measured BC size distributions were scaled to normalize those efficiencies. As discussed in a technical paper16, for the nebulizer used here, these shifts did not result in large changes in distributions, but rather only in changes corresponding to a 20% shift in the estimate of BC mass loadings.

Due to the presence of large BC particle masses in snow melt, the sensitivity of the SP2 to larger particles than it typically quantifies in the ambient atmosphere had to be determined. This was accomplished, as discussed in full detail in the technical paper16, via an extrapolation of the SP2 calibration curve to large mass based on previous results162829. Briefly, other experimentalists28 have measured the non-linearity of SP2 response to BC mass for larger particles up to volume-equivalent diameters of 600 nm; the response is linear over most of the range where BC mass is typically observed in the air at smaller sizes. We calibrated the SP2 in the linear regime and then adopted the published non-linearity dependences to extend the calibration to large sizes, up to 2000 nm, for which calibration information is not available (and will be difficult to obtain). In terms of BC mass size distributions and MAC, discussed in the main text, the uncertainties associated with these manipulations for the specific nebulizer used in the study and the SP2 calibration lead to minor shifts in mass size distribution at large size and only a ~10% contribution to MAC uncertainty. Note that the SP2 is quite insensitive to dust in comparison to its sensitivity to BC. Examination of the vaporization temperatures of the large particles detected here confirmed their BC composition. The dust contribution influences BC determinations primarily at small sizes (<400 nm), and, for the samples used here, had an insignificant effect on the measured BC mass size distributions. Hence the overall uncertainty in the BC mass size distributions examined here are not large enough to substantially affect the underlying conclusions of this manuscript.

The majority of the 31 Arctic samples tested showed enhanced populations of large BC, consistent with the local snow results. However, these samples experienced two thaw/refreeze cycles before first sampling in the laboratory with the SP2, with 1 day or less total spent thawed before sampling by the SP2, and thus are not of optimal quality. To explore the effects of the additional freeze cycles and time spent thawed, we replicated this behavior on a local snow sample, and on laboratory standards of BC in liquid. The results indicated that, although the BC mass concentration in the samples were affected at the level of ~50%, the mass size distribution shifted smaller amounts, ranging from no shift at all (for the laboratory standards), to slightly less than half the amount observed in the arctic samples compared to atmospheric BC (for the local snow). Hence we include these samples as supporting the likelihood that BC size in snow can be larger than in the atmosphere over large regions of the globe.

Author Contributions

JPS performed the snow measurements and MAC calculations, and was the primary author of the manuscript; AEP performed the instrument calibrations and wrote the coagulation model code; JPS, AEP and JRS operated the instrument for the atmospheric measurements; RSG and DWF were part of the intellectual development of the work; all authors contributed to the manuscript.