ARTICLE IN PRESS Atmospheric Environment 40 (2006) 3574–3587 www.elsevier.com/locate/atmosenv Climate-relevant properties of primary particulate emissions from oil and natural gas combustion Tami C. Bonda,, Birgit Wehnerb, Antje Plewkab, Alfred Wiedensohlerb, Jost Heintzenbergb, Robert J. Charlsonc a Department of Civil and Environmental Engineering, University of Illinois at Urbana-Champaign, Urbana, IL, USA b Institute for Tropospheric Research, 15 Permoserstrasse, Leipzig, Germany c Department of Atmospheric Sciences, University of Washington, Seattle, Washington 98195, USA Received 18 August 2005; received in revised form 22 December 2005; accepted 27 December 2005 Abstract We report emissions of mass, light absorption, particle number, chemical composition and size-resolved organic species from an industrial boiler that burned natural gas and residual oil. Organic compounds detected from oil combustion are mainly alkanes; it is not a major source of identifiable polyaromatic hydrocarbons. Elemental carbon (EC) and organic carbon (OC) make up approximately 38% and 15% of the particles from oil burning, respectively. Mass emissions from natural gas were below detection limits. A number peak of ultrafine aerosol (diameters lower than 10 nm) was always associated with oil burning. Burning at full power produced the greatest number of particles in the accumulation mode. Natural gas also produced fine particles, but at a much lower rate. The emission rate of light-absorbing particles from this relatively new boiler is lower than that in current emission inventories. However, real-time measurements show a large contribution to emitted light absorption from boiler warm-up and transients, even those with very short durations. The measured absorption is best explained with a constant absorption cross-section for EC, rather than predictions based on size distribution or mixed aerosol; this finding is consistent with EC in fractal-aggregate form. We compare the emissions with those of a lignite stoker, which this boiler replaced during environmental cleanup in the mid-1990s. Emissions of mass, light absorption and particles are lowest from natural gas, but the oil boiler is also a substantial improvement: emissions of particulate matter are 100 times lower, and emitted absorption is three times lower. However, the oil-burning emissions have a greater net warming effect per mass than those of the lignite plant. r 2006 Elsevier Ltd. All rights reserved. Keywords: Particle size distributions; Light absorption; Elemental carbon; Organic speciation; Technology change 1. Introduction Anthropogenic aerosols from the burning of fossil fuels affect the solar input to the earth–atmoCorresponding author. Tel.: +1 217 244 5277; fax: +1 217 333 6968. E-mail address: [email protected] (T.C. Bond). 1352-2310/$ - see front matter r 2006 Elsevier Ltd. All rights reserved. doi:10.1016/j.atmosenv.2005.12.030 sphere system by scattering and absorbing radiation and by providing sites for the condensation of cloud drops. Emission rates of climate-relevant properties, such as number of particles and absorption crosssection, are used in simulations of net radiative and climatic effects of these particles. Particles that do not absorb light, such as sulfates, produce a cooling effect, especially in industrial regions (Charlson ARTICLE IN PRESS T.C. Bond et al. / Atmospheric Environment 40 (2006) 3574–3587 et al., 1991). Warming by light-absorbing aerosols offsets this cooling; even small amounts of absorption may significantly alter the net radiative forcing by particles (Haywood and Shine, 1995). (‘‘Radiative forcing’’, a measure of climatic impact, is the net change in flux at the Earth’s tropopause.) Black carbon (BC), produced from the incomplete combustion of carbon-based fuels, is thought to be the largest component of light-absorbing atmospheric particles (Rosen et al., 1978). Carbonaceous aerosols that do not absorb light, often called ‘‘organic carbon’’ (OC), may also make up a large fraction of anthropogenic particles (Novakov et al., 1997; Quinn et al., 2004). The relative amount of absorption versus scattering is one factor that determines whether those particles warm or cool the earth system (Jacobson, 2002; Penner et al., 2003). Anthropogenic increases in cloud condensation nuclei (CCN) may contribute to the ‘‘indirect effect’’, or negative forcing due to increases in cloud albedo (Twomey et al., 1984; Chuang et al., 2002). Primary particles appear to have a large impact on CCN (Adams and Seinfeld, 2003). Aerosol emissions from combustion sources also affect human health, and have been the subject of several epidemiological studies (e.g. Pope and Dockery, 1999). Although some researchers suggest that respiratory diseases are linked to the number concentration of particles below 100 nm (Pope, 2000; Oberdörster, 2001; Kreyling et al., 2004), source-characterization measurements typically focus on aerosol mass, which is dominated by larger particles. Direct measurements of the climate- and healthrelevant properties of emissions are rare. In most cases, mass emission measurements must be converted to the quantities of interest using assumptions about optical cross-sections or size distributions. Lack of information on the emissions translates directly to uncertainties in climate models. Measurements of climate-relevant properties at the point of emission provide initial constraints for model assumptions, even if those properties are altered in the atmosphere. In this paper, we report climate-relevant characteristics—particulate light absorption and number size distribution—of particles emitted from an industrial boiler that burns oil and natural gas. We provide size-dependent information on the chemical composition. While our work focuses on identifying the properties needed to estimate climatic impacts, the information may also be useful for studies on human health. 3575 The modern, computer-controlled boiler measured in this study replaced a lignite-burning stoker that was decommissioned in late 1996. The replacement occurred as part of the environmental cleanup following German reunification. We reported the properties of particulate matter emitted from the lignite stoker in earlier work (Bond et al., 1999b; Wehner et al., 1999). In this paper, we include a comparison with the old stoker to assess how emissions change due to the rapid changes in technology associated with economic development. 2. Method 2.1. Source description The boiler measured in this study provided heat and hot water for a research complex of about 20 buildings. The combustion was thought to be clean enough that no particulate control devices are needed, and none were installed. Natural gas was the preferred fuel due to its lower cost, but the heating plant received a discount on its gas prices in exchange for agreeing to use oil during peak load times. A large tank stored the oil for the boiler outside the heating plant. Table 1 contains analyses of the fuel used during testing. The building complex was not heated at night, so the boiler was either off or operating at very low power in the evenings. Our measurements each day began with the startup period and lasted until sufficient steam for the daily operation of the buildings was provided. The boiler typically operated for 7 h when burning natural gas, while the demand was met after 3 h each day with oil fuel. Table 1 Ultimate analysis of residual oil (DIN 51603), and composition of natural gas Residual Oil (DIN 51603) Natural Gas Element Mass fraction (%) Compound Mole fraction (%) Carbon Hydrogen Oxygen Nitrogen Sulfur Other 86 13 0.50 0.24 0.20 0.06 CH4 C2H6 C3H8 C4H10 CO2 N2 Other Data were provided by the City of Leipzig. 88.0 5.9 1.4 0.4 1.1 3.0 0.2 ARTICLE IN PRESS 3576 T.C. Bond et al. / Atmospheric Environment 40 (2006) 3574–3587 The oil boiler warmed up for about 1 h upon startup each day. During that time, the boiler operated at 100% power with 15–20% excess air, and the exhaust gas temperature oscillated between about 120 and 140 1C. After the warmup period, the excess air increased to 40–50%, the power dropped to 60–80% of maximum, and the exhaust temperature became more constant. When burning natural gas, the power behavior was the same, but there was no warmup period when higher excess air was used. 2.2. Measurement description The emission characteristics we measured that are directly relevant to climate modeling include the emission factor of particulate matter [EFPM, g (kg fuel) 1], the absorption cross-section of the emitted particles (MAC, m2 g 1), and the number size distribution of the particles. Values of EFPM and MAC can be combined into a single, easily measured parameter, the absorption emission index [EIabs, (m2 absorption) (kg fuel) 1] (Bond et al., 1998), bypassing the need to determine absorptive properties from chemical composition. An analogous quantity for particle number is the particle emission index [EIN, particles (kg fuel) 1]. Both EIabs and EIN can be obtained from real-time measurements of light absorption, number concentration, and CO2 in the exhaust gas. We sampled the exhaust gas with a porous-tube diluting probe, which reduced coagulation of particles, eliminated condensation of water, and simulated dilution in the atmosphere. The sample was drawn into the probe isoaxially, and was diluted with dried, filtered air in less than 10 ms; dilution ratios ranged from 10:1 to 30:1. Several studies have discussed the influence of dilution on particle number size distributions (Abdul-Khalek et al., 1999; Shi and Harrison, 1999), especially in the ultrafine ‘‘nucleation’’ mode. A recent study (Lyyrännen et al., 2004) suggested that the porous-tube diluter leads to the fewest artifacts in the nucleation mode, compared with other portable dilution systems. We measured particles with diameters below 1 mm, which interact most efficiently with visible light, by removing larger particles with two preimpactors from the diluted gas. Simultaneous measurements of light absorption by all particles smaller than 10 mm and by all particles smaller than 1 mm were statistically identical, so very little coarseparticle light absorption was ignored by choosing a 1-mm cutpoint. Sampling rates were below or near isokinetic. Submicrometer particles are little affected by variations from isokinetic sampling, and biases due to these variations were calculated to be less than 1% for these particles. Continuous measurements in the diluted gas included particulate light absorption, particle size distributions and CO2. We measured the diluted-air CO2 with an NDIR analyzer (Horiba MEXA GE534) and light absorption by particles with a Particle Soot Absorption Photometer (PSAP, Radiance Research, calibrated to a wavelength of 550 nm). We made spot measurements of CO2, O2, CO, and unburned hydrocarbons in the exhaust gas about every 5 min. We obtained dilution ratios by comparing CO2 measurements in the sample stream and exhaust. The power-plant computer provided charts of power and exhaust gas temperature. Particle number size distributions from 3 to 800 nm (particle diameter) were measured using a twin differential mobility particle size (TDMPS)system (Birmili et al., 1999). The TDMPS-system consists of two subsystems measuring different size ranges of dry particles at the same time. Very small particles (3oDpo21 nm) were measured by an ultrafine differential mobility analyzer (UDMA, Hauke-type, center rod length 11 cm) in conjunction with an ultrafine condensation particle counter (UCPC, Model 3025, TSI Inc., St. Paul, MN, USA). For larger particles (Dp421 nm), a differential mobility analyzer (DMA, Hauke-type, center rod length 28 cm) and a CPC (Model 3010, TSI Inc., St. Paul, MN, USA) were used. We collected samples on Nuclepores filters, which were weighed to determine mass emissions and analyzed with a four-wavelength integrating plate for light absorption efficiency. Values of EFPM, MAC, and wavelength dependence of absorption come from these measurements, which necessarily span larger time periods and are less numerous than the real-time measurements. Both PSAP and integrating-plate measurements were corrected according to the recommendations of Bond et al. (1999a). Size-segregated particle samples were collected by a five-stage low-pressure cascade Berner-impactor with cut-off diameters of 0.14, 0.42, 1.2, and 3.5 mm. This sample stream operated separately from the other measurements and therefore was not subjected to the pre-impactors. This sample was also diluted 1:10 by a dilution system similar to that used for the submicrometer particle measurements. Aluminum foils were used to collect the particles and to determine organic and elemental carbon. For ARTICLE IN PRESS T.C. Bond et al. / Atmospheric Environment 40 (2006) 3574–3587 3. Results and discussion 3.1. Particle number size distributions 3.1.1. Natural gas Fig. 1 shows representative size distributions measured during natural gas combustion, with full power decreasing power 60% of full power 107 102 106 100 105 10-2 104 103 ambientair 1 10 100 Dp in nm 1000 dV/ dlog Dp (µm3 cm-3) 108 dN/ dlog Dp (cm-3) analysis of single organic compounds, two pieces of ferromagnetic foils (Fe–Ni alloy, Curie point 590 1C) were placed on each of the four smallest stages. The mass of these samples was determined by weighing the foils before and after exposure (Mettler Toledo UMT 2, 50% RH). The absolute precision of this balance is in the range of 1–2 mg. The fraction of OC and elemental carbon (EC) was determined with a thermographic method (C-mat 5500, Ströhlein). A detailed description of the thermographic method is given by Neusüss et al. (2002). The aluminum foil was placed in a quartz tube and heated rapidly to a specific temperature. In a first step the sample is heated under nitrogen to 500 1C, and compounds which evaporate under these conditions are referred to as OC. In a second step, the sample is heated under oxygen to 650 1C, where all elemental carbon is oxidized to CO2 which is analyzed with an IR detector. No adjustment to account for the associated hydrogen and oxygen is reported here, so the total associated organic matter is expected to be higher than the OC reported by about 40%. For analysis of single organic compounds, we used the Curie-point pyrolysis gas chromatography mass spectrometry apparatus described by Neusüss et al. (2000). In this method, the ferromagnetic foil with the particles is heated very rapidly (about 20 ms) to the Curie point temperature by eddy currents in the foil stimulated by a high-frequency magnetic field. When the alloy reaches its Curie-point temperature the foil becomes non-magnetic and the heating effect ceases. Because the heating is rapid, the organic compounds are evaporated, not pyrolyzed, and introduced by a helium carrier into the GC/MS. Identification is carried out by matching the mass spectra and retention data with those of reference compounds. Deuterated polycyclic aromatic hydrocarbons and aliphatic hydrocarbons were used as reference materials. Only compounds which evaporate easily can be measured; thus, the identified mass is a small fraction of the organics. This method has the advantage that only small sample amounts (a few nanograms) are necessary for analysis, which allows size-segregated measurements of particle composition. 3577 10-4 Fig. 1. Number size distributions measured at the heating plant fired with natural gas (filled symbols) and at ambient air and the corresponding volume size distributions (open symbols) during varying operating conditions. concentration adjusted to the undiluted exhaust. These number size distributions were taken during full-power combustion, which always occurred during the first hour of operation; during 60% power after the warmup period; and during the transition between the two. With decreasing power, the number concentration of ultrafine particles (3oDpo20 nm) strongly decreases. The particle number in the size range 20–100 nm also decreases with decreasing power, but the concentration of particles larger than 100 nm is nearly constant during the entire measurement period. The number size distribution measured during the same time on the top of a building located next to the heating plant is included as the solid line. The concentration of particles larger than 20 nm is relatively low and similar to those values measured in surrounding the urban area. Number concentrations in the ultrafine regions at full-power operation are several orders of magnitude higher than the background concentration, and examination of all size distributions (not shown) suggests that the number concentration in the nuclei mode decreases approximately linearly with burning rate. 3.1.2. Fuel oil Size distributions with fuel oil performed at the same boiler differ significantly from those taken with natural gas, as shown in Fig. 2. Volume size distributions calculated using the assumption of spherical particles are shown as dotted lines with the same symbols. Fig. 2 shows three size distributions: one at full power, one at power decreasing between 100% and 60%, and one at 60% of full power. Like ARTICLE IN PRESS number concentration volume concentration 108 104 103 107 102 106 101 105 full power decreasing power approx. 60% of full power 104 1 10 100 100 10-1 1000 Dp in nm Fig. 2. Selected number size distributions measured at the heating plant fired with heating oil (filled symbols) and the corresponding volume size distributions (open symbols). those of Linak et al. (2000) for residual fuel oil, the volume mode peaks around 100 nm. The full-power scan differs from the others because it has higher concentrations in the accumulation mode range (80–200 nm). At lower power, the ultrafine mode increases slightly and the accumulation mode decreases significantly. During full-power conditions, the higher available surface area associated with the larger accumulation mode may remove some of the ultrafine particles by coagulation. Because the accumulation mode decreases as power decreases, the volume size distributions also depend significantly on the boiler power. 3.1.3. Comparison of three fuels The size distributions measured from residual oil combustion differ in both concentration and shape from those measured with natural gas. For oil, the number concentration maximum (3 108 cm 3) is about a factor of 30 greater than the highest value for natural gas. The maxima in the ultrafine region exceed those from full-power natural gas burning by an order of magnitude, and do not show the dependence on power that appears in the naturalgas combustion. The size distributions from oil also exhibit a much higher accumulation mode, by two to three orders of magnitude, than those from natural gas. The particle number size distributions measured during combustion of natural gas and heating oil show significant differences with respect to operating conditions. Fig. 3 shows how the total number concentration of particles varies with the boiler operating power for both fuels. The number concentration was calculated by integration over Particle number concentration [cm-3] dN/ dlog Dp (cm-3) 109 dV/ dlog Dp (µm3cm-3) T.C. Bond et al. / Atmospheric Environment 40 (2006) 3574–3587 3578 109 100% power 108 Natural Gas Heating Oil 107 106 105 60 - 70% power 104 0.0 0.5 1.0 1.5 2.0 2.5 time since start of operation [h] 3.0 Fig. 3. Change in total number of emitted particles (o800 nm) with decreasing boiler power during combustion of natural gas (dotted lines) and heating oil (solid lines). Boiler power trace is marked with filled squares. the measurement size range (3–800 nm), and is dominated by the ultrafine mode because of the high concentration in that size range. As expected from Fig. 1, the emitted particle number concentration decreases with decreasing power when the boiler is fired with natural gas, although the trend is not linear with power. During the warmup period, the particle number concentration decreases, while power values are constant. For the heating oil measurements, the number concentration remains nearly constant. Because of the slight increase in the ultrafine mode, the number concentration actually increases as the power drops. Fig. 4 shows the time series of the total volume concentration of the emitted particles (o800 nm). The particle volume was calculated from the number size distributions with the assumption of spherical particles. The solid curve (natural gas) decreases very slightly during the 2 h of measurement. The particle volume emitted at heating oil combustion decreases strongly after 1 h of operation when the power decreases and the excess air supplied to the combustion increases. Table 2 summarizes the characteristics of the lognormal modes for the two fuels, and also compares them with our previous study on a lignite-burning heating plant (Wehner et al., 1999). The measured number size distributions from different combustion processes were approximated by three lognormal modes, using an iterative moment-preserving least-squares fitting algorithm (Birmili, 2001). During this procedure, not only the number size distributions were fitted but ARTICLE IN PRESS T.C. Bond et al. / Atmospheric Environment 40 (2006) 3574–3587 12500 60 Natural Gas (left scale) Heating Oil (right scale) 10000 40 7500 5000 20 2500 0.0 0.5 1.0 1.5 2.0 2.5 time since start of operation [h] 3.0 Particle volume concentration [µm3 cm-3] Particle volume concentration [µm3 cm-3] simultaneously the distributions of the moments +2 and +4. This improved the convergence of the least-squares algorithm, and also ensured that the lognormal modes represented upper moments of the distribution, such as the condensational sink, and particle volume concentration at a similar accuracy as particle number concentration. Table 2 shows parameters of the approximated number size distributions at the different combustion stages, averaged over the entire test period. Each individual size distribution was taken over a period of 6–8 min. During the experiment at the old coal-fired heating plant, neither the boiler power nor the size distributions remained as constant as those in the new boiler. Therefore, three different size distributions, representative of most of the combustion, were selected and approximated by the lognormal modes. Fig. 4. Variation of particle volume concentration as a function of operation time for natural gas and oil combustion. The warmup period, with lower excess air, ends after 1 h of operation. 3579 The number concentration from the brown coal plant was nearly constant, about 2 107 cm 3. Oil combustion emits more particles than natural gas, with the number concentration varying between 5 107 and 2 108 cm 3. The emissions of the brown-coal heating plant have number maxima at particle diameters between 50 and 100 nm. The largest number concentration of natural gas and oil particle emissions occurs in the ultrafine range, especially during operation at full power. For oil combustion, the difference between full and halfpower emissions occurs in all modes, but most markedly in the accumulation mode. The number concentration drops and the peak diameter shifts to larger particles as the power decreases. For natural gas, the number concentration of the two modes that are significantly above ambient concentration decreases between full and half-power, but the ultrafine mode (around 10 nm) decreases most markedly. The particles are apparently larger at low power, but the difference is less noticeable than that exhibited by oil combustion. 3.2. Bulk emission properties Measured emission properties are summarized in Table 3. The uncertainties reflect both measurement error and one standard deviation in the observed variability; the latter contributed most of the uncertainty. Particulate emissions from the natural gas burner were quite low. The mass on the Nuclepores filters was below the detection limit, Table 2 Parameters of lognormal modes describing the measured number size distributions measured in the exhaust of different combustion processes at different fuels and different loads Source Ultrafine mode Aitken mode Accumulation mode Nt3 s3 Dm3 Nt2 s2 Dm2 Nt1 s1 Dm1 Natural gas Full power Decrease Steady state 1.7 107 1.1 105 8.4 103 1.3 1.6 1.6 5 8 13 2.5 104 1.7 104 1.4 104 1.8 1.8 1.8 40 52 59 1.5 103 1.4 103 1.3 103 1.7 1.6 1.7 206 206 198 Fuel oil Full power Decrease Steady state 1.7 108 3.4 108 1.2 108 1.3 1.4 1.7 5 6 7 1.5 107 9.5 105 4.1 106 1.8 1.6 1.8 36 60 35 1.2 107 1.2 105 2.6 105 1.6 1.4 1.6 89 200 150 Brown coal Scan1 Scan2 Scan3 1.6 105 1.0 104 2.5 105 1.4 1.4 1.4 3 5 3 6.8 107 4.9 107 7.5 107 1.5 1.5 1.4 48 51 52 1.5 105 4.5 104 1.2 105 1.6 1.6 1.6 400 300 280 Nt means the total particle number concentration in the current mode in part cm 3, s the geometric standard deviation, and Dm the mean diameter in nm. ARTICLE IN PRESS 3580 T.C. Bond et al. / Atmospheric Environment 40 (2006) 3574–3587 Table 3 Characteristics of submicrometer emissions from industrial combustion of natural gas and residual oil Natural gas EFPM (g kg 1) 0.00470.004 MAC(m2 g 1) 435 nm 525 nm 660 nm 800 nm Åap (absorption spectral dependence) n.d. n.d. n.d. n.d. n.d. EIabso1 mm (m2 kg 1) Integratedb, 525 nm Continuousb (Average) (excluding high event) (high event only) EINo1 mm (1015 kg 1) n.d. 0.002 7 0.002 n/a n/a 0.00370.001 Residual oil Warmup Normal Overall pa 0.0670.05 0.0370.02 0.0570.04 0.28 1.471.7 1.171.2 1.071.1 0.870.8 0.870.2 0.970.6 0.770.4 0.670.4 0.470.2 1.170.1 1.271.3 1.071.0 0.970.9 0.770.7 0.970.2 0.76 0.66 0.68 0.62 0.12 0.0470.03 0.2670.29 0.0570.05 0.6170.18 1.270.6 0.0270.01 0.0470.04 0.0370.05 0.0570.01 1.770.3 0.0370.03 0.1270.15 0.0470.05 0.2070.26 1.570.5 0.36 0.00 0.18 0.00 0.00 a The p-statistic is the result of an analysis of variance and represents the probability that the oil-combustion ‘‘warmup’’ and ‘‘normal’’ measurements do not differ. b Integrated values are from particles collected on Nuclepore filters and post-analyzed; continuous values are made with real-time PSAP measurements. and our estimated values of EFPM come from the Berner impactor measurements, which had a higher flow rate and a correspondingly higher mass. The emission of submicrometer particulate matter from natural-gas burning (0.003 g kg 1) is an order of magnitude lower than that for oil burning (0.05 g kg 1). For comparison, values published by the US Environmental Protection Agency (EPA) are 0.21 g kg 1 for natural gas and 0.5 g kg 1 for No. 5 oil (US Environmental Protection Agency (US EPA), 1996); the present measurements are substantially lower, possibly because of the recent installation. The value of MAC for oil emissions, 1.0 m2 g 1, is nearly an order of magnitude lower than that of pure BC (approximately 7.5 m2 g 1). Of course, the emitted particles are not wholly black carbon. Particle size and shape also affect absorption cross-section, and spherical particles that are large relative to the wavelength of light have lower absorption cross-sections. These particles cannot be considered small, but they are probably fractal aggregates rather than spherical; such particles would have absorption cross-sections similar to those of small particles. Thus, chemical composition and not morphology is probably the reason for the low value of MAC. This value is also highly variable, as shown in the table, implying wide variability in chemical composition. The absorption is nearly inversely proportional to the wavelength of light, as one would expect from black particles that are small relative to the wavelength of light (Bohren and Huffman, 1983) or fractal aggregates composed of small spherules. Values of the absorption emission index (EIabs) and particle number emission index (EIN) are calculated from 15-min averages of absorption, CO2, and DMPS measurements. These real-time data are more numerous and span more time than the filter measurements. For natural gas, the value of EIabs (0.002 m2 kg 1) was statistically different from that of filtered air, but not from the background air in the heating plant. For the oil boiler, the average value of EIabs is 0.12 m2 kg 1. As discussed previously, emitted absorption cross-section (EIabs) should be equivalent to the product of EFPM and MAC. Table 3 shows that this is not the case. The continuous measurements that yielded EIabs spanned one event with anomalously high absorption, while the filter measurements used to calculate EFPM and MAC did not include that period. The average value of EIabs in Table 3 is strongly affected by that period, and is not comparable to the product of EFPM and MAC for that reason. This anomalous event will be discussed later. ARTICLE IN PRESS T.C. Bond et al. / Atmospheric Environment 40 (2006) 3574–3587 500 [µg/m3] 400 300 200 100 (a) 3.5-10 1.2-3.5 0.42-1.2 0.42 0.14- 0.14 0.05- 0 particle diameter [µm] 5.0 [µg/m3] 4.0 3.0 2.0 1.0 3.5-10 1.2-3.5 0.42-1.2 0.42 0.14- 0.05- 0.140 0.0 particle diameter [µm] (b) Fig. 5. Particle mass concentration emitted during: (a) oil burning and (b) natural gas burning. in % of particle mass For comparison purposes, emission factors used in recent emission inventories of BC can be translated to values of EIabs by assuming a value of MAC. The value of MAC (7.5 m2 g 1) for pure BC can be used since we are comparing the emission factors for BC alone. EFBC for natural gas was given by Cooke et al. (1999) and Bond et al. (2004) as 0.011 and 0.002 g kg 1, respectively; these translate to EIabs of 0.08 and 0.015 m2 kg 1. Both are higher than our measured value by an order of magnitude or more. However, emissions from natural gas do not have a large impact on global inventories of absorbing aerosols because even the previous emission estimates were very low. For oil burning with no emission controls, 0.07 g kg 1 is given by Cooke et al. (1999) and 0.02 g kg 1 by Bond et al. (2004), equivalent to 0.52 and 0.15 m2 kg 1, respectively. The value in Table 3 lies between the two estimates, but only when the average includes the short period of anomalously high absorption. Excluding that period, the average lies even below the lower estimate. The newness of the boiler may be a cause of these low emissions. Fig. 5 shows the particle mass concentrations from natural gas and oil burning from the impactor measurements. As expected, the concentrations from oil combustion are much higher. For these emissions, most of the mass is found in the smallest particles, although there are still significant concentrations in larger particles. Eighty percent of the mass has diameters below one micrometer, and 90% has diameters below 3.5 mm. Emissions from burning natural gas exhibit very low mass concentrations, as shown in Fig. 5(b). Most of the emitted mass is found in the coarse particles. Fig. 6 shows the fraction of OC and EC in the particle mass from oil burning (natural gas emissions were below the detection limit). EC and OC represent about 40% and 5%, respectively, of the submicron particles. EC forms a substantial component of the fine particle mass. In contrast, the larger particles contain only small amounts of EC, confirming the observation that there was very little absorption by particles above 1 mm. The preferential concentration of EC in small particles agrees with the results of Kleeman et al. (1999). The high burning temperature contributes to the very low OC percentage for all particle sizes. For particles with diameters below 2.5 mm, both EC and OC fractions (0.38 and 0.08, respectively) are similar to those measured in an earlier study (0.29 and 0.13, respectively; Hildemann et al., 1991). Both studies 3581 70 60 50 40 30 20 10 0 OC EC 0.05-0.14 0.14-0.42 0.42-1.2 1.2-3.5 particle diameter [µm] 3.5-10 Fig. 6. OC and EC as percentages of particle mass in each size bin during oil burning. confirm that the particles are not purely carbonaceous. The remaining 54% of the particles probably includes water, mineral matter, ionic compounds, and some hydrogen and oxygen associated with the OC. ARTICLE IN PRESS 3582 T.C. Bond et al. / Atmospheric Environment 40 (2006) 3574–3587 Table 4 Speciation of organic single compounds for the four impactor stages Substance Biphenyle Dibenzofurane Methylphenanthrene or -anthracene Methylphenanthrene Methylphenanthrene Methylphenanthrene Cyclopenta[c,d]-pyrene oxygenated PAH 9H-Fluorenone Phenalen-1-one 9.10-Anthracendion 4H-Cyclopenta[def]phenanthrenone EPA-PAH Phenanthrene Anthracene Fluoranthene Pyrene Chrysene Benz[a]anthracene Benzo[b]+Benzo[k]fluoranthene Benzo[a]pyrene Total EPA-PAHs Alkanes n-Octadecane (C18) n-Nonadecane (C19) n-Eicosane (C20) n-Heneeicosane (C21) n-Docosane (C22) n-Tricosane (C23) n-Tetracosane (C24) n-Pentacosane (C25) n-Hexacosane (C26) n-Heptacosane (C27) n-Octacosane (C28) n-Nonacosane (C29) n-Triacontane (C30) n-Hentriacontane (C31) Total alkanes CPIodd Alkan-2-ones C14 C15 C16 Total Alkan-2-ones 0.05–0.14 mm 0.14-0.42 mm 0.42-1.2 mm 1.2–3.5 mm ng/mg ng/m3 ng/mg ng/m3 ng/mg ng/m3 3.4 5.2 1.5 1.6 2.4 0.7 0.7 2.3 22.1 0.1 0.4 0.1 n.d. 1.6 n.d. n.d. 0.2 n.d. n.d. n.d. n.d. n.d. n.d. n.d. 2.2 1.3 2.0 14.2 1.0 0.6 0.9 6.6 0.3 1.3 0.3 0.3 0.2 0.1 0.1 1.1 n.d. n.d. n.d. n.d. n.d. n.d. n.d. n.d. n.d. n.d. n.d. n.d. n.d. n.d. n.d. n.d. 7.7 0.5 2.1 6.0 3.6 0.2 1.0 2.8 3.2 1.2 n.d. n.d. 0.5 0.2 n.d. n.d. n.d. n.d. n.d. n.d. n.d. n.d. n.d. n.d. n.d. n.d. n.d. n.d. n.d. n.d. n.d. n.d. 8.0 1.5 5.9 15.4 1.7 2.3 2.0 3.4 40.3 3.7 0.7 2.7 7.2 0.8 1.1 0.9 1.6 18.7 7.0 1.0 4.9 22.1 0.7 1.5 1.0 2.5 40.7 1.2 0.2 0.9 3.8 0.1 0.3 0.2 0.4 7.0 2.1 n.d. n.d. 1.8 n.d. n.d. n.d. n.d. 3.9 0.3 n.d. n.d. 0.3 n.d. n.d. n.d. n.d. 0.6 5.2 n.d. n.d. n.d. n.d. n.d. n.d. n.d. 5.2 0.5 n.d. n.d. n.d. n.d. n.d. n.d. n.d. 0.5 3.8 3.8 8.6 10.5 18.6 23.0 32.5 24.7 16.2 8.8 5.5 5.5 1.5 0.3 163.2 1.7 1.8 4.0 4.9 8.6 10.7 15.1 11.5 7.5 4.1 2.6 2.5 0.7 0.1 75.9 0.9 3.4 3.7 9.1 10.7 17.0 23.8 21.8 21.8 17.0 6.8 3.8 3.6 1.6 1.0 145.0 0.6 0.6 1.6 1.8 2.9 4.1 3.8 3.8 2.9 1.2 0.7 0.6 0.3 0.2 25.0 1.0 0.9 0.9 0.4 2.1 0.3 0.3 0.1 0.7 n.d. n.d. n.d. n.d. n.d. n.d. 5.4 5.7 14.5 9.4 11.1 15.4 11.0 11.0 8.1 2.8 1.5 1.6 0.7 0.2 98.4 n.d. n.d. n.d. 0.8 0.8 2.2 1.4 1.6 2.3 1.6 1.6 1.2 0.4 0.2 0.2 0.1 0.0 14.6 0.9 n.d. n.d. n.d. ng/mg 6.2 6.8 16.1 9.4 10.8 14.5 10.1 10.7 10.4 4.4 3.2 2.9 1.5 1.6 108.8 n.d. n.d. n.d. ng/m3 0.6 0.6 1.5 0.9 1.0 1.3 0.9 1.0 1.0 0.4 0.3 0.3 0.1 0.1 10.1 0.9 n.d. n.d. n.d. n.d.: not detected. Values are per particle mass and per m3 of sampled air. Table 4 shows the mass size distribution for the identified polyaromatic hydrocarbons, alkanes and ketones for oil combustion. Less than 1% of the OC was identified for each size fraction, which is typical for this method. The highest concentrations (in ng mg 1and ng m 3 for all identified substances were ARTICLE IN PRESS T.C. Bond et al. / Atmospheric Environment 40 (2006) 3574–3587 found for the fine particles. At the high temperatures present in the boiler, the semi-volatile organic substances are primarily in the gas phase. During cooling of the hot exhaust, the substances condense mainly on the fine particles that provide the highest surface area. Identified compounds with the highest concentrations include pyrene, cyclopenta[c,d]pyrene, and several alkanes. Oil combustion is not a major source of identifiable polyaromatic hydrocarbons. The alkanes form the main part of the identified compounds. They result mainly from unburned fuel, and possibly from thermocracking of compounds during combustion (Kissin, 1990; Garrett et al., 2000). The concentration ratio of odd to even numbered n-alkanes, or the Carbon Preference Index (CPIodd), should be close to 1.0 for the combustion of fossil fuels (Simoneit and Mazurek, 1982; Simoneit, 1989). Our results agree with this prediction at all particle sizes. We are not aware of other organic speciation measurements for external-combustion oil boilers. The closest comparable measurements might be those of diesel trucks reported by Schauer et al. (1999). However, diesel engines burn a lighter distillate fraction of the fuel than reported here, and one which contains less mineral matter; additionally, the internal-combustion process in engines is unsteady and greatly differs from the external combustion in boilers. Nevertheless, we compare our speciation results with this study. The ratio between alkane mass and OC is 0.004 in our study, compared with 0.01 for diesel engines. Given the disparity in combustion and fuel, the agreement is fairly good. However, n-alkanes from oil burning are 0.02% of the total mass below 3.5 mm, compared 3583 with 0.2% for diesel engines. The greater disparity for the total mass comparison results because of the much smaller fraction of OC in our work (5% versus 19%). 3.3. Effect of operating conditions Table 5 separates averages for warmup and nonwarmup periods for the oil boiler. Average values of all characteristics were different under the two types of operation. In particular, the emission of light absorption was about six times higher during warmup. For values of EFPM and MAC, the difference is not statistically significant because of the small number of measurements relative to the high variability. Conventional wisdom states that an increase in excess air or equivalence ratio reduces the emission of ‘‘soot’’ or light-absorbing material; Goldstein and Siegmund (1976) observed this relationship in boilers burning heavy fuel oil. Our observations agree, although the correlation is not very significant. Increasing the burning rate has a slightly negative effect on the emission of absorption, probably because the higher temperatures result in better burnout of the carbon particles. This relationship appears valid on an instantaneous basis only, as the correlation is better for 1-s observations (correlation coefficient 0.44) and much worse for 1-min averages (correlation coefficient 0.10). Transient operation causes large increases in absorption emission. The most marked difference is that between absorption emission on the anomalous day (1 April) and on two other days. Fig. 7 shows a period of 10 min when the power was sharply decreased and then increased. The graph Table 5 Predicted light absorption at 550 nm based on chemical and physical properties, compared with measured light absorption Bap (M m 1) Error (%) All times Measured average Constant absorption cross-section Predicted from EC mass on impactor stages Mie theory applied to DMPS measurement (external mixture) Mie theory applied to DMPS measurement (internal mixture) 2200 1900 1400 3000 3800 16 57 +27 +42 Excluding high absorption period Measured average Mie theory applied to DMPS measurement (external mixture) Mie theory applied to DMPS measurement (internal mixture) 1200 1300 1700 +8 +29 ARTICLE IN PRESS T.C. Bond et al. / Atmospheric Environment 40 (2006) 3574–3587 3584 100 80 40000 60 40 20000 20 0 8.2 8.3 8.4 8.5 8.6 Hour of Day 8.7 Power Potentiometer (%) Bap (exhaust, Mm-1) 60000 0 8.8 Fig. 7. Behavior of absorption (Bap) during power turndown. Bap is a commonly used atmospheric measurement, given in units of meters squared absorption per cubic meter of air, multipled by 106. shows the setting of the power controller from the plant charts; the event is also reflected in the exhaust CO2 concentration. The absorption efficiency, MAC, of the particles from that test is higher than that from other tests, and the accumulation mode is larger. If that result were excluded, the average MAC from all oil-boiler tests would be about 30% lower, and the variability would be reduced by about 50%. An additional unknown factor caused the difference between the anomalous period on 1 April, when EIabs during warmup was 12 times higher than during the other days. The pre-warmup conditions were the same: natural gas burned the day before, with the boiler cooled for 1–2 h prior to initiating oil combustion. Power outputs, exhaust temperatures, and excess air matched closely, although the exhaust temperature on 1 April was cooler by about 5 1C. We also discount errors in our measurement system as a possible explanation, because emissions after the warmup period were comparable to other days. US EPA reports cite poor fuel/air mixing, decreases in combustion intensity or temperature, or wall quenching as causes of incomplete combustion (US Environmental Protection Agency (US EPA), 1996). Our measurements also concur with previous findings of within-source variability, implying that a number of tests, or continuous monitoring, are needed to adequately represent emissions (US Environmental Protection Agency (US EPA), 1997). Extensive realtime monitoring on a few sources may assist in understanding the variability of combustion sources and in identifying the number and type of measurements required to adequately characterize each source. 3.4. Light absorption and chemical composition The chemical composition of submicrometer particles is used to determine the climatic impacts of those particles. We now determine the combination of particle size information and chemical composition required to estimate particulate light absorption. Table 5 shows several methods of estimating light absorption and compares these estimates with the submicrometer absorption that is actually measured. We use the following methods: (1) applying a constant absorption cross-section; (2) applying Mie theory to the impactor measurements; (3) applying Mie theory to the number size distributions measured by the DMA, assuming an external mixture; and (4) applying Mie theory to the DMA measurements assuming an internal mixture. More details and results of these comparisons are given below. We used Matlab implementations of BHMIE and BHCOAT for the Mie theory (Bohren and Huffman, 1983; Mätzler, 2002). (1) We apply a constant absorption cross-section (7.5 m2 g 1, Bond and Bergstrom, 2006) to the measured mass of elemental carbon. The predicted absorption differs from the measurement by only 16%. (2) We apply Mie theory to the size distribution of EC on the submicrometer impactor stages. These estimates of absorption are lower than measurements, probably because the larger particles are predicted to absorb less light but in fact are fractal aggregates which have a constant absorption per mass (Bond and Bergstrom, 2006). (3) We apply Mie theory to the DMA size distributions, assuming EC occurs as pure particles in isolation, and that the EC fraction on the impactor stages represents the fraction of particles composed of EC. Despite the fact that the DMA misses the largest particles and thus a substantial fraction of the mass, this prediction is much greater (27%) than the measured absorption. (4) We apply Mie theory to the DMA size distributions, assuming that all particles of each size are identical to each other and contain cores of EC coated with non-absorbing material. Again, we use the fraction of EC from the impactor stages. This method predicts even higher absorption than method 3, and thus there is a greater discrepancy between prediction and measurement. This is true even when the transient period is excluded, eliminating the possibility that the discrepancy could be ARTICLE IN PRESS T.C. Bond et al. / Atmospheric Environment 40 (2006) 3574–3587 attributed to changing size distributions. The DMA is known to overestimate particle mass for combustion aerosol because fractal aggregates have greater drag than spheres of equivalent mass. Thus, it would also overestimate absorption. Incorporating recent research that would define such a correction will be the subject of future work. 3.5. Environmental impact of boiler replacement Table 6 compares emission characteristics between the old, lignite-burning plant and the new heating plant burning oil and gas fuel. Clearly, natural gas is the preferred fuel for clean burning; emissions of mass, light absorption and particles are all orders of magnitude below the other two fuels. The oil boiler is also a substantial improvement over the lignite plant. Compared with the brown coal stoker, emissions of particulate matter per mass of fuel burned are 100 times lower, and emitted absorption is almost three times lower. More particles per mass of fuel are emitted from oil than from brown coal because of higher concentrations in the ultrafine mode. Furthermore, the heating value of the fuel is four times higher for oil than for brown coal, and over six times higher for natural gas. The net result is an improvement in all measures when normalized to energy produced. Compared with brown coal burning, oil burning has 12 times lower absorption cross-section per unit energy, and slightly lower particle emissions. Natural gas improves over brown coal by three orders of magnitude. These improved emissions result not only from fuel-switching, but from the addition of modern technology, such as computer control of the fuel- and air-flow rates. 3585 The net climatic impact of particulate emissions is determined not only by the amount of light absorption, but also by the amount of light scattering by these particles. We do not have measurements of light scattering, but we can infer that the submicrometer mass has similar scattering properties for both oil and lignite emissions, because the mass is dominated by the accumulation mode. Compared with emissions from the lignite plant, the oil-burning emissions have a higher mass absorption cross-section and thus a greater net warming effect per mass. This presents a dilemma: the new boiler has decreased the amount of emitted absorption, but increased the warming by the particles that are emitted. We suggest that despite the net warming per mass, the reduction in emission of all climate-relevant properties is a positive step. In this work, we have quantified the actual change in health- and climate-relevant emissions due to one instance of targeting a ‘‘dirty’’ boiler for replacement. Continuing such replacement efforts will lead to significant benefits for air quality and health. The impact of these replacements will reduce humans’ impact on climate, but may lead to net warming. Acknowledgements We gratefully acknowledge the technical support by Technoserv-Center Dr. Modes & Partner OHG, Leipzig, in operating the heating plant. TCB acknowledges travel support from a gift from the Ford Motor Company to the University of Washington. References Table 6 Comparison of submicrometer emissions for old and new heating plants Lignite EFPM EIabs EIN g kg g GJ 1 4.7 520 1 m2 kg m2 GJ 1 0.32 36 1 1015 particles kg particles GJ 1 1 0.4 44 Oil Natural gas 0.05 1.2 0.003 0.05 0.12 2.9 0.0015 0.025 1.5 37 0.003 0.05 Values are emissions per kg of fuel burned or per GJ released. Abdul-Khalek, I., Kittelson, D., Brear, F., 1999. The influence of dilution conditions on diesel exhaust particle size distribution measurements, SAE Paper 1999-01-1142. Adams, P.J., Seinfeld, J.H., 2003. Disproportionate impact of particulate emissions on global cloud condensation nuclei concentrations. Geophysical Research Letters 30 (5), 1239. Birmili, W., 2001. A moment-preserving parametrisation of wide particle size distributions. Journal of Aerosol Science 32, S191–S192. Birmili, W., Stratmann, F., Wiedensohler, A., 1999. Design of a DMA-based size spectrometer for a large particle size range and stable operation. Journal of Aerosol Science 30 (4), 549–554. Bohren, C.F., Huffman, D.R., 1983. Absorption and Scattering of Light by Small Particles. Wiley, New York. ARTICLE IN PRESS 3586 T.C. Bond et al. / Atmospheric Environment 40 (2006) 3574–3587 Bond, T.C., Bergstrom, R.W., 2006. Light absorption by carbonaceous particles: an investigative review. Aerosol Science and Technology 41 (1), 27–47. Bond, T.C., Charlson, R.J., Heintzenberg, J., 1998. Quantifying the emission of light-absorbing particles: measurements tailored to climate studies. Geophysical Research Letters 25 (3), 337–340. Bond, T.C., Anderson, T.L., Campbell, D., 1999a. Calibration and intercomparison of filter-based measurements of visible light absorption by aerosols. Aerosol Science and Technology 30, 582–600. Bond, T.C., Bussemer, M., Wehner, B., Keller, S., Charlson, R.J., Heintzenberg, J., 1999b. Light absorption by primary particle emissions from a lignite burning plant. Environmental Science and Technology 33, 3887–3891. Bond, T.C., Streets, D.G., Yarber, K.F., Nelson, S.M., Woo, J.-H., Klimont, Z., 2004. A technology-based global inventory of black and organic carbon emissions from combustion. Journal of Geophysical Research 109, D14203. Charlson, R.J., Langner, J., Rodhe, H., Leovy, C.B., Warren, S.G., 1991. Perturbation of the northern hemisphere radiative balance by backscattering from anthropogenic sulfate aerosols. Tellus 43AB (4), 152–163. Chuang, C.C., Penner, J.E., Prospero, J.M., Grant, K.M., Rau, G.H., Kawamoto, K., 2002. Cloud susceptibility and the first aerosol indirect forcing: sensitivity to black carbon and aerosol concentrations. Journal of Geophysical Research 107 (D21), doi:10.1029/2000JD000215. Cooke, W.F., Liousse, C., Cachier, H., Feichter, J., 1999. Construction of a 11 11 fossil fuel emission data set for carbonaceous aerosol and implementation and radiative impact in the ECHAM4 model. Journal of Geophysical Research 104 (D18), 22137–22162. Garrett, R.M., Guenette, C.C., Haith, C.E., Prince, R.C., 2000. Pyrogenic polycyclic aromatic hydrocarbons in oil burn residues. Environmental Science and Technology 34 (10), 1934–1937. Goldstein, H.L., Siegmund, C.W., 1976. Influence of heavy fuel oil composition and boiler combustion conditions on particulate emissions. Environmental Science and Technology 10 (12), 1109–1114. Haywood, J.M., Shine, K.P., 1995. The effect of anthropogenic sulfate and soot aerosol on the clear sky planetary radiation budget. Geophysical Research Letters 22 (5), 603–606. Hildemann, L.M., Markowski, G.R., Cass, G.R., 1991. Chemical composition of emissions from urban sources of fine organic aerosol. Environmental Science and Technology 25 (4), 744–759. Jacobson, M.Z., 2002. Control of fossil-fuel particulate black carbon and organic matter, possibly the most effective method of slowing global warming. Journal of Geophysical Research 107 (D19), doi:10.1029/2001JD001376. Kissin, Y.V., 1990. Acyclic components in dewaxed heavy distillates. Fuel 69 (10), 1283–1291. Kleeman, M.J., Schauer, J.J., Cass, G.R., 1999. Size and composition distribution of fine particulate matter emitted from motor vehicles. Environmental Science and Technology 34 (7), 1132–1142. Kreyling, W.G., Semmler, M., Möller, W., 2004. Dosimetry and toxicology of ultrafine particles. Journal of Aerosol Medicine 17, 140–152. Linak, W.P., Miller, C.A., Wendt, J.O.L., 2000. Comparison of the particle size distributions and elemental partitioning from the combustion of pulverized coal and residual fuel oil. Journal of the Air and Waste Management Association 50, 1532–1544. Lyyrännen, J., Jokiniemi, J., Kauppinen, E.I., Backmann, U., Vesela, H., 2004. Comparison of different dilution methods for measuring diesel particle emissions. Aerosol Science and Technology 38, 12–23. Mätzler, C., 2002. MATLAB Functions for Mie Scattering and Absorption. Institut für Angewandte Physik, Bern, Switzerland. Neusüss, C., Pelzing, M., Plewka, A., Herrmann, H., 2000. A new analytical approach for size-resolved speciation of organic compounds in atmospheric aerosol particles: methods and first results. Journal of Geophysical Research 105, 4513–4527. Neusüss, C., Gnauk, T., Plewka, A., Herrmann, H., 2002. Carbonaceous aerosol over the Indian Ocean: OC/EC fractions and selected speciations from size-segregated onboard samples. Journal of Geophysical Research 107, D19, 8031, doi:10.1029/2001JD000327. Novakov, T., Hegg, D.A., Hobbs, P.V., 1997. Airborne measurements of carbonaceous aerosols on the East Coast of the United States. Journal of Geophysical Research 102 (D25), 30023–30330. Oberdörster, G., 2001. Pulmonary effects of inhaled ultrafine particles. International Archives of Occupational and Environmental Health 74, 1–8. Penner, J.E., Zhang, S.Y., Chuang, C.C., 2003. Soot and smoke aerosol may not warm climate. Journal of Geophysical Research 108 (D21), 4657. Pope, C.A., 2000. What do epidemiologic findings tell us about health effects of environmental aerosols? Journal of Aerosol Medicine-Deposition Clearance and Effects in the Lung 13 (4), 335–354. Pope, C.A., Dockery, D.W., 1999. Epidemiology of particle effects. Air Pollution and Health 31, 673–705. Quinn, P.K., Coffman, D.J., Bates, T.S., Welton, E.J., Covert, D.S., Miller, T.L., Johnson, J.E., Maria, S., Russell, L., Arimoto, R., Carrico, C.M., Rood, M.J., Anderson, J., 2004. Aerosol optical properties measured on board the Ronald H. Brown during ACE-Asia as a function of aerosol chemical composition and source region. Journal of Geophysical Research 109, D19S01. Rosen, H., Hansen, A.D.A., Gundel, L., Novakov, T., 1978. Identification of the optically absorbing component in urban aerosols. Applied Optics 17 (24), 3859–3861. Schauer, J.J., Kleeman, M.J., Cass, G.R., Simoneit, B.R.T., 1999. Measurement of emissions from air pollution sources. 2. C1 through C30 organic compounds from medium duty diesel trucks. Environmental Science and Technology 33, 1578–1587. Shi, J.P., Harrison, R.M., 1999. Investigation of ultrafine particle formation during diesel exhaust dilution. Environmental Science and Technology 33, 3730–3736. Simoneit, B.R.T., 1989. Organic matter of the troposphere. 5. Application of molecular marker analysis to biogenic emissions into the troposphere for source reconciliations. Journal of Atmospheric Chemistry 8 (3), 251–275. Simoneit, B.R.T., Mazurek, M.A., 1982. Organic matter of the troposphere. 2. Natural background of biogenic lipid matter ARTICLE IN PRESS T.C. Bond et al. / Atmospheric Environment 40 (2006) 3574–3587 in aerosols over the rural western United States. Atmospheric Environment 16, 2139–2159. Twomey, S.A., Piepgrass, M., Wolfe, T.L., 1984. An assessment of the impact of pollution on global cloud albedo. Tellus B 36, 356–366. US Environmental Protection Agency (US EPA), 1996. Compilation of air pollutant emission factors. Report AP-42. 3587 US Environmental Protection Agency (US EPA), 1997. Procedures for Preparing Emission Factor Documents. Report # EPA-454/R95-015. Wehner, B., Bond, T.C., Birmili, W., Heintzenberg, J., Wiedensohler, A., Charlson, R.J., 1999. Climate-relevant particulate emission characteristics of a coal-fired heating plant. Environmental Science and Technology 33, 3881–3886.
© Copyright 2026 Paperzz