Characterization of Non-methane Hydrocarbons Emitted from Various Cookstoves Used in China

Similar documents
Chapter URL:

ANNEX A: Methodology for Estimating Emissions of CO 2 from Fossil Fuel Combustion

Air Pollutant Emissions from a Four-Stroke Motorcycle Engine Influenced by Gasoline Aromatic Content

Prediction of NO X Emissions

Numerical modeling of multi-component fuel combustion using reactivity adjustment of a chemical reaction mechanism

RESEARCHES ON THE INFLUENCE OF PRESSURE WAVE COMPRESSOR ON THE INTAKE AIR TEMPERATURE AT THE SUPERCHARGED ENGINES

631 ALLOWANES AND TOLERANCES

Shanghai University of Electric Power (SUEP)

COTTON VARIETY FIBER CHARACTERISTICS AND YIELD COMPARISON. King Ranch Farms, Kleberg County, 1998

Report of Progress 961

Department of Technology, Yinjian Automobile Repair Co., Ltd., Beijng, China

Emission Factors for Gas-Powered Vehicles Traveling Through Road Tunnels in São Paulo, Brazil

TECHNICAL INSIGHT INTERNAL CLEARANCE - TYPES AND NORMS

1. Energy. Reference approach Sectoral approach Difference Explanation for difference as reported in table 1.A(c) of the CRF

Economic Contribution of Off-Highway Vehicle Recreation in Colorado

Appendix A: Lumped Parameter Assumption Using Biot Number

Unit 2. DC Motors and Generators UNIT OBJECTIVE DISCUSSION OF FUNDAMENTALS. Operating principle of dc motors

As an example of the type of analyses included

Design and Implementation of the Performance Test System for Battery-Powered Motorcycles

Travel Effects of a Suburban Commuter Carsharing Service

Applied Energy 87 (2010) Contents lists available at ScienceDirect. Applied Energy. journal homepage:

Analysis of Recoverable Exhaust Energy from a Light-Duty Gasoline Engine

JUMO MIDAS C18 SW. OEM seawater pressure transmitter. Applications. Brief description. Customer benefits. Special features. Approvals/approval marks

Research and Implementation of Tractor Power Shift Clutch Control System

Spark Plug Fouling 3.0 Liter OptiMax

Limitations of Using Skid Number in Accident Analysis and Pavement Management

AIRCRAFT ENGINES SERVICE INSTRUCTION ENGINE START AT LOW TEMPERATURES AT ROTAX SI R1

Policy Analysis. Environmental Implications on the Oxygenation of Gasoline with Ethanol in the Metropolitan Area of Mexico City.

Safety Relay Unit. Ordering Information. Ultra Slim Safety Relay

APPLICATION OF MASS MEASUREMENT TECHNIQUE IN CALIBRATION AND CONTROL FOR GAS CYLINDERS

11 263/ Radial Piston Pump Model R4 /Series 2X Fixed Displacement. Ordering Code 1 PF 1 R4 2X R Replaces:

Brazed heat exchanger XB

Policy Analysis. Estimate of Mercury Emissions to the Atmosphere from Petroleum. EPA Estimates of Mercury in Crude Oil.

Series 250 Type and Type Pneumatic Control Valves Type 3251 Globe Valve

DPST-NO None 2 channels Auto-reset Inverse 24 VAC/VDC G9SB-2002-A 4 1 channel or 2

Installation and Operational Instructions for ROBA -slip hubs Sizes 0 12 (B.1.0.GB)

TECHNICAL REPORT Kathleen T. Williams, PhD, NCSP Overview Revisions in the New Edition

Study: Draft Final Report

JUMO MIDAS C08. OEM Pressure Transmitter Basic. Applications. Brief description. Special features. Customer benefits

Lincoln Under-cabinet Light Installation Instructions

Power electronic topologies for wind turbines

The Study on Mechanical Reliability Design Method and Its Application

Effectiveness of Alternative Detector Configurations for Option Zone Protection on High-speed Approaches to Traffic Signals

Create a pe ph diagram for the following solid and aqueous iron species:

Low Pressure Domestic Burners

Cable Trolleys for C-Rails Program

Effect of minor components in soy biodiesel on cloud point and fi lterability

330mm (13") MAXIMUM WALL THICKNESS IN AREA OF UNIT. 330mm (13") MAXIMUM WALL THICKNESS IN AREA OF UNIT TYPICAL ACCESS TO LEVELING 15

POWER TRIM. Table of Contents. Section 5D - Auto Trim II

POWER TRIM. Table of Contents. Section 5D - Auto Trim II

(mω) 1. Temperature ( C) Voltage (V) ULRG1 /

Application Examples of Air-piloted Valves

CS 7750 FULL FUNCTION TERMINAL 15" OR 19" DISPLAY - THROUGH THE WALL 330mm (13")

DIEBOLD 5550 CASH DISPENSER THROUGH 330mm (13") MAXIMUM WALL PLAN VIEW FRONT ELEVATION SIDE ELEVATION CONSUMER ACCESS DIMENSIONS

Specialized Lending Rating Model using AHP

Comparison of retarder solvents. used in high-solids automotive refinish clear coatings

WATERTIGHT AIRTIGHT SPECIAL ACCESS DOORS & HATCHES, PART 1 GENERAL

330mm (13") MAXIMUM WALL THICKNESS IN AREA OF UNIT. 330mm (13") MAXIMUM WALL THICKNESS IN AREA OF UNIT (2) EACH SIDE OF SAFE CEN I SAFE ONLY

203mm (8") MAXIMUM WALL THICKNESS IN AREA OF UNIT NOTE: 203mm (8") MAXIMUM WALL THICKNESS IN AREA OF UNIT (2) EACH SIDE OF SAFE CEN I SAFE ONLY

CCS-810 KIT SERIES INSTRUCTION SHEET

Some thoughts on requirements for languages in engineering Requirements for Languages for modelling big systems World Ontology Summit,

4/3, 4/2 and 3/2 explosion-proof solenoid directional valve

World s Leading Supplier of Fastener Inspection Equipment

Accelerated detection of brown-rot decay: Comparison of soil block test, chemical analysis, mechanical properties, and immunodetection

Fourier Transform Infrared Quantitative Analysis of Sugars and Lignin in Pretreated Softwood Solid Residues

Effect of RME Fuel on Combustion in Optical Engine Operating in HCCI Mode

ORIENTAL MOTOR GENERAL CATALOG. Linear Heads A-227. Features A-228 General Specifications A-232 Product Specifications A-233. Standard AC Motors

Energy Management at University Farms

Roller and ball bars for easy and safe die change

MANUAL ON HYDROCARBON ANALYSIS

Case Study NEW LOGIC RESEARCH. SEP Membrane Filtration of Waste Oil. Overview. Used Oil Management System Site Map. Background

!!!!!!!!!!!!!!!!!!!!!!!!!!

DALI electronic ballasts

LTV 29-2 and 39-2 angle nutrunners. Higher speed, lower energy consumption

* * * DIEBOLD 3700 LOBBY CASH RECYCLER- FRONT LOAD 15" DISPLAY OR 19" DISPLAY CONSUMER ACCESS DIMENSIONS PLAN VIEW 19" DISPLAY FRONT ELEVATION

CCD-493 KIT SERIES INSTRUCTION SHEET

Special Valves. Connection size G1/8 to G2. Contents Version Actuation Port size Page Characteristics. Instructions Overview. Electrical G1/8 G1/4

Pneumatic Actuators 1000, 1400, 2800, and 2 x 2800 cm² Type 3271 and Type 3277 for integral positioner attachment Hand-operated Actuator Type 3273

STANDARD AND FULL HEIGHT - REAR LOAD - WITH 40mm (1 16") CEN I, CEN III GAS or CEN IV GAS SAFES

CC-410 KIT SERIES INSTRUCTION SHEET

Product range MEGAlife

Capabilities of airborne infrared remote sensing systems to detect hotspots

DIEBOLD 3750 CASH RECYCLER TERMINAL. 13mm ( 2") UL SAFE - 15" OR 19" DISPLAY THROUGH THE WALL 330mm (13") PLAN VIEW FRONT ELEVATION SECTION

Small Scale Washing Test System with a Magnetic Stirrer as a Model of Pulsator Type Washing Machine

CS 7780 FULL FUNCTION DRIVE-UP TERMINAL

Current methods for enumeration of total aerobic microorganisms

Carl Nielsen Music is Life. SET UP MANUAL Travelling Exhibition from Odense City Museums DK-5000 Denmark

MANN+HUMMEL Oil-bath air cleaner Single-stage air cleaner without spare parts

Metering pumps LK series

Estimation of Biodiesel Properties from Chemical Composition An Artificial Neural Network (ANN) Approach

PIPELINE PLAIN DENT FATIGUE A COMPARISON OF ASSESSMENT METHODOLOGIES

Exhaust System. Table of Contents. Section 6A - Exhaust System

JUMO MIDAS S05 OEM pressure transmitter - universal

Data of Traffic Accidents of Large Trucks in Japan

DEVELOPMENT OF A FAN AIRFLOW STATION FOR AIRFLOW CONTROL IN VAV SYSTEMS

Cooling Efficiency 2017/2018 Cooling Equipment Rebate Application

Selection guide. High performance fluids for precise temperature control

Redacted for Privacy DIT. JAMES MKNUDSIN

Hot Air Thermostats. WTHc series. Special features. Brief description. Switching function. Approval/approval marks (see "Technical data")


Transcription:

Environ. Sci. Technol. 2003, 37, 2869-2877 Chrcteriztion of Non-methne Hydrocrbons Emitted from Vrious Cookstoves Used in Chin STELLA MANCHUN TSAI, JUNFENG (JIM) ZHANG,*,, KIRK R. SMITH, YUQING MA, R. A. RASMUSSEN, AND M. A. K. KHALIL # University of Medicine nd Dentistry of New Jersey - School of Public Helth, Pisctwy, New Jersey, Environmentl nd Occuptionl Helth Sciences Institute, Pisctwy, New Jersey, Environmentl Helth Sciences, University of Cliforni, Berkeley, Cliforni, Institute for Techno-economics nd Energy System Anlysis, Tsinghu University, Beijing, Chin, Deprtment of Environmentl Sciences nd Engineering, Oregon Grdute Institute, Beverton, Oregon, nd Deprtment of Physics, Portlnd Stte University, Portlnd, Oregon Emission contributions from cookstoves to indoor, regionl, nd globl ir pollution lrgely depend on stove nd fuel types. This pper presents dtbse on emission fctors of specited non-methne hydrocrbons (NMHCs) for 16 fuel/stove combintions burning 2 types of crop residue, wood, 4 types of col, kerosene, nd 3 types of gseous fuels. The emission fctors re presented both on fuel mss bsis (compound mss per fuel mss) nd on cooking tsk bsis (compound mss per unit energy delivered to the pot). These fuel/stove combintions cover lrge spectrum of the cookstoves used in both urbn nd rurl households in Chin. Up to 54 hydrocrbons were identified, some of which re rective precursors of photochemicl smog. Bsed on published mximum incrementl rectivity (MIR) vlues for NMHCs, we estimted stove-specific nd fuel-specific ozone forming potentils (OFPs). The results indicte tht rw col powder, wood, nd crop residues hve higher OFP vlues thn the other types of fuels tested. Strikingly, burning the col briquette nd honeycomb col briquette produced OFP vlues more thn 2 orders of mgnitude lower thn burning unprocessed (rw) col, even in the sme vented metl stove, for every 1 MJ delivered to the pot. Introduction In mny prts of developing countries, household cookstoves re the mjor sources of ir pollution exposure (1-3). Elsewhere we hve presented results from studies of cookstove emissions, including pilot study in the Philippines * Corresponding uthor phone: (732)445-0158; fx: (732)445-0116; e-mil: jjzhng@eohsi.rutgers.edu. Corresponding uthor ddress: EOHSI, Room 358, 170 Frelinghuysen Rod, Pisctwy, NJ 08854. University of Medicine nd Dentistry of New Jersey - School of Public Helth. Environmentl nd Occuptionl Helth Sciences Institute. University of Cliforni. Tsinghu University. Oregon Grdute Institute. # Portlnd Stte University. (4), nd two reporting the results of more extensive studies in which emissions from 28 fuel/stove combintions were mesured in Indi nd in Chin (5, 6). In both the pilot nd the more extensive studies, we mesured stove performnce nd emissions of severl importnt greenhouse gses (GHGs) nd other pollutnts of helth concern, including crbon dioxide (CO 2), crbon monoxide (CO), methne (CH 4), totl non-methne hydrocrbon (TNMHC), totl suspended prticles (TSP), oxides of nitrogen (NO x) (for Chinese stoves only), nd sulfur dioxide (for Chinese stoves only). These studies found tht gret frction of fuel crbon ws converted to products of incomplete combustion (PICs), rther thn CO 2, in simple household stoves nd tht burning solid fuels generted substntilly greter emissions of totl PICs thn burning liquid nd gseous fuels to deliver the sme mount of energy to the pot. The studies lso found tht the emission of non-co 2 GHGs (e.g., CH 4, nitrous oxide, nd CO) cn be significnt contributor to totl globl wrming potentils for mny mesured cookstoves (7). In the forementioned extensive studies of cookstove emissions, subset of the fuel/stove combintions were lso mesured for specited NMHCs, the results of which hve not yet been reported. In this pper, we present dtbse of emission fctors of NMHCs for the 16 fuel/stove combintions mesured in Chin. The dtbse for individul NMHCs cn provide better informtion for future study on helth risk ssessment bsed on prticulr toxic compounds emitted from household cookstoves. Respirtory illnesses, including cute respirtory infection, chronic obstructive lung diseses, nd lung cncer, hve been strongly ssocited with exposure to indoor smoke from solid-fuel stoves (8). Other effects, including tuberculosis, dverse pregnncy outcomes, ero-digestive cncer, ctrcts, nd sthm, hve lso been ssocited with such exposures, but the evidence bse is still smll (8-10). Synergistic nd ntgonistic effects of toxic compounds present in the cookstove smoke my contribute to the observed dverse heth effects. Knowing detiled constituents of the smoke will certinly help in understnding its toxicologicl properties. In ddition to some NMHCs tht hve direct helth concerns, some of the compounds identified re rective precursors, in the tmosphere, leding to the formtion of ozone (O 3) nd photochemicl smog. The dtbse on individul NMHCs, thus, cn be useful in modeling locl or regionl O 3. In this pper, we compre stove- nd fuel-specific ozone forming potentils (OFPs) using the mesured emission fctors nd the published rectivity indices for NMHCs. The nlysis of stove/fuel OFPs my provide useful energy policy implictions in Chin from stnd point of O 3 control strtegies. Methods Fuel/Stove Combintion Mesured. The fuels mesured included biomss fuels (crop residues nd wood), severl types of cols, kerosene, liquefied petroleum gs (LPG), col gs, nd nturl gs. These selected fuels covered fuels of common use in rurl nd urbn households in Chin in 1995-1996 when the mesurements were mde. The stove selection ws bsed on populr models for ech fuel type. Since 1996, usge of gs fuels my hve been incresed in some mjor Chinese cities (e.g., Beijing) (11, 12). However, no significnt chnges in household fuel types hve occurred in most res of Chin. The stove/fuel combintions tested in this study re still commonly found in Chin s whole, lthough some stove/fuel combintions found in 1995-1996 hve been replced with others in some cities or res. Detiled 10.1021/es026232 CCC: $25.00 2003 Americn Chemicl Society VOL. 37, NO. 13, 2003 / ENVIRONMENTAL SCIENCE & TECHNOLOGY 9 2869 Published on Web 06/04/2003

TABLE 1. Description of the Fuel/Stove Combintions Tested in Chin description seril no. symbol fuel/stove fuel stove 1 Honey- honeycomb col briquette metl col stove with flue 2 Honey-Metl honeycomb col briquette metl col stove without flue 3 Honey-Imp. honeycomb col briquette improved metl col stove without flue 4 Briq-Metl col briquette metl col stove without flue 5 Wsh wshed col powder metl col stove with flue 6 unprocessed col (col powder) metl col stove with flue 7 unprocessed col (col powder) brick stove with flue 8 Wood- fuel wood brick stove with flue 9 Wood-Imp.-v fuel wood improved brick stove with flue 10 Whet- whet residue brick stove with flue 11 Mize- mize residue brick stove with flue 12 Mize-Imp.-v mize residue improved brick stove with flue 13 Kero-Wick kerosene kerosene wick stove without flue 14 LPG- liquefied petroleum gs LPG trditionl stove without flue 15 col gs trditionl gs stove without flue 16 Nturl nturl gs trditionl gs stove without flue Flue code: v ) vented, i.e., with flue. The fuel ctegory includes: #1-7 for col; #8-9 for wood; #10-12 for crop residue; #13 for kerosene; nd #14-16 for gseous fuel. powder is rw col with lrge pieces being sieved out. power consists of col prticles (pieces) with dimeters typiclly smller thn 1 cm. powder is the most common unprocessed col used in Chinese households. informtion bout the fuel sources, fuel compositions, nd stove chrcteristics cn be found in previous published pper (6). Briefly, for the stoves using piped gseous fuels, the mesurements were conducted in ctul homes. The emissions of ll the other fuel/stove combintions were mesured in simulted villge kitchen locted in the rurl cmpus of Tsinghu University, Beijing, Chin. Among the 28 fuel/stove combintions mesured primrily for mjor GHGs, 16 were lso successfully mesured for specited NMHCs. These 16 fuel/stove combintions, described in Tble 1, represent the most commonly used ones nd covers t lest one fuel/stove combintion within ech brod fuel type (col, wood, crop residue, kerosene, nd gs). With the exception of Honey- combintion, every fuel/stove combintion ws mesured for NMHCs only once for complete burn cycle, due to the budgetry constrint. However, we were ble to mke three repeted mesurements for the Honey- to get sense bout the vrition from one burn cycle to nother burn cycle. This ws evluted by reporting the coefficient of vrition (CV) determined for this prticulr fuel/stove combintion. Uncertinties ssocited with emission fctors for the other tested fuel/stove combintions were discussed bsed on the CV vlues of TNMHC reported in the previous pper nd n nlysis of uncertinty sources (6). Mesurement Methods. Additionl detils on experimentl design nd smple collection procedures cn be found in Zhng et l. (6). Briefly, flue gs smples were collected using smpling configurtion tht included, from upstrem to downstrem of the smpling trin, stinless steel probe, filter holder, smpling pump (SKC Inc., U.S.A.), nd clen Tedlr bg (SKC Inc., U.S.A.). The filter holder contined 37 mm qurtz fiber filter tht ws used to collect TSP. For stoves with flue, the probe ws inserted in the flue. For those hving no flues, the stoves were plced under hood built for the test purpose nd the probe ws plced inside the hood exhust duct. The flue gs smples hd gone through dilution nd cooling before they were collected in the Tedlr bgs. Thus, the temperture of the flue gs smples ws the sme s the mbient temperture. The flow rte of the smpling pump, rnging from 0.2 to 2 L/min, ws djusted to fill one or two 80-L Tedlr bgs throughout whole burn cycle (i.e., from fire strt to fire extinction). If two Tedlr bgs were used, time-weighted frction of ir from ech of the two bgs ws tken nd then mixed in the third Tedlr bg for finl smple nlysis. Hence, the flue smples collected were time-verged emissions covering complete burn cycle consisting of different combustion stges. Two mbient ir smples were collected nerby the simulted kitchen for bckground corrections. Aliquots of ir smples were tken out of the Tedlr bgs for vrious subsequent chemicl nlyses. These included the determintion of flue gs concentrtion of CO 2, CO, CH 4, nd TNMHC, nd crbon in irborne prticultes, ll necessry for the construction of crbon blnce model, long with mesured crbon mss in the fuel combusted nd in combustion residues (for solid fuel only). The crbon blnce model ws used to determine emission fctors of ech mesured species in the flue gs (6). As described in detil by Zhng et l. (6), the crbon blnce pproch only requires the mesurement of emission rtio, i.e., [species concentrtion]/[co 2 concentrtion] of n emitted irborne species. The emission fctor of given species in the flue gs cn be clculted from the CO 2 emission fctor nd the emission rtio of the species. Therefore, here we report emission rtios of ll specited NMHCs, which were used to derive emission fctors. Emission fctors were reported, by convention, on dry fuel mss bsis (i.e., on moisture free bsis). The moisture content of ech fuel ws mesured through stndrd fuel moisture nlysis, s described in detil in the previous pper (6). Since different mounts of fuels re needed for the sme cooking tsk using different fuel/stove combintions, tskbsed emission fctors (pollutnt mss per cooking tsk) rther thn the fuel mss bsed is better performnce index to compre the ir pollution potentil of different fuel/stove combintions (2, 13). For this reson, the emission fctors were lso reported s mss of compound emitted per unit energy, one meg-joule (MJ), delivered to the pot. We chose 1 MJ s unit energy becuse bout 5 MJ of delivered het is roughly wht typicl fmily mel would require, lthough this obviously vry lrgely by mel size, type, nd cooking method (5). The conversion between the mss-bsed nd tsk-bsed emission fctors ws chieved using the fuel energy content (lso known s clorific vlue) (MJ/kg) nd stove therml efficiency (%), both of which were mesured in the study. Therml efficiencies [energy effectively received by the pot]/[totl energy in the fuel burned] of the stoves were determined using wter boiling tests, s described in detil in the previous pper (6). The mesurement of flue gs concentrtions of specited NMHCs involved filling clen evcuted 850-mL stinless 2870 9 ENVIRONMENTAL SCIENCE & TECHNOLOGY / VOL. 37, NO. 13, 2003

steel cnister with the flue gs collected in the lrge Tedlr bgs to 1.5-2 tm pressure using bttery-powered pump. The cnisters were preclened nd provided by the lbortory t Oregon Grdute Institute (OGI), Beverton, OR. The cnisters contining flue gs smples were shipped bck by ir crriers to the OGI lbortory for nlysis. Storge time before nlysis vried from one btch to nother but no more thn 30 dys. There were no field blnk or duplicte smples collected due to budget constrint. The bckground concentrtion ws determined from two mbient smples collected nerby the simulted kitchen nd ws used to correct for net flue gs concentrtions. The hydrocrbon specition nlysis ws mde using the procedure estblished s EPA Compendium Method TO-14A (14). A system of gs chromtogrphy (GC) equipped with flme ioniztion detector (FID) nd mss-selective detector ws used for the nlysis. In ech nlyticl run, smll liquot of ir smple ws collected from the cnister into the nlyticl system through dryer (to remove moisture), chromtogrphic vlve, nd then cryogenic trp. The trp ws heted rpidly, nd the smple ws injected onto n OV-1 cpillry column. High purity helium ws used s the crrier gs. A stndrd mixture contining 74 NMHCs (Scott Specilty Gses, Inc., Sn Bernrdino, CA) ws used to generte clibrtion curves for every btch of smples. All clibrtion curves hd liner regression R 2 vlues > 0.99. The detection limit vlues vried by individul compound but were in sub-ppb rnges. Clcultion of Ozone Forming Potentil. Among the NMHCs identified in the cookstove flue gs, totl of 44 compounds were selected for ozone forming potentil clcultion. These compounds re on trget list, which ws developed by the U.S. EPA Photochemicl Assessment Monitoring Sttions (PAMS) progrm, due to their rectivity in terms of producing tropospheric ozone nd photochemicl smog. We estimted ozone forming potentils (OFPs) of the 16 fuel/stove combintions using the mesured emission fctors nd vlues of the Mximum Incrementl Rectivity (MIR) scle developed by Crter (15) for most of the NMHCs mesured, except o-xylene nd sec-butylbenzne. Since only totl xylenes (o+m+p) were mesured, the published MIR vlue for m,p-xylenes ws pplied to totl xylenes. Similrly, since 1,2,4-trimethylbenzene nd sec-butylbenzene were not seprted for the flue gs smples, we pplied the MIR vlue for 1,2,4-trimethylbenzene to the sum of 1,2,4-trimethylbenzene nd sec-butylbenzene. The MIR of individul NMHC ws defined in terms of grms of ozone formed by photochemicl rections in the tmosphere per grm of crbon (C) in NMHC emitted (g ozone/g C emitted). The totl ozone forming potentil (OFP) for ech fuel/stove combintion is the sum of OFPs of individul NMHCs. The OFP from ech individul NMHC ws clculted s follows. Per kg fuel mss OFP: OFP (g ozone/kilogrm of dry fuel) ) [MIR of given NMHC (g ozone/g C)] [emission fctor of the NMHC (g C/ kilogrm of dry fuel)] Per tsk (1 MJ energy delivered) OFP: OFP (g ozone/mj energy delivered) ) [MIR of given NMHC (g ozone/ g C)] [emission fctor of individul NMHC (g C/MJ energy delivered)] Results Molr Emission Rtios. Non-CO 2 crbon-contining compounds detected in the flue gs re typiclly clled products of incomplete combustion (PICs). The rtio of n individul species to CO 2 in the flue gs, nmely emission rtio, represents reltive bundnce of the species in the flue gs. Tble 2 presents molr emission rtios of 54 specited NMHCs, mong the 74 trgeted species, tht were detected for ny of the 16 fuel/stove combintions tested. The emission rtios were clculted using net concentrtions of individul NMHC nd net concentrtions of CO 2. (The net concentrtion equls flue gs concentrtion minus bckground concentrtion). The vlues for totl specited NMHC were the sums of the molr emission rtios of ll the quntified NMHCs listed in Tble 2. The rnking of reltive bundnce of totl specited NMHC to CO 2 is lso shown in Tble 2. Burning biomss fuels (whet residue, mize residue, nd wood) generlly hd greter emission rtios thn burning gseous fuels. It is striking tht col type hd profound effect on the emission rtio of the sum of ll the quntified NMHCs. Among the 16 tested fuel/stove combintions, the stoves burning unprocessed col powder nd wshed col powder hs highest emission rtios while the vented honeycomb col stove hd the lowest emission rtio. The results from the three repeted mesurements of Honey- showed lrge rnge of CV (8-113%) for the emission fctors of the NMHCs detected in the flue gs, providing rough ide bout run-to-run vritions in NMHC emissions. By nture, the emission cn vry lrgely from one burn cycle to nother cycle. Therefore the uncertinty ssocited with the single mesurement vlues reported in this pper is expected to be lrge, which is discussed lter in the Discussion Section. Emission Fctors. Tble 3 presents emission fctors of ech identified compound, by fuel/stove combintion, on dry-fuel mss bsis. The emission fctors in Tble 3 represent the mount of individul compound in grms (g) generted by burning 1 kg of dry fuel with ech fuel/stove combintion. Tht is, for exmple, burning 1 kg of honeycomb col briquette using metl col stove with flue (Honey-) produced 1.98 mg of benzene. The emissions of totl specited NMHC from burning 1 kg of fuel rnged from 11.9 mg to 6.48 g cross the fuel/stove combintions mesured. Even within fuel type, the rnge of totl specited NMHC emission fctor cn be lrge for different fuel species nd stove types. The col/stove combintions hd rnge from 11.9 mg to 6.48 g; the wood stoves hd rnge from 1.36 nd 2.21 g; the crop residues hd rnge from 0.73 to 2.28 g; nd the stoves burning gseous fuels hd rnge from 0.024 to 0.42 g. Tble 4 presents emission fctors on delivered energy bsis. The dt on stove efficiency nd fuel energy content (mesured s low-het clorific vlue) re lso included in Tble 4. These vlues were used to convert between fuelmss-bsed nd tsk-bsed emission fctors. Emission fctors in Tble 4 represent the mount of ech compound generted from delivering 1 meg-joule (MJ) het to the pot by given fuel/stove combintion. These tsk-bsed emission fctors of totl specited NMHCs lso hd lrge rnge, from 1.11 mg to 2.41 g per 1 MJ het delivered. Within-fuel type rnges were 3.28 mg to 2.41 g for cols, 0.57 g to 0.82 g for wood, 0.37 g to 1.52 g for crop residues, 0.02 g for kerosene, nd 1.11 mg to 17.5 mg for gseous fuels. Ozone Forming Potentils. The estimted ozone forming potentils (OFPs) on dry-fuel mss bsis rnged from 9 mg to 11 g ozone per kilogrm of dry fuel from 16 fuel-stove combintions tested in this study. The estimted men OFP ws 6.03 g ozone per kilogrm of dry fuel for col (col powder), 0.03 g for col briquettes nd honeycomb col briquettes, 1.68 g for wood, 1.75 g for crop residue, 0.74 g for kerosene, 0.14 g for LPG, 0.03 g for col gs, nd 0.06 g for nturl gs. The estimted OFPs on delivered energy bsis rnged from 1.41 mg to 4.05 g per 1 MJ of het delivered to the pot. Figure 1 shows delivering the sme mount of het to the VOL. 37, NO. 13, 2003 / ENVIRONMENTAL SCIENCE & TECHNOLOGY 9 2871

TABLE 2. Molr Emission Rtios to CO 2 (Unit in 10-7 ) Honey- CV fuel/stove b men c (%) Honey- Metl Honey- Imp. Briq- Metl Wsh Wood- Wood- Imp.-v Whet- Mize- Mize- Imp.-v benzene 5.91 42 26.0 29.3 25.8 989 2437 78.8 977 2225 2122 495 949 81.2 147 12.8 81.9 1,3-butdiene nd n nd 0.061 nd 43.0 71.5 nd 5.16 4.31 26.3 nd nd 8.16 0.458 nd nd xylenes (o+m+p) 0.488 60 3.31 nd 8.74 4.28 14.0 7.59 1.53 nd 3.36 nd 0.578 nd 73.3 3.99 21.2 styrene nd n nd nd nd nd nd nd nd nd nd nd nd nd 243 nd 0.850 ethne 4.87 30 34.3 2.92 63.6 6557 4193 94.4 177 137.5 372 513 371 10.5 3.48 31.0 33.8 ethylene 9.30 16 44.6 11.9 86.0 8708 18474 467 2652 6038 8068 3986 4539 933 4.52 58.5 26.1 cetylene 9.78 8 34.6 3.88 17.8 2678 5725 352 5591 9237 10490 3447 7846 529 9.10 5.16 5.83 propne 0.689 14 6.18 0.597 14.7 1614 639 29.0 20.1 34.3 55.4 44.0 113 4.09 6.37 5.72 12.4 propene 1.16 12 14.4 5.27 24.5 2098 1930 73.3 154 205.5 704 273 517 117 18.3 8.10 3.02 i-butne 0.106 n 1.88 nd 1.66 94.1 31.1 3.60 nd nd 1.98 2.93 4.49 3.41 5.97 0.558 2.89 i-butene 0.156 49 2.91 1.06 2.66 124 nd 2.58 5.82 nd 27.2 7.52 26.3 8.99 3.54 0.899 0.316 1-butene 0.073 n 1.21 0.352 1.89 352 261 6.64 13.7 20.8 67.4 31.6 80.9 25.6 5.89 0.514 0.158 n-butne 0.142 80 1.95 nd 3.95 368 134 7.59 3.81 3.55 5.18 5.87 14.0 0.537 1.33 1.22 3.48 trns-2-butene nd n 1.03 1.00 1.50 55.8 83.1 1.95 43.3 46.4 91.6 23.8 55.2 4.72 2.80 0.289 2.17 2,2-dimethylpropne nd n nd nd nd nd nd nd nd nd nd nd nd nd 3.19 nd nd cis-2-butene nd n 0.486 0.235 1.37 10.8 nd 7.42 2.05 nd 8.97 3.09 10.4 2.12 8.25 3.34 12.8 3-methyl-1-butene nd n nd nd 0.206 22.3 5.82 0.623 nd nd 4.94 1.92 3.49 1.70 0.236 0.205 nd i-pentne 0.064 n 2.52 nd 1.98 65.7 19.2 0.914 nd nd 0.620 2.22 0.176 nd nd 0.283 0.716 1-pentene nd n 0.486 0.188 0.481 2.72 nd 2.07 1.83 nd 8.37 nd 7.59 8.94 nd 0.205 2.21 2-methyl-1-butene nd n 0.486 nd 0.275 20.3 3.64 nd nd nd 4.17 nd 1.54 2.00 0.236 0.051 nd n-pentne 0.101 38 1.13 0.183 1.74 167 60.8 2.20 nd nd 1.64 2.93 0.249 nd nd 0.300 0.982 trns-2-pentene nd n nd nd 0.481 2.56 nd 0.131 nd nd 4.15 nd 3.43 1.03 nd nd nd cis-2-pentene nd n nd nd 0.206 1.78 nd nd nd nd 2.71 nd 2.21 0.616 0.236 nd nd 2-methyl-2-butene nd n nd nd 0.618 2.27 nd nd nd nd 0.224 nd 7.34 0.822 0.943 nd nd cyclopentene nd n nd nd nd nd nd nd nd nd 0.767 nd nd 1.48 nd nd nd 4-methyl-1-pentene nd n nd nd nd nd nd nd nd nd 0.828 nd nd nd nd nd nd cyclopentne nd n 0.049 nd 0.034 10.9 nd 0.213 nd nd nd nd nd nd nd 0.077 nd 2,3-dimethylbutne nd n 0.198 nd 0.140 15.6 6.75 nd nd nd nd nd nd nd nd nd nd 2-methylpentne 0.075 21 1.27 nd 0.615 36.9 10.3 0.347 nd nd 0.222 nd nd nd nd 0.125 0.514 3-methylpentne 0.066 88 0.633 nd 0.447 9.98 3.01 0.160 nd nd nd nd nd nd nd nd nd 1-hexene nd n 0.324 nd 0.573 56.6 24.0 2.13 nd nd 7.22 nd 5.85 7.36 0.786 0.428 nd n-hexne 0.151 38 1.11 0.153 1.23 96.1 34.9 1.12 nd nd 0.445 0.921 1.10 0.084 nd 0.125 0.462 methylcyclopentne nd n 0.594 nd 0.191 27.4 9.18 nd nd nd nd nd nd nd nd nd 0.070 2,4-dimethylpentne nd n 0.680 nd 0.144 nd nd nd nd nd nd nd nd nd nd nd nd cyclohexne nd n nd nd nd 0.842 2.95 nd nd nd nd nd nd nd 13.6 nd nd 2-methylhexne 0.052 45 nd nd 0.265 6.65 1.78 nd nd nd nd nd nd nd nd nd 0.066 2,3-dimethylpentne nd n nd nd 0.024 9.87 1.82 nd nd nd nd nd nd nd nd nd 0.243 3-methylhexne 0.113 n nd nd 0.417 6.20 2.03 0.865 nd nd 0.064 nd nd nd nd 0.060 0.780 2,2,4-trimethylpentne 0.671 n nd nd nd 10.1 5.58 nd nd nd 1.45 nd nd 3.82 nd nd nd n-heptne 0.169 87 0.793 0.307 1.14 56.1 20.3 1.00 nd nd 0.203 nd 0.710 0.851 nd nd 0.250 methylcyclohexne 0.094 n 0.879 nd 0.573 11.0 5.75 0.016 nd nd nd nd nd 2.04 nd nd nd 2,3,4-trimethylpentne nd n nd nd nd nd nd nd nd nd nd nd nd 0.252 nd nd nd toluene 1.17 57 6.75 4.24 10.9 129 316 18.5 58.5 86.6 155 10.7 39.4 12.4 51.0 4.90 24.0 2-methylheptne nd n 0.030 nd 0.274 6.06 nd nd nd nd 0.084 nd nd 0.552 nd nd nd 3-ethylhexne nd n 0.239 nd 0.253 1.06 nd nd nd nd nd nd nd 0.978 nd nd nd n-octne 0.157 45 1.11 0.557 0.914 14.9 5.95 0.577 nd nd 0.010 1.36 nd 1.50 nd 0.053 0.258 ethylbenzene 0.218 75 0.855 0.145 1.15 nd 2.25 1.59 1.32 nd 0.891 nd 0.515 nd 5.03 0.317 1.89 n-nonne 0.396 64 1.88 2.22 2.08 1.38 1.41 0.078 nd nd nd nd nd nd 4.43 nd nd i-propylbenzene nd n nd nd 0.581 nd nd nd nd nd nd nd nd nd nd 0.015 nd n-propylbenzene nd n nd nd nd nd nd nd nd nd nd nd nd nd nd 0.659 nd p-ethyltoluene nd n nd nd nd nd nd nd nd nd nd nd nd nd nd 0.389 nd m-ethyltoluene nd n nd nd nd nd nd nd nd nd nd nd nd nd nd 0.060 nd 1,2,4-trimethylbenzene nd sec-butylbenzene 2.43 n nd 2.57 nd nd nd 2.93 nd nd nd nd nd nd nd 0.060 nd n-decne 0.599 113 2.27 4.37 3.40 nd nd nd 529 nd nd nd nd nd 24.7 0.160 nd totl specited NMHC 36.3 21 197 71.5 285 24488 34536 1167 10236 18039 22237 8853 14601 1775 638 141 239 rnking (1 ) highest) 16 13 15 11 2 1 9 6 4 3 7 5 8 10 14 12 Key: nd ) not detected or bckground level g flue gs concentrtion; CV (coefficient of vrition) ) (stndrd devition)/men; n ) vrince clcultion cnnot be mde, only one or no mesurement with detected level. b See Tble 1 for fuel/stove codes. c Averge of 3 burning cycles using Honey- combintion. Kero- Wick LPG- Nturl pot by different col/stove combintions would contribute quite differently to the ozone forming potentils. The results indicted tht burning col powder nd biomss fuels would contribute highest ozone forming potentil mong the 16 tested fuel-stove combintions on delivered energy bsis. Strikingly, burning the col briquette nd honeycomb col briquette produced OFP vlues more thn 2 orders of mgnitude lower thn burning unprocessed col, even in the sme vented metl stove. However, the interprettion of these OFP results is very crude becuse the MIR vlues used in our estimtion were from the dtbse developed for the United Sttes photochemicl sitution which my be different from the Chinese sitution due to the possible difference in reltive mkeup of ozone precursors (NO x nd NMHCs). To put these cookstove OFPs into perspective, we compre them with the published vehiculr OFPs reported by Chng et l. for gsoline-powered nd LPG-powered vehicles representtive of the Tipei fleet (16). The OFPs were estimted bsed on emission fctors of some 56 NMHCs nd their MIR vlues: 1.8 g/km-vehicle for (verge) gsoline-powered vehicles nd 0.97 g/km-vehicle for (verge) LPG-powered vehicles. If we use 5 MJ delivered energy s wht is roughly needed for cooking mel (5), cooking one mel would produce OFPs rnging from 0.007 g (the stove) 2872 9 ENVIRONMENTAL SCIENCE & TECHNOLOGY / VOL. 37, NO. 13, 2003

TABLE 3. Emission Fctors (mg of Compound per kg of Dry Fuel, by Fuel/Stove Combintion) b fuel/stove c Honey Honey- Honey- d Metl Imp. Briq- Metl Wsh Wood- Wood- Imp.-v Whet- Mize- Mize- Imp.-v Kero- Wick LPG- Nturl benzene 2.71 11.0 14.0 7.40 440 1050 25.8 264 629 512 102 194 44.9 80.4 4.19 50.0 1,3-butdiene nd nd 0.020 nd 13.3 21.3 nd 0.963 0.843 4.40 nd nd 3.12 0.173 nd nd xylenes (o+m+p) 0.305 1.91 nd 3.41 2.58 8.22 3.38 0.560 nd 1.10 nd 0.161 nd 54.3 1.77 17.6 styrene nd nd nd nd nd nd nd nd nd nd nd nd nd 177 nd 0.691 ethne 0.857 5.58 0.538 7.02 1121 694 11.9 18.3 14.9 34.6 40.8 29.3 2.23 0.731 3.90 7.92 ethylene 1.52 6.77 2.04 8.87 1390 2856 54.9 257 612 699 296 334 185 0.885 6.87 5.72 cetylene 1.47 4.88 0.618 1.71 397 822 38.4 503 870 844 237 536 97.6 1.65 0.563 1.19 propne 0.177 1.48 0.161 2.37 405 155 5.36 3.06 5.47 7.54 5.12 13.1 1.28 1.96 1.06 4.28 propene 0.283 3.27 1.36 3.79 502 448 12.9 22.3 31.3 91.4 30.3 57.1 34.7 5.39 1.43 0.991 i-butne 0.034 0.590 nd 0.354 31.1 9.97 0.877 nd nd 0.355 0.450 0.683 1.40 2.42 0.136 1.31 i-butene 0.050 0.885 0.363 0.549 39.7 nd 0.607 1.13 nd 4.71 1.12 3.86 3.57 1.39 0.211 0.138 1-butene 0.023 0.369 0.121 0.390 113 80.6 1.56 2.65 4.22 11.7 4.68 11.9 10.2 2.31 0.121 0.069 n-butne 0.046 0.615 nd 0.844 122 42.8 1.85 0.763 0.747 0.929 0.902 2.14 0.221 0.539 0.297 1.58 trns-2-butene nd 0.313 0.343 0.310 17.8 25.7 0.458 8.39 9.41 15.9 3.53 8.12 1.87 1.10 0.068 0.950 2,2-dimethylpropne nd nd nd nd nd nd nd nd nd nd nd nd nd 1.62 nd nd cis-2-butene nd 0.148 0.081 0.283 3.44 nd 1.74 0.397 nd 1.55 0.458 1.53 0.842 3.23 0.784 5.60 3-methyl-1-butene nd nd nd 0.053 8.90 2.25 0.183 nd nd 1.07 0.356 0.642 0.842 0.115 0.060 nd i-pentne 0.027 0.984 nd 0.526 27.0 7.64 0.276 nd nd 0.138 0.424 0.033 nd nd 0.085 0.403 1-pentene nd 0.184 0.081 0.124 1.09 nd 0.607 0.443 nd 1.81 nd 1.40 4.44 nd 0.060 1.21 2-methyl-1-butene nd 0.184 nd 0.071 8.12 1.41 nd nd nd 0.904 nd 0.283 0.995 0.115 0.015 nd n-pentne 0.042 0.443 0.081 0.461 68.5 24.2 0.665 nd nd 0.366 0.558 0.047 nd nd 0.090 0.553 trns-2-pentene nd nd nd 0.124 1.02 nd 0.039 nd nd 0.898 nd 0.630 0.510 nd nd nd cis-2-pentene nd nd nd 0.053 0.709 nd nd nd nd 0.586 nd 0.407 0.306 0.115 nd nd 2-methyl-2-butene nd nd nd 0.159 0.904 nd nd nd nd 0.048 nd 1.35 0.408 0.462 nd nd cyclopentene nd nd nd nd nd nd nd nd nd 0.161 nd nd 0.714 nd nd nd 4-methyl-1-pentene nd nd nd nd nd nd nd nd nd 0.215 nd nd nd nd nd nd cyclopentne nd 0.018 nd 0.009 4.35 nd 0.063 nd nd nd nd nd nd nd 0.023 nd 2,3-dimethylbutne nd 0.092 nd 0.044 7.66 3.21 nd nd nd nd nd nd nd nd nd nd 2-methylpentne 0.038 0.590 nd 0.195 18.1 4.91 0.125 nd nd 0.059 nd nd nd nd 0.045 0.346 3-methylpentne 0.032 0.295 nd 0.142 4.89 1.43 0.058 nd nd nd nd nd nd nd nd nd 1-hexene nd 0.148 nd 0.177 27.1 11.1 0.752 nd nd 1.88 nd 1.29 4.39 0.462 0.151 nd n-hexne 0.076 0.516 0.081 0.390 47.1 16.6 0.405 nd nd 0.118 0.210 0.247 0.051 nd 0.045 0.311 methylcyclopentne nd 0.270 nd 0.059 13.1 4.26 nd nd nd nd nd nd nd nd nd 0.046 2,4-dimethylpentne nd 0.369 nd 0.053 nd nd nd nd nd nd nd nd nd nd nd nd cyclohexne nd nd nd nd 0.403 1.37 nd nd nd nd nd nd nd 7.96 nd nd 2-methylhexne 0.031 nd nd 0.097 3.79 0.983 nd nd nd nd nd nd nd nd nd 0.052 2,3-dimethylpentne nd nd nd 0.009 5.63 1.00 nd nd nd nd nd nd nd nd nd 0.190 3-methylhexne 0.069 nd nd 0.154 3.54 1.12 0.363 nd nd 0.020 nd nd nd nd 0.025 0.610 2,2,4-trimethylpentne 0.446 nd nd nd 6.54 3.51 nd nd nd 0.511 nd nd 3.09 nd nd nd n-heptne 0.102 0.430 0.188 0.419 32.0 11.2 0.421 nd nd 0.063 nd 0.187 0.604 nd nd 0.196 methylcyclohexne 0.056 0.467 nd 0.207 6.13 3.11 0.006 nd nd nd nd nd 1.42 nd nd nd 2,3,4-trimethylpentne nd nd nd nd nd nd nd nd nd nd nd nd 0.204 nd nd nd toluene 0.631 3.37 2.39 3.71 67.5 161 7.15 18.6 28.9 44.2 2.60 9.52 8.12 32.8 1.89 17.3 2-methylheptne nd 0.018 nd 0.115 3.94 nd nd nd nd 0.030 nd nd 0.446 nd nd nd 3-ethylhexne nd 0.148 nd 0.106 0.689 nd nd nd nd nd nd nd 0.791 nd nd nd n-octne 0.104 0.688 0.390 0.384 9.70 3.74 0.276 nd nd 0.004 0.411 nd 1.22 nd 0.025 0.230 ethylbenzene 0.137 0.492 0.094 0.449 nd 1.32 0.707 0.485 nd 0.292 nd 0.143 nd 3.73 0.141 1.57 n-nonne 0.304 1.30 1.75 0.980 1.00 0.995 0.042 nd nd nd nd nd nd 3.96 nd nd i-propylbenzene nd nd nd 0.257 nd nd nd nd nd nd nd nd nd nd 0.008 nd n-propylbenzene nd nd nd nd nd nd nd nd nd nd nd nd nd nd 0.332 nd p-ethyltoluene nd nd nd nd nd nd nd nd nd nd nd nd nd nd 0.196 nd m-ethyltoluene nd nd nd nd nd nd nd nd nd nd nd nd nd nd 0.030 nd 1,2,4-trimethylbenzene 1.77 nd 1.90 nd nd nd 1.47 nd nd nd nd nd nd nd 0.030 nd nd sec-butylbenzene n-decne 0.514 1.75 3.80 1.78 nd nd nd 260 nd nd nd nd nd 24.5 0.096 nd totl specited NMHC 11.9 50.6 30.4 48.6 4977 6480 173 1361 2207 2282 727 1208 416 409 24.7 121 Zhng et l. (6) presented detiled informtion on clcultion for emission fctors using the crbon blnce model. b Key: nd ) not detected or bckground level g flue gs concentrtion. c See Tble 1 for fuel/stove codes. d Averge of 3 burning cycles using Honey- combintion. to 20.3 g (the stove). At the lower bound of the cookstove OFP rnge, cooking 257 mels would produce OFP equivlent to driving gsoline-powered vehicle for 1 km, nd cooking 139 mels would produce OFP equivlent to driving LPG-powered vehicle for 1 km. At the upper bound of the cookstove OFP rnge, in contrst, cooking 1 mel would produce OFP equivlent to driving gsoline-powered vehicle for bout 11 km or equivlent to driving LPG-powered vehicle for bout 21 km. Discussion Uncertinties for Emission Fctors. Due to the budgetry constrint, only one fuel/stove combintion (Honey-Metlv), mong the 16 combintions mesured for NMHC emissions, ws repetedly mesured for three independent burn cycles (experimentl runs). These mesurements showed high run-to-run vritions (8-113% CV) in emission fctors estimted using the crbon blnce model for the Honey- combintion. High CV vlues were lso expected for the emission fctors of specited NMHCs for the other 15 fuel/stove combintions bsed on the nlysis of uncertinty sources of the crbon blnce pproch. There were two mjor sources of uncertinty for this pproch. One ws the vribility in emissions from one burn cycle to nother, lrgely ssocited with vritions in fire tending behvior nd vritions in fuel prmeters (e.g., fuel size, fuel lyout in the stove). This type of vribility ws obviously greter for the solid fuel stoves thn for the LPG nd gs stoves. The other source ws the mesurement errors for ll the prmeters nd concentrtion vlues used in the crbon VOL. 37, NO. 13, 2003 / ENVIRONMENTAL SCIENCE & TECHNOLOGY 9 2873

TABLE 4. Emission Fctors (mg Compound/MJ Delivered to the Pot, by Fuel/Stove Combintion) fuel/stove Honey Honey- Honey- b Metl Imp. Briq- Metl Wsh Wood- Wood- Whet- Mize- Mize- Imp.-v Imp.-v Kero- Wick LPG- Nturl stove efficiency (%) 16 22 48 37 10 10 26 10 24 11 12 19 47 48 51 54 low-het CV(MJ/kg fuel) 19.2 19.2 19.2 13.9 30.1 27.3 27.3 16.3 16.3 14.0 16.1 16.1 43.3 49.0 43.8 51.3 benzene 0.888 2.57 1.51 1.46 146 390 3.71 159 161 342 52.2 62.6 2.20 3.44 0.188 1.81 1,3-butdiene nd nd 0.002 nd 4.39 7.92 nd 0.582 0.216 2.94 nd nd 0.153 0.007 nd nd xylenes (o+m+p) 0.100 0.445 nd 0.674 0.857 3.05 0.486 0.338 nd 0.735 nd 0.052 nd 2.32 0.080 0.635 styrene nd nd nd nd nd nd nd nd nd nd nd nd nd 7.56 nd 0.025 ethne 0.281 1.30 0.058 1.39 372 258 1.71 11.1 3.83 23.1 20.8 9.42 0.109 0.031 0.175 0.286 ethylene 0.494 1.58 0.221 1.75 461 1061 7.90 155 157 467 151 108 9.06 0.038 0.309 0.207 cetylene 0.476 1.14 0.067 0.337 132 305 5.52 304 223 564 121 173 4.77 0.071 0.025 0.043 propne 0.058 0.345 0.017 0.469 134 57.7 0.771 1.85 1.40 5.04 2.62 4.22 0.062 0.084 0.047 0.155 propene 0.092 0.764 0.147 0.747 167 166 1.86 13.5 8.01 61.1 15.5 18.4 1.70 0.231 0.064 0.036 i-butne 0.010 0.138 nd 0.070 10.3 3.70 0.126 nd nd 0.237 0.230 0.220 0.069 0.104 0.006 0.047 i-butene 0.016 0.207 0.039 0.108 13.2 nd 0.087 0.681 nd 3.15 0.570 1.24 0.175 0.059 0.009 0.005 1-butene 0.007 0.086 0.013 0.077 37.3 30.0 0.225 1.60 1.08 7.80 2.39 3.83 0.498 0.099 0.005 0.002 n-butne 0.014 0.144 nd 0.167 40.4 15.9 0.266 0.461 0.191 0.620 0.461 0.688 0.011 0.023 0.013 0.057 trns-2-butene nd 0.073 0.037 0.061 5.90 9.55 0.066 5.07 2.41 10.6 1.80 2.62 0.092 0.047 0.003 0.034 2,2-dimethylpropne nd nd nd nd nd nd nd nd nd nd nd nd nd 0.069 nd nd cis-2-butene nd 0.034 0.009 0.056 1.14 nd 0.251 0.240 nd 1.04 0.234 0.491 0.041 0.138 0.035 0.202 3-methyl-1-butene nd nd nd 0.010 2.95 0.836 0.026 nd nd 0.715 0.182 0.207 0.041 0.005 0.003 nd i-pentne 0.009 0.230 nd 0.104 8.94 2.84 0.040 nd nd 0.092 0.217 0.011 nd nd 0.004 0.015 1-pentene nd 0.043 0.009 0.024 0.360 nd 0.087 0.268 nd 1.21 nd 0.450 0.217 nd 0.003 0.044 2-methyl-1-butene nd 0.043 nd 0.014 2.69 0.522 nd nd nd 0.604 nd 0.091 0.049 0.005 0.001 nd n-pentne 0.013 0.103 0.009 0.091 22.7 8.97 0.096 nd nd 0.244 0.285 0.015 nd nd 0.004 0.020 trns-2-pentene nd nd nd 0.024 0.338 nd 0.006 nd nd 0.600 nd 0.203 0.025 nd nd nd cis-2-pentene nd nd nd 0.010 0.235 nd nd nd nd 0.392 nd 0.131 0.015 0.005 nd nd 2-methyl-2-butene nd nd nd 0.031 0.300 nd nd nd nd 0.032 nd 0.435 0.020 0.020 nd nd cyclopentene nd nd nd nd nd nd nd nd nd 0.108 nd nd 0.035 nd nd nd 4-methyl-1-pentene nd nd nd nd nd nd nd nd nd 0.144 nd nd nd nd nd nd cyclopentne nd 0.004 nd 0.002 1.44 nd 0.009 nd nd nd nd nd nd nd 0.001 nd 2,3-dimethylbutne nd 0.022 nd 0.009 2.54 1.19 nd nd nd nd nd nd nd nd nd nd 2-methylpentne 0.012 0.138 nd 0.038 5.99 1.82 0.018 nd nd 0.040 nd nd nd nd 0.002 0.012 3-methylpentne 0.010 0.069 nd 0.028 1.62 0.532 0.008 nd nd nd nd nd nd nd nd nd 1-hexene nd 0.034 nd 0.035 8.99 4.14 0.108 nd nd 1.25 nd 0.416 0.214 0.020 0.007 nd n-hexne 0.025 0.121 0.009 0.077 15.6 6.16 0.058 nd nd 0.079 0.107 0.080 0.002 nd 0.002 0.011 methylcyclopentne nd 0.063 nd 0.012 4.36 1.58 nd nd nd nd nd nd nd nd nd 0.002 2,4-dimethylpentne nd 0.086 nd 0.010 nd nd nd nd nd nd nd nd nd nd nd nd cyclohexne nd nd nd nd 0.134 0.508 nd nd nd nd nd nd nd 0.341 nd nd 2-methylhexne 0.010 nd nd 0.019 1.26 0.365 nd nd nd nd nd nd nd nd nd 0.002 2,3-dimethylpentne nd nd nd 0.002 1.87 0.372 nd nd nd nd nd nd nd nd nd 0.007 3-methylhexne 0.024 nd nd 0.030 1.17 0.416 0.052 nd nd 0.013 nd nd nd nd 0.001 0.022 2,2,4-trimethylpentne 0.144 nd nd nd 2.17 1.30 nd nd nd 0.341 nd nd 0.151 nd nd nd n-heptne 0.034 0.101 0.020 0.083 10.6 4.17 0.061 nd nd 0.042 nd 0.060 0.030 nd nd 0.007 methylcyclohexne 0.019 0.109 nd 0.041 2.03 1.16 0.001 nd nd nd nd nd 0.069 nd nd nd 2,3,4-trimethylpentne nd nd nd nd nd nd nd nd nd nd nd nd 0.010 nd nd nd toluene 0.207 0.787 0.258 0.732 22.4 59.7 1.03 11.2 7.39 29.5 1.33 3.07 0.397 1.40 0.085 0.625 2-methylheptne nd 0.004 nd 0.023 1.31 nd nd nd nd 0.020 nd nd 0.022 nd nd nd 3-ethylhexne nd 0.034 nd 0.021 0.229 nd nd nd nd nd nd nd 0.039 nd nd nd n-octne 0.033 0.161 0.042 0.076 3.21 1.39 0.040 nd nd 0.002 0.210 nd 0.059 nd 0.001 0.008 ethylbenzene 0.046 0.115 0.010 0.089 nd 0.489 0.102 0.293 nd 0.195 nd 0.046 nd 0.160 0.006 0.057 n-nonne 0.103 0.305 0.189 0.194 0.333 0.370 0.006 nd nd nd nd nd nd 0.170 nd nd i-propylbenzene nd nd nd 0.051 nd nd nd nd nd nd nd nd nd nd 0.000 nd n-propylbenzene nd nd nd nd nd nd nd nd nd nd nd nd nd nd 0.015 nd p-ethyltoluene nd nd nd nd nd nd nd nd nd nd nd nd nd nd 0.009 nd m-ethyltoluene nd nd nd nd nd nd nd nd nd nd nd nd nd nd 0.001 nd 1,2,4-trimethylbenzene nd 0.607 nd 0.205 nd nd nd 0.212 nd nd nd nd nd nd nd 0.001 nd sec-butylbenzene n-decne 0.175 0.408 0.411 0.351 nd nd nd 157 nd nd nd nd nd 1.05 0.004 nd totl specited NMHC 3.91 11.8 3.28 9.60 1650 2407 24.9 822 565 1524 371 389 20.3 17.5 1.11 4.37 Key: nd ) not detected or bckground level g flue gs concentrtion. b Averge of 3 burning cycles using Honey- combintion. blnce equtions. The mesurement errors ppered to be greter for those compounds present t level below or close to the method detection limit. The gs stoves emitted more compounds t low concentrtion levels. The estimted CV vlues of emission fctors reflected the overll contribution of these two uncertinty sources. This nlysis is supported by lrge CV vlues reported for the emission fctors of TNMHC for the 28 Chinese fuel/stove combintions including the 16 ones mesured for specited NMHCs (6). The CV vlues for the TNMHC emission fctors rnged from 13% to 173% cross the 28 fuel/stove combintions with the highest vlue being for the Honey-. The CV vlues for the specited NMHCs ppered to be within the CV vlue rnges for the TNMHC for the Honey-, the only fuel/stove combintions tht were repetedly mesured for both specited NMHCs nd TNMHC. [The sum of mesured individul NMHCs did not ccount for ll the HCs present in the smples. It represents, in theory, frction of TNMHC. However, TNMHC nd specited NMHCs were mesured using different methods (see ref 6), nd thus the sums of specited NMHCs cnnot be directly compred with the concentrtions of TNMHC. Mesurement precisions for TNMHC nd specited NMHCs were, however, expected to be similr.] Bsed on these dt, we expect the emission fctors of specited NMHCs reported in this pper my hve lrge CV vlues, but probbly smller thn 173%. Another source of uncertinty nd reson to do more tests in field settings is tht wood fuels seem to hve higher emission rtes of PICs during the smoldering phse of combustion (17). Since smoldering cn continue in the field fter cooking 2874 9 ENVIRONMENTAL SCIENCE & TECHNOLOGY / VOL. 37, NO. 13, 2003

FIGURE 1. Ozone forming potentil (mg ozone/mj) for the 16 fuel/stove combintions (log-scle). The stndrd devition br for Honey- ws derived from three repeted mesurements. No repeted mesurements were mde for the other 15 fuel/stove combintions. session is completed, the emissions from the smoldering re difficult to mesure in simulted settings. This is only serious issue with wood fuels, however, s our tests of the Chinese fuel/stove combintions did include smoldering phses for col nd crop residues do not smolder significntly fter fueling stops. Gses nd liquids, of course, hve no such smoldering phses. Specited NMHCs were lso mesured in the pilot study in which 6 fuel/stove combintions in the Philippines were tested using the crbon blnce pproch (2, 4, 18). The 6 fuel/stove combintions tested included the stoves burning LPG, kerosene, chrcol, nd wood fuels. As it ws typicl in Mnil, none of the tested stove hd flues or chimneys. Compred to the extensive study (5, 6), the pilot study mesured fewer prmeters tht needed to construct the complete crbon blnce model nd did not use hood to promote mixing of flue gses. No repeted mesurements were mde in the pilot study. Hence, the results from the pilot study cn be expected to be less ccurte thn those from the extensive study. Nevertheless, both the pilot study nd the present study showed tht burning wood fuels yielded greter emissions of benzene (0.2-0.4 g s crbon per kg dry fuel) thn burning other fuels. 1,3-Butdiene emissions mesured in the present study, however, were substntilly lower (by 1 or 2 orders of mgnitude) thn those mesured in the pilot study for the similr fuel/stoves tested in both studies. Reltively lrge differences in other NMHC emissions were lso found when the results from the two studies were compred. A very limited number of other studies of cookstove emissions only mesured certin specific compounds of interest nd usully did not report NMHCs. Therefore, it is necessry to conduct future systemtic studies nd uncertinty nlyses to better understnd the nturl vribility nd to reduce the mesurement errors ssocited with emission fctors for vrious fuel-stove systems. Ozone Concerns. Locl nd regionl pollution of ozone nd photochemicl smog is common worldwide nd occurs in mny urbn nd suburbn res in Chin. Among the specited NMHCs identified in this study, some compounds hd higher ozone forming potentils thn others. Sturted orgnic compounds such s lknes usully hve lower MIR vlues thn lkenes nd romtic species tht re more ctive in terms of forming ozone. When both MIR vlues nd mount of emissions were considered, ethylene, mong the 54 specited NMHCs, ws the mjor contributor to OFP vlues for most fuel-stove combintion tested in this study. In most fuel/stove combintions, emission of ethylene contributes more thn 50% of ozone (rnging from 12% to 85%), followed by cetylene (1-46%), xylenes (0-32%), styrene (0-31%), propene (10-24%), cis-2-butene (0-19%), 1,2,4-trimethylbenzene nd sec-butylbenzene (0-18%), nd toluene (0-10%). The type nd mount of fuels used in Chinese households vries in urbn nd rurl res (11, 12). In urbn res of Chin, the mjor types of fuel used in the residentil sector re col, electricity, LPG, nturl gs, nd col gs. In rurl res, biomss fuels (crop residues nd wood) re the mjor types of fuel used in the residentil sector. Becuse groundlevel ozone ccumultion results from complex chemicl rections, in the tmosphere, between NMHCs nd NO x, the ctul ozone formtion my be limited by NO x concentrtions in rurl res where no significnt sources of NO x re present. Therefore, here we only discuss ozone-forming potentils for urbn res. However, some NMHCs emitted from rurl household stoves cn be trnsported to urbn res, contributing to the urbn ozone pollution. Tble 5 shows the estimted totl ozone forming potentils from domestic cookstoves in the urbn res of Chin bsed on the fuel usge nd the specific OFP vlue for ech type of fuel mesured in this study. The estimtes in Tble 5 were VOL. 37, NO. 13, 2003 / ENVIRONMENTAL SCIENCE & TECHNOLOGY 9 2875

TABLE 5. Estimted Yerly Ozone Forming Potentil (OFP) in the Urbn Are of Chin relted to household cookstoves fuel types v OFP (g ozone /kg fuel) mount of consumption in the urbn re of Chin in 1995 (metric ton 10 6 ) estimted totl OFP per yer in the urbn re (ton ozone) 10 2 rw col 6.03 58.3 b 3515 col briquette 0.03 0.25 0.08 kerosene 0.74 0.07 0.52 LPG 0.144 5.10 7.34 col gs 0.031 2.37 0.74 nturl gs 0.064 1.68 1.07 totl 3525 Reference 11. b Totl mount of rw col nd wshed col consumed. fuel specific rther thn fuel/stove specific, becuse dt on fuel consumption were only fuel specific. We brek down the col consumption dt by col type to seprte rw col nd col briquette (including honeycomb col briquette) consumptions. Our estimtes indicte tht the totl ozone mount ttributble to urbn residentil NMHC emissions in 1995 were pproximtely 353 000 tons per yer. The results indicted tht col-fired cookstoves would contribute > 99% of totl ozone mount ttributble to NMHC emissions from ll the cookstoves used in urbn Chin. This lrge col contribution ws minly driven by the high OFP vlue for rw col nd the lrge mount of rw col consumption. However, relible dt sources of col briquette consumption or production re uncertin becuse most col briquettes were produced in smll workshops where no effective reporting or trcking systems were in plce to record the production dt. We suspect tht the production of col briquettes might be substntilly underestimted in the Chin ntionl energy sttistics nd, hence, the OFP estimte for col briquette in Tble 5 my be lower thn the ctul vlue. Bsed on the sttisticl informtion from the Chin Energy Yerbook (11, 12), the rw col usge hs decresed grdully yerly nd replced by col briquettes, kerosene, nd gseous fuels. However, since the mount of rw col consumption is much lrger nd rw col hs the highest OFP vlue mong ll the fuel tested, the totl OFP from the urbn residentil sector (for cooking) should hve chnged insignificntly over the pst few yers. Our nlysis shows, for exmple, in 1997 rw col still contributes more thn 95% of totl OFP ttributble to household cookstove emissions of NMHCs. This type of OFP nlysis using the mesured emission fctors of NMHCs my provide useful informtion for energy decision-mking, from stndpoint of ozone control strtegy, in Chin or other countries hving similr suite of fuel/stove combintions. By switching the fuel/stove combintions, the NMHC contribution to ozone formtion cn be significntly reduced (see Figure 1). For exmple, on verge, switch of col powder to nturl gs would led to reduction of OFP by 70 times, nd switch of col powder to col gs would led to 40-fold OFP reduction. Among different types of col fuels, burning col briquettes nd honeycomb col briquettes contribute much less OFPs thn burning col powder. However, our estimtes were subject to lrge uncertinties due to the uncertinty ssocited with the emission fctors mesured, the uncertinty ssocited with fuel consumption dt, nd the uncertinty ssocited with MIR vlues (difference in tmospheric NO x/nmhc conditions between the United Sttes nd Chin urbn res). It will probbly be necessry to conduct mesurements under field conditions to reduce these uncertinties sufficiently to mke definitive policy decisions bout ozone control policy in Chin. Other Potentil Applictions of the Dtbse. The specited NMHC dtbse cn be useful in source pportionment studies. Emissions from numerous cooking stoves my be significnt source to prticulte mtter nd voltile orgnic compounds nd thus significntly influence the locl nd regionl ir qulity including the formtion of photochemicl smog. Sorting out the reltive importnce of sources is lwys ttrctive to the scientific community nd policy mkers. Source pportionment techniques hve been widely used to estimte reltive contributions of vrious sources to the concentrtions of given (trget) pollutnts mesured t the receptor. Using VOCs, similr to wht we mesured in the study, to do source pportionment hs been done in previous studies (19, 20). For exmple, Anderson et l. (20) pplied personl exposure mesurements for VOCs to four receptor-oriented source pportionment models. The previous studies used the VOC dt collected in the Totl Exposure Assessment Methodology (TEAM) studies conducted in New Jersey nd Cliforni s well s the dt collected in the Cliforni Indoor Exposure study to evlute severl receptor models. The models evluted include the Chemicl Mss Blnce model, Principl Component Anlysis/Absolute Principl Component Scores, Positive Mtrix Fctoriztion, nd Grphicl Rtio Anlysis for Composition Estimtes/Source Apportionment by Fctors with Explicit Restriction. The source pportionment results from the four models greed resonbly well for the New Jersey dt but were less consistent for the Cliforni dt. Model performnce vries from severl fctors such s vilbility of informtion on emission composition, number of receptor smples, rectivity of compounds used s trcers, mesurement uncertinty, frction of missing dt, frction of belowdetection-limit dt, number of fctors determined, fctors interprettion for sources, whether ll importnt sources re included, nd whether there is collinerity within/mong the source profiles used in the model. Nevertheless, this indictes tht n ccurte source profile is the key to the successful implementtion of source pportionment models. A previous nlysis by Zhng nd Smith (2) evluted lifetime cncer risks ssocited with exposure to benzene, 1,3-butdiene, styrene, nd xylenes relesed from the 6 fuel/ stove combintions tested in the Mnil pilot study. The nlysis found tht estimted cncer risk of benzene nd tht of styrene from the use of the tested wood stove could exceed published risk estimtes from ll sources of irborne benzene or styrene (excluding ctive tobcco smoking) in the United Sttes. These kind of importnt findings bsed on the limited dt of the pilot study should be further evluted using similr or more sophisticted nlysis bsed on the more extensive dt collected in this study. In the future nlysis, cncer risks for ll the tested fuel/stove combintions cn be compred to mke quntittive recommendtions on fuel or stove switching (2). Additionl nlyses my be preformed to lso ssess the noncncer helth risks. For exmple, this dtbse my be used to conduct chemicl mixture risk ssessment using groups (or clsses) of hydrocrbons bsed on crbon number frctions (e.g., C 5-C 8,C 9-C 16), their fte nd trnsport, nd toxicity chrcteristics. These will be discussed in more detiled helth risk nlysis pper to follow. Acknowledgments We gretly pprecite the effective nd timely peer review process of our pper. We pprecite the comments from Dr. Rufus Edwrds of University of Cliforni t Berkeley nd editoril ssistnce from Dr. Steven Miller of the New Jersey Deprtment of Helth nd Senior Services nd Mr. Robert C. Hrrington of Environmentl nd Occuptionl Helth Sciences Institute, New Jersey. We pprecite the ssistnce from Dr. Jonthn Sinton of Lwrence Berkeley Ntionl Lbortory in providing informtion on col consumption in Chin. The originl dt collection ws funded through 2876 9 ENVIRONMENTAL SCIENCE & TECHNOLOGY / VOL. 37, NO. 13, 2003