Attribution of ground-level ozone to anthropogenic and natural sources of nitrogen oxides and reactive carbon in a global chemical transport model
- 1Institute for Advanced Sustainability Studies, Potsdam, Germany
- 2Institut für Meteorologie, Freie Universität Berlin, Germany
Correspondence: Tim Butler (firstname.lastname@example.org)
We perform a source attribution for tropospheric and ground-level ozone using a novel technique that accounts separately for the contributions of the two chemically distinct emitted precursors (reactive carbon and oxides of nitrogen) to the chemical production of ozone in the troposphere. By tagging anthropogenic emissions of these precursors according to the geographical region from which they are emitted, we determine source–receptor relationships for ground-level ozone. Our methodology reproduces earlier results obtained via other techniques for ozone source attribution, and it also delivers additional information about the modelled processes responsible for the intercontinental transport of ozone, which is especially strong during the spring months. The current generation of chemical transport models used to support international negotiations aimed at reducing the intercontinental transport of ozone shows especially strong inter-model differences in simulated springtime ozone. Current models also simulate a large range of different responses of surface ozone to methane, which is one of the major precursors of ground-level ozone. Using our novel source attribution technique, we show that emissions of NOx (oxides of nitrogen) from international shipping over the high seas play a disproportionately strong role in our model system regarding the hemispheric-scale response of surface ozone to changes in methane, as well as to the springtime maximum in intercontinental transport of ozone and its precursors. We recommend a renewed focus on the improvement of the representation of the chemistry of ship NOx emissions in current-generation models. We demonstrate the utility of ozone source attribution as a powerful model diagnostic tool and recommend that similar source attribution techniques become a standard part of future model intercomparison studies.
Tropospheric ozone plays a central role in the chemistry and self-cleansing capacity of the troposphere (Crutzen, 1973; Monks et al., 2015); however, at high concentrations close to the ground, it is harmful to human health (Haagen-Smit, 1952; Fleming et al., 2018) and vegetation (Reich and Amundson, 1985; Mills et al., 2018). As well as being transported into the troposphere through exchange with the stratosphere, ozone can be formed via chemical reactions in the troposphere involving two chemically distinct precursors: oxides of nitrogen (collectively NOx) and reactive carbon species, including carbon monoxide, methane, and volatile organic compounds (Crutzen, 1973; Atkinson, 2000). Increases in tropospheric ozone since preindustrial times have been attributed primarily to increases in anthropogenic emissions of NOx and methane – the latter of which is the most abundant reactive carbon species in the atmosphere (Wang and Jacob, 1998; Stevenson et al., 2013).
Ozone is long-lived enough in the troposphere to circumnavigate the Northern Hemisphere along the prevailing westerly winds (Jacob et al., 1999). Thus, emissions of NOx or reactive carbon in any Northern Hemisphere source region can contribute to the ozone mixing ratio in any other region of the Northern Hemisphere. This long-range contribution to the ozone mixing ratio is often referred to as “baseline” ozone (Parrish et al., 2017; Derwent et al., 2018). Due to seasonal variation in the lifetime of ozone, this effect is strongest in spring and weakest in summer (Fiore et al., 2009). The ambient ozone mixing ratio at any location is a combination of ozone transported from the hemispheric background and in-situ photochemical production. Recent analyses of long-term trends in baseline ozone in western regions of North America (Parrish et al., 2017) and Europe (Derwent et al., 2018) have shown increasing trends since reliable measurements began in the 1980s until approximately 2000–2010, and they indicate that these trends may be beginning to reverse.
Chemical transport models (CTMs) are commonly used to interpret observations of ozone and to synthesize an understanding of the fundamental processes controlling its origin and fate in the atmosphere in order to project future trends (e.g. Young et al., 2018). The range in values of the surface ozone mixing ratio over the Northern Hemisphere simulated by contemporary CTMs is extremely high (see, for example, our Fig. 1 in Sect. 3), requiring the use of a large ensemble of models (e.g. HTAP, 2010; Young et al., 2018). When compared with available measurements of ozone for the Northern Hemisphere (e.g. Schultz et al., 2017), ensembles of global CTMs are generally able to simulate the spatial distribution and seasonal cycles of surface ozone; however, they are consistently biased high in the Northern Hemisphere and have difficulty simulating long-term trends (Young et al., 2018). Potential sources of uncertainty in CTMs include uncertainties in their chemical mechanisms (the representations of the relevant chemical reactions and their rates); their representation of atmospheric transport processes; and their representation of exchange processes between the atmosphere and the surface of the Earth, including emissions of the ozone precursors NOx and reactive carbon.
The most important class of reactions for the formation of ozone in the troposphere is the reaction of NO (nitric oxide) with a peroxy radical, which is itself formed during the oxidation of reactive carbon (Atkinson, 2000). During this process, the NO is converted to NO2 (nitrogen dioxide), which can be rapidly photolysed, ultimately forming ozone and recycling NO. The ozone production efficiency of NOx (the combined concentration of NO and NO2) can vary significantly depending on the location and timing of the NOx emissions. In the polluted boundary layer, NOx is rapidly removed from the atmosphere through the reaction of NO2 with OH, forming HNO3, which is subsequently lost via dry or wet deposition. Under less polluted conditions, NO2 photolysis competes more effectively with HNO3 production, allowing each unit of NOx to react with a higher number of peroxy radicals before eventually being scavenged by OH, thereby leading to a higher ozone production efficiency per unit of NOx. When NOx is lofted into the free troposphere, its ozone productivity increases substantially (Jacob et al., 1996). Thus, emissions of NOx in the tropics are especially effective at producing tropospheric ozone due the fact that they are transported aloft by deep convection (Zhang et al., 2016). NOx emissions from both aircraft and lightning are also highly efficient at producing tropospheric ozone (Beck et al., 1992; Dahlmann et al., 2011). Combustion of fossil fuels is the largest source of NOx in the atmosphere (Galloway et al., 2008).
Lawrence and Crutzen (1999) first pointed out that international shipping can have a disproportionately high influence on tropospheric ozone due to the disperse nature of NOx emissions from this source. Hoor et al. (2009) quantified the sensitivity of tropospheric ozone to NOx emissions from different modes of transport (land, sea, and air), finding that aircraft emissions were most efficient at producing ozone (per molecule of NOx emitted), followed by ships and then land transport. For near-surface ozone, however, NOx emissions from ships were shown to have a higher influence than NOx emissions from aircraft. Nevertheless, ozone production from ship NOx is highly uncertain in current CTMs. Kasibhatla et al. (2000) and von Glasow et al. (2003) showed that global CTMs, due to their coarse resolution (usually in the hundreds of kilometres), do not resolve the chemistry of ship exhaust plumes; instead, there is a tendency for NOx to be removed from the atmosphere more quickly than simulated by the global CTMs, which effectively instantly dilute these emissions into very large volumes. Wild and Prather (2006) also showed that this effect applies more generally to other concentrated emission sources such as urban areas. Vinken et al. (2011) introduced a method for parameterizing ship exhaust plume chemistry using lookup tables in their global CTM, but this method has not been widely adopted by the modelling community. Modelling of ship NOx and the effects of NOx on the atmosphere remains a challenge for global CTMs.
The term “reactive carbon” encompasses a very wide range of atmospheric constituents (e.g. Chameides et al., 1992; Goldstein and Galbally, 2007; Heald and Kroll, 2020). In contrast to NOx, most of the reactive carbon emitted to the atmosphere is not of anthropogenic origin but rather emitted from the biosphere. In this study, we restrict our definition to molecules that yield peroxy radicals (either hydroperoxy radicals HO2 or organic peroxy radicals RO2) during their gas-phase oxidation and, thus, contribute to ozone formation by potentially converting NO to NO2. Therefore, this definition includes carbon monoxide (CO) and the large family of molecules known as volatile organic compounds (VOCs). The simplest VOC is methane, which is often considered separately from non-methane volatile organic compounds (NMVOCs) due to its very long lifetime in the troposphere. The ozone production potential of reactive carbon depends on the rate at which it is oxidized in the atmosphere, usually through reaction with the OH radical (Carter, 1994), as well as the subsequent chemistry of its oxidation products (Butler et al., 2011; Derwent, 2020). Most reactive carbon species have relatively short lifetimes in the troposphere due to their reaction with OH radicals. Methane, due to its exceptionally low reactivity, is well mixed in the troposphere. In contrast to other forms of reactive carbon, emissions of methane can contribute to ozone formation at any location in the troposphere where photochemical conditions are favourable (Fiore et al., 2008). Despite its low reactivity in comparison to other types of reactive carbon, methane is highly abundant and has been shown to make a large contribution to tropospheric ozone (Wang and Jacob, 1998; Fiore et al., 2008; Stevenson et al., 2013; Butler et al., 2018).
PAN (peroxyacetyl nitrate) is an important reservoir species for both peroxy radicals and for NOx (Fischer et al., 2014). Peroxyacetyl radicals are formed during the oxidation of a wide range of different types of NMVOCs from a wide range of different sources. PAN is formed through the reaction of peroxyacetyl radicals with NO2, primarily in the polluted boundary layer where both are abundant (Atkinson, 2000). The lifetime of PAN is strongly temperature dependent. At colder temperatures higher in the troposphere, PAN can be transported over long distances, and it can act as a source of NO2 and peroxyacetyl radicals in remote regions upon subsidence and thermal decomposition (Fischer et al., 2014). The chemical mechanisms and reaction rate constants involved in the formation and decomposition of PAN vary widely between CTMs (Emmerson and Evans, 2009; Knote et al., 2015), leading to large inter-model differences in simulated PAN (Emmons et al., 2015). Fiore et al. (2018) suggested that measurements of PAN at northern midlatitude mountaintop sites in spring could provide a useful constraint on CTMs, although the number of observations available is limited.
With careful interpretation, the results of ensembles of CTMs can be used to diagnose long-range transboundary transport of ozone and to develop intercontinental source–receptor relationships, which relate the effects of precursor emissions from different regions of the Northern Hemisphere to mixing ratios of ground-level ozone in other regions of the Northern Hemisphere. An example of this is the activity of the Task Force on Hemispheric Transport of Air Pollution (TF-HTAP; HTAP, 2010), which reports to the Convention on Long-Range Transboundary Air Pollution (CLRTAP) and, thus, informs international policymaking for the mitigation of air pollution. The TF-HTAP studies used a “perturbation” approach that compared a control simulation with sensitivity simulations in which emissions of particular ozone precursors were reduced by 20 %. Combined 20 % reductions of global average methane and remote anthropogenic emissions of NOx, CO, and NMVOCs were shown to have an approximately equal effect on annual average ozone as a 20 % reduction of local precursor emissions, indicating the strong role of long-range transport in influencing surface ozone in the Northern Hemisphere. Results derived from Phase 1 of the TF-HTAP exercise (and the Phase 2 exercise described by Galmarini et al., 2017) are discussed in more detail by Fiore et al. (2009), Reidmiller et al. (2009) and Huang et al. (2017); Jonson et al. (2018).
An alternative approach to the perturbation technique for source attribution is “tagging” (e.g. Wang et al., 1998; Dunker et al., 2002; Grewe et al., 2010, 2017; Emmons et al., 2012; Derwent et al., 2015; Butler et al., 2018; Bates and Jacob, 2020). When applied to ozone source attribution in a CTM, this technique involves labelling (or “tagging”) modelled ozone with the identity of either the geographical region in which it is chemically produced or with the identity of the emitted precursor(s) that ultimately led to its production. A common challenge faced by all tagging approaches is that the production of one molecule of ozone in the troposphere requires two precursors: one molecule of NO and one peroxy radical (produced during reactive carbon oxidation). This then poses the question of whether the ozone molecule should inherit its tag from the emitted NOx, from the emitted reactive carbon, or in some other way? Butler et al. (2018) provides a detailed review of the different approaches to answering this question, including the trade-offs made in each case. Butler et al. (2018) also describe a novel and unique tagging methodology that allows for the separate attribution of tropospheric ozone to both its NOx and its reactive carbon precursor, at the cost of extra computational expense compared with other tagging methodologies. Recent work from Bates and Jacob (2020) takes the approach of defining an extended odd oxygen family including peroxy radicals, which effectively shifts the production of odd oxygen purely to photolysis reactions. Further comparison of different tagging approaches is beyond the scope of this paper, but it remains an interesting topic for future work.
Tagging and perturbation are complementary approaches (Clappier et al., 2017; Thunis et al., 2019; Mertens et al., 2020). While tagging delivers information about the contribution of different emission sources to a pollutant of interest, perturbation studies deliver information about the sensitivity of pollutants to changes in emissions, including changes in the chemical lifetimes of pollutants in response to variations in emissions. In the absence of non-linear chemical interactions, these two different approaches ultimately yield the same results; however, for tropospheric ozone, which can show highly non-linear interactions between its NOx and reactive carbon precursors under some circumstances, these approaches can sometimes yield very different results (Grewe et al., 2010; Mertens et al., 2018). As air pollution mitigation strategies must involve some change in emissions, perturbation studies will always be necessary for policy-relevant modelling of atmospheric chemistry. However, tagging studies on their own can play a role in helping to identify which emissions to mitigate (Grewe et al., 2010). When combined with perturbation studies, tagging can reveal how the contribution of unmitigated sources to ozone changes in response to mitigation measures (Mertens et al., 2018). Butler et al. (2018) also noted that tagging studies provide useful diagnostic information about model processes and argued for their inclusion in model intercomparison exercises. The method described by Butler et al. (2018) is currently the only available approach that provides the separate attribution of tropospheric ozone to its NOx and reactive carbon precursors. Other schemes take different approaches to the attribution of ozone to these two chemically distinct precursors. A thorough review of several different approaches is presented in Butler et al. (2018).
In this study, we use the ozone tagging methodology previously described by Butler et al. (2018) to perform a source attribution for ground-level ozone to both NOx and reactive carbon. This work builds on the work of Butler et al. (2018) by tagging anthropogenic emissions of NOx and reactive carbon by their geographical source region as well as examining the seasonal cycle of the surface ozone attribution in receptor regions as defined in the HTAP Phase 2 exercise. By performing the separate attribution of ground-level ozone to both NOx and reactive carbon (including methane), we hope to provide more useful information to inform emission mitigation scenarios. We also show how our tagging methodology can be used as a model diagnostic tool to understand the atmospheric budgets of ozone and PAN in more detail than previously possible, potentially informing efforts to reduce the currently high level of inter-model uncertainty. Furthermore, we examine the changing contributions of the different sources of NOx and reactive carbon to a perturbation of the global methane burden, showing how the contribution of emissions from unmitigated sectors would respond to the mitigation of methane emissions.
The tagging approach and model set-up is described in Sect. 2. In Sect. 3, we evaluate our simulations against observations from the Tropospheric Ozone Assessment Report (TOAR) and the ensemble of simulations from HTAP Phase 2, show the intercontinental source attribution for ozone and its precursors, and examine the response of this source attribution to a 20 % perturbation in the global methane burden. Conclusions are drawn in Sect. 4.
Simulations are performed with CAM4-chem (Community Atmosphere Model version 4 with chemistry), which is a component of the CESM (Community Earth System Model) version 1.2.2 (Tilmes et al., 2015; Lamarque et al., 2012), using the same model configuration as Butler et al. (2018). The model is run at a horizontal resolution of 1.9∘ × 2.5∘, with 56 vertical levels using specified dynamics for the year 2010 from the MERRA reanalysis (Rienecker et al., 2011). As in Butler et al. (2018), we have replaced the default chemical mechanism with a tagged mechanism based on an earlier version of the MOZART-4 mechanism (Emmons et al., 2012). Our tagging system allows for the attribution of tropospheric ozone to chemical production by either NOx or reactive carbon precursors (as well as transport from the stratosphere). A complete attribution of tropospheric ozone to both kinds of precursors requires two model runs: one with NOx emissions tagged and another with reactive carbon emissions tagged. Chemical production of ozone in the stratosphere (primarily through the photolysis of molecular oxygen) and other minor production pathways for tropospheric ozone are also tagged, as described in Butler et al. (2018). For both NOx tagging and VOC tagging, the sum of the tagged ozone tracers is equal to the total ozone as simulated by the model.
As in Butler et al. (2018), anthropogenic emissions of NOx, CO, and NMVOCs for 2010 are taken from the EDGAR-HTAPv2 emission inventory (Janssens-Maenhout et al., 2015), biomass burning emissions are from GFEDv3 (van der Werf et al., 2010), and methane is held fixed at the surface to a global average value of 1760 ppb, as in Tilmes et al. (2015). Two simulations (base runs) are performed with this model set-up: one in which all sources of NOx are tagged as described below (the “NOx-tagged” run) and one in which all sources of reactive carbon are tagged as described below (the “VOC-tagged” run). As in Butler et al. (2018), the length of the spin-up period was 1 year for the NOx-tagged run, and it was 2 years for the VOC-tagged run. The model was deemed to be spun up when the maximum difference between the simulated December mean surface ozone attributable to any tagged source was less than 1 % in any 2 subsequent years of simulation.
With the exception of surface-based anthropogenic emissions of NOx, CO, and NMVOCs, the tag identities used in this study are identical to those used in Butler et al. (2018). In this study, all surface-based anthropogenic emissions are tagged with a label representing the geographical location at which the emissions occur. This approach allows for the attribution of simulated ozone to anthropogenic precursor emissions from specific locations. Specifically, anthropogenic emissions of NOx and reactive carbon are tagged according to their Tier 1 source region as defined for the HTAP Phase 2 multi-model ensemble experiment, which is described in more detail in (Galmarini et al., 2017). Due to computational constraints, not all of the HTAP Tier 1 regions are tagged in this study. As the primary focus of this study is on the attribution of ground-level ozone in the Northern Hemisphere, only the major anthropogenic Northern Hemisphere source regions are tagged, whereas other anthropogenic sources are tagged with the label “Rest of the world”. A full list of the tags used in the NOx- and VOC-tagged runs is given in Table 1. The explicitly tagged source regions differ between the NOx-tagged and VOC-tagged runs because VOC tagging is computationally more expensive than the NOx tagging (Butler et al., 2018). One important difference between this study and Butler et al. (2018) is that anthropogenic emissions of CO for each source region are tagged together with emissions of NMVOCs in this study in order to save computational resources. For the emissions tagged as “Oceanic emissions” (Table 1), we note that the only source of NOx from this region in our simulations is from shipping and that the major source of reactive carbon is biogenic emissions of dimethyl sulfide (DMS).
1 NOx from “Oceans” is exclusively from shipping, whereas reactive carbon from this region is predominantly biogenic. 2 For NOx tagging, the stratosphere tag is applied directly to ozone produced in the stratosphere (as for VOC tagging) and also to NO produced from the dissociation of N2O. 3 For NOx tagging, “Extra production” of ozone is due to the self reaction of OH radicals and reactions between HO2 and organic peroxy radicals. 4 For VOC tagging, “Extra production” of ozone is due to the self reaction of OH radicals and reaction of OH with H2O2.
In addition to the NOx- and VOC-tagged base runs described above, we also perform two additional runs in order to investigate the response of tropospheric ozone to a perturbation in the tropospheric burden of methane: one with NOx tagging and another with VOC tagging. In each of these methane perturbation runs, the initial atmospheric methane burden and the methane mixing ratio imposed at the surface as a boundary condition are reduced by 20 %. This translates to a surface methane mixing ratio of 1410 ppb in these methane perturbation runs. In the above-mentioned perturbation runs, all other sources of NOx and reactive carbon are left unchanged. The methane perturbation runs also require 2 years of spin-up for the model to arrive at steady state.
CAM4-chem in version 1.2.2 of the CESM has previously been evaluated by Tilmes et al. (2015), and the modified version used in this study has also been discussed thoroughly by Butler et al. (2018). In Sect. 3, we describe the key differences in methane and tropospheric ozone between our base simulation and the CAM4-chem simulation reported by Tilmes et al. (2015), and we compare our simulated surface ozone with observations from TOAR (Schultz et al., 2017) as well as with the ensemble of CTM simulations from the HTAP Phase 2 multi-model study (Galmarini et al., 2017). The full set of CTMs that participated in the HTAP Phase 2 multi-model ensemble is given in Table 3 of Galmarini et al. (2017). In this study we compare surface ozone from our base simulation with results from a subset of 12 CTMs: CAM-chem (simulations performed by the National Center for Atmospheric Research, NCAR), CHASER_re1, CHASER_t106, C-IFS, C-IFS_v2, EMEP_rv4.5, EMEP_rv48, GEM-MACH, GEOS-Chem Adjoint, GEOS-Chem, OsloCTM3.v2, and RAQMS. Details of the configurations used by each of these models in the HTAP Phase 2 ensemble can be found in Galmarini et al. (2017) and references therein.
All results presented in this study are based on the definition of the troposphere as the model grid cells below the level of 150 ppb of ozone. By design, the ozone simulated in our base model runs is identical to the simulation reported in Butler et al. (2018). Our simulation for 2010 produces a tropospheric ozone burden of 319 Tg(O3), which is within 1 standard deviation of the multi-model mean reported by Young et al. (2013) for the year 2000 (337±23 Tg(O3)). Our simulated tropospheric ozone burden is slightly higher than the burden reported by Tilmes et al. (2015) using a similar model set-up (309 Tg(O3)), which could be due to the use of different emissions datasets. Our simulated tropospheric methane burden (4150 Tg(CH4)) is the same as that reported by Tilmes et al. (2015), but our methane lifetime (due to oxidation in the troposphere by OH), at 7.59 years, is shorter than the 8.82 years reported by Tilmes et al. (2015), which is also likely due to the use of different emission datasets. Our methane lifetime is towards the lower end of the range (7.1–10.6 years) simulated in CTMs, as reported by Saunois et al. (2016).
In Fig. 1, we compare our simulated monthly mean surface ozone mixing ratio for 2010 with data from TOAR and with the other models in the HTAP Phase 2 CTM ensemble. Results are shown averaged over HTAP Tier 2 receptor regions, and they only include grid cells for which TOAR observations are available. In general, most of the HTAP models overestimate the monthly mean surface ozone mixing ratio in regions where observations are available, which is consistent with the high model bias reported by Young et al. (2018). Also apparent from Fig. 1 is the large range in simulated surface ozone between members of the HTAP model ensemble, which is especially high in the northern spring, approaching a spread of approximately 30 ppb between the lowest and highest ensemble members, which is of a similar order to the Northern Hemisphere annual mean surface ozone itself. Our modelled monthly average surface ozone mixing ratio in the HTAP Tier 2 receptor regions is generally close to the HTAP ensemble mean and is usually within 1 standard deviation of the ensemble mean.
3.1 Source attribution of tropospheric ozone
The attribution of annual average tropospheric ozone to emissions of NOx and reactive carbon precursors based on the source tags from Table 1 is shown in Figs. 2 and 3 and is quantified in Tables 2 and 3. Figures showing the attribution of monthly mean ozone are available in the Supplement. For each tagged source, Tables 2 and 3 include the respective emissions of NOx and reactive carbon and, where applicable, the contribution of each source to the 2010 average tropospheric ozone burden, the contribution to the Northern Hemisphere 2010 annual average surface mixing ratio, and the ozone production efficiency of each emission source (defined here as the contribution of each emission source to the tropospheric ozone burden, in units of moles of ozone per mole of N or C emitted). Figures 2 and 3 show the spatial distribution of the annual average surface ozone as attributed to each source of NOx and reactive carbon, respectively. In Fig. 2, ozone attributable to anthropogenic NOx emissions in some source regions (specifically South East Asia, northern Africa, the Middle East, Middle America, and Central Asia) has been added to the “Rest of the world” total in order to unify the definition of this source region with the definition of this region in the VOC-tagged run. The difference in the stratospheric contribution between the NOx- and VOC-tagged runs is due to the role of NOx produced in the stratosphere from the dissociation of N2O. Ozone produced in reactions involving this stratospheric source of NOx are counted in our source attribution as stratospheric ozone, as described in Butler et al. (2018).
3.1.1 Attribution to NOx emissions
Anthropogenic NOx emissions from the three major high-latitude source regions (Europe, East Asia, and North America) contribute to high modelled ozone concentrations both locally and in the Northern Hemisphere background. Lightning NOx, soil NOx, and ozone input from the stratosphere all contribute additionally to modelled global background ozone. Emissions of NOx from shipping contribute significantly to ozone over the major Northern Hemisphere ocean basins, which is also transported over continental regions. South Asia stands out in comparison with the other major Northern Hemisphere source regions in that ozone produced from NOx emitted in South Asia is relatively localized to the South Asian region itself and is not transported into the hemispheric background to the same extent as ozone produced from NOx emissions in the other major Northern Hemisphere source regions.
Table 2 shows that NOx emissions from lightning and aircraft are especially efficient at producing ozone in the free troposphere, consistent with previous work (e.g. Beck et al., 1992; Jacob et al., 1996; Dahlmann et al., 2011). Similarly, surface emissions of NOx from regions closer to the tropics (e.g. South East Asia and Middle America) produce ozone more effectively due to the rapid convective transport of emitted NOx into the free troposphere, as shown by Zhang et al. (2016). Of the major Northern Hemisphere source regions, NOx emissions from South Asia are the most efficient at producing ozone, consistent with a stronger role of vertical transport over this region. In contrast, NOx emissions from the major anthropogenic source regions in the high northern latitudes (Europe, East Asia, and North America) are among the least productive of all global NOx emissions, consistent with a relatively small amount of convective transport, leading to higher rates of NOx removal. Despite their low ozone production efficiency, emissions of NOx in the high northern latitudes contribute significantly to surface ozone across the Northern Hemisphere (Fig. 2, Table 2).
Table 2 also shows that NOx emissions from shipping are also relatively efficient at producing ozone, which is also consistent with previous work (e.g. Lawrence and Crutzen, 1999; Hoor et al., 2009). The high ozone production efficiency of ship emissions is due to their location in relatively pristine regions with few other sources of NOx. Owing to the high ozone productivity of ship emissions and the fact that they are emitted at relatively high latitudes, they contribute significantly to the Northern Hemispheric background (Fig. 2, Table 2). As noted above, the ozone production from ship NOx is likely to be overestimated due to the artificial dilution of emissions into relatively coarse model grid cells.
Mertens et al. (2018) reported a contribution of shipping to the tropospheric ozone burden of 18 Tg(O3) using their tagging technique and based on a model simulation with ship NOx emissions of 6 Tg(N) yr−1. In our study, we calculate a contribution of ship NOx to tropospheric ozone of 19.9 Tg(O3) based on ship NOx emissions of 4.28 Tg(N) yr−1, implying a much higher ozone production efficiency for ship NOx. As the tagging technique used by Mertens et al. (2018) is based on the method described by Grewe et al. (2017), which combines the effects of tagged NOx and reactive carbon precursors into a single tagged ozone molecule during ozone production, we do not expect our results to be directly comparable. As shipping emits significantly more NOx than reactive carbon, we would expect the combinatorial tagging approach of Mertens et al. (2018) to attribute less ozone to shipping than our method, because the ozone produced from ship NOx would also be partially attributed to the reactive carbon precursor involved in the ozone production. Indeed, Mertens et al. (2018) reported maximum contributions of shipping to surface ozone of about 10 ppb in summer over major Northern Hemisphere ocean basins. In our study, surface ozone attributable to shipping over these regions can exceed 20 ppb (see the Supplement).
3.1.2 Attribution to reactive carbon emissions
Methane and biogenic emissions clearly stand out as major reactive carbon precursors to tropospheric ozone, contributing 35 % and 24 % to the tropospheric ozone burden in our simulation, respectively. Anthropogenic emissions of reactive carbon (excluding biomass burning) together contribute about 14 % to the tropospheric ozone burden. The relatively low influence of anthropogenic reactive carbon emissions on ground-level ozone has been noted elsewhere (e.g. HTAP, 2010; Butler et al., 2018); however, despite this low overall ozone productivity, anthropogenic reactive carbon emissions from source regions in higher northern latitudes still have disproportionately high contributions to surface ozone in the Northern Hemisphere (Table 3).
Due to the emissions of CO being tagged together with emitted VOCs in this study, the contribution of each tagged source to the tropospheric ozone burden (and, therefore, also the ozone production efficiency of each tagged source) is a mixture of ozone production due to emitted CO and emitted NMVOCs. The ozone attributed to methane oxidation in Table 3 is due to all stages of methane oxidation in the MOZART-4 chemical mechanism, including the final step in which CO from earlier stages of methane oxidation is itself oxidized to CO2. The oxidation of CO can produce a maximum of one peroxy radical (HO2). Therefore, the maximum ozone production potential of CO is 1 mol of ozone per mole of emitted CO. VOCs (including methane) can produce significantly more ozone per mole of emitted carbon, when the subsequent oxidation of the initial oxidation products is taken into account (Bowman and Seinfeld, 1994; Atkinson, 2000; Butler et al., 2011). Future studies using this tagging methodology should consider tagging CO emissions separately from NMVOC emissions if they aim to determine the ozone production efficiency of anthropogenic NMVOC emissions from different world regions. Butler et al. (2018) tagged NMVOC emissions separately from CO emissions, but they did not tag anthropogenic emissions separately according to their geographical region. We re-examined the output of the otherwise identical VOC-tagged run described by Butler et al. (2018) in order to determine the ozone production efficiency of NMVOC emissions from anthropogenic, biomass burning, and biogenic sources: they are 0.0580, 0.0354, and 0.0268 mol(O3) mol (C)−1 for anthropogenic, biomass burning, and biogenic sources, respectively. The ozone production efficiency of biogenic NMVOCs recalculated from Butler et al. (2018) is not significantly different from the value reported here in Table 3, reflecting the relatively minor contribution of CO to the total amount of emitted biogenic reactive carbon. For biomass burning and anthropogenic sources, however, the ozone production efficiency of NMVOCs emitted from these sources is greater than the corresponding value from Table 3, reflecting the fact that the numbers from Table 3 also include emissions of CO.
We note that methane has a higher ozone production efficiency per unit of reactive carbon (0.0689 mol mol (C)−1; Table 3) than any of the NMVOCs in our runs. The low ozone production efficiency of biogenic NMVOCs is consistent with large amounts of isoprene being emitted in remote regions under low-NOx conditions, where loss of peroxy radicals through reaction with other peroxy radicals could be expected to dominate (Atkinson, 2000). It might, however, be expected that anthropogenic NMVOCs would have a higher ozone production efficiency, due to the fact that they are co-emitted with anthropogenic NOx, favouring the conversion of NO to NO2 via reaction with peroxy radicals and, thus, the production of ozone. The relatively low production efficiency of anthropogenic NMVOCs in our model runs could be due to the fact that the relatively simple chemistry of methane oxidation is well described in the version of the MOZART-4 chemical mechanism used here, in which the relatively complex chemistry of the higher NMVOCs has been simplified. Coates and Butler (2015) noted that the ozone production potential of NMVOCs in simplified chemical mechanisms tended to be lower than the more comprehensive Master Chemical Mechanism (Saunders et al., 2003). Utembe et al. (2010) previously noted increased tropospheric ozone in a CTM when using a more explicit oxidation mechanism for NMVOCs. The extremely low ozone production efficiency of reactive carbon from oceanic sources in Table 3 is due to the lack of any ozone-forming pathways in the oxidation of dimethyl sulfide (DMS) in the MOZART-4 chemical mechanism as used in this study. DMS is the dominant source of reactive carbon over the oceans in our model simulations.
To our knowledge, the only other study to perform source attribution of global tropospheric ozone specifically to reactive carbon precursors is Butler et al. (2018), on which the present study builds. Here, we attribute 113 Tg(O3) to methane oxidation (Table 3). Grewe et al. (2017) attribute 45 Tg(O3) to methane using their tagging approach, which combines the effects of tagged NOx and reactive carbon precursors into a single tagged ozone molecule during ozone production. Ozone production due to methane oxidation under their combinatorial tagging approach would be expected to also include the attribution to the source of NOx involved in the ozone production. Thus, we would expect Grewe et al. (2017) to attribute approximately half of the amount of ozone to methane that we would. Doubling the value reported by Grewe et al. (2017) yields 90 Tg(O3) attributable to methane oxidation, which is much closer to our value of 113 Tg.
Widespread implementation of tagging techniques for the separate tagging of NOx and reactive carbon emissions in other CTMs as well as systematic intercomparisons of their results could help to understand differences in the simulated budgets of tropospheric ozone. Due to the variety of approaches taken by different tagging techniques (summarized in Butler et al., 2018), an intercomparison of a range of different tagging techniques with other methods for ozone source attribution could be informative.
3.2 Source–receptor relationships for ozone
3.2.1 Annual average surface ozone
Figure 4 shows the modelled annual average surface ozone concentration in the five major HTAP Tier 1 regions in the Northern Hemisphere (Europe; Russia, Belarus, and the Ukraine; South Asia; East Asia; and North America), including a full attribution of ozone in each of these regions to all sources, such as transport from the stratosphere and emitted precursors of both NOx and reactive carbon. Annual average ozone in most of the regions shown in Fig. 4 is close to the Northern Hemisphere annual average of 30 ppb (Table 2), except in South Asia and East Asia where the annual average surface ozone mixing ratio is closer to 40 ppb. The difference in each case is primarily due to a larger source of ozone produced from locally emitted precursors. Transport from the stratosphere contributes approximately 2–4 ppb to annual average surface ozone depending on the receptor region, which is consistent with the 2.91 ppb contribution of stratospheric ozone to the annual average surface ozone in the Northern Hemisphere average surface ozone (Table 3). As already shown in the previous section, anthropogenic sources of NOx dominate other NOx sources as ozone precursors, while the major reactive carbon precursors are methane and biogenic volatile organic compounds (BVOCs).
In each of the five regions shown in Fig. 4, natural sources and long-range transport of ozone produced from extra-regional anthropogenic precursors together contribute more to the annual average surface ozone than anthropogenic emissions within the region itself. In each region, the local anthropogenic NOx emissions produce more ozone than can be attributed to anthropogenic NOx emissions in any other Tier 1 regions; moreover, with the exception of (South Asia), the combined contribution of external anthropogenic NOx emissions to annual average surface ozone is greater than the local contribution. The importance of long-range transboundary transport of ozone has also been noted elsewhere (HTAP, 2010).
While anthropogenic precursor emissions from South Asia contribute significantly to surface ozone within the South Asian region, they contribute relatively little to surface ozone in the other four regions shown in Fig. 4. This is also consistent with the surface ozone maps in Fig. 2 and the higher ozone productivity of NOx emissions from South Asia when compared with the other major Northern Hemisphere source regions (Table 2). Emissions of ozone precursors (particularly NOx) from South Asia are transported efficiently into the free troposphere, where they contribute disproportionately to the global tropospheric ozone burden (as also noted by Zhang et al., 2016), but the contribution of South Asian emissions to surface ozone in other parts of the Northern Hemisphere is disproportionately smaller than emissions from the other HTAP Tier 1 regions.
3.2.2 Seasonal cycles of surface ozone
Figures 5, 6, and 7 show the seasonal cycles of surface ozone in the three selected Tier 2 regions “North West Europe”, ”North East China”, and “North West United States” respectively. These regions are selected in order to compare two regions on the western side of their respective continents (where long-range transport is expected to be important) with a region on the eastern side of its continent that is also a major source region. A set of figures for other HTAP Tier 2 regions with a complete attribution of surface ozone to all tagged HTAP Tier 1 source regions is available in the Supplement. In Figs. 5, 6, and 7, results are shown for both NOx and VOC tagging. In each receptor region, the contribution of long-range transport due to extra-regional anthropogenic emissions from HTAP Tier 1 regions is shown both in aggregate (panels a and c) and by individual Tier 1 source region (panels b and d). The definition of the “Rest of the world” ozone tracer has been harmonized between the NOx- and VOC-tagged runs in these figures. Consistent with the annual averages from Fig. 4, anthropogenic NOx sources also dominate the seasonal cycle of modelled ozone, while the major reactive carbon precursors of ozone are methane and BVOCs.
All three receptor regions show a seasonal cycle of ozone with a spring–summer ozone maximum superimposed on a year-round ozone baseline. The summertime maximum in ozone is clearly due to local photochemical production from the combination of locally emitted anthropogenic NOx and BVOCs. The strong role of locally emitted precursors in the production of ozone in summer is consistent with earlier work (e.g. Reidmiller et al., 2009; Huang et al., 2017; Jonson et al., 2018; Han et al., 2019), while the importance of BVOC emissions, especially isoprene, for ozone production in summer has also been noted elsewhere (e.g. Chameides et al., 1992; Andersson and Engardt, 2010; Han et al., 2019). Biogenic emissions of NOx (from soils) also contribute to this summertime maximum in the local photochemical ozone production in all three of the regions shown in Figs. 5, 6, and 7, although to a much smaller extent than anthropogenic NOx emissions.
The year-round baseline ozone in our model simulations in all three receptor regions can be primarily explained by slower photochemistry involving methane as the reactive carbon precursor, in combination with extra-regional anthropogenic NOx (Figs. 5, 6, and 7). The contribution of methane to surface ozone is slightly higher in summer, coinciding with the peak in local anthropogenic NOx emissions, which is consistent with local photochemical ozone production from enhanced local methane oxidation. The contribution of extra-regional anthropogenic NOx to surface ozone is largest in spring, coinciding with the peak in the contribution of extra-regional anthropogenic reactive carbon, which is consistent with long-range transboundary transport of ozone produced elsewhere.
Maxima in springtime ozone have previously been linked to long-range transboundary transport in all major receptor regions (HTAP, 2010; Lin et al., 2012; Jonson et al., 2018; Ni et al., 2018). This transported ozone can be attributed to input from the stratosphere as well as extra-regional anthropogenic emissions of NOx and reactive carbon. In our simulations, the contribution of stratospheric ozone peaks around March, whereas the contribution of extra-regional anthropogenic emissions tends to peak around April, when it contributes more strongly to the monthly average surface ozone in each region than local anthropogenic NOx (Figs. 5, 6, and 7). In all Northern Hemisphere regions, the springtime peak in the contribution of extra-regional anthropogenic reactive carbon is smaller than the corresponding springtime peak in the contribution of extra-regional anthropogenic NOx. Previous work has identified uncertainties in the treatment of ozone production from NMVOC oxidation as a potential source of inter-model differences (e.g. Emmerson and Evans, 2009; Utembe et al., 2010; Coates and Butler, 2015). The relatively large influence of anthropogenic NMVOCs on springtime ozone (compared with their influence during other times of the year) could be a contributing factor to the large spread in springtime ozone simulated by current-generation CTMs (Fig. 1).
NOx from shipping is the largest single contributor to springtime transboundary ozone transport in all three receptor regions shown here. We note, however, that the coarse resolution of our model (2∘) would be expected to exaggerate the effects of ship NOx on ozone production due to rapid dilution of the emissions (von Glasow et al., 2003) and exaggerate the transport of NOx and ozone near coastlines due to unrealistically high diffusion between adjacent land and ocean grid cells. Thus, the contribution of shipping emissions to surface ozone in our simulations should be considered an upper bound, especially in coastal regions. The high contribution of ship NOx to summertime ozone in Europe (Fig. 5) may be an artefact of the coarse model resolution combined with the high shipping volume near coasts in the eastern North Atlantic Ocean as well as the North Sea. Early work by Lawrence and Crutzen (1999) showed a stronger influence of ship NOx on surface ozone in northwestern Europe than on any other continental region. Jonson et al. (2018) showed that the only other CTM in the HTAP Phase 2 ensemble to report the results of a perturbation of shipping emissions (the EMEP_rv48 CTM with a resolution of 0.5∘ × 0.5∘) shows a similar magnitude for the influence of ship NOx on summertime ozone in Europe as for springtime ozone. Using a regional model at 50 km ×50 km resolution and a similar ozone tagging system as that used in the present study, Lupaşcu and Butler (2019) showed that the contribution of ship NOx to ozone in coastal regions of Europe reaches a maximum level in summer. Jonson et al. (2020) used a global model with a resolution of 0.5∘×0.5∘ to show that shipping near coastal areas strongly influences ozone over northwestern Europe in both spring and summer, while NOx emissions from shipping on the high seas have a stronger influence on European ozone in spring than in summer. In contrast, using a higher-resolution regional model (20 km ×20 km) for Europe, Aksoyoglu et al. (2016) did not indicate that ship NOx played such a strong role in summertime ozone over Europe.
In other receptor regions, the influence of ship NOx emissions on surface ozone is largest in spring (Figs. 6, 7), suggesting a stronger influence of NOx emissions over the high seas on springtime ozone in our simulations in these regions. Model-dependent inconsistencies in the treatment of ship NOx emissions may play a role in the large spread of simulated ozone between models in springtime (Fig. 1). Future work should examine the contribution of ship NOx emissions to ozone in both spring and summer, using model systems that include the better representations of plume dilution over the major shipping routes, more refined attribution to shipping emissions from coastal regions and the high seas, and higher resolution over receptor regions.
Previous work has indicated a strong influence of anthropogenic emissions from both North America and East Asia on springtime ozone in Europe (Jonson et al., 2018), a strong influence of East Asian emissions on springtime ozone in North America (Lin et al., 2012), and a diverse range of intercontinental influences on springtime ozone in East Asia (Ni et al., 2018). The direct numerical comparison of our results with these previous studies is difficult due to the different methodologies used. These previous studies all employed the perturbation technique to determine the influence of all anthropogenic emissions (both NOx and reactive carbon) from each source region, whereas this work uses a tagging approach that separately attributes ozone to emitted NOx and reactive carbon. However, qualitatively, our results (as shown in Figs. 5, 6, and 7) do appear to be consistent with this earlier work. Comparison of the NOx-tagged and VOC-tagged results in these figures shows that anthropogenic NOx emissions from most source regions have a stronger influence on springtime ozone in any given receptor region than anthropogenic emissions of reactive carbon. The only exception to this is East Asia, where reactive carbon emissions are substantially higher than in other HTAP Tier 1 source regions (Table 3). Reactive carbon emissions from East Asia contribute approximately equally to springtime ozone in North America when compared to the East Asian NOx emissions, whereas East Asian reactive carbon contributes more to springtime ozone than East Asian NOx in Europe.
3.3 Long-range transport of ozone precursors
Fiore et al. (2018) suggested that measurements of the abundance of PAN at mountaintop sites in spring may be useful as an indicator of the intercontinental transport of ozone and its precursors, as well as being a diagnostic for uncertainties in CTM simulations, which show large inter-model differences in simulated PAN (Emmons et al., 2015). The column integrated density of PAN in the lower troposphere (defined here as the model layers between 500 and 800 hPa) in the three HTAP Tier 2 receptor regions “North West Europe”, ”North East China”, and “North West United States” are shown in Figs. 8, 9, 10 respectively. A set of figures for other HTAP Tier 2 regions, with the complete attribution of surface ozone to all tagged HTAP Tier 1 source regions is available in the Supplement. Simulated PAN is highest in late winter to early spring, which is consistent with earlier work (Fischer et al., 2014; Fiore et al., 2018). The extra-regional contribution to PAN is also highest in spring; this is primarily due to anthropogenic NMVOCs, which is also consistent with Fischer et al. (2014).
Our model simulations with NOx and VOC tagging provide a unique opportunity to examine the origin and fate of PAN as simulated in our model, as this allows the simultaneous attribution of simulated PAN to both its NOx precursor and its reactive carbon precursor. Comparison of panels b and d in Figs. 8, 9, and 10 consistently shows that, for any given land-based HTAP Tier 1 source region, the anthropogenic NMVOC emissions contribute more to PAN formation than the anthropogenic NOx emissions from that region to the PAN modelled in all HTAP Tier 2 receptor regions. In all cases, the balance of extra-regional PAN is due to NOx emissions from shipping. In our simulations, significant amounts of PAN are formed downwind of the regions in which the anthropogenic NMVOC precursors are emitted, often via reaction with NOx emitted from shipping. A strong influence of anthropogenic NOx emissions on PAN in the northern midlatitudes is consistent with the results shown by Fischer et al. (2014; their Fig. 7, which does not distinguish between different sources of anthropogenic NOx).
Figures 8, 9, and 10 also show that the reactive carbon component of PAN is generally more persistent than the NOx component. For example, the contribution of anthropogenic NMVOCs from North America to springtime PAN over East Asia is only slightly lower than its contribution to springtime PAN over Europe, which is much closer to North America considering the prevailing westerly winds (panel d of Figs. 8 and 9). In contrast, the contribution of anthropogenic NOx from North America to springtime PAN in East Asia is substantially less than its contribution to springtime PAN over Europe (panels b of Figs. 8 and 9).
3.4 Attribution of Northern Hemisphere total organic reactivity
We examine the Northern Hemisphere budget of reactive carbon in more detail in Fig. 11. This figure shows the seasonal cycle of the Northern Hemisphere column-integrated total reactivity with respect to the OH radical of all reactive carbon-containing species in our simulation attributed to their emission source. The total OH reactivity of reactive carbon species of an air mass is often linked to its ozone production potential (Chameides et al., 1992; Kleinman et al., 2002). In each case, the OH reactivities shown in Fig. 11 include the OH reactivity of the primary emitted species as well as the OH reactivity of each carbon-containing oxidation product. These were calculated using the monthly averaged output of the modelled concentration of each carbon-containing species (including its associated tags) and the temperature- and pressure-dependent rate coefficients for their reaction with the OH radical; they were then averaged over all Northern Hemisphere grid cells and weighted by air density.
The total Northern Hemisphere OH reactivity of reactive carbon remains fairly constant year-round at about 0.6–0.7 s−1, but the seasonal cycles of the OH reactivity attributable to different reactive carbon sources show more variability. Methane (and its oxidation products) contribute about 0.2–0.3 s−1 (almost half of the total hemispheric reactivity), with a slight maximum in the summer, which is consistent with enhanced oxidation (and, thus, the enhanced availability of more reactive methane oxidation products) due to higher OH in summer. The contributions of anthropogenic and biogenic reactive carbon sources to total hemispheric reactivity are similar, ranging between about 0.1 and 0.3 s−1, but with distinct seasonal cycles. The reactivity of biogenic carbon is highest in summer–autumn (which is consistent with the Northern Hemisphere growing season), whereas the reactivity of anthropogenic carbon is highest in winter–spring (which is consistent with constant year-round anthropogenic emissions and a build-up of reactive carbon over winter due to lower hemispheric OH). The build-up of anthropogenic reactive carbon throughout the Northern Hemisphere over winter combined with the resumption of OH chemistry in spring is consistent with the disproportionate effect of extra-regional anthropogenic reactive carbon on springtime ozone seen in Figs. 5, 6, and 7. Thus, uncertainties in the model chemical mechanisms associated with the oxidation of anthropogenic NMVOCs (e.g. Emmerson and Evans, 2009; Utembe et al., 2010; Coates and Butler, 2015) may also contribute to the large spread in simulated ozone seen in the HTAP ensemble during spring (Fig. 1).
Figure 11 also shows the geographical origin of Northern Hemisphere anthropogenic carbon reactivity. Emissions of reactive carbon from East Asia stand out as the single major source of enhanced anthropogenic carbon reactivity in winter and spring in our simulations. This is consistent with the high emissions of reactive carbon from this region in 2010 noted earlier (Table 3). Growth in NMVOC emissions from East Asia may have continued since this time (Li et al., 2019), while NOx emissions have been decreasing (Liu et al., 2017). Increasing trends in the local production of ozone during summer over East Asia (e.g. Li et al., 2018) should be associated with the increased oxidation of reactive carbon and, thus, potentially less export of reactive carbon into the Northern Hemisphere background during summer. We expect, however, that increasing emissions of reactive carbon in East Asia could lead to an increased build-up of East Asian reactive carbon in the Northern Hemisphere over winter and, in turn, to an increased East Asian contribution to extra-regional springtime ozone in other parts of the Northern Hemisphere.
Our tagging technique is currently the only one we know of that is capable of examining the budget of reactive carbon in the level of detail presented in this study. The separate tracking of the carbon-containing and nitrogen-containing components of PAN is particularly informative, suggesting that significant amounts of PAN are formed downwind of source regions in our model, especially during winter and spring, due to a build-up of anthropogenic reactive carbon over winter when photochemistry is relatively slow. Given the large variety in model representations of NMVOC chemistry, including PAN formation and decomposition processes (Emmerson and Evans, 2009; Knote et al., 2015) and the large inter-model differences in simulated PAN (Emmons et al., 2015), the widespread implementation of similar tagging diagnostics in other CTMs may help to provide additional information about the origin and fate of simulated PAN and, more generally, about the influence of reactive carbon on atmospheric composition. In combination with routine mountaintop observations of springtime PAN, this may aid understanding of the global PAN budget (Fiore et al., 2018) and other processes responsible for the intercontinental transport of air pollution. Better constraints on these chemical and transport processes should also help to reduce inter-model differences in simulated springtime ozone (Fig. 1).
3.5 Tropospheric ozone sensitivity to methane
We performed an additional set of model runs with both NOx and VOC tagging with the methane surface boundary condition reduced from 1760 to 1410 ppb, which amounts to a reduction of 350 ppb (or 20 %). This perturbation can also be expressed as an increase of 25 %. Here we interpret the methane perturbation run in terms of the atmospheric response to a 25 % increase in the methane surface mixing ratio at steady state.
In response to the 25 % increase in the imposed surface mixing ratio of methane, the total tropospheric burden increased by 776 Tg(CH4), which amounts to an increase of 23 %. The strength of the annual tropospheric chemical sink of methane due to OH increased by 72.5 Tg(CH4) (or 15.2 %). The corresponding increase in the methane lifetime was 0.48 years (or 6.75 %). The relatively small growth in the chemical methane sink compared with the magnitude of the perturbation in methane itself is consistent with the feedback of methane on its own lifetime due to the depletion of OH (Prather, 1996). Table 4 shows the response of the tropospheric ozone burden (and the contributions of different reactive carbon precursors) to the 25 % increase in the imposed surface mixing ratio of methane. The 1 ppb simulated increase in Northern Hemisphere surface ozone in response to a 25 % increase in methane burden is consistent with previous work (HTAP, 2010). The 9.22 Tg increase in tropospheric ozone burden is also consistent with the review of Fiore et al. (2008), who derived a sensitivity of 0.11–0.16 Tg(O3) per Tg(CH4) yr−1 emitted based on an analysis of the literature. We calculate 0.13 Tg(O3) per Tg(CH4) yr−1 based on our results.
The relative increase in tropospheric ozone attributed to methane (13.0 %) by our tagging scheme is comparable to, but slightly smaller than, the increase in the magnitude of the chemical methane sink due to OH (15.2 %), which is consistent with the troposphere as a whole becoming slightly more NOx limited with increasing methane. The absolute increase in the total ozone burden (9.22 Tg(O3)) is, however, significantly lower than the increase in the burden of ozone attributed to methane by our tagging scheme (13.0 Tg(O3)). When the methane burden is increased, the contribution of every other reactive carbon source to the tropospheric ozone burden decreases (each by approximately 1 %–2 %) to partially offset the increased ozone production from methane oxidation. This is also consistent with a slightly more NOx-limited atmosphere with increasing methane. In a future with an increased methane burden, control of NMVOC emissions could be expected to be less effective at the large-scale reduction of annual average ground-level ozone.
Table 5 shows the change in the contributions of different NOx sources to tropospheric ozone in response to the 25 % increase in the methane burden. As expected, all NOx sources become more productive when the total atmospheric burden of reactive carbon is increased (which is consistent with the troposphere as a whole becoming more NOx limited). The increase in the productivity of the different NOx sources under an increased burden of methane is, however, not uniform. Ozone production due to NOx from shipping stands out as highly sensitive to the global methane burden in our simulations. Ship NOx accounts for almost 30 % of the 1 ppb increase in Northern Hemisphere average surface ozone when the methane burden is increased by 25 % (Table 5), despite being a much smaller percentage of total global NOx emissions (Table 2).
The spatial distribution of the increase in annual average surface ozone from ship NOx in response to the 25 % increase in methane is similar to the spatial distribution of surface ozone due to ship NOx in our base run (Fig. 2). Figures showing the response of the attributed surface ozone are available in the Supplement. The response is largest over the major Northern Hemisphere ocean basins, but it also extends over continental regions. The seasonal cycle of the increase in annual average surface ozone from ship NOx in the three HTAP Tier 2 regions examined here in response to the 25 % increase in methane is similar to the seasonal cycle of surface ozone due to ship NOx in our base run (Figs. 5, 6, and 7). The maximum response of surface ozone from ship NOx to rising methane is simulated over the major Northern Hemisphere ocean basins in summer (which influences surface ozone in northwestern Europe in our simulations; Fig. 5), while the influence of this response over most Northern Hemisphere continental regions is generally higher in winter–spring (as seen in northeastern China; Fig. 6).
Previous work (Lawrence and Crutzen, 1999) has noted the disproportionate influence of ship NOx on tropospheric ozone due to the diffuse and widespread nature of this source over regions that would otherwise have very low NOx mixing ratios. Fiore et al. (2008) noted that the response of surface ozone to increased methane was especially strong in ship tracks. Myhre et al. (2011) also showed that ship NOx emissions reduce the global methane lifetime much more than terrestrial NOx emissions. We note again that the contribution of ship NOx to ozone in our simulations (as in most current-generation CTMs) is likely to be an overestimate due to the unrealistic dilution of these emissions into coarse model grid cells (von Glasow et al., 2003) and the lack of explicit plume chemistry (Vinken et al., 2011). However, we do expect that the interaction between ship NOx and methane for ozone production would persist in our model even with a more realistic treatment of ship emissions, as this interaction is likely due to the location rather than the magnitude of ship emissions. We are not aware of any previous work linking the combined influence of these two sources to a potentially disproportionate influence on background ozone in the Northern Hemisphere or on modelled surface ozone air quality in inhabited regions of the Northern Hemisphere, especially in spring. Given the current uncertainty in the attribution of recent trends in methane (Turner et al., 2019) and the potential for future increases in methane emissions, combined with slower reductions of NOx emissions from international shipping than from other sectors (e.g. the SSP5 future emission scenario Rao et al., 2017), we expect that model simulations of future background ozone in the Northern Hemisphere, especially during spring, may come to be increasingly influenced by ozone produced via the interaction of methane and ship NOx. Future work should investigate the ozone production via the interaction of these two sources in more detail.
We have performed a source attribution for tropospheric ozone in a chemical transport model using a novel technique that separately accounts for the influence of both the emitted NOx and the emitted reactive carbon precursors on simulated tropospheric ozone. By tagging anthropogenic emissions of NOx and reactive carbon according to their geographical region, we have calculated source–receptor relationships for the Northern Hemisphere. The results of our study are consistent with previous work, and they provide a number of important new insights that are relevant to both the mitigation of intercontinental transboundary air pollution and ongoing efforts to reduce the uncertainty in the current generation of chemical transport models.
Consistent with previous work, the annual average ground-level ozone in all major Northern Hemisphere regions is primarily influenced by extra-regional emissions of both NOx and reactive carbon. In all cases, local anthropogenic emissions of ozone precursors have a smaller influence on annual average ozone than the combined effect of precursor emissions from the rest of the world. As a reactive carbon precursor, methane contributes 35 % of the tropospheric ozone burden and 41 % of the Northern Hemisphere annual average surface mixing ratio, which is more than any other source of reactive carbon. Our novel tagging methodology also reproduces the well-known dependence of summer ozone maxima on local emissions of anthropogenic NOx and biogenic reactive carbon as well as the enhanced importance of the intercontinental transport of ozone from remote anthropogenic sources in spring. Consistent with previous work, we find that emissions of NOx at low latitudes produce free-tropospheric ozone more effectively due to more efficient vertical transport. We show, however, that NOx sources at higher northern latitudes have a stronger influence on ground-level ozone, which is known to have a lower radiative forcing but a higher influence on human health and ecosystems.
The current generation of chemical transport models has particular difficulty simulating the intercontinental transport of ozone, as shown by the large spread in ensemble simulations of ground-level ozone during the spring months. We show that our tagging methodology can deliver detailed diagnostic information about the origin and budget of springtime ozone in our model, along with information about the springtime budget of peroxyacetyl nitrate (PAN), which is also associated with springtime long-range transport and ozone production. We show that a substantial proportion of the free-tropospheric PAN simulated by our model in spring is not produced in the polluted boundary layer over the major anthropogenic source regions; instead, it is produced downwind of these regions in our model via the interaction of transported anthropogenic reactive carbon and NOx emitted from international shipping. Reactive carbon of anthropogenic origin (and its oxidation products, including PAN) builds up in our model across the entire Northern Hemisphere during the winter months and then contributes to a short burst of hemispheric-scale ozone production during spring in our simulations. In all but the most polluted source regions, anthropogenic NMVOCs do not make a significant contribution to simulated ground level ozone in any other season but spring.
In this study, we showed that the export of anthropogenic reactive carbon from East Asia may be playing a dominant role in contributing to the build-up of reactive carbon in the Northern Hemisphere over winter and, in turn, to the hemispheric-scale production of ground-level ozone in spring. Given the likely lack of recent mitigation in reactive carbon emissions from East Asia, we expect this effect to be ongoing, and we recommend that future work continue to investigate this possibility using updated emission inventories.
In addition to a contribution from the stratosphere, the springtime peak in transported ozone in our model is influenced by the interaction of two processes known to be especially poorly represented in current models: the chemistry of the intermediate oxidation products of NMVOCs and the emissions of NOx from international shipping. Furthermore, the response of ground-level ozone to changes in methane also appears highly sensitive to the treatment of ship NOx, especially in spring. We believe that our tagging technique could deliver useful information about the large differences in simulated springtime ozone between current-generation models if it were to be implemented in a larger number of models and used systematically in model intercomparison exercises. This could potentially point the way to improved representations of the processes responsible for the intercontinental transport of ozone.
Improved global CTMs are required to inform effective policies aimed at reducing the intercontinental transport of ground-level ozone – a problem which is most urgent in the springtime. In particular, we recommend that developers of emission inventories and CTMs revisit their representations of anthropogenic NMVOC emissions and the associated oxidation chemistry in order to reduce the uncertainties in the modelled springtime ozone. Additionally, more explicit representations of the NOx chemistry of ship exhaust plumes should be prioritized in order to improve the suitability of the current models for simulating both the intercontinental transport of ozone and the response of ozone to changing atmospheric methane.
The CESM is maintained by NCAR and is provided free to the community. The specific modifications made to the model to enable the tagging of ozone production have been described by Butler et al. (2018) and are available in the online supplement to that work, which is archived and accessible via the following DOI: https://doi.org/10.5194/gmd-11-2825-2018 (Butler et al., 2018).
The supplement related to this article is available online at: https://doi.org/10.5194/acp-20-10707-2020-supplement.
TB designed the study. Model runs were performed by AL. The analysis of model runs was performed by AL, AN, and TB. TB wrote the paper with input from AL and AN.
The authors declare that they have no competing interests.
The authors would like to thank Mark Lawrence, Louisa Emmons, Simone Tilmes, and Terry Keating for numerous helpful discussions during the preparation of this paper.
This work was hosted by IASS Potsdam, and financial support was provided by the Federal Ministry of Education and Research of Germany (BMBF) and the Ministry for Science, Research and Culture of the State of Brandenburg (MWFK).
This paper was edited by Frank Dentener and reviewed by Volker Grewe and one anonymous referee.
Aksoyoglu, S., Baltensperger, U., and Prévôt, A. S. H.: Contribution of ship emissions to the concentration and deposition of air pollutants in Europe, Atmos. Chem. Phys., 16, 1895–1906, https://doi.org/10.5194/acp-16-1895-2016, 2016. a
Andersson, C. and Engardt, M.: European ozone in a future climate: Importance of changes in dry deposition and isoprene emissions, J. Geophys. Res., 115, D02303, https://doi.org/10.1029/2008jd011690, 2010. a
Bates, K. H. and Jacob, D. J.: An Expanded Definition of the Odd Oxygen Family for Tropospheric Ozone Budgets: Implications for Ozone Lifetime and Stratospheric Influence, Geophys. Res. Lett., 47, e2019GL084486, https://doi.org/10.1029/2019gl084486, 2020. a, b
Beck, J. P., Reeves, C. E., de Leeuw, F. A., and Penkett, S. A.: The effect of aircraft emissions on tropospheric ozone in the northern hemisphere, Atmos. Environ. Part A.-Gen., 26, 17–29, https://doi.org/10.1016/0960-1686(92)90257-l, 1992. a, b
Bowman, F. M. and Seinfeld, J. H.: Ozone productivity of atmospheric organics, J. Geophys. Res., 99, 5309–5324, 1994. a
Butler, T., Lawrence, M., Taraborrelli, D., and Lelieveld, J.: Multi-day ozone production potential of volatile organic compounds calculated with a tagging approach, Atmos. Environ., 45, 4082–4090, https://doi.org/10.1016/j.atmosenv.2011.03.040, 2011. a, b
Butler, T., Lupascu, A., Coates, J., and Zhu, S.: TOAST 1.0: Tropospheric Ozone Attribution of Sources with Tagging for CESM 1.2.2, Geosci. Model Dev., 11, 2825–2840, https://doi.org/10.5194/gmd-11-2825-2018, 2018. a, b, c, d, e, f, g, h, i, j, k, l, m, n, o, p, q, r, s, t, u, v, w, x, y, z, aa, ab, ac
Carter, W.: Development of ozone reactivity scales for volatile organic compounds, J. Air Waste Manage, 44, 881–899, 1994. a
Chameides, W. L., Fehsenfeld, F., Rodgers, M. O., Cardelino, C., Martinez, J., Parrish, D., Lonneman, W., Lawson, D. R., Rasmussen, R. A., Zimmerman, P., Greenberg, J., Mlddleton, P., and Wang, T.: Ozone precursor relationships in the ambient atmosphere, J. Geophys. Res., 97, 6037, https://doi.org/10.1029/91jd03014, 1992. a, b, c
Clappier, A., Belis, C. A., Pernigotti, D., and Thunis, P.: Source apportionment and sensitivity analysis: two methodologies with two different purposes, Geosci. Model Dev., 10, 4245–4256, https://doi.org/10.5194/gmd-10-4245-2017, 2017. a
Coates, J. and Butler, T. M.: A comparison of chemical mechanisms using tagged ozone production potential (TOPP) analysis, Atmos. Chem. Phys., 15, 8795–8808, https://doi.org/10.5194/acp-15-8795-2015, 2015. a, b, c
Dahlmann, K., Grewe, V., Ponater, M., and Matthes, S.: Quantifying the contributions of individual NOx sources to the trend in ozone radiative forcing, Atmos. Environ., 45, 2860–2868, https://doi.org/10.1016/j.atmosenv.2011.02.071, 2011. a, b
Derwent, R. G.: Representing Organic Compound Oxidation in Chemical Mechanisms for Policy-Relevant Air Quality Models under Background Troposphere Conditions, Atmosphere, 11, 171, https://doi.org/10.3390/atmos11020171, 2020. a
Derwent, R. G., Utembe, S. R., Jenkin, M. E., and Shallcross, D. E.: Tropospheric ozone production regions and the intercontinental origins of surface ozone over Europe, Atmos. Environ., 112, 216–224, https://doi.org/10.1016/j.atmosenv.2015.04.049, 2015. a
Derwent, R. G., Manning, A. J., Simmonds, P. G., Spain, T. G., and O'Doherty, S.: Long-term trends in ozone in baseline and European regionally-polluted air at Mace Head, Ireland over a 30-year period, Atmos. Environ., 179, 279–287, https://doi.org/10.1016/j.atmosenv.2018.02.024, 2018. a, b
Dunker, A., Yarwood, G., Ortmann, J., and Wilson, G.: Comparison of source apportionment and source sensitivity of ozone in a three-dimensional air quality model, Environ. Sci. Technol., 36, 2953–2964, https://doi.org/10.1021/es011418f, 2002. a
Emmerson, K. M. and Evans, M. J.: Comparison of tropospheric gas-phase chemistry schemes for use within global models, Atmos. Chem. Phys., 9, 1831–1845, https://doi.org/10.5194/acp-9-1831-2009, 2009. a, b, c, d
Emmons, L. K., Hess, P. G., Lamarque, J.-F., and Pfister, G. G.: Tagged ozone mechanism for MOZART-4, CAM-chem and other chemical transport models, Geosci. Model Dev., 5, 1531–1542, https://doi.org/10.5194/gmd-5-1531-2012, 2012. a, b
Emmons, L. K., Arnold, S. R., Monks, S. A., Huijnen, V., Tilmes, S., Law, K. S., Thomas, J. L., Raut, J.-C., Bouarar, I., Turquety, S., Long, Y., Duncan, B., Steenrod, S., Strode, S., Flemming, J., Mao, J., Langner, J., Thompson, A. M., Tarasick, D., Apel, E. C., Blake, D. R., Cohen, R. C., Dibb, J., Diskin, G. S., Fried, A., Hall, S. R., Huey, L. G., Weinheimer, A. J., Wisthaler, A., Mikoviny, T., Nowak, J., Peischl, J., Roberts, J. M., Ryerson, T., Warneke, C., and Helmig, D.: The POLARCAT Model Intercomparison Project (POLMIP): overview and evaluation with observations, Atmos. Chem. Phys., 15, 6721–6744, https://doi.org/10.5194/acp-15-6721-2015, 2015. a, b, c
Fiore, A. M., Fischer, E. V., Milly, G. P., Pandey Deolal, S., Wild, O., Jaffe, D. A., Staehelin, J., Clifton, O. E., Bergmann, D., Collins, W., Dentener, F., Doherty, R. M., Duncan, B. N., Fischer, B., Gilge, S., Hess, P. G., Horowitz, L. W., Lupu, A., MacKenzie, I. A., Park, R., Ries, L., Sanderson, M. G., Schultz, M. G., Shindell, D. T., Steinbacher, M., Stevenson, D. S., Szopa, S., Zellweger, C., and Zeng, G.: Peroxy acetyl nitrate (PAN) measurements at northern midlatitude mountain sites in April: a constraint on continental source–receptor relationships, Atmos. Chem. Phys., 18, 15345–15361, https://doi.org/10.5194/acp-18-15345-2018, 2018. a, b, c, d
Fiore, A. M., West, J. J., Horowitz, L. W., Naik, V., and Schwarzkopf, M. D.: Characterizing the tropospheric ozone response to methane emission controls and the benefits to climate and air quality, J. Geophys. Res., 113, D08307, https://doi.org/10.1029/2007JD009162, 2008. a, b, c, d
Fiore, A. M., Dentener, F. J., Wild, O., Cuvelier, C., Schultz, M. G., Hess, P., Textor, C., Schulz, M., Doherty, R. M., Horowitz, L. W., MacKenzie, I. A., Sanderson, M. G., Shindell, D. T., Stevenson, D. S., Szopa, S., Van Dingenen, R., Zeng, G., Atherton, C., Bergmann, D., Bey, I., Carmichael, G., Collins, W. J., Duncan, B. N., Faluvegi, G., Folberth, G., Gauss, M., Gong, S., Hauglustaine, D., Holloway, T., Isaksen, I. S. A., Jacob, D. J., Jonson, J. E., Kaminski, J. W., Keating, T. J., Lupu, A., Marmer, E., Montanaro, V., Park, R. J., Pitari, G., Pringle, K. J., Pyle, J. A., Schroeder, S., Vivanco, M. G., Wind, P., Wojcik, G., Wu, S., and Zuber, A.: Multimodel estimates of intercontinental source-receptor relationships for ozone pollution, J. Geophys. Res., 114, D04301, https://doi.org/10.1029/2008JD010816, 2009. a, b
Fischer, E. V., Jacob, D. J., Yantosca, R. M., Sulprizio, M. P., Millet, D. B., Mao, J., Paulot, F., Singh, H. B., Roiger, A., Ries, L., Talbot, R. W., Dzepina, K., and Pandey Deolal, S.: Atmospheric peroxyacetyl nitrate (PAN): a global budget and source attribution, Atmos. Chem. Phys., 14, 2679–2698, https://doi.org/10.5194/acp-14-2679-2014, 2014. a, b, c, d
Fleming, Z. L., Doherty, R. M., Schneidemesser, E. V., Malley, C. S., Cooper, O. R., Pinto, J. P., Colette, A., Xu, X., Simpson, D., Schultz, M. G., Lefohn, A. S., Hamad, S., Moolla, R., Solberg, S., and Feng, Z.: Tropospheric Ozone Assessment Report: Present-day ozone distribution and trends relevant to human health, Elem. Sci. Anth., 6, p. 12, https://doi.org/10.1525/elementa.273, 2018. a
Galloway, J. N., Townsend, A. R., Erisman, J. W., Bekunda, M., Cai, Z., Freney, J. R., Martinelli, L. A., Seitzinger, S. P., and Sutton, M. A.: Transformation of the Nitrogen Cycle: Recent Trends, Questions, and Potential Solutions, Science, 320, 889–892, https://doi.org/10.1126/science.1136674, 2008. a
Galmarini, S., Koffi, B., Solazzo, E., Keating, T., Hogrefe, C., Schulz, M., Benedictow, A., Griesfeller, J. J., Janssens-Maenhout, G., Carmichael, G., Fu, J., and Dentener, F.: Technical note: Coordination and harmonization of the multi-scale, multi-model activities HTAP2, AQMEII3, and MICS-Asia3: simulations, emission inventories, boundary conditions, and model output formats, Atmos. Chem. Phys., 17, 1543–1555, https://doi.org/10.5194/acp-17-1543-2017, 2017. a, b, c, d, e
Goldstein, A. H. and Galbally, I. E.: Known and unexplored organic constituents in the earth's atmosphere, Environ. Sci. Technol., 41, 1514–1521, 2007. a
Grewe, V., Tsati, E., and Hoor, P.: On the attribution of contributions of atmospheric trace gases to emissions in atmospheric model applications, Geosci. Model Dev., 3, 487–499, https://doi.org/10.5194/gmd-3-487-2010, 2010. a, b, c
Grewe, V., Tsati, E., Mertens, M., Frömming, C., and Jöckel, P.: Contribution of emissions to concentrations: the TAGGING 1.0 submodel based on the Modular Earth Submodel System (MESSy 2.52), Geosci. Model Dev., 10, 2615–2633, https://doi.org/10.5194/gmd-10-2615-2017, 2017. a, b, c, d, e
Haagen-Smit, A. J.: Chemistry and physiology of Los-Angeles smog, Ind. Eng. Chem., 44, 1342–1346, 1952. a
Han, H., Liu, J., Yuan, H., Wang, T., Zhuang, B., and Zhang, X.: Foreign influences on tropospheric ozone over East Asia through global atmospheric transport, Atmos. Chem. Phys., 19, 12495–12514, https://doi.org/10.5194/acp-19-12495-2019, 2019. a, b
Heald, C. L. and Kroll, J. H.: The fuel of atmospheric chemistry: Toward a complete description of reactive organic carbon, Science Advances, 6, eaay8967, https://doi.org/10.1126/sciadv.aay8967, 2020. a
Hoor, P., Borken-Kleefeld, J., Caro, D., Dessens, O., Endresen, O., Gauss, M., Grewe, V., Hauglustaine, D., Isaksen, I. S. A., Jöckel, P., Lelieveld, J., Myhre, G., Meijer, E., Olivie, D., Prather, M., Schnadt Poberaj, C., Shine, K. P., Staehelin, J., Tang, Q., van Aardenne, J., van Velthoven, P., and Sausen, R.: The impact of traffic emissions on atmospheric ozone and OH: results from QUANTIFY, Atmos. Chem. Phys., 9, 3113–3136, https://doi.org/10.5194/acp-9-3113-2009, 2009. a, b
HTAP: Hemispheric Transport of Air Pollution 2010, Part A: Ozone and Particulate Matter, United Nations Publication ECE/EB.AIR/100, edited by: Dentener, F., Keating, T., Akimoto, H., 2010. a, b, c, d, e, f
Huang, M., Carmichael, G. R., Pierce, R. B., Jo, D. S., Park, R. J., Flemming, J., Emmons, L. K., Bowman, K. W., Henze, D. K., Davila, Y., Sudo, K., Jonson, J. E., Tronstad Lund, M., Janssens-Maenhout, G., Dentener, F. J., Keating, T. J., Oetjen, H., and Payne, V. H.: Impact of intercontinental pollution transport on North American ozone air pollution: an HTAP phase 2 multi-model study, Atmos. Chem. Phys., 17, 5721–5750, https://doi.org/10.5194/acp-17-5721-2017, 2017. a, b
Jacob, D., Heikes, B., Fan, S., Logan, J., Mauzerall, D., Bradshaw, J., Singh, H., Gregory, G., Talbot, R., Blake, D., and Sachse, G.: Origin of ozone and NOx in the tropical troposphere: A photochemical analysis of aircraft observations over the South Atlantic basin, J. Geophys. Res.-Oceans, 101, 24235–24250, 1996. a, b
Janssens-Maenhout, G., Crippa, M., Guizzardi, D., Dentener, F., Muntean, M., Pouliot, G., Keating, T., Zhang, Q., Kurokawa, J., Wankmüller, R., Denier van der Gon, H., Kuenen, J. J. P., Klimont, Z., Frost, G., Darras, S., Koffi, B., and Li, M.: HTAP_v2.2: a mosaic of regional and global emission grid maps for 2008 and 2010 to study hemispheric transport of air pollution, Atmos. Chem. Phys., 15, 11411–11432, https://doi.org/10.5194/acp-15-11411-2015, 2015. a
Jonson, J. E., Schulz, M., Emmons, L., Flemming, J., Henze, D., Sudo, K., Tronstad Lund, M., Lin, M., Benedictow, A., Koffi, B., Dentener, F., Keating, T., Kivi, R., and Davila, Y.: The effects of intercontinental emission sources on European air pollution levels, Atmos. Chem. Phys., 18, 13655–13672, https://doi.org/10.5194/acp-18-13655-2018, 2018. a, b, c, d, e
Jonson, J. E., Gauss, M., Schulz, M., Jalkanen, J.-P., and Fagerli, H.: Effects of global ship emissions on European air pollution levels, Atmos. Chem. Phys. Discuss., https://doi.org/10.5194/acp-2020-293, in review, 2020. a
Kasibhatla, P., Levy, H., Moxim, W. J., Pandis, S. N., Corbett, J. J., Peterson, M. C., Honrath, R. E., Frost, G. J., Knapp, K., Parrish, D. D., and Ryerson, T. B.: Do emissions from ships have a significant impact on concentrations of nitrogen oxides in the marine boundary layer?, Geophys. Res. Lett., 27, 2229–2232, https://doi.org/10.1029/2000gl011387, 2000. a
Kleinman, L. I., Daum, P. H., Imre, D., Lee, Y.-N., Nunnermacker, L. J., Springston, S. R., Weinstein-Lloyd, J., and Rudolph, J.: Ozone production rate and hydrocarbon reactivity in 5 urban areas: A cause of high ozone concentration in Houston, Geophys. Res. Lett., 29, 105-1–105-4, https://doi.org/10.1029/2001gl014569, 2002. a
Knote, C., Tuccella, P., Curci, G., Emmons, L., Orlando, J. J., Madronich, S., Baró, R., Jiménez-Guerrero, P., Luecken, D., Hogrefe, C., Forkel, R., Werhahn, J., Hirtl, M., Pérez, J. L., José, R. S., Giordano, L., Brunner, D., Yahya, K., and Zhang, Y.: Influence of the choice of gas-phase mechanism on predictions of key gaseous pollutants during the AQMEII phase-2 intercomparison, Atmos. Environ., 115, 553–568, https://doi.org/10.1016/j.atmosenv.2014.11.066, 2015. a, b
Lamarque, J.-F., Emmons, L. K., Hess, P. G., Kinnison, D. E., Tilmes, S., Vitt, F., Heald, C. L., Holland, E. A., Lauritzen, P. H., Neu, J., Orlando, J. J., Rasch, P. J., and Tyndall, G. K.: CAM-chem: description and evaluation of interactive atmospheric chemistry in the Community Earth System Model, Geosci. Model Dev., 5, 369–411, https://doi.org/10.5194/gmd-5-369-2012, 2012. a
Li, K., Jacob, D. J., Liao, H., Shen, L., Zhang, Q., and Bates, K. H.: Anthropogenic drivers of 2013–2017 trends in summer surface ozone in China, P. Natl. Acad. Sci. USA, 116, 422–427, https://doi.org/10.1073/pnas.1812168116, 2018. a
Li, M., Zhang, Q., Zheng, B., Tong, D., Lei, Y., Liu, F., Hong, C., Kang, S., Yan, L., Zhang, Y., Bo, Y., Su, H., Cheng, Y., and He, K.: Persistent growth of anthropogenic non-methane volatile organic compound (NMVOC) emissions in China during 1990–2017: drivers, speciation and ozone formation potential, Atmos. Chem. Phys., 19, 8897–8913, https://doi.org/10.5194/acp-19-8897-2019, 2019. a
Lin, M., Fiore, A. M., Horowitz, L. W., Cooper, O. R., Naik, V., Holloway, J., Johnson, B. J., Middlebrook, A. M., Oltmans, S. J., Pollack, I. B., Ryerson, T. B., Warner, J. X., Wiedinmyer, C., Wilson, J., and Wyman, B.: Transport of Asian ozone pollution into surface air over the western United States in spring, J. Geophys. Res.-Atmos., 117, D00V07, https://doi.org/10.1029/2011jd016961, 2012. a, b
Liu, F., Beirle, S., Zhang, Q., van der A, R. J., Zheng, B., Tong, D., and He, K.: NOx emission trends over Chinese cities estimated from OMI observations during 2005 to 2015, Atmos. Chem. Phys., 17, 9261–9275, https://doi.org/10.5194/acp-17-9261-2017, 2017. a
Mertens, M., Grewe, V., Rieger, V. S., and Jöckel, P.: Revisiting the contribution of land transport and shipping emissions to tropospheric ozone, Atmos. Chem. Phys., 18, 5567–5588, https://doi.org/10.5194/acp-18-5567-2018, 2018. a, b, c, d, e, f
Mertens, M., Kerkweg, A., Grewe, V., Jöckel, P., and Sausen, R.: Attributing ozone and its precursors to land transport emissions in Europe and Germany, Atmos. Chem. Phys., 20, 7843–7873, https://doi.org/10.5194/acp-20-7843-2020, 2020. a
Mills, G., Pleijel, H., Malley, C. S., Sinha, B., Cooper, O. R., Schultz, M. G., Neufeld, H. S., Simpson, D., Sharps, K., Feng, Z., Gerosa, G., Harmens, H., Kobayashi, K., Saxena, P., Paoletti, E., Sinha, V., and Xu, X.: Tropospheric Ozone Assessment Report: Present-day tropospheric ozone distribution and trends relevant to vegetation, Elem. Sci. Anth., 6, p. 47, https://doi.org/10.1525/elementa.302, 2018. a
Monks, P. S., Archibald, A. T., Colette, A., Cooper, O., Coyle, M., Derwent, R., Fowler, D., Granier, C., Law, K. S., Mills, G. E., Stevenson, D. S., Tarasova, O., Thouret, V., von Schneidemesser, E., Sommariva, R., Wild, O., and Williams, M. L.: Tropospheric ozone and its precursors from the urban to the global scale from air quality to short-lived climate forcer, Atmos. Chem. Phys., 15, 8889–8973, https://doi.org/10.5194/acp-15-8889-2015, 2015. a
Myhre, G., Shine, K., Rädel, G., Gauss, M., Isaksen, I., Tang, Q., Prather, M., Williams, J., van Velthoven, P., Dessens, O., Koffi, B., Szopa, S., Hoor, P., Grewe, V., Borken-Kleefeld, J., Berntsen, T., and Fuglestvedt, J.: Radiative forcing due to changes in ozone and methane caused by the transport sector, Atmos. Environ., 45, 387–394, https://doi.org/10.1016/j.atmosenv.2010.10.001, 2011. a
Parrish, D. D., Petropavlovskikh, I., and Oltmans, S. J.: Reversal of Long-Term Trend in Baseline Ozone Concentrations at the North American West Coast, Geophys. Res. Lett., 44, 10675–10681, https://doi.org/10.1002/2017gl074960, 2017. a, b
Prather, M. J.: Time scales in atmospheric chemistry: Theory, GWPs for CH4 and CO, and runaway growth, Geophys. Res. Lett., 23, 2597–2600, 1996. a
Rao, S., Klimont, Z., Smith, S. J., Dingenen, R. V., Dentener, F., Bouwman, L., Riahi, K., Amann, M., Bodirsky, B. L., van Vuuren, D. P., Reis, L. A., Calvin, K., Drouet, L., Fricko, O., Fujimori, S., Gernaat, D., Havlik, P., Harmsen, M., Hasegawa, T., Heyes, C., Hilaire, J., Luderer, G., Masui, T., Stehfest, E., Strefler, J., van der Sluis, S., and Tavoni, M.: Future air pollution in the Shared Socio-economic Pathways, Global Environ. Chang., 42, 346–358, https://doi.org/10.1016/j.gloenvcha.2016.05.012, 2017. a
Reidmiller, D. R., Fiore, A. M., Jaffe, D. A., Bergmann, D., Cuvelier, C., Dentener, F. J., Duncan, B. N., Folberth, G., Gauss, M., Gong, S., Hess, P., Jonson, J. E., Keating, T., Lupu, A., Marmer, E., Park, R., Schultz, M. G., Shindell, D. T., Szopa, S., Vivanco, M. G., Wild, O., and Zuber, A.: The influence of foreign vs. North American emissions on surface ozone in the US, Atmos. Chem. Phys., 9, 5027–5042, https://doi.org/10.5194/acp-9-5027-2009, 2009. a, b
Rienecker, M. M., Suarez, M. J., Gelaro, R., Todling, R., Bacmeister, J., Liu, E., Bosilovich, M. G., Schubert, S. D., Takacs, L., Kim, G.-K., Bloom, S., Chen, J., Collins, D., Conaty, A., da Silva, A., Gu, W., Joiner, J., Koster, R. D., Lucchesi, R., Molod, A., Owens, T., Pawson, S., Pegion, P., Redder, C. R., Reichle, R., Robertson, F. R., Ruddick, A. G., Sienkiewicz, M., and Woollen, J.: MERRA: NASA's Modern-Era Retrospective Analysis for Research and Applications, J. Climate, 24, 3624–3648, https://doi.org/10.1175/jcli-d-11-00015.1, 2011. a
Saunders, S. M., Jenkin, M. E., Derwent, R. G., and Pilling, M. J.: Protocol for the development of the Master Chemical Mechanism, MCM v3 (Part A): tropospheric degradation of non-aromatic volatile organic compounds, Atmos. Chem. Phys., 3, 161–180, https://doi.org/10.5194/acp-3-161-2003, 2003. a
Saunois, M., Bousquet, P., Poulter, B., Peregon, A., Ciais, P., Canadell, J. G., Dlugokencky, E. J., Etiope, G., Bastviken, D., Houweling, S., Janssens-Maenhout, G., Tubiello, F. N., Castaldi, S., Jackson, R. B., Alexe, M., Arora, V. K., Beerling, D. J., Bergamaschi, P., Blake, D. R., Brailsford, G., Brovkin, V., Bruhwiler, L., Crevoisier, C., Crill, P., Covey, K., Curry, C., Frankenberg, C., Gedney, N., Höglund-Isaksson, L., Ishizawa, M., Ito, A., Joos, F., Kim, H.-S., Kleinen, T., Krummel, P., Lamarque, J.-F., Langenfelds, R., Locatelli, R., Machida, T., Maksyutov, S., McDonald, K. C., Marshall, J., Melton, J. R., Morino, I., Naik, V., O'Doherty, S., Parmentier, F.-J. W., Patra, P. K., Peng, C., Peng, S., Peters, G. P., Pison, I., Prigent, C., Prinn, R., Ramonet, M., Riley, W. J., Saito, M., Santini, M., Schroeder, R., Simpson, I. J., Spahni, R., Steele, P., Takizawa, A., Thornton, B. F., Tian, H., Tohjima, Y., Viovy, N., Voulgarakis, A., van Weele, M., van der Werf, G. R., Weiss, R., Wiedinmyer, C., Wilton, D. J., Wiltshire, A., Worthy, D., Wunch, D., Xu, X., Yoshida, Y., Zhang, B., Zhang, Z., and Zhu, Q.: The global methane budget 2000–2012, Earth Syst. Sci. Data, 8, 697–751, https://doi.org/10.5194/essd-8-697-2016, 2016. a
Schultz, M. G., Schröder, S., Lyapina, O., Cooper, O., Galbally, I., Petropavlovskikh, I., Schneidemesser, E. V., Tanimoto, H., Elshorbany, Y., Naja, M., Seguel, R., Dauert, U., Eckhardt, P., Feigenspahn, S., Fiebig, M., Hjellbrekke, A.-G., Hong, Y.-D., Kjeld, P. C., Koide, H., Lear, G., Tarasick, D., Ueno, M., Wallasch, M., Baumgardner, D., Chuang, M.-T., Gillett, R., Lee, M., Molloy, S., Moolla, R., Wang, T., Sharps, K., Adame, J. A., Ancellet, G., Apadula, F., Artaxo, P., Barlasina, M., Bogucka, M., Bonasoni, P., Chang, L., Colomb, A., Cuevas, E., Cupeiro, M., Degorska, A., Ding, A., Fröhlich, M., Frolova, M., Gadhavi, H., Gheusi, F., Gilge, S., Gonzalez, M. Y., Gros, V., Hamad, S. H., Helmig, D., Henriques, D., Hermansen, O., Holla, R., Huber, J., Im, U., Jaffe, D. A., Komala, N., Kubistin, D., Lam, K.-S., Laurila, T., Lee, H., Levy, I., Mazzoleni, C., Mazzoleni, L., McClure-Begley, A., Mohamad, M., Murovic, M., Navarro-Comas, M., Nicodim, F., Parrish, D., Read, K. A., Reid, N., Ries, L., Saxena, P., Schwab, J. J., Scorgie, Y., Senik, I., Simmonds, P., Sinha, V., Skorokhod, A., Spain, G., Spangl, W., Spoor, R., Springston, S. R., Steer, K., Steinbacher, M., Suharguniyawan, E., Torre, P., Trickl, T., Weili, L., Weller, R., Xu, X., Xue, L., and Zhiqiang, M.: Tropospheric Ozone Assessment Report: Database and Metrics Data of Global Surface Ozone Observations, Elem. Sci. Anth., 5, p. 58, https://doi.org/10.1525/elementa.244, 2017. a, b
Stevenson, D. S., Young, P. J., Naik, V., Lamarque, J.-F., Shindell, D. T., Voulgarakis, A., Skeie, R. B., Dalsoren, S. B., Myhre, G., Berntsen, T. K., Folberth, G. A., Rumbold, S. T., Collins, W. J., MacKenzie, I. A., Doherty, R. M., Zeng, G., van Noije, T. P. C., Strunk, A., Bergmann, D., Cameron-Smith, P., Plummer, D. A., Strode, S. A., Horowitz, L., Lee, Y. H., Szopa, S., Sudo, K., Nagashima, T., Josse, B., Cionni, I., Righi, M., Eyring, V., Conley, A., Bowman, K. W., Wild, O., and Archibald, A.: Tropospheric ozone changes, radiative forcing and attribution to emissions in the Atmospheric Chemistry and Climate Model Intercomparison Project (ACCMIP), Atmos. Chem. Phys., 13, 3063–3085, https://doi.org/10.5194/acp-13-3063-2013, 2013. a, b
Thunis, P., Clappier, A., Tarrason, L., Cuvelier, C., Monteiro, A., Pisoni, E., Wesseling, J., Belis, C., Pirovano, G., Janssen, S., Guerreiro, C., and Peduzzi, E.: Source apportionment to support air quality planning: Strengths and weaknesses of existing approaches, Environ. Int., 130, 104825, https://doi.org/10.1016/j.envint.2019.05.019, 2019. a
Tilmes, S., Lamarque, J.-F., Emmons, L. K., Kinnison, D. E., Ma, P.-L., Liu, X., Ghan, S., Bardeen, C., Arnold, S., Deeter, M., Vitt, F., Ryerson, T., Elkins, J. W., Moore, F., Spackman, J. R., and Val Martin, M.: Description and evaluation of tropospheric chemistry and aerosols in the Community Earth System Model (CESM1.2), Geosci. Model Dev., 8, 1395–1426, https://doi.org/10.5194/gmd-8-1395-2015, 2015. a, b, c, d, e, f, g
Utembe, S., Cooke, M., Archibald, A., Jenkin, M., Derwent, R., and Shallcross, D.: Using a reduced Common Representative Intermediates (CRIv2-R5) mechanism to simulate tropospheric ozone in a 3-D Lagrangian chemistry transport model, Atmos. Environ., 44, 1609–1622, 2010. a, b, c
van der Werf, G. R., Randerson, J. T., Giglio, L., Collatz, G. J., Mu, M., Kasibhatla, P. S., Morton, D. C., DeFries, R. S., Jin, Y., and van Leeuwen, T. T.: Global fire emissions and the contribution of deforestation, savanna, forest, agricultural, and peat fires (1997–2009), Atmos. Chem. Phys., 10, 11707–11735, https://doi.org/10.5194/acp-10-11707-2010, 2010. a
Vinken, G. C. M., Boersma, K. F., Jacob, D. J., and Meijer, E. W.: Accounting for non-linear chemistry of ship plumes in the GEOS-Chem global chemistry transport model, Atmos. Chem. Phys., 11, 11707–11722, https://doi.org/10.5194/acp-11-11707-2011, 2011. a, b
von Glasow, R., Lawrence, M. G., Sander, R., and Crutzen, P. J.: Modeling the chemical effects of ship exhaust in the cloud-free marine boundary layer, Atmos. Chem. Phys., 3, 233–250, https://doi.org/10.5194/acp-3-233-2003, 2003. a, b, c
Wang, Y., Jacob, D. J., and Logan, J. A.: Global simulation of tropospheric O3-NOx-hydrocarbon chemistry: 3. Origin of tropospheric ozone and effects of nonmethane hydrocarbons, J. Geophys. Res.-Atmos., 103, 10757–10767, https://doi.org/10.1029/98jd00156, 1998. a
Young, P. J., Archibald, A. T., Bowman, K. W., Lamarque, J.-F., Naik, V., Stevenson, D. S., Tilmes, S., Voulgarakis, A., Wild, O., Bergmann, D., Cameron-Smith, P., Cionni, I., Collins, W. J., Dalsøren, S. B., Doherty, R. M., Eyring, V., Faluvegi, G., Horowitz, L. W., Josse, B., Lee, Y. H., MacKenzie, I. A., Nagashima, T., Plummer, D. A., Righi, M., Rumbold, S. T., Skeie, R. B., Shindell, D. T., Strode, S. A., Sudo, K., Szopa, S., and Zeng, G.: Corrigendum to “Pre-industrial to end 21st century projections of tropospheric ozone from the Atmospheric Chemistry and Climate Model Intercomparison Project (ACCMIP)” published in Atmos. Chem. Phys., 13, 2063–2090, 2013, Atmos. Chem. Phys., 13, 5401–5402, https://doi.org/10.5194/acp-13-5401-2013, 2013. a
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. R., 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, Elem. Sci. Anth., 6, p.10, https://doi.org/10.1525/elementa.265, 2018. a, b, c, d
Zhang, Y., Cooper, O. R., Gaudel, A., Thompson, A. M., Nédélec, P., Ogino, S.-Y., and West, J. J.: Tropospheric ozone change from 1980 to 2010 dominated by equatorward redistribution of emissions, Nat. Geosci., 9, 875–879, https://doi.org/10.1038/ngeo2827, 2016. a, b, c