Effects of heterogeneous reactions on tropospheric chemistry: a global simulation with the chemistry climate model CHASER V4.0

. This study uses a chemistry-climate model CHASER (MIROC) to explore the roles of heterogeneous reactions (HRs) in global tropospheric chemistry. Three distinct HRs of N 2 O 5 , HO 2 , and RO 2 are considered for surfaces of aerosols and cloud particles. The model simulation is verified with EANET and EMEP stationary observations, R/V MIRAI ship-based 10 data, ATom1 aircraft measurements, satellite observations by OMI, ISCCP, and CALIPSO-GOCCP, and reanalysis data JRA55. The heterogeneous chemistry facilitates improvement of model performance with respect to observations for NO 2 , OH, CO, and O 3 , especially in the lower troposphere. The calculated effects of heterogeneous reactions cause marked changes in global abundances of O 3 (-2.96%), NO x (-2.19%), CO (+3.28%), and global mean CH 4 lifetime (+5.91%). These global effects were contributed mostly by N 2 O 5 uptake onto aerosols in the middle troposphere. At the surface, HO 2 uptake gives the 15 largest contributions, with a particularly significant effect in the North Pacific region (-24% O 3 , +68% NO x , +8% CO, and - 70% OH), mainly attributable to its uptake onto clouds. The RO 2 reaction has a small contribution, but its global-mean negative effects on O 3 and CO are not negligible. In general, the uptakes onto ice crystals and cloud droplets that occur mainly by HO 2 and RO 2 radicals cause smaller global effects than the aerosol-uptake effects by N 2 O 5 radicals (+1.34% CH 4 lifetime, +1.71% NO x , -0.56% O 3 , +0.63% CO abundances). Nonlinear responses of tropospheric O 3 , NO x , and OH to the N 2 O 5 and HO 2 uptakes 20 are found in the same modelling framework of this study ( R > 0.93). Although all HRs showed negative tendencies for OH and O 3 levels, the effects of HR(HO 2 ) on the tropospheric abundance of O 3 showed a small increment with an increasing loss rate. However, this positive tendency turns to reduction at higher rates (>5 times). Our results demonstrate that the HRs affect not only polluted areas but also remote areas such as the mid-latitude sea boundary layer and upper troposphere. Furthermore, HR(HO 2 ) can bring challenges to pollution reduction efforts because it causes opposite effects between NO x (increase) and 25 surface O 3 (decrease).


Introduction
Heterogeneous reactions (HRs) on the surfaces of atmospheric aerosols and cloud droplets are regarded as playing crucial roles in atmospheric chemistry.They affect ozone (O 3 ) concentrations in various pathways via the cycle of odd hydrogen (HO x ) and nitrogen oxides (NO x ) (Jacob, 2000).Tropospheric ozone, an important greenhouse gas, causes damage to human health, crops, and ecosystem productivity (Monks et al., 2015).Although tropospheric O 3 was recognized as a critical oxidant species, its global distribution has not been adequately captured to date because of the limited number of observations.Whereas many sites in the heavily polluted regions of eastern Asia show ozone increases since 2000 (Liu and Wang, 2020), many sites in other regions show decreases (Gaudel et al., 2018).Moreover, O 3 responds to changes of multiple pollutants such as NO x and volatile organic compounds (VOCs) in different ways, which challenge the local pollutant control policy.For instance, since the Chinese government released the Air Pollution Prevention and Control Action Plan in 2010 (Zheng et al., 2018), the targets of SO 2 , NO x , and particulate matter (PM) decreased drastically, but urban ozone pollution has been worsening (Liu and Wang, 2020).Indeed, the O 3 responses are controlled by several mechanisms, including heterogeneous effects of HO 2 and N 2 O 5 onto aerosols (Kanaya et al., 2009;Li et al., 2019;Liu and Wang, 2020;Taketani et al., 2012).
Stationary observations and laboratory experiments are important for enhancing the understanding of the tropospheric chemistry of O 3 and other essential components (NO x , HO x ).However, direct observation of vertical O 3 distribution, including upper tropospheric O 3 , was not available before 1970.It has been deployed only at limited sites across the globe.Global atmospheric modeling is a useful method to reanalyze or forecast the past and future changes in O 3 and their effects on human health and the climate.To serve this task, atmospheric models use both laboratory and observational data to help achieve accurate simulations of O 3 and its precursors (HO x , NO x , hydrocarbons).To date, many modeling studies have suggested that heterogeneous chemistry be included in a standard model for tropospheric chemistry (Jacob, 2000;Macintyre andEvans, 2010, 2011;de Reus et al., 2005).
One fundamentally important HR in the troposphere is the uptake of N 2 O 5 onto aqueous aerosols, known as a removal pathway for NO x at night (Platt et al., 1984).Actually, NO x plays crucially important roles in the troposphere because it controls the cycle of HO x and the production rate of tropospheric O 3 (Logan et al., 1981;Riemer et al., 2003).The morning photochemistry can be affected by NO 3 and N 2 O 5 , which are important nocturnal oxidants.Since the early 1980s, the role of urban NO x chemistry in Los Angeles pollution (National Research Council, 1991) has been acknowledged, but the proclamation of nighttime radicals remained sparse.It was only recognized in the past decade that N 2 O 5 radical chemistry could have a much more perceptible effect stemming from reasons including a refined understanding of heterogeneous processes occurring at night (Brown and Stutz, 2012).The HR of N 2 O 5 was revealed under different meteorological conditions in the US, Europe, and China (photosmog, high relative humidity (RH), or seasonal variation) for particles of various types: ice, aqueous aerosols with organic-coating, urban aerosols, dust, and soot (Apodaca et al., 2008;Lowe et al., 2015;Qu et al., 2019;Riemer et al., 2003Riemer et al., , 2009;;Wang et al., 2018Wang et al., , 2017;;Xia et al., 2019).The uptake of N 2 O 5 can markedly enhance nitrate concentration in nocturnal chemistry or PM 2.5 explosive growth events in summer, decrease NO x , and either increase or decrease O 3 concentrations in different NO x conditions (Dentener and Crutzen, 1993;Qu et al., 2019;Riemer et al., 2003;Wang et al., 2017).Even during daytime, N 2 O 5 in the marine boundary layers can enhance the NO x to HNO 3 conversion, and chemical destruction of O 3 (Osthoff et al., 2006).A 10-20 ppbv reduction of O 3 because of N 2 O 5 uptake in the polluted regions of China has also been reported (Li et al., 2018).At mid-to high latitudes, N 2 O 5 uptakes on sulfate aerosols could engender 80 % and 10 % NO x reduction, respectively, in winter and summer, leading to an ap-proximate 10 % reduction of O 3 in both seasons (Li et al., 2018).
Another vital process taking place on particles is the HRs of peroxy radicals (HO 2 and RO 2 ).Peroxy radicals are the primary chain carriers driving O 3 production in the troposphere.Moreover, it can drive the hydrocarbon and NO x concentrations, which are important for nocturnal radical chemistry (Geyer et al., 2003;Richard, 2000;Salisbury et al., 2001).In the past, the HR(HO 2 ) effects have been well considered in the laboratory (Macintyre and Evans, 2011) and field observations (Kanaya et al., 2001(Kanaya et al., , 2002a(Kanaya et al., , b, 2003(Kanaya et al., , 2007;;Taketani et al., 2012), but many technical problems (e.g., detecting HO 2 ) have created difficulties that challenge its reported importance in the troposphere, as asserted from recent studies (Liao and Seinfeld, 2005;Martin et al., 2003;Tie et al., 2001).More recently, global modeling reports have described that the inclusion of HO 2 uptake can affect atmospheric constituents strongly by the increment in tropospheric abundances for carbon monoxide (CO) and other trace gases because of reduced oxidation capacity (Lin et al., 2012;Macintyre and Evans, 2011).The HO x loss on aerosols can reduce O 3 concentrations by up to 33 % in remote areas and up to 10 % in a smog episode (Saathoff et al., 2001;Taketani et al., 2012).The HO x loss on sea salt, sulfate, and organic carbon in various environments can decrease HO 2 levels by 6 %-13 %, 10 %-40 %, and 40 %-70 %, respectively (Martin et al., 2003;Taketani et al., 2008Taketani et al., , 2009;;Tie et al., 2001).For RO 2 with a typical representative of CH 3 CO.O 2 (peroxyacetyl radical, PA), it plays a big role in the longrange transport of pollution (VOC, NO x ) (Richard, 2000;Villalta et al., 1996).It can bring NO x from polluted domains as peroxyacyl nitrates (PAN) to remote regions in the ocean and higher altitudes (Qin et al., 2018;Richard, 2000).The concentrations of HO 2 and RO 2 at nighttime in the marine boundary layer were measured and confirmed (Geyer et al., 2003;Salisbury et al., 2001).Moreover, some evidence suggests uptake of HO 2 and PA on clouds, aqueous aerosols, and other surfaces in high-humidity conditions, although the mechanism is uncertain (Geyer et al., 2003;Jacob, 2000;Kanaya et al., 2002b;Liao and Seinfeld, 2005;Lin et al., 2012;Richard, 2000;Salisbury et al., 2001).The predominance of peroxy uptake to clouds results from the ubiquitous existence and larger SAD maxima of cloud droplets in the atmosphere.Indeed, aqueous-phase chemistry might represent an important sink for O 3 (Lelieveld and Crutzen, 1990).In addition, PA loss on aqueous particles can mediate the loss of PAN (CH 3 CO.O 2 NO) in fog (Villalta et al., 1996).Some modeling studies indicate that HO x loss (including HO 2 loss) on aqueous aerosols reduces OH by 2 % , increases CO by 7 % and increases O 3 by 0.5 % in the annual mean global burden (Huijnen et al., 2014).However, in a coastal environment in the Northern Hemisphere it increases OH by 15 % and reduces HO 2 by 30 % (Sommariva et al., 2006;Thornton et al., 2008).
Although the contributions of each uptake category to tropospheric chemistry differ and must be considered both separately and as a whole, few studies have provided a global overview of heterogeneous chemistry the comprehensively examines the uptakes of N 2 O 5 , HO 2 , and RO 2 on widely various particles.For instance, uptakes of both N 2 O 5 and HO 2 tend to reduce O 3 in particular environments (Li et al., 2018;Saathoff et al., 2001;Taketani et al., 2012), but the HO 2 loss on clouds can increase the tropospheric O 3 burden (Huijnen et al., 2014).The latter trend is not widely suggested yet because the cloud chemistry is still neglected in many O 3 models (Stadtler et al., 2018;Thornton et al., 2008).The predominant effects of HO 2 uptake on aerosols compared to the effect by N 2 O 5 were reported during the summer smog condition (Saathoff et al., 2001) but with lack of confirmation on a global scale.Moreover, the heterogeneous effects of RO 2 have been investigated only insufficiently (Jacob, 2000).In this study, we examine these uncertainties using the global model CHASER to perceive the respective and total effects of the HRs of N 2 O 5 , HO 2 , and RO 2 on the tropospheric chemistry.For the interface of HRs in the atmosphere, we tentatively consider surfaces of cloud particles and those of aerosols and discuss details of its effects in this study.In the following text, the research method, including model description and configuration, is described in Sect. 2. In Sect.3.1, our model is verified with available observations including ground stations and ship, aircraft, and satellite measurements, particularly addressing the roles of the HRs.The global effects of N 2 O 5 , HO 2 , and RO 2 uptake are discussed in Sect.3.2 to elucidate cloud particles and aerosol effects.Section 3.3 will discuss sensitivities of tropospheric chemistry to the magnitudes of HRs.Section 4 presents a summary and concluding remarks.

Global chemistry model
The global chemistry model used for this study is CHASER (MIROC-ESM) (Sudo et al., 2002;Sudo and Akimoto, 2007;Watanabe et al., 2011), which considers detailed photochemistry in the troposphere and stratosphere.The chemistry component of the model, based on CHASER-V4.0,calculates the concentrations of 92 chemical species and 262 chemical reactions (58 photolytic, 183 kinetic, and 21 heterogeneous reactions including reactions on polar stratospheric clouds); more details on CHASER can be found in an earlier report of the literature (Morgenstern et al., 2017).Its tropospheric chemistry considers the fundamental chemical cycle of O x -NO x -HO x -CH 4 -CO, along with oxidation of nonmethane volatile organic compounds (NMVOCs).Its stratospheric chemistry simulates chlorine and bromine-containing compounds, CFCs, HFCs, carbonyl sulfide (OCS), NO 2 , and the formation of polar stratospheric clouds (PSCs) and het-erogeneous reactions on PSC surfaces.In the framework of MIROC-Chem, CHASER is coupled with the MIROC-AGCM atmospheric general circulation model (version 4; Watanabe et al., 2011).The meteorological fields simulated by MIROC-AGCM were nudged toward the 6-hourly NCEP FNL data (https://rda.ucar.edu/datasets/ds083.2/, last access: 30 October 2018).For this study, the spatial resolution of the model was set as T42 (about 2.8 • × 2.8 • grid spacing) in horizontal and L36 (surface to approx.50 km) in vertical.Anthropogenic emissions for O 3 and aerosol precursors like NO x , CO, VOCs, and SO 2 are specified using the HTAP-II inventory (Janssens-Maenhout et al., 2015), with biomass burning emissions derived from the MACC reanalysis system (Inness et al., 2013).
In the model, the aerosol concentrations for black carbon (BC) / organic carbon (OC), sea salt, and soil dust are handled by the SPRINTAR module, which is also based on the CCSR/NIES AGCM (Takemura et al., 2000).The bulk thermodynamics for aerosols are applied, including SO 2− 4 chemistry (SO 2 oxidation with OH, O 3 /H 2 O 2 , which is cloud-pH dependent) SO 2− 4 -NO − 3 -NH + 4 and SO 2− 4 -dust interaction.

Heterogeneous reactions in the chemistry-climate model (CHASER)
The CHASER-V4 model considers HRs in both the troposphere and stratosphere.In this work, we particularly examine HRs in the troposphere.In the current version of CHASER, tropospheric HRs are considered for N 2 O 5 , HO 2 , and RO 2 , using uptake coefficients for the distinct surfaces of aerosols (sulfate, sea salt, dust, and organic carbons) and cloud particles (liquid/ice) as listed in Table 2.Although some other views incorporate the catalysis of transition metal ions (TMIs) Cu(I)/Cu(II) and Fe(II)/Fe(III) for the HO 2 conversion on aqueous aerosols (Li et al., 2018;Mao et al., 2013;Taketani et al., 2012), this mechanism remains uncertain (Jacob, 2000).The TMI mechanism might lead to either H 2 O 2 (Jacob, 2000) or H 2 O products (Mao et al., 2013).However, this may not cause any significant difference, since recycling HO 2 from H 2 O 2 is ineffective (Li et al., 2018).For this study, the uptake of HO 2 is affirmed with H 2 O 2 as the product (Loukhovitskaya et al., 2009;Taketani et al., 2009), as it is generally used in many atmospheric models such that this is not counted as a terminal sink for HO 2 (Jacob, 2000;Lelieveld and Crutzen, 1990;Morita et al., 2004;Thornton et al., 2008).The RO 2 uptakes are assumed with inert products, as suggested by Jacob (2000).The heterogeneous pseudofirst-order loss rate β for the species i is given using the theory of Schwartz (Dentener and Crutzen, 1993;Jacob, 2000;Schwartz, 1986), in which it is simply treated with the mass transfer limitations operating two conductances representing free molecular and continuum regimes for tropospheric clouds and aerosols, in addition to using reactive uptake coefficient (γ ) instead of the mass accommodation coefficient https://doi.org/10.5194/gmd-14-3813-2021Geosci.Model Dev., 14, 3813-3841, 2021 as follows: where ν i stands for the mean molecular speed (cm s −1 ) of species i, D ij is the gaseous mass transfer (diffusion) coefficient (cm 2 s −1 ) of species i for particle type j , and A j expresses the surface area density (cm 2 cm −3 ) for particle type j .In the model, the particle size and effective radius R j for aerosols are calculated as a function of RH (Takemura et al., 2000).The aerosol concentrations are based on SPRINTAR for BC/OC, sea salt, and dust (Takemura et al., 2000).The surface area density (SAD) for aerosols (A j ) is estimated using lognormal distributions of particle size (SF j ) with mode radii variable with the RH (Sudo et al., 2002) as follows where C N represents number density (cm −3 ) and R j signifies the effective radii (cm) of particle type j .To calculate SAD for cloud particles, the liquid water content (LWC) and ice water content (IWC) in the AGCM are converted using the cloud droplet distribution of Battan and Reitan (1957) and the relation between IWC and the surface area density for ice clouds (Lawrence and Crutzen, 1998;McFarquhar and Heymsfield, 1996).
In these equations, A c represents the cross-section area for ice crystals (cm 2 cm −3 ).For liquid clouds, the following holds: The uptake coefficient parameter (γ ) is defined as the net probability that a molecule X undergoing a gas-kinetic collision with a surface is actually taken up onto the surface.
Although several recent model studies that consider dependency of γ on RH and/or T, the majority of the earlier studies use constant γ values that only vary with aerosol particle compositions (Chen et al., 2018;Evans and Jacob, 2005;Macintyre andEvans, 2010, 2011).For one study, γ HO 2 for the uptake onto aqueous aerosols is considered with pH dependence (Thornton et al., 2008).However, another study demonstrated that the uptake is large, irrespective of the solubility in cloud water or pH (Morita et al., 2004).Therefore, we instead choose γ HO 2 as fixed values depending on the type of particle.Indeed, from Eq. (1) it is apparent that uptake coefficients should be unimportant for uptake onto large particles such as cloud droplets.In this study, γ for cloud particles of liquid and ice phases are given based on suggestions from earlier reports (Dentener and Crutzen, 1993;Jacob, 2000).One study (Dentener and Crutzen, 1993) used a constant γ N 2 O 5 of 0.1 for uptake on sea salt, sulfate, and cloud particles.They also revealed that a smaller γ N 2 O 5 of 0.01, which had been reported as laboratory measurements, is insensitive to effects on tropospheric oxidant components.Results of another study (Jacob, 2000) indicated constants γ N 2 O 5 = 0.1 and γ HO 2 = 0.2 for uptakes on both liquid clouds and aerosols, the latter aiming to involve HO 2 scavenging by clouds without accounting for details of aqueous-phase chemistry.For ice crystals, Jacob (2000) suggested γ HO 2 = 0.025 based on a report by Cooper and Abbatt (1996).Jacob (2000)

Experiment setup
In this study, simulations of two types were conducted to isolate the distinct effects of each HR for the surface types considered in the model (Tables 3 and S1).Whereas a control simulation standard (STD) run considers all HRs, cases with no HRs (noHR) cases intentionally ignore one or all of the HRs to calculate effects of individual HRs.The sensitivity runs turned off the separate HRs onto clouds (liquid and ice), and aerosols were also added to exploit the separate aerosol-heterogeneous and cloud-heterogeneous effects, as suggested in many earlier studies (Apodaca et al., 2008;Jacob, 2000;Lelieveld andCrutzen, 1990, 1991;Morita et al., 2004).All simulations were run in the 2009-2017 timeframe, with 2009 being treated as a spin-up year.The HR effects are determined as the differences between noHR cases and an STD simulation as in Eq. ( 5): where STD i stands for the concentration of investigated atmospheric component i in the STD run and noHR(j ) i denotes the concentration of component i in the sensitivity run in which the HRs of/onto j was ignored (j could be N 2 O 5 , HO 2 , RO 2 , clouds, aerosols).An additional sensitivity test was run to examine the sensitivity of the troposphere's responses with the amplified HRs magnitudes (Table S1).These simulations only apply for HR(N 2 O 5 ) and HR(HO 2 ) to verify some uncertainties that have been argued among earlier studies (Chen et al., 2018;Evans and Jacob, 2005;Macintyre andEvans, 2010, 2011).
product 0 0.001 0.004 0.004 0.004 0.004 References and the details of these adjustments are given in the main text.The RO 2 uptakes are assumed with inert products, as suggested by Jacob (2000).ISO 2 is denoted for peroxy radicals from C 5 H 8 + OH, and MACRO 2 stands for peroxy radicals from methacrolein (CH 2 = C(CH 3 )CHO).

Observation data for model evaluation
Model simulations with and without HRs are evaluated distinctively with stationary, ship-based, aircraft-based, and satellite-based measurements.The observational information, locations of the surface sites, and ship or aircraft tracks for the observations used for this study are summarized in Table 4 and Fig. 1.
EANET is well known as the Acid Deposition Monitoring Network in eastern Asia.The monthly data from 45 stations over 13 countries during 2010-2016 were used to verify surface concentrations of aerosols (sulfate, nitrate) and trace gases (HNO 3 , NO x , O 3 ) in eastern Asia.We also used data of the European Monitoring and Evaluation Programme (EMEP), which compiles observations over 245 European stations.
Additionally, we exploited ship-based observational data from R/V Mirai cruise (http://www.godac.jamstec.go.jp/ darwin/e, last access: 30 June 2020) undertaken by the Japan Agency for Marine-Earth Science and Technology (JAMSTEC).This study used data for surface CO and O 3 concentrations in summer 2015-2017 along the Japan-Alaska and Japan-Indonesia-Australia routes (Kanaya et al., 2019).The model data were compiled in hourly time steps and were interpolated corresponding with the Mirai time step and coordinates.For verification of the vertical tropospheric profiles, we used Atmospheric Tomography (ATom1) aircraft measurements (https://espo.nasa.gov/atom/content/ATom, last access: 30 June 2020) for NO 2 , OH, CO, and O 3 .
The simulated tropospheric ozone was also evaluated using the tropospheric column O 3 (TCO) derived from the OMI satellite data (https://daac.gsfc.nasa.gov/,https://doi.org/10.5194/gmd-14-3813-2021 Geosci.Model Dev., 14, 3813-3841, 2021  For this study, we tentatively consider HRs on the cloud particle surface.Given the great uncertainties related to the reaction coefficient (γ ) (Macintyre andEvans, 2010, 2011), the cloud distributions must be examined adequately in the model to the greatest extent possible.The model-calculated cloud distributions were verified using satellite observation data ISCCP D2, CALIPSO-GOCCP, and reanalysis data JRA55.
For the entire troposphere, the calculated cloud fraction was generally underestimated against the satellite observations and reanalysis data (Fig. 2, the first row).In the North Pacific region in JJA (Fig. 2, the second row), when the cloud fraction peaked in the region, the model was able to reproduce the satellite observations (ISCCP and CALIPSO).How-ever, for the lower troposphere over the region, the cloud fraction calculated using CHASER in JJA appears to be overestimated (Fig. 2, the fourth row), suggesting that the resulting HR effects would also be exaggerated to some extent.

Verification with stationary observations
Verifications with EANET and EMEP stationary observations were conducted to assess the model performance on land domains of eastern Asia and Europe, particularly addressing the roles of the heterogeneous reactions considered for this study.
The mass concentrations of particulate matter (PM 2.5 ); sulfate (SO 2− 4 ); nitrate (NO − 3 ); and gaseous HNO 3 , NO x , O 3 , and CO (CO only for EMEP) of 2010-2016 were evaluated (see Figs. S1 to S8 for monthly concentrations and Fig. S9 for correlations).In general, the model can moderately reproduce the PM 2.5 , SO 2− 4 , and NO − 3 aerosol concentrations at these locations (R = 0.3-0.7,Table 5), although PM 2.5 was underestimated, sulfate was overestimated slightly.Nitrate was underestimated for EANET and overestimated for EMEP.It is noteworthy that the model performance for EMEP stations was better than that for EANET.The PM 2.5 concentration was better estimated with the inclusion of N 2 O 5 and HO 2 uptakes (bias reduction in Table 5).Figure 3a-g present the median values of NO x , O 3 , and CO for grouped stations in the Chinese and South Korean regions (stations in China: Jinyunshan; stations in South Korea: Kanghwa, Imsil, Cheju), remote stations with low NO x levels of EANET, and all EMEP stations.Figure 3h-n show changes in NO x , O 3 , and CO for these stations.The model's positive bias for NO − 3 at Kanghwa as a remote area is different from the model underestimates at other EANET stations (e.g., Bangkok, Hanoi, and Hongwen in Fig. S1).These high negative biases for NO − 3 can be associated with undervaluation for NO x and can thereby lessen the effects of N 2 O 5 uptake.
Nitric acid in both regions was overestimated.The correlations, biases, and normalized root-mean-square error (NRMSE) of the model for SO 2− 4 , NO − 3 , and HNO 3 are in the ranges reported in a multi-model study by Bian et al. (2017) (Table 6).In ISCCP and CALIPSO data, the pressure boundary layer of the low troposphere is > 680 hPa.In JRA55, the low troposphere was defined as 850-1100 hPa of pressure.
The NO x concentration for eastern Asia and Europe was underestimated, with significant bias for polluted Asian locations (Bangkok, Metro Manila, Nai Muaeng, Samutprakarn, Si Phum, Ulaanbaatar, not shown).In Fig. 3a and c, simulated NO x levels still underestimated the observed values for Chinese and European regions, at which the observed NO x could reach 16 and 7 ppb, respectively.For the low-NO x EANET region, excluding the abovementioned sites (Fig. 3b), simulated NO x levels turned to overestimate the observed levels in January, February, September, and October.The increasing effects of NO x attributable to heterogeneous reactions, although minor, mitigated these underestimations (Fig. 3a-c).Although NO x was partly reduced via uptake of N 2 O 5 , the NO x level was mostly increased because of HO 2 and RO 2 uptakes (Fig. 3h-j).
In this comparison, the low correlations of the model with EANET and EMEP sites for HNO 3 and NO x are still a problem.The high biases for nitrogen species could be ascribed to the low horizontal resolution in this study (∼ 2.8 • ).Higher resolutions could improve the model reproduction for sur-face NO x as previously investigated by Sekiya et al. (2018).Moreover, the low reproducibility of the model for NO x is probably caused by lacking mechanisms that reduce HNO 3 and enhance NO x in our model.One possible mechanism is the heterogeneous reaction of HNO 3 on soot surfaces (Reaction R8) (Akimoto et al., 2019, and references therein): The additional Reaction (R4) followed by NO 2 uptakes onto soot (Jacob, 2000), NO 2 + particles → 0.5 HONO + 0.5 HNO 3 (Reaction R9), can be expected to increase NO and decrease O 3 via the consequent titration reaction.These changes could reduce the model overestimates for HNO 3 and O 3 and the model underestimates for NO x with EANET and EMEP stations.Further tests for this issue shall be discussed in a future report.
CO for EMEP was partly underestimated by the model, especially during January-March (Fig. 3g).This underestimate was mitigated by increasing effects because of HRs of N 2 O 5 and HO 2 .The uptakes of RO 2 , in contrast, minorly reduced CO levels (Fig. 3n) so that the model bias was wors-ened slightly.For O 3 , whereas the model tends to overestimate this tracer for both regions (Fig. 3d-f), O 3 reduction effects of all HRs (Fig. 3k-m) alleviated the model overestimates from April to December, although advanced reduction is still needed.In January-March, the model tended to underestimate O 3 levels (Fig. 3d-f), which was exaggerated by reduction effects for O 3 .In general, the STD simulation with coupled HRs partly improved the agreement related to the particulate and gaseous species, showing less bias than that of simulations without HRs (Table 5).

Verification with ship-based measurements
The model simulations were also verified with O 3 and CO observations from the Research Vessel (R/V) Mirai for the Pacific Ocean region.This study specifically examines data from the four cruises of R/V Mirai for the Japan-Alaska region (40-75 • N, 140 • E-150 • W) in summer, designated as MR15-03 leg 1 and leg 2 (28 August-21 October 2015, labeled as Track 1 in this study), MR16-06 (22 August-3 October 2016 as Track 4), MR1704 leg 1 (11 July-2 August 2017 as Track 5), and MR1705C (24 August-29 September 2017 as Track 6).Two other cruises during DJF for the Indonesia-Australia region (5-25 • S, 105-115 • E) and Indonesia-Japan region (10-35 • N, 129-140 • E) are also explored in this study, respectively designated as MR15-05 leg 1 (23 December 2015-10 January 2016 as Track 2) and MR15-05 leg 2 (17-24 January 2016 as Track 3).All measuring data for CO and O 3 from the six cruises are respectively plotted in Fig. 4a and b as grey dots.The ship-based data used in this study was partly reported (T1-4) in the work of Kanaya et al. (2019), including the extraordinary peak of CO on 26 September 2016 exceeded 500 ppbv (off the scale in Fig. 4a-T4) associated with heavy fires in Russia (Kanaya et al., 2019).Data for the North Pacific region (40-60 • N) are addressed in light-blue shades in Fig. 4 (T1, T4-6) for analysis in Sect.3.2.
Table 7 shows correlation coefficients (plotted in Fig. S10), indicating that the CHASER simulations for CO and O 3 are in good agreement with Mirai observations (R = approx.0.6).However, the model still shows some discrepancies for both CO and O 3 concentrations.In general, the estimated CO and O 3 are both reduced for T1 and T4-6 as compared to observations, whereas they are superior for the data located in 20 • S-20 • N during T2-3.Overestimations for CO and O 3 occurring in the region with considerably low levels of these species might be attributed to the lack of halogen chemistry in the model, as also discussed for the nearby region in a past report (Kanaya et al., 2019).Underestimates for O 3 levels up to 70 ppbv in the higher latitudes (Fig. 4b: T1, T4-6) are ascribable to the insufficient downward mixing process of stratosphere O 3 in the model (Kanaya et al., 2019).Except for the CO's peak on 26 September 2016 mentioned above, the mild reductions for CO (< 30 ppbv) in the model during T1 and T4-6 as compared to the observations could be attributed to the insufficient emission from territories, in-Table 5. Model correlations and biases with EANET/EMEP observations: three-sigma-rule outlier detection is applied for each station before calculating all data.For NO x , all data were filtered once more using the two-sigma-rule.Bias of the sensitivity run is shown in bold if it is higher than the bias of the STD run.R has no unit; the units in brackets are for biases.ternational shipping, and aviation, as we used the HTAP-II emission inventory.These reductions for CO in Kanaya's work are minor, e.g., for T1 (MR15-03 leg 1 and 2 in their work), due to reanalysis data with finer horizontal resolution (1.1 • ) utilized for the CO emission rate (Kanaya et al., 2019).

EANET
The negative biases in the noHR simulations for CO are lower in the STD run for all cruises, as they are for the North Pacific region (second versus third/fourth/fifth data rows for CO, Table 7).The CO-increasing effects by N 2 O 5 and HO 2 uptakes (Fig. 4c) are consistent with the comparison for EMEP.This is also true for CO-reduction effects because of HR(RO 2 ).Whereas the effects by N 2 O 5 and HO 2 reduce the model bias, the CO-reducing effects by HR(RO 2 ) exaggerated the CO bias (second versus sixth data rows for CO in Table 7), which is already apparent for comparison with EMEP (last column, Table 5).
For O 3 level, the model underestimates (Table 7) are in the opposite direction to the O 3 overestimates for EANET and EMEP stations (Table 5).The lower panels presented in Fig. 4b show marked O 3 reduction with all HRs (gaps between red and black lines), mostly contributed from the HO 2 uptake onto cloud particles (Fig. 4d: green and hatched fields).This marked reduction of the O 3 level is evident at some points during the cruises, especially in the North Pacific region (the shaded areas) and for T5.Unlike comparisons for land domain data (Table 5), O 3 reduction because of HRs worsens the model underestimates during the Mirai cruises.It is noteworthy that one cannot necessarily confirm whether the STD run simulates these species better than the noHR does because tropospheric CO and O 3 levels are controlled by a complicated chemical mechanism and an interplay of emission, transport, deposition, and local mixing in the boundary layers.As discussed later in Sect.3.2, the surface aerosols concentration in the western Pacific Ocean is mostly dominated by liquid clouds (exceed-ing 50 000 µm 2 cm −3 during JJA) and sulfate aerosols (approximately 75 µm 2 cm −3 in JJA).The model improvements in reproducing CO by adding N 2 O 5 and HO 2 uptake indicate that the appropriate mechanisms of these processes onto cloud droplets and sulfate aerosols are well established in the model.For HR(RO 2 ), which induces the smallest and opposite effects on CO compared with the effects of N 2 O 5 and HO 2 uptakes, it can be stated in general for the total HR effects that including all three HRs partially improves the model during Mirai cruises.

Verification using aircraft measurements
To verify the model performance in the free troposphere, we used ATom1 aircraft measurements in August 2016 (for NO 2 , OH, CO, and O 3 ).Table 8 lists the model's correlation coefficients and biases of each sensitivity run against ATom1.Figure 5 shows observed and simulated concentrations of CO and O 3 during flight no. 2 in the North Pacific (NP) region.In general, the model simulations for NO 2 , OH, CO, and O 3 adequately agree with aircraft measurements with R > 0.5 (Fig. S12).However, NO 2 and CO still tend to be underestimated by the model, which is consistent with comparisons for EANET/EMEP and Mirai observations.In Fig. 5, the CO-increasing effects, mostly due to the uptake of N 2 O 5 and HO 2 (Fig. 6o), mitigated the negative bias of the model.This CO bias reduction was visible for all flight altitudes, the lower troposphere, and the North Pacific region (Table 8; Fig. 6d and h).Both N 2 O 5 and HO 2 uptakes show improvements for CO reproduction of the model.However, RO 2 uptake seems to worsen the model's CO bias due to its reducing https://doi.org/10.5194/gmd-14-3813-2021 Geosci.Model Dev., 14, 3813-3841, 2021 effect for CO (Fig. 6o and s), which is consistent with the Mirai comparison.
For the O 3 level, the model generally overestimates O 3 when calculating for all altitudes or lower troposphere, which is similar to the EANET/EMEP observations.In the North Pacific region with P > 600 hPa (40-60 • N, 198-210 • W), the model bias for O 3 in STD run turns to underestimate (second data row and second column from the right, Table 8), which might be similar with Mirai data verification for the western North Pacific (143 • E-193 • W).However, the underlayers (> 700 hPa) again show overestimation (second data row -last column, Table 8).As Mirai and ATom1 data show, the underestimates for O 3 exist at the marine boundary layer in the western North Pacific and extend to the upper troposphere (< 700 hPa) of the east side, which might be ascribed to the insufficient downward mixing process of stratosphere O 3 in the model as discussed previously.
The HR effects on O 3 are generally negative effects (Fig. 6m), although they are small and barely recognizable in Fig. 5, which mitigates the model bias in the noHR run.This model improvement is consistent for all flight altitudes, the low troposphere, and the North Pacific region (second versus third data rows in Table 8).Both HR(N 2 O 5 ) and HR(RO 2 ) typically contribute to this improvement (Fig. 6m and q).In contrast, HR(HO 2 ) seems to only reduce the model bias in a thin layer: from the ground up to 800 hPa for all flights and 700 hPa for the North Pacific region (Fig. 6m and q).At the bottommost layers in this region, the model's overestimates for O 3 are reduced by the negative effects of HO 2 uptake (Fig. 6f and q).The extension of model bias because of HO 2 uptake above 800 hPa is attributable to its increasing effect on O 3 level (Fig. 6m).We recognize that this O 3 increase effect above 800 hPa is opposite to the effects for EANET/EMEP and Mirai comparisons, which is discussed in Sect.3.2 for HO 2 uptake effects.
The vertical means of model biases and changes for all four species (NO 2 , OH, CO, O 3 ) are presented in Fig. 6.In general, the STD run reduces model biases for all four species, with better performance for broader regions (all flights) than for the smaller region (North Pacific).In the North Pacific region, the negative bias for O 3 is observed only for the 500-900 hPa layers (Fig. 6f).The model bias is apparently extended in this region.However, the inclusion of HR(HO 2 ) leading to O 3 increment (Fig. 6f) reduces O 3 bias in this region, which might indicate that the O 3 increase effect by HR(HO 2 ) is verified, particularly in 500-900 hPa layers during ATom1.
We also verify the total uptake of N 2 O 5 , HO 2 , and RO 2 onto ice and liquid clouds using data obtained from ATom1 flights within the free troposphere.As Table 8 shows, the inclusion of HRs onto clouds reduces the model biases for CO and O 3 in all calculations.In Fig. 6l-s, because the HRs(cloud) effects occupy the major part of total HRs effects for NO 2 , O 3 , and OH, especially for NP region and low troposphere (> 800 hPa), cloud uptakes could also contribute to the overall reduction in model bias against ATom1.For O 3 , HRs(cloud) mostly induce negative effects (Fig. 6m  and q).At the layer of 600-800 hPa in the NP region, this O 3 reduction due to HR(RO 2 ) onto clouds (Fig. 6q: yellow and hatched patterns) might account for the model worsening in O 3 levels as described above.This result might prove that cloud overestimation for the North Pacific, as revealed at the beginning of this section, affects the model bias in this region.

Verification with OMI satellite observation for TCO
We also tested STD and noHR simulations using the tropospheric column ozone (TCO) derived from the OMI satellite instrument (Figs.S13 and 7).In a large area of the Northern Hemisphere, the inclusion of HRs (STD run) generally improved the consistency with the OMI TCO, as seen in Fig. 7b, d, f with less overestimation than Fig. 7a, c, e, particularly enhancing the winter minima (first and second panels in Fig. S13).This improvement in DJF is attributed mostly to the reductive effects of HR(N 2 O 5 ) and HR(RO 2 ) in the lower (800 hPa) and middle troposphere (500 hPa), respectively (see Fig. 10 for vertical profiles of HR(N 2 O 5 ) on O 3 and Fig. 13 for vertical profiles of HR(RO 2 )).In the North Pacific, HRs appeared to exaggerate O 3 underestimates, especially for latitudes higher than 40 • N (Fig. 7) during the first half of the year (third panel, Fig. S13).However, such a discrepancy, which was also observed from comparison for R/V Mirai observations (Fig. 4), might result from other factors such as deposition or vertical mixing rather than by HRs.

HR effects
This section presents a discussion of the global effects of HRs calculated using CHASER with their spatial distributions in the troposphere using standard (STD) and sensitivity simulations (noHR_n2o5, noHR_ho2, noHRs_ro2, noHR) for the meteorological year of 2011.Aside from the main simulations described in Table 3, additional runs that separately turned off the uptakes onto clouds or aerosols for each HR are also conducted to exploit the contributions of effects to the troposphere.

Distribution of clouds and aerosols surface area density (SAD)
To obtain the parameters for uptake to clouds and aerosols, SAD estimations are used together with cloud fraction and https://doi.org/10.5194/gmd-14-3813-2021Geosci.Model Dev., 14, 3813-3841, 2021  aerosols concentration.Hereinafter, we discuss SAD distributions for total aerosol, ice clouds, and cloud droplets, which are estimated for the model using Eqs.( 2), (3), and (4), respectively.In Fig. 8, total surface area concentrations of liquid clouds and aerosols are both much lower aloft than at the surface (as counted on the dry and wet depositions).The liquid cloud SAD results are 2 orders of magnitude larger than ice cloud SAD and total aerosol SAD.The ice cloud SAD, distributed at the middle and upper troposphere, is enhanced for N/S middle latitudes in wintertime.Liquid cloud SAD concentrates mainly at the surface with distributions extending to 500 hPa and maximized at approx.800 hPa over the midlatitude storm tracks and in tropical convective systems, especially at 60 • N in JJA.Total aerosol SAD was derived mainly from pollution sources at 40 • N during both seasons, with higher concentrations apparent for DJF and a greater spatial spread observed for JJA.Sulfate aerosols are becoming the dominant source of aerosol surface area in the model above 600 hPa (approx.20 µm 2 cm −3 ) in addition to organic carbons and soil dust (both are approx.10 µm 2 cm −3 in JJA) for the Northern Hemisphere.
In Fig. S14, showing the SAD distribution at the surface, the SAD for liquid clouds is dominant in JJA, reaching approx.50 000 µm 2 cm −3 for mid-latitude and high-latitude ocean regions.Liquid clouds make the greatest contribution to the SAD at the surface.Our model performance for aerosol SAD shows agreement with that presented in an earlier report (Thornton et al., 2008).Sulfate aerosols are prevalent in the northern mid-latitudes near industrial bases; maximize at the surface in DJF for the Chinese region (exceeding 1000 µm 2 cm −3 ), NE U.S. (approx.500 µm 2 cm −3 ), and western Europe; and transport to the North Pacific region in JJA (approx.250 µm 2 cm −3 ).Soil dust aerosol SAD dominate in the regions of the Sahara and Gobi deserts, reaching annual average values exceeding 100 µm 2 cm −3 .Organic carbon (OC) is a dominant source of aerosol SAD over biomass burning regions in China (up to 1000 µm 2 cm −3 in DJF), South Africa (up to 800 µm 2 cm −3 in JJA), western Europe, and South America.The black carbon (BC) surface area can reach values exceeding 600 µm 2 cm −3 in DJF for the region of China or other significant industrial areas (India, which reaches 75 µm 2 cm −3 , NE U.S., and Europe) or over tropical forests, primarily in Africa.Sea salt aerosols are most important in high-latitude oceans during winter.However, the maximum contributions only reach 2 µm 2 cm −3 in our model, which is greatly underestimated compared to Thornton's work (75 µm 2 cm −3 ) (Thornton et al., 2008).In brief, SAD for aerosols of all types contributes the most during DJF, whereas during JJA, the SAD for liquid clouds and sulfate aerosols are dominant, particularly for the northern highlatitude and mid-latitude oceans.The total aerosol SAD in  this region is approx.75 µm 2 cm −3 , which is consistent with the estimation by Thornton et al. (2008).

Effects of N 2 O 5 heterogeneous reaction (HR(N 2 O 5 ))
The inclusion of HR(N 2 O 5 ) in the model increases global methane lifetime by +4.48 % and changes NO x , O 3 , and CO abundances by −5.51 %, −2.12 %, and +3.42 %, respectively (Table 9).In Fig. 9a-h, the changes in OH, NO x , O 3 , and CO are most significant in the middle troposphere (400-600 hPa).These changes are attributed mostly to uptakes of N 2 O 5 onto aerosols, rather than onto clouds.Marked negative effects on NO x concentration are apparent for DJF in the middle troposphere (600-700 hPa) of the 60 • N and the Arctic region (> −20 % at 700 hPa) (Fig. 9c).The effects are probably associated with high concentrations of sulfate aerosols, organic carbons, or soil dusts in the middle troposphere (see the paragraph above) and are also related to a long chemical lifetime of NO y in the middle and upper troposphere in winter.When it comes to JJA, these negative effects become significant at higher altitudes around 30 • N/S (> −10 % at 400 hPa).At the surface (Fig. 9i-p), HR(N 2 O 5 ) causes negative effects on NO x , O 3 , and OH concentrations (up to about −24 %, −5 % and −8 % respectively) and positive effects on CO concen-tration (up to about +4 %), which are also mainly attributable to the N 2 O 5 uptake on aerosols.
In Fig. 10, the latitude-longitude (lat-long) means of HR(N 2 O 5 ) effects are calculated for each pressure range (pressure ranges are defined as in Fig. 6).The global NO x decrease is up to −9 % at 300-400 hPa.This decrease causes correspondent reductions in O 3 and OH, which are calculated as about −3 % and −7 % at 400-600 hPa, respectively, for lat-long mean O 3 and OH.About 4 % lat-long mean CO increment throughout the entire troposphere responds to decreased OH.
The small effects of HR(N 2 O 5 ) on O 3 in the lower troposphere are consistent with findings from an earlier study (Riemer et al., 2003).Reductions in O 3 and NO x concentrations also agree well with the collective knowledge summarized in work reported by Brown and Stutz (2012).Despite a considerable HR(N 2 O 5 ) effect calculated in the middle troposphere, its effect in the whole troposphere is apparently not as great as reported to date.Another study assessed HR(N 2 O 5 ) effects on annual burdens of NO x , O 3 , and OH as −11 %, −5 %, and −7 %, respectively, when using a similar γ N 2 O 5 value (0.1) (Macintyre and Evans, 2010).Although the effects of magnitude estimated in our work (Table 9) are almost half those of this earlier study (probably because of differences in NO x emissions, estimation of SAD, and chemical mechanism), the effect tendencies are similar.A strong increase of ozone attributed to N 2 O 5 uptake under high-NO x  conditions calculated using box models was reported from an earlier study (Riemer et al., 2003), but this is only slightly apparent in our global model.Our results revealed that the HR(N 2 O 5 ) effect might help clean up NO x pollutant.However, it increases the concentration of other pollutants (such as CO) because of the effects of reducing oxidizing agents in the atmosphere.
https://doi.org/10.5194/gmd-14-3813-2021Geosci.Model Dev., 14, 3813-3841, 2021 As Fig. 11a-f shows, the zonal-mean effects of HR(HO 2 ) on NO x , OH, and O 3 are more widespread in DJF but are more concentrated at the surface in JJA because of the high level of HO 2 .The most substantial effects by HR(HO 2 ) are calculated in JJA at the surface of North Pacific (140-240 • E, 40-60 • N) by as much as about +69 % (NO x ), +7 % (CO), −70 % (OH), and −21 % (O 3 ), which are more significant than those of HR(N 2 O 5 ) at the surface.In the lower troposphere, HR(HO 2 ) suppresses the NO oxidation (Reaction R10) and thus preserves a high NO/NO x ratio and gen-erally restricts OH and O 3 formations.
These effects are primarily attributable to HR(HO 2 ) in clouds rather than to aerosols (which is opposite to N 2 O 5 uptake).These OH and O 3 reduction effects go along with past studies in which approx.50 % OH and approx.10 % O 3 of reductions are calculated for the low troposphere of northern mid-latitude region ascribed to aqueous-phase HO x sink in clouds (Lelieveld andCrutzen, 1990, 1991).The efficient scavenge of HO 2 radical by cloud droplets might associate with acid-base dissociation HO 2 /O − 2 and electron transfer of O − 2 to HO 2 to produce H 2 O 2 (Jacob, 2000).Furthermore, cloud droplets SAD in our model are 2 orders of magnitude higher than total aerosol SAD (Fig. 8), which also contributes to the preference of the aqueous-phase HO 2 sink.Our large calculated effects for the North Pacific region are new findings from other models, which have considered only aqueous aerosols (Stadtler et al., 2018;Thornton et al., 2008) because cloud particles are dominant at remote marine areas in addition to sulfate and aqueous sea salt particles (discussed at the beginning of Sect.3.2).The HO 2 uptake onto aerosols is minor; it is observed only in DJF in the Arctic region and polluted areas (China and US), with apparent changes of up to +17 % for NO x , −40 % for OH, and −14 % for O 3 at the local surface (Fig. 11i, k, m).The aerosol negative effect of HR(HO 2 ) on surface O 3 concentration is significant in the Chinese area, which might be in line with other studies of the Chinese O 3 trend (Kanaya et al., 2009;Li et al., 2019;Liu and Wang, 2020;Taketani et al., 2012), which suggests that the observed recent O 3 increases can be attributed in part to reduced HO 2 uptake under aerosol (PM) decreases brought about by the new Chinese air pollution policy.
In Fig. 10, vertical profiles show that the lat-long averaged effect of HR(HO 2 ) on OH is −9 % in the lower troposphere.As a result, the lat-long mean CO level increases by 2 % at the surface.Additionally, the daytime NO x oxidation by OH is suppressed.Also, NO x might be preserved in clouds (Dentener, 1993), which increases the lat-long averaged NO x level by +6 % at 900 hPa.The lat-long mean O 3 is reduced by −1 % at the surface, but it is increased at higher altitudes (about +0.2 % at 300 hPa).The reduction of O 3 is associated with HO 2 depletion in clouds and aqueous aerosols as described above, coupled with the NO x preservation in clouds, which enhances the NO/NO x ratio.The preserved NO x in clouds might remain available for O 3 production after the cloud evaporates (Dentener, 1993), along with the low SAD for both liquid clouds and aerosols at higher altitudes (Fig. 8), thereby increasing O 3 in places other than the aqueous phase.The O 3 increment might be trivial in DJF but enhanced in JJA.As a result, the Northern Hemisphere mean O 3 in JJA exhibits only positive effects.In contrast, for the North Pacific region in JJA, due to its large cloud fraction, an O 3 reduction effect is apparent in this region.The effects in JJA for this region show changes of −25 % OH, +35 % NO x , −12 % O 3 , and +5 % CO at 900-100 hPa as the most remarkable HR(HO 2 ) effects, as described above.In general, the regional mean effects of HR(HO 2 ) in the North Pacific region are enhanced in JJA, but the mean global effects of HR(HO 2 ) are slightly favored in DJF because of the additional effects of aerosols during this season.Macintyre and Evans (2011) also found a similar contrast between the behaviors of HR(N 2 O 5 ) and HR(HO 2 ): the uptake of N 2 O 5 produces both regional and global effects on O 3 , whereas the uptake of HO 2 affects O 3 at regional scales more strongly than on a global scale (Macintyre and Evans, 2011).Such features are generally consistent with our results.

Effects of RO 2 heterogeneous reactions (HR(RO 2 ))
Effects of HR(RO 2 ) increase the global mean methane lifetime by +0.15 % and change tropospheric abundances of NO x (+0.52 %), O 3 (−0.93%), and CO (−1.78 %) (Table 9).In Fig. 12k-l, significant latitudinal contrasts exist in the NO x changes: large NO x increases at high latitudes and decreases at lower latitudes.These NO x changes probably reflect the reduced formation of PANs, which decreases NO x transport from source regions to remote areas and from the surface to the upper troposphere (Villalta et al., 1996).The model calculated especially large NO x increases (> 50 %) for high latitudes around the Arctic Ocean in JJA, indicating a reduction in the formation of PANs (NO 2 +RO 2 → PANs), which is linked tightly to the enhanced biogenic emissions of VOCs such as isoprene and terpenes in summer.The PANreducing effect of HR(RO 2 ), associated with the suppression of NO oxidation, as seen in the case of HR(HO 2 ), causes a double increase for NO x at the surface as compared to HR(HO 2 ) (144 % in Fig. 12l versus 66 % in Fig. 11l).The increases in NO x in DJF at high latitudes offshore of southern oceans could also reflect reduced transport of NO x under reduced PAN formation caused by HR(RO 2 ) since these coastal areas are located downwind of the major biogenic volatile organic compounds (BVOC) sources of South America, South Africa, and Australia.Moreover, the areas with pronounced NO x increases in Fig. 12k-l are all associated with high-cloud SAD, as seen in Fig. S14 (left panels).
In the troposphere, the NO oxidation carried out by RO 2 via Reaction (R11) produces another HO x molecule (Reaction R12).Hence, the uptakes of RO 2 onto particles are also expected to confine the OH and O 3 formations, similar to what the uptake of HO 2 does.
However, the different catalytic role of NO x in the oxidizing mechanisms between polluted and remote regions results in different tendencies for OH and O 3 (Jacob, 1999).Accordingly, O 3 and OH molecules are produced through the oxidation of hydrocarbons in presence of high NO x , whereas HO x molecules are consumed and no O 3 is produced in the same process that occurred devoid of NO x .For that reason, less RO 2 participating in the hydrocarbon oxidation only reduces OH and O 3 levels in polluted regions, e.g., the Chinese region in DJF (Fig. 12i and m), while enhancing OH level and leaving no significant effect on O 3 at remote regions, e.g., the NP region (Fig. 12i, j, m, n).The corresponding changes in OH concentration at the surface are in the range of about −5 % to +20 % in JJA (Fig. 12j).O 3 levels are slightly reduced throughout the troposphere, at most about −5 % at the surface of the Arctic region above Canada, due to HR(RO 2 ), with only about a 3 % increase, corresponding with the NO x https://doi.org/10.5194/gmd-14-3813-2021Geosci.Model Dev., 14, 3813-3841, 2021   increases at the tropical coasts (Fig. 12n).For CO, different from the increasing impacts by HR(N 2 O 5 ) and HR(HO 2 ) (+3.42 % and +1.95 %, respectively), a reduction effect occurs through the free troposphere (−1.78 %).CO decrease might be due to reduced secondary CO production from RO 2 oxidation (RO 2 → HCHO/RCHO or ROOH → CO) such as isoprene (Bates and Jacob, 2019) when functionalized RO 2 species uptake onto aerosols and clouds particles instead.The effects of HR(RO 2 ) are primarily attributed to the heterogeneous reaction on clouds rather than on aerosols in terms of changes in NO x and CH 4 lifetime, although this cloud effect is far smaller than the cloud effect on the HO 2 uptake.The areas with pronounced NO x increases in Fig. 12b are all associated with high-cloud SAD, as seen in Fig. S14 (left panels).Although it is proper to expect the high solubility of RO 2 (e.g., CH 3 O 2 ) from its peroxy substituent (Betterton, 1992;Shepson et al., 1996), it is much less soluble than HO 2 because of its lower polarity and thus lower Henry law constant (Jacob, 2000).Consequently, the possible accumulation of CH 3 O 2 in the cloud is attributable to suppression of its gas-phase sink with HO 2 (Jacob, 1986).
Figure 13a-c show lat-long means of HR(RO 2 ) effects calculated for the respective pressure ranges: the lat-long values are constrained for the entire globe, the Northern https://doi.org/10.5194/gmd-14-3813-2021 Geosci.Model Dev., 14, 3813-3841, 2021 Hemisphere, and North Pacific region.For the entire globe, the contrast effects of HR(RO 2 ) between the lower and higher troposphere on NO x and OH are shown clearly (+3.50 % NO x and +0.55 % OH at 900 hPa, but −2.50 % NO x and −0.75 % OH at 400-500 hPa annually).As a result, the annual mean O 3 and CO levels decreased throughout the troposphere, reaching their lowest at −1.60 % O 3 and −1.50 % CO at the surface.In JJA, the global effects by HR(RO 2 ) are more concentrated in the lower troposphere, especially in the North Pacific (+3 % OH, +10 % NO x , −3 % O 3 , −2 % CO at 900-1000 hPa).In DJF, the HR(RO 2 ) effects are observed mostly in the middle and higher troposphere, especially when considering the Northern Hemisphere (−1.25 % OH, −4 % NO x , −2 % O 3 at 500-800 hPa).

Total effects of all HRs
As discussed above, different heterogeneous reactions affect tropospheric chemistry differently.However, their effects can either augment or negate others in performing for the atmospheric chemistry.HR(N 2 O 5 ) is the greatest contributor to reduction of tropospheric OH, O 3 , and NO x abundances and is more active in the middle troposphere.HR(HO 2 ) reduces OH but increases the abundances of O 3 and NO x globally, whereas it exposes a negative effect on O 3 level at the surface of the North Pacific region.HR(RO 2 ) similarly has a smaller distribution to the total heterogeneous effects, but its global mean negative effects for O 3 are not negligible.The uptake of N 2 O 5 mainly takes place for aerosols, whereas the uptakes of HO 2 and RO 2 occur more for liquid and ice clouds.Overall, the total effects of all HRs for the whole troposphere are +5.91 % for global mean CH 4 lifetime, −2.19 % for NO x (tropospheric abundance), −2.96 % for O 3 , and +3.28 % for CO (Table 9).At the surface, the annual effects ranged from about −53 % to +2 % for OH, −13 % to +51 % for NO x , −13 % to −2 % for O 3 , and −0.3 % to +6 % for CO (Fig. 14).As Fig. 13d-f show for the vertical profiles of HR effects, the change of OH largely concentrated in the lower troposphere (−10 % OH at 900 hPa, calculated for the entire globe) is associated with the HO 2 uptake.By contrast, the NO x change is more intensive at higher altitudes (−9 % NO x at 400 hPa, calculated for the entire globe), associated with N 2 O 5 and RO 2 uptake.The global mean HR effects on O 3 and CO are vertically even, with the highest effects reaching −4 % O 3 and +4 % CO at the surface.Globally, HR effects on atmospheric oxidants (OH and O 3 ) are enhanced in DJF because of the higher pollution in the Northern Hemi- sphere.However, the largest HR effects are apparent for JJA at the surface of the North Pacific (−25 % OH, +38 % NO x , −14 % O 3 , and +6 % CO as calculated for the 950-1000 hPa layer).These effects are mostly ascribed to HO 2 uptake onto clouds.This finding is also apparent from Fig. S15b: these effects reach −66 % for OH, +206 % for NO x , −23 % for O 3 , and +4.4 % for CO at the surface.They were able to extend up to 400 hPa in the atmosphere.These substantial effects are readily apparent for the large reduction of O 3 level during Mirai observation (hatched field in Fig. 4d, T5).However, the major contribution of HR(HO 2 ) to these effects is only partially verified by the ATom1 measurements in this study (Fig. 6q, hatched bars at 500-700 hPa).Because of model overestimates of cloud fraction in JJA for the North Pacific region, these effects of HR(HO 2 ) should have existed at a smaller magnitude.For HR effects in the middle to upper troposphere, the N 2 O 5 uptake on aerosols is dominant in these layers and intense in both DJF and JJA.

Sensitivities of tropospheric chemistry respond to heterogeneous reactions
From the discussion presented above, marked effects of HRs on global tropospheric chemistry are apparent.Here we examine how the tropospheric chemistry responds to the mag-nitude of HRs' loss rates.To do this, we introduced a factor F for application to the first-order loss rate shown in Eq. ( 1) for artificially manipulating the HR magnitude.
For this sensitivity test, we only specifically examine HR(HO 2 ) and HR(N 2 O 5 ) and consider factors of 0-10 for the STD run (Table S1).This test might help to show the effective oxidation sensitivity of the troposphere because future pollution and climate change might enhance the activities of these HRs.
For both effects, we performed nonlinear function fitting with their uptake loss rates, which yielded correlation coefficients higher than 0.93 (Fig. 15).Although both HRs showed negative tendencies for OH and O 3 levels, the effect of HR(HO 2 ) on the tropospheric abundance of O 3 showed only a small increment with an increasing loss rate (maxima at around F = 3) and turned to reduction at higher rates (F > 5).As discussed alongside HR(HO 2 ) effects, the O 3 level is expected to be primarily reduced only in JJA at the surface of the North Pacific region.At the same time, O 3 will be increased gradually elsewhere because of the persistent NO x increment.This behavior produces a positive global https://doi.org/10.5194/gmd-14-3813-2021 Geosci.Model Dev., 14, 3813-3841, 2021 mean effect.Figure 10 (dashed lines) shows that manipulation of the HR(HO 2 ) loss rate that is 10 factors higher will effectively increase the negative HR(HO 2 ) effects on O 3 in DJF (dashed blue versus solid blue lines, third row, fourth column), which results in a higher tendency of negative values for global mean effects.This sensitivity in DJF might be attributable to the HO 2 uptake to aerosols rather than to clouds during this polluted period, which is apparent through comparison of Figs.11 and S16 for notable events.In DJF, by amplifying the HO 2 uptake loss rate by a factor of 10, the effects for the polluted Chinese area (because of HO 2 uptake onto aerosols) significantly magnify from −18 % (Fig. 11m) to −47 % (third row, first column in Fig. S16b).In contrast, effects at the surface O 3 level in JJA for the North Pacific region (because of HO 2 uptake onto clouds) only enhance from −21 % (Fig. 11n) to −29 % (third row, second column in Fig. S16b).
When amplifying HR(N 2 O 5 ) by a factor of 10, the sensitivities of global effects show no seasonal variation.The HR(N 2 O 5 ) effects are more sensitive in DJF for the North Pacific region, which link to the higher concentration of aerosol in this season.Otherwise, the HR(N 2 O 5 ) effects for the generic Northern Hemisphere tend to be more sensitive in JJA as a result of pollutant transportation to the higher troposphere.
Consequently, we suggest that the sensitivity of tropospheric chemistry to HR(N 2 O 5 ) and HR(HO 2 ) might be attributable to loss activities of aerosols rather than to clouds.The sharply curved effect on O 3 because of amplification of HR(HO 2 ) makes sense in plans for ozone pollution control when increased pollution or climate change factors cause the rate of HRs to increase by 3-5 times or more in the future.

Conclusion
The "CHASER" chemistry-climate model was used to investigate global effects of N 2 O 5 , HO 2 , and RO 2 uptake.Verification of the model with observations from inland and ocean domains showed adequate agreement for PM 2.5 , SO 2− 4 , and NO − 3 particles and gaseous HNO 3 , NO x , OH, CO, and O 3 concentrations.R, bias, and NRMSE values for SO 2− 4 , NO − 3 , and HNO 3 at EANET and EMEP stations are comparable with other models.Inclusion of HRs reduced model bias for OH, NO 2 , CO, and O 3 , especially in the low troposphere.However, verification with satellite and reanalysis data showed deterioration by HRs for TCO and an overestimate for cloud fraction in the North Pacific region.
The total effects of HRs are important for the tropospheric chemistry that might change +5.91 % CH 4 lifetime and −2.19 % NO x , −2.96 % O 3 , and +3.28 % CO abun-dances.Global effects are −9 % NO x at 400 hPa, −10 % OH at 900 hPa, and −4 % O 3 and +4 % CO at the surface.Global HR effects tend to be enhanced in DJF because of greater amounts of pollution in the Northern Hemisphere.
Total HR effects are contributed mainly by HR(N 2 O 5 ) onto aerosols in the middle troposphere.At the surface, HR(HO 2 ) is more active and leaves a remarkable disturbance in JJA at the North Pacific region with changes of −70 % for OH, −24 % for O 3 , +68 % for NO x , and +8 % for CO.These effects were attributed to the uptake of HO 2 on cloud particles, which were partially verified with ATom1 observations.HR(RO 2 ), which also favors cloud particles, minorly contributes to the tropospheric chemistry, but it has enormous impacts on PAN and NO x transportation (+144 % NO x for the North Pacific and Atlantic regions in JJA) and the negative changes in CO (−1.78 %), as compared to positive effects by HR(N 2 O 5 ) and HR(HO 2 ), that can not be neglected.However, the effect magnitude requires further investigation because of model overestimates for cloud fractions in this region.
The sensitivity of tropospheric chemistry to HR magnitude was determined as a nonlinear function.The increasing effect on the global O 3 abundance from HR(HO 2 ) will sharply change to a decreasing effect when the uptake rate is amplified by more than 3 times.This turning is ascribed to the uptake onto aerosols in DJF.In general, uptake to aerosols is more responsive to the heterogeneous loss rate than uptake to clouds.
Overall, the N 2 O 5 and HO 2 uptakes will sweep away atmospheric oxidants, thereby enhancing concentrations of pollutants.Our results reveal that although HRs are reported to be associated with polluted regions, the global effects of HRs reach further remote regions such as the marine boundary layer at mid-latitudes and the upper troposphere.For ground-based studies of polluted regions such as China, it should be considered that HR(HO 2 ) and HR(RO 2 ) were able to contribute to the NO x increment in DJF and JJA, respectively.Moreover, the HR(HO 2 ) effect might hinder efforts at reducing environmental pollution in urban areas because it increases NO x but decreases O 3 at the surface.Therefore, if this reaction is minimized because of a decrease in particulate matter, then the surface ozone level might increase.

Figure 3 .
Figure 3. Monthly averaged concentrations at EANET and EMEP from 2010-2016 (a-g) and corresponding HR effects (h-n) for NO x : (a, h) Jinyunshan (China), (b, i) other low-NO x EANET stations, and (c, j) EMEP stations; for O 3 : (d, k) Korean stations, (e, l) other EANET stations, and (f, m) EMEP stations; for CO: (g, n) EMEP stations.In (a)-(g), grey lines are observations (OBS), red lines are the STD simulation, and black lines are the noHR simulation.All lines showed the median of monthly mean concentrations for each group of stations, except (a).In (b)-(g), vertical thin lines with markers show the 25th-75th percentiles of monthly mean concentration at the particular group of stations.In (h)-(n), blue fields are changes caused by HR(N 2 O 5 ), green fields are changes in HR(HO 2 ), and yellow fields are changes caused by HR(RO 2 ).

Figure 6
exhibits vertical biases of the model and computed HR effects caused by each HR for all ATom1 flights and NP region.The spatial and temporal concentrations are available in Fig. S11.Correlations are shown in Fig. S12.

Figure 4 .
Figure 4. Observed and simulated concentrations for CO (a) and O 3 (b) and their changes caused by each HR (c, d) during Mirai cruises.The left axis shows concentrations.Dashed lines show latitudes of cruises scaled with the right axis.The horizontal axis shows cruise travel times.Light-blue areas show data in the North Pacific region (140-240 • E, 40-60 • N).In (a) and (b), grey dots are for observations, red lines are for STD simulations, and black lines are for noHR simulations.In (c) and (d), blue fields are for changes in HR(N 2 O 5 ), green fields for changes in HR(HO 2 ), yellow fields for changes in HR(RO 2 ), and hatched fields for changes in HRs (cloud).Note that in (d) all three fields are stacked but in (c) only blue and green fields are stacked for better illustration of negative and positive changes.

Figure 6 .
Figure 6.Vertical bias against ATom1 (a-h) and vertical HR effects (l-s).Data for each pressure level P are calculated within the range of P ±50 hPa, with the applied three-sigma-rule for outlier detection.All rows show calculations for all flight domains the and the North Pacific region.The horizontal axis shows model bias and absolute changes with units written in each panel.The vertical axis shows pressure (hPa).The red numbers in (a)-(h) represent relative reductions (%) of the STD run's bias compared to that of the noHR run.

Figure 7 .
Figure 7. Differences in tropospheric column ozone (TCO) in CHASER's noHR (a, c, e) and STD (b, d, f) runs from OMI: ANN is annual; DJF is December, January, February; and JJA is June, July, and August.

Figure 8 .
Figure 8. Zonal, seasonal mean (upper and middle), and annual mean (lower) distribution of cloud droplet surface aerosol density (SAD), cloud ice SAD, total aerosol SAD, sulfate aerosol SAD, and organic carbon SAD (from left to right).

Figure 9 .
Figure 9. Effects by N 2 O 5 uptake onto both clouds and aerosols in zonal mean (a-h) and at the surface (i-p).Note that the color scale for (a)-(h) is different from that for (i)-(p).

Figure 10 .
Figure 10.Averaged effects of N 2 O 5 (a-c) and HO 2 uptakes (d-f) for each air pressure range.Calculations are for global data (a, d), Northern Hemisphere data (b, e), and North Pacific region data (c, f).Dashed lines show effects of the FCTHR_10 run (Sect.3.3).The legend is in panel (f) for December-January-February (DJF), June-July-August (JJA), and annual (ANN) data.

Figure 11 .
Figure11.Effects of HR(HO 2 ) on both clouds and aerosols in zonal mean (a-h) and at the surface layer (i-p).Note that the color scale for (a-h) is different from that for (i-p).

Figure 12 .
Figure 12.Effects of HR(RO 2 ) on both clouds and aerosols in the zonal mean (a-h) and at the surface layer (i-p).Note that the color scale for (a)-(h) is different from that for (i)-(p).

Figure 13 .
Figure 13.Averaged effects of HR(RO 2 ) (a-c) and all HRs (d-f) on OH (first row), NO x (second row), O 3 (third row), and CO (fourth row) for each air pressure range.Calculations are for global (a,d), Northern Hemisphere (b, e), and North Pacific region (c, f).Legend is on the first panel for December-January-February (DJF), June-July-August (JJA), and annual (ANN).

Figure 14 .
Figure 14.Annual zonal mean (a) and surface (b) total HR effects.

Figure 15 .
Figure 15.Trend lines for the sensitivity of HR(N 2 O 5 ) (a, b) and HR(HO 2 ) effects (c, d), with uptake rates shown.Panels (a) and (c) show the CH 4 lifetime (blue) and tropospheric abundance of CO (red).Panels (b) and (d) show tropospheric abundances of NO x (blue) and O 3 (red).

Table 1 .
Computation packages in the chemistry-climate model CHASER.

Table 4 .
Datasets used for verification in this study.

Table 6 .
Comparisons between EANET and EMEP observations with atmospheric models.Outlier detection follows that in Table5.The model result is shown in bold if it is better than or agreed with Bian's report.R has no unit; units in brackets are for biases and NRMSE.

Table 7 .
Model correlations and biases for Mirai.No outlier filtration is applied.The bias of the sensitivity run is shown in bold if it is higher than the bias of STD run.R has no unit; the units in brackets are for biases.

Table 8 .
Model correlations and biases with ATom1: three-sigma-rule was applied for CO and O 3 .NP denotes North Pacific region (140-240 • E, 40-60 • N).The bias of the sensitivity run is presented in bold when it is higher than the bias of STD run.R has no unit; the units in brackets are for biases.hPa) (> 600 hPa) (> 800 hPa) (NP, > 600 hPa) (NP, > 600 hPa) (NP, > 700 hPa)