Abstract. We evaluate monoterpene-derived peroxy radical (MT-RO2) unimolecular autoxidation and self and cross reactions with other RO2 in the GEOS-Chem global chemical transport model. Formation of associated highly oxygenated organic molecule (HOM) and accretion products are tracked in competition with other bimolecular reactions. Autoxidation is the dominant fate up to 6–8 km for first-generation MT-RO2 which can undergo unimolecular H-shifts. Reaction with NO can be a more common fate for H-shift rate constants < 0.1 s−1 or at altitudes higher than 8 km due to the imposed Arrhenius temperature dependence of unimolecular H-shifts. For MT-derived HOM-RO2, generated by multi-step autoxidation of first-generation MT-RO2, reaction with other RO2 is predicted to be the major fate throughout most of the boreal and tropical forested regions, while reaction with NO dominates in temperate and subtropical forests of the Northern Hemisphere. The newly added reactions result in ~4 % global average decrease of HO2 and RO2 mainly due to faster self-/cross-reactions of MT-RO2, but the impact upon HO2/OH/NOx abundances is only important in the planetary boundary layer (PBL) over portions of tropical forests. Within the bounds of formation kinetics and HOM photochemical lifetime constraints from laboratory studies, predicted HOM concentrations in MT-rich regions and seasons reach 10 % or even exceed total organic aerosol as predicted by the standard GEOS-Chem model. Comparisons to observations reveal large uncertainties remain for key reaction parameters and processes, especially the photochemical lifetime of HOM and associated accretion products. Using the highest reported yields and H-shift rate constants of MT-RO2 that undergo autoxidation, HOM concentrations tend to exceed the limited set of observations. Similarly, we infer that RO2 cross reactions rate constants near the gas-kinetic limit with accretion product branching greater than ~0.25 are inconsistent with total organic aerosol unless there is rapid decomposition of accretion products, the accretion products have saturation vapor concentrations > > 1 μg m−3, or modeled MT emission rates are overestimated. This work suggests further observations and laboratory studies related to MT-RO2 derived HOM and gas-phase accretion product formation kinetics, and especially their atmospheric fate, such as gas-particle partitioning, multi-phase chemistry, and net SOA formation, are needed.
Special Issue on “Recent Advances in Population Balance Modeling”
