Sensitivity of tropospheric ozone to halogen chemistry in the chemistry–climate model LMDZ-INCA vNMHC

. The atmospheric chemistry of halogenated species (Cl, Br, I) participates in the global chemical sink of tropo-spheric ozone and perturbs the oxidising capacity of the tro-posphere, notably by inﬂuencing the atmospheric lifetime of methane. Global chemistry–climate models are commonly used to assess the global budget of ozone and its sensitivity to emissions of its precursors, as well as to project its long-term evolution. Here, we report on the implementation of tropospheric sources and chemistry of halogens in the chemistry–climate model LMDZ-INCA (Laboratoire de Météorologie Dynamique general circulation model, LMDZ, and Interactions with


Introduction
The atmospheric chemistry of halogenated species was first studied for its importance for the composition of the stratosphere and in particular the ozone layer (Molina and Rowland, 1974;Stolarski and Cicerone, 1974), but its importance for tropospheric chemistry was also suggested as soon as the early 1980s (Chameides and Davis, 1980). With the progress of measurement techniques, marine aerosols were thought to harbour and release reactive halogens into the gasphase (Cicerone, 1981). Later on, halogenated species were found to be responsible for ozone destruction in the Arctic troposphere (Bottenheim et al., 1990;Oltmans and Komhyr, 1986). Meanwhile, natural and anthropogenic sources of halogens were starting to be identified and characterised (Graedel and Keene, 1995;Lobert et al., 1999;McCulloch et al., 1999;Carpenter, 2003;Platt and Hönninger, 2003;Saiz-Lopez and Plane, 2004). In recent years, the addition of halogenated chemistry into global chemistry transport models has revealed its significant influence on the photo-oxidising chemistry of the troposphere on a global scale (Saiz-Lopez et al., 2012b;Simpson et al., 2015). This influence arises from the high reactivity of atomic halogens (e.g. Cl, Br, I) and halogen oxide radicals (e.g. ClO, BrO, IO and higher oxides), which are produced in the atmosphere by the decomposition Published by Copernicus Publications on behalf of the European Geosciences Union.
of very short-lived halogenated substances (VSL) and by heterogeneous reactions on sea salt aerosols and which, in the case of iodine, are also released via ozone deposition to the ocean surface. They can also be generated by polar-specific halogen activation mechanisms on sea ice surfaces, blowing snow and the snowpack (as reviewed by Abbatt et al., 2012).
Details on the tropospheric chemistry of chlorine, bromine and iodine can be found in review articles (e.g. Saiz-Lopez and von Glasow, 2012;Simpson et al., 2015). Reactive halogen species mainly affect tropospheric chemistry by (1) modifying the partitioning of HO x (HO 2 + OH) and nitrogen oxides (NO x = NO 2 + NO) and reacting with NO 2 to form species that remove NO x from the system via deposition or other loss processes, (2) being involved in the oxidation of non-methane volatile organic compounds (NMVOCs) and mercury, (3) participating in the formation of new particles (Giorgi and Chameides, 1985;O'Dowd et al., 2002;von Glasow et al., 2004;Gómez Martín et al., 2020;Saiz-Lopez and Von Glasow, 2012), and (4) constituting a sink for tropospheric ozone at the global scale. Halogenated species play a role in the ozone and hydroxyl radical (OH) budgets and therefore influence the oxidative capacity of the troposphere (Iglesias-Suarez et al., 2020;Badia et al., 2021). In addition, by influencing the abundances of tropospheric ozone, which is an important greenhouse gas, halogens can indirectly affect the climate (Saiz-Lopez et al., 2012a;Hossaini et al., 2015;Sherwen et al., 2017).
The loss of tropospheric ozone by halogens is mainly induced by iodine and bromine species with a lower contribution of chlorine species (Sherwen et al., 2016b). It can be caused by direct catalytic ozone destruction and reduction in ozone production through the removal of NO 2 by halogen radicals. Halogen chemistry is thought to be especially active in pristine atmospheres that host 60 % to 70 % of the global burden of tropospheric ozone (Holmes et al., 2013). Based on model calculations, the halogen-driven ozone loss is estimated to be around 30 % in the marine boundary layer (von Glasow et al., 2002(von Glasow et al., , 2004Saiz-Lopez et al., 2014;Sherwen et al., 2016b) and 5 % to 20 % globally (Yang et al., 2005;Saiz-Lopez et al., 2012aSherwen et al., 2016b), going up to 25 % to 30 % in Badia et al. (2019). Consequently, it might be important to consider halogens in chemistryclimate modelling. It could also help to reconstruct more reliably past atmospheric compositions .
Despite its potential importance, only three global chemistry models have so far implemented detailed representations of halogenated (Cl, Br and I) tropospheric chemistry: CAM-Chem , GEOS-Chem (Sherwen et al., 2016b) and TOMCAT (Hossaini et al., 2015). Halogen chemistry is complex with a high computational cost (Fernandez et al., 2021), Earth system models used to assess how the evolution of atmospheric composition influences climate need to find the right compromise between the level of details (e.g. number of species and processes) here for halogen chemistry representations and the computing efficiency (box 6.1 in Szopa et al., 2021). Indeed, the computational cost has to be sufficiently low for multiple, long chemistry-climate simulations. So far, none of the chemistry-climate models involved in recent large intercomparison projects included tropospheric chemistry of very short-lived halogens and hence the impacts of their chemistry on key short-lived climate forcers (such as ozone and methane). Such multi-model projects are important. They are used to assess the budget of key chemical species, which can impact directly or indirectly climate or global air quality, and their sensitivity to changes in anthropogenic and natural emissions. The results are also used to evaluate the effects of international policies (e.g. the Gothenburg Protocol to Abate Acidification, Eutrophication and Ground-level Ozone) by quantifying source-receptor relations at continental scales (in the case of the Hemispheric Transport of Air Pollution project, HTAP; Fiore et al., 2009), to explore the sensitivity of short-lived climate forcers to changes in emissions of their precursors since the pre-industrial period, or to project their possible future evolutions under a range of scenario assumptions: in the case of the Atmospheric Chemistry and Climate Model Intercomparison Project (ACCMIP;  or, more recently, the Aerosol Chemistry Model Intercomparison Project (AerChemMIP; Stevenson et al., 2020).
The primary purpose of this work is to assess whether iodine, bromine and chlorine chemistry can, beyond their already established effect on global tropospheric ozone budget, affect the sensitivity of ozone to the changes in its precursors. In order to do so, we implemented a simplified representation of tropospheric halogen chemistry in a global chemistry-climate model, the LMDZ-INCA model (Laboratoire de Météorologie Dynamique (LMD) general circulation model, LMDZ, and Interactions with Chemistry and Aerosols, INCA) with a focus on reproducing the overall effect of halogen chemistry on global tropospheric ozone. Results of the three global tropospheric chemistry models that already include detailed representations of this chemistry are used to guide the choices of the chemical reactions and emissions in our simple halogen scheme and to assess the LMDZ-INCA simulations. First, we evaluate the ability of the model to simulate the global tropospheric budget of halogens and their effect on ozone by using the results of the CAM-Chem and GEOS-Chem models as references. The limitations of our model are also discussed, notably how some simplifications in our model affect the simulation of the spatial distribution of some specific halogenated species. Second, we investigate how the sensitivities of ozone, methane and hydroxyl radicals to perturbations in anthropogenic emissions of photo-oxidant precursors change when halogen chemistry is accounted for. For instance, we explore changes in ozone burden and net destruction in simulations considering present-day (PD) and pre-industrial (PI) emissions with and without halogen chemistry. Finally, the last section re-calls key findings, and we conclude with implications for chemistry-climate simulations of tropospheric ozone evolution.  Sadourny and Laval (1984). The primitive equations of the GCM are solved with a time step of 3 min, a large-scale transport of tracers is carried out every 15 min, and the physical processes are calculated at a time interval of 30 min. The LMDZ model provides at each time step the physical quantities (temperature, pressure, water vapour, turbulence, boundary layer, etc.) needed to compute the transport and chemistry in INCA.
INCA is an atmospheric chemistry model developed at LSCE (Laboratoire des Sciences du Climat et de l'Environnement). The model takes into account the primary emissions from natural sources or by anthropogenic activities. It computes physical processes such as dry and wet deposition and photochemical reactions (Hauglustaine et al., 2004). The atmospheric concentration fields are integrated in time by solving differential equations when called by LMDZ, with a time step of 30 min. For this development, we used a version of INCA (version v_5.2.1) simulating interactively the gaseous photochemistry, whereas aerosol distributions obtained with a more complex version of INCA are prescribed through climatologies. This version represents in particular the chemistry of the inorganic species O x , NO x and HO x , as well as CH 4 , CO, non-methane volatile organic compounds (NMVOCs) and their photo-oxidised products (Folberth et al., 2006).
Biogenic emissions like those of isoprene are provided by the global vegetation model ORCHIDEE (Organising Carbon and Hydrology In Dynamic Ecosystems) (Folberth et al., 2006). Anthropogenic emissions are based on the Community Emissions Data System (CEDS) v2016-07-2 emission inventory for the year 2010 (Hoesly et al., 2018).
The LMDZ-INCA model has been used and assessed in numerous previous model intercomparison exercises aimed at simulating the pre-industrial, present and future tropospheric composition (e.g. Szopa et al., 2013); the long range transport of air pollution in the HTAP exercise (Fiore et al., 2009); or the chemistry-climate interactions in ACCMIP .

Halogen chemistry
The chemistry of halogenated compounds implemented in LMDZ-INCA is essentially based on the scheme developed for CAM-Chem (Ordóñez et al., 2012;Prados-Roman et al., 2015a;Saiz-Lopez et al., 2012a and GEOS-Chem (Sherwen et al., 2016a), with a priority given to the reactions of importance for the ozone sinks ( Fig. 1).
A total of 42 halogenated species have been added to the INCA chemical scheme (Table 1). The list of the 160 chemical reactions added to the INCA chemical scheme can be found in Tables 2, 3, 4 and 5. Photolysis reactions of halogens as well as relative quantum yields and absorption crosssections are listed in Table 2. Uptake coefficients for heterogeneous reactions are based on experimental values and data selected for other models and are reported in Table 5. Recycling through heterogeneous chemistry on sea salt and sulfate aerosols is included, while the heterogeneous bromine reactions in water droplets and on ice crystals in clouds are not. Dry deposition is calculated with the Wesely scheme (Wesely, 1989), while washout of gases by precipitation is simulated by using Henry's law constants detailed in Table S1 in the Supplement. Some VSL halocarbons, CH 2 BrCl, CHBr 2 Cl and CHBrCl 2 , represent a few percent of the emissions of halocarbons (Ordóñez et al., 2012). Considering the high uncertainty in their emissions and their limited abundances, we chose not to consider their chemistry in this version.

Simulations
For these developments, LMDZ-INCA is used with a 3.75 • longitude × 1.9 • latitude resolution over 39 vertical layers (up to 70 km with 15 levels in the stratosphere). This resolution, coarser than the resolution now commonly used in the CMIP6 model exercise, was the standard LMDZ resolution for CMIP5 (Hourdin et al., 2013) and allows faster simulations necessary in model development. The LMDZ wind fields and sea surface temperatures are nudged on ECMWF reanalysis for the year 2010 in all simulations. The aerosol distribution and monthly variations used for heterogeneous reaction calculations are specified from climatologies precomputed with another version of the same model but with fully interactive calculations of aerosols. In this version, the stratospheric composition is prescribed. All simulations cover 1 year (2010) and follow a spin-up period of 1 year. Three types of simulations have been performed to simulate (i) present-day conditions, (ii) pre-industrial conditions and (iii) sensitivity to ozone precursors. The surface emissions of NO x , CO and NMVOCs and CH 4 concentrations considered in the present-day and pre-industrial simulations are presented in Table S2. For sensitivity tests, simulations with individual emission perturbations of photo-oxidant precursors (−20 %) of NO x , CO and NMVOCs and concentration reduction (−20 %) of CH 4 were performed. For each type of simulation, two versions of the model are applied: one accounting for the halogen chemistry (called Halo simulations hereafter) and a reference simulation without halogen chemistry called "NoHalo".

Present-day set-up
Primary halogenated compounds can originate from both inorganic and organic sources. Inorganic species are reactive species emitted by different sources (e.g. oceans, sea salts) but also result from atmospheric reactions (Simpson et al., 2015). Organic species are only of primary origin and are photo-oxidised in the atmosphere to form halogen radicals (Simpson et al., 2015). Table 6 shows the emissions of halogenated compounds as considered in INCA. Organobromines (e.g. CHBr 3 , CH 2 Br 2 ), organo-iodines (e.g. CH 3 I, CH 2 I 2 ) and inter-halogen species (e.g. CH 2 IBr, CH 2 ICl) are those used by the CAM-Chem model (Ordóñez et al., 2012).
Sources of organic chlorine are diverse, mainly originating from biological activity on the surface of the oceans and from forest fires. The emissions of these chlorocarbons are detailed in Table S2. In INCA, inorganic chlorine, HCl -the major reservoir of chlorine in the atmosphere -is considered to be emitted by sea salts (90 000 Gg Cl yr −1 ) and by forest fires (6400 Gg Cl yr −1 ) with total based on Lobert et al. (1999) and Hossaini et al. (2016). Note that more comprehensive but hard to implement approaches relying on pHdependent parameterisation of the dechlorination of sea salts have been recently implemented in GEOS-Chem. HCl emissions from industrial sources (McCulloch et al., 1999) are neglected as recommended by Wang et al. (2019) to avoid overestimating HCl in urban areas. ClNO 2 is emitted by the heterogeneous reaction of N 2 O 5 on sea salts and produces 2217 Gg Cl yr −1 .
CH 3 Br emissions are from oceanic origin (40 %) and from agriculture (60 %) (Barker et al., 2016) and represent 91 Gg of Br per year (Schmidt et al., 2016). IBr and ICl originate from the heterogeneous reaction of HOI, INO 2 and INO 3 on sea salt aerosols (Table 5) as considered in other models Sherwen et al., 2016a) but without accounting for pH dependence (McFiggans et al., 2000).
The inorganic iodine distribution of emissions (HOI and I 2 ) are prescribed from a climatology computed by the CAM-Chem model considering sea-air fluxes as a function of the oxidation of aqueous iodide by atmospheric ozone on the ocean surface parameterised by Carpenter et al. (2013) and MacDonald et al. (2014) and which depends on the ozone deposition to the ocean surface, the sea surface temperature and the wind speed (Prados-Roman et al., 2015b).

Pre-industrial set-up
The aim of this simulation is to explore how the sensitivity of ozone to anthropogenic emissions perturbations have changed since the pre-industrial period (a present-day climate is thus considered). Only emissions from natural sources are included for the simulation of the pre-industrial troposphere, following the approach used by other modelling teams to quantify this perturbation. Emissions of halogenated compounds from biomass burning are reduced to 10 % of their current values (Wang and Jacob, 1998). The methane concentration is set at 700 ppbv (IPCC, 2013). Emissions of biogenic volatile organic compounds are kept at their present-day level due to uncertainties in the effect of their drivers over the past century and the resulting high uncertainty regarding the quantification of change in their emissions (Szopa et al., 2021). Organo-iodine emissions are almost unchanged between the PI and PD. For organo-bromines, pre-industrial emissions of CHBr 3 and CH 2 Br 2 are equivalent to the PD ones, following the methodology of Parrella et al. (2012). The concentration of CH 3 Br is set to ∼ 5 ppbv and is considered to be in agreement with the measurements from ice cores (Saltzman et al., 2004). Pre-industrial emissions of CH 3 Cl, CHCl 3 and CH 2 Cl 2 are from , Hu (2012), Montzka et al. (2010), and Worton et al. (2006). HOI and I 2 emissions, pre-computed by CAM-Chem, are ∼ 40 % lower than today because they are sensitive to atmospheric ozone content which was lower in the pre-industrial era (Prados-Roman et al., 2015a;Sherwen et al., 2017). Global emissions for each primary halogenated compound are reported in Table 6.
3 Evaluation of the representation of the halogen-driven ozone sinks We first assess the ozone sensitivity to halogen chemistry in LMDZ-INCA, and then we evaluate the performance of the model in simulating the climatology of halogen species which are important for the global ozone loss. In order to do  The reactions come from GEOS-Chem (Parrella et al., 2012;Eastham et al., 2014;Schmidt et al., 2016;Sherwen et al., 2016a, b) and TOMCAT (Hossaini et al., 2016) and from Thomas et al. (2011). Lumped species are detailed in Folberth et al. (2006).  The reactions come from previous updates to halogen chemistry in GEOS-Chem (Parrella et al., 2012;Eastham et al., 2014;Schmidt et al., 2016;Sherwen et al., 2016a, b). The lower pressure limit rate (k 0 ) is given by A 0 ( 300 T ) x . The high pressure limit rate is given by k ∞ . Fc characterises the fall-off curve of the reaction as described by Atkinson et al. (2007). so, model outputs for compounds of interest in this chemistry are compared to a range of halogen observations and simulations from other models with halogen chemistry. An analysis of the model tropospheric budgets of key halogen species is also performed.   Sources of organic species (CH 3 X, CH 2 X 2 , CHX 3 ) are shown in terms of emissions. I 2 and HOI are the inorganic oceanic source that depends on surface ozone (Carpenter et al., 2013). IBr and ICl are released following gaseous iodine uptake on sea salt. ClNO 2 is emitted after absorption of N 2 O 5 on sea salt. Figure 2 shows the absolute changes in surface ozone and relative change in zonally averaged ozone in a simulation with halogen chemistry relative to a simulation without it. Significant relative losses are found in the Southern Hemisphere, which are due to the strong emissions of halogenated species from the ocean in this region (Long et al., 2014;Saiz-Lopez et al., 2012a;Schmidt et al., 2016;Sherwen et al., 2016a, b) and the smaller absolute abundance of ozone. Most of the change in total tropospheric ozone concentration occurs in the free troposphere between 350 and 900 hPa (63 %). This result is similar in magnitude to the relative change reports in GEOS-Chem (65 %; Sherwen et al., 2016a) and CAM-Chem (65 %; Iglesias-Suarez et al., 2020).

Impact on ozone loss in the present day
The annual production of ozone decreases by 4.7 % when the chemistry of halogens is considered. This decrease is due to a reduction in NO x concentrations of 32.8 % (Table 7) due to the hydrolysis of XNO 3 (X = Cl, Br, I) leading to deposition losses, in accordance with previous studies (Saiz-Lopez et al., 2012a;Long et al., 2014;von Glasow et al., 2004;Parrella et al., 2012;Schmidt et al., 2016;Sherwen et al., 2016a).
The vertical profile of the different contributions to chemical ozone loss is shown in Fig. 3. The processes leading to ozone loss differ considerably between simulations with and without halogenated compounds (Table 7). Halogens account for nearly ∼ 22 % of the ozone loss in the boundary layer in the halogen simulation, as shown in Fig. 3. Halogens represent 55 % of the Ox (Ox = odd oxygen) loss in the upper troposphere (350 hPa > pressure > tropopause), while this value is about 33 % and 40 % in GEOS-Chem with the iodine plus bromine and the iodine plus bromine plus chlorine chemistry, respectively. Such differences could be partly due to differences in convection schemes which induce different uplifts of IO in the upper troposphere resulting in different IO concentrations (see Sect. 3.1.3 below) but also to different representations of the stratospheric composition or chemistry between models.
When chemistry of halogenated compounds is considered, the total annual chemical ozone loss decreases by 2.4 % (Table 7). The relative contributions of different halogens to ozone loss are very similar to those of GEOS-Chem in Sherwen et al. (2016b). The loss of ozone by chlorine, bromine and iodine represents 0.5 %, 4.7 % and 18 % of the total losses, respectively; thus, halogens represent 23.2 % of the total loss of ozone. The sum of the ozone loss caused by the halogens is 984 Tg yr −1 . Iodine is responsible for 78 % of that halogen-mediated loss, while bromine accounts for 19 %, which is close to the values reported by Iglesias-Suarez et al. (2020) (80 %, 16 %, respectively). Iodine is responsible for an overall tropospheric ozone loss of 762 Tg yr −1 (comparable to the 748 Tg yr −1 reported by Sherwen et al., 2016b). Much of this loss arises from the photolysis of HOI after its production from the reaction of IO with HO 2 and the photolysis of OIO after its production from reactions implicating IO and I x O y . Bromine is responsible for an overall tropospheric loss of 187 Tg yr −1 which mostly arises from the photolysis of HOBr after its production from the reaction of BrO with HO 2 .
The inclusion of halogens leads to an overall decrease in the net production of tropospheric ozone (PO 3 -LO 3 ) of 20 % (126 Tg yr −1 ). Consequently, the lifetime of ozone (= tropospheric burden / losses) is significantly decreased by 3.2 d and reaches 17.9 d.

Present-day budgets for halogenated species in LMDZ-INCA
In order to reflect the role played by the halogenated species in their respective cycle, inorganic halogens are presented hereafter under the name X y which can be divided into HX  (the acidic forms of halogenated compounds) and X* (composed of the nocturnal reservoirs of reactive halogenated species and the reactive halogenated species). HX in LMDZ-INCA takes part in heterogeneous bimolecular reactions. In CAM-Chem and GEOS-Chem, the representation of this chemistry is more detailed: HX can be incorporated in droplets or particles and dissociate (Ordóñez et al., 2012;Prados-Roman et al., 2015a;Saiz-Lopez et al., 2012aSherwen et al., 2016a). X* is subdivided into the "reactive halogen" XO x ; grouping X + XO*; and the reservoir species X 2 , HOX, XNO 2 and XNO 3 , which have a relatively longer lifetime. These last species are formed mostly at night and serve as reservoirs of reactive halogens far from their primary sources.
In the following subsections, we present, for each halogen, the species considered in each reservoir, their respective tropospheric mass burden and the main fluxes between these reservoirs. The purpose is to evaluate the representation of halogenated species of interest for ozone and to discuss the limitations of the model in the computation of the budget of each halogen family.

Iodine
Iodine chemistry is the most impactful with regard to the ozone loss. The simplified iodine cycle is shown schematically in Fig. 4. The burden of inorganic iodine excluding IBr or ICl and iodine in aerosols is ∼ 22.0 Gg in the troposphere in LMDZ-INCA, comparable to the value of ∼ 27.9 Gg reported by Sherwen et al. (2016a, Fig . 7). HOI and OIO species, whose photolysis result in atomic iodine production, are the main intermediates for direct ozone loss. HOI is the most abundant iodine species, and its wet deposition controls the total atmospheric iodine burden. Its tropospheric burden of 12.8 Gg I in INCA is very close to the 13 Gg simulated by Sherwen et al. (2016b). The tropospheric burden of OIO in INCA, 0.43 Gg I, is comparable to the 0.55 Gg I also reported in Sherwen et al. (2016b). Overall, the global levels of HOI and OIO simulated by INCA are similar to their simulation by other models described in the literature. The tropospheric burdens of HI, INO 2 and I 2 sum up to 1 Gg I, lower than the burdens reported for GEOS-Chem and totalling 2.48 Gg I (Sherwen et al., 2016b). The global burden of I 2 in INCA is lower than that in GEOS-Chem (Sherwen et al., 2016b) because in Sherwen et al. (2016b) heterogeneous reactions transform INO 3 , INO 2 and HOI into I 2 , whereas in INCA, these reactions produce IBr and ICl, similarly to the heterogeneous reactions adopted by Sherwen et al. (2016a). Other quantifications in the literature are lacking to explain the differences in INO 2 and HI; however, Saiz-Lopez et al. (2014) indicate that minor iodine species (I min = HI + OIO + INO 2 + INO) represent less than 5 % of the mass in tropospheric inorganic iodine (2.3 % in our case). In INCA, the tropospheric I y mean is 1.3 pptv (parts per trillion volume), whereas previous model studies found a range of ∼ 0.5 to 1 pptv Sherwen et al., 2016a, b). A large set of observations has been assembled by Sherwen et al. (2016b) for IO. Measurements of IO during oceanic campaigns show a strong diurnal variation (IO being present only during daytime) with climatological concentrations (local observations averaged over 20 • bins) comprised between 0.3 and 1 pptv. These IO concentrations show a weak latitudinal variability and a rapid decrease with altitude to reach values of about 0.2 pptv from 2 km to the tropopause (observations gathered and reported by Sherwen et al., 2016b). Surface values simulated by INCA agree within a factor of 2 with these observations. Over oceans, the annual mean surface concentrations of IO is approximately 0.2 pptv with concentrations of between 0 and 2 pptv. This is similar to the range simulated by Sherwen et al. (2016b) (0.25-1 pptv) and to the observations reported by Prados-Roman et al. (2015b). In the tropics (20 • S-20 • N), levels of surface IO in INCA (0.35 pptv) are in agreement with model results reported for CAM-Chem (∼ 0.4 pptv) . In the free troposphere, IO levels slightly increase with altitude to reach diurnal values of ∼ 0.5 pptv at 10 km in INCA. This overestimation of IO in the upper troposphere is common to all current model-simulated iodine and remains unclear (Sherwen et al., 2016b;Badia et al., 2019).
The IO x family influences the oxidative capacity of the troposphere through the catalytic cycles of ozone depletion (Brasseur, and Solomon, 2005), modifies the HO x (HO 2 /OH) and NO x (NO 2 /NO) ratios (Bloss et al., 2010), produces photolabile higher-order iodine oxides (I x O y ) Lewis et al., 2020), and generates different forms of inorganic iodine (Saiz-Lopez et al., 2012b). For that reason, a correct representation of the branching ratios and the generation and loss of IO x is crucial. Globally, IO x production in INCA is dominated by inorganic photolysis of HOI (80.5 %), OIO (10.2 %), INO 3 (5.4 %) and I 2 O x (3.9 %). This partition is comparable to the one reported in Sherwen et al. (2016a), who found a 76 % HOI photolysis contribution and a 11 % OIO photolysis contribution. The main IO x loss pathway is the production of HOI by IO reaction with HO 2 (78.5 %), with additional loss pathways by self-reaction and reactions with NO x and BrO contributing 11.3 %, 5.7 % and 4.5 %, respectively -numbers that are very close to the branching ratios reported in Sherwen et al. (2016a) (77 %, 10 %, 7.7 % and 4.6 %, respectively).

Bromine
The tropospheric bromine budget and simplified recycling cycle are illustrated in Fig. 5. The tropospheric inorganic bromine burden is ∼ 27 Gg in LMDZ-INCA, close to the 28 Gg reported for the most recent version of GEOS-Chem (Wang et al., 2021). On an annual average, INCA simulates daytime BrO values in the free troposphere of ∼ 0.3 pptv within the tropics, close to the results in CAM-Chem (∼ 0.2 pptv) . At the tropical tropopause layer, the modelled tropical BrO concentrations in INCA are ∼ 3 pptv, close to the observed concentrations of 2.0 ± 1.5 pptv and of 3.2 ± 1.6 pptv at two different altitudes performed by Dorf et al. (2008) and reported by Fernandez et al. (2014) for the evaluation of CAM-Chem. Overall, the global tropospheric mean of BrO simulated by INCA is 0.6 pptv (in agreement with model results from Parrella et al., 2012) but is in the lower end of the 0.5 to 2 pptv range inferred from observations (Dorf et al., 2008;Saiz-Lopez and von Glasow, 2012;Stone et al., 2012).
In this first implementation of halogen chemistry in a version of LMDZ-INCA, the focus is on gaseous chemistry. Aqueous chemistry in droplets or on ice crystals in clouds is not accounted for in the first version of the scheme because it is complex and uncertain and, more importantly, might not be a dominant mechanism for global ozone loss. Nonetheless, the scheme includes recycling through heterogeneous chemistry on sea salt and sulfate aerosols. A study by Badia et al. (2019) using WRF-Chem shows that heterogeneous halogen chemistry (on all types of aerosols) allows a better vertical Br y partitioning but contributes little (∼ 6 %) to ozone loss compared to the gas-phase halogen chemistry (25 %) on a global scale. Our choice to include, in this first step, reactions on sea salt and sulfate aerosols only is suited to reproduce bromine-mediated ozone loss in the marine boundary layer but probably tends to slightly underestimate the halogen-mediated ozone loss in the global troposphere.
Due to this lack of in-cloud halogen chemistry in cloud droplets and ice crystals, whose importance increases with altitude, the bromine chemistry in INCA is less dynamic and the recycling of bromine is less efficient in LMDZ-INCA compared to those reported for some versions of GEOS-Chem (Parrella et al., 2012;Schmidt et al., 2016;Sherwen et al., 2016a;Zhu et al., 2019). Note, however, that the recently updated HOBr heterogeneous chemistry implemented in the most recent GEOS-Chem version (Wang et al., 2021) results in much slower recycling of HOBr in clouds and aqueous aerosols than previously calculated. This results in fluxes from heterogeneous reactions (specifically HOBr + HBr → Br 2 + H 2 O and HOBr + HCl → BrCl + H 2 O) that sum up to 32 500 Gg Br yr −1 , which is smaller for example than our estimates of 86 000 Gg Br yr −1 but remains of the same order of magnitude. Finally, our branching ratios are found to be consistent with those in Schmidt et al. (2016) and Wang et al. (2021). This leads us to conclude that, although LMDZ-INCA does not integrate the same degree of details in tropospheric halogen heterogeneous chemistry compared to its peers, global fluxes in the bromine cycle of INCA appear to be more comparable in terms of the bromine cycle with those in the newest version of GEOS-Chem (see the detailed bromine budget and cycle in Fig. S1) and CAM-Chem (Iglesias-Suarez et al., 2020) than fluxes in the older versions of GEOS-Chem (Parrella et al., 2012;Schmidt et al., 2016;Sherwen et al., 2016a;Zhu et al., 2019). All these estimates remain of the same order of magnitude, but the range of fluxes illustrates the level of uncertainties about the strength of the halogen heterogeneous recycling.

Chlorine
The total tropospheric burden of Cl y (inorganic chlorine) is 290 Gg of Cl, largely dominated by HCl. The total load of HCl is smaller than the one reported for GEOS-Chem  essentially because of the difference in the representation of the emissions. However, the total burden of reservoir species has a similar magnitude (10.4 Gg Cl). The halogen chemistry in INCA leads to a 25 % higher quantity of available reactive chlorine compounds: Cl and ClO* (ClO* being ClO + ClOO + OClO + 2Cl 2 O 2 ).
Atomic chlorine has two main chemical sinks. The first and most important in terms of reaction flux is the transformation of atomic Cl into ClO*. This conversion takes place mainly by reaction with ozone (99.5 %), while the conversion of ClO x to Cl is dominated by reaction with NO (93.9 %). The second path is represented by the chemical loss of reactive chlorine Cl to inactive chlorine HCl (Fig. 6), which is due to the reactions of Cl with methane and several NMVOCs. The conversion of Cl to ClO* is in strong competition with the second route. By dividing the fluxes of the conversion Cl to ClO* by the flux of conversion of Cl to HCl, we find that the chain length is short and equal to 1.3 (4.8×10 4 /3.6×10 4 ). This is close to the chain length of 1.6 calculated by Wang et al. (2019) and shows that the branching ratios are of the same order of magnitude in both models despite differences in NMVOC and CH 4 concentrations between the two models. Overall, we find that the partition between Cl and ClO* is correct in the model as well as the balance between radical and non-radical species, although the latter is a little more shifted towards reactive species in INCA compared to GEOS-Chem and CAM-Chem.
The average of climatological observations (mean over several months and representative of large areas) over six oceanic and coastal remote locations shows a mean surface HCl mixing ratio of 196±114 pptv, whereas INCA simulates an average of 193 ± 98 pptv, compared to the 246 ± 90 pptv in GEOS-Chem ) (see Fig. S2). However, for continental regions, LMDZ-INCA underestimates the HCl concentrations compared to observations mainly because it does not include industrial emissions. ClNO 2 over continental areas does not exceed a few tens of pptvs in INCA, whereas observations at several continental locations, as gathered by Wang et al. (2019) (their Table 5), show concentrations of a few hundreds to thousands of pptvs over such locations for hourly maxima. This difference can arise from the coarse resolution of the model which does not resolve the spatial heterogeneity of these local observations (Sherwen et al., 2016a) and the absence of ClNO 2 pro-duction in sulfate-nitrate-ammonium aerosols . The chlorine Cl* measurements to which we compare LMDZ-INCA results represent the sum of Cl 2 , HOCl, ClNO 2 , ClNO 3 and other minor components as explained in Pszenny et al. (1993). LMDZ-INCA Cl* values are of the same order of magnitude as measurements (see Table S3).

Perturbation of the ozone sensitivity to its precursors by halogen chemistry
In this section, the evolution of the model-calculated ozone sensitivity from the pre-industrial period to the present day is assessed, with and without halogen chemistry, and compared with results from other models. Then the impact of halogen chemistry on the sensitivity of ozone, OH radical and CH 4 concentrations to the emissions of anthropogenic ozone precursors is discussed.

Pre-industrial to present-day changes
With and without halogen chemistry, the model simulates significantly lower ozone concentrations in the pre-industrial compared to present due to the decrease in NO x concentrations. Taking into account halogen chemistry reduces the average surface ozone concentrations for pre-industrial condi- Globally, halogens reduce the tropospheric ozone burden by 67.1 Tg in the present day and 48.2 Tg in pre-industrial conditions (Table 7). Between PI and PD, ozone increases by 98.1 Tg (113 Tg in Sherwen et al., 2017) without halogen chemistry and by 79.2 Tg (90 Tg in Sherwen et al., 2017) with halogen chemistry. Both in PD and PI, halogens are responsible for about 23 %-28 % of the ozone destruction. Tropospheric ozone lifetimes drop from 21.1 to 17.9 d in the present with the inclusion of halogens and from 25.6 to 21.2 d in the pre-industrial period. The sensitivity of the burden of ground-level ozone to pre-industrial-present changes is ∼ 20 % lower when the chemistry of the halogens is considered. The changes observed in LMDZ-INCA and schematised in Fig. 7 are in agreement with the changes in Sherwen et al. (2017).

Sensitivity runs
In order to study the effect of halogen chemistry on the sensitivity of the tropospheric oxidative system, we explore, in Fig. 8, the relative changes (%) in ozone and OH when present-day emissions of photo-oxidant precursors (NO x , CO and NMVOCs and of CH 4 concentrations) are perturbed by a 20 % reduction, as done in the HTAP protocol. Our results show that the importance of sensitivities of ozone and OH to individual photo-oxidant precursor emissions increases as follows in terms of importance: NMVOCs, CO, NO x and CH 4 (Fig. 8). It is also clear that the OH radical is more sensitive to perturbations in CH 4 , NO x and CO than ozone. Notably, when accounting for halogens, ozone becomes more sensitive to changes in CH 4 and NO x and to a lesser extent to changes in NMVOCs. On the other hand, in the presence of halogen chemistry, the OH radical becomes more resilient in the presence of halogens to CH 4 , NO x and NMVOCs perturbations. Usually, solely the CH 4 consumption by OH is considered when the chemical lifetime of CH 4 in the troposphere is reported in the literature since it represents the main CH 4 sink. It thus provides a useful proxy for global tropospheric oxidising capacity (Wild et al., 2020). The simulations examined here show that the CH 4 lifetime with respect to OH remains a useful proxy for global tropospheric oxidising capacity for all perturbations. However, in the presence of halogen chemistry, when NO x emissions are reduced, atomic chlorine, a secondary oxidant of methane, becomes more efficient and increases the sensitivity of CH 4 to NO x emission change. This exercise shows for the first  time the influence of halogens on the sensitivity of the tropospheric oxidative system to changes in the emissions of photo-oxidant precursor (NO x , CO and NMVOCs and CH 4 ).

Impact on methane lifetime
The inclusion of tropospheric halogens affects the atmospheric concentrations of a large number of compounds in the model. A marked decrease in oxidants (O 3 , OH, HO 2 , H 2 O 2 ) is then computed. This leads to an increase in the con-centrations of CO (+8.5 %) and some NMVOCs, including isoprene (+7.1 %). The lifetime of methane with respect to chlorine and the OH radical, as well as the share of chlorine in the losses of methane in the literature, is shown in Table 8. In the reference simulation, the tropospheric chemical lifetime of methane due to its consumption by OH is 7.9 years. With the inclusion of halogen chemistry, the OH concentration decreases, increasing this chemical lifetime of methane to 9.3 years. The chemical lifetime of methane linked to oxidation by Cl is 333 years, less than the 384 years suggested by Hossaini et al. (2016) but closer to IPCC Fifth Assessment Report (AR5) estimates. This oxidation by Cl represents 2.7 % of the total loss (close to the values reported in the literature and in Table 8). Thus, in LMDZ-INCA, accounting for both OH and the chlorinated radical consumption, methane oxidation leads to a total chemical lifetime of CH 4 of 9.0 years.

Conclusions
Global chemistry climate models allow us to simulate the evolution of the atmospheric composition and are used for future projections in the framework of international model intercomparisons such as AerChemMIP. Most of these models do not account for halogen chemistry. We have implemented a simplified representation of the halogenated chemistry in the three-dimensional climate-chemistry model LMDZ-  INCA to account for its effects on the tropospheric chemistry, notably ozone. First, the halogenated species, their sources and sinks, their reactions in the gas and heterogeneous phases, and their physicochemical characteristics have been integrated into INCA. In order to simplify the integration of this chemistry, a chemical scheme already tested and evaluated in the CAM-Chem model has been used as a basis for our developments.
The representation of this chemistry has been shown to be sufficient in a first approach to simulate properly the impact of halogens on the photo-oxidising system in the troposphere and particularly in the atmospheric boundary layer. The response of the tropospheric chemical system O x , HO x , NO x , CH 4 and NMVOCs to the incorporation of halogen chemistry has been quantified. On a global tropospheric scale, this chemistry induces decreases in the ozone burden by 22 %, in OH by 8 % and in NO x by 33 %. An increase in CO and NMVOC concentrations is simulated, and the lifetime of methane increases by 1 year. These changes are consistent with those reported in the literature. Comparisons between simulations considering pre-industrial and present-day emissions show that the sensitivity of the ozone burden to changes in emissions over the pre-industrial era is ∼ 20 % lower in the presence of halogen chemistry (mainly because of the chemistry of iodine).
Sensitivity tests, consisting of reducing in turn the emissions of ozone precursors and methane concentrations, show for the first time that ozone becomes more sensitive to NO x , CO and NMVOC perturbations in the presence of halogen chemistry.
Although still rarely represented in chemistry climate models, our model results confirm that the chemistry of tropospheric halogens (Cl, Br, I) plays a significant role in the global loss of tropospheric ozone. Since halogen chemistry affects substantially the ozone burden both in pre-industrial and present-day chemical conditions and the ozone sensitivity to changes in anthropogenic emissions, international exercises examining the evolution of ozone between preindustrial, present and future periods should consider taking into account halogen chemistry.
Code and data availability. The main model codes and the chemical scheme with halogenated compounds as well as the data outputs are publicly accessible through https://doi.org/10.14768/56b1b070deae-4357-a2dd-d75d27659f63 (Szopa et al., 2023).
Author contributions. The implementation of the halogens chemical scheme has been done by CC, SS and AC. The design of the numerical simulations was conducted by CC, SS, ASL and CAC. CC and SS prepared the figures. CC, SS, SB, ASL and CAC analysed the results and drafted the paper.
Competing interests. At least one of the (co-)authors is a member of the editorial board of Geoscientific Model Development. The peer-review process was guided by an independent editor, and the authors also have no other competing interests to declare.
Disclaimer. Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.