Iodine chemistry in the chemistry–climate model SOCOL-AERv2-I

Abstract

In this paper, we present a new version of the chemistry–climate model SOCOL-AERv2 supplemented by an iodine chemistry module. We perform three 20-year ensemble experiments to assess the validity of the modeled iodine and to quantify the effects of iodine on ozone. The iodine distributions obtained with SOCOL-AERv2-I agree well with AMAX-DOAS observations and with CAM-chem model simulations. For the present-day atmosphere, the model suggests that the iodine-induced chemistry leads to a 3 %–4 % reduction in the ozone column, which is greatest at high latitudes. The model indicates the strongest influence of iodine in the lower stratosphere with 30 ppbv less ozone at low latitudes and up to 100 ppbv less at high latitudes. In the troposphere, the account of the iodine chemistry reduces the tropospheric ozone concentration by 5 %–10 % depending on geographical location. In the lower troposphere, 75 % of the modeled ozone reduction originates from inorganic sources of iodine, 25 % from organic sources of iodine. At 50 hPa, the results show that the impacts of iodine from both sources are comparable. Finally, we determine the sensitivity of ozone to iodine by applying a 2-fold increase in iodine emissions, as it might be representative for iodine by the end of this century. This reduces the ozone column globally by an additional 1.5 %–2.5 %. Our results demonstrate the sensitivity of atmospheric ozone to iodine chemistry for present and future conditions, but uncertainties remain high due to the paucity of observational data of iodine species.

1 Introduction

Emissions of chlorine- and bromine-containing halogen compounds have long been the subject of scientific investigation because they play an important role in the catalytic destruction cycles of stratospheric ozone . Recent studies demonstrate the success of the Montreal Protocol and its amendments in phasing out the emissions of chlorine- and bromine-containing substances and also point to early signs of stratospheric ozone recovery e.g.,. At the same time, our understanding of iodine-induced ozone depletion is incomplete because of the extremely low stratospheric iodine concentrations, typically in the parts-per-trillion (pptv) range . However, the ozone depletion efficiency of iodine-containing species on a per-molecule basis is hundreds of times higher than that of chlorine , because iodine reservoir species readily convert to free radicals, I and IO, in the presence of sunlight. Therefore, iodine chemistry remains suspected of affecting stratospheric ozone   and slowing the recovery of the ozone layer . Still, iodine-induced ozone depletion is not regulated under the Montreal Protocol as the mixing ratio of total inorganic iodine in the atmosphere is extremely low ($\sim \mathrm{1}$ pptv) and because there is almost no direct anthropogenic production of iodine-containing species .

The study of iodine source gases (ISGs) began in the late 1970s. The main organic source of atmospheric iodine is found to be iodomethane (CH${}_{\mathrm{3}}$I) as suggested by and . Photochemical production in ocean surface water involving phytoplankton and algae is the main source of organic and inorganic iodine that enters the atmosphere ( and references therein). Iodocarbons are further emitted by wetland plants , biomass burning , rice paddies , and volcanoes . ISGs are produced mostly in the tropics as the sea surface temperature (SST) is an important driver of biological activity. The 5–6 d lifetime of CH${}_{\mathrm{3}}$I is long enough to allow a small fraction of it to reach the upper troposphere via tropical convection with an estimated mixing ratio $\sim \mathrm{0.1}$ pptv in the tropical tropopause layer . Other organic ISGs, such as CH${}_{\mathrm{2}}$I${}_{\mathrm{2}}$, CH${}_{\mathrm{2}}$ICl, CH${}_{\mathrm{2}}$IBr, and C${}_{\mathrm{2}}$H${}_{\mathrm{5}}$I, are only minor contributors to iodine abundance in the upper troposphere and stratosphere since their lifetimes after volatilization are just a few hours on average.

Knowledge of the sources of atmospheric iodine has changed when the importance of inorganic iodine emissions was recognized, such as molecular iodine, I${}_{\mathrm{2}}$ , and hypoiodous acid, HOI . These are formed through reactions of near-surface ozone with oceanic iodide (I${}^{-}$), i.e., O${}_{\mathrm{3}}$ $+$ H${}^{+}$ $+$ I${}^{-}$ $\to$ HOI $+$ O${}_{\mathrm{2}}$ and HOI $+$ H${}^{+}$ $+$ I${}^{-}$ $\to$ I${}_{\mathrm{2}}$ $+$ H${}_{\mathrm{2}}$O, which are possibly responsible for up to 75 % of total iodine emissions . In conjunction with these findings, used a numerical parametrization of sea-surface iodide in numerical modeling of HOI and I${}_{\mathrm{2}}$ fluxes. Caused by the growing anthropogenic air pollution and the concomitant increase in surface ozone levels, iodine emissions have risen continuously during the past decades. and argued that the atmospheric iodine loading must have increased by at least a factor of 3 since the 1950s, because of increasing anthropogenic NO${}_x$ emissions resulting in an increase in near-surface ozone. This is in addition to a simultaneous increase in SST due to global warming with correspondingly enhanced metabolic rates of oceanic biota. While the future surface ozone evolution has a large spread in model projections covering the 21st century , which results in large uncertainty in future iodine emissions, the continuous increase in surface temperatures is predicted to raise tropospheric iodine levels throughout the 21st century based on Representative Concentration Pathway (RCP) scenarios . In addition, the effectiveness of iodine for ozone depletion (O${}_{\mathrm{3}}$ destruction by one I atom relative to O${}_{\mathrm{3}}$ destruction by one Cl atom) is predicted to change, with the share of halogen-induced ozone loss due to reactions of iodine likely to grow in the future stratosphere .

Using a two-dimensional model of atmospheric chemistry and dynamics, demonstrated a high impact of iodine chemistry on ozone even when only 1 pptv of total stratospheric iodine is present. These model results were later put in perspective by , who presented an upper limit of concentration of reactive IO of less than 0.1–0.2 pptv in the upper troposphere–lower stratosphere measured during nine balloon flights over Scandinavia and France. This suggested only a small influence of iodine chemistry on stratospheric ozone. In contrast, using differential optical absorption spectroscopy (DOAS) found significant concentrations of IO between 0.65 and 0.8 pptv ($\pm \mathrm{0.2}$ pptv) over Spitsbergen Island. Recent aircraft campaigns, also applying DOAS over the West Pacific during CONTRAST (Convective Transport of Active Species in the Tropics) in January–February 2014 and over the East Pacific during TORERO (Tropical Ocean Troposphere Exchange of Reactive Halogen Species and Oxygenated VOC) in January–February 2012

The GV AMAX-DOAS Data can be found here: https://data.eol.ucar.edu/ (last access: 20 October 2021).

, measured a second maximum in IO mixing ratio of more than 0.1 pptv in the tropical tropopause layer that had not been observed before. Yet, globally representative quantitative measurements of iodine species are still not available. Quasi-global IO has been inferred from the SCIAMACHY satellite measurements , though with a very high level of uncertainty.

Modeling of atmospheric iodine chemistry started with proposing a detailed iodine chemistry scheme with gas and aqueous phase reactions in a box model, allowing iodine species to recycle via the uptake of HOI onto sea-salt particles resulting in the production of ICl and IBr (i.e., HOI $+$ Cl${}^{-}$ $+$ H${}^{+}$ $\to$ ICl $+$ H${}_{\mathrm{2}}$O, HOI $+$ Br${}^{-}$ $+$ H${}^{+}$ $\to$ IBr $+$ H${}_{\mathrm{2}}$O). Note that these recycling reactions constitute a net source of bromine and chlorine into the atmosphere but represent only a change in partitioning for the iodine species, which avoids iodine uptake and washout (e.g., as H${}^{+}$ $+$ I${}^{-}$). hypothesized that the measured increase in IO mixing ratio in the Pacific free troposphere might hint at heterogeneous recycling processes from aerosols back to the gas phase in the upper troposphere. were the first to use a global chemistry–climate model, namely CAM-chem, to incorporate a comprehensive bromine and iodine chemistry scheme. The global distribution of organic iodocarbons from CAM-chem showed reasonable agreement with observations in the marine boundary layer when including global oceanic ISGs of approximately 1.9 Tg (I)/year . Later, an updated version of the GEOS-Chem model demonstrated 3.83 Tg (I)/year of total ISGs emissions, when using the sea-surface iodide field from . used the sea-surface values from , which possibly underestimates the sea-surface iodide emissions in comparison to observations , suggesting only 2.7 Tg (I)/year of total emissions. In order to reproduce the IO observations of , introduced the heterogeneous recycling of iodine reservoirs (HOI and IONO${}_{\mathrm{2}}$) on ice crystals and implemented it in CAM-chem, assuming it to be analogy to bromine recycling on ice crystals in the upper troposphere . Accordingly, the total iodine injected into the stratosphere is currently thought to be up to 0.8 pptv , i.e., about 5 times larger than that from the previously published assessment . Also, CAM-chem results suggest that stratospheric iodine might be responsible for up to 30 % of halogen-mediated ozone depletion in the lower tropical stratosphere. At the same time, the results of obtained using the 3-D chemistry transport model TOMCAT suggest that iodine lowers stratospheric ozone only by $\sim \mathrm{3}$ % and total column ozone by 0.5 %. However, this model did not take the recycling mechanism described in into consideration, which increases the stratospheric iodine levels significantly beyond what was assumed by . It remains difficult to judge which results are closer to reality because of the paucity of observations.

This paper introduces the new version of the chemistry–climate model (CCM) SOCOL-AERv2 (Solar Climate Ozone Links coupled to a size-resolving sulfate aerosol module), which has been extended to include an iodine chemistry scheme. The main objective of this study is to verify the model and further constrain the impact of iodine chemistry on stratospheric ozone depletion, first based on present-day emissions, and second when applying a 2-fold increase in iodine emissions. Section 2 provides a short description of SOCOL-AERv2 (Sect. 2.1), improvements which have been made to obtain the new model version subsequently referred to as SOCOL-AERv2-I (Sect. 2.2), and the numerical experiments conducted with the new model version (Sect. 2.3). Section 3 presents the model results, i.e., the comparison of simulated and observed iodine: first, different aspects of iodine effect on present-day stratospheric ozone climatology are considered (Sect. 3.1), and, second, the sensitivity of ozone to iodine is presented (Sect. 3.2). Finally, the discussion and summary of the present study are provided in Sect. 4.

2 Model description and conducted experiments

2.1 The SOCOL-AERv2 chemistry–climate model

The chemistry–climate model SOCOL-AERv2 is the CCM SOCOLv3 coupled to a size-resolving sulfate aerosol module (AER) along with other important modifications for chemistry and deposition. The AER module of SOCOL was established by (CCM SOCOL-AERv1). The CCM SOCOL-AERv1 was substantially updated by (CCM SOCOL-AERv2) with an interactive deposition scheme, expanded tropospheric chemistry scheme, and improved sulfate mass and particle number conservation (less susceptible to numerical diffusion). The SOCOL-AERv2 consists of a dynamical core that is the middle atmosphere version of the spectral transform general circulation model MA-ECHAM5.4 (the Middle Atmosphere version of the European Centre/Hamburg Model version 5.4) , which has been interactively coupled to the MEZON atmospheric chemistry-transport module (Model for the Evaluation of OZONe Trends) . The coupling takes account of radiative forcing caused by ozone, H${}_{\mathrm{2}}$O, N${}_{\mathrm{2}}$O, CH${}_{\mathrm{4}}$, and CFCs. The MA-ECHAM developed at the Max Planck Institute for Meteorology (Hamburg, Germany) is based on primitive prognostic equations for meteorological parameters such as logarithm of surface pressure, temperature, humidity, vorticity, etc. The advection in MA-ECHAM is regulated by a flux-transform semi-Lagrangian scheme based on mass conservation and shape retention  . The standard SOCOL-AERv2 utilizes the Gaussian transform horizontal grid with T42 triangular truncation (64 latitudes $\times$ 128 longitudes) splitting the model space into grid cells of $\sim$2.5 $\times$ 2.5${}^{\circ }$ each. The vertical direction model grid consists of 39 levels in the hybrid sigma–pressure coordinate system covering the altitudes ranging from the ground surface to about 80 km (0.01 hPa). The model time step is 15 min for dynamical and physical processes whereas it is 2 h for atmospheric full radiation and chemical calculations. The CCM SOCOL-AERv2 uses the prescribed monthly fields of the sea surface temperature (SST) and ice coverage acquired from the Hadley Centre dataset . MEZON shares the horizontal and vertical spatial resolutions with MA-ECHAM5 and includes 95 chemical species and 215 gas-phase, 16 heterogeneous, and 75 photolysis reactions. For more details on CCM SOCOL-AERv2, see .

2.2 SOCOL-AERv2 updates to SOCOL-AERv2-I

Here, we introduce the SOCOL-AERv2-I that is the SOCOL-AERv2 extended with the iodine chemistry module. This module includes 61 gas-phase and 4 heterogeneous chemical reactions involving iodine, boundary conditions for prescribed iodocarbon emissions, interactive inorganic iodine emissions, wet/dry depositions of iodine species augmented with the deposition on sea-salt and sulfate aerosol particles, and effective uptake (removal)/reactive uptake (recycling) on tropospheric cloud ice. The following section describes all components of the iodine chemistry module in SOCOL-AERv2-I. The prescribed organic and interactive inorganic iodine source gases (ISGs) are presented in Table .

Table 1

Iodine source gases (ISGs) incorporated in SOCOL-AERv2-I.

ISG	Type	Resolution	Reference
I${}_{\mathrm{2}}$	interactive	model time step
HOI	interactive	model time step
CH${}_{\mathrm{3}}$I	prescribed	monthly
CH${}_{\mathrm{2}}$I${}_{\mathrm{2}}$	prescribed	monthly
CH${}_{\mathrm{2}}$ICl	prescribed	monthly
CH${}_{\mathrm{2}}$IBr	prescribed	monthly

Following , the organic iodocarbons have been obtained from the inventory of for CH${}_{\mathrm{3}}$I and from the 1-D model estimates of for CH${}_{\mathrm{2}}$I${}_{\mathrm{2}}$, CH${}_{\mathrm{2}}$ICl, and CH${}_{\mathrm{2}}$IBr. Organic ISGs have been parameterized in by a biogenic chlorophyll-$a$ (chl-$a$) dependent source in the tropical oceans . In our scheme, organic emissions are prescribed on a monthly basis. To simulate the long-term period, we repeat organic emissions fluxes at the beginning of every model year, so any interannual variability is not included.

Iodocarbon source fluxes were directly extracted from the GEOS-Chem model (v10-01) at a resolution of a 2 $\times$ 2.5${}^{\circ }$, including updates to both iodine and bromine chemistry from and , and interpolated on the SOCOL-AERv2-I horizontal grid (T42). The inorganic HOI $/$ I${}_{\mathrm{2}}$ fluxes are interactively calculated in the model using the numerical parametrization of and global sea surface iodide concentration calculated following . To derive HOI $/$ I${}_{\mathrm{2}}$ fluxes this parametrization utilizes the model fields of near-surface ozone (at the model level closest to the surface), surface wind speed, and SST.

The parametrization for the emission fluxes of $F_{\mathrm{HOI}}$ and $F_{I_{\mathrm{2}}}$ (in nmol m${}^{-\mathrm{2}}$ d${}^{-\mathrm{1}}$) as a function of the surface ozone mixing ratio ${\mathit{\chi }}_{{\mathrm{O}}_{\mathrm{3},\mathrm{surf}}}$ (in ppbv) are given by $$\begin{array}{}\text{1}& {\displaystyle }\begin{array}{rl}{F}_{\mathrm{HOI}}& ={\mathit{\chi }}_{{\mathrm{O}}_{\mathrm{3},\mathrm{surf}}}\times (\mathrm{4.15}\times {\mathrm{10}}^{\mathrm{5}}\times (M_{{\mathrm{I}}_{\mathrm{aq}}^{-}}^{\mathrm{1}/\mathrm{2}}/W_{\mathrm{s}})\\ & -(\mathrm{20.6}/W_{\mathrm{s}})-\mathrm{23}\mathrm{600}\times M_{{\mathrm{I}}_{\mathrm{aq}}^{-}}^{\mathrm{1}/\mathrm{2}}),\end{array}\text{2}& {\displaystyle }{F}_{{\mathrm{I}}_{\mathrm{2}}}={\mathit{\chi }}_{{\mathrm{O}}_{\mathrm{3},\mathrm{surf}}}\times M_{{\mathrm{I}}_{\mathrm{aq}}^{-}}^{\mathrm{1.3}}\times (\mathrm{1.74}\times {\mathrm{10}}^{\mathrm{9}}-\mathrm{6.54}\times {\mathrm{10}}^{\mathrm{8}}\times {\mathrm{lnW}}_{\mathrm{s}}),\end{array}$$ where $W_{\mathrm{s}}$ is the surface wind speed (in m s${}^{\mathrm{-}\mathrm{1}}$) and $M_{{\mathrm{I}}_{\mathrm{aq}}^{-}}$ is the sea surface iodide concentration expressed as molarity (mol dm${}^{\mathrm{-}\mathrm{3}}$)

3$$M_{I_{\mathrm{aq}}^{-}}=\mathrm{1.46}\times {\mathrm{10}}^{\mathrm{6}}\exp (-\mathrm{9134}/T_{\mathrm{ss}})$$ with sea surface temperature $T_{\mathrm{ss}}$ in K.

It is important to note that this parameterization is intended for calculating I${}_{\mathrm{2}}$ and HOI fluxes at wind speeds below 14 m s${}^{\mathrm{-}\mathrm{1}}$ (higher wind speeds lead to negative values of fluxes) since it is based on the approximation of measurements and there are no measurements when a storm occurs. Also, it should be mentioned that turbulent mixing of the interfacial layer with bulk seawater reduces the proportion of I${}_{\mathrm{2}}$ and HOI escaping into the atmosphere; therefore fluxes decrease with the wind speed . The obtained inorganic ISGs fluxes in the SOCOL-AERv2-I exceed GEOS-Chem and CAM-Chem emissions because SOCOL-AERv2-I overestimates near-surface O${}_{\mathrm{3}}$ . When we compared the ground-level ozone fields from SOCOL and GEOS-chem, we found that for most regions where HOI $/$ I${}_{\mathrm{2}}$ emissions are produced, the overestimation of ground-level ozone in SOCOL is $\sim$100 %. To ameliorate this bias, the surface ozone used within the parameterization was scaled by 0.5. This halved the emission, making HOI $/$ I${}_{\mathrm{2}}$ fluxes comparable: $\sim$2.4 Tg (I)/year in SOCOL-AERv2-I and $\sim$2.2 Tg (I)/year in GEOS-chem . Since there is an uncertainty in inorganic iodine emissions , such a difference in emissions between models is admissible. Iodocarbon source fluxes in both SOCOL-AERv2-I and GEOS-chem are identical and correspond to 0.5 Tg (I)/year. Iodine emission fluxes used in CCM SOCOL-AERv2-I are demonstrated in Fig. .

Figure 1

Annual-mean emission fluxes of organic (CH${}_{\mathrm{3}}$I, a; CH${}_{\mathrm{2}}$I${}_{\mathrm{2}}$, b; CH${}_{\mathrm{2}}$ICl, c; and CH${}_{\mathrm{2}}$IBr, d) and inorganic (HOI e and I${}_{\mathrm{2}}$ f) iodine source gases (in nmol m${}^{-\mathrm{2}}$ d${}^{-\mathrm{1}}$) used in SOCOL-AERv2-I.

[Figure omitted. See PDF]

As it is seen in Fig. , the sum inorganic fluxes of ISGs [HOI $+$ I${}_{\mathrm{2}}$] is roughly 3 times higher than that of organic fluxes [CH${}_{\mathrm{3}}$I $+$ CH${}_{\mathrm{2}}$I${}_{\mathrm{2}}$ $+$ CH${}_{\mathrm{2}}$ICl $+$ CH${}_{\mathrm{2}}$IBr] of ISGs that is corresponding well with .

Chemical compounds involved in the iodine chemistry scheme in SOCOL-AERv2-I are presented in Table .

Table 2

The iodine species considered in SOCOL-AERv2-I.

Species	Molar mass (g mol${}^{-\mathrm{1}}$)	Henry’s law constant (H) (mol m${}^{-\mathrm{3}}$ atm${}^{-\mathrm{1}}$)	Temperature dependence of H (K)	Reference
I	126.90	7.9 $\times$ 10${}^{\mathrm{1}}$	0
I${}_{\mathrm{2}}$	253.81	2.8 $\times$ 10${}^{\mathrm{3}}$	3900
IO	142.90	–	–	–
OIO	158.90	–	–	–
INO${}_{\mathrm{2}}$	172.91	3.0 $\times$ 10${}^{\mathrm{2}}$	0
INO	156.91	–	–	–
IONO${}_{\mathrm{2}}$	188.91	$\mathrm{\infty }$	0
HOI	143.91	4.1 $\times$ 10${}^{\mathrm{5}}$	0
HI	127.91	7.3 $\times$ 10${}^{\mathrm{16}}$	3190
IBr	206.81	2.43 $\times$ 10${}^{\mathrm{4}}$	0
ICl	162.36	1.1 $\times$ 10${}^{\mathrm{5}}$	0
I${}_{\mathrm{2}}$O${}_{\mathrm{2}}$	285.81	$\mathrm{\infty }$	0
I${}_{\mathrm{2}}$O${}_{\mathrm{3}}$	301.81	$\mathrm{\infty }$	0
I${}_{\mathrm{2}}$O${}_{\mathrm{4}}$	317.81	$\mathrm{\infty }$	0
CH${}_{\mathrm{3}}$I	141.94	2.0 $\times$ 10${}^{\mathrm{2}}$	3600
CH${}_{\mathrm{2}}$I${}_{\mathrm{2}}$	267.84	7.3 $\times$ 10${}^{\mathrm{3}}$	0
CH${}_{\mathrm{2}}$I	140.93	–	–	–
CH${}_{\mathrm{2}}$IBr	220.83	2.0 $\times$ 10${}^{\mathrm{3}}$	0
CH${}_{\mathrm{2}}$ICl	176.38	2.0 $\times$ 10${}^{\mathrm{3}}$	0

A dash means that Henry's law constant is unavailable. For IONO${}_{\mathrm{2}}$ and I${}_x$O${}_y$ Henry's law constants are assumed to be infinity and represented by a very large number.

Overall, 19 iodine species are included in the current version of the iodine chemistry scheme of SOCOL-AERv2-I: 6 iodine source gases (ISGs) and 13 product gases (PGs) including two species (HOI and I${}_{\mathrm{2}}$), which are both being emitted and chemically produced. Photolysis rates and reaction cross sections for iodine species were taken from . An exception is made only for the photolysis rates of higher iodine oxide species, such as I${}_{\mathrm{2}}$O${}_{\mathrm{2}}$, I${}_{\mathrm{2}}$O${}_{\mathrm{3}}$, and I${}_{\mathrm{2}}$O${}_{\mathrm{4}}$, where we followed the recommendations of and , who suggested using photolysis rates of iodine oxides 9 times higher than those of Cl${}_{\mathrm{2}}$O${}_{\mathrm{2}}$. This is a more simplified approach than used in CAM-chem, where cross sections for I${}_{\mathrm{2}}$O${}_{\mathrm{2}}$ and I${}_{\mathrm{2}}$O${}_{\mathrm{3}}$ are adopted from , and for I${}_{\mathrm{2}}$O${}_{\mathrm{4}}$ the used spectrum is measured at the University of Leeds . In GEOS-chem, the cross sections for higher iodine oxides are equal to those of IONO${}_{\mathrm{2}}$ based on the assumption made by . In our scheme, we checked this possibility too and found that there is no noticeable difference between using IONO${}_{\mathrm{2}}$ photolysis rates or 9*Cl${}_{\mathrm{2}}$O${}_{\mathrm{2}}$ photolysis rates as substitutes for unknown photolysis rates of higher iodine oxides. The Henry's law constant of IONO${}_{\mathrm{2}}$ is taken in SOCOL-AERv2-I to be infinity by analogous to BrONO${}_{\mathrm{2}}$ and ClONO${}_{\mathrm{2}}$ (1 $\times$ 10${}^{\mathrm{30}}$ mol m${}^{-\mathrm{3}}$ atm${}^{-\mathrm{1}}$). The Henry's law constant of I${}_x$O${}_y$ is infinity and represented by 2.65 $\times$ 10${}^{\mathrm{18}}$ mol m${}^{-\mathrm{3}}$ atm${}^{-\mathrm{1}}$ . The Henry's law constant of INO${}_{\mathrm{2}}$ is equal to that of BrNO${}_{\mathrm{2}}$ presented in .

The full iodine reaction scheme, composed of 61 gas-phase chemical reactions including 17 photolysis reactions and 4 heterogeneous reactions proceeding on tropospheric cloud ice, is presented in Table .

Table 3

The list of chemical reactions with iodine included to SOCOL-AERv2-I.

Part 1: Gas-phase chemical reactions.
Chemical reaction	A factor [cm${}^{-\mathrm{3}}$	$E_{\mathrm{a}}/$	Reference
	molecule${}^{-\mathrm{1}}$ s${}^{-\mathrm{1}}$]	$R$ (K)
I${}_{\mathrm{2}}$ $+$ O${}_{\mathrm{3}}$ $\to$ IO $+$ I $+$ O${}_{\mathrm{2}}$	3.8 $\times$ 10${}^{-\mathrm{18}}$	–
I${}_{\mathrm{2}}$ $+$ O${}_{\mathrm{3}}$ $\to$ OIO $+$ IO	3.8 $\times$ 10${}^{-\mathrm{18}}$	–
I $+$ O${}_{\mathrm{3}}$ $\to$ IO $+$ O${}_{\mathrm{2}}$	2.3 $\times$ 10${}^{-\mathrm{11}}$	870
I $+$ HO${}_{\mathrm{2}}$ $\to$ HI $+$ O${}_{\mathrm{2}}$	1.5 $\times$ 10${}^{-\mathrm{11}}$	1090
IO $+$ NO $\to$ I $+$ NO${}_{\mathrm{2}}$	9.1 $\times$ 10${}^{-\mathrm{12}}$	$-$240
IO $+$ HO${}_{\mathrm{2}}$ $\to$ HOI $+$ O${}_{\mathrm{2}}$	1.3 $\times$ 10${}^{-\mathrm{11}}$	$-$570
IO $+$ IO $\to$ OIO $+$ I	2.13 $\times$ 10${}^{-\mathrm{11}}$	$-$180
IO $+$ IO $\to$ I${}_{\mathrm{2}}$O${}_{\mathrm{2}}$	3.24 $\times$ 10${}^{-\mathrm{11}}$	$-$180
I${}_{\mathrm{2}}$ $+$ OH $\to$ HOI $+$ I	1.8 $\times$ 10${}^{-\mathrm{10}}$	–
I${}_{\mathrm{2}}$ $+$ NO${}_{\mathrm{3}}$ $\to$ I $+$ IONO${}_{\mathrm{2}}$	1.5 $\times$ 10${}^{-\mathrm{12}}$	–
I $+$ NO${}_{\mathrm{3}}$ $\to$ IO $+$ NO${}_{\mathrm{2}}$	1.0 $\times$ 10${}^{-\mathrm{10}}$	–
OH $+$ HI $\to$ I $+$ H${}_{\mathrm{2}}$O	1.6 $\times$ 10${}^{-\mathrm{11}}$	440
HOI $+$ OH $\to$ IO $+$ H${}_{\mathrm{2}}$O	2.0 $\times$ 10${}^{-\mathrm{13}}$	–
INO $+$ INO $\to$ I${}_{\mathrm{2}}$ $+$ NO $+$ NO	8.4 $\times$ 10${}^{-\mathrm{11}}$	2620
INO${}_{\mathrm{2}}$ $+$ INO${}_{\mathrm{2}}$ $\to$ I${}_{\mathrm{2}}$ $+$ NO${}_{\mathrm{2}}$ $+$ NO${}_{\mathrm{2}}$	2.9 $\times$ 10${}^{-\mathrm{11}}$	2600
IO $+$ BrO $\to$ Br $+$ I $+$ O${}_{\mathrm{2}}$	3.0 $\times$ 10${}^{-\mathrm{12}}$	$-$510
IO $+$ BrO $\to$ Br $+$ OIO	1.2 $\times$ 10${}^{-\mathrm{11}}$	$-$510
I $+$ BrO $\to$ IO $+$ Br	1.44 $\times$ 10${}^{-\mathrm{11}}$	–
IO $+$ ClO $\to$ I $+$ Cl $+$ O${}_{\mathrm{2}}$	1.175 $\times$ 10${}^{-\mathrm{12}}$	$-$280
IO $+$ O $\to$ I $+$ O${}_{\mathrm{2}}$	1.4 $\times$ 10${}^{-\mathrm{10}}$	–
O $+$ I${}_{\mathrm{2}}$ $\to$ IO $+$ I	1.4 $\times$ 10${}^{-\mathrm{10}}$	0
OH $+$ CH${}_{\mathrm{3}}$I $\to$ CH${}_{\mathrm{2}}$I $+$ H${}_{\mathrm{2}}$O	2.9 $\times$ 10${}^{-\mathrm{12}}$	1100
Cl $+$ CH${}_{\mathrm{3}}$I $\to$ CH${}_{\mathrm{2}}$I $+$ HCl	2.9 $\times$ 10${}^{-\mathrm{11}}$	1000
IO $+$ DMS $\to$ SO${}_{\mathrm{2}}$* $+$ I	3.2 $\times$ 10${}^{-\mathrm{13}}$	925
CH${}_{\mathrm{2}}$I $+$ O${}_{\mathrm{2}}$ $\to$ CH${}_{\mathrm{2}}$O $+$ IO	4.0 $\times$ 10${}^{-\mathrm{13}}$	–
IO $+$ NO${}_{\mathrm{3}}$ $\to$ OIO $+$ NO${}_{\mathrm{2}}$	9.0 $\times$ 10${}^{-\mathrm{12}}$	–
OIO $+$ NO $\to$ IO $+$ NO${}_{\mathrm{2}}$	1.1 $\times$ 10${}^{-\mathrm{12}}$	$-$542
IO $+$ OIO $\to$ I${}_{\mathrm{2}}$O${}_{\mathrm{3}}$	1.5 $\times$ 10${}^{-\mathrm{10}}$	–
OIO $+$ OIO $\to$ I${}_{\mathrm{2}}$O${}_{\mathrm{4}}$	1.5 $\times$ 10${}^{-\mathrm{10}}$	–
IO $+$ Br $\to$ I $+$ BrO	2.49 $\times$ 10${}^{-\mathrm{11}}$	–
I $+$ IONO${}_{\mathrm{2}}$ $\to$ I${}_{\mathrm{2}}$ $+$ NO${}_{\mathrm{3}}$	9.1 $\times$ 10${}^{-\mathrm{11}}$	$-$146
IO $+$ CH${}_{\mathrm{3}}$O${}_{\mathrm{2}}$ $\to$ CH${}_{\mathrm{2}}$O $+$ I $+$ HO${}_{\mathrm{2}}$	2.0 $\times$ 10${}^{-\mathrm{12}}$	–
IO $+$ ClO $\to$ ICl $+$ O${}_{\mathrm{2}}$	9.4 $\times$ 10${}^{-\mathrm{13}}$	$-$280
IO $+$ O${}_{\mathrm{3}}$ $\to$ OIO $+$ O${}_{\mathrm{2}}$	3.6 $\times$ 10${}^{-\mathrm{16}}$	–
IO $+$ OH $\to$ HO${}_{\mathrm{2}}$ $+$ I	1.0 $\times$ 10${}^{-\mathrm{10}}$	–
HI $+$ NO${}_{\mathrm{3}}$ $\to$ I $+$ HNO${}_{\mathrm{3}}$	1.3 $\times$ 10${}^{-\mathrm{12}}$	1830
I${}_{\mathrm{2}}$O${}_{\mathrm{2}}$ $+$ M $\to$ OIO $+$ I $+$ M	2.5 $\times$ 10${}^{\mathrm{14}}$	9770
I${}_{\mathrm{2}}$O${}_{\mathrm{2}}$ $+$ M $\to$ IO $+$ IO $+$ M	1.0 $\times$ 10${}^{\mathrm{12}}$	9770
I${}_{\mathrm{2}}$O${}_{\mathrm{4}}$ $+$ M $\to$ OIO $+$ OIO $+$ M	3.8 $\times$ 10${}^{-\mathrm{2}}$	–
I $+$ NO${}_{\mathrm{2}}$ $+$ M $\to$ INO${}_{\mathrm{2}}$ $+$ M	$k_{\mathrm{0}}=\mathrm{3.0}$ $\times$ 10${}^{-\mathrm{31}}$	1.0
	$k_{\mathrm{\infty }}=\mathrm{6.6}$ $\times$ 10${}^{-\mathrm{11}}$	0
INO${}_{\mathrm{2}}$ $+$ M $\to$ I $+$ NO${}_{\mathrm{2}}$ $+$ M	9.94 $\times$ 10${}^{\mathrm{17}}$	11859
IO $+$ NO${}_{\mathrm{2}}$ $+$ M $\to$ IONO${}_{\mathrm{2}}$ $+$ M	$k_{\mathrm{0}}=\mathrm{7.5}$ $\times$ 10${}^{-\mathrm{31}}$	3.5
	$k_{\mathrm{\infty }}=\mathrm{7.6}$ $\times$ 10${}^{-\mathrm{12}}$	1.5
IONO${}_{\mathrm{2}}$ $+$ M $\to$ IO $+$ NO${}_{\mathrm{2}}$ $+$ M	$k_{\mathrm{0}}=\mathrm{5.0}$ $\times$ 10${}^{-\mathrm{28}}$	14120
	$k_{\mathrm{\infty }}=\mathrm{1.9}$ $\times$ 10${}^{-\mathrm{7}}$	2.5
I $+$ NO $+$ M $\to$ INO $+$ M	$k_{\mathrm{0}}=\mathrm{1.8}$ $\times$ 10${}^{-\mathrm{32}}$	1
	$k_{\mathrm{\infty }}=\mathrm{1.7}$ $\times$ 10${}^{-\mathrm{11}}$	0

Table 3

Continued.

Part 2: Photochemical reactions.
Chemical reaction	Method	Reference
CH${}_{\mathrm{3}}$I $+$ h$\mathit{\nu }$ $\to$ CH${}_{\mathrm{3}}$ $+$ I	Look-up table
CH${}_{\mathrm{2}}$I${}_{\mathrm{2}}$ $+$ h$\mathit{\nu }$ $\to$ CH${}_{\mathrm{2}}$I $+$ I	Look-up table
I${}_{\mathrm{2}}$ $+$ h$\mathit{\nu }$ $\to$ I $+$ I	Look-up table
IO $+$ h$\mathit{\nu }$ $\to$ I $+$ O	Look-up table
OIO $+$ h$\mathit{\nu }$ $\to$ I $+$ O${}_{\mathrm{2}}$	Look-up table
INO $+$ h$\mathit{\nu }$ $\to$ I $+$ NO	Look-up table
INO${}_{\mathrm{2}}$ $+$ h$\mathit{\nu }$ $\to$ I $+$ NO${}_{\mathrm{2}}$	Look-up table
IONO${}_{\mathrm{2}}$ $+$ h$\mathit{\nu }$ $\to$ I $+$ NO${}_{\mathrm{3}}$	Look-up table
HOI $+$ h$\mathit{\nu }$ $\to$ I $+$ OH	Look-up table
HI $+$ h$\mathit{\nu }$ $\to$ I $+$ H	Look-up table
I${}_{\mathrm{2}}$O${}_{\mathrm{2}}$ $+$ h$\mathit{\nu }$ $\to$ I $+$ OIO	Look-up table
I${}_{\mathrm{2}}$O${}_{\mathrm{3}}$ $+$ h$\mathit{\nu }$ $\to$ IO $+$ OIO	Look-up table
I${}_{\mathrm{2}}$O${}_{\mathrm{4}}$ $+$ h$\mathit{\nu }$ $\to$ OIO $+$ OIO	Look-up table
IBr $+$ h$\mathit{\nu }$ $\to$ I $+$ Br	Look-up table
ICl $+$ h$\mathit{\nu }$ $\to$ I $+$ Cl	Look-up table
CH${}_{\mathrm{2}}$ICl $+$ h$\mathit{\nu }$ $\to$ I $+$ Cl	Look-up table
CH${}_{\mathrm{2}}$IBr $+$ h$\mathit{\nu }$ $\to$ I $+$ Br	Look-up table
Part 3: Heterogeneous chemical reactions on tropospheric cloud ice.
Chemical reaction	Reactive uptake	Reference
	coefficient [unitless]
HOI $+$ HI $\to$ I${}_{\mathrm{2}}$ $+$ H${}_{\mathrm{2}}$O	$\mathit{\gamma }=\mathrm{0.12}$
HOI $+$ HCl $\to$ ICl $+$ H${}_{\mathrm{2}}$O	$\mathit{\gamma }=\mathrm{0.12}$
HOI $+$ HBr $\to$ IBr $+$ H${}_{\mathrm{2}}$O	$\mathit{\gamma }=\mathrm{0.12}$
IONO${}_{\mathrm{2}}$ $+$ (H${}_{\mathrm{2}}$O) $\to$ HOI $+$ HNO${}_{\mathrm{3}}$	$\mathit{\gamma }=\mathrm{0.1}$

A factor: the pre-exponential factor. $E_{\mathrm{a}}$: the activation energy. $R$: the universal gas constant. SO${}_{\mathrm{2}}$ instead of DMSO is used (DMSO is absent in SOCOL-AERv2-I). $k_{\mathrm{0}}$: low-pressure limit (cm${}^{\mathrm{6}}$ molecule${}^{-\mathrm{2}}$ s${}^{-\mathrm{1}}$). $k_{\mathrm{\infty }}$: high-pressure limit (cm${}^{\mathrm{3}}$ molecule${}^{-\mathrm{1}}$ s${}^{-\mathrm{1}}$).

We briefly review SOCOL's chlorine and bromine reactions because of their interaction with iodine chemistry. Apart from reactions involving iodine, there are about 100 gas-phase and 11 heterogeneous reactions on sulfates and polar stratospheric clouds (PSCs) for chlorine species and about 50 gas-phase and 4 heterogeneous reactions on sulfates and polar stratospheric clouds (PSCs) for bromine species. The total gas-phase chlorine (Cl${}_{\mathrm{tot}}$) and bromine (Br${}_{\mathrm{tot}}$) in the current version of the model are $$\begin{array}{rl}& {\mathrm{Cl}}_{\mathrm{tot}}=\mathrm{ClO}+\mathrm{Cl}+\mathrm{2}\times {\mathrm{Cl}}_{\mathrm{2}}+{\mathrm{ClNO}}_{\mathrm{3}}+\mathrm{HOCl}+\mathrm{HCl}+\mathrm{2}\\ & \times {\mathrm{Cl}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{2}}+\mathrm{BrCl}+\mathrm{3}\times \mathrm{CFC}-\mathrm{11}+\mathrm{2}\times \mathrm{CFC}-\mathrm{12}+\mathrm{3}\\ & \times \mathrm{CFC}-\mathrm{113}+\mathrm{2}\times \mathrm{CFC}-\mathrm{114}+\mathrm{CFC}-\mathrm{115}+\mathrm{4}\times {\mathrm{CCl}}_{\mathrm{4}}\\ & +\mathrm{3}\times {\mathrm{CH}}_{\mathrm{3}}{\mathrm{CCl}}_{\mathrm{3}}+\mathrm{HCFC}-\mathrm{22}+\mathrm{2}\times \mathrm{HCFC}-\mathrm{141}\mathrm{B}\\ & +\mathrm{HCFC}-\mathrm{142}\mathrm{B}+\mathrm{H}-\mathrm{1211}+{\mathrm{CH}}_{\mathrm{3}}\mathrm{Cl}+\mathrm{2}\times \mathrm{HCFC}-\mathrm{21}\\ & +\mathrm{2}\times \mathrm{HCFC}-\mathrm{123}+\mathrm{HCFC}-\mathrm{31}+\mathrm{ICl}+{\mathrm{CH}}_{\mathrm{2}}\mathrm{ICl}\\ & \mathrm{and}\\ & {\mathrm{Br}}_{\mathrm{tot}}=\mathrm{Br}+\mathrm{BrO}+\mathrm{2}\times {\mathrm{Br}}_{\mathrm{2}}+\mathrm{BrCl}+{\mathrm{BrNO}}_{\mathrm{3}}+\mathrm{HBr}+\mathrm{2}\\ & \times {\mathrm{CH}}_{\mathrm{2}}{\mathrm{Br}}_{\mathrm{2}}+{\mathrm{CH}}_{\mathrm{3}}\mathrm{Br}+{\mathrm{CBrF}}_{\mathrm{3}}+\mathrm{3}\times {\mathrm{CHBr}}_{\mathrm{3}}+\mathrm{HOBr}+\mathrm{H}\\ & -\mathrm{1301}+\mathrm{H}-\mathrm{1211}+\mathrm{2}\times \mathrm{H}-\mathrm{2402}+\mathrm{IBr}+{\mathrm{CH}}_{\mathrm{2}}\mathrm{IBr}.\end{array}$$ Similar to the iodine schemes in CAM-chem and GEOS-Chem , we implemented the free molecular transfer approximation of . This allows the introduction of the iodine scavenging and deposition on sea-salt and sulfate aerosols as well as effective ice uptake (removal)/reactive ice uptake (recycling) on a surface of tropospheric cloud ice crystals . The transfer coefficient (s${}^{-\mathrm{1}}$) is calculated as follows:

4$${\displaystyle }k={\displaystyle \frac{\mathrm{1}}{\mathrm{4}}}\mathit{\gamma }⟨c⟩A,$$ where $\mathit{\gamma }$ is the effective/reactive uptake coefficient (see Table ), $A$ is the surface area density (in cm${}^{\mathrm{2}}$ cm${}^{-\mathrm{3}}$) of the particles on which the deposition occurs, $⟨c⟩=(\mathrm{8}RT/\mathit{\pi }M{)}^{\mathrm{1}/\mathrm{2}}$ is the mean thermal molecular speed (in cm s${}^{-\mathrm{1}}$) of molecules with molar mass $M$ (in kg mol${}^{-\mathrm{1}}$) at absolute temperature $T$ (in K), and $R=\mathrm{8.3145}$ J mol${}^{-\mathrm{1}}$ K${}^{-\mathrm{1}}$.

In SOCOL-AERv2-I, sea-salt aerosols are prescribed by monthly means from observational data and aqueous sulfuric acid aerosols are calculated interactively; from both, surface area densities (SADs) are available. However, there is no readily available SAD for cloud ice in SOCOL-AERv2-I. Therefore, we calculate the effective radius $R_{\mathrm{eff}}$ (in $\mathrm{\mu }\mathrm{m}$) of ice crystals following :

5$${\displaystyle }{R}_{\mathrm{eff}}=\mathit{\alpha }\exp \left(\mathit{\beta }{T}_c\right),$$ where $T_c$ is the temperature (in ${}^{\circ }$C); $\mathit{\alpha }=\mathrm{154.2}$ and $\mathit{\beta }=\mathrm{0.0152}$ for $-\mathrm{56}$ ${}^{\circ }$C $<T<$ 0 ${}^{\circ }$C; $\mathit{\alpha }=\mathrm{4.5872}$ $\times$ 10${}^{\mathrm{4}}$ and $\mathit{\beta }=\mathrm{0.117}$ for $-\mathrm{71}$ ${}^{\circ }$C $<T<$ $-$56 ${}^{\circ }$C; and $\mathit{\alpha }=\mathrm{41.65}$ and $\mathit{\beta }=\mathrm{0.0184}$ for $-$85 ${}^{\circ }$C $<T<$ $-$71${}^{\circ }$C. Following , the SAD for ice particles (in cm${}^{\mathrm{2}}$ cm${}^{-\mathrm{3}}$) is calculated from ($R_{\mathrm{eff}}$) as follows:

6$${\displaystyle }{\mathrm{SAD}}_{\mathrm{ice}}=\mathrm{6.75}\times \mathrm{IWC}/\left(\mathit{\rho }{R}_{\mathrm{eff}}\right),$$ where $\mathit{\rho }$ is the density of ice (9.167 $\times$ 10${}^{-\mathrm{4}}$ kg cm${}^{-\mathrm{3}}$) and IWC is the ice water content, i.e., the mass of ice per volume of air in the cloud (in kg cm${}^{-\mathrm{3}}$).

Since the timescale of the physical process of removal/recycling of iodine is shorter than the model time step (2 h), using an explicit integration scheme may result in excessive removal/recycling of iodine species leading to errors (such as negative concentrations). To avoid this, we decided to implement a simple implicit scheme:

7$${\displaystyle }{C}_{\mathrm{1}}=C_{\mathrm{0}}/(\mathrm{1}+K\times \mathrm{\Delta }t),$$ where $C_{\mathrm{0}}$ is the initial concentration, $K$ is the transfer coefficient, $C_{\mathrm{1}}$ is the final concentration, and $\mathrm{\Delta }t$ is the model time step for chemistry. This scheme avoids producing negative $C_{\mathrm{1}}$.

As effective uptake coefficients ($\mathit{\gamma }$), we applied ${\mathit{\gamma }}_{{\mathrm{IONO}}_{\mathrm{2}}}=\mathrm{0.01}$, ${\mathit{\gamma }}_{{\mathrm{INO}}_{\mathrm{2}}}=\mathrm{0.02}$, ${\mathit{\gamma }}_{\mathrm{HOI}}=\mathrm{0.06}$, ${\mathit{\gamma }}_{{\mathrm{I}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{2}}}={\mathit{\gamma }}_{{\mathrm{I}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{3}}}={\mathit{\gamma }}_{{\mathrm{I}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{4}}}=\mathrm{0.01}$ to model an effective uptake and removal of iodine species on sea-salt aerosols .

Since $\mathit{\gamma }$ values for sulfate aerosols are currently unknown and if we consider effective uptake on sulfate aerosols with the same gamma as for sea-salt aerosols, we would obtain modeled I${}_y$ values with a large bias ($X$ %) against observations for sulfate particles. $\mathit{\gamma }$ values for sea-salt aerosols in SOCOL-AERv2-I were divided by 100 (as the amount of iodine to be removed on tropospheric sulfate aerosols is assumed to be 100 less than for sea salt), which brought the iodine from SOCOL-AERv2-I to be closer to available observations. For effective ice uptake of iodine, $\mathit{\gamma }$ values are taken to be the same as in CAM-chem : ${\mathit{\gamma }}_{\mathrm{HOI}}=\mathrm{0.0003}$; ${\mathit{\gamma }}_{{\mathrm{IONO}}_{\mathrm{2}}}=\mathrm{0.005}$; ${\mathit{\gamma }}_{\mathrm{HI}}=\mathrm{0.02}$. The removal on all presented surfaces is operated only within the troposphere (it is confined by the tropopause height). It must be mentioned that the values of effective uptake coefficients for iodine species for different types of surfaces are highly uncertain . It is worth saying that we only consider the deposition on sea-salt and sulfate aerosols and did not implement any heterogeneous recycling of iodine species on sea-salt aerosols like those described in , , , and . It should be noted that the lack of recycling on sea-salt aerosols may affect the I${}_y$ burden and the amount of stratospheric iodine injection. Yet, it requires additional study to estimate it quantitatively. In the heterogeneous recycling mechanism on ice, reactive uptake coefficients for cloud ice are taken as follows: ${\mathit{\gamma }}_{\mathrm{HOI}}=\mathrm{0.12}$; ${\mathit{\gamma }}_{{\mathrm{IONO}}_{\mathrm{2}}}=\mathrm{0.1}$ (see supplements). Transfer coefficients for heterogeneous reactions (see Table ) are also calculated by free molecular transfer approximation but using reactive uptake coefficients. The reactive ice uptake and recycling of HOI and IONO${}_{\mathrm{2}}$ is applied in SOCOL-AERv2-I after effective ice uptake and removal of HOI, IONO${}_{\mathrm{2}}$, and HI. The sequence of removal/recycling processes is unclear but the chosen sequence shows reasonable results. Also, it should be noted here that there is no evidence of reactive uptake and recycling of iodine species in liquid clouds, and in our scheme, we use it only for ice crystals as was done by .

2.3 Conducted experiments

To evaluate the iodine chemistry scheme of SOCOL-AERv2-I as well as to estimate the influence of iodine chemistry on ozone we designed and carried out three transient numerical experiments. The first one is the control experiment where iodine emissions are set to zero. For the second experiment (1 $\times$ iodine), we applied a basic configuration with the present-day iodine emissions. These experiments are used to evaluate the veracity of iodine in SOCOL-AERv2-I and to estimate the influence of iodine chemistry on present-day ozone climatology.

To assess whether the potential intensification of iodine emissions in the future will have a tangible effect on the ozone layer, we designed a sensitivity experiment (to verify the sensitivity of ozone to iodine) in which all iodine emissions are doubled to present-day emissions (2 $\times$ iodine). In essence, it could be assumed as a worst-case scenario compared to the present-day emissions that might become closer to reality by the middle of this century if the rise of emissions presented by and continued despite no dramatic forecast of iodine emission’s evolution being made by . The iodine content in the future is also difficult to predict due to a huge discrepancy between scenarios for the future evolution of iodine precursors like tropospheric ozone and SST . Albeit, it should be said that in our study this assumption was used only to verify the sensitivity of ozone to iodine chemistry. The sensitivity of ozone to the increase in the iodine emissions was characterized by comparing experiments with 2 $\times$ and 1 $\times$ loading of iodine. All experiments were run for the 1990–2009 period including the 10-year spinup (1990–1999) from the initial conditions that is necessary for iodine to reach the quasi-equilibrium state. The spinup period was excluded from further analysis. Each experiment consists of 10 ensemble members with a 1-month perturbation of initial CO${}_{\mathrm{2}}$ concentration to get 10 different atmospheric realizations and to calculate the statistical significance of the iodine influence on ozone using Student's $t$ test. The summary of the experimental setup can be found in Table .

Table 4

The experiments with SOCOL-AERv2-I.

Name of experiment	Experiment description	Period of simulation and spinup
0 $\times$ iodine	Control run (10 ensemble members)	1990–2009 (1990–1999 spinup)
1 $\times$ iodine	Present-day emissions (10 ensemble members)	1990–2009 (1990–1999 spinup)
2 $\times$ iodine	Doubled emissions (10 ensemble members)	1990–2009 (1990–1999 spinup)

3 Results of simulation

3.1 Evaluation of the iodine from SOCOL-AERv2-I against CAM-chem and observations

The total gas-phase inorganic I${}_y$ in both 2 $\times$ iodine and 1 $\times$ iodine experiments of SOCOL-AERv2-I averaged over the tropics (20${}^{\circ }$ N–20${}^{\circ }$ S), for the 2000–2009 period and 10 ensemble members, is presented in Fig. .

Figure 2

Modeled vertical distribution of total organic (I${}_{\mathrm{org}}$) and inorganic (I${}_y$) gas-phase iodine simulated with SOCOL-AERv2-I averaged over the tropics (20${}^{\circ }$ N–20${}^{\circ }$ S), for the 2000–2009 period and 10 ensemble members. Red curve: I${}_y$ from the experiment 2 $\times$ iodine. Blue curve: I${}_y$ from the experiment 1 $\times$ iodine. Light red and blue curves: I${}_{\mathrm{org}}$ from 2 $\times$ iodine and 1 $\times$ iodine experiments, correspondingly. Shadings represent a standard deviation of tropical iodine (20${}^{\circ }$ N–20${}^{\circ }$ S).

[Figure omitted. See PDF]

The I${}_y$ was calculated as a sum of inorganic iodine compounds presented in the SOCOL-AERv2-I’s iodine scheme: I${}_y$ $=$ IO $+$ I $+$ OIO $+$ HOI $+$ 2 $\times$ I${}_{\mathrm{2}}$ $+$ HI $+$ INO $+$ INO${}_{\mathrm{2}}$ $+$ 2 $\times$ (I${}_{\mathrm{2}}$O${}_{\mathrm{2}}$ $+$ I${}_{\mathrm{2}}$O${}_{\mathrm{3}}$ $+$ I${}_{\mathrm{2}}$O${}_{\mathrm{4}}$) $+$ ICl$+$ IBr $+$ IONO${}_{\mathrm{2}}$ $+$ CH${}_{\mathrm{2}}$I. I${}_{\mathrm{org}}$ is calculated as a sum of all organic compounds in SOCOL-AERv2-I: I${}_{\mathrm{org}}$ $=$ CH${}_{\mathrm{3}}$I $+$ 2 $\times$ CH${}_{\mathrm{2}}$I${}_{\mathrm{2}}$ $+$ CH${}_{\mathrm{2}}$ICl $+$ CH${}_{\mathrm{2}}$IBr. The standard deviation of iodine between ensemble members is found to be less than 1 %. In the lower troposphere, the I${}_y$ is rapidly dropping with altitude until about 600 hPa following the washout that is most effective at this layer. In the upper troposphere (above 200 hPa), I${}_y$ increases after the initiation of the recycling heterogeneous mechanism on ice that works at these altitudes in the tropics because the availability of cloud ice is essential. The recycling on cloud ice mostly defines the amount of inorganic iodine injected into the stratosphere because it competes with the washout of reservoir species by converting reservoirs into species with lower washout rates and thus increasing the residence lifetime. The stratospheric I${}_y$ shows an increase until about 50 hPa and then stays constant because there is no deposition of iodine above this layer. The simulated stratospheric I${}_y$ by SOCOL-AER2-I agrees well with the results of the CAM-chem model in the lowermost stratosphere showing about 0.75–0.8 pptv of I${}_y$; however, it becomes higher in the middle stratosphere than in . We suggest that it might have resulted from model non-conservative transport scheme and dynamics, for example, the deep tropical convection cells over the area of iodine production that overcomes the deposition velocity enhancing the stratospheric iodine loading. Higher values of I${}_y$ in the stratosphere relative to the upper and middle troposphere can be also explained by the rapid tropospheric sinks and relatively slow transport of I${}_y$ into the troposphere. A hard upper border in SOCOL-AER2-I that prevents chemical species from going through it also might impact the profile. Nevertheless, the stratospheric I${}_y$ abundance calculated with SOCOL-AER2-I does not exceed 1.05 pptv and corresponds well to the estimation given by , although it is slightly larger than the most recent assessment from (0.8 pptv I${}_y$).

Figure 3

Modeled vertical distribution and global map of total inorganic gas-phase iodine (I${}_y$) and iodine monoxide (IO) simulated with SOCOL-AERv2-I from the experiment with present-day iodine emissions ($\mathrm{1}\times$ iodine), averaged for the 2000–2009 period and 10 ensemble members. (a, b) Zonal-mean vertical distributions of I${}_y$ and IO. (c, d) Global maps of I${}_y$ and IO averaged for 100–70 hPa region. Blue solid line in (a) and (b): annual mean tropopause height.

[Figure omitted. See PDF]

Figure 4

January–February averages of modeled and observed IO in the tropical troposphere [15${}^{\circ }$ N–15${}^{\circ }$ S] for (a) the TORERO campaign from Costa Rica (January–February 2012, 10${}^{\circ }$ N–40${}^{\circ }$ S, 250–285${}^{\circ }$ E), and (b) the CONTRAST campaign from Guam (January–February 2014, 40${}^{\circ }$ N–15${}^{\circ }$ S, 115–175${}^{\circ }$ E). Red line: IO from SOCOL-AERv2-I. Green line: IO from CAM-chem. Blue dots: IO observed by AMAX-DOAS. Shadings: IO standard deviations of all modeled and measured IO during the January–February period.

[Figure omitted. See PDF]

Figure  presents the vertical distributions of I${}_y$ and IO at all latitudes and global maps of I${}_y$ and IO averaged for the lower stratosphere (100–70 hPa). I${}_y$ and IO have the highest mixing ratio at the polar regions of the stratosphere due to transport carried out by the Brewer–Dobson circulation (BDC). The larger stratospheric values of IO at high latitudes are also related to the higher O${}_{\mathrm{3}}$ abundance at those locations. At the tropics, there are two pronounced areas with a higher (over the Indian ocean) and with a lower (over South America) I${}_y$ mixing ratio as depicted in Fig. c. Their formation might have resulted from different convection patterns (weaker/stronger) over these areas. Also, the area with a higher stratospheric I${}_y$ burden is located right over the region with higher surface emissions of HOI $/$ I${}_{\mathrm{2}}$ (see Fig. ). In the troposphere, the iodine level is decreasing toward the poles and far from the iodine source regions. The iodine distribution demonstrates the highest mixing ratio of iodine in the lower stratosphere over middle-to-high latitudes with a maximum I${}_y$ of more than 1.15 pptv in the polar region of the Northern Hemisphere and about 1 pptv in the Southern Hemisphere. The IO has two maximums of about 0.4 pptv in the lower-to-middle stratosphere and high latitudes. Also, note that IO decreases at higher levels of the stratosphere exhibiting more than 5 times less abundance than in the lower stratosphere, which might result from decreasing the efficiency of O${}_{\mathrm{3}}$–I reactions. We compared the mixing ratio of iodine compounds with some of the local measurements mentioned above. The modeled reactive IO over Scandinavia (70${}^{\circ }$ N; 20${}^{\circ }$ E) in March is $>\mathrm{0.45}$ pptv at 17 km (monthly-mean value), which is in agreement with IO simulated with a box model that was partly initialized with the IO retrieved by balloon flights (a daytime concentration of $\sim$0.65 pptv at 17 km) . Yet, the balloon-borne measured upper limit of IO mixing ratio is 0.2 pptv , which is more than 2 times lower than the IO mixing ratio simulated with SOCOL-AER2-I. SOCOL-AER2-I also captures well the I${}_y$ estimated by the box model showing a mixing ratio of 1–1.1 pptv. Also, IO simulated with SOCOL-AERv2-I corresponds well with DOAS measurements over Spitsbergen island (79${}^{\circ }$ N; 12${}^{\circ }$ E) in March , showing $>\mathrm{0.48}$ pptv in the lower stratosphere. Further, we compare the iodine monoxide (IO) obtained by SOCOL-AER2-I with the one from the CAM-chem model and the recent aircraft observations with AMAX-DOAS conducted during the TORERO and CONTRAST campaign . The results are shown in Fig. .

We compare the observed and measured IO only over the tropics (15${}^{\circ }$ N–15${}^{\circ }$ S) as extratropical measurements from both campaigns are scarcer and the comparison will be biased. Modeled IO in Fig.  has been obtained by doubling monthly averaged IO, which is fair to do for tropics as during nighttime the IO concentration is negligible. Also, the observations are given as so-called footprints to get profiles for each day; they were averaged for certain altitudes designated in Fig.  as blue points. The averaged values for Fig.  are obtained by averaging all measurements taken in the nearby altitude ranges. The IO from both models is sampled at the same altitudes as measurements to conduct an equitable comparison. It should be noted that from the CONTRAST campaign, tropical IO measurements in the upper troposphere are only available. Both SOCOL-AERv2-I and CAM-chem overestimate observations from the TORERO campaign at near-surface levels. The sharp decrease in observed and modeled IO from both models is similar until 2 km, showing compatible IO mixing ratios. At 5 km, CAM-chem shows a maximum IO mixing ratio of about 0.2 pptv that is not seen in SOCOL-AERv2-I and observations. CAM-chem generally overestimates SOCOL-AERv2-I and observations during the TORERO campaign up to about 11–12 km. After 11 km, both models overestimate the mean of observed IO from the TORERO campaign but in general stays within the uncertainty of observations. The increase in modeled IO after 12 km might follow the recycling activity on ice, but it is not seen over the region of the CONTRAST campaign. IO from both models is generally compatible with IO measurements during the CONTRAST campaign showing the IO mixing ratio of about 0.1–0.15 pptv.

Besides, we made a comparison of modeled I${}_y$ from SOCOL-AERv2-I and CAM-chem models with the I${}_y$ derived from IO $/$ I${}_y$ ratio modeled with the University of Colorado (CU) chemical box model constrained by measured temperature, pressure, chemical concentrations, particle size distributions, and photolysis frequencies. The uncertainty in derived I${}_y$ is estimated as 30 % of the IO $/$ I${}_y$ ratio including errors in the calibration of in situ and remote sensing data and accounting for differences in the spatial scales . IO was taken as an average of IO fields measured with AMAX-DOAS during both TORERO and CONTRAST campaigns . Thus, the inferred I${}_y$ is based on the measured IO and modeled IO $/$ I${}_y$. The result of this comparison is illustrated in Fig. .

Figure 5

Total gas-phase I${}_y$ modeled with SOCOL-AERv2-I and CAM-chem as well as estimated by CU chemical box model based on AMAX-DOAS IO observations in the upper troposphere and lowermost stratosphere. All model data and observations are taken only for the region of TORERO (10${}^{\circ }$ N–40${}^{\circ }$ S, 250–285${}^{\circ }$ E) and CONTRAST measurements (40${}^{\circ }$ N–15${}^{\circ }$ S, 115–175${}^{\circ }$ E). Red and green solid and dashed lines: gas-phase I${}_y$ simulated with SOCOL-AERv2-I and CAM-chem averaged for the January–February period. Light red and light green dashed lines: assumed gas-phase I${}_y$ simulated with SOCOL-AERv2-I and CAM-chem excluding approx. iodine in the particulate phase. Blue dots: gas-phase I${}_y$ modeled by the University of Colorado (CU) chemical box model. Black error bars: the uncertainty in I${}_y$ modeled by the CU box model. Blue, purple, and orange solid lines represent the empirical fit for gas-phase I${}_y$ (gas), particulate I${}_y$ (part), and total I${}_y$ (gas $+$ part), correspondingly. Shadings represent the uncertainty in the empirical fit.

[Figure omitted. See PDF]

Data were plotted with respect to the H${}_{\mathrm{2}}$O $/$ O${}_{\mathrm{3}}$ ratio as it is a good indicator to distinguish the upper troposphere air enriched in water vapor from the dehydrated air in the lower stratosphere as proposed by . To make CU box-model data clearer they were thinned out by averaging the data every 50 units of the H${}_{\mathrm{2}}$O $/$ O${}_{\mathrm{3}}$ ratio between 1000 and 100 of the H${}_{\mathrm{2}}$O $/$ O${}_{\mathrm{3}}$ ratio; every 5 units between 10 and 1; and every 1 unit between 10 and 1. The blue empirical fit is to a subset of the blue dots, so that the orange line is based on I${}_y$ (gas) and I${}_y$ (part) data in the upper troposphere, and the purple line is derived under the assumption of conversion. Equations to calculate empirical fit lines for experimental data can be found here (see Eqs. 2, 4, and 5 in , supplements).

Below the H${}_{\mathrm{2}}$O $/$ O${}_{\mathrm{3}}$ ratio of $\sim$70, both models capture well the I${}_y$ estimated from observations. Higher up, the red and green dashed lines represent the gas-phase I${}_y$ that is simulated with both models. Because the particulate iodine is not considered in these models, the total inorganic I${}_y$ in SOCOL-AERv2-I and CAM-chem presented here is only in the gas phase. The light red and green lines are the assumed gas-phase I${}_y$ if the estimated particulate iodine is excluded from the modeled I${}_y$. To exclude the unknown particulate iodine, the modeled gas-phase I${}_y$ was subtracted from the 0.87 pptv of iodine that is empirical total I${}_y$ plus uncertainty of observations (0.1 pptv) (see the orange solid line and its uncertainty in Fig. ). Such an approach can show what would be the approximate modeled gas-phase I${}_y$ if the iodine in the particulate phase was presented in models. It is seen that after the level of H${}_{\mathrm{2}}$O $/$ O${}_{\mathrm{3}}$ ratio of $\sim \mathrm{70}$, the gas-phase I${}_y$ estimated from measurements rapidly decreases and reaches the value of about 0.1 pptv at the approximate altitude of the lowermost stratosphere (H${}_{\mathrm{2}}$O $/$ O${}_{\mathrm{3}}$ ratio $<\mathrm{10}$). After excluding the estimated particulate iodine, the assumed modeled gas-phase I${}_y$ becomes very close to the one from the CU box model.

It bears mentioning that there is evidence that a certain part of gas-phase iodine undergoes partitioning to aerosol in the stratosphere. This mechanism is not fully understood due to a lack of measurements. Note that the approach used here is different from the simplified parameterization of IPART in CAM-Chem . The assumption that the overestimation of modeled gas-phase I${}_y$ against observations is because of an absence of iodine in the particulate phase is reasonable as is seen in Fig. . In this work, we do not consider particulate iodine for the analysis of ozone loss as it is out of the scope of the paper. Here, we assume that the total stratospheric I${}_y$ is only in the gas phase. Nevertheless, the total amount of iodine obtained with SOCOL-AERv2-I is in good agreement with other estimates and observations and can be used in further analysis of its effect on ozone.

3.2 Iodine chemistry effect on present-day ozone climatology

We estimate the iodine effect on present-day ozone climatology by comparing the experiment with a single loading of iodine (1 $\times$ iodine), and the control experiment neglecting iodine chemistry (0 $\times$ iodine). The contribution of iodine chemistry to present-day ozone climatology estimated by the SOCOL-AERv2-I is presented in Fig. .

Figure 6

Modeled effect of iodine chemistry on annual-mean ozone climatology averaged for the 2000–2009 period and 10 ensemble members. (a, c) Ozone changes of the case with present-day iodine emissions (1 $\times$ iodine) relative to the control run (0 $\times$ iodine), presented as zonal mean ozone mixing ratios and total ozone columns (TOC) in percent. (b, d) Corresponding absolute ozone changes in parts per billion by volume (ppbv) and Dobson units (DU), respectively. Blue solid line in (a) and (b): annual mean tropopause height. Hatching marks regions with iodine effect on ozone having a confidence level less than 95 %.

[Figure omitted. See PDF]

The relative and absolute responses of ozone to iodine chemistry activation are shown for both the ozone mixing ratio and total ozone column (TOC). The relative iodine effect on ozone is calculated as follows: $\left(\right(\mathrm{EXP}-\mathrm{REF})/\mathrm{REF})\times \mathrm{100}$, where EXP is the ozone from the experiment with 1 $\times$ loading of iodine and REF is the control run without iodine being included. The absolute difference is simply defined as EXP-REF. The crisp iodine signal in the ozone mixing ratio is observed in the lower stratosphere and intensifies over the polar regions where the effect of halogens is estimated to be higher . The peak of ozone loss resides in the lower stratosphere over the high southern latitudes, where the ozone loss reaches about 10 % or 100 ppbv. The hemispheric asymmetry of the effect over high latitudes might be caused by a difference not only in the IO${}_x$ mixing ratio but also in the mixing ratio of chlorine and bromine in active form (ClO and BrO) that in turn react with iodine. The effectiveness of iodine also depends upon the coupling with chlorine and bromine chemistry intensifying their ozone depletion cycles in the lower stratosphere (see chemical reactions in Table ). There is evidence that cross-cycles between bromine and chlorine are the dominant contribution to ozone reduction in the mid-latitudes and polar regions . The active iodine in the form of IO is also an effective reaction partner for $\mathrm{BrO}+\mathrm{ClO}$ as established by . assumed that the interaction of IO with BrO and ClO might be more effective over high latitudes because of their higher concentration in that region .

To verify whether or not cross-reactions of iodine with bromine and chlorine largely contribute to the total iodine effect on ozone, we performed an additional experiment where we set the reaction coefficients of all chemical reactions of iodine with chlorine (I-Cl) and bromine (I-Br) to zero. The contribution of iodine chemistry to present-day ozone climatology excluding I-Cl and I-Br cross-reactions estimated by the SOCOL-AERv2-I is presented in Fig. .

Figure 7

Modeled effect of iodine chemistry on annual-mean ozone climatology excluding I-Cl and I-Br cross-halogen reactions, averaged for the 2000–2009 period. Presented ozone changes of the case with present-day iodine emissions (1 $\times$ iodine) relative to the control run (0 $\times$ iodine). Blue solid line: annual mean tropopause height.

[Figure omitted. See PDF]

Comparing Figs. a and , one can say that the inter-halogen cross-reactions are minor contributors to tropospheric ozone responsible for only about 20 % of total iodine-induced ozone loss. However, in the lower stratosphere, the contribution of these cross-cycles is much higher, covering the 60 %–70 % of total iodine effect on ozone in the Northern Hemisphere. In the Southern Hemisphere, the contribution is found to be less than $\sim \mathrm{30}$ %. Thus, our model suggests that cross-cycles of I with Cl and Br might largely affect ozone and be responsible for more than a half of total iodine-induced ozone reduction in the lower stratosphere, which agrees with and .

The iodine-mediated ozone loss in the tropical lower stratosphere ranged between 4 %–5 % or 20 ppbv below 20 km. The ozone destruction caused by iodine chemistry is not so pronounced over the tropics because the extremely low temperature in the cold trap on the tropical tropopause that resulted from adiabatic cooling is even lower than the temperature over high latitudes at the same height making catalytic cycles less effective.

It is important to mention that the effects of iodine chemistry are confined mostly to the lower stratosphere and are greatly decreasing from the middle to the upper stratosphere. The “glasses-like” iodine effect in the upper stratosphere similar to the one of ClO${}_x$ cycle is not observed. The IO${}_x$-catalytic cycle is ineffective in the upper stratosphere similarly to BrO${}_x$ cycle despite the presence of atomic iodine up there in reasonable concentrations. It might be because iodine reservoirs are much more unstable than those of chlorine and, similarly to bromine, the iodine species in lower altitudes will be more likely in active form (IO) than that of chlorine will be in the form of ClO . So, the probability of the terminal reaction $\mathrm{ClO}+\mathrm{O}$ is higher than those of $\mathrm{IO}+\mathrm{O}$ or $\mathrm{BrO}+\mathrm{O}$ in the upper stratosphere. It should be noted that SOCOL-AERv2-I demonstrates the decrease in the IO mixing ratio in the upper stratosphere (see Fig. b). Hence, the effect of iodine chemistry on upper stratospheric ozone loss is negligibly small.

The tropospheric effect is $\sim$4–5 ppbv and maximizes over the tropics where iodine sources are mostly emitted from the ocean thanks to its higher temperature. About 6 %–8 % of tropospheric ozone loss is comparable to what was reported by , but surface iodine emissions were a bit higher. It is also in agreement with an estimation made by . Tropospheric ozone loss in SOCOL-AERv2-I is a bit higher than in CAM-chem where the iodine-induced ozone loss does not exceed 2–3 ppbv .

The total ozone column (TOC) is affected by iodine mostly over high latitudes (see Fig. c and d). The highest impact of iodine on climatological TOC is over high latitudes of the Southern Hemisphere showing the TOC loss of about 4 % or 11–12 DU. Over the Northern Hemisphere, the iodine effect on TOC does not exceed 3 %.

3.3 Ozone response to the increased iodine emissions

To estimate the consequences for ozone from the continuous increase in iodine emissions, we compare ozone from the experiment with doubled emissions to that from the experiment with observed (or present-day) emissions. The results are shown in Fig. 8.

Figure 8

Modeled effect of 2-fold iodine chemistry on annual-mean ozone climatology averaged for the 2000–2009 period and 10 ensemble members. (a, c) Ozone changes of the case with 2-fold iodine chemistry (2 $\times$ iodine) relative to present-day iodine emissions (1 $\times$ iodine) in percent. (b, d) Corresponding absolute ozone changes in parts per billion by volume (ppbv) and Dobson units (DU), respectively. Blue solid line in (a) and (b): annual mean tropopause height. Hatching marks regions with iodine effect on ozone having a confidence level less than 95 %.

[Figure omitted. See PDF]

In this case, the EXP was taken to be the experiment with 2 $\times$ loading of iodine whereas the REF is the one with observed (1 $\times$) iodine emissions. In the troposphere, a 2-fold increase in emissions leads to the ozone loss of about 6 %–8 %, which is similar to what is seen in Fig. . In the stratosphere, the contribution of additional iodine is different. The southern hemispheric maximum is weakened and shifted to the middle latitude showing the ozone loss of up to 7 % or 50 ppbv. The iodine contribution to ozone loss in the lowermost northern polar stratosphere also shows a similar pattern and magnitude of the effect on the climatology in Fig. , showing about 50 ppb of ozone loss. For the Northern Hemisphere, the effect might be characterized as a linear kind too and the intensification of ozone loss by a factor of 2 can be expected. The iodine-induced total column ozone loss would enhance by 2 %–3 % following the 2-fold increase in iodine emissions.

Thus, we can expect that a 2-fold increase in iodine injection into the atmosphere would lead to a mostly linear increase in the ozone loss over the troposphere and lower stratosphere. The reason why we do not see the linearly changed effect over the Southern Hemisphere is possibly related to the saturation effect when the ozone was almost destroyed even for smaller iodine loading. Hence, it can be predicted that the iodine effect on ozone in the lower atmosphere, if the assumed negative iodine scenario plays out in the future, would simply hinge on a factor of the increasing iodine injection into the atmosphere.

3.4 Ozone response to organic and inorganic iodine sources

We also address the effect of organic vs inorganic iodine sources on tropical ozone. To do this, we repeated the experiment with observed (or present-day) emissions twice, setting either organic or inorganic surface emissions to zero. The results are shown in Fig. 9.

Figure 9

Modeled vertical profiles of O${}_{\mathrm{3}}$ reduction by iodine from organic, inorganic, and total surface emissions averaged over the tropics (20${}^{\circ }$ N–20${}^{\circ }$ S), for the 2000–2009 period and 10 ensemble members. Red curve: difference in O${}_{\mathrm{3}}$ relative to the control run (0 $\times$ iodine), if only organic iodine source gases are considered; blue curve: a percentage difference in O${}_{\mathrm{3}}$ relative to the control run (0 $\times$ iodine), if inorganic iodine source gases are only considered; green curve: a percentage difference in O${}_{\mathrm{3}}$ relative to the control run (0 $\times$ iodine), if both organic and inorganic iodine source gases are considered. Shadings represent a standard deviation of ensemble members. The results have a confidence level more than or equal to 95 %.

[Figure omitted. See PDF]

In the lower troposphere, the iodine from inorganic emissions is responsible for $\sim \mathrm{75}$ % of total iodine effect on ozone, and the contribution of iodine from organic sources is $\sim \mathrm{25}$ %, as expected . In the upper troposphere and lower stratosphere, the ozone responses to iodine from both kinds of sources are closer, but still the impact of inorganic sources prevails, i.e., $\sim \mathrm{60}$ % and $\sim \mathrm{40}$ %, correspondingly. At 50 hPa, the contributions of iodine originating from organic and inorganic sources to total iodine-induced ozone reduction are compatible.

4 Discussion and conclusions

In this study, we describe the new version of the chemistry–climate model SOCOL-AERv2 improved with the addition of the iodine chemistry module. The iodine chemistry scheme in SOCOL-AERv2-I was developed based on up-to-date knowledge about atmospheric iodine. We performed a set of numerical experiments to test the fidelity of the developed iodine chemistry scheme and to estimate iodine contribution to ozone depletion. The model results show about 0.75 pptv of iodine in the lowermost tropical stratosphere in agreement with previous estimations . A gradual increase in I${}_y$ up to 1.05 pptv in the stratosphere might be related to the dynamical features, the effect of tropospheric sinks and/or recycling, or the rigid upper border in the model atmosphere. The comparison of modeled and observed IO within the tropical troposphere showed that IO from SOCOL-AERv2-I is in a generally good agreement with AMAX-DOAS observations. The additional comparison of total gas-phase inorganic I${}_y$ with the values determined by the CU box model based on AMAX-DOAS observations of IO showed that the model reproduces the I${}_y$ in the upper troposphere very well, while in the stratosphere, I${}_y$ is much overestimated due to the absence of simulated iodine in the particulate phase. If the assumed particulate iodine was excluded from the modeled stratospheric I${}_y$, SOCOL-AERv2-I corresponds well with I${}_y$ from the CU box model. The simulation of particulate iodine is a subject for future studies as the mechanism of its formation is still not fully understood due to a lack of measurements.

The estimated contribution of iodine chemistry on the lower stratospheric ozone is higher than those discussed in , showing up to 10 % of the lower stratospheric ozone loss driven by iodine chemistry. It should be noted that reported only 0.15 pptv of I${}_y$ injected into the stratosphere, which is more than 5 times less than in SOCOL-AERv2-I and CAM-chem. In the lower troposphere, the share of ozone loss induced by iodine originating only from inorganic sources is estimated to be 75 %; it is estimated to be 25 % if considering only organic sources. At 50 hPa, contributions of iodine from organic and inorganic sources to total iodine-induced ozone loss are found to be compatible. We also verified that even if the concentration of iodine is much less than other halogens it could play a noticeable role in the lower stratospheric ozone depletion especially over high latitudes. Nevertheless, negative lower stratospheric changes recently found by might be driven by iodine chemistry only in the lowermost stratosphere as the iodine effect in the extratropical lower stratosphere is not supposed to produce such a broad signal. At the same time, SOCOL-AERv2-I suggests that cross-cycles of iodine with bromine and chlorine might sufficiently affect ozone covering the 50 %–70 % of total iodine contribution to ozone reduction in the lower stratosphere at mid-to-high latitudes with the strongest contribution in the Northern Hemisphere. It should be noted that the $\mathrm{I}-\mathrm{Br}/\mathrm{I}-\mathrm{Cl}$ ratio in the total effect of iodine on ozone requires additional studying. The upper stratospheric ozone is not affected by iodine chemistry, similar to the impact of bromine chemistry, as presumably the iodine species in the upper stratosphere are hardly in more active form (IO) than that of chlorine due to longer lifetime of chlorine precursors . Our results also demonstrate a lowering of the IO mixing ratio in the upper stratosphere. Hence, the net effect is that iodine is relatively more important in the lower stratosphere although the abundance of total inorganic gas-phase iodine is expected to be similar throughout the stratosphere (see Fig. ).

Also, we would like to address here some shortcomings of the current iodine scheme and further updates that are anticipated to increase the accuracy of iodine simulations in chemistry–climate models. A strong local decrease in free tropospheric ozone can be expected from iodine collected inside iodized aerosol particles from deserts and oceans which can reach the above-cloud troposphere where they release iodine . However, such complex aerosol iodine chemistry cannot be properly simulated with current chemistry–climate models . The recent studies of , , and suggested that iodine species can make new aerosol particles that are big enough to be cloud condensation nuclei (CCNs). Polar ice-melting in turn leads to the increase in the atmospheric amount of iodine which may enhance the formation of CCNs. These findings are worth studying further using global climate models with advanced aerosol chemistry and cloud microphysics. Additionally, higher iodine oxide species presented in our study by I${}_{\mathrm{2}}$O${}_{\mathrm{2}}$, I${}_{\mathrm{2}}$O${}_{\mathrm{3}}$, and I${}_{\mathrm{2}}$O${}_{\mathrm{4}}$ might decompose to form iodine oxoacids that grow further to become CCNs . However, the formation of iodine oxoacids is still an open question that needs further study . Here, we did not consider any iodine in aerosol form, and these species are represented only in the gas phase. Also, as was mentioned above, in this work we used a simplified approach for photolysis of higher-order iodine oxides. However, the cross sections for photolysis of these species are recently measured and can be used to increase the accuracy of simulation . Recent field evidence indicates that the recycling of iodine on sea-salt aerosol, and perhaps other types of aerosols, may be much faster than currently represented . Also, we did not include CF${}_{\mathrm{3}}$I that also could modify the total concentration of iodine in the atmosphere. However, based on the recent studies, it will not substantially impact the stratospheric ozone loss, showing even less impact than that of CH${}_{\mathrm{3}}$I , and will mostly affect the tropical and northern mid-latitude tropospheric ozone because of higher concentrations of pollutants . All organic iodine emissions in our scheme are prescribed. However, they can be interactively calculated in the model utilizing ocean biogenic sources . It can also be embodied in the next-generation Earth system model where the ocean biosphere is interactively calculated. As was mentioned above, the sea-surface iodide that is the precursor for HOI $/$ I${}_{\mathrm{2}}$ fluxes has some uncertainty in the models and observations . However, argued that the iodide field does not necessarily solve any uncertainty in the HOI $/$ I${}_{\mathrm{2}}$ flux calculation or other impacts on these fluxes. Nonetheless, there is a prediction for increased iodide that will potentially impact the iodine abundance in the future . Also, there are region-specific parameterizations for sea surface iodide concentration that can be implemented in the next version of the iodine scheme in SOCOL-AERv2-I to increase the reliability of abiotic iodine emissions . Recent investigations of the chemical basis of the HOI $/$ I${}_{\mathrm{2}}$ source can help improve the generalization of empirical source functions .

Our sensitivity study showed that the contribution of increased iodine to ozone is almost linear compared to the present-day iodine. To simulate the reliable future impact of iodine on ozone, the recent estimations on future iodine emissions based on RCP scenarios can be used . Tropospheric ozone content in SOCOL-AERv2-I is overestimated compared to other models , which affects the sea-surface deposition of O${}_{\mathrm{3}}$ and its concentration inside the marine boundary layer; hence the accuracy of simulated iodine emissions. There are plans to fix this problem in future versions of SOCOL.

One of the most controversial parts of atmospheric iodine studies is the scrutiny of the role of volcanic iodine in stratospheric chemistry and its effect on ozone. The volcanic iodine impact is worth studying since some of the powerful volcanoes are supposedly capable of directly injecting the iodine into the stratosphere and supposing to result in negative and long-lasting implications for the ozone layer . However, there is no solid evidence that volcanoes can inject a sufficient amount of iodine into the atmosphere , and it is needed to organize the measurement campaigns to make estimations of emitted volcanic iodine in a more precise way. The results of this work showed the highest impact of iodine on ozone in the lowermost stratosphere at high latitudes. This finding indicates the necessity of having broad measurements of iodine species in this region.

We also stress that the iodine can presumably be more noteworthy in the future. The intensified Brewer–Dobson circulation might bring more iodine into the stratosphere in the future than today. Also, the further increase in sea surface temperature due to global warming, near-surface ozone, and sea surface iodide concentrations could be the reason for the intensification of iodine emissions in the future, making the atmospheric amount of iodine vastly higher . The effectiveness of iodine for ozone destruction is found to be stable in future warming scenarios and therefore its relative importance increases relative to the other halogens .

All of this inspires further efforts to better characterize the iodine in the atmosphere and its impact on ozone loss.

Alongside the further improvements of iodine chemistry simulations in chemistry–climate models, it is necessary to overcome the scarcity of global measurements of iodine chemistry especially in the upper troposphere and lower stratosphere to increase the accuracy of estimations for iodine impact on ozone loss and to make better predictions of the future ozone evolution.

Code and data availability

The SOCOL-AERv2-I code is available here: 10.5281/zenodo.4844994 upon request to the corresponding author. However, since the ECHAM5 model is a part of SOCOL-AERv2-I, the license agreement must be signed before using the code (http://www.mpimet.mpg.de/en/science/models/license/, last access: 11 September 2021). The SOCOL-AERv2-I simulation data can be accessed here: 10.5281/zenodo.4820523 . CU-AMAX-DOAS CONTRAST IO data are available at https://data.eol.ucar.edu/dataset/383.023 (last access: 11 September 2021; 10.5065/D6F769MF,). CU-AMAX-DOAS TORERO IO data are available at https://data.eol.ucar.edu/dataset/352.082 . The CU box model data are available here: 10.5281/zenodo.4916787 .

Author contributions

AKD conducted all simulations, visualized the results, and drafted the paper. ER, TiS, TE, and TG analyzed the simulations. ASL, CAC, and RPF provided the CAM-chem data and assisted in the analysis. ToS provided iodine organic fluxes for boundary conditions. RV and TKK provided CU box model data. This study was conceptualized and supervised by ER and TP. All authors participated in the model development, discussed the results, and contributed to writing and editing the paper.

Competing interests

The contact author has declared that neither they nor their co-authors have any competing interests.

Disclaimer

Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Acknowledgements

Arseniy Karagodin-Doyennel, Eugene Rozanov, Timofei Sukhodolov, and Tatiana Egorova express gratitude to the Swiss National Science Foundation for supporting this research through the no. 200020-182239 project POLE (Past and Future Ozone Layer Evolution). Rainer Volkamer acknowledges funding by the US National Science Foundation (awards AGS-2027252, AGS-1261740, AGS-1104104). Rainer Volkamer is currently an ETH guest professor and recipient of a Swiss National Science Foundation fellowship (award 199407). The authors thank the Center for Climate Systems Modeling (C2SM) for their support and ETH’s High Performance Computing Center (ID SIS) for the possibility to use the Euler Linux cluster to conduct our numerical experiments. The work of Eugene Rozanov and Timofei Sukhodolov on the analysis of the results was performed in the SPbSU Ozone Layer and Upper Atmosphere Research Laboratory supported by the Ministry of Science and Higher Education of the Russian Federation under grant 075-15-2021-583. We thank the editor and two anonymous referees for their constructive comments, which helped us to improve the paper.

Financial support

This research has been supported by the Schweizerischer Nationalfonds zur Förderung der Wissenschaftlichen Forschung (grant no. 200020_182239 (POLE)). Rainer Volkamer was supported by the US National Science Foundation (awards AGS-2027252, AGS-1261740, AGS-1104104). The research has been partly supported by the Ministry of Science and Higher Education of the Russian Federation under grant 075-15-2021-583.

Review statement

This paper was edited by Patrick Jöckel and reviewed by two anonymous referees.