Articles | Volume 14, issue 10
Model description paper
29 Oct 2021
Model description paper |  | 29 Oct 2021

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

Arseniy Karagodin-Doyennel, Eugene Rozanov, Timofei Sukhodolov, Tatiana Egorova, Alfonso Saiz-Lopez, Carlos A. Cuevas, Rafael P. Fernandez, Tomás Sherwen, Rainer Volkamer, Theodore K. Koenig, Tanguy Giroud, and Thomas Peter

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 (Solomon1999). 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., Newman et al.2009; Egorova et al.2013; WMO2018). 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 (Solomon et al.1994; Saiz-Lopez et al.2012). However, the ozone depletion efficiency of iodine-containing species on a per-molecule basis is hundreds of times higher than that of chlorine (Solomon et al.1994; Gutmann et al.2018; Koenig et al.2020), 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  (Saiz-Lopez et al.2015) and slowing the recovery of the ozone layer (Koenig et al.2020). 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 (∼1 pptv) and because there is almost no direct anthropogenic production of iodine-containing species (Fuge and Johnson2015).

The study of iodine source gases (ISGs) began in the late 1970s. The main organic source of atmospheric iodine is found to be iodomethane (CH3I) as suggested by Lovelock and Maggs (1973) and Chameides and Davis (1980). Photochemical production in ocean surface water involving phytoplankton and algae is the main source of organic and inorganic iodine that enters the atmosphere (Saiz-Lopez et al.2012 and references therein). Iodocarbons are further emitted by wetland plants (Manley et al.2007), biomass burning (Akagi et al.2011), rice paddies (Lee-Taylor and Redeker2005), and volcanoes (Bureau et al.2000; Aiuppa et al.2005). 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 CH3I is long enough to allow a small fraction of it to reach the upper troposphere via tropical convection with an estimated mixing ratio ∼0.1 pptv in the tropical tropopause layer (WMO2018). Other organic ISGs, such as CH2I2, CH2ICl, CH2IBr, and C2H5I, 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, I2 (Saiz-Lopez and Plane2004), and hypoiodous acid, HOI (Carpenter et al.2013). These are formed through reactions of near-surface ozone with oceanic iodide (I), i.e., O3+ H++ I HOI + O2 and HOI + H++ I I2+ H2O, which are possibly responsible for up to 75 % of total iodine emissions (Carpenter et al.2013; Prados-Roman et al.2015b). In conjunction with these findings, (MacDonald et al.2014) used a numerical parametrization of sea-surface iodide in numerical modeling of HOI and I2 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. Cuevas et al. (2018) and Legrand et al. (2018) argued that the atmospheric iodine loading must have increased by at least a factor of 3 since the 1950s, because of increasing anthropogenic NOx 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 (Archibald et al.2020), 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 (Iglesias-Suarez et al.2020). In addition, the effectiveness of iodine for ozone depletion (O3 destruction by one I atom relative to O3 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 (Klobas et al.2021).

Using a two-dimensional model of atmospheric chemistry and dynamics, Solomon et al. (1994) 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 Pundt et al. (1998), 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) Wittrock et al. (2000) found significant concentrations of IO between 0.65 and 0.8 pptv (±0.2 pptv) over Spitsbergen Island. Recent aircraft campaigns, also applying DOAS (Baidar et al.2013) over the West Pacific during CONTRAST (Convective Transport of Active Species in the Tropics) in January–February 2014 (Pan et al.2017) and over the East Pacific during TORERO (Tropical Ocean Troposphere Exchange of Reactive Halogen Species and Oxygenated VOC) in January–February 2012 (Volkamer et al.2015)1, 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 (Schönhardt et al.2008), though with a very high level of uncertainty.

Modeling of atmospheric iodine chemistry started with Vogt et al. (1996, 1999) 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+ ICl + H2O, HOI + Br+ H+ IBr + H2O). 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). Dix et al. (2013) 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. Ordóñez et al. (2012) 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 (Prados-Roman et al.2015a, b). Later, an updated version of the GEOS-Chem model (Sherwen et al.2016a) demonstrated 3.83 Tg (I)/year of total ISGs emissions, when using the sea-surface iodide field from Chance et al. (2014). Sherwen et al. (2016c) used the sea-surface values from MacDonald et al. (2014), which possibly underestimates the sea-surface iodide emissions in comparison to observations (Sherwen et al.2019), suggesting only 2.7 Tg (I)/year of total emissions. In order to reproduce the IO observations of Volkamer et al. (2015), Saiz-Lopez et al. (2015) introduced the heterogeneous recycling of iodine reservoirs (HOI and IONO2) on ice crystals and implemented it in CAM-chem, assuming it to be analogy to bromine recycling on ice crystals in the upper troposphere (Aschmann and Sinnhuber2013). Accordingly, the total iodine injected into the stratosphere is currently thought to be up to 0.8 pptv (WMO2018), i.e., about 5 times larger than that from the previously published assessment (WMO2014). 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 Hossaini et al. (2015) obtained using the 3-D chemistry transport model TOMCAT suggest that iodine lowers stratospheric ozone only by ∼3 % and total column ozone by 0.5 %. However, this model did not take the recycling mechanism described in Saiz-Lopez et al. (2015) into consideration, which increases the stratospheric iodine levels significantly beyond what was assumed by Hossaini et al. (2015). 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 (Stenke et al.2013; Revell et al.2015) coupled to a size-resolving sulfate aerosol module (AER) (Weisenstein et al.1997) along with other important modifications for chemistry and deposition. The AER module of SOCOL was established by Sheng et al. (2015) (CCM SOCOL-AERv1). The CCM SOCOL-AERv1 was substantially updated by Feinberg et al. (2019) (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) (Manzini et al.2006), which has been interactively coupled to the MEZON atmospheric chemistry-transport module (Model for the Evaluation of OZONe Trends) (Egorova et al.2003). The coupling takes account of radiative forcing caused by ozone, H2O, N2O, CH4, 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  (Lin and Rood1996). The standard SOCOL-AERv2 utilizes the Gaussian transform horizontal grid with T42 triangular truncation (64 latitudes × 128 longitudes) splitting the model space into grid cells of 2.5 × 2.5 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 (Rayner et al.2003). 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 Feinberg et al. (2019).

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 1.

Carpenter et al. (2013)Carpenter et al. (2013)Ordóñez et al. (2012)Ordóñez et al. (2012)Ordóñez et al. (2012)Ordóñez et al. (2012)

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

Download Print Version | Download XLSX

Following Ordóñez et al. (2012), the organic iodocarbons have been obtained from the inventory of Bell et al. (2002) for CH3I and from the 1-D model estimates of Jones et al. (2010) for CH2I2, CH2ICl, and CH2IBr. Organic ISGs have been parameterized in Ordóñez et al. (2012) by a biogenic chlorophyll-a (chl-a) dependent source in the tropical oceans (McClain et al.2004). 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 × 2.5, including updates to both iodine and bromine chemistry from Sherwen et al. (2016c) and Schmidt et al. (2016), and interpolated on the SOCOL-AERv2-I horizontal grid (T42). The inorganic HOI / I2 fluxes are interactively calculated in the model using the numerical parametrization of Carpenter et al. (2013) and global sea surface iodide concentration calculated following MacDonald et al. (2014). To derive HOI / I2 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 FHOI and FI2 (in nmol m−2 d−1) as a function of the surface ozone mixing ratio χO3,surf (in ppbv) are given by


where Ws is the surface wind speed (in m s−1) and MIaq- is the sea surface iodide concentration expressed as molarity (mol dm−3)

(3) M I aq - = 1.46 × 10 6 exp ( - 9134 / T ss )

with sea surface temperature Tss in K.

It is important to note that this parameterization is intended for calculating I2 and HOI fluxes at wind speeds below 14 m s−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 I2 and HOI escaping into the atmosphere; therefore fluxes decrease with the wind speed (MacDonald et al.2014). The obtained inorganic ISGs fluxes in the SOCOL-AERv2-I exceed GEOS-Chem and CAM-Chem emissions because SOCOL-AERv2-I overestimates near-surface O3 (Revell et al.2018). When we compared the ground-level ozone fields from SOCOL and GEOS-chem, we found that for most regions where HOI / I2 emissions are produced, the overestimation of ground-level ozone in SOCOL is 100 %. To ameliorate this bias, the surface ozone used within the parameterization was scaled by 0.5. This halved the emission, making HOI / I2 fluxes comparable: 2.4 Tg (I)/year in SOCOL-AERv2-I and 2.2 Tg (I)/year in GEOS-chem (Sherwen et al.2016c). Since there is an uncertainty in inorganic iodine emissions (Chance et al.2014; Sherwen et al.2016c), 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. 1.

Figure 1Annual-mean emission fluxes of organic (CH3I, a; CH2I2, b; CH2ICl, c; and CH2IBr, d) and inorganic (HOI e and I2 f) iodine source gases (in nmol m−2 d−1) used in SOCOL-AERv2-I.

As it is seen in Fig. 1, the sum inorganic fluxes of ISGs [HOI + I2] is roughly 3 times higher than that of organic fluxes [CH3I + CH2I2+ CH2ICl + CH2IBr] of ISGs that is corresponding well with Carpenter et al. (2013).

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

Burkholder et al. (2015)Burkholder et al. (2015)Ordóñez et al. (2012)Burkholder et al. (2015)Badia et al. (2019)Badia et al. (2019)Wagman et al. (1989)Wagman et al. (1989)Badia et al. (2019)Badia et al. (2019)Badia et al. (2019)Burkholder et al. (2015)Burkholder et al. (2015)Hilal et al. (2008)Hilal et al. (2008)

Table 2The iodine species considered in SOCOL-AERv2-I.

A dash means that Henry's law constant is unavailable. For IONO2 and IxOy Henry's law constants are assumed to be infinity and represented by a very large number.

Download Print Version | Download XLSX

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 I2), which are both being emitted and chemically produced. Photolysis rates and reaction cross sections for iodine species were taken from Burkholder et al. (2015). An exception is made only for the photolysis rates of higher iodine oxide species, such as I2O2, I2O3, and I2O4, where we followed the recommendations of Davis et al. (1996) and Vogt et al. (1999), who suggested using photolysis rates of iodine oxides 9 times higher than those of Cl2O2. This is a more simplified approach than used in CAM-chem, where cross sections for I2O2 and I2O3 are adopted from Gómez-Martín et al. (2005), and for I2O4 the used spectrum is measured at the University of Leeds (Saiz-Lopez et al.2014). In GEOS-chem, the cross sections for higher iodine oxides are equal to those of IONO2 (Sherwen et al.2016a) based on the assumption made by Bloss et al. (2010). In our scheme, we checked this possibility too and found that there is no noticeable difference between using IONO2 photolysis rates or 9*Cl2O2 photolysis rates as substitutes for unknown photolysis rates of higher iodine oxides. The Henry's law constant of IONO2 is taken in SOCOL-AERv2-I to be infinity by analogous to BrONO2 and ClONO2 (1 × 1030 mol m−3 atm−1). The Henry's law constant of IxOy is infinity and represented by 2.65 × 1018 mol m−3 atm−1 (Badia et al.2019). The Henry's law constant of INO2 is equal to that of BrNO2 presented in Ordóñez et al. (2012).

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 3.

Ordóñez et al. (2012)Ordóñez et al. (2012)Burkholder et al. (2015)Ordóñez et al. (2012)Burkholder et al. (2015)Burkholder et al. (2015)Saiz-Lopez et al. (2014)Badia et al. (2019)Burkholder et al. (2015)Saiz-Lopez et al. (2014)Badia et al. (2019)Badia et al. (2019)Saiz-Lopez et al. (2014)Burkholder et al. (2015)Burkholder et al. (2015)Badia et al. (2019)Badia et al. (2019)Saiz-Lopez et al. (2014)Saiz-Lopez et al. (2014)Burkholder et al. (2015)Burkholder et al. (2015)Burkholder et al. (2015)Burkholder et al. (2015)Saiz-Lopez et al. (2015)Gravestock et al. (2010)Saiz-Lopez et al. (2014)Saiz-Lopez et al. (2014)Badia et al. (2019)Badia et al. (2019)Saiz-Lopez et al. (2014)Saiz-Lopez et al. (2014)Saiz-Lopez et al. (2014)Saiz-Lopez et al. (2014)Saiz-Lopez et al. (2014)Saiz-Lopez et al. (2014)Saiz-Lopez et al. (2014)Ordóñez et al. (2012)Ordóñez et al. (2012)Badia et al. (2019)Burkholder et al. (2015)Badia et al. (2019)Burkholder et al. (2015)Burkholder et al. (2015)Burkholder et al. (2015)Burkholder et al. (2015)Burkholder et al. (2015)Burkholder et al. (2015)Burkholder et al. (2015)Burkholder et al. (2015)Burkholder et al. (2015)Burkholder et al. (2015)Burkholder et al. (2015)Burkholder et al. (2015)Burkholder et al. (2015)Davis et al. (1996)Davis et al. (1996)Davis et al. (1996)Burkholder et al. (2015)Burkholder et al. (2015)Burkholder et al. (2015)Burkholder et al. (2015)Saiz-Lopez et al. (2015)Saiz-Lopez et al. (2015)Saiz-Lopez et al. (2015)Saiz-Lopez et al. (2015)

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

A factor: the pre-exponential factor. Ea: the activation energy. R: the universal gas constant. SO2 instead of DMSO is used (DMSO is absent in SOCOL-AERv2-I). k0: low-pressure limit (cm6 molecule−2 s−1). k: high-pressure limit (cm3 molecule−1 s−1).

Download XLSX

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 (Cltot) and bromine (Brtot) in the current version of the model are


Similar to the iodine schemes in CAM-chem (Ordóñez et al.2012) and GEOS-Chem (Sherwen et al.2016a), we implemented the free molecular transfer approximation of McFiggans et al. (2000). 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 (Fernandez et al.2014; Saiz-Lopez et al.2014, 2015). The transfer coefficient (s−1) is calculated as follows:

(4) k = 1 4 γ c A ,

where γ is the effective/reactive uptake coefficient (see Table 3), A is the surface area density (in cm2 cm−3) of the particles on which the deposition occurs, c=(8RT/πM)1/2 is the mean thermal molecular speed (in cm s−1) of molecules with molar mass M (in kg mol−1) at absolute temperature T (in K), and R=8.3145 J mol−1 K−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 Reff (in µm) of ice crystals following Heymsfield et al. (2014):

(5) R eff = α exp ( β T c ) ,

where Tc is the temperature (in C); α=154.2 and β=0.0152 for −56C <T< 0 C; α=4.5872× 104 and β=0.117 for −71C <T<56 C; and α=41.65 and β=0.0184 for 85 C <T<71C. Following Holmes et al. (2019), the SAD for ice particles (in cm2 cm−3) is calculated from (Reff) as follows:

(6) SAD ice = 6.75 × IWC / ( ρ R eff ) ,

where ρ is the density of ice (9.167 × 10−4 kg cm−3) and IWC is the ice water content, i.e., the mass of ice per volume of air in the cloud (in kg cm−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) C 1 = C 0 / ( 1 + K × Δ t ) ,

where C0 is the initial concentration, K is the transfer coefficient, C1 is the final concentration, and Δt is the model time step for chemistry. This scheme avoids producing negative C1.

As effective uptake coefficients (γ), we applied γIONO2=0.01, γINO2=0.02, γHOI=0.06, γI2O2=γI2O3=γI2O4=0.01 to model an effective uptake and removal of iodine species on sea-salt aerosols (Ordóñez et al.2012).

Since γ 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 Iy values with a large bias (X %) against observations for sulfate particles. γ 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, γ values are taken to be the same as in CAM-chem (Saiz-Lopez et al.2014): γHOI=0.0003; γIONO2=0.005; γHI=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 (Saiz-Lopez et al.2014). 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 Vogt et al. (1996, 1999), Ordóñez et al. (2012), Saiz-Lopez et al. (2014), and Tham et al. (2021). It should be noted that the lack of recycling on sea-salt aerosols may affect the Iy 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: γHOI=0.12; γIONO2=0.1 (see Saiz-Lopez et al.2015 supplements). Transfer coefficients for heterogeneous reactions (see Table 3) are also calculated by free molecular transfer approximation (McFiggans et al.2000) but using reactive uptake coefficients. The reactive ice uptake and recycling of HOI and IONO2 is applied in SOCOL-AERv2-I after effective ice uptake and removal of HOI, IONO2, 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 Saiz-Lopez et al. (2015).

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 × 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 × 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 Cuevas et al. (2018) and Legrand et al. (2018) continued despite no dramatic forecast of iodine emission’s evolution being made by Iglesias-Suarez et al. (2020). 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 (Archibald et al.2020) and SST (Taylor et al.2012). 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 ×  and 1 × 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 CO2 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 4.

Table 4The experiments with SOCOL-AERv2-I.

Download Print Version | Download XLSX

3 Results of simulation

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

The total gas-phase inorganic Iy in both 2 × iodine and 1 × iodine experiments of SOCOL-AERv2-I averaged over the tropics (20 N–20 S), for the 2000–2009 period and 10 ensemble members, is presented in Fig. 2.

Figure 2Modeled vertical distribution of total organic (Iorg) and inorganic (Iy) gas-phase iodine simulated with SOCOL-AERv2-I averaged over the tropics (20 N–20 S), for the 2000–2009 period and 10 ensemble members. Red curve: Iy from the experiment 2 × iodine. Blue curve: Iy from the experiment 1 × iodine. Light red and blue curves: Iorg from 2 × iodine and 1 × iodine experiments, correspondingly. Shadings represent a standard deviation of tropical iodine (20 N–20 S).


The Iy was calculated as a sum of inorganic iodine compounds presented in the SOCOL-AERv2-I’s iodine scheme: Iy= IO + I + OIO + HOI + 2 ×  I2+ HI + INO + INO2+ 2 ×  (I2O2+ I2O3+ I2O4) + ICl+ IBr + IONO2+ CH2I. Iorg is calculated as a sum of all organic compounds in SOCOL-AERv2-I: Iorg= CH3I + 2 ×  CH2I2+ CH2ICl + CH2IBr. The standard deviation of iodine between ensemble members is found to be less than 1 %. In the lower troposphere, the Iy 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), Iy 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 Iy shows an increase until about 50 hPa and then stays constant because there is no deposition of iodine above this layer. The simulated stratospheric Iy 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 Iy; however, it becomes higher in the middle stratosphere than in Saiz-Lopez et al. (2015). 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 Iy in the stratosphere relative to the upper and middle troposphere can be also explained by the rapid tropospheric sinks and relatively slow transport of Iy 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 Iy abundance calculated with SOCOL-AER2-I does not exceed 1.05 pptv and corresponds well to the estimation given by Solomon et al. (1994), although it is slightly larger than the most recent assessment from WMO (2018) (0.8 pptv Iy).

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

Figure 4January–February averages of modeled and observed IO in the tropical troposphere [15 N–15 S] for (a) the TORERO campaign from Costa Rica (January–February 2012, 10 N–40 S, 250–285 E), and (b) the CONTRAST campaign from Guam (January–February 2014, 40 N–15 S, 115–175 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 3 presents the vertical distributions of Iy and IO at all latitudes and global maps of Iy and IO averaged for the lower stratosphere (100–70 hPa). Iy 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 O3 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) Iy mixing ratio as depicted in Fig. 3c. Their formation might have resulted from different convection patterns (weaker/stronger) over these areas. Also, the area with a higher stratospheric Iy burden is located right over the region with higher surface emissions of HOI / I2 (see Fig. 1). 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 Iy 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 O3–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 N; 20 E) in March is >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 0.65 pptv at 17 km) (Pundt et al.1998). Yet, the balloon-borne measured upper limit of IO mixing ratio is 0.2 pptv (Pundt et al.1998), which is more than 2 times lower than the IO mixing ratio simulated with SOCOL-AER2-I. SOCOL-AER2-I also captures well the Iy estimated by the box model (Pundt et al.1998) 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 N; 12 E) in March (Wittrock et al.2000), showing >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 (Volkamer et al.2015; Pan et al.2017). The results are shown in Fig. 4.

We compare the observed and measured IO only over the tropics (15 N–15 S) as extratropical measurements from both campaigns are scarcer and the comparison will be biased. Modeled IO in Fig. 4 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. 4 as blue points. The averaged values for Fig. 4 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 Iy from SOCOL-AERv2-I and CAM-chem models with the Iy derived from IO / Iy 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 Iy is estimated as 30 % of the IO / Iy ratio including errors in the calibration of in situ and remote sensing data and accounting for differences in the spatial scales (Wang et al.2015; Koenig et al.2017). IO was taken as an average of IO fields measured with AMAX-DOAS during both TORERO and CONTRAST campaigns (Volkamer et al.2015; Pan et al.2017; Koenig et al.2020). Thus, the inferred Iy is based on the measured IO and modeled IO / Iy. The result of this comparison is illustrated in Fig. 5.

Figure 5Total gas-phase Iy 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 N–40 S, 250–285 E) and CONTRAST measurements (40 N–15 S, 115–175 E). Red and green solid and dashed lines: gas-phase Iy simulated with SOCOL-AERv2-I and CAM-chem averaged for the January–February period. Light red and light green dashed lines: assumed gas-phase Iy simulated with SOCOL-AERv2-I and CAM-chem excluding approx. iodine in the particulate phase. Blue dots: gas-phase Iy modeled by the University of Colorado (CU) chemical box model. Black error bars: the uncertainty in Iy modeled by the CU box model. Blue, purple, and orange solid lines represent the empirical fit for gas-phase Iy (gas), particulate Iy (part), and total Iy (gas + part), correspondingly. Shadings represent the uncertainty in the empirical fit.


Data were plotted with respect to the H2O / O3 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 Koenig et al. (2020). To make CU box-model data clearer they were thinned out by averaging the data every 50 units of the H2O / O3 ratio between 1000 and 100 of the H2O / O3 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 Iy (gas) and Iy (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 Koenig et al.2020, supplements).

Below the H2O / O3 ratio of 70, both models capture well the Iy estimated from observations. Higher up, the red and green dashed lines represent the gas-phase Iy that is simulated with both models. Because the particulate iodine is not considered in these models, the total inorganic Iy 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 Iy if the estimated particulate iodine is excluded from the modeled Iy. To exclude the unknown particulate iodine, the modeled gas-phase Iy was subtracted from the 0.87 pptv of iodine that is empirical total Iy plus uncertainty of observations (0.1 pptv) (see the orange solid line and its uncertainty in Fig. 5). Such an approach can show what would be the approximate modeled gas-phase Iy if the iodine in the particulate phase was presented in models. It is seen that after the level of H2O / O3 ratio of ∼70, the gas-phase Iy estimated from measurements rapidly decreases and reaches the value of about 0.1 pptv at the approximate altitude of the lowermost stratosphere (H2O / O3 ratio <10). After excluding the estimated particulate iodine, the assumed modeled gas-phase Iy 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 (Koenig et al.2020). The assumption that the overestimation of modeled gas-phase Iy against observations is because of an absence of iodine in the particulate phase is reasonable as is seen in Fig. 5. 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 Iy 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 × iodine), and the control experiment neglecting iodine chemistry (0 × iodine). The contribution of iodine chemistry to present-day ozone climatology estimated by the SOCOL-AERv2-I is presented in Fig. 6.

Figure 6Modeled 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 × iodine) relative to the control run (0 × 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 %.

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: ((EXP-REF)/REF)×100, where EXP is the ozone from the experiment with 1 × 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 (Chipperfield et al.2018). 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 IOx 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 3). There is evidence that cross-cycles between bromine and chlorine are the dominant contribution to ozone reduction in the mid-latitudes and polar regions (Fernandez et al.2017; Barrera et al.2020). The active iodine in the form of IO is also an effective reaction partner for BrO+ClO as established by Read et al. (2008). Thompson et al. (2015) 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 (Sioris et al.2006).

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. 7.

Figure 7Modeled 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 × iodine) relative to the control run (0 × iodine). Blue solid line: annual mean tropopause height.


Comparing Figs. 6a and 7, 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 ∼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 Fernandez et al. (2017) and Barrera et al. (2020).

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 ClOx cycle (Zubov et al.2013) is not observed. The IOx-catalytic cycle is ineffective in the upper stratosphere similarly to BrOx 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 (Daniel et al.1999). So, the probability of the terminal reaction ClO+O is higher than those of IO+O or BrO+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. 3b). Hence, the effect of iodine chemistry on upper stratospheric ozone loss is negligibly small.

The tropospheric effect is 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 Sherwen et al. (2016a), but surface iodine emissions were a bit higher. It is also in agreement with an estimation made by Davis et al. (1996). 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 (Saiz-Lopez et al.2014).

The total ozone column (TOC) is affected by iodine mostly over high latitudes (see Fig. 6c 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 8Modeled 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 × iodine) relative to present-day iodine emissions (1 × 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 %.


In this case, the EXP was taken to be the experiment with 2 × loading of iodine whereas the REF is the one with observed (1 ×) 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. 6. 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. 6, 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 9Modeled vertical profiles of O3 reduction by iodine from organic, inorganic, and total surface emissions averaged over the tropics (20 N–20 S), for the 2000–2009 period and 10 ensemble members. Red curve: difference in O3 relative to the control run (0 × iodine), if only organic iodine source gases are considered; blue curve: a percentage difference in O3 relative to the control run (0 × iodine), if inorganic iodine source gases are only considered; green curve: a percentage difference in O3 relative to the control run (0 × 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 %.


In the lower troposphere, the iodine from inorganic emissions is responsible for ∼75 % of total iodine effect on ozone, and the contribution of iodine from organic sources is ∼25 %, as expected (Koenig et al.2020). 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., ∼60 % and ∼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 (Saiz-Lopez et al.2015). A gradual increase in Iy 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 Iy with the values determined by the CU box model based on AMAX-DOAS observations of IO showed that the model reproduces the Iy in the upper troposphere very well, while in the stratosphere, Iy 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 Iy, SOCOL-AERv2-I corresponds well with Iy 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 Hossaini et al. (2015), showing up to 10 % of the lower stratospheric ozone loss driven by iodine chemistry. It should be noted that Hossaini et al. (2015) reported only 0.15 pptv of Iy 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 Ball et al. (2018) 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 I-Br/I-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 (Daniel et al.1999). 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. 3).

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 (Volkamer et al.2021). However, such complex aerosol iodine chemistry cannot be properly simulated with current chemistry–climate models (Baker2004; Sherwen et al.2016b). The recent studies of Gómez-Martín et al. (2020), Baccarini et al. (2020), and He et al. (2021) 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 (Cuevas et al.2018) 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 I2O2, I2O3, and I2O4 might decompose to form iodine oxoacids that grow further to become CCNs (McFiggans et al.2004; Burkholder et al.2004; Saunders et al.2010). However, the formation of iodine oxoacids is still an open question that needs further study (He et al.2021). 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 (Lewis et al.2020). 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 (Tham et al.2021). Also, we did not include CF3I 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 CH3I (Zhang et al.2020), and will mostly affect the tropical and northern mid-latitude tropospheric ozone because of higher concentrations of pollutants (Youn et al.2010). All organic iodine emissions in our scheme are prescribed. However, they can be interactively calculated in the model utilizing ocean biogenic sources (Ordóñez et al.2012). 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 / I2 fluxes has some uncertainty in the models and observations (Chance et al.2014; Sherwen et al.2019). However, Shaw and Carpenter (2013) argued that the iodide field does not necessarily solve any uncertainty in the HOI / I2 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 (Carpenter et al.2021). 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 (Inamdar et al.2020). Recent investigations of the chemical basis of the HOI / I2 source can help improve the generalization of empirical source functions (Moreno et al.2020).

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 (Iglesias-Suarez et al.2020). Tropospheric ozone content in SOCOL-AERv2-I is overestimated compared to other models (Revell et al.2018), which affects the sea-surface deposition of O3 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 (Bureau et al.2000; Aiuppa et al.2005; Balcone-Boissard et al.2010; Cadoux et al.2015). However, there is no solid evidence that volcanoes can inject a sufficient amount of iodine into the atmosphere (Schönhardt et al.2017), 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 (Cuevas et al.2018; Legrand et al.2018; Koenig et al.2020; Iglesias-Suarez et al.2020; Carpenter et al.2021). 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 (Klobas et al.2021).

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: (Karagodin-Doyennel et al.2021) 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 (, last access: 11 September 2021). The SOCOL-AERv2-I simulation data can be accessed here: (Karagodin-Doyennel2021). CU-AMAX-DOAS CONTRAST IO data are available at (last access: 11 September 2021;,Volkamer et al.2020). CU-AMAX-DOAS TORERO IO data are available at (Volkamer and Dix2017). The CU box model data are available here: (Volkamer and Koenig2021).

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.


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


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.


Aiuppa, A., Federico, C., Franco, A., Giudice, G., Gurrieri, S., Inguaggiato, S., Liuzzo, M., McGonigle, A. J. S., and Valenza, M.: Emission of bromine and iodine from Mount Etna volcano, Geochem. Geophy. Geosy., 6, Q08008,, 2005. a, b

Akagi, S. K., Yokelson, R. J., Wiedinmyer, C., Alvarado, M. J., Reid, J. S., Karl, T., Crounse, J. D., and Wennberg, P. O.: Emission factors for open and domestic biomass burning for use in atmospheric models, Atmos. Chem. Phys., 11, 4039–4072,, 2011. a

Archibald, A., Turnock, S., Griffiths, P., Cox, T., Derwent, R. G., Knote, C., and Shin, M.: On the changes in surface ozone over the twenty-first century: sensitivity to changes in surface temperature and chemical mechanisms, Philos. T. Roy. Soc. A, 378, 20190329,, 2020. a, b

Aschmann, J. and Sinnhuber, B.-M.: Contribution of very short-lived substances to stratospheric bromine loading: uncertainties and constraints, Atmos. Chem. Phys., 13, 1203–1219,, 2013. a

Baccarini, A., Karlsson, L., Dommen, J., Duplessis, P., Vüllers, J., Brooks, I. M., Saiz-Lopez, A., Salter, M., Tjernström, M., Baltensperger, U., Zieger, P., and Schmale, J.: Frequent new particle formation over the high Arctic pack ice by enhanced iodine emissions, Nat. Commun., 11, 4924,, 2020. a

Badia, A., Reeves, C. E., Baker, A. R., Saiz-Lopez, A., Volkamer, R., Koenig, T. K., Apel, E. C., Hornbrook, R. S., Carpenter, L. J., Andrews, S. J., Sherwen, T., and von Glasow, R.: Importance of reactive halogens in the tropical marine atmosphere: a regional modelling study using WRF-Chem, Atmos. Chem. Phys., 19, 3161–3189,, 2019. a, b, c, d, e, f, g, h, i, j, k, l, m, n, o

Baidar, S., Oetjen, H., Coburn, S., Dix, B., Ortega, I., Sinreich, R., and Volkamer, R.: The CU Airborne MAX-DOAS instrument: vertical profiling of aerosol extinction and trace gases, Atmos. Meas. Tech., 6, 719–739,, 2013. a

Baker, A. R.: Inorganic iodine speciation in tropical Atlantic aerosol, Geophys. Res. Lett., 31, L23S02,, 2004. a

Balcone-Boissard, H., Villemant, B., and Boudon, G.: Behavior of halogens during the degassing of felsic magmas, Geochem. Geophy. Geosy., 11, Q09005,, 2010. a

Ball, W. T., Alsing, J., Mortlock, D. J., Staehelin, J., Haigh, J. D., Peter, T., Tummon, F., Stübi, R., Stenke, A., Anderson, J., Bourassa, A., Davis, S. M., Degenstein, D., Frith, S., Froidevaux, L., Roth, C., Sofieva, V., Wang, R., Wild, J., Yu, P., Ziemke, J. R., and Rozanov, E. V.: Evidence for a continuous decline in lower stratospheric ozone offsetting ozone layer recovery, Atmos. Chem. Phys., 18, 1379–1394,, 2018. a

Barrera, J. A., Fernandez, R. P., Iglesias-Suarez, F., Cuevas, C. A., Lamarque, J.-F., and Saiz-Lopez, A.: Seasonal impact of biogenic very short-lived bromocarbons on lowermost stratospheric ozone between 60 N and 60 S during the 21st century, Atmos. Chem. Phys., 20, 8083–8102,, 2020. a, b

Bell, N., Hsu, L., Jacob, D. J., Schultz, M. G., Blake, D. R., Butler, J. H., King, D. B., Lobert, J. M., and Maier-Reimer, E.: Methyl iodide: Atmospheric budget and use as a tracer of marine convection in global models, J. Geophys. Res.-Atmos., 107, 4340,, 2002. a

Bloss, W. J., Camredon, M., Lee, J. D., Heard, D. E., Plane, J. M. C., Saiz-Lopez, A., Bauguitte, S. J.-B., Salmon, R. A., and Jones, A. E.: Coupling of HOx, NOx and halogen chemistry in the antarctic boundary layer, Atmos. Chem. Phys., 10, 10187–10209,, 2010. a

Bureau, H., Keppler, H., and Métrich, N.: Volcanic degassing of bromine and iodine: experimental fluid/melt partitioning data and applications to stratospheric chemistry, Eart Planet Sc. Lett., 183, 51–60,, 2000. a, b

Burkholder, J. B., Curtius, J., Ravishankara, A. R., and Lovejoy, E. R.: Laboratory studies of the homogeneous nucleation of iodine oxides, Atmos. Chem. Phys., 4, 19–34,, 2004. a

Burkholder, J. B., Sander, S. P., Abbatt, J., Barker, J. R., Huie, R. E., Kolb, C. E., Kurylo, M. J., Orkin, V. L., Wilmouth, D. M., and Wine, P. H.: Chemical Kinetics and Photochemical Data for Use in Atmospheric Studies, Jet Propulsion Laboratory, available at: (last access: 20 October 2021), 2015. a, b, c, d, e, f, g, h, i, j, k, l, m, n, o, p, q, r, s, t, u, v, w, x, y, z, aa, ab, ac, ad, ae, af, ag, ah

Cadoux, A., Scaillet, B., Bekki, S., Oppenheimer, C., and Druitt, T. H.: Stratospheric Ozone destruction by the Bronze-Age Minoan eruption (Santorini Volcano, Greece), Sci. Rep., 5, 12243,, 2015. a

Carpenter, L. J., MacDonald, S. M., Shaw, M. D., Kumar, R., Saunders, R. W., Parthipan, R., Wilson, J., and Plane, J. M. C.: Atmospheric iodine levels influenced by sea surface emissions of inorganic iodine, Nat. Geosci., 6, 108–111,, 2013. a, b, c, d, e, f

Carpenter, L. J., Chance, R. J., Sherwen, T., Adams, T. J., Ball, S. M., Evans, M. J., Hepach, H., Hollis, L. D. J., Hughes, C., Jickells, T. D.,Mahajan, A., Stevens, D. P., Tinel, L., and Wadley, M. R.: Marine iodine emissions in a changing world, Proc. Roy. Soc. Lond. A, 477, 20200824,, 2021. a, b

Chameides, W. L. and Davis, D. D.: Iodine: Its possible role in tropospheric photochemistry, J. Geophys. Res.-Oceans, 85, 7383–7398,, 1980. a

Chance, R., Baker, A. R., Carpenter, L., and Jickells, T. D.: The distribution of iodide at the sea surface, Environ. Sci.: Processes Impacts, 16, 1841–1859,, 2014. a, b, c

Chipperfield, M. P., Dhomse, S., Hossaini, R., Feng, W., Santee, M. L., Weber, M., Burrows, J. P., Wild, J. D., Loyola, D., and Coldewey-Egbers, M.: On the Cause of Recent Variations in Lower Stratospheric Ozone, Geophys. Res. Lett., 45, 5718–5726,, 2018. a

Cuevas, C., Maffezzoli, N., and Corella, J. e. a.: Rapid increase in atmospheric iodine levels in the North Atlantic since the mid-20th century, Nat. Commun., 9, 1452,, 2018. a, b, c, d

Daniel, J. S., Solomon, S., Portmann, R. W., and Garcia, R. R.: Stratospheric ozone destruction: The importance of bromine relative to chlorine, J. Geophys. Res.-Atmos., 104, 23871–23880,, 1999. a, b

Davis, D., Crawford, J., Liu, S., McKeen, S., Band y, A., Thornton, D., Rowland, F., and Blake, D.: Potential impact of iodine on tropospheric levels of ozone and other critical oxidants, J. Geophys. Res.-Atmos., 101, 2135–2147,, 1996. a, b, c, d, e

Dix, B., Baidar, S., Bresch, J. F., Hall, S. R., Schmidt, K. S., Wang, S., and Volkamer, R.: Detection of iodine monoxide in the tropical free troposphere, P. Natl. Acad. Sci. USA, 110, 2035–2040,, 2013. a

Egorova, T., Rozanov, E., Zubov, V., and Karol, I.: Model for investigating ozone trends (MEZON), IZV Atmos. Ocean Phy+, 39, 277–292, 2003. a

Egorova, T., Rozanov, E., Gröbner, J., Hauser, M., and Schmutz, W.: Montreal Protocol Benefits simulated with CCM SOCOL, Atmos. Chem. Phys., 13, 3811–3823,, 2013. a

Feinberg, A., Sukhodolov, T., Luo, B.-P., Rozanov, E., Winkel, L. H. E., Peter, T., and Stenke, A.: Improved tropospheric and stratospheric sulfur cycle in the aerosol–chemistry–climate model SOCOL-AERv2, Geosci. Model Dev., 12, 3863–3887,, 2019. a, b

Fernandez, R. P., Salawitch, R. J., Kinnison, D. E., Lamarque, J.-F., and Saiz-Lopez, A.: Bromine partitioning in the tropical tropopause layer: implications for stratospheric injection, Atmos. Chem. Phys., 14, 13391–13410,, 2014. a

Fernandez, R. P., Kinnison, D. E., Lamarque, J.-F., Tilmes, S., and Saiz-Lopez, A.: Impact of biogenic very short-lived bromine on the Antarctic ozone hole during the 21st century, Atmos. Chem. Phys., 17, 1673–1688,, 2017. a, b

Fuge, R. and Johnson, C. C.: Iodine and human health, the role of environmental geochemistry and diet, a review, Appl. Geochem., 63, 282–302,, 2015. a

Gómez-Martín, J. C., Spietz, P., and Burrows, J. P.: Spectroscopic studies of the I2/O3 photochemistry: Part 1: Determination of the absolute absorption cross sections of iodine oxides of atmospheric relevance, J. Photoch. Photobio. A, 176, 15–38,, 2005. a

Gómez-Martín, J. C., Lewis, T. R., Blitz, M. A., Plane, J. M. C., Kumar, M., Francisco, J. S., and Saiz-Lopez, A.: A gas-to-particle conversion mechanism helps to explain atmospheric particle formation through clustering of iodine oxides, Nat. Commun., 11, 4521,, 2020. a

Gravestock, T. J., Blitz, W. J., and Heard, D. E.: A multidimensional study of the reaction CH2I + O2: products and atmospheric implications, Chemphyschem, 11, 3928–3941,, 2010. a

Gutmann, A., Bobrowski, N., Roberts, T. J., Rüdiger, J., and Hoffmann, T.: Advances in bromine speciation in volcanic plumes, Front. Earth Sci., 6, 213,, 2018. a

He, X.-C., Tham, Y. J., Dada, L., Wang, M., Finkenzeller, H., Stolzenburg, D., Iyer, S., Simon, M., Kürten, A., Shen, J., Rörup, B., Rissanen, M., Schobesberger, S., Baalbaki, R., Wang, D. S., Koenig, T. K., Jokinen, T., Sarnela, N., Beck, L. J., Almeida, J., Amanatidis, S., Amorim, A., Ataei, F., Baccarini, A., Bertozzi, B., Bianchi, F., Brilke, S., Caudillo, L., Chen, D., Chiu, R., Chu, B., Dias, A., Ding, A., Dommen, J., Duplissy, J., El Haddad, I., Gonzalez Carracedo, L., Granzin, M., Hansel, A., Heinritzi, M., Hofbauer, V., Junninen, H., Kangasluoma, J., Kemppainen, D., Kim, C., Kong, W., Krechmer, J. E., Kvashin, A., Laitinen, T., Lamkaddam, H., Lee, C. P., Lehtipalo, K., Leiminger, M., Li, Z., Makhmutov, V., Manninen, H. E., Marie, G., Marten, R., Mathot, S., Mauldin, R. L., Mentler, B., Möhler, O., Müller, T., Nie, W., Onnela, A., Petäjä, T., Pfeifer, J., Philippov, M., Ranjithkumar, A., Saiz-Lopez, A., Salma, I., Scholz, W., Schuchmann, S., Schulze, B., Steiner, G., Stozhkov, Y., Tauber, C., Tomé, A., Thakur, R. C., Väisänen, O., Vazquez-Pufleau, M., Wagner, A. C., Wang, Y., Weber, S. K., Winkler, P. M., Wu, Y., Xiao, M., Yan, C., Ye, Q., Ylisirniö, A., Zauner-Wieczorek, M., Zha, Q., Zhou, P., Flagan, R. C., Curtius, J., Baltensperger, U., Kulmala, M., Kerminen, V.-M., Kurtén, T., Donahue, N. M., Volkamer, R., Kirkby, J., Worsnop, D. R., and Sipilä, M.: Role of iodine oxoacids in atmospheric aerosol nucleation, Science, 371, 589–595,, 2021. a, b

Heymsfield, A., Winker, D., Avery, M., Vaughan, M., Diskin, G., Deng, M., Mitev, V., and Matthey, R.: Relationships between Ice Water Content and Volume Extinction Coefficient from In Situ Observations for Temperatures from 0 to 86 C: Implications for Spaceborne Lidar Retrievals, J. Appl. Meteorol. Clim., 53, 479–505,, 2014. a

Hilal, S. H., Ayyampalayam, S. N., and Carreira, L. A.: Air-Liquid Partition Coefficient for a Diverse Set of Organic Compounds: Henry's Law Constant in Water and Hexadecane, Environ. Sci. Technol., 42, 9231–9236,, 2008. a, b

Holmes, C. D., Bertram, T. H., Confer, K. L., Graham, K. A., Ronan, A. C., Wirks, C. K., and Shah, V.: The Role of Clouds in the Tropospheric NOx Cycle: A New Modeling Approach for Cloud Chemistry and Its Global Implications, Geophys. Res. Lett., 46, 4980–4990,, 2019. a

Hossaini, R., Chipperfield, M. P., Montzka, S. A., Rap, A., Dhomse, S., and Feng, W.: Efficiency of short-lived halogens at influencing climate through depletion of stratospheric ozone, Nat. Geosci., 8, 186–190,, 2015. a, b, c, d

Iglesias-Suarez, F., Badia, A., Fernandez, R. P., Cuevas, C. A., Kinnison, D. E., Tilmes, S., Lamarque, J.-F., Long, M. C., Hossaini, R., and Saiz-Lopez, A.: Natural halogens buffer tropospheric ozone in a changing climate, Nat. Clim. Change, 10, 147–154,, 2020. a, b, c, d

Inamdar, S., Tinel, L., Chance, R., Carpenter, L. J., Sabu, P., Chacko, R., Tripathy, S. C., Kerkar, A. U., Sinha, A. K., Bhaskar, P. V., Sarkar, A., Roy, R., Sherwen, T., Cuevas, C., Saiz-Lopez, A., Ram, K., and Mahajan, A. S.: Estimation of reactive inorganic iodine fluxes in the Indian and Southern Ocean marine boundary layer, Atmos. Chem. Phys., 20, 12093–12114,, 2020. a

Jones, C. E., Hornsby, K. E., Sommariva, R., Dunk, R. M., von Glasow, R., McFiggans, G., and Carpenter, L. J.: Quantifying the contribution of marine organic gases to atmospheric iodine, Geophys. Res. Lett., 37, L18804,, 2010. a

Karagodin-Doyennel, A.: The iodine chemistry and ozone data simulated with SOCOL AERv2-iodine, Zenodo [data set],, 2021. a

Karagodin-Doyennel, A., Rozanov, E., and Sukhodolov, T.: Aerosol-Chemistry-Climate Model SOCOL-AERv2-iodine code (Version 1.0), Zenodo [code],, 2021. a

Klobas, J. E., Hansen, J., Weisenstein, D. K., Kennedy, R. P., and Wilmouth, D. M.: Sensitivity of Iodine-Mediated Stratospheric Ozone Loss Chemistry to Future Chemistry-Climate Scenarios, Frontiers, 9, 268,, 2021. a, b

Koenig, T. K., Volkamer, R., Baidar, S., Dix, B., Wang, S., Anderson, D. C., Salawitch, R. J., Wales, P. A., Cuevas, C. A., Fernandez, R. P., Saiz-Lopez, A., Evans, M. J., Sherwen, T., Jacob, D. J., Schmidt, J., Kinnison, D., Lamarque, J.-F., Apel, E. C., Bresch, J. C., Campos, T., Flocke, F. M., Hall, S. R., Honomichl, S. B., Hornbrook, R., Jensen, J. B., Lueb, R., Montzka, D. D., Pan, L. L., Reeves, J. M., Schauffler, S. M., Ullmann, K., Weinheimer, A. J., Atlas, E. L., Donets, V., Navarro, M. A., Riemer, D., Blake, N. J., Chen, D., Huey, L. G., Tanner, D. J., Hanisco, T. F., and Wolfe, G. M.: BrO and inferred Bry profiles over the western Pacific: relevance of inorganic bromine sources and a Bry minimum in the aged tropical tropopause layer, Atmos. Chem. Phys., 17, 15245–15270,, 2017. a

Koenig, T. K., Baidar, S., Campuzano-Jost, P., Cuevas, C. A., Dix, B., Fernandez, R. P., Guo, H., Hall, S. R., Kinnison, D., Nault, B. A., Ullmann, K., Jimenez, J. L., Saiz-Lopez, A., and Volkamer, R.: Quantitative detection of iodine in the stratosphere, P. Natl. Acad. Sci. USA, 117, 1860–1866,, 2020. a, b, c, d, e, f, g, h

Lee-Taylor, J. and Redeker, K. R.: Reevaluation of global emissions from rice paddies of methyl iodide and other species, Geophys. Res. Lett., 32, L15801,, 2005. a

Legrand, M., McConnell, J. R., Preunkert, S., Arienzo, M., Chellman, N., Gleason, K., Sherwen, T., Evans, M. J., and Carpenter, L. J.: Alpine ice evidence of a three-fold increase in atmospheric iodine deposition since 1950 in Europe due to increasing oceanic emissions, P. Natl. Acad. Sci. USA, 115, 12136–12141,, 2018. a, b, c

Lewis, T. R., Gómez Martín, J. C., Blitz, M. A., Cuevas, C. A., Plane, J. M. C., and Saiz-Lopez, A.: Determination of the absorption cross sections of higher-order iodine oxides at 355 and 532 nm, Atmos. Chem. Phys., 20, 10865–10887,, 2020. a

Lin, S.-J. and Rood, R. B.: Multidimensional Flux-Form Semi-Lagrangian Transport Schemes, Mon. Weather Rev., 124, 2046,<2046:MFFSLT>2.0.CO;2, 1996. a

Lovelock, J. E. and Maggs, R. J.: Halogenated Hydrocarbons in and over the Atlantic, Nature, 241, 194–196,, 1973. a

MacDonald, S. M., Gómez Martín, J. C., Chance, R., Warriner, S., Saiz-Lopez, A., Carpenter, L. J., and Plane, J. M. C.: A laboratory characterisation of inorganic iodine emissions from the sea surface: dependence on oceanic variables and parameterisation for global modelling, Atmos. Chem. Phys., 14, 5841–5852,, 2014. a, b, c, d

Manley, S. L., Wang, N.-Y., Walser, M. L., and Cicerone, R. J.: Methyl halide emissions from greenhouse-grown mangroves, Geophys. Res. Lett., 34, L01806,, 2007. a

Manzini, E., Giorgetta, M. A., Esch, M., Kornblueh, L., and Roeckner, E.: The Influence of Sea Surface Temperatures on the Northern Winter Stratosphere: Ensemble Simulations with the MAECHAM5 Model, J. Climate, 19, 3863,, 2006. a

McClain, C. R., Feldman, G. C., and Hooker, S. B.: An overview of the SeaWiFS project and strategies for producing a climate research quality global ocean bio-optical time series, Deep-Sea Res. Pt. II, 51, 5–42,, 2004. a

McFiggans, G., Plane, J. M. C., Allan, B. J., Carpenter, L. J., Coe, H., and O'Dowd, C.: A modeling study of iodine chemistry in the marine boundary layer, J. Geophys. Res.-Atmos., 105, 14371–14385,, 2000. a, b

McFiggans, G., Coe, H., Burgess, R., Allan, J., Cubison, M., Alfarra, M. R., Saunders, R., Saiz-Lopez, A., Plane, J. M. C., Wevill, D., Carpenter, L., Rickard, A. R., and Monks, P. S.: Direct evidence for coastal iodine particles from Laminaria macroalgae – linkage to emissions of molecular iodine, Atmos. Chem. Phys., 4, 701–713,, 2004. a

Moreno, C., Baeza-Romero, M.-T., Sanz, M., Gálvez, Ó., López Arza, V., Ianni, J. C., and Espíldora, E.: Iodide conversion to iodate in aqueous and solid aerosols exposed to ozone, Phys. Chem. Chem. Phys., 22, 5625–5637,, 2020. a

Newman, P. A., Oman, L. D., Douglass, A. R., Fleming, E. L., Frith, S. M., Hurwitz, M. M., Kawa, S. R., Jackman, C. H., Krotkov, N. A., Nash, E. R., Nielsen, J. E., Pawson, S., Stolarski, R. S., and Velders, G. J. M.: What would have happened to the ozone layer if chlorofluorocarbons (CFCs) had not been regulated?, Atmos. Chem. Phys., 9, 2113–2128,, 2009. a

Ordóñez, C., Lamarque, J.-F., Tilmes, S., Kinnison, D. E., Atlas, E. L., Blake, D. R., Sousa Santos, G., Brasseur, G., and Saiz-Lopez, A.: Bromine and iodine chemistry in a global chemistry-climate model: description and evaluation of very short-lived oceanic sources, Atmos. Chem. Phys., 12, 1423–1447,, 2012. a, b, c, d, e, f, g, h, i, j, k, l, m, n, o, p, q, r

Pan, L. L., Atlas, E. L., Salawitch, R. J., Honomichl, S. B., Bresch, J. F., Randel, W. J., Apel, E. C., Hornbrook, R. S., Weinheimer, A. J., Anderson, D. C., Andrews, S. J., Baidar, S., Beaton, S. P., Campos, T. L., Carpenter, L. J., Chen, D., Dix, B., Donets, V., Hall, S. R., Hanisco, T. F., Homeyer, C. R., Huey, L. G., Jensen, J. B., Kaser, L., Kinnison, D. E., Koenig, T. K., Lamarque, J. F., Liu, C., Luo, J., Luo, Z. J., Montzka, D. D., Nicely, J. M., Pierce, R. B., Riemer, D. D., Robinson, T., Romashkin, P., Saiz-Lopez, A., Schauffler, S., Shieh, O., Stell, M. H., Ullmann, K., Vaughan, G., Volkamer, R., and Wolfe, G.: The Convective Transport of Active Species in the Tropics (CONTRAST) Experiment, B. Am. Meteorol. Soc., 98, 106–128,, 2017. a, b, c

Prados-Roman, C., Cuevas, C. A., Fernandez, R. P., Kinnison, D. E., Lamarque, J.-F., and Saiz-Lopez, A.: A negative feedback between anthropogenic ozone pollution and enhanced ocean emissions of iodine, Atmos. Chem. Phys., 15, 2215–2224,, 2015a. a

Prados-Roman, C., Cuevas, C. A., Hay, T., Fernandez, R. P., Mahajan, A. S., Royer, S.-J., Galí, M., Simó, R., Dachs, J., Großmann, K., Kinnison, D. E., Lamarque, J.-F., and Saiz-Lopez, A.: Iodine oxide in the global marine boundary layer, Atmos. Chem. Phys., 15, 583–593,, 2015b. a, b

Pundt, I., Pommereau, J.-P., Phillips, C., and Lateltin, E.: Upper Limit of Iodine Oxide in the Lower Stratosphere, J. Atmos. Chem., 30, 173–185,, 1998. a, b, c, d

Rayner, N. A., Parker, D. E., Horton, E. B., Folland , C. K., Alexander, L. V., Rowell, D. P., Kent, E. C., and Kaplan, A.: Global analyses of sea surface temperature, sea ice, and night marine air temperature since the late nineteenth century, J. Geophys. Res.-Atmos., 108, 4407,, 2003. a

Read, K. A., Mahajan, A. S., Carpenter, L. J., Evans, M. J., Faria, B. V. E., Heard, D. E., Hopkins, J. R., Lee, J. D., Moller, S. J., Lewis, A. C., Mendes, L., McQuaid, J. B., Oetjen, H., Saiz-Lopez, A., Pilling, M. J., and Plane, J. M. C.: Extensive halogen-mediated ozone destruction over the tropical Atlantic Ocean, Nature, 453, 1232–1235,, 2008. a

Revell, L. E., Tummon, F., Stenke, A., Sukhodolov, T., Coulon, A., Rozanov, E., Garny, H., Grewe, V., and Peter, T.: Drivers of the tropospheric ozone budget throughout the 21st century under the medium-high climate scenario RCP 6.0, Atmos. Chem. Phys., 15, 5887–5902,, 2015. a

Revell, L. E., Stenke, A., Tummon, F., Feinberg, A., Rozanov, E., Peter, T., Abraham, N. L., Akiyoshi, H., Archibald, A. T., Butchart, N., Deushi, M., Jöckel, P., Kinnison, D., Michou, M., Morgenstern, O., O'Connor, F. M., Oman, L. D., Pitari, G., Plummer, D. A., Schofield, R., Stone, K., Tilmes, S., Visioni, D., Yamashita, Y., and Zeng, G.: Tropospheric ozone in CCMI models and Gaussian process emulation to understand biases in the SOCOLv3 chemistry–climate model, Atmos. Chem. Phys., 18, 16155–16172,, 2018. a, b

Saiz-Lopez, A. and Plane, J. M. C.: Novel iodine chemistry in the marine boundary layer, Geophys. Res. Lett., 31, L04112,, 2004. a

Saiz-Lopez, A., Plane, J. M. C., Baker, A. R., Carpenter, L. J., Von Glasow, R., Gómez Martín, J. C., McFiggans, G., and Saunders, R. W.: Atmospheric Chemistry of Iodine, Chem. Rev., 112, 1773–1804,, 2012. a, b

Saiz-Lopez, A., Fernandez, R. P., Ordóñez, C., Kinnison, D. E., Gómez Martín, J. C., Lamarque, J.-F., and Tilmes, S.: Iodine chemistry in the troposphere and its effect on ozone, Atmos. Chem. Phys., 14, 13119–13143,, 2014. a, b, c, d, e, f, g, h, i, j, k, l, m, n, o, p, q, r, s, t

Saiz-Lopez, A., Baidar, S., Cuevas, C. A., Koenig, T. K., Fernandez, R. P., Dix, B., Kinnison, D. E., Lamarque, J. F., Rodriguez-Lloveras, X., Campos, T. L., and Volkamer, R.: Injection of iodine to the stratosphere, Geophys. Res. Lett., 42, 6852–6859,, 2015. a, b, c, d, e, f, g, h, i, j, k, l, m

Saunders, R. W., Kumar, R., Martín, J. C. G., Mahajan, A. S., Murray, B. J., and Plane, J. M. C.: Studies of the Formation and Growth of Aerosol from Molecular Iodine Precursor, Z. Phys. Chem., 224, 1095–1117,, 2010. a

Schmidt, J. A., Jacob, D. J., Horowitz, H. M., Hu, L., Sherwen, T., Evans, M. J., Liang, Q., Suleiman, R. M., Oram, D. E., Le Breton, M., Percival, C. J., Wang, S., Dix, B., and Volkamer, R.: Modeling the observed tropospheric BrO background: Importance of multiphase chemistry and implications for ozone, OH, and mercury, J. Geophys. Res.-Atmos., 121, 11819–11835,, 2016. a

Schönhardt, A., Richter, A., Wittrock, F., Kirk, H., Oetjen, H., Roscoe, H. K., and Burrows, J. P.: Observations of iodine monoxide columns from satellite, Atmos. Chem. Phys., 8, 637–653,, 2008. a

Schönhardt, A., Richter, A., Theys, N., and Burrows, J. P.: Space-based observation of volcanic iodine monoxide, Atmos. Chem. Phys., 17, 4857–4870,, 2017. a

Shaw, M. D. and Carpenter, L. J.: Modification of Ozone Deposition and I2Emissions at the Air-Aqueous Interface by Dissolved Organic Carbon of Marine Origin, Environ. Sci. Technol., 47, 10947–10954,, 2013. a

Sheng, J.-X., Weisenstein, D. K., Luo, B.-P., Rozanov, E., Stenke, A., Anet, J., Bingemer, H., and Peter, T.: Global atmospheric sulfur budget under volcanically quiescent conditions: Aerosol-chemistry-climate model predictions and validation, J. Geophys. Res.-Atmos., 120, 256–276,, 2015. a

Sherwen, T., Evans, M. J., Carpenter, L. J., Andrews, S. J., Lidster, R. T., Dix, B., Koenig, T. K., Sinreich, R., Ortega, I., Volkamer, R., Saiz-Lopez, A., Prados-Roman, C., Mahajan, A. S., and Ordóñez, C.: Iodine's impact on tropospheric oxidants: a global model study in GEOS-Chem, Atmos. Chem. Phys., 16, 1161–1186,, 2016a. a, b, c, d

Sherwen, T. M., Evans, M. J., Spracklen, D. V., Carpenter, L. J., Chance, R., Baker, A. R., Schmidt, J. A., and Breider, T. J.: Global modeling of tropospheric iodine aerosol, Geophys. Res. Lett., 43, 10012–10019,, 2016b. a

Sherwen, T., Schmidt, J. A., Evans, M. J., Carpenter, L. J., Großmann, K., Eastham, S. D., Jacob, D. J., Dix, B., Koenig, T. K., Sinreich, R., Ortega, I., Volkamer, R., Saiz-Lopez, A., Prados-Roman, C., Mahajan, A. S., and Ordóñez, C.: Global impacts of tropospheric halogens (Cl, Br, I) on oxidants and composition in GEOS-Chem, Atmos. Chem. Phys., 16, 12239–12271,, 2016c. a, b, c, d

Sherwen, T., Chance, R. J., Tinel, L., Ellis, D., Evans, M. J., and Carpenter, L. J.: A machine-learning-based global sea-surface iodide distribution, Earth Syst. Sci. Data, 11, 1239–1262,, 2019. a, b

Sioris, C. E., Kovalenko, L. J., McLinden, C. A., Salawitch, R. J., van Roozendael, M., Goutail, F., Dorf, M., Pfeilsticker, K., Chance, K., von Savigny, C., Liu, X., Kurosu, T. P., Pommereau, J. P., BöSch, H., and Frerick, J.: Latitudinal and vertical distribution of bromine monoxide in the lower stratosphere from Scanning Imaging Absorption Spectrometer for Atmospheric Chartography limb scattering measurements, J. Geophys. Res.-Atmos., 111, D14301,, 2006. a

Solomon, S.: Stratospheric ozone depletion: A review of concepts and history, Rev. Geophys., 37, 275–316,, 1999. a

Solomon, S., Garcia, R. R., and Ravishankara, A. R.: On the role of iodine in ozone depletion, J. Geophys. Res.-Atmos., 99, 20491–20499,, 1994. a, b, c, d

Stenke, A., Schraner, M., Rozanov, E., Egorova, T., Luo, B., and Peter, T.: The SOCOL version 3.0 chemistry–climate model: description, evaluation, and implications from an advanced transport algorithm, Geosci. Model Dev., 6, 1407–1427,, 2013. a

Taylor, K. E., Stouffer, R. J., and Meehl, G. A.: An Overview of CMIP5 and the Experiment Design, B. Am. Meteorol. Soc., 93, 485–498,, 2012. a

Tham, Y. J., He, X.-C., Li, Q., Cuevas, C. A., Shen, J., Kalliokoski, J., Yan, C., Iyer, S., Lehmusjärvi, T., Jang, S., Thakur, R. C., Beck, L., Kemppainen, D., Olin, M., Sarnela, N., Mikkilä, J., Hakala, J., Marbouti, M., Yao, L., Li, H., Huang, W., Wang, Y., Wimmer, D., Zha, Q., Virkanen, J., Spain, T. G., O'Doherty, S., Jokinen, T., Bianchi, F., Petäjä, T., Worsnop, D. R., Mauldin, R. L., Ovadnevaite, J., Ceburnis, D., Maier, N. M., Kulmala, M., O'Dowd, C., Dal Maso, M., Saiz-Lopez, A., and Sipilä, M.: Direct field evidence of autocatalytic iodine release from atmospheric aerosol, P. Natl. Acad. Sci. USA, 118, e2009951118,, 2021. a, b

Thompson, C. R., Shepson, P. B., Liao, J., Huey, L. G., Apel, E. C., Cantrell, C. A., Flocke, F., Orlando, J., Fried, A., Hall, S. R., Hornbrook, R. S., Knapp, D. J., Mauldin III, R. L., Montzka, D. D., Sive, B. C., Ullmann, K., Weibring, P., and Weinheimer, A.: Interactions of bromine, chlorine, and iodine photochemistry during ozone depletions in Barrow, Alaska, Atmos. Chem. Phys., 15, 9651–9679,, 2015. a

Vogt, R., Crutzen, P. J., and Sander, R.: A mechanism for halogen release from sea-salt aerosol in the remote marine boundary layer, Nature, 383, 327–330,, 1996. a, b

Vogt, R., Sander, R., von Glasow, R., and Crutzen, P. J.: Iodine Chemistry and its Role in Halogen Activation and Ozone Loss in the Marine Boundary Layer: A Model Study, J. Atmos. Chem., 32, 375–395,, 1999. a, b, c

Volkamer, R., Dix, B.: GV AMAX-DOAS Data. Version 5.0, UCAR/NCAR – Earth Observing Laboratory [data set], available at: (last access: 20 October 2021), 2017. a

Volkamer, R. and Koenig, K. T.: The CU box-model data (Version 1.0). Zenodo [data set],, 2021. a

Volkamer, R., Baidar, S., Campos, T. L., Coburn, S., DiGangi, J. P., Dix, B., Eloranta, E. W., Koenig, T. K., Morley, B., Ortega, I., Pierce, B. R., Reeves, M., Sinreich, R., Wang, S., Zondlo, M. A., and Romashkin, P. A.: Aircraft measurements of BrO, IO, glyoxal, NO2, H2O, O2–O2 and aerosol extinction profiles in the tropics: comparison with aircraft-/ship-based in situ and lidar measurements, Atmos. Meas. Tech., 8, 2121–2148,, 2015. a, b, c, d

Volkamer, R., Koenig, T., Baidar, S., and Dix, B.: Airborne Multi-AXis Differential Optical Absorption Spectroscopy (AMAX-DOAS) Data, Version 2.0, UCAR/NCAR – Earth Observing Laboratory [data set],, 2020. a

Volkamer, R., Koenig, T., Apel, E., Bresch, J., Cuevas, C., Dix, B., Eloranta, E., Fernandez, R., Hall, S., Hornbrook, R., Pierce, B., Reeves, M., Saiz-Lopez, A., Spuhler, S., and Ullman, K.: Mini ozone holes due to dust release of iodine, EGU General Assembly 2021, online, 19–30 April 2021, EGU21-13841,, 2021. a

Wagman, D. D., Evans, W. H., Parker, V. B., Schumm, R. H., Halow, I., Bailey, S. M., Churney, K. L., and Nuttall, R. L.: Erratum: The NBS tables of chemical thermodynamic properties. Selected values for inorganic and C1 and C2 organic substances in SI units [J. Phys. Chem. Ref. Data 11, Suppl. 2 (1982)], J. Phys. Chem. Ref. Data, 18, 1807–1812,, 1989. a, b

Wang, S., Schmidt, J. A., Baidar, S., Coburn, S., Dix, B., Koenig, T. K., Apel, E., Bowdalo, D., Campos, T. L., Eloranta, E., Evans, M. J., DiGangi, J. P., Zondlo, M. A., Gao, R.-S., Haggerty, J. A., Hall, S. R., Hornbrook, R. S., Jacob, D., Morley, B., Pierce, B., Reeves, M., Romashkin, P., ter Schure, A., and Volkamer, R.: Active and widespread halogen chemistry in the tropical and subtropical free troposphere, P. Natl. Acad. Sci. USA, 112, 9281–9286,, 2015.  a

Weisenstein, D. K., Yue, G. K., Ko, M. K. W., Sze, N.-D., Rodriguez, J. M., and Scott, C. J.: A two-dimensional model of sulfur species and aerosols, J. Geophys. Res.-Atmos., 102, 13019–13035,, 1997. a

Wittrock, F., Müller, R., Richter, A., Bovensmann, H., and Burrows, J. P.: Measurements of iodine monoxide (IO) above Spitsbergen, Geophys. Res. Lett., 27, 1471–1474,, 2000. a, b

WMO: Scientific assessment of ozone depletion: 2014, WMO, World Meteorological Organization Geneva, Switzerland, Global Ozone Research and Monitoring Project, 2014. a

WMO: Scientific assessment of ozone depletion: 2018. Publisher: WMO, World Meteorological Organization Geneva, Switzerland, GlobalOzone Research and Monitoring Project, 86, 2018. a, b, c, d

Youn, D., Patten, K. O., Wuebbles, D. J., Lee, H., and So, C.-W.: Potential impact of iodinated replacement compounds CF3I and CH3I on atmospheric ozone: a three-dimensional modeling study, Atmos. Chem. Phys., 10, 10129–10144,, 2010. a

Zhang, J., Wuebbles, D. J., Kinnison, D. E., and Saiz-Lopez, A.: Revising the Ozone Depletion Potentials Metric for Short-Lived Chemicals Such as CF3I and CH3I, J. Geophys. Res.-Atmos., 125, e32414,, 2020. a

Zubov, V., Rozanov, E., Egorova, T., Karol, I., and Schmutz, W.: Role of external factors in the evolution of the ozone layer and stratospheric circulation in 21st century, Atmos. Chem. Phys., 13, 4697–4706,, 2013. a


The GV AMAX-DOAS Data can be found here: (last access: 20 October 2021).

Short summary
Here, we present the iodine chemistry module in the SOCOL-AERv2 model. The obtained iodine distribution demonstrated a good agreement when validated against other simulations and available observations. We also estimated the iodine influence on ozone in the case of present-day iodine emissions, the sensitivity of ozone to doubled iodine emissions, and when considering only organic or inorganic iodine sources. The new model can be used as a tool for further studies of iodine effects on ozone.