DEVELOPMENT OF THE SAPRC-07 CHEMICAL MECHANISM AND UPDATED OZONE REACTIVITY SCALES

                 Report to the California Air Resources Board
                   Contracts No. 03-318, 06-408, and 07-730

                                      By

                             William P. L. Carter

                           Revised January 27, 2010

               Center for Environmental Research and Technology
                            College of Engineering
                           University of California
                          Riverside, California 92521

ABSTRACT
A completely updated version of the SAPRC-99 chemical mechanism, designated SAPRC-07, has been developed and is documented in this report. This includes a complete update of the rate constants and reactions based on current data and evaluations, reformulated and less parameterized aromatics mechanisms, a representation of chlorine chemistry, a reformulated method to represent peroxy reactions that is more appropriate for modeling secondary organic aerosol formation, and improved representations for many types of VOCs. This mechanism was evaluated against the result of ~2400 environmental chamber experiments carried out in 11 different environmental chambers, including experiments to test mechanisms for over 120 types of VOCs. The performance of the mechanism in simulating the chamber data was comparable to SAPRC-99, with generally satisfactory results for most types of VOCs but some increases in biases in simulations of some mixture experiments. The mechanism was used to derive an update to the MIR and other ozone reactivity scales for over 1100 types of VOCs. The average changes in relative MIR values was about 10%, with >90% of the VOCs having changes less than 30%, but with larger changes for some types of VOCs, including halogenated compounds. Recommendations are given for future mechanism development research.
The mechanism documentation includes some large tabulations that are being provided only in electronic form. Links to downloading these tabulations are available at http://www.cert.ucr.edu/~carter /SAPRC.
Note: This version of the report and the associated supplementary materials and files superceded those available previously at the above-referenced web site. The changes are summarized in Appendix E to this report.
ACKNOWLEDGEMENTS AND DISCLAIMERS
This work was carried out at the College of Engineering Center for Environmental Research and Technology (CE-CERT) at the University of California at Riverside (UCR), with funding by the California Air Resources Board (CARB) primarily through contract number 03-318, with additional work carried out through Contracts 06-408 and 07-730. However, the contents of this report reflect only the opinions and conclusions of the author, and not CE-CERT, UCR, or the CARB.
The author wishes to thank Dr. Roger Atkinson of the Air Pollution Research Center at UCR for many helpful discussions concerning available kinetic and mechanistic data and atmospheric chemistry in general, and for partial support during the period of this project. The author also appreciates the patience of the staff of the Research Division of the CARB, particularly Dr. Dongmin Luo, regarding the delays in completing this project and their continuing support. Experiments for this project that were used in this mechanism evaluation were carried out primarily by Ms. Irina Malkina with assistance from Mr. Kurt Bumiller and other members of the CE-CERT Atmospheric Processes Laboratory staff, and their contributions are gratefully acknowledged. Finally, the author wishes to acknowledge the University of California Retirement System, which provided major support to the author during the completion of mechanism development and the preparation of this report.
When discussing materials whose mechanisms were evaluated, this report makes references to "Texanol", which is a trademark of Eastman Chemical Company, and "Exxol" and "Isopar", which are trademarks of ExxonMobil Chemical Company. Mention of trade names and commercial products does not constitute endorsement or recommendation for use.

TABLE OF CONTENTS
EXECUTIVE SUMMARY	1
Background	1
Accomplishments	1
Results	2
Recommendations	3
Introduction	5
Mechanism Documentation	8
Base Mechanism	8
Listing of Standard Base Mechanism and Summary of Changes	8
Representation of Peroxy Radical Operators	13
Base Chlorine Mechanism	15
Representation of Organic Products	16
Mechanisms for Individual VOCs	17
Mechanism Generation and Estimation System	17
Enhanced Capabilities for VOCs and Reactions	18
Support for Chlorine Reactions	19
Updates to Rate Constants and Mechanisms	26
Adjusted Product Mechanisms	28
Updated Aromatics Mechanisms	31
Reactions of Alkylbenzenes with OH	32
Representation of Unsaturated Dicarbonyl Products	41
Aromatic Ring Retaining Products	49
Non-Alkylbenzene Aromatics	50
Reactions with Chlorine	50
Mechanisms for Amines	54
Miscellaneous Assigned Mechanisms	56
Lumped Molecule Representations	59
Lumped Mechanisms for Airshed Models	60
Adjustable Parameter Mechanisms	60
Fixed Parameter Mechanism	60
Emissions Assignments	64
Mechanism Evaluation	66
Methods	70
Data Presented	76
Results	78
Results for Chamber Characterization Experiments	78
Results for Mechanism Evaluation Experiments	79
Overall Performance in Simulating Entire Dataset	79
Evaluation Experiments for Individual VOCs	79
Evaluation Experiments for Mixtures	91
Mechanism Performance for Direct Reactivity Experiments	99
Updated Reactivity Scales	100
Methods	100
Results	102
Predicted Ozone Formation in the Scenarios	102
Updated Ozone Reactivity Scales	104
Updated Uncertainty Estimates	105
Summary and recommendations	108
Summary of Project Accomplishments	108
Recommendations	113
References	120
Appendix A. General Mechanism Listing Tables	136
Appendix B. Mechanism and Reactivity Listings for Individual VOCs	177
Appendix C. Mechanism Evaluation Tabulations and Plots	271
Appendix D. Supplementary Material available	377
Appendix E. Change Log	378

LIST OF TABLES

Table 1.	Reactions where the rate constants or photolysis rates changed by more than 5% or that were added or removed in the standard base mechanism relative to SAPRC-99.	11
Table 2.	Products formed in the reactions of the chemical operators used to represent peroxy radical reactions in the SAPRC-07 mechanism.	14
Table 3.	Group additivity parameters used for estimating rate constants for abstraction reactions by OH and NO3 radicals and by Cl atoms.	20
Table 4.	Rate constants for reactions with chlorine atoms for organic compounds used to derive parameters for group-additivity estimates.	24
Table 5.	Group additivity rate constants and factors used for estimating rates of addition of Cl atoms to double bonds.	26
Table 6.	Branching ratio assignments for chlorine-substituted alkoxy radicals made to generate mechanisms for the reactions of chlorine with representative alkenes.	27
Table 7.	Assigned mechanisms for Crigiee Biradicals used for deriving mechanisms for reactions of O3 with alkenes.	29
Table 8.	Compounds whose Averaged Conditions MIR reactivities changed by more than 20% when adjustable product mechanisms were employed.	31
Table 9.	Measured or estimated yields of known products for benzene and methylbenzene reactions that were incorporated in the SAPRC-07 aromatics mechanisms	33
Table 10.	Rate constants and assumptions used to estimate fractions of reaction of on and off the aromatic ring for the reactions of OH with aromatic compounds or groups of compounds for which mechanistic parameters have been derived.	39
Table 11.	Model species used to represent non-nitrate products from the reactions of OH radicals at aromatic rings for various types of ring structures. Yields are normalized to OH ring addition only, excluding nitrate formation.	42
Table 12.	Estimated yields for unsaturated dicarbonyl products in the OH reactions of the methylbenzenes for which -dicarbonyl yield data are available.	43
Table 13.	Mechanism estimation assignments made for representative unsaturated dicarbonyl compounds for the purpose of generating mechanisms for the aromatic fragmentation product model species AFG1-3.	45
Table 14.	Discussion of mechanisms for the non-alkylbenzene aromatic compounds that for which mechanistic assignments have been derived. The rate constants and additional documentation are given in Table B-4 and the reactions are given in Table B-2 in Appendix B.	51
Table 15.	Assignments made for alkoxy radical intermediates predicted to be formed in the photooxidation mechanisms of the amines.	55
Table 16.	Discussion of mechanisms for miscellaneous compounds that were not derived using the procedures discussed previously in this report. The rate constants, reactions, and additional documentation are given in Table B-2, and Table B-4 through Table B-6 in Appendix B.	57
Table 17.	List of model species used in the lumped mechanism for airshed models for representing the various types of emitted VOCs.	61
Table 18.	Compounds and weighting factors used to derive the parameters for the ALK and ARO model species in the fixed parameter mechanism for airshed models	63
Table 19.	Compounds and weighting factors used to derive the parameters for the OLE and TERP model species in the fixed parameter mechanism for airshed models	64
Table 20.	Summary of environmental chambers whose data were used for mechanism evaluation	67
Table 21.	Summary of types of experiments used for mechanism evaluation.	69
Table 22.	Chamber wall effect and background characterization parameters used in the environmental chamber model simulations for mechanism evaluation.	72
Table 23.	Summary of mechanism evaluation results using chamber experiments for individual compounds and mixtures.	83
Table 24.	Scenarios used for reactivity assessment, with updated calculated maximum O3, Integrated OH, and MIR, MOIR, and EBIR NOx inputs.	101
Table 25.	List of compounds and mixtures whose reactivities in the MIR scale changed by more than 30%.	106

Table A-1.	List of model species used in the standard base SAPRC-07 mechanism. The lumped model species added for the fixed parameter lumped mechanism for airshed models are also listed.	136
Table A-2.	Listing of reactions and rate parameters in the base SAPRC-07 mechanism.	140
Table A-3.	Absorption cross-sections and quantum yields for the all photolysis reactions in the base mechanism. (Available in electronic form only)	163
Table A-4.	List of model species used in the base chlorine SAPRC-07 mechanism.	163
Table A-5.	Listing of reactions and rate parameters added to represent the reactions of chlorine species in the SAPRC-07 mechanism.	164
Table A-6.	Absorption cross-sections and quantum yields for the all photolysis reactions added for the chlorine chemistry mechanism.	170
Table A-7.	Listing of reactions and rate parameters used for the lumped model species in the fixed parameter version of the lumped SAPRC-07 mechanism.	173

Table B-1.	Listing of detailed model species, their representation in the model, and calculated atmospheric reactivity values in various reactivity scales	178
Table B-2.	Listing of mechanisms for all VOCs for which mechanism assignments have been derived. (Available in electronic form only)	210
Table B-3.	Listing of adjusted product mechanisms for all VOCs for which such mechanisms have been derived. (Available in electronic form only)	210
Table B-4.	Listing of compounds for which mechanisms have been derived. 1. OH radical rate constants at 300[o]K, mechanism types, and structures used for compounds with generated mechanisms.	211
Table B-5.	Listing of compounds for which mechanisms have been derived. 2. Rate constants for reactions with O3, NO3, and O[3]P at 300[o]K, where applicable.	236
Table B-6.	Listing of compounds for which mechanisms have been derived. 3. Absorption cross-sections, quantum yields, and mechanisms for photolysis reactions, where applicable.	242
Table B-7.	Listing of compounds for which mechanisms have been derived. 4. Rate constants and mechanisms for reactions with chlorine atoms, where used.	245
Table B-8.	Absorption cross sections used for the photolysis reactions for represented photoreactive VOCs that are not in the base mechanism.	253
Table B-9.	Lumped molecule representations used in the SAPRC-07 mechanism.	257
Table B-10.	Listing of components of complex mixtures for which atmospheric reactivities have been calculated. (Available in electronic form only)	270
Table B-11.	Assignments of explicitly represented compounds to lumped model species in the fixed parameter mechanism. (Available in electronic form only)	270

Table C-1.	List of environmental chamber experiments used in the mechanism evaluation. . (Available in electronic form only)	271

LIST OF FIGURES

Figure 1.	Effect of using the adjustable product mechanisms vs. the standard fixed product mechanisms on incremental reactivities in the "Averaged Conditions" MIR scale for those VOCs for which adjustable product mechanisms could be derived.	30
Figure 2.	Reactions of the OH - aromatic ring adduct that are assumed in the current mechanism. The example shown is for benzene, but analogous reactions are assumed for the alkylbenzenes.	36
Figure 3.	Nitrate yield factors used for estimating aromatic oxidation mechanisms	38
Figure 4.	Action spectra for photodecompositions of 2-buene-1,4-dial, 4-oxo-2-penteneal and used for the model species AFG1 and AFG2.	49
Figure 5.	Distribution of ([O3]-[NO]) model errors for the radical source characterization experiments.	80
Figure 6.	Distribution of O3 model errors for the NOx offgasing characterization experiments	81
Figure 7.	Distributions of ([O3]-[NO]) model errors for all of the single VOC - NOx and mixture - NOx experiments used in the mechanism evaluation.	82
Figure 8.	Experimental and calculated effects of CO addition in representative toluene and m-xylene NOx experiments carried out in the UCR EPA chamber.	90
Figure 9.	Plots of model errors in simulations of miscellaneous mixture - NOx experiments carried out in various chambers.	93
Figure 10.	Plots of model errors in simulations of miscellaneous surrogate - NOx experiments carried out in various UCR chambers and in the TVA chamber.	94
Figure 11.	Plots of model errors in simulations of surrogate - NOx experiments carried out in the UNC chamber.	95
Figure 12.	Plots of model errors in simulations of the standard base case surrogate - NOx irradiations carried out in conjunction with the incremental reactivity experiments in the ETC and DTC (blacklight), and CTC and XTC (arc light)  chambers.	96
Figure 13.	Plots of model errors in simulations of the surrogate - NOx experiments carried out in the UCR EPA chamber.	97
Figure 14.	Plots of model errors in simulations of 6-hour ([O3]-[NO]) against initial ROG surrogate /NOx ratio.	98
Figure 15.	Model performance in simulating the direct reactivity data of Carter and Malkina (2002).	99
Figure 16.	Plots of maximum O3 concentrations calculated for the various scenarios used for reactivity assessment against the initial ROG levels in the scenarios.	103
Figure 17.	Change in maximum O3 concentrations in the scenarios used for reactivity assessment, calculated using SAPRC-07 compared to SAPRC-99.	104
Figure 18.	Plots of incremental reactivities in the MIR scale computed using the updated mechanism against those reported previously by Carter (2003a).	105

Figure C-1.	Plots of experimental and calculated environmental chamber reactivity results for carbon monoxide.	272
Figure C-2.	Plots of experimental and calculated environmental chamber reactivity results for ethane and propane	273
Figure C-3.	Plots of experimental and calculated environmental chamber reactivity results for n-butane	274
Figure C-4.	Plots of experimental and calculated environmental chamber reactivity results for n-hexane and the Surg-3 MIR1 experiments for n-octane.	275
Figure C-5.	Plots of experimental and calculated environmental chamber reactivity results for n-octane (excluding the Surg-3 MIR1 experiments, which are on Figure C-4).	276
Figure C-6.	Plots of experimental and calculated environmental chamber reactivity results for n-dodecane.	277
Figure C-7.	Plots of experimental and calculated environmental chamber reactivity results for n-tetradecane	278
Figure C-8.	Plots of experimental and calculated environmental chamber reactivity results for n-pentadecane and n-hexadecane.	279
Figure C-9.	Plots of experimental and calculated environmental chamber reactivity results for isobutane, 2,2,4-trimethyl pentane, and 2,6-dimethyl octane.	280
Figure C-10.	Plots of experimental and calculated environmental chamber reactivity results for 2-methyl nonane and 3,4-diethyl hexane.	281
Figure C-11.	Plots of experimental and calculated environmental chamber reactivity results for cyclohexane and hexyl cyclohexane.	282
Figure C-12.	Plots of experimental and calculated environmental chamber reactivity results for octyl cyclohexane.	283
Figure C-13.	Plots of model errors in simulations of the ethene - NOx environmental chamber experiments.	284
Figure C-14.	Plots of experimental and calculated environmental chamber reactivity results for ethene.	285
Figure C-15.	Plots of model errors in simulations of the propene - NOx environmental chamber experiments.	286
Figure C-16.	Plots of experimental and calculated environmental chamber reactivity results for propene.	287
Figure C-17.	Plots of model errors in simulations of the 1-buteme and 1-hexene - NOx environmental chamber experiments.	288
Figure C-18.	Plots of model errors in simulations of the trans-2-butene and isobutene - NOx environmental chamber experiments.	289
Figure C-19.	Plots of experimental and calculated environmental chamber reactivity results for isobutene and trans-2-butene.	290
Figure C-20.	Plots of model errors in simulations of the isoprene - NOx environmental chamber experiments.	291
Figure C-21.	Plots of experimental and calculated environmental chamber reactivity results for isoprene	292
Figure C-22.	Plots of model errors in simulations of the terpene - NOx environmental chamber experiments.	293
Figure C-23.	Plots of experimental and calculated environmental chamber reactivity results for - and -pinene and d-limonene.	294
Figure C-24.	Plots of experimental and calculated environmental chamber reactivity results for cyclohexene.	295
Figure C-25.	Plots of experimental and calculated environmental chamber reactivity results for styrene.	295
Figure C-26.	Plots of model errors in simulations of the benzene and benzene + CO - NOx environmental chamber experiments.	296
Figure C-27.	Plots of experimental and calculated environmental chamber reactivity results for benzene.	297
Figure C-28.	Plots of model errors in simulations of the toluene and toluene + CO - NOx environmental chamber experiments.	298
Figure C-29.	Plots of model errors in simulations of the ethylbenzene - NOx environmental chamber experiments.	299
Figure C-30.	Plots of experimental and calculated environmental chamber reactivity results for toluene and ethylbenzene.	300
Figure C-31.	Plots of model errors in simulations of the m-xylene and m-xylene + CO - NOx environmental chamber experiments.	301
Figure C-32.	Plots of model errors in simulations of the o-xylene - NOx environmental chamber experiments.	302
Figure C-33.	Plots of model errors in simulations of the p-xylene - NOx environmental chamber experiments.	303
Figure C-34.	Plots of experimental and calculated environmental chamber reactivity results for m-xylene (additional experiments shown on Figure C-35).	304
Figure C-35.	Plots of experimental and calculated environmental chamber reactivity results for m-xylene (with variable surrogate conditions), o-xylene, and p-xylene.	305
Figure C-36.	Plots of model errors in simulations of the 1,2,3-trimethylbenzene - NOx environmental chamber experiments.	306
Figure C-37.	Plots of model errors in simulations of the 1,2,4-trimethylbenzene - NOx environmental chamber experiments.	307
Figure C-38.	Plots of model errors in simulations of the 1,3,5-trimethylbenzene and 1,3,5-trimethylbenzene + CO - NOx environmental chamber experiments.	308
Figure C-39.	Plots of experimental and calculated environmental chamber reactivity results for 1,2,3-, 1,2,4- and 1,3,4-trimethylbenzenes.	309
Figure C-40.	Plots of model errors in simulations of the naphthalene, 2,3-dimethylnaphthalene, and tetralin - NOx environmental chamber experiments.	310
Figure C-41.	Plots of model errors in simulations of the acetylene - NOx environmental chamber experiments.	311
Figure C-42.	Plots of experimental and calculated environmental chamber reactivity results for acetylene.	312
Figure C-43.	Plots of experimental and calculated environmental chamber reactivity results for methanol and ethanol.	313
Figure C-44.	Plots of experimental and calculated environmental chamber reactivity results for isopropyl alcohol.	314
Figure C-45.	Plots of experimental and calculated environmental chamber reactivity results for t-butyl alcohol.	315
Figure C-46.	Plots of experimental and calculated environmental chamber reactivity results for 1-, 2-, and 3-octanol.	316
Figure C-47.	Plots of experimental and calculated environmental chamber reactivity results for ethylene glycol.	317
Figure C-48.	Plots of experimental and calculated environmental chamber reactivity results for propylene glycol.	318
Figure C-49.	Plots of experimental and calculated environmental chamber reactivity results for dimethyl ether and diethyl ether.	319
Figure C-50.	Plots of experimental and calculated environmental chamber reactivity results for methyl t-butyl ether.	320
Figure C-51.	Plots of experimental and calculated environmental chamber reactivity results for 1-methoxy-2-propanol.	321
Figure C-52.	Plots of experimental and calculated environmental chamber reactivity results for 2-ethoxyethanol and 2-(2-ethoxyethoxy) ethanol (DGEE).	322
Figure C-53.	Plots of experimental and calculated environmental chamber reactivity results for 2-butoxyethanol.	323
Figure C-54.	Plots of experimental and calculated environmental chamber reactivity results for 2-(2-butoxyethoxy)-ethanol.	324
Figure C-55.	Plots of experimental and calculated environmental chamber reactivity results for methyl acetate.	325
Figure C-56.	Plots of experimental and calculated environmental chamber reactivity results for ethyl acetate.	326
Figure C-57.	Plots of experimental and calculated environmental chamber reactivity results for isopropyl and t-butyl acetates.	327
Figure C-58.	Plots of experimental and calculated environmental chamber reactivity results for methyl isobutyrate.	328
Figure C-59.	Plots of experimental and calculated environmental chamber reactivity results for methyl pivalate.	329
Figure C-60.	Plots of experimental and calculated environmental chamber reactivity results for n-butyl acetate.	330
Figure C-61.	Plots of experimental and calculated environmental chamber reactivity results for dimethyl carbonate and methyl isopropyl carbonate.	331
Figure C-62.	Plots of experimental and calculated environmental chamber reactivity results for propylene carbonate.	332
Figure C-63.	Plots of experimental and calculated environmental chamber reactivity results for 1-methoxy-2-propyl acetate.	333
Figure C-64.	Plots of experimental and calculated environmental chamber reactivity results for dimethyl succinate and dimethyl glutarate.	334
Figure C-65.	Plots of experimental and calculated environmental chamber reactivity results for the Texanol(R) isomers.	335
Figure C-66.	Plots of model errors in simulations of the furan - NOx environmental chamber experiments.	336
Figure C-67.	Plots of model errors in simulations of the methyl furan - NOx environmental chamber experiments.	337
Figure C-68.	Plots of model errors in simulations of the benzyl alcohol and benzyl alcohol + CO - NOx environmental chamber experiments.	338
Figure C-69.	Plots of experimental and calculated environmental chamber reactivity results for benzyl alcohol.	339
Figure C-70.	Plots of model errors in simulations of the formaldehyde - NOx environmental chamber experiments.	340
Figure C-71.	Plots of experimental and calculated environmental chamber reactivity results for formaldehyde.	341
Figure C-72.	Plots of model errors in simulations of the acetaldehyde - NOx environmental chamber experiments.	342
Figure C-73.	Plots of experimental and calculated environmental chamber reactivity results for acetaldehyde.	343
Figure C-74.	Plots of model errors in simulations of the acrolein and methacrolein (MA) - NOx environmental chamber experiments.	344
Figure C-75.	Plots of experimental and calculated environmental chamber reactivity results for benzaldehyde.	345
Figure C-76.	Plots of model errors in simulations of the acetone - NOx environmental chamber experiments.	346
Figure C-77.	Plots of experimental and calculated environmental chamber reactivity results for acetone.	347
Figure C-78.	Plots of model errors in simulations of the methyl ethyl ketone (MEK), 2-pentanone (MPK) and 2-heptanone - NOx environmental chamber experiments.	348
Figure C-79.	Plots of experimental and calculated environmental chamber reactivity results for methyl ethyl ketone.	349
Figure C-80.	Plots of experimental and calculated environmental chamber reactivity results for 2-pentanone.	350
Figure C-81.	Plots of experimental and calculated environmental chamber reactivity results for 4-methyl-2-pentanone and 2-heptanone.	351
Figure C-82.	Plots of experimental and calculated environmental chamber reactivity results for cyclohexanone.	352
Figure C-83.	Plots of model errors in simulations of the methyl vinyl ketone - NOx environmental chamber experiments.	353
Figure C-84.	Plots of experimental and calculated ([O3]-[NO]), cresol, and PAN concentrations in the cresol - NOx environmental chamber experiments.	354
Figure C-85.	Plots of experimental and calculated environmental chamber reactivity results for m-cresol.	354
Figure C-86.	Plots of experimental and calculated environmental chamber reactivity results for para toluene isocyanate.	355
Figure C-87.	Plots of experimental and calculated environmental chamber reactivity results for 2,4-toluene diisocyanate (TDI1) and 2,6-toluene diisocyanate (TDI2).	356
Figure C-88.	Plots of experimental and calculated environmental chamber reactivity results for n-methyl-2-pyrrolidone.	357
Figure C-89.	Plots of experimental and calculated environmental chamber reactivity results for n-propyl bromide and n-butyl bromide.	358
Figure C-90.	Plots of experimental and calculated environmental chamber reactivity results for trichloroethylene.	359
Figure C-91.	Plots of model errors in simulations of the 1,3-dichloropropene and 1,3-dichloropropene + n-butane - NOx environmental chamber experiments.	360
Figure C-92.	Plots of experimental and calculated environmental chamber reactivity results for the 1,3-dichloropropenes.	361
Figure C-93.	Plots of experimental and calculated environmental chamber reactivity results for 2-(chloromethyl)-3-chloro-propene.	361
Figure C-94	Plots of model errors in simulations of the chloropicrin (CP)  -  alkane - NOx and the chlorine + alkane - NOx experiments.	362
Figure C-95.	Plots of experimental and calculated environmental chamber reactivity results for chloropicrin.	363
Figure C-96.	Plots of experimental and calculated environmental chamber reactivity results for chlorine. Note that chlorine was added to the "base case" experiment after four hours of irradiation.	363
Figure C-97.	Plots of experimental and calculated environmental chamber reactivity results for hexamethyldisiloxane, d4 cyclosiloxane, and hydroxymethyldisiloxane.	364
Figure C-98.	Plots of experimental and calculated environmental chamber reactivity results for carbon disulfide and methyl isothiocyanate.	365
Figure C-99.	Plots of experimental and calculated O3, NO, DMSO and formaldehyde in the dimethyl sulfoxide - NOx environmental chamber experiments.	366
Figure C-100.	Plots of experimental and calculated environmental chamber reactivity results for dimethyl sulfoxide.	366
Figure C-101.	Plots of experimental and calculated environmental chamber reactivity results for s-ethyl dipropylthiocarbamate (EPTC).	367
Figure C-102.	Plots of experimental and calculated environmental chamber reactivity results for the mineral spirits samples studied for Safety-Kleen (Carter et al, 1997e).	368
Figure C-103.	Plots of experimental and calculated environmental chamber reactivity results for Exxon Exxol(R) D95 Fluid and Exxon Isopar(R) M Fluid studied for ExxonMobil (Carter et al, 2000e).	369
Figure C-104.	Plots of experimental and calculated environmental chamber reactivity results for oxo-decyl Acetate fluid studied for ExxonMobil (Carter et al, 2000e).	370
Figure C-105.	Plots of experimental and calculated environmental chamber reactivity results for VMP Naphtha sample studied by Carter and Malkina (2005) and the Kerosene sample studied by Carter and Malkina (2007).	371
Figure C-106.	Plots of experimental and calculated environmental chamber reactivity results for the dearomatized alkanes mixed, predominately C10-C12 (ASTM-1C), synthetic isoparaffinic alkane mixture, predominately C10-C12 (ASTM-3C1), Reduced Aromatics Mineral Spirits (ASTM-1B), and Regular mineral spirits (ASTM-1A) solvents studied by Carter and Malkina (2005).	372
Figure C-107.	Plots of experimental and calculated environmental chamber reactivity results for the aromatic-100 solvent studied by Carter and Malkina (2005).	373
Figure C-108.	Plots of experimental and calculated environmental chamber reactivity results for ethanolamine and isopropyl amine. In order to obtain the fits shown, the initial amine concentrations had to be adjusted to optimize the model fits to the reactivity data.	373
Figure C-109.	Plots of experimental and calculated environmental chamber reactivity results for t-butyl amine and 2-amino-2-methyl-1-propanol. In order to obtain the fits shown, the initial amine concentrations had to be adjusted to optimize the model fits to the reactivity data.	374
Figure C-110.	Plots of selected experimental and calculated concentration-time plots for the methyl iodide - NOx, methyl iodide - CO - NOx, and methyl iodide - O3 experiments.	375
Figure C-111.	Plots of experimental and calculated environmental chamber reactivity results for the methyl iodide experiments of Carter (2007).	376

EXECUTIVE SUMMARY
Background
Airshed models are essential for the development of effective control strategies for reducing photochemical air pollution because they provide the only available scientific basis for making quantitative estimates of changes in air quality resulting from changes in emissions. The chemical mechanism is an important component of the model that represents the processes by which emitted volatile organic compound (VOC) pollutants and oxides of nitrogen (NOx) react to form secondary pollutants such as ozone (O3) and other oxidants. If the mechanism is incorrect or incomplete in significant respects, then the model's predictions of secondary pollutant formation may also be incorrect, and its use might result in implementation of inappropriate or even counter-productive air pollution control strategies.
One airshed model application where the accuracy of the chemical mechanism is important is the calculation of reactivity scales that measure relative impacts of different types of VOCs on ozone formation. VOCs differ significantly in their impacts on O3 formation, and regulations that take this into account are potentially much more cost-effective than those that regulate all VOCs equally. In view of this, several VOC regulations implemented (or being considered) in California take reactivity into account. The current California regulations use the Maximum Incremental Reactivity (MIR) scale calculated using the SAPRC-99 chemical mechanism.
The SAPRC-99 mechanism includes representations of atmospheric reactions of almost 780 types of VOCs for reactivity assessment, and is widely used in other airshed modeling applications for research and regulatory applications. Although this represented the state of the art at the time it was developed, since then there has been continued progress in basic atmospheric chemistry, and new information has become available concerning the reactions and O3 impacts of many individual VOCs. In addition, the California Air Resources Board (CARB) is obligated to update the reactivity scales used in its regulations approximately every three years so they reflect the current state of the science. Since the mechanism was developed in 1999, updates to the mechanism and the reactivity scale are now due.
Another reason for updating the SAPRC mechanism is to make it more suitable for prediction of secondary particulate matter (PM), which is another air quality issue of concern. SAPRC-99, like most other mechanisms used in current airshed models, incorporates simplifications and approximations that may be appropriate for O3 modeling, but that restricts its capability to represent how secondary organic aerosol (SOA) formation is affected by chemical conditions. This needs to be addressed.
In view of this, the CARB funded us to provide an update to the SAPRC-99 mechanism used for modeling and VOC reactivity assessment, and to provide some needed improvements and enhancements. This report documents this mechanism, its evaluation against available environmental chamber data, and the updated reactivity scales that were developed.
Accomplishments
The major accomplishment of this project is the development of the SAPRC-07 chemical mechanism and its associated reactivity scales, which are documented in this report. Specific accomplishments are summarized as follows.
Base mechanism updated. All the reactions and rate constants in the mechanism have been reviewed based on results of current evaluations, and updated as needed. Most of the rate constant changes were relatively small, but a few errors were found and corrected and some potentially significant changes occurred. These have been assessed in the evaluations against chamber data.
Aromatics mechanisms reformulated. The mechanisms for the aromatic ring fragmentation reactions were reformulated to be more consistent with estimated explicit mechanisms, and to give predictions that are somewhat more consistent with available data. However, although an improvement over that used in SAPRC-99, the updated mechanism is still simplified in many respects, and is still not completely consistent with all of the available data
Chlorine chemistry added. A representation of chlorine chemistry has been added to the mechanism as an optional capability. In addition to improving the ability of airshed models to simulate air quality in regions impacted by chlorine emissions, the representation of chlorine chemistry has resulted in reduced uncertainties in reactivity estimates for chlorinated VOCs.
Capability for adaptation to SOA predictions improved. The method that the mechanism used to represent the reactions of peroxy radicals was reformulated so that effects of changes in NOx conditions on organic product formation can be more accurately represented. Because development of SOA mechanisms was beyond the scope of this project, the current mechanism does not fully take advantage of this capability, but it provides a framework upon which improved SOA mechanisms can be developed.
Mechanisms for many types of VOCs added or improved. The number of types of VOCs for which reactivity estimates have been made has been increased by over 20%, and the methods used to estimate mechanisms for a number of compounds were improved. This has involved enhancements of the capabilities of the mechanism estimation and generation system that is used to derive many of the mechanisms, and deriving estimated mechanisms for new classes of VOCs. A few errors found in the SAPRC-99 mechanism for some VOCs were corrected.
Updated mechanism evaluated against chamber experiments. The updated mechanism was comprehensively evaluated by comparing predictions with results of all environmental chamber experiments used for SAPRC-99 evaluation, plus the results of more recent UCR experiments, and experiments in other chambers. The results are summarized below.
Reactivity scales updated. The updated mechanism used to calculate MIR and other reactivity scales for all the ~1100 types of VOCs that are currently represented. Uncertainty classifications were also updated as part of this work. It is recommended that these be used to supercede the reactivity values distributed previously.
The mechanism developed in this project was implemented for the box model calculations used for reactivity scale calculations, and the data files used in this implementation can serve as the basis for implementing in more comprehensive airshed models such as CMAQ or CAMx. The data files are being made available at the project web site at http://www.cert.ucr.edu/~carter/SAPRC.
Results
Evaluation Results. In general, the performance of the updated mechanism in simulating the available environmental chamber data for individual compounds was comparable to SAPRC-99, though there were some differences. Some uncertain parameters for some compounds that were adjusted to fit the data for SAPRC-99 had to be re-adjusted for this mechanism. The updated aromatic mechanism simulated most of the experiments about as well as SAPRC-99, but some discrepancies observed with the previous version were reduced. Model performance was improved in simulating data for some compounds whose mechanisms were not changed, but biases were slightly increased with others. This is attributed to changes in the base mechanism, but the specific causes have not been determined.
One area of potential concern is that the mechanism update caused a slight increase in overall biases in model simulations of experiments with mixtures of VOCs. These changes are small compared to the 30% variability of the fits overall  -  which is less than the average biases  -  but because of the large number of such experiments (>1500 total) it may be statistically significant. The mechanism update also did not solve the problem, noted previously for other mechanisms, of underpredicting O3 formation and NO oxidation in ambient surrogate - NOx experiments carried out at relatively low NOx levels in the new UCR EPA chamber. These biases cannot be attributed to problems with individual VOCs, whose experiments are generally reasonably well simulated.
Preliminary assessment of impacts of updates on predictions of ambient ozone. Results of box model simulations of 1-day urban scenarios used for the reactivity scales indicate that the mechanism update caused changes in maximum ozone concentrations ranging from a ~10% decrease to a ~5% increase, with the predicted O3 decreasing by about 5% on average.  The largest increase appears at the lower ROG/NOx ratios, but other factors appear to be equally important. However, the scenarios used in the reactivity assessment calculations are highly simplified representations of ambient conditions, and comprehensive models are needed to fully assess the impacts of this update on ambient ozone and control strategy predictions.
Reactivity Scale Update.  For most VOCs, the reactivity scale update did significantly affect the reactivity values, with the average change in relative MIR values for the VOCs on the previous scales, being on the order of 10%. However, the MIR value changed by more than 5% for 56% of the VOCs given in the previous tabulation. The MIR change was less than 30% for approximately 93% of these VOCs, but 35 VOCs had changes greater than 35% and 5 VOCs had changes greater than a factor of 2. The latter consisted of halogenated compounds and one compound where a rate constant error was corrected.
Recommendations
Although the accomplishments of this project were significant, there were some objectives of this project that could not be met within the available time and resource, and there are other areas where additional work is recommended. These are summarized below.
Aromatics Mechanisms. Although the aromatics mechanisms developed in this work represent an improvement, problems remain and additional work is needed. Work was begun during this project in developing a more explicit aromatics mechanism, and new environmental chamber data useful for this effort was obtained. Work in this area needs to continue.
Chlorine and Halogen Mechanisms. Uncertainties remain in the mechanisms of halogenated compounds. Estimation methods need to be developed for reactions of halogen-containing radicals, and reactivities of halogenated oxidation products need to be assessed. Reactivity data are available for only a limited number of such compounds, and the available data indicates problems that need to be addressed. Mechanisms for bromine chemistry need to be developed.
Mechanism Generation System. An important component of the current mechanism is the mechanism generation and estimation system that is used to derive the mechanisms for most of the non-aromatic VOCs that are represented. However, because of time constraints we were unable to update most of the estimation methods incorporated in the system, except for the initial VOC rate constants and those related to chlorine chemistry. This needs to be done. The capability of the system needs to be enhanced to more reliably estimate mechanisms for additional classes of VOCs and intermediates, such as halogenated radicals, aromatics, and aromatic products. This will make it valuable for future mechanism improvements and reactivity scale updates. Finally, the system needs to be undergo peer review and be documented in the peer-reviewed literature.
Adaptation to SOA Predictions. Although the capability of the mechanism for improved SOA predictions has been enhanced, the potential of this capability has not been exploited. Recommendations in this regard include adding new species to the mechanism to represent low volatility products, implementing methods in the mechanism generation system to estimate volatility and incorporate them in the mechanisms so derived, and evaluating the predictions against SOA formation measured in chamber experiments.
Mechanism Performance Issues. The reason that the mechanism has biases in simulations of mixture experiments, while simulating single compound experiments reasonably well, needs to be investigated. This is necessary to assess the implications on ambient simulations of the biases in the simulations of the mixture experiments that were found, and for developing methods to reduce these biases. The reason for the relatively poor performance of the mechanisms in simulating the University of North Carolina outdoor chamber database also needs to be investigated. Although existing tools involved in sensitivity, uncertainty, and process analysis may be useful, new analysis methods probably need to be developed. This is an area where original research is needed.
Mechanism Evaluation Database. Although the database of chamber experiments useful for mechanism evaluation is very extensive and comprehensive in some respects, there are gaps and problems that need to be addressed. There are a number of classes of compounds where reactivity chamber data are needed to reduce mechanism uncertainties. Incremental reactivity experiments need to be developed that are more sensitive to reactions of organic oxidation products, which affect predicted reactivities in ambient scenarios much more than in chamber experiments. The current chamber dataset is not adequate for evaluating effects of temperature on mechanism performance. Experiments are needed for testing mechanisms for predictions of SOA formation, particularly under lower pollutant conditions more representative of ambient conditions than are most PM chamber data.
Mechanism Condensation. One of the objectives of this project that was not accomplished was to develop a condensed version of the mechanism to serve as an alternative to CB4/CB05 for modeling applications where computer speed is more important than chemical detail. This still needs to be done.
Reactivity Scenarios. The scenarios used for deriving the reactivity scales developed in this work are poorly documented, oversimplified, and do not represent current ambient conditions. Evaluations carried out by the Reactivity Research Working Group (RRWG) indicate that this methodology could be improved in a number of respects, particularly the scenarios and modeling methods. This update is way overdue, and may result in changes in relative reactivity values that are greater than those resulting from updates to the mechanism.
Next Mechanism Update. The CARB is still committed to updating its regulatory reactivity scale on a periodic basis. This will obviously need to include updated reactivity scenarios, as indicated above. Another problem is that the developer of the SAPRC mechanisms and the MIR scale is now semi-retired, and may be completely retired by the time the next update is needed, and it is unclear who will be carrying out this work in the future. Funding agencies need to show an interest in providing support for this type of mechanism development on a sufficiently consistent basis that it will attract younger researchers into this field. As it is now, support for mechanism development is relatively limited, and not of the type needed for the long-term commitment that this type or research requires.
Introduction
Airshed models are essential for the development of effective control strategies for reducing photochemical air pollution because they provide the only available scientific basis for making quantitative estimates of changes in air quality resulting from changes in emissions. The chemical mechanism is the portion of the model that represents the processes by which emitted primary pollutants, such as volatile organic compounds (VOCs) and oxides of nitrogen (NOx), interact in the gas phase to form secondary pollutants such as ozone (O3) and other oxidants. This is an important component of airshed models because if the mechanism is incorrect or incomplete in significant respects, then the model's predictions of secondary pollutant formation may also be incorrect, and its use might result in implementation of inappropriate or even counter-productive air pollution control strategies.
One airshed model application where the accuracy of the chemical mechanism is particularly important is the assessment or implementation of control strategies to encourage use of VOCs that have lower impacts on ozone or other secondary pollutant formation than VOCs that are currently emitted. Such strategies require a means to quantify the impacts, or "reactivities" of the VOCs with respect to O3 or other measures of air quality. There are several examples of control strategies where accurate O3 reactivity estimates are important. In the California Air Resources Board (CARB)'s "Low Emissions Vehicle/Clean Fuels" regulations, "reactivity adjustment factors" are used to place exhaust emissions standards for alternatively-fueled vehicles on an equal ozone impact basis as those for vehicles using conventional gasoline (CARB, 1993). More recently, the CARB implemented reactivity-based regulations for aerosol coatings (CARB, 2000), and is considering expanding such regulations to other types of stationary sources (e.g., CARB, 2006, 2007). In addition, the EPA has used O3 impacts of VOCs calculated for various environments among the factors they consider when evaluating proposals to exempt various compounds from controls as ozone precursors (Dimitriades, 1999).
 The MIR scale initially adopted in the CARB vehicle regulation was calculated using the SAPRC-90 chemical mechanism (Carter, 1990), but it has since been recalculated using an updated version of this mechanism, designated SAPRC-99 (Carter, 2000, 2003a). This mechanism has assigned or estimated mechanisms for over 500 types of VOCs. Although other state-of-the-art mechanisms are available for airshed model applications (e.g., Stockwell et al, 1997; Yarwood et al, 2005, Jenkin et al, 2003; Saunders et al, 2003), the SAPRC mechanisms were used for this purpose because they are the only mechanisms that represent a large number of VOCs that was comprehensively evaluated against environmental chamber data. However, although the SAPRC-99 mechanism represented the state of the art at the time it was developed, since then there has been continued progress in basic atmospheric chemistry, and new information has become available concerning the reactions and O3 impacts of many individual VOCs. In addition, the CARB is obligated to update the reactivity scales used in its regulations approximately every three years so they reflect the current state of the science. Since the last update was made in 2003 (Carter, 2003a), updates to the scale are now due.
In addition to calculation of reactivity scales for regulatory applications, a condensed version of the SAPRC-99 mechanism (Carter, 2000b) is widely used in comprehensive airshed models for prediction of effects of emissions on secondary pollutant formation in regional and urban atmospheres. Such regional models are required for many research and regulatory applications, so they should represent the current state of the science. In the United States, the two alternative mechanisms most generally used in comprehensive airshed models are the Carbon Bond 4 (CB4) (Gery et al, 1988) and the SAPRC-99 mechanism, and the two mechanisms have been found to give quite different predictions in some cases. The use of CB4 is preferred for many model applications because of it compact nature and because it is more computationally efficient, but SAPRC-99 is preferred in applications where chemical accuracy is a priority because it is more chemically detailed and more comprehensively evaluated against available environmental chamber data. CB4 has recently been updated to "Carbon Bond 05" (CB05) (Yarwood et al, 2005; Sarwar et al, 2008), and this is now being implemented in such models. Therefore, SAPRC-99 needs to be updated so it can continue to be a viable alternative to the Carbon Bond mechanisms for applications where chemical accuracy is a priority.
Another reason for updating the SAPRC mechanism is to make it more suitable for prediction of secondary particulate matter (PM). Fine particulate matter pollution an important issue because of the major health impacts it is believed to cause, and secondary PM is a major contribution to this problem. The formation of secondary PM is even more complex and incompletely understood than the formation of ozone. Models for prediction of secondary PM have appended aerosol models to gas-phase mechanisms such as SAPRC-99 and CB4, but the treatment is necessarily parameterized and approximate because of the condensation approaches incorporated in these mechanisms. These approaches may be appropriate for ozone modeling, but are not necessarily appropriate for prediction of secondary organic aerosol. For example, such mechanism use lumped species to represent oxidation products based on their gas-phase reactivity, but without considerations of volatility, which is the major factor in secondary organic aerosol (SOA) formation. In addition, SAPRC-99, like CB4 and CB05, uses a condensed method to represent the many peroxy + peroxy radical reactions that does not readily permit models to represent how organic product distributions may change under very low NOx conditions where these reactions are important. Recent data (e.g., Song et al, 2005) indicate that, at least for aromatics, secondary PM formation may be much more dependent on NOx conditions than represented in current models. The condensation methods used in current mechanisms need to be modified to be more suitably adapted for models for PM prediction
Based on these considerations, the California Air Resources Board funded us to develop and document an updated version of the SAPRC mechanism, and use it to derive updated MIR and other VOC ozone reactivity scales. The specific objectives of the project were as follows:
   * Update rate constants and reactions to current state of science.
   * Improve mechanisms for aromatics to incorporate new data and improve performance in simulating available chamber data.
   * Conduct environmental chamber experiments as appropriate to support this effort.
   * Add chlorine chemistry, to support modeling areas impacted by chlorine emissions and also calculating reactivities of chlorinated organics.
   * Update and enhance the mechanism generation system used to derive the mechanisms for most of the VOCs (Carter, 2000a). This includes updating the estimation methods and assignments as needed, and also enhancing the capabilities of the system, e.g., to support generating explicit mechanisms for aromatics and chlorine atom reactions;
   * Improve capability of the mechanism to be adapted to secondary PM models
   * Increase the number of VOC mechanisms to include more compounds present in emissions inventories or otherwise of interest to the CARB.
   * Develop new condensed mechanisms from the detailed version (including a highly condensed version as an alternative to CB4).
   * Make the mechanism available for implementation in airshed models
   * Calculate updated reactivity scales and update associated uncertainty classifications.
The major accomplishments of this project are the development of a completely updated version of SAPRC-99, which is designated SAPRC-07, and the update of the associated MIR and other reactivity scales, and their uncertainty classifications. The bulk of this report consists of the documentation of this mechanism and its evaluation, and the reactivity scales. The final section of this report includes a summary of the accomplishments of this project, the objectives that were not fully addressed, recommendations for future research.
Mechanism Documentation
The major components of the SAPRC mechanisms are the base mechanism used to represent the reactions of the inorganic reactants and the common organic oxidation products and radicals, the representation of the reactions of the individual VOCs, the lumped mechanism used for complex mixtures in airshed models, and the emissions assignments necessary for implementing the mechanism in airshed models. The general structure of this mechanism is essentially the same as that for SAPRC-99 mechanism documented by Carter (2000a), except that in this case chlorine chemistry is added and can be incorporated in the base and lumped mechanisms as an option, and is incorporated in representations and reactivity calculations for halogenated VOCs as needed. Each of these components, and their updates relative to SAPRC-99 are discussed in this section of this report.
Reference is made to the comprehensive documentation of the SAPRC-99 mechanism (Carter, 2000a) for documentation of the features of the mechanism that were retained in this version, so that document must be considered an integral part of the documentation of this mechanism. Appendices to the report include the tabulations or data plots that are too large to include in the main body of the text, and additional information, such as complete listings of the mechanisms of the over 700 types of VOCs that are explicitly represented, is available in electronic form at the SAPRC mechanism website at http://www.cert.ucr.edu/~carter/SAPRC.
Base Mechanism
The base mechanism is the portion of the mechanism that represents the reactions of the inorganic species, the common organic products, and the intermediate radicals leading to these products, including those formed from the initial reactions of the represented VOCs not in the base mechanism. Most of the VOCs that can be separately represented are not in the base mechanism, but can be added to the mechanism, either as explicit reactions for individual VOCs or as lumped model species whose parameters are derived from the mixture of detailed model species they represent, as needed in the model application. However, a few VOCs are represented explicitly and incorporated in the base mechanism, either because they are also common organic oxidation products that are represented explicitly, or are sufficiently important in emissions inventories and have sufficiently different mechanisms than most other VOCs that representation with lumped model species is inappropriate. These include formaldehyde, acetaldehyde, acetone, ethylene, isoprene, and (new to this version) benzene and acetylene. This portion of the mechanism is discussed in this section.
Listing of Standard Base Mechanism and Summary of Changes
The "Standard Base Mechanism" is the portion of the base mechanism excluding the optional additional reactions used to represent chlorine chemistry, and is the portion that is directly comparable to the SAPRC-99 base mechanism. Table A-1 in Appendix A gives the list of model species used in the standard base SAPRC-07 mechanism. These include the inorganic reactants, common reactive organic product species, chemical operators used to represent peroxy radical reactions (discussed below), and explicitly represented primary organics. Except as discussed below, the species used are the same as in the base SAPRC-99, though some model species have been renamed to be more compatible with some airshed model software systems. The changes, and additions, are as follows:
The set of lumped peroxy radical species and chemical operators was changed because as discussed in the "Project Summary" section, above, a different method was used to represent peroxy radical reactions. SAPRC-99 used three active chemical operators (RO2-R., RO2-N., and R2O2.) to represent the effects of peroxy radical reactions on NOx and organic nitrates. The updated mechanism uses two active chemical operators (RO2C and RO2XC) to represent the effects of peroxy radical reactions on NOx and 34 steady-state operators to represent the effects on formation of organic products (not counting the 3 used in the added chlorine mechanism).
       Two additional model species were added to represent the reactions of the higher hydroperoxides that were previously lumped with one generic hydroperoxide species in SAPRC-99. These higher hydroperoxides are expected to have different reactivities and SOA formation potentials than propyl hydroperoxide, which was used as the basis for the single generic hydroperoxide in SAPRC-99.
Although three active model species are still used to represent the reactions of the unsaturated aromatic ring fragmentation products, the types of compounds they represent have been changed. This is discussed in conjunction with the discussion of the revised aromatic mechanisms, below.
The reactions of formic acid, acetic acid, and the lumped higher organic acids were added to the mechanism. These species were in SAPRC-99 as inert tracers, but their reactions with OH radicals were added to this mechanism because they may be non-negligible loss processes in some regional modeling applications.
Acetylene and benzene are now represented explicitly in this mechanism. This is because both compounds are relatively important in emissions inventories, their reactivities are quite different than the other compounds with which they previously were lumped (monoalkylbenzenes and low reactivity alkanes, respectively), and are not well represented by other species used in the lumped mechanism. In addition, explicit simulations of benzene are of interest for toxics modeling, and acetylene can provide a useful tracer for vehicle emissions.
Although the objective of this update was to provide a mechanism with similar or, if appropriate, greater level of detail as SAPRC-99, a few SAPRC-99 model species were judged to be unnecessary and were removed from this version. These are as follows.
Phenol was removed from the mechanism because it is important only in the oxidation of benzene, and representing it by the lumped cresol species (CRES) did not significantly affect results of simulations of benzene, whose mechanism is very uncertain in any case.
The SAPRC-99 model species BZNO2-O, used in the mechanism for the reaction of NO3 with lumped nitrophenols, was removed. It was found that representing it with the phenoxy model species used in the cresol mechanism gave model simulations that were essentially the same, especially considering the large uncertainty in the nitrophenol and cresol mechanisms.
The SAPRC-99 model species CCO-OOH and RCO-OOH, used to represent various peroxy acids formed in the reactions of acyl peroxy radicals with HO2, were removed. These are represented by the reactions of the corresponding acid model species. It was judged that separate representation of these species was not necessary, though this could be changed in future versions of the mechanism if desired.
The reactions and rate parameters used in the base mechanism are given in Table A-2 in Appendix A, and Table A-3 gives the absorption cross sections and quantum yields used for the photolysis reactions listed in this table. Footnotes to Table A-2 indicate the source of the rate constants and mechanisms used. As indicated there, most of the updated rate constants are based on results of the IUPAC (2006) and NASA (2006) evaluations, though a number of other sources were also used as the basis for the updates. The major changes to the inorganic and common organic radical and product mechanisms are discussed further below. The changes to the mechanisms for the explicitly represented species are discussed later in conjunction with the mechanisms for the other individual VOCs.
Table 1 shows the changes in rate constants or atmospheric photolysis rates for this version of the base mechanism compared to SAPRC-99, for the cases where the rate constants or photolysis rates were changed by more than 5%. Reactions used to represent unsaturated aromatic ring fragmentation products, or reactions of chemical operators used to represent peroxy radical reactions are not shown because the representations are not comparable on a reaction-by-reaction basis; these are discussed below or later in this report. Specific changes of potential interest are as follows.
The most important single change to the base mechanism may be the ~19% increase in the rate constant for the OH + NO2 reaction, based on the results of the recent NASA (2006) evaluation. This affected the mechanism evaluation against the chamber experiments because it required re-deriving some chamber effects parameters and also some uncertain mechanistic parameters derived to fit chamber data. It is also expected to result in somewhat lower O3 predictions in ambient simulations, though this may be offset somewhat by the changes in the parameters adjusted to fit chamber data.
The ~7% increase in the NO2 photolysis rate under atmospheric conditions has no effect on NO2 photolysis rates used in chamber simulations because all photolysis rates are normalized to measured NO2 photolysis rates. However, rates of other photolysis reactions in chamber simulations will decrease accordingly.
There was a relatively large increase in the calculated atmospheric photolysis rates for the -dicarbonyl aromatic ring fragmentation products. However, this will not result in increases in calculated reactivities of aromatics because of changes that were made to the representation of the other reactive aromatic fragmentation products, discussed later in this report.
The ~30% increase in the calculated photolysis rates for methacrolein and the model species used to represent C5 aldehyde products formed from isoprene did not seem to have a significant effect on simulations of isoprene chamber experiments or calculations of the atmospheric reactivity of this important biogenic compound. Despite the relatively large decrease in the peroxy + HO2 rate constant, this reaction is still calculated to be the major loss process for peroxy radicals competing with reaction with NO under low NOx conditions.
 The photolysis of PAN was added to the mechanism at the request of Deborah Luecken of the EPA. Although not important in urban simulations, this reaction may be important under low temperature conditions such as occur at higher altitudes.
New information available concerning the reactions of nitrophenols indicates that the major atmospheric loss processes are photolysis and reactions with OH (see footnotes for these reactions on Table A-2). The speculative reaction of nitrophenols with NO3 was deleted, though it may still occur to some extent.
As discussed below, the mechanism for PROD2 is based on mechanism for various C5-C9 ketones. The updated mechanism gives better simulations of incremental reactivity environmental chamber experiments for higher ketones if lower photolysis rates for higher ketones are assumed (see discussion of mechanisms of individual VOCs and the evaluation against chamber experiments later in this report). The overall quantum yield for photolyses of C7+ ketones are set to a sufficiently low value that photolysis is unimportant. This is reflected in the average overall quantum yield used for PROD2.
As indicated above, several model species were removed and the base mechanism now explicitly represents the reactions of some additional compounds. These changes should have relatively small effects on most simulations except that atmospheric simulations will now simulate these species explicitly, rather than ignoring their reactions (as with the acids) or lumping them with other compounds. Note that this change does not, by itself, affect reactivity calculations for these compounds because the compound of interest (or the compound used to represent them if the "lumped molecule" approach is used) is always represented explicitly in reactivity calculations.

Table 1.	Reactions where the rate constants or photolysis rates changed by more than 5% or that were added or removed in the standard base mechanism relative to SAPRC-99.
                     Rate constant or photolysis rate [a]
Reaction [b]
                                   Notes [c]
                                   SAPRC-07
                                   SAPRC-99
                                    Change

                                    7.23e-1
                                    6.69e-1
                                      8%
NO2 + HV = NO + O3P
                                       1
                                   1.64e-12
                                   2.45e-12
                                     -33%
O3P + NO = NO2

                                   1.03e-11
                                   9.70e-12
                                      6%
O3P + NO2 = NO + O2

                                   3.24e-12
                                   1.79e-12
                                      81%
O3P + NO2 = NO3
                                       2
                                   2.02e-14
                                   1.87e-14
                                      8%
O3 + NO = NO2 + O2

                                   1.24e-12
                                   1.53e-12
                                     -19%
NO2 + NO3 = N2O5

                                    5.69e-2
                                    6.74e-2
                                     -16%
N2O5 = NO2 + NO3

                                   1.80e-39

N2O5 + H2O + H2O = #2 HNO3 + H2O
                                       3
                                    1.91e+0
                                    1.59e+0
                                      20%
NO3 + HV = NO + O2

                                   1.99e-10
                                   2.20e-10
                                     -10%
O1D + H2O = #2 OH

                                   3.28e-11
                                   2.87e-11
                                      14%
O1D + M = O3P + M

                                    1.14e-1
                                    1.27e-1
                                     -10%
HONO + HV = OH + NO

                                       -
                                    1.60e-2

HONO + HV = HO2 + NO2

                                   5.95e-12
                                   6.42e-12
                                      -7%
OH + HONO = H2O + NO2

                                   1.05e-11
                                   8.81e-12
                                      19%
OH + NO2 = HNO3
                                       1
                                   2.28e-13
                                   2.08e-13
                                      10%
OH + CO = HO2 + CO2

                                   7.41e-14
                                   6.78e-14
                                      9%
OH + O3 = HO2 + O2

                                   8.85e-12
                                   8.36e-12
                                      6%
HO2 + NO = OH + NO2

                                   1.12e-12
                                   1.36e-12
                                     -18%
HO2 + NO2 = HNO4

                                    1.07e-1
                                    9.61e-2
                                      11%
HNO4 = HO2 + NO2

                                    5.42e-4
                                    4.69e-4
                                      16%
HNO4 + HV = Products

                                   4.61e-12
                                   4.98e-12
                                      -7%
HNO4 + OH = H2O + NO2 + O2

                                   1.69e-15
                                   1.89e-15
                                     -11%
HO2 + O3 = OH + #2 O2

                                   1.80e-12
                                   1.70e-12
                                      6%
HO2H + OH = HO2 + H2O

                                   4.65e-12
                                   5.12e-12
                                      -9%
MEO2 + HO2 = COOH + O2

                                   4.50e-13

MEO2 + HO2 = HCHO + O2 + H2O

                                   2.16e-13
                                   2.61e-13
                                     -17%
MEO2 + MEO2 = MEOH + HCHO + O2

                                   1.31e-13
                                   1.08e-13
                                      21%
MEO2 + MEO2 = #2 {HCHO + HO2}

                                   7.63e-12
                                   1.45e-11
                                     -47%
Peroxy + HO2 = Products
                                       4
                                   9.37e-12
                                   1.04e-11
                                     -10%
MECO3 + NO2 = PAN

                                    6.27e-4
                                    7.04e-4
                                     -11%
PAN = MECO3 + NO2

                                    6.12e-5

PAN + HV = Products
                                       1
                                   1.97e-11
                                   2.12e-11
                                      -7%
MECO3 + NO = MEO2 + CO2 + NO2

                                   2.30e-12
                                   4.00e-12
                                     -43%
MECO3 + NO3 = Products
                                      2,5
                                   1.06e-11
                                   9.53e-12
                                      11%
MECO3 + MEO2 = Products
                                       5
                                   1.56e-11
                                   7.50e-12
                                     108%
MECO3 + Peroxy = Products
                                     2,4,5
                                    5.48e-4
                                    5.90e-4
                                      -7%
PAN2 = RCO3 + NO2

                                   2.08e-11
                                   2.78e-11
                                     -25%
RCO3 + NO = Products

                                   7.63e-12
                                   1.45e-11
                                     -47%
BZO + HO2 = CRES

                                   3.80e-11

BZNO2-O. + NO2 = inert products
                                       6

                                   1.49e-11

BZNO2-O. + HO2 = NPHE
                                       6

                                    1.00e-3

BZNO2-O. = NPHE
                                       6
                                    2.76e-3
                                    2.32e-3
                                      19%
HCHO + HV = #2 HO2 + CO

                                   8.47e-12
                                   9.19e-12
                                      -8%
HCHO + OH = HO2 + CO + H2O

                                   1.49e-11
                                   1.57e-11
                                      -6%
CCHO + OH = MECO3 + H2O

                                   6.74e-15
                                   3.82e-15
                                      76%
RCHO + NO3 = HNO3 + RCO3

                                    3.24e-5
                                    4.16e-5
                                     -22%
ACET + HV = Products

                                    1.69e-4
                                    1.42e-4
                                      19%
MEK + HV = Products

                                   6.78e-12
                                   1.10e-11
                                     -38%
ROOH + OH = Products

                                    3.18e-3
                                    1.08e-3
                                     193%
GLY + HV = HCHO + CO

                                    1.56e-2
                                    1.10e-2
                                      42%
MGLY + HV = HO2 + CO + MECO3

                                    2.67e-2
                                    1.90e-2
                                      41%
BACL + HV = #2 MECO3

                                   2.63e-11

PHEN + OH = products
                                       6

                                   3.78e-12

PHEN + NO3 = products
                                       6
                                   3.50e-12
                                       -

NPHE + OH = BZO + XN
                                       1

                                   3.78e-12

NPHE + NO3 = products
                                       1
                                    1.08e-3

NPHE + HV = HONO + #6 XC
                                       1
                                    1.08e-2

NPHE + HV = #6 XC + XN
                                       1
                                   1.20e-11
                                   1.29e-11
                                      -7%
BALD + OH = BZCO3

                                   2.84e-11
                                   3.35e-11
                                     -15%
MACR + OH = Products

                                   1.28e-18
                                   1.19e-18
                                      8%
MACR + O3 = Products

                                   3.54e-15
                                   4.76e-15
                                     -26%
MACR + NO3 = Products

                                    1.97e-4
                                    1.36e-4
                                      45%
MACR + HV = Products

                                   1.99e-11
                                   1.87e-11
                                      6%
MVK + OH = Products

                                   5.36e-18
                                   4.74e-18
                                      13%
MVK + O3 = Products

                                    7.50e-5
                                    6.97e-5
                                      8%
MVK + HV = Products

                                    1.97e-4
                                    1.36e-4
                                      45%
IPRD + HV = Products

                                    4.69e-6
                                    1.90e-5
                                     -75%
PROD2 + HV = Products
                                       1
                                   7.20e-12
                                   7.80e-12
                                      -8%
RNO3 + OH = Products

                                   4.50e-13

HCOOH + OH = HO2 + CO2
                                       7
                                   7.26e-13

CCOOH + OH =products
                                       7
                                   1.20e-12

RCOOH + OH = products
                                       7
                                   1.64e-11
                                       -

R6OOH + OH = products
                                       7
                                    3.94e-4
                                       -

R6OOH + HV = products
                                       7
                                   1.08e-10
                                       -

RAOOH + OH = products
                                       7
                                    3.94e-4

RAOOH + HV = products
                                       7
                                   7.56e-13

ACETYLEN + OH = products
                                       7
                                   1.16e-20

ACETYLEN + O3 = products
                                       7
                                   1.22e-12

BENZENE + OH = products
                                       7
 [a]	Rate constants are for T=300[o]K and are in molecule, cm[3] sec-1 units. Photolysis rates are calculated for direct overhead sunlight based on the actinic fluxes used in the reactivity scale calculations (Carter, 1994a,b).
[b]	Reactions used to represent unsaturated aromatic ring fragmentation products, and reactions of chemical operators used to represent peroxy radical reactions, are not directly comparable in the two mechanisms. They are discussed separately later in this report.
[c]	Notes concerning changes for reaction. See also text and footnotes to Table A-2
   1	See text for a discussion of this change.
   2	This is a relatively minor process under most conditions of interest, so the change should not have a significant effect on most model predictions.
   3	This reaction is added to represent the expected humidity dependence of this process.
   4	"Peroxy" refers to the various model species used to represent various types of peroxy radicals or peroxy radical operators.
   5	The same rate constants are used for all the higher acyl peroxy radical model species.
   6	This model species was deleted. See text.
   7	The reactions of these compounds were added to the base mechanism. See text.

Note that the discussion given above of expected effects of these changes on model simulations is based on expectations that have not, in most cases, been verified by actual model sensitivity calculations. Such sensitivity calculations would be useful to assess the effects of various changes, and help focus on areas where basic research may be useful.
Representation of Peroxy Radical Operators
Because of the large number of peroxy radicals that are involved even in condensed atmospheric chemistry mechanisms, it is generally not practical to represent the many possible peroxy + peroxy reactions explicitly, especially considering that under most conditions, especially conditions favorable for O3 formation, most of these reactions are relatively unimportant. Even highly explicit mechanisms such as the MCM (Jenkin et al, 2003; Saunders et al, 2003) use an approximate method to represent the many peroxy + peroxy cross reactions involving a chemical operator representing the total peroxy radical concentration. The RADM-2 and RACM mechanisms (Stockwell et al, 1990, 1997) have separate peroxy radical species for each VOC or VOC product model species whose reactions form peroxy radicals, but neglect peroxy + peroxy reactions except for those involving HO2 and methyl peroxy radicals. Because of the large number of reactions and model species involved even with the more approximate RADM-2 representation, SAPRC-99 represents methyl peroxy radicals explicitly, but uses a limited number of "chemical operators" to represent effects of peroxy radical reactions on NOx and radicals, and represents the organic products formed when higher peroxy radicals react with other peroxy radicals by those formed when they react with NOx. The CB4/05 mechanisms use a similar, though somewhat more condensed approach.
Use of these condensed representations of peroxy radical reactions has been shown to have relatively little effects on predictions of O3 formation and overall gas-phase reactivity because they involve no approximation when the major fate of peroxy radicals is reaction with NO, as is the case when O3 formation occurs, and because they give reasonably good representations of how NO to NO2 conversions, organic nitrate formation, and regeneration of radicals change when NOx levels are reduced to the point where the competing peroxy + peroxy become non-negligible. However the representations incorporated in SAPRC-99 and CB4/05 do not represent the changes in organic oxidation products that occur when these peroxy + peroxy reactions become non-negligible, since they use the set of products formed in the peroxy + NO reaction as the surrogate for the generally different products formed in the competing reactions. The inability of this representation to represent the formation of hydroperoxides formed in the peroxy + HO2 reaction is of particular concern, because these are predicted to be the major competing products formed under low NOx conditions (e.g., see Carter, 2004), and as discussed above are believed to be important precursors to secondary organic aerosol (SOA) formation.
Therefore, the peroxy radical representation used in SAPRC-99 is not satisfactory for use of the mechanism for prediction of SOA formation, and for this reason was changed as part of this update. The peroxy radical operator method implemented in SAPRC-99 is described in the SAPRC-99 documentation (Carter, 2000a). Briefly, it involves the operator RO2-R. to represent the reactions converting NO to NO2 forming HO2, RO2-N. to represent reactions with NO forming organic nitrates, and R2O2. to represent extra NO to NO2 conversions involved in multi-step peroxy radical reactions. Since it uses the set of products formed in the NO reaction to represent the products in the competing peroxy + peroxy reactions, it uses no additional model species for this purpose except for the inclusion of a generic "ROOH" species to represent reactions at the hydroperoxide group formed in the peroxy + HO2 reaction (with the set of products formed in the NO reaction also being formed in the HO2 reaction.)
The SAPRC-07 representation is similar in that operators are also used to represent the effects of peroxy radicals on NOx, with "RO2C" representing NO to NO2 conversions and "RO2XC" representing NO consumption that occurs in conjunction with nitrate formation, but uses separate chemical operator model species to represent the formation of radicals and oxidation products, and how they depend on which peroxy radical reactions are occurring. Three sets of chemical operators are used for this purpose: "xPROD" species are used to represent the formation of the alkoxy radical products resulting when the peroxy radical react with NO, NO3, and (in part) RO2; "yPROD" species are used to represent the formation of hydroperoxides formed when peroxy radicals react with HO2 or H-shift disproportion products formed when peroxy radicals react with acyl peroxy radicals or (in part) with RO2; and "zRNO3" species are used to represent the formation of organic nitrates when peroxy radicals react with NO or the alkoxy radical formed (in part) in RO2 + RO2 reactions. The products formed in the reactions of these operators are summarized on Table 2. Table A-2 and footnotes to Table A-2 indicate how they are implemented in the model.

Table 2.	Products formed in the reactions of the chemical operators used to represent peroxy radical reactions in the SAPRC-07 mechanism.
Reaction
RO2C
RO2XC
xPROD
yROOH
zRNO3
NO
- NO + NO2 [a]
- NO [a]
PROD

RNO3
NO3
- NO3 + NO2 [a]

PROD

PROD2 + HO2 [b]
HO2

ROOH

RO2·
 RO·.

PROD

PROD2 + HO2 [b]

H-Shift

MEK or PROD2 [c]

RC(O)O2·

PROD

PROD2 + HO2 [b]
[a]	"-NO" or "-NO3" means that the reaction consumes these reactions. If this is not shown, it means that the reaction should not consume any reactant other than the peroxy radical operator. For example, the reaction of RO2C with HO2 would be simulated as "RO2C + HO2  HO2".
[b]	This represents products formed from alkoxy radicals formed in the absence of NOx that are not represented by xPROD model species because they are not used for the portions of the reactions that form organic nitrates in the presence of NOx. The present version of the mechanism has only a single zRNO3 species, and uses PROD2 for this purpose. HO2 is used to represent the radicals formed.
[c]	The model species used to represents the H-shift disproportion products depends on the size of molecules being represented by yROOH. For the operators forming the ROOH model species used for 4 or fewer carbons, MEK is used for these processes. For operators forming larger hydroperoxides (R6OOH or RAOOH), PROD2 is used.
Note that in this representation the total yield of RO2C or RO2C + RO2XC may be greater than unity in multi-step processes involving formation of secondary peroxy radicals that cause additional NO to NO2 conversions or nitrate formation. However, radical conservation requires that the sum of RO2XC and all xPROD operators where PROD is a radical be equal to 1 in VOC + radical (e.g., VOC + OH) reactions, and the sum of all yROOH species must be equal to the total yield of all peroxy radicals formed in the initial reaction.
The main approximation involved in this representation concerns the treatment of the cases of multi-step mechanisms where the peroxy radicals undergo reactions to form other peroxy radicals, which can react with NO or other peroxy radicals, etc. This would result in different (generally intermediate) branching ratios in the competitions of the peroxy radical reactions in terms of the overall products. In this representation it is assumed that the overall branching ratios, e.g., the extent to which hydroperoxides are formed from peroxy + HO2 reactions, are the same in multi-step as in single step mechanism. Removing this approximation would require a much more complex mechanism, with more operators or model species, than is appropriate given its level of importance. Note that this approximation is incorporated in all mechanisms that use the more approximate peroxy radical representations, such as SAPRC-99 or CB4/05.
Although this representation can otherwise potentially give the same predictions as fully explicit mechanisms, it is not particularly intuitive from a chemical perspective. It involves using separate model species for each type of product that is formed, and its "reactions" do not have a straightforward correspondence to explicit reactions. An alternative method, that is more straightforward to understand in terms of actual chemical processes (and gives the same predictions), is to use separate model species for each group of peroxy radicals formed from the reactions of each of the various VOCs, and represent the competing overall reactions of these lumped groups of peroxy radicals with NO, NO3, HO2, RCO3, and other peroxy radicals. Chemical operators would still be needed to represent NO to NO2 conversions in multi-step processes and to determine a total RO2 concentration for calculating peroxy + peroxy rates, but the result would be closer to the actual chemical processes that occur. This is the approach used in the RADM2/RACM mechanisms (Stockwell et al, 1990, 1997), with the problem of representing peroxy + peroxy reactions dealt with by ignoring all but reactions with methyl peroxy; which is not a bad approximation (Carter and Lurmann, 1990).
Although we considered use of the RADM2 approach for this mechanism, we had to abandon it because it was incompatible with the "lumped parameter" approach incorporated in the SAPRC mechanisms for flexible representation of the hundreds of individual VOCs or deriving parameters for lumped model species based on the specific compounds they represent. This representation involves having a numerical parameter represent each of the product or radical model species involved in their overall reactions, including NO to NO2 conversions and NO consumptions. This is not readily adaptable to the RADM2 peroxy radical representation because the overall products are not associated directly with the reacting VOCs but with their peroxy radicals, and the yield of peroxy radicals from VOCs can vary. On the other hand, it is readily adaptable to the peroxy radical representation incorporated in SAPRC-07, where each overall product is still directly related to each reacting VOC.
Base Chlorine Mechanism
The model species added to the base mechanism to represent the atmospheric reactions of chlorine species are listed on Table A-4 in Appendix A. These include 8 model species to represent active inorganic reactants and radicals, 2 to represent chlorine-containing oxidation products, 3 steady state chemical operators to represent formation of chlorine species in peroxy radical reactions, and 2 chlorine-containing inert tracer species. Table A-5 gives the reactions and rate parameters used in the base chlorine mechanism, and Table A-6 gives the absorption cross sections and quantum yields for the photolysis reactions listed in this table. Footnotes to Table A-5 and indicate the source of the rate constants and mechanisms used. As indicated there, most of the rate constants are based on results of the IUPAC (2006) and NASA (2006) evaluations, though several other sources were also used in some cases.
Except for the updated rate constants and photolysis data, the base inorganic chlorine mechanism is very similar to that developed previously by Carter et al (1996a, 1997b). The major difference is the addition of separate model species to represent the reactions of chlorinated aldehydes and ketones that may be formed in the reactions of chlorinated VOCs. Carter and Malkina (2007) found that the reactivities of the 1,3-dichloropropenes are significantly underpredicted if the standard aldehyde model species are used to represent the reactions of the chloroacetaldehyde predicted to be formed, but satisfactory simulations are obtained if a separate model species is used, using the NASA (2006)-recommended absorption cross sections and quantum yields for this compound. This is because this compound is calculated to photolyze ~15 times faster than acetaldehyde and ~4.5 times faster than the lumped higher aldehyde model species used in the standard base mechanism. This model species (CLCCHO) is also used to represent the reactions of other -chloroaldehydes, which are assumed to be similarly photoreactive. Because chloroacetone is calculated to photolyze ~7 times faster than MEK and even more for higher ketones, a chloroacetone (CLACET) model species is added to represent the reactions of -chloroketones. However, -chloro- aldehydes and ketones and other chlorinated aldehydes are still represented by the generic higher aldehyde or ketone species (RCHO, MEK, or PROD2) in the standard mechanism.
Chlorine atoms react rapidly with most reactive VOCs and any complete chlorine mechanism must include a representation of their reactions, at least for the VOCs present in the chamber experiments with chlorine-containing species, and for the explicit and lumped VOC species used in atmospheric models. The mechanisms for these reactions are discussed in conjunction with the mechanisms for individual VOCs and lumped mechanisms for airshed models, below.
Representation of Organic Products
The set of model species to represent the reactions of the organic oxidation products are given in Table A-1 in Appendix A. Some of these are compounds represented explicitly (e.g. formaldehyde) and some are lumped species whose mechanisms are derived based on that for a compound or group of compounds chosen as representative. In the latter case, the compound or compounds used to derive the mechanisms for the lumped species are given in Table A-1 and in applicable footnotes to Table A-2. Except for the additional species added to represent higher hydroperoxides and the removal of phenol and the organic peroxyacids (now lumped with cresols or organic acids, respectively), this mechanism uses essentially the same set of model species to represent the organic products, and the same set of compounds to derive their mechanisms in the case of lumped species, as does SAPRC-99. Therefore the discussion of these species, and the choice of representative compounds used to derive the mechanisms of the lumped species, given in the SAPRC-99 documentation (Carter, 2000a) are still applicable. However, the rate constants and in some cases the products used for the reactions of these species or representative compounds were updated, and the reactions of these species with chlorine atoms were added and incorporated in the base chlorine mechanism discussed above.
The optional added chlorine mechanism given in Table A-5 includes the reactions of Cl atoms with most of the model species used reactive organic products in the base mechanism. Although reactions of Cl with some organic product species have been omitted, the processes represented should be sufficient for most ambient and chamber simulations, where the major sinks for the chlorine atoms or the product species are other reactions. As indicated above the added chlorine mechanism also includes model species for chlorinated aldehydes and ketones that may be formed from the reactions of chlorinated VOCs. The mechanism also includes the reactions of chlorine atoms with these compounds.
Mechanisms for Individual VOCs
An important distinguishing feature of the current SAPRC mechanisms is their ability to optionally include separate representations of the atmospheric reactions of many hundreds of different types of VOCs. Generally these are not all included in the airshed model at the same time, but selected individual compounds can be represented explicitly for the purpose of calculating reactivity scales or for toxics modeling, or mechanisms for groups of compounds can be used for deriving parameters for lumped model species used for representing complex mixtures such as ambient emissions. The current mechanism has separate explicit representations for 734 types of VOCs, making it the most comprehensive in terms of types of VOCs than any current gas-phase atmospheric mechanism. In addition, for deriving reactivity scales or representations of complex mixtures, 304 additional types of VOCs are represented using the "lumped molecule approach", where their impacts are estimated by assuming they are the same, on a per molecule basis, of an explicitly represented VOC. This yields a total of 1038 types of VOCs that are represented by a "detailed model species" in the mechanism.
A listing of all of the detailed model species in the SAPRC-07 in the mechanism, their molecular weights and general representation method (explicit or lumped molecule) is given in Table B-1 in Appendix B. That table also gives codes for availability of mechanism evaluation data, calculated ozone reactivity values, and reactivity uncertainty classifications. Table B-2 and Table B-4 through Table B-6 give the reactions and rate parameters for all the explicitly represented VOCs, Table B-3 gives the mechanisms for those VOCs where "adjusted product" mechanisms (discussed below) were used in the reactivity calculations, Table B-7 gives the mechanisms for the reactions of chlorine with the VOCs that were used in developing the fixed parameter lumped mechanism discussed in the following section. In addition to giving the OH radical rate constants and references for all VOCs that are explicitly represented in the mechanism, Table B-4 gives codes indicating the general type of mechanism used for each, and gives the structures used for compounds whose mechanisms were derived using the mechanism generation system. Footnotes to Table B-2 and Table B-4 through Table B-7 document the sources of the rate constants used, and Table B-8 gives the absorption cross sections and quantum yields for those that are photoreactive and have different absorption cross sections and quantum yields than those used in the base mechanism. Table B-9 gives the lumped molecule assignments used for the types of VOCs that are represented using this approximation. The derivations of these mechanisms are discussed further in the remainder of this section.
Mechanism Generation and Estimation System
A major tool used in the derivation of the SAPRC-99 mechanisms for individual VOCs was the automated mechanism generation and estimation software system that is described in the SAPRC-99 documentation (Carter, 2000a). Given the structure of the molecule, available information concerning applicable rate constants and branching ratios, and various estimation methods, this system can derive fully explicit mechanisms for the atmospheric reactions of many classes of VOCs in the presence of NOx, and then use these mechanisms to derive the appropriate representations of the VOCs in terms of the model species in the base mechanism. This was used to derive the SAPRC-99 mechanisms for the acyclic and monocyclic alkanes, acyclic and monocyclic monoalkenes, many classes of oxygenates including alcohols, ethers, glycols, esters, aldehydes, ketones, glycol ethers, carbonates, etc, and the organic nitrates. Although many of the estimated rate constants and rate constant ratios are highly uncertain, this procedure provided a consistent basis for deriving "best estimate" mechanisms for chemical systems which are too complex to be examined in detail in a reasonable amount of time. The system allows for assigning or adjusting rate constants or branching ratios in cases where data are available, or where adjustments are necessary for model simulations to fit chamber data. Therefore, it could be used for deriving fully detailed mechanisms for VOCs that fully incorporate whatever relevant data are available, relying on various estimation methods only when information is not otherwise available. The program also outputs documentation for the generated mechanism, indicating the source of the estimates or assumptions or explicit assignments that were used.
This SAPRC-99 mechanism generation system, with some enhancements and updates as described below, was also used for deriving mechanisms for many of the VOCs in the SAPRC-07 mechanism. Since this system is comprehensively documented by Carter (2000a) and most of the estimation methods were not updated as part of this project, in this section we will restrict our discussion to the changes, updates, and enhancements that were made as part of this project and utilized when deriving the updated mechanisms for the individual VOCs. These are described below
Enhanced Capabilities for VOCs and Reactions
During the period of this project, the capabilities of the mechanism generation system were enhanced to support generating mechanisms for VOCs and radicals that could not previously be processed. These are briefly summarized below. These enhancements, and the associated estimates and branching ratio assignments, will be described in more detail in a future report when the current system is more comprehensively documented.
Reactions of VOCs with chlorine atoms can now be generated. This is discussed in the following subsection.
The ability to represent compounds with more than one ring has been added. This required revising the way structures were identified. Because of this, mechanisms for terpenes can now be generated. This is also useful in generating estimated mechanisms for aromatics, discussed below.
The ability to represent compounds with more than one double bond has been added. This also required revising the way structures were identified. However, estimation methods for reactions of compounds with conjugated or cumulated double bonds were not developed, and automatically generated rate constants and branching ratios for OH and NO3 reactions of such compounds are not necessarily chemically reasonable, and reactions of O3 with compounds with more than one double bond are still not supported. Explicit assignments have to be made for such reactions, as was the case with, for example, isoprene and 1,3-butadiene.
The ability to estimate reactions of unsaturated radicals has been added. The system recognizes allylic resonance and the probability of reaction at various positions of alkyl radicals with allylic resonance is estimated by assuming that the probability of reaction at a particular structure is proportional to exp(Hf/RT), where Hf is the estimated heat of formation of the structure, R is the gas constant, and T is the temperature. If the heat of formation cannot be estimated, equal probability for reaction at all positions (or the positions with the most conjugated double bonds if there are more than one double bond) is assumed.
Reactions involving H-atom abstractions are now considered when generating mechanisms for reactions of OH, NO3, and Cl with alkenes. These previously had to be assumed to be negligible because of lack of support for generating mechanisms for unsaturated radicals. These are still neglected if they are estimated to occur less than 1% of the time, as is generally the case in OH and NO3 reactions. However, abstraction reactions are generally not negligible in chlorine + alkene reactions, as discussed below.
The ability to generate mechanisms for amines (compounds with >N-, -NH-, or -NH2 groups) was added. The new amine mechanisms, and their evaluation against chamber data, are documented by Carter (2008a). The available kinetic data indicates that reaction with OH occurs by H-atom abstraction from the amino group and that the rate of H-atom abstraction from the carbon next to the amino group is significantly enhanced, and group-additivity estimates for rates of reaction at these positions were derived. The data of Aschmann and Atkinson (1998) on N-methyl pyrrolidinone (NMP) indicate that reaction of amines with NO3 radicals may also occur to a non-negligible extent, and estimates (albeit highly uncertain) of NO3 reactions for amines in general were derived. Methods for estimating rate constants and mechanisms for reactions of amines with O3 were also developed. A more complete discussion and documentation of the amine mechanisms are given by Carter (2008a).
The ability to represent alkynes was added. The mechanisms for OH reactions were assumed to be analogous to that derived for acetylene (Carter et al, 1997c, Carter, 2000a).
The ability to represent alkylbenzenes and their initial reactions has also been added. However, our ability to estimate these reactions has not evolved to the point where this can be used to generate reliable mechanisms useful for modeling, and this has been used primarily as a research tool in working towards eventually developing more explicit mechanisms for aromatics. This work, which is still underway, was useful in deriving the updated condensed mechanisms for aromatics discussed later in this section.
The enhanced capabilities required group additivity estimates deriving rate constants for additional types of abstraction reactions that were not previously supported, such as abstractions from compounds with double bonds, abstractions from aromatics, and abstraction reactions by NO3 radicals and Cl atoms. As discussed previously (e.g., see Kwok and Atkinson, 1995 or Carter, 2000a), rate constants for abstraction processes can be estimated using

where "group" refers to various parts of the molecule for which group additivity assignments are made. Group additivity parameters used for estimating abstraction reactions are given in Table 3. The assignments for abstractions by OH are the same as used previously (Carter, 2000a), with assignments for reactions of aromatics or amines added as indicated by footnotes to the table. The assignments for the NO3 reactions are primarily from Atkinson (1991). The assignments for abstractions by chlorine atoms are discussed in the following section.
Support for Chlorine Reactions
As indicated above, in order to support the extension of the mechanism to represent chlorine chemistry, the mechanism generation system was enhanced to generate mechanisms for chlorine + VOC reactions. Since the reactions of Cl radicals with VOCs are similar to those of OH. This general procedure can be readily adapted to generate mechanisms for Cl + VOC reactions. This requires (1) making group-additivity estimates for all the possible initial reactions of Cl with VOCs, (2) generating mechanisms for the Cl-containing radicals that can be formed when Cl adds to double bonds, and (3) generating mechanisms for the radicals formed. The general procedures are discussed in detail by Carter (2000a). The specific adaptations to chlorine atom reactions are discussed below.
Chlorine can react with VOCs either by abstracting a hydrogen atom to form HCl and the corresponding alkyl radical, or by adding to a double bond. For abstraction reactions, the rate constants can be estimated using group additivity methods, with the rate constant being determined by the sum of the abstraction rate constant assigned to the group, multiplied by substituent correction factors for each non-hydrogen substituent on the group, summed over all groups with hydrogen atoms (e.g., see Kwok and Atkinson, 1995 or Carter, 2000a). Note that the correction factor for methyl substitution is arbitrarily set at unity, with the factors for the other substituents being determined based on differences in rate constants at groups that are only methyl substituted.

Table 3.	Group additivity parameters used for estimating rate constants for abstraction reactions by OH and NO3 radicals and by Cl atoms.
            Table 3a. Abstraction rate constants assigned to groups
Group
                                    OH [a]

                                      NO3

                                      Cl

                                    k(300)
                                       A
                                       B
                                       D
                                   Notes [b]

                                    k(300)
                                     Notes

                                    k(300)
                                     Notes
CH3
                                   1.39e-13
                                   4.49e-18
                                       2
                                      320
                                      1,2

                                   7.00e-19
                                       9

                                   3.43e-11
                                      13
CH2
                                   9.41e-13
                                   4.50e-18
                                       2
                                     -253
                                      1,2

                                   1.50e-17
                                       9

                                   6.77e-11
                                      13
CH
                                   1.94e-12
                                   2.12e-18
                                       2
                                     -696
                                      1,2

                                   8.20e-17
                                       9

                                   4.46e-11
                                      13
OH
                                   1.42e-13
                                   2.10e-18
                                       2
                                      85
                                      1,2

                                       0
                                      10

                                       0
                                      10
CHO
                                   1.56e-11
                                   5.55e-12
                                       0
                                     -311
                                      1,3

                                 2.84e-15 [c]
                                      11

                                   6.64e-11
                                      13
HCO(O)
                                       0

                                      1,4

                                       0
                                      10

                                       0
                                      12
OH(O)
                                   5.30e-12
                                   7.82e-17
                                       2
                                      85
                                       5

                                       0
                                      10

                                       0
                                      10
CH3(Bz) [d]
                                   4.92e-13

                                       6

                                   7.00e-19
                                      12

                                   3.43e-11
                                      12
CH2(Bz) [d]
                                   1.88e-12

                                       7

                                   1.50e-17
                                      12

                                   6.77e-11
                                      12
CH(Bz) [d]
                                   1.33e-12

                                       8

                                   8.20e-17
                                      12

                                   4.46e-11
                                      12
CH3(NHx) [e]
                                   1.67e-11

                                      22

                                   ~6.2e-14
                                      24

                                      [f]

CH2(NHx) [e]
                                   1.84e-11

                                      22

                                   ~4.4e-14
                                      24

                                      [f]

CH(NHx) [e]
                                   2.01e-11

                                      23

                                   ~9.9e-14
                                      24

                                      [f]

NH2
                                   1.41e-11

                                      22

                                   ~5.2e-14
                                      24

                                      [f]

NH
                                   3.17e-11

                                      22

                                   ~5.7e-14
                                      24

                                      [f]

[a]	Temperature dependences for OH rate constants given by k(T) = A B[T] exp(-D/T), where T is in [o]K.
[b]	Notes for derivations of the group rate constants and substituent correction factors are given below (Note [a] with Table 3b).
[c]	The temperature dependence of this group rate constant is given by 1.40e-12 x exp(-1860/T).
[d]	"Bz" refers to any aromatic carbon.
[e]	Applicable for x=0,1,2.
[f]	Group rate constants for reactions of chlorine with amines were not derived

         Table 3b. Group correction factors for abstraction reactions
Group
                         Substituent Correction Factor

                                      OH

                                      NO3

                                      Cl

                                    Factor
                                   Note [a]

                                    Factor
                                     Note

                                    Factor
                                     Note

CH3
                                       1
                                    1,2,14

                                       1
                                      17

                                       0
                                      14
CHx (x<3)
                                     1.23
                                      1,2

                                     1.34
                                      11

                                     0.95
                                      13
CHx(CO)
                                     3.90
                                      1,2

                                     1.34
                                      17

                                     0.95
                                      17
CHx(CO-O)
                                     1.23
                                      1,2

                                     1.34
                                      17

                                     0.95
                                      17
CHx(Cl)
                                     0.36
                                      1,2

                                     1.34
                                      17

                                     0.19
                                      19
CHx(Br)
                                     0.46
                                      1,2

                                     1.34
                                      17

                                     0.95
                                      17
CHx(F)
                                     0.61
                                      1,2

                                     1.34
                                      17

                                     0.95
                                      17
OH
                                     3.50
                                      1,2

                                       0
                                      18

                                     1.07
                                      13
CHO
                                     0.75
                                      1,2

                                     0.18
                                      11

                                     0.40
                                      13
CO
                                     0.75
                                      1,2

                                     0.89
                                      11

                                     0.04
                                      13
CO(O)
                                     0.31
                                      15

                                       0
                                      12

                                     0.04
                                      17
CO(OH)
                                     0.74
                                      1,2

                                       0
                                      12

                                     0.04
                                      17
O
                                     8.40
                                      1,2

                                       0
                                      18

                                     1.07
                                      20
O(CO)
                                     1.60
                                      1,2

                                       0
                                      18

                                     1.07
                                      17
O(HCO)
                                     0.90
                                     1,15

                                       0
                                      18

                                     1.07
                                      17
O(NO2)
                                     0.04
                                      1,2

                                       0
                                      18

                                     1.07
                                      17
O(OH)
                                     16.26
                                       5

                                       1
                                      18

                                     1.07
                                      17
NO2
                                     0.00
                                      1,2

                                       0
                                      12

                                      n/a

F
                                     0.09
                                      1,2

                                       0
                                      12

                                     0.01
                                      21
Cl
                                     0.38
                                      1,2

                                       0
                                      12

                                     0.01
                                      13
Br
                                     0.28
                                      1,2

                                       0
                                      12

                                      n/a

I
                                     0.53
                                      1,2

                                       0
                                      12

                                      n/a

C=C
                                     1.00
                                      1,2

                                       1
                                      17

                                     0.95
                                      13
ONO2
                                     0.04
                                      1,2

                                       0
                                      12

                                     0.12

Bz
                                       1
                                      14