Seasonal impact of biogenic very short-lived bromocarbons on lowermost stratospheric ozone between 60 N and 60 S during the 21st century

Biogenic very short-lived bromocarbons (VSLBr) currently represent∼ 25 % of the total stratospheric bromine loading. Owing to their much shorter lifetime compared to anthropogenic long-lived bromine (e.g. halons) and chlorine (e.g. chlorofluorocarbons), the impact of VSLBr on ozone peaks in the lowermost stratosphere, which is a key climatic and radiative atmospheric region. Here we present a modelling study of the evolution of stratospheric ozone and its chemical loss within the tropics and at mid-latitudes during the 21st century. Two different experiments are explored: considering and neglecting the additional stratospheric injection of 5 ppt biogenic bromine naturally released from the ocean. Our analysis shows that the inclusion of VSLBr results in a realistic stratospheric bromine loading and improves the agreement between the model and satellite observations of the total ozone column (TOC) for the 1980–2015 period at mid-latitudes. We show that the overall ozone response to VSLBr at mid-latitudes follows the stratospheric evolution of long-lived inorganic chlorine and bromine throughout the 21st century. Additional ozone loss due to VSLBr is maximized during the present-day period (1990–2010), with TOC differences of −8 DU (−3 %) and −5.5 DU (−2 %) for the Southern Hemisphere and Northern Hemisphere midlatitudes (SH-MLs and NH-MLs), respectively. Moreover, the projected TOC differences at the end of the 21st century are ∼ 50 % lower than the values found for the present-day period. We find that seasonal VSLBr impact on lowermost stratospheric ozone at mid-latitude is influenced by the seasonality of the heterogeneous inorganic-chlorine reactivation processes on ice crystals. Indeed, due to the more efficient reactivation of chlorine reservoirs (mainly ClONO2 and HCl) within the colder SH-ML lowermost stratosphere, the seasonal VSLBr impact shows a small but persistent hemispheric asymmetry through the whole modelled period. Our results indicate that, although the overall VSLBr-driven ozone destruction is greatest during spring, the halogen-mediated (Halogx-Loss) ozone loss cycle in the mid-latitude lowermost stratosphere during winter is comparatively more efficient than the HOx cycle with respect to other seasons. Indeed, when VSLBr are considered, Halogx-Loss dominates wintertime lowermost stratospheric ozone loss at SH-MLs between 1985 and 2020, with a contribution of inter-halogen ClOx–BrOx cycles to Halogx-Loss of ∼ 50 %. Within the tropics, a small (<−2.5 DU) and relatively constant (∼−1 %) ozone depletion mediated by VSLBr is closely related to their fixed emissions throughout the modelled period. By including the VSLBr sources, the seasonal Published by Copernicus Publications on behalf of the European Geosciences Union. 8084 J. A. Barrera et al.: Seasonal impact of VSL bromocarbons on lowermost stratospheric ozone Halogx-Loss contribution to lowermost stratospheric ozone loss is practically dominated by the BrOx cycle, reflecting the low sensitivity of very short-lived (VSL) bromine to background halogen abundances to drive tropical stratospheric ozone depletion. We conclude that the link between biogenic bromine sources and seasonal changes in heterogeneous chlorine reactivation is a key feature for future projections of mid-latitude lowermost stratospheric ozone during the 21st century.

Abstract. Biogenic very short-lived bromocarbons (VSL Br ) currently represent ∼ 25 % of the total stratospheric bromine loading. Owing to their much shorter lifetime compared to anthropogenic long-lived bromine (e.g. halons) and chlorine (e.g. chlorofluorocarbons), the impact of VSL Br on ozone peaks in the lowermost stratosphere, which is a key climatic and radiative atmospheric region. Here we present a modelling study of the evolution of stratospheric ozone and its chemical loss within the tropics and at mid-latitudes during the 21st century. Two different experiments are explored: considering and neglecting the additional stratospheric injection of 5 ppt biogenic bromine naturally released from the ocean. Our analysis shows that the inclusion of VSL Br results in a realistic stratospheric bromine loading and improves the agreement between the model and satellite observations of the total ozone column (TOC) for the 1980-2015 period at mid-latitudes. We show that the overall ozone response to VSL Br at mid-latitudes follows the stratospheric evolution of long-lived inorganic chlorine and bromine throughout the 21st century. Additional ozone loss due to VSL Br is maximized during the present-day period , with TOC differences of −8 DU (−3 %) and −5.5 DU (−2 %) for the Southern Hemisphere and Northern Hemisphere midlatitudes (SH-MLs and NH-MLs), respectively. Moreover, the projected TOC differences at the end of the 21st century are ∼ 50 % lower than the values found for the present-day period.
We find that seasonal VSL Br impact on lowermost stratospheric ozone at mid-latitude is influenced by the seasonality of the heterogeneous inorganic-chlorine reactivation processes on ice crystals. Indeed, due to the more efficient reactivation of chlorine reservoirs (mainly ClONO 2 and HCl) within the colder SH-ML lowermost stratosphere, the seasonal VSL Br impact shows a small but persistent hemispheric asymmetry through the whole modelled period. Our results indicate that, although the overall VSL Br -driven ozone destruction is greatest during spring, the halogen-mediated (Halog x-Loss ) ozone loss cycle in the mid-latitude lowermost stratosphere during winter is comparatively more efficient than the HO x cycle with respect to other seasons. Indeed, when VSL Br are considered, Halog x-Loss dominates wintertime lowermost stratospheric ozone loss at SH-MLs between 1985 and 2020, with a contribution of inter-halogen ClO x -BrO x cycles to Halog x-Loss of ∼ 50 %.
Within the tropics, a small (< −2.5 DU) and relatively constant (∼ −1 %) ozone depletion mediated by VSL Br is closely related to their fixed emissions throughout the modelled period. By including the VSL Br sources, the seasonal Halog x-Loss contribution to lowermost stratospheric ozone loss is practically dominated by the BrO x cycle, reflecting the low sensitivity of very short-lived (VSL) bromine to background halogen abundances to drive tropical stratospheric ozone depletion. We conclude that the link between biogenic bromine sources and seasonal changes in heterogeneous chlorine reactivation is a key feature for future projections of mid-latitude lowermost stratospheric ozone during the 21st century.

Introduction
The role of bromine in stratospheric ozone depletion has been discussed in several studies (Wofsy et al., 1975;Prather and Watson, 1990;Daniel et al., 1999;Dvortsov et al., 1999;Solomon, 1999). Although bromine is much less abundant than chlorine in the atmosphere, it is known to deplete stratospheric ozone 45 to 69 times more efficiently on a per atom basis (Daniel et al., 1999;Sinnhuber et al., 2009). Moreover, the ozone depletion efficiency of active bromine (Br and BrO) is strongly related to the available amount of activated chlorine (mainly Cl and ClO radicals) in the atmosphere (McElroy et al., 1986;Solomon et al., 1999, and references therein) as well as to enhanced sulfate aerosol loading via heterogeneous reactions (Salawitch et al., 2005). Consequently, even low concentrations of bromine have a relatively large impact on stratospheric ozone. In addition to anthropogenic long-lived chlorine (LL Cl ) and bromine (LL Br ) substances such as chlorofluorocarbons (CFCs), chlorocarbons, methyl bromide, and bromofluorocarbons (halon fire suppressants), other substances with photochemical lifetimes shorter than 6 months, often referred to as very short-lived (VSL) substances, have the potential to transport a significant amount of reactive halogens into the stratosphere. Owing to their short lifetimes, the impact of VSL substances on stratospheric ozone peaks in the lowermost stratosphere (Salawitch et al., 2005;Feng et al., 2007;Sinnhuber et al., 2009;Sinnhuber and Meul, 2015;Yang et al., 2014;Hossaini et al., 2015a;Falk et al., 2017), which is an important atmospheric region because surface temperature and climate are most sensitive to ozone perturbations (Riese et al., 2012;Hossaini et al., 2015b;Iglesias-Suarez et al., 2018). In fact, Hossaini et al. (2015b) found that very short-lived bromocarbons (VSL Br ) exert an ozone radiative effect (normalized by halogen content) that is 3.6 times larger than that arising from long-lived substances.
Production of anthropogenic LL Cl and LL Br substances has been restricted by the Montreal Protocol and its subsequent amendments and adjustments (Solomon, 1999;WMO, 2018). After application, the stratospheric LL Cl and LL Br load peaked at the end of the 20th century and has been decreasing at a rate that depends on their respective atmospheric lifetimes. This implies, for example, that hydrogen chloride (HCl) shows a long-term decrease at a rate of about 0.5 % yr −1 in the middle stratosphere between 60 • N and 60 • S, while total stratospheric bromine decreased at a rate of about 0.75 % yr −1 from 2004 to 2014 (WMO, 2018). Accordingly, stratospheric ozone is expected to recover from the effects of anthropogenic halogen-induced loss on a similar timescale, which is already detectable in the Antarctic and the upper stratosphere (e.g. Solomon et al., 2016a;Chipperfield et al., 2017;Dhomse et al., 2018;Strahan and Douglass, 2008, and references therein). Hossaini et al. (2015a) quantified the stratospheric injection of organic and inorganic chlorine from anthropogenic VSL Cl sources such as chloroform (CHCl 3 ), dichloromethane (CH 2 Cl 2 ), tetrachloroethene (C 2 Cl 4 ), trichloroethene (C 2 HCl 3 ), and 1,2dichloroethane (CH 2 ClCH 2 Cl), which are not controlled by the Montreal Protocol. From their results, the total stratospheric chlorine load from VSL Cl inferred for 2013 is 123 ppt, with a stratospheric injection dominated by source gases (∼ 83 %). Moreover, the stratospheric injection of VSL Cl increased by ∼ 52 % between 2005 and 2013, mainly due to recent and ongoing growth in anthropogenic CH 2 Cl 2 emissions. In fact, Hossaini et al. (2017) showed that the impact of CH 2 Cl 2 on stratospheric ozone has increased markedly in recent years, and if these increases continue into the future, the recovery of the Earth's ozone layer could be delayed even further, offsetting some of the gains achieved by the Montreal Protocol.
The total stratospheric chlorine and bromine budgets derived for 2016 were 3.29 ppb Cl and 19.60 ppt Br, respectively (WMO, 2018), and are expected to return to their 1980 values, an arbitrary reference date before the discovery of the Antarctic ozone hole, around the middle of this century (Dhomse et al., 2018). Although we still lack scientific consensus with respect to the future evolution of VSL Br oceanic source strength and stratospheric injection (WMO, 2014;Lennartz et al., 2015;Ziska et al., 2017), the continuous decrease in LL Br leads to an increase in the relative VSL Br contribution to total stratospheric bromine into the future. Thus, understanding the role of natural VSL Br sources is a key issue for chemistry-climate projections.
The impact of additional bromine is closely related to the heterogeneous chemistry of chlorine (Solomon, 1999;Salawitch et al., 2005;Müller, 2012). The influence of temperature and water vapour in the heterogeneous conversion processes of inorganic-chlorine reservoirs (mostly HCl and ClONO 2 ) into active radicals on sulfate aerosols has been reviewed elsewhere (e.g. Anderson et al., 2012Anderson et al., , 2017Drdla and Müller, 2012;Anderson and Clapp, 2018). In particular, Drdla and Müller (2012) highlighted that an increase in water vapour above background values would allow chlorine activation at higher temperatures than those observed in polar regions, which led to the hypothesis that chlorine reactivation and subsequent ozone loss could occur in the lower stratosphere at mid-latitudes during summer (Robrecht et al., 2019). Indeed, the spatial and seasonal distributions of chlorine monoxide and chlorine nitrate in the monsoon regions provide indicators of heterogeneous chlorine processing in the tropical and subtropical lowermost stratosphere of the Northern Hemisphere (Solomon et al., 2016b). Moreover, heterogeneous chorine reactivation has also been observed on ice particles in cirrus clouds located near the tropopause (Borrman et al., 1996(Borrman et al., , 1997Solomon et al., 1997;Bregman et al., 2002;Thornton et al., 2003). The occurrence of cirrus clouds was reported in the lowermost stratosphere of northern high-mid-latitudes (40-75 • N; Spang et al., 2015), with summertime clouds located up to ∼ 3 km (or 40-50 K of potential temperature) above the tropopause (Dessler, 2009). Furthermore, von Hobe et al. (2011 suggested that cirrus ice particles and/or liquid aerosol at low temperatures may promote a significant activation of heterogeneous chlorine in the tropical upper troposphere and lower stratosphere. Finally, satellite observations indicate that chlorine activation also occurs in the Antarctic and Arctic sub-vortex regions, where processed air dispersing from the decaying vortex in spring induces rapid changes in extra-vortex trace gas abundances (Santee et al., 2011).
Previous chemistry-climate modelling studies considering VSL Br chemistry have mainly focused on improving the model vs. observed mid-latitude ozone trends at the end of the 20th century (Feng et al., 2007;Sinnhuber et al., 2009). More recently, Sinnhuber and Meul (2015) have shown that the model simulations including VSL Br resulted in a much greater ozone decrease in the lowermost stratosphere during the 1979-1995 period and a faster ozone increase during the 1996-2005 period, which is in better agreement with observed ozone trends. Moreover, Falk et al. (2017) have reported that the projected additional VSL Br impact on ozone at the end of the 21st century is reduced compared to present day. These results are in agreement with those of Yang et al. (2014), who performed time slice simulations to address the sensitivity of stratospheric ozone to a speculative doubling of VSL Br sources under different LL Cl scenarios.
However, neither of the previous studies had distinguished the contribution to ozone depletion -including its seasonal variability -that arises from VSL Br with respect to LL Cl and LL Br throughout the 21st century.
In this work, using the CAM-chem model, we present a coupled chemistry-climate modelling study from 1960 to 2100 with and without the contribution from VSL Br sources. We focus on natural VSL Br -driven impacts on the temporal evolution of stratospheric ozone and the total ozone column (TOC) in the tropics and at mid-latitudes, distinguishing both the long-term seasonal change and the resulting hemispheric asymmetries. Additionally, we present a timeline assessment of the individual contribution from anthropogenic (LL Cl and LL Br ) and biogenic (VSL Br ) sources to stratospheric halogen loading throughout the 21st century while also recognizing the VSL Br contribution to overall halogen-catalysed ozone loss in the lowermost stratosphere. The layout of the paper is as follows: Sect. 2 describes the main characteristics and configuration of the CAM-chem model as well as the observation data sets that are used for the evaluation of the CAM-chem performance. Section 3 presents a quantitative comparison between the model and observations of stratospheric chlorine and bromine loading (Sect. 3.1) and the TOC (Sect. 3.2) as well as discusses long-term seasonal impacts mediated by VSL Br on the TOC (Sect. 3.2) and the lowermost stratospheric ozone budget (Sect. 3.3). Changes in the seasonal evolution of halogendriven ozone loss rates due to VSL Br are shown in Sect. 3.4. The final concluding remarks are summarized in Sect. 4.

Methods: model description
The Community Earth System Model (CESM1) with the Community Atmospheric Model including interactive chemistry (CAM4-chem version; Lamarque et al., 2012;Tilmes et al., 2016) has been used to explore VSL Br -driven stratospheric ozone loss in three latitude bands: Southern Hemisphere mid-latitudes (SH-MLs; 60-35 • S), tropics (25 • N-25 • S), and Northern Hemisphere mid-latitudes (NH-MLs; 35-60 • N). The model set-up follows the Chemistry-Climate Model Initiative (CCMI) REFC2 configuration described in detail by Tilmes et al. (2016), with the exception that in this work we consider a full halogen chemistry mechanism and oceanic emission -seasonally dependent and geographically distributed -of six bromocarbons (VSL Br = CHBr 3 , CH 2 Br 2 , CH 2 BrCl, CHBrCl 2 , CHBr 2 Cl, and CH 2 IBr) from the Ordóñez et al. (2012) emission inventory. This emission inventory is based on a monthly varying satellite chl a climatology, which allows the introduction of a complete seasonal cycle to the emission strength and spatial distribution of oceanic bromocarbons. A full description of the VSL Br mechanism implemented in CAM-chem -including both biogenic and anthropogenic sources, heterogeneous recycling reactions, dry and wet deposition, convective uplift, and large-scale transport -has been given elsewhere (Ordóñez et al., 2012;Fernandez et al., 2014). Monthly and seasonally varying lower boundary conditions were considered for long-lived chlorine (LL Cl = CH 3 Cl, CH 3 CCl 3 , CCl 4 , CFC-11,  and long-lived bromine (LL Br = CH 3 Br, H-1301, H-1211, H-1202, and H-2402) following the A1 halogenated ozone-depleting substance emissions scenario (WMO, 2011). The surface concentrations of CO 2 , CH 4 , H 2 , and N 2 O were specified following the moderate Representation Concentration Pathway 6.0 (RCP6.0) scenario (Meinshausen et al., 2011;Eyring et al., 2013).
The CAM-chem configuration used in this work extends from the surface to approximately 40 km (3.5 hPa in the upper stratosphere), with 26 vertical levels, and includes a horizontal resolution of 1.9 • latitude by 2.5 • longitude. The number of stratospheric levels changes depending on the latitudinal region: within the tropics, there are eight levels above the tropopause (∼ 100 hPa), with a mean thickness of 1.25 km (15.5 hPa) for the lower stratospheric levels and 5.2 km (3.8 hPa) between the two highest levels, while for the mid-latitudes, the tropopause is located at approximately ∼ 200 hPa, and the stratosphere contains up to 12 levels. To ensure a consistent dynamical description of the stratosphere within the relatively low model top of CAM-chem (i.e. gravity wave dragging and the Brewer-Dobson circulation), the integrated momentum that would have been deposited above the model top is specified by an upper boundary condition (Lamarque et al., 2008). A similar procedure is applied to the altitude-dependent photolysis rate computations, which include an upper boundary condition that considers the ozone column fraction prevailing above the model top. The CAM4chem version used in this work includes a non-orographic gravity wave scheme based on the inertia-gravity wave parameterization and an observation-based implementation of stratospheric aerosol and surface area density. The quasibiennial oscillation is imposed in the model by relaxing equatorial zonal winds between 90 and 3 hPa to the observed interannual variability (see Tilmes et al., 2016, for more details).
Equivalent CAM-chem configurations, such as the one implemented in this work, have been used for the CCMVal-2 (Chemistry-Climate Model Validation 2) and CMIP5 (Coupled Model Intercomparison Project, Phase 5) activities in order to represent trends in the ozone evolution and to estimate the return dates, which lie in the middle of the multi-model range (Eyring et al., 2010(Eyring et al., , 2013. Moreover, our model configuration uses a fully coupled Earth system model approach, i.e. the ocean and sea ice are explicitly computed (Neale et al., 2013). This implies that instead of isolating the chemical impact mediated by VSL Br on the climatological ozone budget (as would be the case using a specified dynamic approach), the chemical interaction between VSL bromine and ozone in the lowermost stratosphere is affected by dynamic feedbacks (i.e. temperature, radiation, etc.).
Even though the interactive ocean coupling would have allowed us to compute the evolution of VSL Br ocean source emissions throughout the modelled period, we still lack a considerable understanding of the seasonal processes that dominate VSL Br ocean emissions in both the present and the future, which is combined with a limited observation data set of VSL bromocarbons (WMO, 2014(WMO, , 2018. Therefore, in this work, we forced the model with annually cycled VSL Br fluxes, replicating the Ordoñez et al. emission inventory for all years between 1950 and 2100. This procedure aims to reduce the uncertainties associated with the unknown evolution of both biogeochemical production and oceanic emissions of VSL Br into the future (Lennartz et al., 2015;Ziska et al., 2017).
Two ensembles of independent modelled experiments were performed, each with three individual simulations from 1960 to 2100, considering approximately 10 years of spin-up to allow stabilization of the stratospheric circulation. The individual simulations 001, 002, and 003 started in 1950, 1951, and 1952, respectively, with initial conditions taken from the equivalent years of the CESM1 WACCM (Whole Atmosphere Community Climate Model) 20th-century ensemble for CMIP5 (Marsh et al., 2013). The baseline simulation set-up of the first experiment considered only the halogen LL Cl and LL Br contribution from anthropogenic chlorofluorocarbons, hydrochlorofluorocarbons, halons, methyl chloride, and methyl bromide, while the simulation set-up of the second experiment included the background biogenic contribution from VSL Br oceanic sources in addition to the anthropogenic sources above. The mean ensemble of each experiment (i.e. with and without the VSL Br sources) was calculated as the average of the three individual simulations. Differences between these two mean ensembles allow quantification of the overall impact of VSL Br on stratospheric ozone as well as the determination of their relative contributions to the halogen-mediated catalytic ozone-depleting families. In addition, from this analysis it is also possible to distinguish the contributions of both LL Br and VSL Br to the inorganic fraction of stratospheric bromine (Br LL+VSL y = Br LL y +Br VSL y = Br+BrO+HBr+BrONO 2 +BrCl+HOBr+ 2Br 2 + BrNO 2 ), whereas the inorganic fraction of chlorine (Cl LL y = ClO+Cl+HOCl+2Cl 2 O 2 +2Cl 2 +OClO+HCl+ ClONO 2 +BrCl+ClONO 2 ) arises only from the degradation of LL Cl , which is identical for both experiments.
For the case of vertical distributions and latitudinal variations, the zonal mean of each mean ensemble was computed from the monthly output before processing the data, while a locally weighted scatter plot smoothing (LOWESS) filter with a 0.2 fraction was applied to all long-term time series to smooth the data. Most of the figures and values in the Results section include geographically averaged quantities for the present day (defined here as the 1990-2010 mean, nominally year 2000) and the end of the 21st century (defined here as the 2080-2100 mean, nominally year 2100) over the latitu-dinal bands defined above. To highlight seasonal changes, we computed the average for the months of December-January-February (DJF), March-April-May (MAM), June-July-August (JJA), and September-October-November (SON). Ozone loss trends due to VSL Br between 2000 and 2100 were calculated on the basis of the least-squares linear trends, while the trend errors were estimated as twice the statistical deviation stemming from the least-squares fit.
To provide a basic evaluation of the model performance, the modelled results of the TOC and the inorganic-halogen abundances were compared with the following selected observation data sets: An important point to highlight about the SBUV MOD (and its uncertainties) is that it was constructed from ozone profiles measured by individual SBUV instruments and retrieved using the Version 8.6 algorithm (Bhartia et al., 2013;Frith et al., 2014Frith et al., , 2017. To maintain consistency over the entire time series, the individual instrument records were analysed with respect to each other, and absolute calibration adjustments were applied as needed based on comparison of radiance measurements during periods of instrument overlap (DeLand et al., 2012;Weber et al., 2018). During the overlap periods, the records of the involved SBUV instruments were combined using a simple average. Frith et al. (2014) have estimated the potential errors in the SBUV MOD through Monte Carlo model simulations, taking into account the range of offsets and drifts observed during the overlap periods between individual SBUV instruments and between SBUV and Dobson ground-based measurements. The authors also assessed the extent to which these potential errors might affect the long-term variability of ozone in a multiple-regression model. See Frith et al. (2014) for a detailed description of the data used and their associated uncertainties in the SBUV MOD.
3 Results and discussion 3.1 Contribution of LL Cl , LL Br , and VSL Br to stratospheric inorganic-halogen loading The evolution (1960-2100) of the modelled annual mean abundances of inorganic chlorine (Cl y ) and bromine (Br y ) as well as the VSL Br in the global upper stratosphere (3.5 hPa) and lower stratosphere (50 hPa) of the tropics and midlatitudes is shown in Fig. 1. In addition, Fig. 1 includes the total Br y trends (WMO, 2014(WMO, , 2018 and the MLS observations (HCl + ClO). The values of bias, normalized mean bias (NMB), and normalized mean error (NME) used to evaluate the agreement between the modelled Cl y and observations for each of the regions under study are shown in Table S1 in the Supplement. Clearly, the temporal evolution of Cl LL y and Br LL y shows a pronounced peak at the end of the 20th century and beginning of the 21st century, respectively, after which the abundance of both halogens declines. Within the A1 halogen emission scenario considered in this work, the Cl LL y abundance in our model returns to its past 1980 levels just before ∼ 2060 in the global upper stratosphere, whereas the prevailing Br LL y abundance for the year 1980 is not recovered even by the end of the 21st century (see Fig. 1a). In comparison, the evolution of Br VSL y abundances remains constant over time, resulting in an additional and time-independent fixed contribution of ∼ 5 ppt total bromine injection, which leads to maximum global upperstratosphere Br LL+VSL y abundances of up to ∼ 20.5 ppt at the beginning of the 21st century. This is in agreement with previous studies performed only for present-day conditions . Furthermore, note that the modelled global Br LL+VSL y abundances are in good agreement with the observations within reported errors, highlighting the importance of considering the additional contribution from VSL Br to determine the overall evolution of stratospheric ozone during the 21st century.
In the lower stratosphere, at 50 hPa, Br LL y abundances return to 1980 levels before ∼ 2080 and ∼ 2050 for the midlatitudes and tropics, respectively, while the Cl LL y return date remains unaltered (see Fig. 1b-d). Furthermore, the contribution from VSL Br to mid-latitude total inorganic bromine also reaches ∼ 5 ppt throughout the modelled period, while in the tropics there prevails a small carbon-bonded organic fraction that has not been converted to the reactive inorganic form. Falk et al. (2017) showed an increase in VSL Br of about 0.5 ppt in the tropical lowermost stratosphere by the end of the century compared to present day under the RCP6.0 scenario and assuming constant VSL Br fluxes, following scenario five of Warwick et al. (2006). This increase in the stratospheric VSL Br attributed to enhanced vertical transport in the tropics is counteracted by a decrease in the Br VSL y injected into the stratosphere so that the total stratospheric VSL bromine remains unchanged between the two periods. In this work, the stratospheric injection of VSL Br and Br VSL y between 2000 and 2100 remains unchanged. Based on the limited observational data, it is uncertain whether any trend in the stratospheric injection of organic and inorganic VSL bromine would be expected, and thus no conclusive statement can be made. Nevertheless, our presentday stratospheric injection of VSL Br (∼ 2 ppt) and Br VSL y (∼ 3 ppt) is in perfect accordance with that reported in the latest WMO Ozone Assessment Report (WMO, 2018). Further studies are needed to evaluate the uncertain evolution of the stratospheric injection of organic and inorganic species from VSL Br throughout the 21st century.
The modelled inorganic-chlorine abundance shows a good agreement with the (HCl + ClO) MLS observations for the 2005-2015 period mainly in the upper stratosphere, where most chlorine has already been photochemically converted to Cl LL y . For example, the normal mean error of the compar-ison between modelled Cl LL y abundances and MLS observations is about 3.5 % and 19 % at 50 hPa for the mid-latitudes and tropics, respectively, while the NME is about 2 % for the global mean at 3.5 hPa (see Table S1). This relatively good agreement between the model and MLS observation occurs even without consideration of the recent contribution of anthropogenic VSL Cl sources (i.e. Hossaini et al., 2015bHossaini et al., , 2017. Moreover, the inter-annual variation in MLS observations at 50 hPa occurs mainly due to heterogeneous chlorine reactivation, which can be influenced by the polar vortex dynamics at mid-latitudes (Dhomse et al., 2018). Note that although we processed output at 120 hPa in the mid-latitude lowermost stratosphere in Sect. 3.4, we did not compare Cl LL y at this altitude level. This is because below 50 hPa, the fractional conversion of organic chlorine to Cl LL y is very small, and thus the inferred HCl + ClO MLS data have a very large uncertainty that prevents a reliable model-observation comparison (Santee et al., 2011).
The stratospheric Br VSL y injection remains constant during the whole modelled period, representing ∼ 25 % of Br LL+VSL y in the global upper stratosphere for the year 2000 and increasing up to ∼ 40 % by the end of the 21st century. These results are in agreement with those reported by Fernandez et al. (2017). Due to the much shorter lifetimes of VSL Br with respect to LL Br , the contribution percentage is larger in the lower stratosphere. For example, the Br VSL y relative contribution to Br LL+VSL y for the year 2000 at 50 hPa represents ∼ 30 % and ∼ 45 % for the mid-latitudes and tropics, respectively. Furthermore, at mid-latitudes, Br VSL y represents up to 45 % of Br LL+VSL y both at the end of the 21st century and in the years prior to 1980, while in the tropics Br VSL y represents up to 65 % for these time periods. Consequently, these changes in the relative contribution of VSL Br to total inorganic bromine with respect to time, region, and altitude in the presence of different inorganic-chlorine abundances lead to relatively different ozone loss within those regions and altitude regimes, as described in detail below. Figure 2 shows the annual zonal mean distribution of inorganic-chlorine and inorganic-bromine abundances as well as the corresponding ClO x /Cl y and BrO x /Br y percentage ratios (black contour lines) for the present day. Note that the ClO x /Cl y and BrO x /Br y ratios reflect the changes in the contribution of the active-chlorine (ClO x = ClO + Cl + HOCl+2Cl 2 O 2 +OClO+2Cl 2 ) and active-bromine (BrO x = BrO + Br) fractions relative to their total inorganic abundances. Equivalent results are presented for the end of the 21st century in the Supplement (see Fig. S1). In the lowermost stratosphere, the available inorganic fraction of both chlorine and bromine rapidly increases with altitude as halogen atoms are released by the photolysis and oxidation from all LL Cl , LL Br , and VSL Br species. In contrast, within the upper stratosphere no major changes in inorganic-halogen abundances are observed since the conversion from their organic sources is nearly complete. Although our model shows a symmetric hemispheric distribution in stratospheric Cl y and Br y abundances, there is a marked difference in the ClO x /Cl y ratio. In fact, the experiments capture a local increase in the ClO x /Cl y ratio at 15 • N and ∼ 100 hPa as a result of a large increase in heterogeneous reactivation of Cl y species due to enhanced vertical transport during the Indian summer monsoon (Solomon et al., 2016b). Along the same lines, the higher ClO x /Cl y values modelled in the lowermost stratosphere at SH-MLs compared to NH-MLs are mainly due to the enhancement in the heterogeneous chlorine reactivation processes on ice crystals (i.e. HCl and ClONO 2 reactivation; Santee et al., 2011;Solomon, 1999, and references therein), influenced mainly by the proximity of the Antarctic polar vortex edge and the lower temperatures prevailing at high-and mid-latitudes of the Southern Hemisphere (see Fig. S2). Additionally, the favourable conditions for isentropic exchange of young subtropical air with extremely old winter polar vortex air masses at SH-MLs (Spang et al., 2015;Rolf et al., 2015) could drive additional heterogeneous chlorine reactivation and thus enhance hemispheric asymmetry. Moreover, the increase in the ClO x abundance and ClO x /Cl y ratio is highest during winter and spring for both SH-MLs and NH-MLs, which is consistent with the seasonal changes of the lowermost stratospheric temperatures that affect chlorine reactivation on ice crystals (see Fig. S3). In contrast, a symmetric hemispheric distribution of the BrO x /Br y ratio is modelled in the mid-latitude lowermost stratosphere and is associated with heterogeneous bromine reactivation on sulfate aerosols instead of ice crystals (i.e.  (Frith et al., 2014(Frith et al., , 2017 within the same spatial and temporal mask as the model output. HBr and BrONO 2 reactivation; Solomon, 1999, and references therein), which show neither any hemispheric asymmetry nor seasonal changes in the reactivation processes (see Fig. S4). Thus, even when the inclusion of VSL Br sources increases the total stratospheric Br y abundance and therefore its active fraction (BrO x ), the BrO x /Br y ratio remains unchanged in the stratosphere during the whole modelled period (see Figs. 2b, c and S1b, c), highlighting that this ratio is nearly independent of the total inorganic bromine abundance.

Impact of VSL Br on the seasonal evolution of the total ozone column (TOC)
The temporal evolution (1960-2100) of the modelled annual mean TOC at mid-latitudes and in the tropics, along with SBUV MOD observations, is illustrated in Fig. 3. Equivalent results for the temporal evolution of the individual simulations of each experiment are shown in Fig. S5. Table 1 shows the annual and seasonal mean values of the absolute and relative TOC differences ( TOC) between the experiments for the present day and the end of the 21st century as well as the TOC trends (% per decade) over the century. In addition, Table S2 shows the values of bias, normalized mean bias (NMB), and normalized mean error (NME) to evaluate the agreement between the modelled TOC and observations for each of the regions under study. Based on the comparison of the TOC results over the different latitudinal bands, the following features can be described: i. The constant emission of VSL Br introduces a continuous reduction in the TOC that exceeds the mean ensemble variability between the experiments (see Fig. S5). Moreover, the inclusion of VSL Br sources improves the overall agreement between the model and observations, with a reduction of 1 % and 0.9 % in NME for SH-MLs and NH-MLs, respectively (see Table S2). However, the inclusion of VSL Br sources leads to a slight increase of 0.2 % in NME for the tropics. Note that the statistical analysis presented in Table S2 does not take into account the uncertainties that may arise from the merging of the individual SBUV instruments in the SBUV MOD (see Sect. 2).
for the present day, which is in line with the projected TOC trends of 0.15 % and 0.10 % per decade for SH-MLs and NH-MLs, respectively.
iii. Within the tropics, the inclusion of VSL Br leads to small impacts on the TOC compared to impacts observed at mid-latitudes. The TOC differences due to VSL Br are close to −1.5 DU (< −1 %) during practically the entire modelled period, with a maximum TOC 2000 reaching −2.1 DU (< −1 %). This is in line with a projected TOC trend of 0.03 % per decade. Moreover, even though VSL Br slightly worsen the agreement between the modelled TOC LL and observations, the minimum TOC LL+VSL is shifted to ∼ 5 years earlier compared to TOC LL (i.e. towards where the inorganic-chlorine peak is located in the year 1995), which is in agreement with observations. iv. Overall, our model results show much higher ozone depletion due to VSL Br at mid-latitudes compared to the tropics.
Our results of the mid-latitude TOC due to VSL Br lie within the lower range of previous modelling studies over the 1960(Sinnhuber and Meul, 2015) and 1980(Feng et al., 2007 periods. This is expected as the detailed treatment of VSL Br sources in Sinnhuber and Meul (2015) results in an additional Br VSL  Feng et al. (2007), the model configurations included the direct supply of tropospheric Br y as a lower condition in addition to a VSL substance tracer. Furthermore, our results are in accord with the future projection trends in Falk et al. (2017), although their analysis on annual mean ozone loss due to VSL Br only focused on the late 21st century (2075-2100). Figure 4 shows the temporal evolution (1960-2100) of the annual and seasonal TOC between the experiments for SH-MLs, the tropics, and NH-MLs. In addition, an analysis of TOC as a function of latitude for the present day and the end of the 21st century is shown in Fig. 5. The inclusion of VSL Br leads to a continuous reduction in the TOC at all latitudes, with a larger ozone destruction efficiency moving from the tropics to the high latitudes. Moreover, note that when comparing the seasonal relative TOC between both periods, a statistically significant difference (here defined as p < 0.05) is only observed at mid-latitudes, as shown by the horizontal lines at the bottom of Fig. 5d. This highlights that the VSL Br -driven ozone depletion efficiency changes significantly over the century even though the contribution from VSL Br to bromine stratospheric injection remains constant. Since gas-phase and heterogeneous inter-halogen reactions involving both bromine and chlorine play a fundamental role in stratospheric ozone loss (McElroy et al., 1986;Solomon, 1999;Salawitch et al., 2005), the VSL bromine efficiency in ozone depletion is primarily linked to the background inorganic-halogen abundance, which shows a continuous decline over the course of the 21st century due to the A1 halogen emission scenario considered in this work. This is in line with the findings of Yang et al. (2014). Therefore, an additional stratospheric chlorine and bromine load from natural or anthropogenic substances not regulated by the Montreal Protocol will induce an increase in the VSL Br impacts mod-elled in this work on projected ozone loss trends throughout the current century.
The largest mid-latitude TOC differences between experiments occur during spring, when the maximum seasonal TOC levels are found, followed by the winter. The maximum springtime TOC 2000 reaches −10 DU (∼ −3 %) and −7.7 DU (∼ −2 %) for SH-MLs and NH-MLs, respectively (see Table 1 and Fig. 4). In contrast, during summer and autumn, the maximum TOC 2000 values remain below −6.8 (∼ −2.4 %) and −4 DU (∼ −4.1 %) for SH-MLs and NH-MLs, respectively. This is in line with the ozone loss mediated by VSL Br within the southern polar cap reported by Fernandez et al. (2017), who showed a maximum ozone depletion of up to −15 DU (∼ −14 % of TOC) during October at the beginning of the century. Moreover, as mentioned above, the projected seasonal VSL Br impact on the TOC significantly decreases towards the end of the century, with springtime TOC trends of up to 0.20 % and 0.14 % per decade for SH-MLs and NH-MLs, respectively (see Table 1). In comparison, the projected TOC trends for summer and autumn reach approximately 0.10 % per decade at mid-latitudes. Note also that the changes in the magnitudes of these seasonal TOC trends are mainly attributed to a marked seasonal difference of ozone loss mediated by VSL Br during the period of highest background halogen abundance (i.e. present day). These seasonal differences are reduced towards the end of the 21st century due to the decline in longlived inorganic-halogen abundance.
Within the tropics, no significant changes in the TOC loss due to VSL Br are projected, though the model captures a slightly major TOC depletion during the 1990s. Therefore, VSL Br -driven ozone depletion is less sensitive to background halogen abundances, introducing an almost constant impact that is in line with its fixed emissions considered throughout the whole modelled period.

VSL Br impact on vertical ozone distributions
Timeline analysis of the annual mean stratospheric ozone difference between the experiments as a function of altitude at mid-latitudes and in the tropics is presented in Fig. 6. The inclusion of VSL Br produces a reduction in ozone concentrations (i.e. O 3 (z) number densities) throughout most of the stratosphere during the whole modelled period. The largest O 3 (z) reductions occur during the present-day period in the mid-latitude lowermost stratosphere, reaching ozone differences ( O 3 (z)) of up to −8 % (−0.25 × 10 12 molecules cm −3 ) and −5 % (−0.14 × 10 12 molecules cm −3 ) for SH-MLs and NH-MLs, respectively. Above 50 hPa (∼ 20 km), O 3 (z) values are practically constant at around −1 % for all latitudes, which is linked to the lesser role played by VSL bromine compared to catalytic ozone loss by hydrogen (HO x = OH + HO 2 ) and nitrogen (NO x = NO + NO 2 ) oxides as well as by long-lived reactive halogens (Brasseur and Solomon, 2005;Salawitch et al., 2005;Müller, 2012).
Our modelled O 3 (z) latitudinal distribution is in agreement with those reported by Sinnhuber and Meul et al. (2015) during the 1980-2005 period, with a larger percentage of ozone loss deepening around the mid-latitude tropopause. In addition, O 3 (z) reduction due to VSL Br in the midlatitude lowermost stratosphere is about 50 % smaller by the end of the 21st century compared to the present-day period (see Fig. 6d and e). This is in line with projected local O 3 (z) trends of up to 1 % and 0.5 % per decade at SH-MLs and NH-MLs, respectively (see Fig. 6f). The largest O 3 (z) reductions are located precisely at the same altitudes and latitudes where the ClO x /Cl y ratio peak is modelled within the lowermost stratosphere at SH-MLs (see Fig. 2), which explains the hemispheric asymmetry in the VSL Br impacts. This highlights that the enhancement in local ozone depletion mediated by VSL Br is largest in the context of maximum background active-chlorine abundance since they provide an additional partner for chlorine species involved in the ozone-depleting inter-halogen chemical reactions (e.g. ClO x -BrO x-Loss ; see Sect. 3.4). These results are in agreement with Yang et al. (2014), who determined that VSL Br have a significantly larger ozone impact in the presence of a high stratospheric chlorine background than a low chlorine background using time slice sensitivity simulations with a variable and speculative VSL Br contribution on different background stratospheric Cl y abundances. Figure 7 shows the seasonal zonal mean distribution of the O 3 (z) during the present-day period. The corresponding O 3 (z) trends (% per decade) over the century are shown in Fig. S6. The largest contributions to the annual mean ozone loss in the mid-latitude lowermost stratosphere correspond to spring, followed by the preceding winter and the posterior summer in each hemisphere, with springtime O 3 (z) reaching up to −10 % and −7 % for SH-MLs and NH-MLs, respectively (see Fig. 7b and d). This seasonal enhancement in VSL Br -mediated local ozone depletion is influenced by the seasonal changes in the heterogeneous chlorine reactivation processes on lowermost-stratospheric ice crystals since heterogeneous bromine reactivation occurs mainly in sulfate aerosols and therefore lacks seasonal changes (see Figs. S3 and S4). Consequently, since the seasonal changes in ozone loss are maximized during the present-day period, the largest local O 3 (z) trends are also projected for spring and winter, with springtime O 3 (z) trends reaching up to 1.6 % and 0.7 % per decade for SH-MLs and NH-MLs, respectively (see Fig. S6b and d).
Compared to the annual mean, springtime absolute O 3 (z) reductions along mid-latitudes are deepened over the first ∼ 5 km above the tropopause during the present-day period, whereas these reductions are located at higher altitudes (between 2 and 5 km above the tropopause) and over a narrower geographic area, mainly at NH-MLs, during summer and autumn. Therefore, VSL bromine introduces seasonal changes not only in overall O 3 (z) reduction but also in vertical O 3 (z) distribution. Moreover, since seasonal VSL Br impacts persist until the end of the century, they could lead to significantly different seasonal perturbations in the radiative effects mediated by ozone. Future studies are needed to explore the potential VSL Br -mediated seasonal effects on the atmosphere's radiative balance.
Near the tropical tropopause, the O 3 (z) reductions due to VSL Br are practically constant during the whole modelled period, with O 3 (z) of approximately −2 % (−0.06 × 10 12 molecules cm −3 ) at 50 hPa (see Fig. 6b). Even with an increase in the contribution of Br VSL y to Br LL+VSL y of up to 65 % by the end of century -due to the implementation of the Montreal Protocol (WMO, 2014(WMO, , 2018 and the fixed stratospheric Br VSL y injection -O 3 (z) changes in the tropical lowermost stratosphere are not significant over the century (see Fig. 6f). This is in line with the modelled impacts on the TOC. In contrast, at mid-latitudes, the seasonal VSL Br impacts on ozone significantly decrease over the cen- tury between 50 and 300 hPa (see Figs. 6f and S6a-d) even though the contribution percentage from Br VSL y to Br LL+VSL y increases towards the end of the century. Moreover, the projected seasonal and annual O 3 (z) trends peak precisely at the same altitudes where the modelled relative O 3 (z) values are significant (see . These results highlight the clear dependence of VSL Br impacts on the temporal evolution of background halogen abundances and the processes of heterogeneous halogen reactivation, which in turn depend on temperature as well as the formation and distribution of ice crystals and stratospheric sulfate aerosols. In fact, the projected impacts on vertical ozone distribution at the end of the century are similar to those prevailing before 1980 for both SH-MLs and NH-MLs since the background halogen abundances are similar between these periods. Thus, potential future changes in the lowermost stratospheric temperature (therefore influencing ice crystal formation), geoengineering of climate via injection of stratospheric sulfate, and the formation and distribution of stratospheric aerosol (e.g. future volcanic eruptions) will directly influence the processes of heterogeneous halogen reactivation and consequently VSL Br -driven ozone depletion efficiency (Tilmes et al., 2008Banerjee et al., 2016;Klobas et al., 2017;Heckendorn et al., 2009).

Changes in the long-term seasonal evolution of halogen-driven ozone loss rate due to VSL Br
In order to determine the drivers controlling the VSL Brdriven seasonal impact on ozone in the lowermost stratosphere at mid-latitudes and in the tropics, we compute the odd-oxygen chemical loss from each ozone-depleting family considering the independent contributions of oxygen (O x-Loss ), hydrogen (HO x-Loss ), nitrogen (NO x-Loss ), and halogen (Halog x-Loss ; Brasseur and Solomon, 2005;Saiz-Lopez et al., 2014). Moreover, we discriminate between the individual contributions of pure bromine (BrO x-Loss ) and chlorine (ClO x-Loss ) cycles and the inter-halogen (ClO x -BrO x-Loss ) cycle to the halogen family (Halog x-Loss = BrO x-Loss + ClO x-Loss + ClO x − BrO x-Loss ). The odd-oxygen loss rate equations for the ozone-depleting families considered in this work are presented in Table S3. Figure 8 shows the temporal evolution of the annual contribution percentage of each ozone-depleting family to the total ozone loss rate in the lowermost stratosphere at SH-MLs (∼ 120 hPa) for both experiments as well as the seasonal evolution of HO x-Loss and Halog x-Loss contributions in each experiment. The corresponding figures for the lowermost stratosphere at NH-MLs (∼ 120 hPa) and in the tropics (∼ 50 hPa) are shown in Figs. S7 and S8, respectively. Halog x-Loss represents the second most important contribution to the total ozone loss rate after HO x-Loss (∼ 80 %) at mid-latitudes in both experiments (see Figs. 8a and S7a), while within the tropics it represents the third family after HO x-Loss and NO x-Loss (see Fig. S8a). This partially explains the smaller modelled VSL bromine impact on ozone within the tropics as Halog LL+VSL x-Loss represents at most 10 % of the total ozone destruction over this region. Moreover, the enhancement in the Halog LL+VSL x-Loss contribution with respect to Halog LL x-Loss by including VSL Br is mostly compensated by a decrease in the HO LL+VSL x-Loss contribution at midlatitudes. For example, for the present day, Halog LL+VSL

x-Loss
shows an increase of ∼ 7 % with respect to Halog LL x-Loss at SH-MLs, which is compensated by a reduction of ∼ 6 % in

HO LL+VSL
x-Loss . In particular, the interaction between different catalytic ozone-depleting cycles involving OH and halogen radicals as well as the influence of reactive halogens in altering the OH/HO 2 ratio has been described elsewhere (Bloss et al., 2005;Saiz-Lopez and von Glasow, 2012).
Even though the additional seasonal ozone depletion due to VSL Br peaks during spring followed by winter at midlatitudes, the largest Halog x-Loss contribution to the total ozone loss rate occurs during winter in both experiments, presenting a much smaller contribution during the rest of the year (see Figs. 8b,c and S7b,c). Indeed, this marked seasonal behaviour prevails both during boreal winter at NH-MLs and during the austral winter at SH-MLs. This is explained by considering the strong seasonal increase in active-chlorine abundance (driven by heterogeneous reactivation on ice crystals) as well as the winter decline in HO x abundance (driven by the decrease in solar radiation with respect to summer-time). In fact, the enhancement in the Halog LL+VSL x-Loss contribution due to VSL Br makes it the dominant ozone-depleting family during winter between 1985 and 2020 at SH-MLs, surpassing in importance the role played by the otherwise dominant HO x-Loss . In contrast, no seasonal variability of Halog x-Loss contribution is modelled within the tropics in either experiment (see Fig. S8b and c), and thus the seasonal VSL Br impact on ozone remains approximately constant during the whole modelled period around the annual mean. Figure 9 shows the long-term evolution of the contribution percentage of BrO x-Loss , ClO x-Loss , and ClO x -BrO x-Loss to the Halog x-Loss as a function of months and years in the lowermost stratosphere at SH-MLs (120 hPa) for both experiments. Equivalent results are presented for NH-MLs and the tropics in Figs. S9 and S10, respectively. For any fixed year in both experiments, the BrO x-Loss contribution peaks during late summer and early autumn, while its minimum contribution occurs during the austral winter. In contrast, both the ClO x-Loss and ClO x -BrO x-Loss reach their maximum contribution during winter and early spring. Consequently, the inclusion of VSL Br increases the BrO LL+VSL x-Loss contribution in all seasons while also increasing the ClO x -BrO LL+VSL x-Loss contribution mainly in winter. This is at the expense of a decrease in the ClO LL+VSL x-Loss contribution. In addition, the wintertime ClO x -BrO LL+VSL x-Loss contribution reaches up to 50 % between 1990 and 2020, decreasing towards the end of the century as Cl LL y and Br LL y abundances decrease within the A1 emission scenario. The decrease in both ClO LL+VSL

x-Loss
and ClO x -BrO LL+VSL x-Loss over time is offset by an increase in the BrO LL+VSL x-Loss contribution. Accordingly, the inclusion of VSL Br produces similar changes in the evolution of the contributions of each of the cycles that compose Halog x-Loss at NH-MLs, although the seasonal BrO x-Loss contribution is greater in both experiments compared to SH-MLs during the entire modelled period (see Fig. S9).
In summary, the Halog LL+VSL x-Loss contribution to the total ozone loss rate in summer and autumn is almost entirely dominated by BrO LL+VSL x-Loss . In the winter months, when the seasonal Halog LL+VSL x−Loss contribution to the total ozone loss rate is comparatively higher, a transition from the ClO x BrO LL+VSL x-Loss contribution (∼ 50 %) at the beginning of the century towards BrO LL+VSL x-Loss at the end of the century is modelled. Thus, the halogen-driven ozone loss projected towards the end of the 21st century is mostly dominated by BrO x-Loss , leading to a less marked seasonal VSL Br impact on both vertical ozone distribution and the TOC, in accordance with the modelled results. On the other hand, within the tropics, BrO LL+VSL x-Loss dominates the seasonal Halog LL+VSL x-Loss contribution for practically the entire modelled period (see Fig. S10), while ClO LL+VSL x-Loss represents a maximum contribution of up to ∼ 35 % during the end of the 20th century. The seasonal increase from BrO LL x-Loss to BrO LL+VSL x-Loss by including VSL Br is well marked and is prac- tically offset by a reduction in the ClO LL+VSL x-Loss contribution. Furthermore, unlike mid-latitudes, ClO x -BrO LL+VSL x-Loss has a lower contribution to the tropical Halog LL+VSL x-Loss , showing that VSL Br -driven topical ozone loss is independent of changes in the background halogen abundance throughout the modelled period.

Conclusions
This study has explored the impact of VSL Br on the temporal evolution of stratospheric ozone both in the late 20th century and in the course of the 21st century at mid-latitudes and in the tropics. The analysis compares two independent experiments, one of which considers only the anthropogenic LL Cl and LL Br source contributions, while the other considers -in addition to long-lived sources -the contributions of VSL Br oceanic sources. We have evaluated annual and seasonal mean changes in vertical stratospheric ozone distribution and the TOC as well as examined the projected ozone depletion ( O 3 and TOC) trends over the century under a scenario in which both the relative contribution of VSL Br to total inorganic-bromine and the background inorganic-chlorine abundances shift towards the future. Finally, we have also assessed the long-term seasonal contribution changes of halogen-mediated ozone loss (Halog x-Loss ) due to VSL bromine, focusing on the lowermost stratosphere and distinguishing the role of pure (BrO x-Loss , ClO x-Loss ) and inter-halogen (ClO x -BrO x-Loss ) cycles that compose Halog x-Loss . The results presented here highlight the importance of considering natural sources of halocarbons to de-termine the overall evolution of stratospheric ozone in the tropics and at mid-latitudes during the 21st century.
Our analysis shows that the inclusion of VSL Br results in a realistic stratospheric bromine loading that improves the agreement between the model and observations of the TOC at mid-latitudes, highlighting the need to consider these additional natural emissions to determine the overall evolution of ozone during the 21st century. Our results are in accordance with previous modelling studies that included natural bromocarbons and explored stratospheric ozone in the recent past (Salawitch et al., 2005;Feng et al., 2007; as well as by the end of the 21st century (Falk et al., 2017). However, unlike these previous works, we have explored seasonal ozone depletion arising from VSL Br with respect to the background stratospheric LL Cl and LL Br abundances between 1960 and 2100.
At mid-latitudes, VSL Br introduce a continuous reduction in the TOC during the entire modelled period, with the largest TOC depletion occurring during the period of highest background halogen abundance ( TOC 2000 = −8 DU, −2.5 %, for SH-MLs and −5.5 DU, −1.6 %, for NH-MLs). In turn, a significant decrease in the impact towards the end of the 21st century at mid-latitudes is projected, with TOC trends of 0.15 % and 0.10 % per decade for SH-MLs and NH-MLs, respectively. Even though the contribution of VSL bromine to total stratospheric bromine reaches almost 50 % by the end of the 21st century, our simulations show a statistically significant smaller VSL Br impact on lowermost stratospheric ozone as we move into the future. This result highlights that VSL Br -driven ozone loss is closely linked to the temporal evolution of LL Cl and LL Br abundances regulated by the Montreal Protocol, as suggested by Yang et al. (2014). Along these same lines, Hossaini et al. (2015a) demonstrated that the chlorine load in the lowermost stratosphere has increased as a consequence of a recent and ongoing growth in emissions from anthropogenic VSL Cl sources (mainly CH 2 Cl 2 ) not controlled by the Montreal Protocol. If this VSL Cl emissions trend continues into the future, additional studies will be required to quantify the increase in VSL Br impact modelled in this work on the projected ozone loss trends over the century considering the additional inorganic chlorine from VSL Cl .
The largest TOC depletion associated with VSL Br corresponds to the spring months followed by winter, with maximum springtime TOC 2000 reaching −10 DU (∼ −3 %) and −7.7 DU (∼ −2 %) for SH-MLs and NH-MLs, respectively. We find that the inclusion of VSL Br leads to seasonal changes in both local O 3 (z) reduction and its vertical O 3 (z) distribution within the lowermost stratosphere, which could result in different seasonal perturbations in the radiative effects mediated by ozone. This seasonal enhancement in VSL Brmediated local ozone loss is strongly influenced by the seasonal changes in the processes of heterogeneous chlorine activation on ice crystals (mainly by ClONO 2 and HCl recycling). In fact, the modelled additional increase in het-erogeneous chlorine reactivation efficiency in the lowermost stratosphere at SH-MLs compared to NH-MLs leads to hemispheric asymmetry in the VSL Br -driven ozone reductions.
Even though the additional seasonal ozone depletion due to VSL Br peaks during spring followed by winter, the largest Halog LL+VSL x-Loss contribution to the total ozone loss rate occurs during winter, with a much smaller contribution during all other seasons. In fact, the wintertime Halog LL+VSL x-Loss contribution dominates the total ozone loss rate between 1985 and 2020 in the lowermost stratosphere at SH-MLs, with a ClO x -BrO LL+VSL x-Loss contribution to Halog LL+VSL x-Loss reaching up to 50%. Due to the expected decrease in the LL Cl and LL Br abundances, a transition from ClO x -BrO LL+VSL x-Loss to BrO x-Loss during the course of the 21st century is projected, leading to a less marked seasonal VSL Br impact towards the end of the century.
Within the tropics, although our model captures a slightly greater TOC depletion associated with VSL Br during the 1990s, no significant change is projected by the end of the 21st century. Halog x-Loss represents the third family in the contribution to the total ozone loss rate after HO x-Loss and NO x-Loss in the tropical lowermost stratosphere, with an almost constant seasonal contribution throughout the modelled period. BrO LL+VSL x-Loss practically dominates the contribution to Halog LL+VSL x-Loss , with a maximum ClO x -BrO LL+VSL x-Loss contribution of up to ∼ 30 % by the end of the 20th century, reflecting the low sensitivity of VSL bromine to background halogen abundances to drive stratospheric ozone depletion.
Data availability. Computing resources, support, and data storage are provided and maintained by the Computational and Information System Laboratory from the National Center of Atmospheric Research (CISL, 2017). The code of the CAM-chem model can be downloaded from the website (https://www2.acom.ucar.edu/gcm/ cam-chem (last access: 10 March 2018) (Atmospheric Chemistry Observations and Modelling Division repository, 2018). Data that support the findings of this study can be downloaded from the AC2 web page (https://ac2.iqfr.csic.es/en/publications, last access: 25 June 2020) (Atmospheric Chemistry and Climate Department, 2020).
Author contributions. RPF and ASL designed the experiments. RPF and JAB configured and ran all the simulations. JAB processed the results of all simulations, led the manuscript, and prepared the figures. RPF, ASL, FIS, and CAC provided additional feedback on geophysical processing and representation. JAB, RPF, FIS, CAC, JFL, and ASL were involved in the discussion and iterations of the manuscript.