Aromatics contribute a significant fraction to organic compounds in the
troposphere and are mainly emitted by anthropogenic activities and biomass burning. Their oxidation in lab experiments is known to lead to the formation of ozone and aerosol precursors. However, their overall impact on
tropospheric composition is uncertain as it depends on transport,
multiphase chemistry, and removal processes of the oxidation intermediates.
Representation of aromatics in global atmospheric models has been either
neglected or highly simplified. Here, we present an assessment of their impact on gas-phase chemistry, using the general circulation model EMAC (ECHAM5/MESSy Atmospheric Chemistry). We employ a comprehensive kinetic model to represent the oxidation of the following monocyclic aromatics: benzene, toluene, xylenes, phenol, styrene, ethylbenzene, trimethylbenzenes, benzaldehyde, and lumped higher aromatics that contain more than nine C atoms.
Significant regional changes are identified for several species. For instance, glyoxal increases by 130 % in Europe and 260 % in East
Asia, respectively. Large increases in HCHO are also predicted in these regions. In general, the influence of aromatics is particularly evident in areas with high concentrations of NOx, with increases up to 12 % in O3 and 17 % in OH.
On a global scale, the estimated net changes of trace gas levels are minor when aromatic compounds are included in our model. For instance, the tropospheric burden of CO increases by about 6 %, while the burdens of OH, O3, and NOx (NO+NO2) decrease between 3 % and 9 %. The global mean changes are small, partially because of compensating effects between high- and low-NOx regions. The largest change is predicted for the important aerosol precursor glyoxal, which increases globally by 36 %. In contrast to other studies, the net change in tropospheric ozone is predicted to be negative, -3 % globally. This change is larger in the Northern Hemisphere where global models usually show positive biases. We find that the reaction with phenoxy radicals is a significant loss for
ozone, on the order of 200–300 Tg yr-1, which is similar to the estimated ozone loss due to bromine chemistry.
Although the net global impact of aromatics is limited, our results indicate that aromatics can strongly influence tropospheric chemistry on a regional scale, most significantly in East Asia. An analysis of the main model uncertainties related to oxidation and emissions suggests that the impact of aromatics may even be significantly larger.
Introduction
Volatile organic compounds (VOCs) comprise a large variety of species
that influence the tropospheric chemistry at local, regional, and global scales. VOCs react mainly with the hydroxyl radical (OH), ozone (O3), and the nitrate radical (NO3), or they are
photolyzed. Their oxidation affects many key atmospheric species,
including OH, O3, and nitrogen oxides
(NOx=NO+NO2). Production and destruction of ozone is controlled by the ratio of VOCs to NOx. In the low-NOx regime, the net effect of VOC oxidation is ozone destruction. Under high-NOx conditions, e.g., in urban areas, O3 is generated by the oxidation of VOCs .
Aromatics are unsaturated planar cyclic organic compounds with enhanced stability due to a strong electron delocalization. Several of them are present in the atmosphere, e.g., benzene, toluene, ethylbenzene, xylenes, styrene, and trimethylbenzenes. In general, aromatic compounds are found in continental areas, especially in industrialized urban and semiurban regions where their emissions are highest. They are responsible for a considerable fraction of ozone and secondary organic aerosol (SOA) formation . In addition, many aromatics are toxic .
Emissions of aromatics are primarily anthropogenic, related to fuel combustion and leakage from fuels and solvents . Emissions from biomass burning play a secondary role but can be important on a regional scale . Biogenic emissions are only relevant for toluene , although recent studies suggest that other aromatics from biogenic sources may rival those from fossil fuel use .
As shown by , aromatic compounds are removed from the atmosphere mainly via chemical oxidation. Dry deposition is a minor sink and wet deposition is almost negligible. The gas-phase chemistry of aromatics has been the subject of many studies e.g.,. Aromatics have relatively short atmospheric lifetimes ranging from a few hours, e.g., for trimethylbenzene, to about 10 days, e.g., for benzene . Their oxidation is mainly controlled by the OH radical but they also react with NO3 and O3. The reaction with OH can proceed along two principal pathways. The first starts with H abstraction from an aliphatic substituent. The following reactions are similar to those of aliphatic compounds and involve the addition of O2, yielding a peroxy radical as an intermediate. Toluene, for example, can be oxidized in this way to benzaldehyde and references therein. The second, which is the dominant path, is OH addition to the aromatic ring. Secondary reactions can lead to ring opening and complex further reactions, eventually generating HCHO, glyoxal, and other smaller organic molecules and references therein. The products from the oxidation of aromatic compounds have a reduced volatility and allow for the formation of SOA , which in turn can significantly reduce the gas-phase concentrations of the aromatic oxidation products.
Numerical models are essential to understand the highly complex chemical
degradation of aromatics and to quantify the impact of these compounds
in atmospheric chemistry. A very detailed modeling of aromatics is possible with the reactions contained in the Master Chemical Mechanism MCM,. However, due to its complexity, the full mechanism
is mainly suitable for box model calculations. For global studies,
simplified reaction schemes are usually used e.g.,.
The main objective of this study is to investigate how tropospheric OH, O3, NOx, and several VOC concentrations are affected by the oxidation of several monocyclic aromatics. The paper is
organized as follows. In Sect. , the numerical
model and the set-up of the simulations are described. Section analyzes the calculated impact on selected chemical
species both on the global and regional scales.
Model description
We used the ECHAM5/MESSy Atmospheric Chemistry (EMAC) model, which is a numerical chemistry and climate simulation system that includes submodels describing tropospheric and middle atmosphere processes . EMAC uses the second version of the Modular Earth Submodel System (MESSy2) to link multi-institutional computer codes. The core atmospheric model is the fifth-generation European Centre Hamburg general circulation model ECHAM5,.
For the present study we performed simulations with EMAC (ECHAM5 version 5.3.02, MESSy version 2.53) in the T63L31ECMWF resolution, which corresponds to a grid with a horizontal cell size of approximately 1.875∘×1.875∘ and 31 vertical hybrid pressure levels, extending from the surface up to about 10 hPa.
Emission rates of the individual aromatics are shown in Table . The sum of all sources is 29.4 Tg C yr-1.
For anthropogenic emissions, we used EDGAR 4.3.2 distributed vertically as in . The MESSy submodel MEGAN calculates biogenic emissions . For biomass burning, the submodel BIOBURN was used, which integrates the Global Fire Assimilation System (GFAS) inventory .
Global annual emission rates of aromatic compounds included in the model simulations and their relative contributions.
Atmospheric chemistry was calculated with the MECCA submodel, which has been evaluated by and . The most recent model version has been described by . The mechanism for aromatic species is a reduced version of the MCM , as described in detail by . In short, the MCM schemes for benzene and toluene were taken. Following the approach of , short-lived intermediates were replaced with their stable products and isomeric peroxy radicals were lumped, preserving the yield of stable products. The initial oxidation steps of aromatics other than benzene and toluene are considered and products replaced by the analogous toluene oxidation products. This approximation inherently introduces an error with respect to the formation of larger and low volatile products. The carbon mass that is not accounted for with this approximation is however tracked by introducing the counter LCARBON for the difference of carbon atoms between the oxidation products of larger aromatics and toluene. In this study, we consider several additions to the MCM reactions:
For several nitrophenols (MCM names: HOC6H4NO2, DNPHEN, TOL1OHNO2,
MNCATECH, DNCRES), their photolytic production of HONO were added
, e.g.,
.In JVAL the cross sections for 2-nitrophenol and 3-methyl-2-nitrophenol and the quantum yield for 2-nitrophenol by are used to calculate the j values.
For the photolysis of benzaldehyde, the MCM uses the rate constant
(j value) of methacrolein as a proxy. Our model calculates with JVAL
the j value based on the UV/VIS spectrum of benzaldehyde
recommended by . In our code, the photolysis of benzaldehyde
produces C6H5O2, HO2, and CO.
For several phenyl peroxy compounds (MCM names: C6H5O2, CATEC1O2,
OXYL1O2, MCATEC1O2, NCRES1O2), their reactions
with NO2 were added , e.g.,.
For the reaction of HO2 with the peroxy radical C6H5CO3
(resulting from the oxidation of benzaldehyde), we use the yields
provided by .
Alkyl nitrate yields are calculated as a function of temperature
and pressure, as described by .
Bicyclic peroxy radicals in the oxidation mechanism of toluene
yield 60 % glyoxal and 40 % methyl glyoxal from the nonradical terminating reactions with NO and HO2 as suggested by . Benzene is treated analogously but yields 100 % glyoxal from the above-mentioned reactions.
The aerosol calculations follow the approach of , with the notable difference of the inclusion of the explicit organic aerosol submodel ORACLEv1.0 by . Although, similar to , low and intermediate volatiles are parameterized as lumped species; the equilibrium with their equivalent aerosol phase is explicitly calculated for ≃600 volatile organic carbon tracers via ORACLE. The volatility and the enthalpy of vaporization of each tracer is estimated with the approaches of and , respectively.
The simulated period covers the years 2009–2010, with the first year as spin-up and year 2010 being used for the analysis. The feedback between radiation and chemistry was decoupled to avoid any influence of chemistry on the dynamics (QCTM mode by ). As a consequence, every simulation discussed here has the same meteorology, i.e., binary identical transport.
To analyze the influence of the aromatic compounds on atmospheric chemistry and composition, we performed three model simulations, as listed in Table . The AROM simulation includes all chemical reactions and emissions of the following monocyclic aromatic compounds: benzene, toluene, xylenes (lumped), phenol, styrene, ethylbenzene, trimethylbenzenes (lumped), benzaldehydes, and higher aromatics (as representative of aromatics with more than nine carbon atoms). The reference simulation (NOAROM) is identical to AROM except that it excludes aromatic compounds. In the OnlyMCM run, we reverted the additions and changes to the MCM that have been described above. Our focus is to compare AROM with NOAROM. The results of OnlyMCM are mainly interesting for benzaldehyde and HONO. As EMAC uses terrain-following vertical hybrid pressure coordinates, we will refer to “surface” as the lowest model level, with an average thickness of roughly 60 m.
Annual mean mixing ratios of the sum of aromatics at the
surface (a) and the zonal mean (b) in the AROM simulation. The solid line between 100 and 300 hPa depicts the mean tropopause level.
Sensitivity studies.
SimulationDescriptionAROMAromatics are fully includedNOAROM (reference)No aromatics (emissions switched off)OnlyMCMOnly MCM reactionsResults and discussion
Globally averaged surface mixing ratios obtained from all model simulations (AROM, NOAROM, and OnlyMCM) are listed in Table . Figure shows the annual average mixing ratios of the sum of all aromatic compounds included in the simulation AROM. They are higher in continental areas and close to the surface. The highest values are predicted in the Northern Hemisphere (NH), in particular, in East and South Asia, as well as in parts of Europe, Africa, and the US, reaching up to about 1 nmol mol-1. The background mean mixing ratios in oceanic areas of the Southern Hemisphere (SH) are on the order of a few pmol mol-1. For a more detailed analysis, we have selected the following five regions, as defined in Fig. : Amazon area (AMA), central Africa (CAF), eastern Asia (EAS), Europe (EUR), and eastern US (EUS). The budgets of selected chemical species were calculated within these regions (Table ).
Globally averaged area-weighted mixing ratios at the surface (annual
averages for 2010). “ABSDIFF” denotes the absolute difference, (e.g., AROM-NOAROM), and “RELDIFF” the relative difference, (e.g., AROM/NOAROM-1).
AROM vs. OnlyMCM AROM vs. NOAROM NOAROMOnlyMCMAROMABSDIFFRELDIFFABSDIFFRELDIFFmol mol -1mol mol -1mol mol -1mol mol -1%mol mol -1%OH4.630×10-144.472×10-144.487×10-141.557×10-160.3482-1.425×10-15-3.078O33.269×10-83.220×10-83.190×10-8-2.964×10-10-0.9204-7.888×10-10-2.413NO3.029×10-112.793×10-112.609×10-11-1.843×10-12-6.599-4.203×10-12-13.87NO23.389×10-103.314×10-103.191×10-10-1.228×10-11-3.706-1.977×10-11-5.834NO31.004×10-129.462×10-131.080×10-121.339×10-1314.157.599×10-147.568HONO7.393×10-127.260×10-127.315×10-125.538×10-140.7628-7.754×10-14-1.049HNO31.420×10-101.393×10-101.426×10-103.352×10-122.4076.607×10-130.4653HCHO5.993×10-105.992×10-106.002×10-109.484×10-130.15838.414×10-130.1404Glyoxal1.040×10-111.444×10-111.505×10-116.117×10-134.2374.646×10-1244.67Methyl glyoxal3.847×10-114.005×10-114.015×10-111.051×10-130.26251.682×10-124.372Benzaldehyde6.798×10-124.479×10-12-2.319×10-12-34.114.479×10-12CO97.6×10-9103.3×10-9103.3×10-9-6.5×10-11-0.062785.7×10-95.847
Selected regions: AMA – Amazon area, CAF – central Africa,
EAS – eastern Asia, EUR – Europe, and EUS – eastern US.
Hydroxyl radical (OH)
Figure shows the model-calculated surface OH in
the AROM and NOAROM simulations. When aromatics are introduced to the model, the global average concentration of OH decreases for two reasons: first, the direct reaction with aromatics consumes OH, and second, additional CO resulting from the degradation of aromatics represents an increased sink for OH. However, in eastern Asia, Europe, and the east coast of the US, where NOx concentrations are high, an increase of OH can be seen. Although the aromatics decrease NOx in these areas (see below), the chemical system remains in the high-NOx regime.
Annual average OH mixing ratios at the surface. Middle rows: seasonal means. Left column: mixing ratios in the NOAROM simulation. Middle column: absolute difference AROM–NOAROM. Right column: relative difference
AROM/NOAROM-1 in % (shown only where OH is above 0.01 pmol mol-1).
We find that inclusion of aromatics emissions leads to an increase of OH in these regions but to a decrease in the low-NOx CAF region.
The increased OH in the high-NOx regions is mainly caused by the reaction of NO with HO2. The production of OH from this important reaction is enhanced by the significant HO2 formation in aromatics oxidation. Compared to OnlyMCM, the AROM simulation has additional HO2 production from the photolysis of ortho-nitrophenols (Reaction ) and benzaldehyde (Sect. ). The enhanced HO2 levels (not shown) overcompensates the negative changes in NO (see Sect. ).
Simulated tropospheric integrals of OH, O3 and NOx, and the lifetime τ of CH4. Tropospheric burdens were reckoned using six different tropopause definitions (provided by the TROPOP submodel, see for details): (1, 2) surfaces of O3 mixing ratio of 125 and 150 nmol mol-1, respectively,
(3) WMO definition , (4) dynamic PV-based (3.5 PVU potential vorticity surface, sought within 50–800 hPa), (5) climatological (invariable zonal profile, i.e., 300–215 ⋅ (cos(latitude))2 hPa) and (6) the combined definition (WMO tropopause within 30∘ N–30∘ S, otherwise dynamic PV-based tropopause). The latter definition is used by default in EMAC and in this work. Estimated changes to tropospheric O3 burden are identical within 0.05 % between the available definitions.
Figure shows the seasonal cycle of the OH mixing ratio in the planetary boundary layer for the NH and SH. Inclusion of the aromatics leads to a relative decrease between 2.5 % and 5.5 %. Higher
OH concentrations are identified over continental areas during the NH autumn, winter, and spring than in summer (Fig. ). In summer, OH concentrations increase only at a few locations when aromatics are included. In general, enhancements are predicted for regions where radical production is not NOx-limited. In the NH there are obviously more such regions compared to the SH. However, the largest decrease in the planetary boundary OH is computed for the NH where most of the emissions of aromatics are located.
(a) Seasonal cycles of OH daily (24 h) mixing ratio means (in 10-14 mol mol-1) in the planetary boundary layer (PBL) for AROM (solid line) and NOAROM (dashed line). (b) Relative difference
(expressed in %) between AROM and NOAROM. In blue, values for the NH; in red, values for the SH. The PBL diagnosis is described in . The PBL is calculated in the model based on the work of . An interactive calculation is performed following the approach of , using the Richardson number, the horizontal velocity components, the buoyancy parameters, and the virtual temperature .
Figure shows the annual zonal mean changes of the OH mixing ratio. The changes are most pronounced in the NH upper troposphere where reductions range from 7 % to 20 %. These predicted changes are associated with similar reductions in NOx. In fact, the upper troposphere is in general NOx-limited and the oxidation of aromatics enhances the formation N2O5 and HNO3, which are lost heterogeneously. This leads to an effective removal of NOx from the gas phase and lowers the radical production. The change in hemispheric burdens of OH are consistent with this picture (Table ). This moderately helps bring the model-simulated interhemispheric OH asymmetry closer to
that derived from observations . Globally, aromatics oxidation reduces OH by 7.7 % and consequently increases methane lifetime by about 5.5 %. The changes are more pronounced in the Northern Hemisphere where aromatics are mostly emitted (Table ). However, in the EMAC model, methane lifetime remains significantly lower than the ACCMIP multimodel mean and the observational-based estimates . Coarse model spatial resolutions (about 200 km) are known to result in an overestimation (underestimation) of global mean OH (methane lifetime) of at least 5 % . This is due to a less efficient conversion of NOx to NOy when strong pollutant emissions are artificially diluted in the model grid boxes. This aspect certainly has a larger impact on the interhemispheric OH asymmetry in atmospheric models that is in contrast to observational estimates .
Annual average zonal mean OH mixing ratios. (a) Mixing ratios in the NOAROM simulation. (b) Absolute difference AROM–NOAROM. (c) Relative difference AROM/NOAROM-1 in %. The solid line between 100 and 300 hPa depicts the mean tropopause level.
Annual average O3 mixing ratios at the surface. (a) Mixing ratios in the NOAROM simulation. (b) Absolute
difference AROM–NOAROM. (c) Relative difference AROM/NOAROM-1 in %.
Regionally averaged mixing ratios of selected species (annual averages for 2010).
Differences for OH between the AROM and OnlyMCM simulations are shown in Figs. and of the Appendix .
Ozone (O3)
In most areas of the globe, surface ozone is slightly lower in AROM than in NOAROM (Fig. ). The O3 reduction is due to (i) the decrease in NOx concentrations (limiting ozone formation) and (ii) increasing radical production (HOx and RO2) in ozone-depleting regimes, which enhances reaction of O3 with HO2. Only a few high-NOx regions, where hydrocarbons are the limiting factor for ozone formation, show increased ozone concentrations: mainly eastern Asia (EAS), but also the eastern US (EUS) and Europe (EUR). The increases in these areas is associated with anthropogenic emissions of aromatics, which have significant ozone formation potentials. We find that anthropogenic emissions of aromatics leads to an increase of O3 in the EAS and EUR regions but to a decrease in the low-NOx CAF region.
The seasonal cycles of the relative differences show lower amplitude than for OH, but similar patterns (Fig. ). The impact of aromatics is smallest in summer. Like for the OH levels, the interhemispheric asymmetry in the emission of aromatics determines the higher O3 decrease in the NH compared to the SH.
Same as in Fig. for ozone (the unit in the left plot is 10-8 mol mol-1.
Annual average zonal mean O3 mixing ratios. (a) Mixing ratios in the NOAROM simulation. (b) Absolute
difference AROM–NOAROM. (c) Relative difference AROM/NOAROM-1 in %. The solid line between 100 and 300 hPa depicts the mean tropopause level.
The zonal mean changes of O3 mixing ratio in the troposphere are uniformly negative (Fig. ). Similar to surface ozone, the annual mean changes for OnlyMCM and AROM are -2.3 % and -3.0 %, respectively. The hemispheric changes are shown in Table . It is well known that MCM for aromatics overestimates ozone production in chamber experiments . This issue has been analyzed in the companion paper , where the best mechanism improvement was found to be an early OH source during oxidation. introduced enhanced HOx sources by photolysis of benzaldehyde and nitrophenols. These modifications consistently result in less ozone produced with respect to MCM. These results deviate from the results by who suggested a global increase of 0.4 % due to aromatics. However, they only considered benzene, toluene, and xylenes. Our results, obtained with a more comprehensive setup, suggest that aromatics could slightly ameliorate the model overestimate in the NH . The overall tropospheric ozone burden decreases from 381 to 369 Tg for the AROM simulation. These estimated changes are robust against the tropopause definition and are about -3.5 % and -2.3 % for the Northern Hemisphere and Southern Hemisphere, respectively (Table ).
The changes in ozone are caused by perturbations of the radical production in different NOx regimes but also by the direct ozone loss in reactions with organic compounds. It is widely acknowledged that this direct loss is only due to the ozonolysis of unsaturated VOCs and is estimated to be about 100 Tg yr-1, less than 2 % of the tropospheric ozone budget (e.g., in ). However, with aromatics a new direct ozone loss process involving organic radicals comes into play. In Fig. the change in tropospheric ozone burden is shown against the change in ozone loss with organic compounds. This change is estimated to be globally in the 200–300 Tg yr-1 range depending on the mechanism used and is comparable to the loss by bromine chemistry in the troposphere . Ozone is known to react with organic radicals like methyl peroxy radical although this loss is an insignificant sink . We find that (substituted) phenoxy radicals from aromatics are a significant sink term of ozone (>200 Tg yr-1). These radicals are unique to aromatics oxidation and they also react with NO and NO2. When the concentrations of NOx are relatively low, C6H5O has a sufficiently long lifetime to react with O3. This ozone loss is modeled based on the results by for phenoxy radical:
.
Although the known rate constant for Reaction () is about one order of magnitude lower than the others, the high abundance in the atmosphere makes ozone the major sink of (substituted) phenoxy radicals.
This direct ozone loss in Reaction () is enhanced by phenoxy radical production in Reaction () and the concurrent loss of odd oxygen by NO3 photolysis and N2O5 heterogeneous loss:
4NO3+hν→NO+O25NO3+NO2→N2O5N2O5+H2O→2HNO3(aq).
In our chemical kinetics mechanism (also in MCM) the reaction system just
described constitutes an effective catalytic destruction cycle of odd oxygen.
The strength of this cycle has not been diagnosed in this study. Nevertheless, we observe that it depends on the (substituted) phenoxy radical levels and is significantly reduced in AROM compared to OnlyMCM (Fig. ). We ascribe this difference to MCM not accounting for the photolysis of nitrophenols (Reaction ) as determined by . In fact, in MCM, the first nitrophenols from benzene (HOC6H4NO2) and toluene (TOL1OHNO2) solely form nitrophenoxy radicals with the same reactivity of the unsubstituted phenoxy radical (C6H5O). Thus, the photolysis of nitrophenols decreases the amount of ozone lost by reaction with nitrophenoxy radicals. The impact of all the additions and modifications to the MCM on the predicted O3 levels is shown in Figs. and of Appendix . Uncertainties on the reactions mentioned in this paragraph are discussed in Sect. .
Change in tropospheric ozone burden versus change in ozone loss for all reactions in the VOC chemistry (G4 category of the MECCA mechanism, see the Supplement of ). The change in ozone loss is due to the reactions with (substituted) phenoxy radicals. Global and hemispheric results for OnlyMCM (blue) and AROM (orange) simulations are shown.
Annual average NO mixing ratios at the surface. (a) Mixing ratios in the NOAROM simulation. (b) Absolute difference AROM–NOAROM. (c) Relative difference AROM/NOAROM-1 in % (shown only where NO is above 10 pmol mol-1).
Annual average NO2 mixing ratios at the surface.
(a) Mixing ratios in the NOAROM simulation. (b) Absolute difference AROM–NOAROM. (c) Relative difference AROM/NOAROM-1 in % (shown only where NO2 is above 100 pmol mol-1).
Our results for ozone differ both in magnitude and sign compared to the global study by . However, the latter used the SAPRC-11
oxidation mechanism , which does not account for the reaction of phenoxy radicals with ozone (Reaction ) and phenylperoxy radicals with NO2 (Reaction ).
Inorganic nitrogen
The simulated annual mean NOx concentrations at the surface are significantly lower in AROM than in NOAROM
(Figs. and ). One
reason is the formation of aromatic species containing nitrogen (e.g.,
nitrophenols) in AROM, thereby transferring part of the NOx burden to the nitrogenated species. The largest decreases (both absolute and relative) are found in regions with high NOx concentrations. Since the ozone chemistry is not NOx-limited in these regions, the impact on ozone is small. This holds for the free troposphere for which zonal average decreases in NOx can be larger than 20 % (not shown), which in turn significantly influences OH (Fig. ).
On the one hand, the reaction with aromatics is a sink for NO3.
On the other hand, NO3 is produced in the phenylperoxy reaction
with NO2 (Reaction ). However the latter seems to dominate and cause a significant and widespread increase in the predicted NO3 levels. Relative to NOAROM, in AROM the global average of the nighttime species NO3 increases by more than 7 % (Table ). In contrast to the global mean tendency, NO3 has modest decreases in several regions in Africa, South America, and India (Fig. ). These decreases correlate well with emissions from biomass burning. Differences for NO3 between the AROM and OnlyMCM simulations are shown in Fig. of Appendix .
Annual average NO3 mixing ratios at the surface.
(a) Mixing ratios in the NOAROM simulation. (b) Absolute difference AROM–NOAROM. (c) Relative difference AROM/NOAROM-1 in % (shown only where NO3 is above 1 pmol mol-1).
Although the net change of global HONO is small (about 3 % less
in AROM than in NOAROM, see Fig. and Table ), the regional differences can be large (Table ). A decrease of HONO is seen mainly in polluted areas (EAS, EUR, EUS) in the winter. In contrast, HONO increases in the regions with emissions from biomass burning (AMA, CAF). Here, HONO is formed by the photolysis of nitrophenols (Reaction ). Since these reactions are not included in the MCM, we do not see any HONO increase in the OnlyMCM simulation (Fig. ).
Annual average HONO mixing ratios at the surface.
(a) Mixing ratios in the NOAROM simulation. (b) Absolute difference AROM–NOAROM. (c) Relative difference AROM/NOAROM-1 in % (shown only where HONO is above 1 pmol mol-1).
Annual average HONO mixing ratios at the surface. (a) Mixing ratios in the NOAROM simulation. (b) Absolute difference OnlyMCM–NOAROM. (c) Relative difference OnlyMCM/NOAROM-1 in % (shown only where HONO is above 1 pmol mol-1).
Annual average HNO3 mixing ratios at the surface.
(a) Mixing ratios in the NOAROM simulation. (b) Absolute difference AROM–NOAROM. (c) Relative difference AROM/NOAROM-1 in % (shown only where HNO3 is above 10 pmol mol-1).
Annual average HCHO mixing ratios at the surface.
(a) Mixing ratios in the NOAROM simulation. (b) Absolute difference AROM–NOAROM. (c) Relative difference AROM/NOAROM-1 in % (shown only where HCHO is above 100 pmol mol-1).
Annual average glyoxal mixing ratios at the surface. (a) Mixing ratios in the NOAROM simulation. (b) Absolute
difference AROM–NOAROM. (c) Relative difference AROM/NOAROM-1 in % (shown only where glyoxal is above 10 pmol mol-1).
Annual average methyl glyoxal mixing ratios at the surface.
(a) Mixing ratios in the NOAROM simulation. (b) Absolute difference AROM–NOAROM. (c) Relative difference AROM/NOAROM-1 in % (shown only where methyl glyoxal is above 10 pmol mol-1).
On a global average level, HNO3 is not affected much by aromatics.
However, an increase can be seen in the regions where ozone increases (EAS) or where biomass burning decreases NO3 and N2O5 (CAF) (see Fig. and Table ). An average zonal mean change of up to 5 % throughout the upper troposphere–lower stratosphere (UTLS) is linked to the enhanced NO3 production by Reaction ().
Selected oxygenated compounds
Globally, HCHO is not affected much by aromatics. There are,
however, regional differences that are moderate because of the concurrent enhancement of the HCHO sink by reaction with OH. We find maximum absolute depletions in the AMA region, where concentrations are typically high (Fig. ). Increased values of HCHO are mainly seen in EAS and EUR (Table ).
Annual average zonal mean CO mixing ratios. (a) Mixing ratios in the NOAROM simulation. (b) Absolute difference AROM–NOAROM. (c) Relative difference AROM/NOAROM-1 in %. The solid line between 100 and 300 hPa depicts the mean tropopause level.
α-dicarbonyls like glyoxal and methyl glyoxal are primarily produced from the bicycloalkyl-radical pathway leading in the case of benzene to BZBIPERO2 (MCM) . A minor secondary formation pathway from conjugated unsaturated dicarbonyls, e.g., MALDIAL (MCM), is also known and taken into account . As expected, the model predicts a very large increase of glyoxal in almost all continental areas (Figs. and ). The global burden is 36 % higher than in the NOAROM model simulation. The largest regional increases are in the EAS and EUR regions (Table ). An exception to the global trend is the AMA region, where OH is too low to produce either glyoxal or methyl glyoxal. Annual mean increases exceed 50 % over the continents close to the surface. In the lower troposphere, zonal mean increases are in the 10 %–20 % range. These changes are of significance for the model SOA budget since these two dicarbonyls are estimated to produce a large fraction of SOA by cloud
processing yielding low-volatile oligomers . However,
a model assessment of SOA formation from α-dicarbonyls is beyond the scope of this study. The reason is that although the simulations were performed with a VBS-based approach to model condensation of organic vapors, the EMAC model version used in this study has no representation of oligomer formation from (methyl)glyoxal. This has been recently implemented explicitly for cloud droplets and its effect is planned to be assessed in a subsequent study together with the contribution of reactive uptake of epoxides from isoprene and aromatics.
Comparing AROM to OnlyMCM, benzaldehyde decreases by more than 50 % when the photolysis rate constant (j value) from the MCM (based on methacrolein) is replaced by our value (based on the UV/VIS spectrum of benzaldehyde). The more realistic photolysis rate enhances the production of radicals like HO2.
Since additional reactive carbon compounds have been introduced in the
model, the oxidation of aromatics produces more CO, which has a lifetime of about 1–3 months . CO can travel long distances from its source, although its lifetime is not long enough to allow it to cross hemispheres . CO concentrations generally increase on the global scale, indicating a small addition to the carbon budget. When comparing AROM to NOAROM, we find an increase of about 6 % in the atmospheric
burden of CO. Interestingly, maximum zonal average increases of
10 %–20 % are found for the NH UTLS region (Fig. ).
Model uncertainties
The model calculations presented in this work are associated with some uncertainties related to the oxidation kinetic model, emissions, and model resolution.
Gas-phase oxidation of aromatics is complex and the kinetic mechanism used in this study reflects the state of knowledge, advancements, and limitations in the mechanism that have recently been discussed . Recent progress has focused in particular on the source strength of aerosol precursors and not on the overall radical production, which also affects ozone. Nevertheless, our kinetic model makes use of only one rate constant for the Reaction () of phenoxy radicals with ozone . It also assigns this rate constant to the substituted phenoxy radicals other than C6H5O. Unfortunately, there is only one study of the rate constant of Reaction () at 298 K. Although the 2σ reported uncertainty is slightly larger than 10 %, the rate constant of 2.86×10-13 cm3 molecules-1 s-1 has to be regarded as a lower limit. On the other hand, experimental evidence for the product of Reaction (), being phenyl peroxy radical (C6H5O2), has not been found although it was expected. If the products are different, then the catalytic O3-destruction cycle illustrated in Sect. would not be in place. However, a significant amount of ozone loss via Reaction () and analogous reactions is to be expected. Moreover, the ozone loss is likely underestimated because of the model not accounting for the photolysis of nitrophenols forming nitrosophenoxy radicals. Different from the HONO-formation channel, which destroys the aromatic ring, channels yielding substituted phenoxy radicals may dominate
and thus enhance ozone loss.
Another source of uncertainty is the direct formation of epoxide upon addition of OH and subsequently by O2 as implemented in the MCM ranging from 11.8 %, for benzene, to 24 %, for trimethylbenzene . There is in fact consistent theoretical evidence that the epoxide formation pathway passes through a second O2 addition. This implies that the epoxide yield likely depends on the abundance of NO, HO2, and RO2and references therein. This uncertainty limits the reliability of the predicted SOA formation from reactive uptake of epoxides by aerosols .
Cloud chemistry of organic compounds is known to suppress gas-phase HOx production and directly consume ozone . The overall effect on ozone depends on the local chemical regime. In our study water-soluble products are set to only undergo wet deposition (dissolution and removal by precipitation). Their aqueous-phase chemistry might however have a nonnegligible effect on ozone and other oxidants. For instance, phenol is known to react very quickly with OH in the aqueous phase . Moreover, phenoxide anions from phenols react quickly with ozone . In particular, nitrophenols might be efficient ozone scavengers as they are stronger acids than unsubstituted phenols. A global assessment of cloud chemistry involving aromatics oxidation products is possible with the modeling system used here . However, considering the complexity of aqueous-phase oxidation of organic compounds, such an assessment is outside the scope of this study and deserves a dedicated model study.
In our study, biomass burning emissions of aromatics are potentially underestimated. In fact, based on the recent update by , we estimate that emissions might be up to 5 Tg yr-1 (65 %) higher than what is implemented in our model. Moreover, emissions from peat fires in 2010 (the simulated year) were up to a factor 15 lower than in the subsequent years . In general, the interannual variability of biomass burning is large and difficult to capture in a study such as the present one. However, it appears that the two major contributions to this variability are the peat fires in Indonesia and boreal forest fires, which are strongly favored by El Niño and heat waves, respectively. An early estimate of anthropogenic emissions of aromatics gave 16 Tg C yr-1. Two relatively recent datasets yield about 50 % higher emissions of 23 Tg C yr-1 for RCP and 22 Tg C yr-1 for EDGAR 4.3.2 . The latter is used in this study and lacks the biofuel burning emissions of phenol, benzaldehyde, and styrene. Interannual variability of anthropogenic emissions of aromatics is not well known but the decadal trends are known to be negative since the 1980s . Aromatics emissions from terrestrial vegetation have been long neglected or considered very low. However, suggested that aromatics emissions from biogenic sources may rival those from anthropogenic ones. In this study we used the same emission algorithm used in but get much lower emissions for toluene (about 0.3 Tg C yr-1 vs. 1.5 Tg C yr-1). However, suggest that emissions of aromatics and benzenoid compounds may be in the 1.4–15 Tg C yr-1 range. The major contributors are toluene and some benzenoids (oxygenated aromatics). The latter are mainly emitted during blossoming and stress-induced reactions by plants. The variability of their emissions is not very well quantified. For instance, the MEGAN model calculates their emission strengths based on the ones for carbon monoxide .
The spatial resolution of atmospheric models has a significant influence on the predicted levels of oxidants and nitrogen oxides. Generally, in polluted regions, the coarser the resolution the larger the ozone production per molecule of NOx will be . This is due to the artificial dilution of strong NOx emissions which, in reality, is efficiently converted to NOy by reacting with HOx. For instance, reducing the spatial resolution over polluted North America, Europe, and East Asia with a two-way nested regional model led to a 9.5 % reduction in the global tropospheric ozone burden . We have shown that at our model resolution of 1.875∘×1.875∘ aromatics are estimated to induce important increases in HOx (Fig. ) and decreases in NOx (Figs. and ) over continental polluted regions. Therefore, at much higher spatial resolutions the predicted enhancement of surface ozone by aromatics in those regions (Fig. ) might be reduced. Based on the results by we expect this effect to translate into a significant enhancement of the tropospheric ozone reduction reported in this study (Sect. ). A quantification of the model resolution effect on chemical regimes is at the moment computationally prohibitive with our very large chemical scheme running in the global EMAC model.
Finally, atmospheric levels of benzene and toluene simulated by our model were shown to underestimate many observations by at least 20 % . It is worth noting that in the total emissions of aromatics were even slightly higher (2.6 Tg C yr-1) than in the AROM
simulation. This underestimate could be explained by an overestimate of the chemical sink in the troposphere by reaction with hydroxyl radical. However, the annual global mean concentration of hydroxyl radicals is potentially 10 % too high , which cannot account for model concentration biases that are larger than 20 %. Therefore, we surmise that the impact of aromatics on the trace gas composition may be larger than estimated in this study.
Summary
This study investigated the effects of several monocyclic aromatics on the tropospheric gas-phase composition by means of the chemistry–climate model EMAC. When aromatics are introduced into our model calculations, large changes are seen for glyoxal and methyl glyoxal. For other species, our results show a relatively small importance of aromatics on the global scale. This is consistent with recent results by who used a simpler chemistry mechanism in the GEOS-Chem model. However, different from that study, we found a negative impact on global ozone. Our results also indicate that by including aromatics chemistry, free tropospheric OH is reduced, especially in the Northern Hemisphere. On a regional scale, the concentrations of several species change significantly, with the relatively largest impacts in East Asia where emissions are higher. Regions with high NOx concentrations show increases of OH and O3. However, since these increases are counteracted by decreases downwind, i.e., in remote areas where NOx concentrations are much lower, the net effects on large scales are small. Of the nitrogen compounds, mainly NO3 and HONO are affected by the aromatics chemistry.
We conclude that although the impact of aromatics is relatively minor on the global scale, it is important on regional scales, notably in the anthropogenic source regions and especially in those where NOx emissions are strongest. Given the uncertainties in the oxidation mechanisms and emissions, the results of our model calculations may underestimate the impact of aromatics on the tropospheric gas-phase composition.
AROM vs. OnlyMCM
In this appendix the impact of the modifications to the MCM chemistry (listed in Sect. ) on the model results are shown for the main atmospheric oxidants.
Hydroxyl radical (OH)
The differences at the surface are shown in Fig. .
Much of the increase in Fig. can be ascribed to the
enhanced HOx production by photolysis of benzaldehyde and HONO from Reaction (). That latter, which results from benzene chemistry, explain the significant enhancement across the UTLS (see Fig. ).
Ozone (O3)
The differences at the surface are shown in Fig. . It can be seen that a large part of the enhancement in surface ozone mixing ratio in Fig. is due to enhanced HOx production in regions that are not NOx-limited. The zonal mean change in ozone is minimal and slightly positive at the tropical UTLS (Fig. ).
Nitrate radical (NO3)
The differences at the surface are shown in Fig. .
It can be seen that the widespread enhancement of NO3 in Fig. is largely to be ascribed to the effect of phenylperoxy reaction with NO2 (Reaction ).
Annual average OH mixing ratios at the surface. (a) Mixing ratios in the OnlyMCM simulation. (b) Absolute difference AROM–OnlyMCM. (c) Relative difference AROM/OnlyMCM-1 in % (shown only where OH is above 0.01 pmol mol-1).
Annual average zonal mean OH mixing ratios. (a) Mixing ratios in the OnlyMCM simulation. (b) Absolute difference AROM–OnlyMCM. (c) Relative difference AROM/OnlyMCM-1 in %. The solid line between 100 and 300 hPa depicts the mean tropopause level.
Annual average O3 mixing ratios at the surface. (a) Mixing ratios in the OnlyMCM simulation. (b) Absolute
difference AROM–OnlyMCM. (c) Relative difference AROM/OnlyMCM-1 in %.
Annual average zonal mean O3 mixing ratios. (a) Mixing ratios in the OnlyMCM simulation. (b) Absolute
difference AROM–OnlyMCM. (c) Relative difference AROM/OnlyMCM-1 in %. The solid line between 100 and 300 hPa depicts the mean tropopause level.
Annual average NO3 mixing ratios at the surface.
(a) Mixing ratios in the OnlyMCM simulation. (b) Absolute difference AROM–OnlyMCM. (c) Relative difference AROM/OnlyMCM-1 in % (shown only where NO3 is above 1 pmol mol-1).
Code availability
The Modular Earth Submodel System (MESSy) is continuously further developed and applied by a consortium of institutions. The usage of MESSy and access to the source code is licensed to all affiliates of institutions that are members of the MESSy Consortium. Institutions can be a member of the MESSy Consortium by signing the MESSy Memorandum of Understanding. More information can be found on the MESSy Consortium web page (http://www.messy-interface.org, last access: 21 February 2021) .
Author contributions
DT, RS, AP, and DC wrote the manuscript. AP and DC performed the model simulations. DT and RS developed and analyzed the chemical mechanism. SB visualized the model results. SG performed extended budgeting of the species' chemical turnover. All coauthors contributed to the analysis of results and the writing of the paper.
Competing interests
The authors declare that they have no conflict of interest.
Special issue statement
This article is part of the special issue “The Modular Earth Submodel System (MESSy) (ACP/GMD inter-journal SI)”. It is not associated with a conference.
Acknowledgements
The authors want to acknowledge the use of the Ferret program for analysis and graphics in this paper. Ferret is a product of NOAA's Pacific Marine Environmental Laboratory (information is available at
https://ferret.pmel.noaa.gov/Ferret/,
last access: 21 February 2021). The work described in this paper has received funding from the Initiative and Networking Fund of the Helmholtz Association through the project “Advanced Earth System Modelling Capacity (ESM)” (information is available at https://www.esm-project.net/, last access: 21 February 2021).
Financial support
The article processing charges for this open-access publication were covered by a Research Centre of the Helmholtz Association.
Review statement
This paper was edited by Anne Perring and reviewed by Alexander Archibald and one anonymous referee.
ReferencesAndreae, M. O.: Emission of trace gases and aerosols from biomass burning – an updated assessment, Atmos. Chem. Phys., 19, 8523–8546,
10.5194/acp-19-8523-2019, 2019.Atkinson, R. and Arey, J.: Atmospheric degradation of volatile organic
compounds, Chem. Rev., 103, 4605–4638, 10.1021/cr0206420, 2003.
Atkinson, R., Aschmann, S. M., Arey, J., and Carter, W. P.: Formation of
ring-retaining products from the OH radical-initiated reactions of benzene and toluene, Int. J. Chem. Kinet., 21, 801–827, 1989.
Barletta, B., Meinardi, S., Sherwood Rowland, F., Chan, C.-Y., Wang, X., Zou,
S., Yin Chan, L., and Blake, D. R.: Volatile organic compounds in 43 Chinese
cities, Atmos. Environ., 39, 5979–5990, 2005.
Bejan, I., El Aal, Y. A., Barnes, I., Benter, T., Bohn, B., Wiesen, P., and
Kleffmann, J.: The photolysis of ortho-nitrophenols: a new gas phase source
of HONO, Phys. Chem. Chem. Phys., 8, 2028–2035, 2006.Birdsall, A. W., Andreoni, J. F., and Elrod, M. J.: Investigation of the role
of bicyclic peroxy radicals in the oxidation mechanism of toluene, J. Phys.
Chem. A, 114, 10655–10663, 10.1021/jp105467e, 2010.Bloss, C., Wagner, V., Bonzanini, A., Jenkin, M. E., Wirtz, K., Martin-Reviejo, M., and Pilling, M. J.: Evaluation of detailed aromatic mechanisms (MCMv3 and MCMv3.1) against environmental chamber data, Atmos. Chem. Phys., 5, 623–639, 10.5194/acp-5-623-2005, 2005a.Bloss, C., Wagner, V., Jenkin, M. E., Volkamer, R., Bloss, W. J., Lee, J. D.,
Heard, D. E., Wirtz, K., Martin-Reviejo, M., Rea, G., Wenger, J. C., and
Pilling, M. J.: Development of a detailed chemical mechanism (MCMv3.1) for the atmospheric oxidation of aromatic hydrocarbons, Atmos. Chem. Phys., 5,
641–664, 10.5194/acp-5-641-2005, 2005b.Cabrera-Perez, D., Taraborrelli, D., Sander, R., and Pozzer, A.: Global
atmospheric budget of simple monocyclic aromatic compounds, Atmos. Chem.
Phys., 16, 6931–6947, 10.5194/acp-16-6931-2016, 2016.Carter, W. P. and Heo, G.: Development of revised SAPRC aromatics mechanisms,
Atmos. Environ., 77, 404–414, 10.1016/j.atmosenv.2013.05.021, 2013.Chen, J., Wenger, J. C., and Venables, D. S.: Near-ultraviolet absorption cross sections of nitrophenols and their potential influence on tropospheric
oxidation capacity, J. Phys. Chem. A, 115, 12235–12242, 10.1021/jp206929r, 2011.Cheng, S.-B., Zhou, C.-H., Yin, H.-M., Sun, J.-L., and Han, K.-L.: OH produced from o-nitrophenol photolysis: a combined experimental and theoretical investigation, J. Chem. Phys., 130, 234311–234311,
10.1063/1.3152635, 2009.
Daniel, J. S. and Solomon, S.: On the climate forcing of carbon monoxide, J.
Geophys. Res., 103, 13249–13260, 1998.Deckert, R., Jöckel, P., Grewe, V., Gottschaldt, K.-D., and Hoor, P.: A quasi chemistry-transport model mode for EMAC, Geosci. Model Dev., 4, 195–206, 10.5194/gmd-4-195-2011, 2011.Emmons, L. K., Walters, S., Hess, P. G., Lamarque, J.-F., Pfister, G. G.,
Fillmore, D., Granier, C., Guenther, A., Kinnison, D., Laepple, T., Orlando,
J., Tie, X., Tyndall, G., Wiedinmyer, C., Baughcum, S. L., and Kloster, S.:
Description and evaluation of the Model for Ozone and Related chemical
Tracers, version 4 (MOZART-4), Geosci. Model Dev., 3, 43–67,
10.5194/gmd-3-43-2010, 2010.Epstein, S. A., Riipinen, I., and Donahue, N. M.: A semiempirical correlation
between enthalpy of vaporization and saturation concentration for organic
aerosol, Environ. Sci. Technol., 44, 743–748, 10.1021/es902497z, 2010.Field, R. J., Raghavan, N. V., and Brummer, J. G.: A pulse radiolysis
investigation of the reactions of bromine dioxide radical (BrO2.cntdot.) with
hexacyanoferrate(II), manganese(II), phenoxide ion, and phenol, J. Phys. Chem., 86, 2443–2449, 10.1021/j100210a040, 1982.Fu, T.-M., Jacob, D. J., Wittrock, F., Burrows, J. P., Vrekoussis, M., and
Henze, D. K.: Global budgets of atmospheric glyoxal and methylglyoxal, and
implications for formation of secondary organic aerosols, J. Geophys. Res.,
113, D15303, 10.1029/2007JD009505, 2008.Guenther, A. B., Jiang, X., Heald, C. L., Sakulyanontvittaya, T., Duhl, T.,
Emmons, L. K., and Wang, X.: The Model of Emissions of Gases and Aerosols
from Nature version 2.1 (MEGAN2.1): an extended and updated framework for
modeling biogenic emissions, Geosci. Model Dev., 5, 1471–1492,
10.5194/gmd-5-1471-2012, 2012.Heiden, A. C., Kobel, K., Komenda, M., Koppmann, R., Shao, M., and Wildt, J.:
Toluene emissions from plants, Geophys. Res. Lett., 26, 1283–1286,
10.1029/1999GL900220, 1999.Henze, D. K., Seinfeld, J. H., Ng, N. L., Kroll, J. H., Fu, T.-M., Jacob, D. J., and Heald, C. L.: Global modeling of secondary organic aerosol formation from aromatic hydrocarbons: high- vs. low-yield pathways, Atmos. Chem. Phys., 8, 2405–2420, 10.5194/acp-8-2405-2008, 2008.Hoigné, J. and Bader, H.: Rate constants of reactions of ozone with organic and inorganic compounds in water – II: Dissociating organic compounds, Water Res., 17, 185–194, 10.1016/0043-1354(83)90099-4, 1983.
Holtslag, A. and Boville, B.: Local versus nonlocal boundary-layer diffusion in a global climate model, J. Climate, 6, 1825–1842, 1993.
Holtslag, A., De Bruijn, E., and Pan, H.: A high resolution air mass
transformation model for short-range weather forecasting, Mon. Weather Rev., 118, 1561–1575, 1990.Hu, L., Millet, D. B., Baasandorj, M., Griffis, T. J., Travis, K. R., Tessum,
C. W., Marshall, J. D., Reinhart, W. F., Mikoviny, T., Müller, M., Wisthaler, A., Graus, M., Warneke, C., and de Gouw, J.: Emissions of
C6–C8 aromatic compounds in the United States: Constraints from tall tower and aircraft measurements, J. Geophys. Res., 120,
826–842, 10.1002/2014JD022627, 2015.Huang, G., Brook, R., Crippa, M., Janssens-Maenhout, G., Schieberle, C., Dore, C., Guizzardi, D., Muntean, M., Schaaf, E., and Friedrich, R.: Speciation of anthropogenic emissions of non-methane volatile organic compounds: a global gridded data set for 1970–2012, Atmos. Chem. Phys., 17, 7683–7701, 10.5194/acp-17-7683-2017, 2017.Jagiella, S. and Zabel, F.: Reaction of phenylperoxy radicals with NO2 at 298 K, Phys. Chem. Chem. Phys., 9, 5036–5051, 10.1039/B705193J, 2007.Jenkin, M. E., Saunders, S. M., Wagner, V., and Pilling, M. J.: Protocol for
the development of the Master Chemical Mechanism, MCM v3 (part B):
tropospheric degradation of aromatic volatile organic compounds, Atmos. Chem.
Phys., 3, 181–193, 10.5194/acp-3-181-2003, 2003.Jöckel, P. and the MESSy Consortium: The highly structured Modular Earth Submodel System (MESSy), availabe at: http://www.messy-interface.org, last access: 21 February 2021.Jöckel, P., Kerkweg, A., Pozzer, A., Sander, R., Tost, H., Riede, H.,
Baumgaertner, A., Gromov, S., and Kern, B.: Development cycle 2 of the
Modular Earth Submodel System (MESSy2), Geosci. Model Dev., 3, 717–752,
10.5194/GMD-3-717-2010, 2010.Jöckel, P., Tost, H., Pozzer, A., Kunze, M., Kirner, O., Brenninkmeijer, C. A. M., Brinkop, S., Cai, D. S., Dyroff, C., Eckstein, J., Frank, F., Garny, H., Gottschaldt, K.-D., Graf, P., Grewe, V., Kerkweg, A., Kern, B., Matthes, S., Mertens, M., Meul, S., Neumaier, M., Nützel, M., Oberländer-Hayn, S., Ruhnke, R., Runde, T., Sander, R., Scharffe, D., and Zahn, A.: Earth System Chemistry integrated Modelling (ESCiMo) with the Modular Earth Submodel System (MESSy), version 2.51, Geosci. Model Dev., 9, 1153–1200, 10.5194/gmd-9-1153-2016, 2016.Kaiser, J. W., Heil, A., Andreae, M. O., Benedetti, A., Chubarova, N., Jones,
L., Morcrette, J.-J., Razinger, M., Schultz, M. G., Suttie, M., and van der
Werf, G. R.: Biomass burning emissions estimated with a global fire
assimilation system based on observed fire radiative power, Biogeoscience, 9,
527–554, 10.5194/bg-9-527-2012, 2012.Koppmann, R.: Volatile organic compounds in the atmosphere, Blackwell
Publishing Ltd, 10.1002/9780470988657, 2007.Lamarque, J.-F., Bond, T. C., Eyring, V., Granier, C., Heil, A., Klimont, Z.,
Lee, D., Liousse, C., Mieville, A., Owen, B., Schultz, M. G., Shindell, D.,
Smith, S. J., Stehfest, E., Van Aardenne, J., Cooper, O. R., Kainuma, M.,
Mahowald, N., McConnell, J. R., Naik, V., Riahi, K., and van Vuuren, D. P.:
Historical (1850–2000) gridded anthropogenic and biomass burning emissions
of reactive gases and aerosols: methodology and application, Atmos. Chem. Phys., 10, 7017–7039, 10.5194/acp-10-7017-2010, 2010.
Lee, S., Chiu, M., Ho, K., Zou, S., and Wang, X.: Volatile organic compounds (VOCs) in urban atmosphere of Hong Kong, Chemosphere, 48, 375–382, 2002.Lelieveld, J. and Crutzen, P. J.: Influences of cloud photochemical processes
on tropospheric ozone, Nature, 343, 227–233, 10.1038/343227A0, 1990.Lelieveld, J., Gromov, S., Pozzer, A., and Taraborrelli, D.: Global
tropospheric hydroxyl distribution, budget and reactivity, Atmos. Chem.
Phys., 16, 12477–12493, 10.5194/acp-16-12477-2016, 2016.Li, Y., Pöschl, U., and Shiraiwa, M.: Molecular corridors and
parameterizations of volatility in the chemical evolution of organic
aerosols, Atmos. Chem. Phys., 16, 3327–3344, 10.5194/acp-16-3327-2016,
2016.Lin, G., Penner, J. E., Sillman, S., Taraborrelli, D., and Lelieveld, J.:
Global modeling of SOA formation from dicarbonyls, epoxides, organic nitrates
and peroxides, Atmos. Chem. Phys., 12, 4743–4774,
10.5194/acp-12-4743-2012, 2012.Misztal, P., Hewitt, C., Wildt, J., Blande, J., Eller, A., Fares, S., Gentner, D., Gilman, J., Graus, M., Greenberg, J., Guenther, A., Hansel, A., Harley, P., Huang, M., Jardine, K., Karl, T., Kaser, L., Keutsch, F.,
Kiendler-Scharr, A., Kleist, E., Lerner, B., Li, T., Mak, J., Nölscher, A., Schnitzhofer, R., Sinha, V., Thornton, B., Warneke, C., Wegener, F., Werner, C., Williams, J., Worton, D., Yassaa, N., and Goldstein, A.: Atmospheric benzenoid emissions from plants rival those from fossil fuels, Scient. Rep., 5, 12064, 10.1038/srep12064, 2015.Naik, V., Voulgarakis, A., Fiore, A. M., Horowitz, L. W., Lamarque, J.-F., Lin, M., Prather, M. J., Young, P. J., Bergmann, D., Cameron-Smith, P. J., Cionni, I., Collins, W. J., Dalsøren, S. B., Doherty, R., Eyring, V., Faluvegi, G., Folberth, G. A., Josse, B., Lee, Y. H., MacKenzie, I. A., Nagashima, T., van Noije, T. P. C., Plummer, D. A., Righi, M., Rumbold, S. T., Skeie, R., Shindell, D. T., Stevenson, D. S., Strode, S., Sudo, K., Szopa, S., and Zeng, G.: Preindustrial to present-day changes in tropospheric hydroxyl radical and methane lifetime from the Atmospheric Chemistry and Climate Model Intercomparison Project (ACCMIP), Atmos. Chem. Phys., 13,
5277–5298, 10.5194/acp-13-5277-2013, 2013.Ng, N. L., Kroll, J. H., Chan, A. W. H., Chhabra, P. S., Flagan, R. C., and
Seinfeld, J. H.: Secondary organic aerosol formation from m-xylene, toluene, and benzene, Atmos. Chem. Phys., 7, 3909–3922,
10.5194/acp-7-3909-2007, 2007.Patra, P. K., Krol, M. C., Montzka, S. A., Arnold, T., Atlas, E. L., Lintner,
B. R., Stephens, B. B., Xiang, B., Elkins, J. W., Fraser, P. J., Ghosh, A.,
Hintsa, E. J., Hurst, D. F., Ishijima, K., Krummel, P. B., Miller, B. R.,
Miyazaki, K., Moore, F. L., Mühle, J., O'Doherty, S., Prinn, R. G.,
Steele, L. P., Takigawa, M., Wang, H. J., Weiss, R. F., Wofsy, S. C., and
Young, D.: Observational evidence for interhemispheric hydroxyl-radical
parity, Nature, 513, 219–223, 10.1038/nature13721, 2014.Paulot, F., Crounse, J. D., Kjaergaard, H. G., Kürten, A., St. Clair, J. M., Seinfeld, J. H., and Wennberg, P. O.: Unexpected epoxide formation in the gas-phase photooxidation of isoprene, Science, 325, 730–733,
10.1126/science.1172910, 2009.Pozzer, A., Jöckel, P., Tost, H., Sander, R., Ganzeveld, L., Kerkweg, A., and Lelieveld, J.: Simulating organic species with the global atmospheric
chemistry general circulation model ECHAM5/MESSy1: a comparison of model
results with observations, Atmos. Chem. Phys., 7, 2527–2550,
10.5194/ACP-7-2527-2007, 2007.Pozzer, A., Jöckel, P., and van Aardenne, J.: The influence of the vertical distribution of emissions on tropospheric chemistry, Atmos. Chem. Phys., 9, 9417–9432, 10.5194/ACP-9-9417-2009, 2009.Pozzer, A., Pollmann, J., Taraborrelli, D., Jöckel, P., Helmig, D., Tans, P., Hueber, J., and Lelieveld, J.: Observed and simulated global distribution and budget of atmospheric C2–C5 alkanes, Atmos. Chem. Phys., 10, 4403–4422, 10.5194/ACP-10-4403-2010, 2010.Pringle, K. J., Tost, H., Message, S., Steil, B., Giannadaki, D., Nenes, A.,
Fountoukis, C., Stier, P., Vignati, E., and Lelieveld, J.: Description and
evaluation of GMXe: a new aerosol submodel for global simulations (v1), Geosci. Model Dev., 3, 391–412, 10.5194/gmd-3-391-2010, 2010.Ran, L., Zhao, C., Geng, F., Tie, X., Tang, X., Peng, L., Zhou, G., Yu, Q., Xu, J., and Guenther, A.: Ozone photochemical production in urban Shanghai,
China: Analysis based on ground level observations, J. Geophys. Res., 114,
D15301, 10.1029/2008JD010752, 2009.Roeckner, E., Brokopf, R., Esch, M., Giorgetta, M., Hagemann, S., Kornblueh,
L., Manzini, E., Schlese, U., and Schulzweida, U.: Sensitivity of simulated
climate to horizontal and vertical resolution in the ECHAM5 atmosphere model,
J. Climate, 19, 3771–3791, 10.1175/JCLI3824.1, 2006.Rosanka, S., Sander, R., Wahner, A., and Taraborrelli, D.: Oxidation of low-molecular weight organic compounds in cloud droplets: development of the JAMOC chemical mechanism in CAABA/MECCA (version 4.5.0gmdd), Geosci. Model Dev. Discuss. [preprint], 10.5194/gmd-2020-337, in review, 2020.Roth, E., Chakir, A., and Ferhati, A.: Study of a benzoylperoxy radical in the gas phase: ultraviolet spectrum and C6H5C(O)O2+HO2 reaction between 295 and 357 K, J. Phys. Chem. A, 114, 10367–10379, 10.1021/jp1021467, 2010.
Sack, T. M., Steele, D. H., Hammerstrom, K., and Remmers, J.: A survey of
household products for volatile organic compounds, Atmos. Environ., 26A,
1063–1070, 1992.Sander, R., Jöckel, P., Kirner, O., Kunert, A. T., Landgraf, J., and Pozzer, A.: The photolysis module JVAL-14, compatible with the MESSy standard, and the JVal PreProcessor (JVPP), Geosci. Model Dev., 7, 2653–2662, 10.5194/GMD-7-2653-2014, 2014.Sander, R., Baumgaertner, A., Cabrera-Perez, D., Frank, F., Gromov, S.,
Grooß, J.-U., Harder, H., Huijnen, V., Jöckel, P., Karydis, V. A.,
Niemeyer, K. E., Pozzer, A., Riede, H., Schultz, M. G., Taraborrelli, D., and
Tauer, S.: The community atmospheric chemistry box model CAABA/MECCA-4.0,
Geosci. Model Dev., 12, 1365—1385, 10.5194/gmd-12-1365-2019, 2019.Sherwen, T., Schmidt, J. A., Evans, M. J., Carpenter, L. J., Großmann, K., Eastham, S. D., Jacob, D. J., Dix, B., Koenig, T. K., Sinreich, R., Ortega, I., Volkamer, R., Saiz-Lopez, A., Prados-Roman, C., Mahajan, A. S., and Ordóñez, C.: Global impacts of tropospheric halogens (Cl, Br, I) on oxidants and composition in GEOS-Chem, Atmos. Chem. Phys., 16, 12239–12271, 10.5194/acp-16-12239-2016, 2016.Sillman, S., Logan, J. A., and Wofsy, S. C.: The sensitivity of ozone to
nitrogen oxides and hydrocarbons in regional ozone episodes, J. Geophys. Res., 95, 1837–1851, 10.1029/JD095iD02p01837, 1990.Tao, Z. and Li, Z.: A kinetics study on reactions of C6H5O with
C6H5 and O3 at 298 K, Int. J. Chem. Kinet., 31, 65–72, 10.1002/(SICI)1097-4601(1999)31:1<65::AID-KIN8>3.0.CO;2-J, 1999.Taraborrelli, D., Lawrence, M. G., Butler, T. M., Sander, R., and Lelieveld,
J.: Mainz Isoprene Mechanism 2 (MIM2): an isoprene oxidation mechanism for
regional and global atmospheric modelling, Atmos. Chem. Phys., 9, 2751–2777,
10.5194/ACP-9-2751-2009, 2009.Tarr, M. A., Miller, W. L., and Zepp, R. G.: Direct carbon monoxide
photoproduction from plant matter, J. Geophys. Res.-Atmos., 100, 11403–11413, 10.1029/94JD03324, 1995.Tilmes, S., Lamarque, J.-F., Emmons, L. K., Kinnison, D. E., Marsh, D., Garcia, R. R., Smith, A. K., Neely, R. R., Conley, A., Vitt, F., Val Martin, M., Tanimoto, H., Simpson, I., Blake, D. R., and Blake, N.: Representation of the Community Earth System Model (CESM1) CAM4-chem within the Chemistry-Climate Model Initiative (CCMI), Geosci. Model Dev., 9, 1853–1890, 10.5194/gmd-9-1853-2016, 2016.Tost, H., Jöckel, P., Kerkweg, A., Sander, R., and Lelieveld, J.: Technical note: A new comprehensive SCAVenging submodel for global atmospheric chemistry modelling, Atmos. Chem. Phys., 6, 565–574,
10.5194/ACP-6-565-2006, 2006.Tost, H., Lawrence, M. G., Brühl, C., Jöckel, P., The GABRIEL Team, and The SCOUT-O3-DARWIN/ACTIVE Team: Uncertainties in atmospheric chemistry
modelling due to convection parameterisations and subsequent scavenging,
Atmos. Chem. Phys., 10, 1931–1951, 10.5194/ACP-10-1931-2010, 2010.
Troen, I. and Mahrt, L.: A simple model of the atmospheric boundary layer;
sensitivity to surface evaporation, Bound.-Lay. Meteorol., 37, 129–148, 1986.Tsimpidi, A. P., Karydis, V. A., Pozzer, A., Pandis, S. N., and Lelieveld, J.: ORACLE (v1.0): module to simulate the organic aerosol composition and
evolution in the atmosphere, Geosci. Model Dev., 7, 3153–3172,
10.5194/gmd-7-3153-2014, 2014.Tyndall, G. S., Wallington, T. J., and Ball, J. C.: FTIR Product Study of the
Reactions CH3O2+CH3O2 and CH3O2+O3, J. Phys. Chem. A, 102, 2547–2554, 10.1021/jp972784h, 1998.van der Werf, G. R., Randerson, J. T., Giglio, L., van Leeuwen, T. T., Chen,
Y., Rogers, B. M., Mu, M., van Marle, M. J. E., Morton, D. C., Collatz, G. J., Yokelson, R. J., and Kasibhatla, P. S.: Global fire emissions estimates during 1997–2016, Earth Syst. Sci. Data, 9, 697–720, 10.5194/essd-9-697-2017, 2017.Vereecken, L.: Reaction Mechanisms for the Atmospheric Oxidation of Monocyclic Aromatic Compounds, in: chap. 6, World Scientific Publishing, 377–527, 10.1142/9789813271838_0006, 2019.Vereecken, L., Chakravarty, H. K., Bohn, B., and Lelieveld, J.: Theoretical
Study on the Formation of H- and O-Atoms, HONO, OH, NO, and NO2 from the Lowest Lying Singlet and Triplet States in Ortho-Nitrophenol Photolysis, Int. J. Chem. Kinet., 48, 785–795, 10.1002/kin.21033, 2016.Volkamer, R., Platt, U., and Wirtz, K.: Primary and secondary glyoxal formation from aromatics: Experimental evidence for the bicycloalkyl-radical pathway from benzene, toluene, and p-xylene, J. Phys. Chem. A, 105, 7865–7874, 10.1021/jp010152w, 2001.
Wallington, T. J., Ammann, M., Cox, R. A., Crowley, J. N., Herrmann, H., Jenkin, M. E., McNeill, V., Mellouki, A., Rossi, M. J., and Troe, J.: IUPAC
Task group on atmospheric chemical kinetic data evaluation: Evaluated kinetic
data, available at: http://iupac.pole-ether.fr (last access: 31 July 2015), 2018.Warneck, P.: Chemistry of the Natural Atmosphere, 2nd Edn., Academic
Press, available at: https://www.sciencedirect.com/bookseries/international-geophysics/vol/71/suppl/C (last access: 22 February 2021), 1999.WMO: Definition of the tropopause and of significant levels, available at:
https://library.wmo.int/doc_num.php?explnum_id=6960 (last access: 15 November 2020), 1957.
WMO: Air quality guidelines for Europe, WHO Regional Office for
Europe, Copenhagen, 2000.Yan, Y., Lin, J., Chen, J., and Hu, L.: Improved simulation of tropospheric
ozone by a global-multi-regional two-way coupling model system, Atmos. Chem. Phys., 16, 2381–2400, 10.5194/acp-16-2381-2016, 2016.Yan, Y., Cabrera-Perez, D., Lin, J., Pozzer, A., Hu, L., Millet, D. B., Porter, W. C., and Lelieveld, J.: Global tropospheric effects of aromatic chemistry with the SAPRC-11 mechanism implemented in GEOS-Chem version 9-02, Geosci. Model Dev., 12, 111–130, 10.5194/gmd-12-111-2019, 2019.Young, P. J., Naik, V., Fiore, A. M., Gaudel, A., Guo, J., Lin, M. Y., Neu, J. L., Parrish, D. D., Rieder, H. E., Schnell, J. L., Tilmes, S., Wild, O., Zhang, L., Ziemke, J., Brandt, J., Delcloo, A., Doherty, R. M., Geels, C., Hegglin, M. I., Hu, L., Im, U., Kumar, R., Luhar, A., Murray, L., Plummer, D., Rodriguez, J., Saiz-Lopez, A., Schultz, M. G., Woodhouse, M. T., and Zeng, G.: Tropospheric Ozone Assessment Report: Assessment of global-scale model performance for global and regional ozone distributions, variability, and trends, Elementa: Sci. Anthrop., 6, 10, 10.1525/elementa.265, 2018.