Journal cover Journal topic
Geoscientific Model Development An interactive open-access journal of the European Geosciences Union
Journal topic
GMD | Articles | Volume 12, issue 3
Geosci. Model Dev., 12, 1227-1239, 2019
© Author(s) 2019. This work is distributed under
the Creative Commons Attribution 4.0 License.
Geosci. Model Dev., 12, 1227-1239, 2019
© Author(s) 2019. This work is distributed under
the Creative Commons Attribution 4.0 License.

Development and technical paper 29 Mar 2019

Development and technical paper | 29 Mar 2019

Improvements to stratospheric chemistry scheme in the UM-UKCA (v10.7) model: solar cycle and heterogeneous reactions

UKCA stratospheric chemistry
Fraser Dennison1,a, James Keeble2,3, Olaf Morgenstern1, Guang Zeng1, N. Luke Abraham2,3, and Xin Yang4 Fraser Dennison et al.
  • 1National Institute of Water and Atmospheric Research, Wellington, New Zealand
  • 2Centre for Atmospheric Science, Department of Chemistry, University of Cambridge, Cambridge, UK
  • 3National Centre for Atmospheric Science, UK
  • 4British Antarctic Survey, Cambridge, UK
  • anow at: CSIRO Oceans and Atmosphere, Aspendale, Australia
Back to toptop

Improvements are made to two areas of the United Kingdom Chemistry and Aerosol (UKCA) module, which forms part of the Met Office Unified Model (UM) used for weather and climate applications. Firstly, a solar cycle is added to the photolysis scheme. The effect on total column ozone of this addition was found to be around 1 %–2 % in midlatitude and equatorial regions, in phase with the solar cycle. Secondly, reactions occurring on the surfaces of polar stratospheric clouds and sulfate aerosol are updated and extended by modification of the uptake coefficients of five existing reactions and the addition of a further eight reactions involving bromine species. These modifications are shown to reduce the overabundance of modelled total column ozone in the Arctic during October to February, southern midlatitudes during August and the Antarctic during September. Antarctic springtime ozone depletion is shown to be enhanced by 25 DU on average, which now causes the ozone hole to be somewhat too deep compared to observations. We show that this is in part due to a cold bias of the Antarctic polar vortex in the model.

1 Introduction
Back to toptop

Stratospheric chemistry is a crucial aspect of chemistry–climate models primarily due to the coupling of ozone with atmospheric dynamics. Ozone strongly absorbs ultraviolet (UV) radiation, thus controlling the temperature of the stratosphere and hence the speed and structure of large-scale stratospheric circulation (e.g. McLandress et al.2010; Braesicke et al.2013; Keeble et al.2014). Ozone depletion in southern high latitudes has been shown to drive trends in the Southern Annular Mode (SAM) which is the leading mode of climate variability in the Southern Hemisphere (SH) (e.g. McLandress et al.2011; Thompson et al.2011; Dennison et al.2015). In turn, the SAM affects numerous other climate features, e.g. Antarctic surface temperatures (Thompson and Solomon2002; Marshall2007; Gillett et al.2006), the Southern Ocean storm track (Yin2005), the atmospheric blocking frequency (Dennison et al.2016) and sea ice (Hall and Visbeck2002; Sen Gupta and England2006). In northern high latitudes, chemical ozone loss is less pronounced (e.g. Solomon et al.2014), although there is some evidence of the effect of strong ozone-depletion events on the large-scale circulation, with impacts on midlatitude weather (e.g. Shindell et al.2001; Ivy et al.2017).

The abundance of ozone is affected by solar radiation via photochemical reactions. The dominant mode of ozone production is via the photolysis of oxygen; ozone destruction is also initiated by photolysis. Ozone-depleting reactions involving only oxygen compounds are known as the Chapman cycle (Chapman1930). Additionally, there are photolysis reactions which are part of various ozone-depleting catalytic cycles, for example, the photolysis of chlorine peroxide (Cl2O2) (Crutzen1974; Molina and Rowland1974; McGrath et al.1990). One of the main modes of variation in solar radiation is the 11-year cycle (Solanki et al.2013, and references therein). This cycle is known to have a significant effect on stratospheric ozone (e.g. Zerefos et al.1997; Calisesi and Matthes2006; Tourpali et al.2007; Kuroda et al.2008; Gruzdev2014). Specifically, studies such as Dameris et al. (2006), Steinbrecht et al. (2004) and Keeble et al. (2018) make the point that ozone increases in the early 2000s were caused by the Sun going through its solar maximum and did not constitute evidence of ozone recovery due to reductions in stratospheric halogen. More recent increases in ozone have been shown to be above those expected by the next solar maximum and are therefore indicative of the effect of declining stratospheric chlorine (Keeble et al.2018; Chipperfield et al.2017).

The importance of simulating the solar cycle in chemistry–climate models (CCMs) has been noted by a number of studies (e.g. Egorova et al.2005; Langematz et al.2005; Tourpali et al.2003; Labitzke et al.2002). Of the 20 CCMs participating in the first phase of the Chemistry-Climate Model Initiative (CCMI-1), only three did not consider solar variability (Morgenstern et al.2017). In this paper, we describe the addition of a solar cycle to the photolysis scheme in a newer version of that model and show the impact this has on modelled stratospheric ozone. In the following, we will refer to the model as “UM-UKCA” (Unified Model – United Kingdom Chemistry and Aerosols) which forms part of the United Kingdom Earth System Model (UKESM) and other configurations of the Met Office Unified Model (UM).

In addition to photolysis reactions, heterogeneous reactions are also important for stratospheric ozone. Halogen species such as hydrogen chloride (HCl), chlorine nitrate (ClONO2), hydrogen bromide (HBr) and bromine nitrate (BrONO2) react on the surfaces of stratospheric aerosols, e.g. polar stratospheric clouds (PSCs), forming species which are subsequently photolysed into ozone-depleting radicals. In this study, two configurations of UKCA are used which differ in the representation of stratospheric chemistry. In previous versions of the UKCA chemistry scheme, only five heterogeneous reactions are considered (Morgenstern et al.2009), none of which involve the activation of bromine. While bromine is much less abundant than chlorine in the stratosphere, on a per-atom basis, it is more effective at depleting ozone (Daniel et al.2007; Newman et al.2007; Sinnhuber et al.2009). As such, it has been shown to have a significant effect in some instances. For example, Sinnhuber et al. (2009) have estimated that approximately half of northern midlatitudes ozone loss over the period 1980–2005 may be attributable to anthropogenic bromine emissions. To account for this, in order to improve the simulation of stratospheric ozone in the UKCA model, eight further heterogeneous reactions involving bromine species are added. In addition, rate constants for the heterogeneous reactions are updated to match latest literature values. Various heterogeneous bromine reactions have been included in other chemistry–climate models (e.g. Jöckel et al.2006; Wegner et al.2013).

This study details the effects on simulated stratospheric ozone of adding solar variability to the UKCA photolysis scheme and of expanding heterogeneous chemistry to include reactions involving bromine compounds. These two additions will be examined separately. Section 2 gives a brief description of the UM-UKCA model. Section 3 deals with the photolysis and Sect. 4 deals with the heterogeneous chemistry. Each section details the changes made to the model, the experiments run to test the impact of these changes, the effect of the changes on the stratospheric ozone and discussion of the results. Section 5 provides a summary.

2 Model description
Back to toptop

The United Kingdom Chemistry and Aerosol (UKCA) module is part of the Met Office Unified Model (UM). Here, we use a configuration with a free-running atmosphere and prescribed sea ice and sea surface temperatures from the Hadley Centre Sea Ice and Sea Surface Temperature (HadISST) dataset (Rayner et al.2003) and operate the model at a resolution of 1.875 longitude by 1.25 latitude with 85 levels extending to 85 km. The UKCA model contains a number of chemistry configurations; here, we use the combined stratosphere and troposphere chemistry (CheST) option (a combination of Morgenstern et al.2009 and O'Connor et al.2014) in a GA7.1 (Walters et al.2017) atmosphere configuration identical to that used in Esentürk et al. (2018) but at Unified Model version 10.7. This configuration includes 75 chemical species and 283 reactions. Further details on the relevant model chemistry are provided in Sects. 3 and 4.

3 Photolysis and the solar cycle
Back to toptop

3.1 Model

Photolysis rates in the configuration of the UKCA model used here are calculated using a combination of the FAST-JX scheme (Wild et al.2000; Bian and Prather2002; Neu et al.2007) and look-up tables. FAST-JX covers wavelengths from 177 to 850 nm over 18 wavelength bins, thus making it suitable for the simulation of chemistry in the stratosphere. FAST-JX calculates scattering for all wavelength bands. The implementation of FAST-JX in UKCA is described by Telford et al. (2012). Above about 60 km, wavelengths shorter then 177 nm become important. In this region, UKCA additionally uses a look-up table of photolysis rates (Lary and Pyle1991; Morgenstern et al.2009) for the contribution to photolysis occurring at these wavelengths. Rates from the two schemes are added together to form the full photolysis rate.

The implementation of the solar cycle in the model uses solar irradiance data from Lean et al. (2005). Using a singular-value decomposition analysis, solar irradiance I (W m−2, which comes as a function of time and wavelength) is expressed as


i.e. a constant term I0 (the temporal average over the time series) and the product of a spectrally varying term I1 and a time series T at monthly resolution. I1 and T are then normalized such that the standard deviation of T equals 1 and I1 is predominantly positive. These terms are illustrated in Fig. 1. The decomposition captures the full variability in the Lean et al. (2005) data extremely well. The spectral term I1 is processed by the FAST-JX binning algorithm into an 18-element array (Bian and Prather2002) and for the look-up table part of photolysis it is processed into a 46-element array covering the wavelength range of 116.7 to 176.2 nm (the Lean et al.2005 data cover wavelengths of 120 nm and above, so the first five elements are zero). The top-of-the-atmosphere incoming shortwave flux is then modified by these time series and wavelength factors. The time series is extended into the future by repeating an average of the last five solar cycles.

To test this modification, two 20-year runs were produced: one with the solar cycle switched on and a control run with the solar cycle switched off. In order to examine only the direct effect of the modification, feedbacks between the chemistry (ozone, methane, nitrous oxide and water vapour) and radiation have been switched off for these runs (i.e. these runs are identical in terms of their meteorological variables). These runs use historic greenhouse gas (GHG) and ozone-depleting substance (ODS) forcing. The effect of the solar cycle was tested by comparing the oxygen photolysis rates in the two runs. This reaction is a component of the cycle which is the main mode of ozone production in the stratosphere and requires radiation with wavelengths shorter than 240 nm (Chapman1930). Figure 1 shows that at these wavelengths the solar irradiance will vary with an amplitude of up to 10 % for a 1-standard-deviation anomaly of solar output I, so we expect an effect on the oxygen loss rate of a similar magnitude.

Figure 1Singular-value decomposition of the solar irradiance data into the (a) time series and (b) spectral components.


3.2 Results

Figure 2a shows the anomalous globally averaged oxygen loss rate in the solar cycle run relative to the control run. The 11-year cycle is clearly evident with the amplitude of approximately ±10 % in the mesosphere. This magnitude is what is expected given the spectral component shown in Fig. 1 is of the order of 10 % for the shorter wavelengths at which this reaction occurs. The amplitude of the anomaly decreases in the stratosphere as this short-wavelength radiation gets attenuated.

The effect of the solar cycle implementation on ozone is shown in Fig. 2b. The percentage difference in the global mean total column ozone (TCO) in the solar cycle run relative to the control run is indicated by the red curve. It can be seen that ozone varies in phase with the solar cycle, illustrated here by the black curve (taken from Fig. 1a), with an amplitude of around 1 % or 2–3 Dobson units (DU).

Figure 2(a) The global mean anomaly in the chemical flux of the oxygen photolysis reaction and (b) the mean total column ozone anomaly 60 S–60 N (cyan, and 3-year running mean in blue), both expressed as percentage difference between the control and solar cycle runs. Overlaid on panel (b) is the solar coefficient time series, T (black).


In addition to ozone production via oxygen photolysis, variations in solar flux also affect ozone loss pathways. These include the photolysis of ozone, a component of the Chapman cycle, and various catalytic loss cycles. We illustrate the effect of the solar cycle in Fig. 3a, which shows the difference in mean chemical fluxes of the ozone production (blue) and various loss pathways as functions of altitude for the year 1989 (i.e. near a solar maximum). The figure also shows the net effect from the changes in the production and loss (black curve). The photolysis of ozone extends to longer wavelengths than that of oxygen (Chapman1930) and hence is less influenced by the solar cycle than the photolysis of oxygen (at solar maximum, the photolysis rate of ozone increases by around 1 %). However, the various changes in the fluxes of the catalytic loss cycles offset the increase in production to a large extent. While the effect on ozone production peaks around 44 km, Fig. 3b shows that the impact on the ozone mixing ratio peaks at a somewhat lower altitude, around 39 km. This is due, in large part, to the HO2 loss cycle that increases strongly above 40 km. The effect on the ClO cycle and the loss branch of the Chapman cycle (i.e. the O3 + O3P reaction) also peak slightly above 40 km.

Figure 3Difference in the (a) chemical flux of the ozone production and loss reactions and (b) global mean ozone mixing ratio between the control and solar cycle runs for 1989.


3.3 Discussion

The magnitude of the solar-cycle-driven ozone effect shown here, i.e. 1 %–2 % between solar maximum and solar minimum, compares well to other observational and climate model studies. A review of a number of observational studies by Calisesi and Matthes (2006) found a magnitude around 2 % for stratospheric ozone. Reinsel et al. (2002), using total column ozone data from the Total Ozone Mapping Spectrometer (TOMS) satellite instruments, found the effect to be 1 %–2 %. More recently, Maycock et al. (2016) found the effect in the lower stratosphere to be around 1 % using the Stratospheric Aerosol and Gas Experiment (SAGE) and Solar Backscatter Ultraviolet Merged Ozone Dataset (SBUVMOD) observations. Maycock et al. (2018) examined a number of CCMI models and found them to be largely consistent with these observational findings.

Figure 2 shows the variation in ozone to be strongly in phase with the solar cycle. Previous work has found there to be phase shift between the solar cycle and its effect on ozone. For example, Gruzdev (2014) found that in the middle and lower stratosphere the oscillation of ozone leads the solar cycle, while in the upper stratosphere it lags the solar cycle, with the difference being as much as a quarter cycle. It is proposed that these results follow from the variation in the strength of the Brewer–Dobson circulation with the solar cycle (Kodera2006). Angell (1989) also found that ozone leads the solar cycle but with a much smaller phase difference of around half a year. Because the model runs used here have chemistry–radiation feedback switched off, we would not expect to see this effect.

Bednarz et al. (2018) also examine an implementation of a solar cycle in the UKCA model. They take a slightly different approach wherein the amplitude of the solar cycle is defined by the difference between the 1981 solar maximum and the 1986 solar minimum. Observational data are used to partition this variation across the 18 FAST-JX wavelength bins as well as the UM shortwave radiation scheme's six bins. No spectral variation in the solar irradiance is accounted for in the look-up table. The results presented here are broadly similar to those shown by Bednarz et al. (2018); they report a solar-driven effect on total column ozone in the 60 S–60 N region of 4.1 (±1.9) and 5.6 (±1.4) DU/W m−2 using two different methods, where the difference in total irradiance between solar maximum and minimum was 1.06 W m−2. The difference between solar maximum and minimum shown in Fig. 2 expressed as an absolute difference is around 4–5 DU. Bednarz et al. (2018) also find the ozone response to increase away from the Equator but do not report values for the polar regions.

The approach we have taken here has the advantage of using historical solar data. Figure 1a shows that the amplitude of the cycle is not constant over the historical record, and as such historical simulations should be more accurate than if using a cycle of constant amplitude. Also, our more comprehensive treatment of the look-up table component, which is used here for pressure levels above 20 Pa, is likely to produce better results in the upper model levels, for example, in the photolysis of oxygen shown in Fig. 2a. Our solar cycle implementation is solely in the UKCA module, unlike Bednarz et al. (2018) who also implemented solar variability in the UM shortwave radiation scheme (which is separate from the photolysis scheme). It may be possible to extend our treatment of solar variability to the UM radiation scheme although it is unlikely to have a large effect given the smaller amplitude of the solar cycle at the relevant wavelengths. Specifically, Bednarz et al. (2018) find that the difference between solar maximum and minimum in the shortest UM wavelength band (covering 200–320 nm) is only 0.56 % and is less than 0.1 % in the longer-wavelength bands.

4 Heterogeneous chemistry
Back to toptop

4.1 Model

The pre-existing UKCA stratospheric scheme included five heterogeneous reactions, the first five listed in Table 1. Reaction rates are calculated based on uptake on three classes of aerosols: ice, nitric acid trihydrate (NAT) and sulfate aerosol (SA). Ice abundance is taken from the UM microphysics scheme (Wilson and Ballard1999) (i.e. it is consistent with cloud physics used elsewhere in the model); the calculation of NAT abundances follows Hanson and Mauersberger (1988), and a climatology of sulfate aerosol based on the year 2000 is prescribed (Morgenstern et al.2010). Table 1 lists the uptake coefficient on each aerosol type for each reaction; the uptake coefficient is the probability of the gas molecule undergoing an irreversible reaction upon collision with the specified surface. Uptake coefficients for these reactions are those recommended by the Jet Propulsion Laboratory (JPL) (Burkholder et al.2015); in many cases, this represents a change from the coefficients in a pre-existing scheme (which are listed under “old” in Table 1).

Table 1Heterogeneous reactions included in the updated UKCA module and the uptake coefficients on ice, NAT and SA.

Uptake coefficients are denoted f where they are not universal constants; see the references for the full formulation. a Burkholder et al. (2015). b Shi et al. (2001) (as recommended by both JPL and IUPAC). c This was retained from the original because the JPL recommendation, which is a function of temperature, evaluates to approximately 0.1 for typical stratospheric values. d Crowley et al. (2010). e Steele and Hamill (1981), Zhang et al. (1994), Cox et al. (1994). f Zhang et al. (1994).

Download Print Version | Download XLSX

We extend the heterogeneous chemistry scheme by adding eight reactions involving bromine species; these are also listed in Table 1. The uptake coefficients for the additional reactions follow International Union of Pure and Applied Chemistry (IUPAC) (Crowley et al.2010) recommendations, except for the BrONO2+HCl reaction on SA surfaces which follows Burkholder et al. (2015).

In Fig. 4, the changes to the heterogeneous chemistry are initially tested in separate parts to assess their relative importance, namely the updating of the NAT and ice uptake coefficients for the original five reactions (labelled “new coeff.”), the new formulation for the uptake coefficients on sulfate aerosol (“Shi 2001”) and the addition of the eight new bromine reactions (“new br”), as well as these modifications in combination (“all”). The 15-month simulations are produced for each and compared to a control simulation with the original heterogeneous chemistry. Similar to the solar cycle runs, chemistry–radiation feedbacks are switched off for these runs in order to examine only the direct effect of the modifications.

Additionally, a 21-year run (1989–2009) for both the original and new configurations is produced, with chemistry–radiation feedbacks switched on, to provide a more reliable indication of the impact on stratospheric chemistry. These runs are forced by the IPCC “historic” and RCP6.0 greenhouse gas scenario (Masui et al.2011; Meinshausen et al.2011; van Vuuren et al.2011) and the A1 scenario for ozone-depleting substances (WMO2011). We compare these runs to the National Institute of Water and Atmospheric Research – Bodeker Scientific (NIWA-BS) total column ozone dataset (version 3.4; see, last access: 12 October 2018). The time span we use for comparison is 1996–2009, excluding 2002. The beginning of the run is not used to ensure any transient effects of the chemistry changes are discounted. The year 2002 is not included because the stratospheric sudden warming that occurred in that year over Antarctica resulted in anomalously low polar ozone in the observations.

4.2 Results

The largest effect of the changes to the heterogeneous chemistry is found at southern polar latitudes. Figure 4a shows the differences in total column ozone, averaged over the region south of 65 S, between the various sensitivity runs and the control run. Updating the NAT and ice uptake coefficients of the original five reactions (yellow curve) increases ozone during spring. This is expected given that the coefficients were either unchanged or decreased from the original version, effectively decreasing chlorine activation. The magnitude of this change is up to 9 DU and occurs in late September. This change in the NAT and ice coefficients is outweighed by adopting the Shi et al. (2001) formulation for the uptake coefficients on sulfate aerosol (blue curve). The effect of this change is an enhancement of the springtime ozone depletion, peaking at 20 DU in early October. In contrast, the effect of the additional bromine reactions (green curve) is evident from the beginning of the run and amounts to around 10 DU in late winter, before the pronounced springtime ozone depletion begins. The additional bromine reactions are at maximum effect in early spring, when they amount to 20 DU of additional ozone depletion. These three changes are combined in the run labelled “all” (red curve), which shows the maximum effect on total column ozone to be a decrease in ozone of 30 DU in late September/early October, but due to the inclusion of the bromine reactions is important in all seasons. The effect from the combination of the three modifications is slightly more than the sum of the individual runs. This effect is confined mostly to September and is up to about 4 DU. Figure 4b shows the equivalent for the Northern Hemisphere. The overall impact on ozone is smaller (at maximum 15 DU during March) and it is notable that the addition of bromine is producing the majority of the effect relative to the either of the modifications to the original reactions.

Figure 4The difference in mean total column ozone (a) south of 65 S and (b) north of 65 N, between the various sensitivity runs and the control run.


For the rest of this section, we focus on the fully interactive simulations, comparing the new heterogeneous chemistry (“new het.”) to the control simulation and to the NIWA-Bodeker total column ozone dataset. Figure 5 shows the zonal mean total column ozone climatology in each of the model runs. The figure uses crosshatching and stippling to illustrate areas where the model underestimates or overestimates, respectively, the NIWA-BS dataset by more than 30 DU. There are no NIWA-BS data over the winter poles, these areas are marked by the dashed contours. Comparing the two runs, it can be seen that the new heterogeneous chemistry reduces biases in a number of areas: the northern high latitudes during October to February, the equatorial region during August to October, the southern midlatitudes during August and the southern polar region during the September onset of the ozone hole. The one area where the changes to the heterogeneous chemistry make the model notably worse is the southern polar region during austral summer, when the existing underestimation of ozone by the model is exacerbated. This can be seen more clearly in Fig. 6, which shows the polar mean (south of 65 S) total column ozone climatology. During springtime, peak ozone depletion is enhanced by around 25 DU (note: this is slightly smaller than that indicated in Fig. 4, likely reflecting that the 1999 forcings during the sensitivity runs were more conducive to ozone depletion than the 1996–2009 average). This increased depletion persists with a similar magnitude through summer. The effect of the new configuration is that the extent of peak ozone depletion is now overestimated relative to the NIWA-BS dataset, whereas for the control run it was underestimated. An issue with the control run illustrated here is the slow summertime replenishing of ozone. Although “new het.” replenishes ozone at a similar rate, the increase in the depth of the ozone hole means that the summertime discrepancy in ozone increases significantly relative to NIWA-BS.

Figure 5Zonal mean total column ozone climatology (1996–2009, excluding 2002) of the (a) control and (b) “new het.” runs. Contours show the NIWA-BS data; crosshatching/stippling indicates where the model under-/overestimates NIWA-BS data by more than 30 DU. Dashed contours mark the absence of NIWA-BS data.


Figure 6Southern polar mean total column ozone climatology of the “new het.” (green) and control (red) runs as well as the NIWA-BS data (black). The shaded bands show ±1 standard deviation; uncertainty in the control run is not shown here for clarity. The average is over the region south of 65 S and for the years 1996–2009 (excluding 2002).


The overestimation of springtime ozone depletion shown in Fig. 5 is linked to model temperature biases which influence polar ozone depletion. We illustrate this in Fig. 7, which shows the relationship between polar total column ozone and polar lower stratospheric mean temperature (south of 65 S and pressure levels between 100 and 10 hPa) for October in the model runs and the NIWA-BS dataset (with temperatures from the ERA-Interim reanalysis; Dee et al.2011). The model has a substantial cold bias relative to the ERA-Interim temperature, which is not solely due to differences in ozone. Based on the gradient of the ozone–temperature relationship shown here, a correction of the temperature bias would reveal a superior simulation of the ozone depletion by the “new het.” configuration.

Figure 7Southern polar mean total column ozone versus the mean temperature in the polar lower stratosphere for the “new het.” (green) and control (red) runs as well as NIWA-BS/ERA-Interim (black). The average is taken over the region south of 65 S; temperature is additionally averaged over the 100–10 hPa pressure range. Data cover the years 1996–2009 (excluding 2002).


Figure 8 shows the difference in zonal mean ozone mixing ratio during October. Student's t test is used to identify regions where the model runs differ significantly. The enhanced polar ozone depletion is primarily at the lower reaches of the ozone layer between around 12 and 24 km in altitude. Figure 9 shows the change in chemical flux of the various ozone production/destruction pathways in this southern polar lower stratosphere region. The largest change is via the Cl2O2 photolysis cycle. This change represents an 8 % increase over the flux in the control run.

Figure 8Difference in zonal mean ozone between the “new het.” and control runs during October 1996–2009 (shading). Contours illustrate the control run climatology. Stippling indicates the difference is significant according to Student's t test (p<0.01).


Figure 9Annual difference in ozone production (red)/destruction (blue) flux between “new het.” and control runs, in the lower stratosphere (12–24 km) south of 65 S (1996–2009).


Table 2 lists the differences in flux of the various heterogeneous reactions between the control and “new het.” runs. The changes in the original five reactions include an increase in the rate of the ClONO2+HCl reaction at the expense of the hydrolysis reaction. This is almost entirely due to the change in the calculation of the uptake coefficient on sulfate aerosol. Also, the rate of N2O5 hydrolysis decreases; in this case, the uptake rates on both ice and NAT have been lowered by one-third (while the reaction on sulfate is unchanged). The most important of the new bromine reactions are BrONO2+HCl and BrONO2+H2O, the two bromine reactions for which reaction on sulfate aerosols is considered.

Table 2Difference in heterogeneous reaction flux between “new het.” and control runs (1996–2009).

Download Print Version | Download XLSX

Figure 10a shows the October SH polar mean vertical profile of the bromine species: BrO, HBr and BrONO2 in the “new het.” (solid curves) and control (dashed) runs. The effect on BrONO2, suggested by the fluxes in Table 2, is shown to be quite substantial here. The lower peak in BrONO2 at around 12 km is strongly reduced, with a co-located decrease in BrO and increase in HBr localized to the same altitude range. Figure 10b shows the equivalent to Fig. 10a for the chlorine species: ClO, Cl2O2, HCl and ClONO2. The impact on the chlorine species varies over a larger altitude range. ClONO2 increases broadly over both peaks around 19 and 26 km; this is primarily due to the decrease uptake on NAT. Below about 15 km, where SA is more prominent, ClONO2 is reduced due to the new loss pathway of reaction with HCl on SA.

Figure 10Vertical profile of (a) the species: BrO, HBr and BrONO2 and (b) ClO, Cl2O2, HCl and ClONO2 for October, averaged over 1996–2009 and south of 65 S. “New het.” is shown with solid curves and the control run with dashed curves. Shading indicates ±1 standard deviation around the “new het.” mean.


4.3 Discussion

The changes to the heterogeneous chemistry result in increased SH springtime polar ozone depletion. This has the advantage of more accurately simulating the timing and rate of decline of ozone in early spring. During October, the peak ozone depletion is slightly overestimated relative to the NIWA-BS data; however, this is likely due to the model's cold bias (see Fig. 7). The downside of this change occurs during summer, when the model has a tendency to replenish the SH polar ozone too slowly. It is possible the increased ozone depletion itself exacerbates this problem, as it has been shown that increased ozone depletion acts to delay the break-up of the polar vortex (Haigh and Roscoe2009; McLandress et al.2010; Keeble et al.2014). This summertime ozone problem could be both dynamical, for example, a lack of planetary wave activity which delays the vortex break-up, or chemical, for example, in the conversion of radicals back into reservoir species. This remains an important area for model development.

Another area of possible model development is the treatment of PSCs in the UKCA. In the case of NAT, particles are assumed to form when HNO3 reaches saturation (following Hanson and Mauersberger1988). More recent work has found that NAT formation does not occur at the equilibrium temperature predicted by this assumption but at temperatures around 3 K lower (Wegner et al.2013). This would suggest that a HNO3 supersaturation of around 10 ought to be required for NAT formation. It is possible that this may not be very important given that the role of NAT is somewhat de-emphasized in our new heterogeneous chemistry scheme, with the corresponding uptake coefficients on NAT lowered in four of the five original reactions and not considered for six of the additional eight reactions. Perhaps of more consequence is the treatment of sulfate aerosol. UKCA considers sulfate aerosol; however, sulfate aerosol absorbs HNO3 at lower temperatures, forming what is known as a supercooled ternary solution (STS) (Hamill et al.1996). Data regarding reactions on STS are more limited, but it has been shown to be important in some cases. For example, Zhang et al. (1995) found that the rate of N2O5 hydrolysis on STS to decrease as the amount of HNO3 in solution is increased, while the rate of ClONO2 hydrolysis is unaffected by the presence of HNO3.

A hazard of modelling stratospheric chemistry is the limited experimental data on which parameterizations are based. Of the eight reactions added to the heterogeneous chemistry scheme, the most active of these are the reactions of BrONO2 with HCl and H2O on sulfate aerosol. The experimental work on these reactions (Hanson and Ravishankara1995; Hanson et al.1996; Hanson2003) focused mostly on the hydrolysis reaction. The uptake coefficient of 0.9 for BrONO2+HCl, recommend by JPL, is based upon just two experiments (Hanson et al.1996). One tests 60 wt % H2SO4 with a small amount of HCl (approximately 10−3 M), and the other 48 wt % with 0.3 M HCl, both at 229 K, finding uptake coefficients of 0.9±0.2 and 1.0±0.2, respectively. Hence, there may be substantial uncertainty associated with this value. Perhaps because of this uncertainty, this reaction is not included in some chemistry–climate models; for example, the Whole Atmosphere Community Climate Model (WACCM) excludes this reaction while including the BrONO2+H2O reaction (Solomon et al.2015). However, a coefficient as large as 0.9 is not necessarily surprising given that the parameterization of BrONO2+H2O (which is based on much more experimental data) gives an uptake coefficient of 0.8 for solutions with H2SO4 at less than 65 wt %. As Table 2 shows BrONO2+HCl to be relatively important, we choose to include it despite the associated uncertainty.

Compared to the BrONO2 reactions with H2O and HCl, the other additional reactions are shown to have relatively little impact. However, we include these in the UKCA heterogeneous chemistry framework for completeness, as it is possible that future work may reveal them to be more relevant. In particular, further investigation of these reactions on sulfate aerosols would be useful given that the uptake of BrONO2 is known to be so large.

The sulfate aerosol used in this work is a climatology representing the “background” level, so the effect on ozone demonstrated here could be perhaps thought of as a lower limit. Further work looking at the sensitivity of this heterogeneous chemistry scheme to the effect of elevated aerosol levels (e.g. from volcanic eruptions) would also be of interest.

Given the importance of bromine activation on sulfate aerosols, another interesting aspect is the sensitivity of ozone depletion to the abundance of bromine. Yang et al. (2014) examined this, also using the UM-UKCA model, finding a 5 % decrease in springtime SH polar TCO (∼10 DU) in response to an increase in stratospheric inorganic bromine resulting from a doubling of the very short-lived species (VSLS) source. The uptake coefficients used by Yang et al. (2014) differ from those used here; this configuration shows a larger impact on ozone. It is therefore possible that the ozone-depletion sensitivity demonstrated by Yang et al. (2014) is somewhat of an underestimation. Yang et al. (2014) also found the ozone-depletion sensitivity to bromine concentration to be dependent on the amount of chlorine present. Thus, a deeper study of the role of the chlorine/bromine cross reactions would be an interesting area for further study with this new configuration of heterogeneous chemistry.

5 Summary
Back to toptop

Improvements to the UKCA model have been made in two areas: the photolysis and the heterogeneous chemistry schemes. The photolysis scheme was improved by the addition of solar variability which modifies the solar flux entering the photolysis rate calculations. The solar cycle is based on observations and is spectrally resolved. At short wavelengths, the amplitude of the solar cycle is of the order of 10 %. Consequently, reactions (such as the photolysis of oxygen) that require this short-wavelength radiation vary proportionately. We show that the effect of the added solar cycle on total column ozone in the extratropical regions is of the order of 1 %, in phase with the solar cycle. The effect is likely to be larger in polar regions although the model run in this study is not of sufficient length to provide reliable estimates given the increased variability of ozone in these regions.

The second improvement is to the heterogeneous chemistry scheme. This consisted of updating the uptake coefficients of the five pre-existing reactions as well as adding a further eight reactions. The impact on ozone from these changes was substantial, amounting to around 25 DU of additional springtime SH polar ozone depletion, with significant effect coming from both the updated original reactions and the new reactions. Comparison with the NIWA-BS ozone dataset shows that the new heterogeneous chemistry scheme simulates springtime ozone better, especially when accounting for the model temperature bias, although the increase in springtime ozone depletion leads to an increase in the low-ozone bias during summer.

Code and data availability
Back to toptop
Code and data availability. 

The Met Office Unified Model is available for use under licence. For information on how to apply for a licence, see (Met Office2018). Both developments are available in the UM truck from vn11.2. NIWA-BS total column ozone data (vn3.4) were obtained from (retrieved 12 October 2018). ERA-Interim data were obtained from (retrieved 16 May 2018). Data from the UM-UKCA model runs and the code used to produce the figures in this paper can be found at (Dennison et al.2018).

Author contributions
Back to toptop
Author contributions. 

FD wrote the paper with input from all authors. The solar cycle code was originally written by OM, with modifications by FD. The heterogeneous chemistry code was written by JK and FD with assistance from all authors. Model runs and analysis were produced by FD with advice from OM and GZ.

Competing interests
Back to toptop
Competing interests. 

The authors declare that they have no conflict of interest.

Back to toptop

This work has been funded by the New Zealand Government Ministry for Business, Innovation, and Employment (MBIE) through the Deep South National Science Challenge. This work has also been supported by NIWA as part of its government-funded core research. Also, funding was received from European Community's Seventh Framework Programme (FP7/2007-2013) under grant agreement no. 603557 (StratoClim). Xin Yang acknowledges financial support from NERC/UK through the project NE/M015394/1. The authors wish to acknowledge the contribution of NeSI high-performance computing facilities to the results of this research. New Zealand's national facilities are provided by the New Zealand eScience Infrastructure (NeSI) and funded jointly by NeSI's collaborator institutions and through MBIE's Research Infrastructure programme ( We would also like to thank Bodeker Scientific, funded by the New Zealand Deep South National Science Challenge, for providing the combined NIWA-BS total column ozone database.

Review statement
Back to toptop
Review statement. 

This paper was edited by Volker Grewe and reviewed by two anonymous referees.

Back to toptop

Angell, J. K.: On the Relation between Atmospheric Ozone and Sunspot Number, J. Climate, 2, 1404–1416,<1404:OTRBAO>2.0.CO;2, 1989. a

Bednarz, E. M., Maycock, A. C., Telford, P. J., Braesicke, P., Abraham, N. L., and Pyle, J. A.: Simulating the atmospheric response to the 11-year solar cycle forcing with the UM-UKCA model: the role of detection method and natural variability, Atmos. Chem. Phys. Discuss.,, in review, 2018. a, b, c, d, e

Bian, H. and Prather, M. J.: Fast-J2: Accurate Simulation of Stratospheric Photolysis in Global Chemical Models, J. Atmos. Chem., 41, 281–296, 2002. a, b

Braesicke, P., Keeble, J., Yang, X., Stiller, G., Kellmann, S., Abraham, N. L., Archibald, A., Telford, P., and Pyle, J. A.: Circulation anomalies in the Southern Hemisphere and ozone changes, Atmos. Chem. Phys., 13, 10677–10688,, 2013. a

Burkholder, J. B., Sander, S. P., Abbatt, J. P. D., 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 Evaluation Number 18, Tech. rep., Jet Propulsion Laboratory, Pasadena, California, available at: (last access: 19 April 2018), 2015. a, b, c

Calisesi, Y. and Matthes, K.: The middle atmospheric ozone response to the 11-year solar cycle, Space Sci. Rev., 125, 273–286,, 2006. a, b

Chapman, S.: On Ozone and Atomic Oxygen in the Upper Atmosphere, Philos. Mag., 10, 369–383, 1930. a, b, c

Chipperfield, M. P., Bekki, S., Dhomse, S., Harris, N. R. P., Hassler, B., Hossaini, R., Steinbrecht, W., Thiéblemont, R., and Weber, M.: Detecting recovery of the stratospheric ozone layer, Nature, 549, 211–218,, 2017. a

Cox, R. A., MacKenzie, A. R., Muller, R. H., Peter, T., and Crutzen, P. J.: Activation of stratospheric chlorine by reactions in liquid sulphuric acid, Geophys. Res. Lett., 22, 1439–1442,, 1994. a

Crowley, J. N., Ammann, M., Cox, R. A., Hynes, R. G., Jenkin, M. E., Mellouki, A., Rossi, M. J., Troe, J., and Wallington, T. J.: Evaluated kinetic and photochemical data for atmospheric chemistry: Volume V – heterogeneous reactions on solid substrates, Atmos. Chem. Phys., 10, 9059–9223,, 2010. a, b

Crutzen, P. J.: Estimates of possible future ozone reductions from continued use of fluoro-chloro-methanes (CF2Cl2, CFCl3), Geophys. Res. Lett., 1, 205–208,, 1974. a

Dameris, M., Matthes, S., Deckert, R., Grewe, V., and Ponater, M.: Solar cycle effect delays onset of ozone recovery, Geophys. Res. Lett., 33, L03806,, 2006. a

Daniel, J. S., Velders, G. J. M., Douglass, A., Forster, P., Hauglustaine, D., Isaksen, I., Kuijpers, L., McCulloch, A., and Wallington, T.: Halocarbon Scenarios, Ozone Depletion Potentials, and Global Warming Potentials, chap. 8 in: Scientific assessment of ozone depletion: 2006, Global Ozone Research and Monitoring Project – Report No. 50, Tech. rep., World Meteorological Organization, Geneva, 2007. a

Dee, D. P., Uppala, S. M., Simmons, A. J., Berrisford, P., Poli, P., Kobayashi, S., Andrae, U., Balmaseda, M. A., Balsamo, G., Bauer, P., Bechtold, P., Beljaars, A. C. M., van de Berg, L., Bidlot, J., Bormann, N., Delsol, C., Dragani, R., Fuentes, M., Geer, A. J., Haimberger, L., Healy, S. B., Hersbach, H., Hólm, E. V., Isaksen, L., Kållberg, P., Köhler, M., Matricardi, M., Mcnally, A. P., Monge-Sanz, B. M., Morcrette, J. J., Park, B. K., Peubey, C., de Rosnay, P., Tavolato, C., Thépaut, J. N., and Vitart, F.: The ERA-Interim reanalysis: Configuration and performance of the data assimilation system, Q. J. Roy. Meteor. Soc., 137, 553–597,, 2011. a

Dennison, F., Keeble, J., Morgenstern, O., Zeng, G., Abraham, N. L., and Yang, X.: Improvements to stratospheric chemistry scheme in the UM-UKCA (v10.7) model: solar cycle and heterogeneous reactions [Data set], Zenodo,, 2018. a

Dennison, F. W., McDonald, A., and Morgenstern, O.: The Influence of Ozone Forcing on Blocking in the Southern Hemisphere, J. Geophys. Res.-Atmos., 121, 14358–14371,, 2016. a

Dennison, F. W., McDonald, A. J., and Morgenstern, O.: The effect of ozone depletion on the Southern Annular Mode and stratosphere–troposphere coupling, J. Geophys. Res.-Atmos., 120, 6305–6312,, 2015. a

Egorova, T., Rozanov, E., Zubov, V., Schmutz, W., and Peter, T.: Influence of solar 11-year variability on chemical composition of the stratosphere and mesosphere simulated with a chemistry-climate model, Adv. Space Res., 35, 451–457,, 2005. a

Esentürk, E., Abraham, N. L., Archer-Nicholls, S., Mitsakou, C., Griffiths, P., Archibald, A., and Pyle, J.: Quasi-Newton methods for atmospheric chemistry simulations: implementation in UKCA UM vn10.8, Geosci. Model Dev., 11, 3089–3108,, 2018. a

Gillett, N. P., Kell, T. D., and Jones, P. D.: Regional climate impacts of the Southern Annular Mode, Geophys. Res. Lett., 33, L23704,, 2006. a

Gruzdev, A. N.: Estimate of the effect of the 11-year solar activity cycle on the ozone content in the stratosphere, Geomagn. Aeronomy, 54, 633–639, 2014. a, b

Haigh, J. D. and Roscoe, H. K.: The final warming date of the Antarctic polar vortex and influences on its interannual variability, J. Climate, 22, 5809–5819,, 2009. a

Hall, A. and Visbeck, M.: Synchronous variability in the Southern Hemisphere atmosphere, sea ice, and ocean resulting from the annular mode, J. Climate, 15, 3043–3057,<3043:SVITSH>2.0.CO;2, 2002. a

Hamill, P., Tabazadeh, A., Kinne, S., Toon, O. B., and Turco, R. P.: On the growth of ternary system HNO3/H2SO4/H2O aerosol particles in the stratosphere, Geophys. Res. Lett., 23, 753–756, 1996. a

Hanson, D. R.: Reactivity of BrONO2 and HOBr on sulfuric acid solutions at low temperatures, J. Geophys. Res., 108, 4239,, 2003. a

Hanson, D. R. and Mauersberger, K.: Laboratory studies of the nitric acid trihydrate: Implications for the south polar stratosphere, Geophys. Res. Lett., 15, 855–858,, 1988. a, b

Hanson, D. R. and Ravishankara, A. R.: Heterogeneous chemistry of bromine species in sulfuric acid under stratospheric conditions, Geophys. Res. Lett., 22, 385–388, 1995. a

Hanson, D. R., Ravishankara, A. R., and Lovejoy, E. R.: Reaction of BrONO2 with H2O on submicron sulfuric acid aerosol and the implications for the lower stratosphere, J. Geophys. Res., 101, 9063–9069,, 1996. a, b

Ivy, D. J., Solomon, S., Calvo, N., and Thompson, D. W.: Observed connections of Arctic stratospheric ozone extremes to Northern Hemisphere surface climate, Environ. Res. Lett., 12, 024004,, 2017. a

Jöckel, P., Tost, H., Pozzer, A., Brühl, C., Buchholz, J., Ganzeveld, L., Hoor, P., Kerkweg, A., Lawrence, M. G., Sander, R., Steil, B., Stiller, G., Tanarhte, M., Taraborrelli, D., van Aardenne, J., and Lelieveld, J.: The atmospheric chemistry general circulation model ECHAM5/MESSy1: consistent simulation of ozone from the surface to the mesosphere, Atmos. Chem. Phys., 6, 5067–5104,, 2006. a

Keeble, J., Braesicke, P., Abraham, N. L., Roscoe, H. K., and Pyle, J. A.: The impact of polar stratospheric ozone loss on Southern Hemisphere stratospheric circulation and climate, Atmos. Chem. Phys., 14, 13705–13717,, 2014. a, b

Keeble, J., Brown, H., Abraham, N. L., Harris, N. R. P., and Pyle, J. A.: On ozone trend detection: using coupled chemistry–climate simulations to investigate early signs of total column ozone recovery, Atmos. Chem. Phys., 18, 7625–7637,, 2018. a, b

Kodera, K.: The role of dynamics in solar forcing, Space Science Reviews, 125, 319–330,, 2006. a

Kuroda, Y., Yamazaki, K., and Shibata, K.: Role of ozone in the solar cycle modulation of the North Atlantic Oscillation, J. Geophys. Res., 113, D14122,, 2008. a

Labitzke, K., Austin, J., Butchart, N., Knight, J., Takahashi, M., Nakamoto, M., Nagashima, T., Haigh, J., and Williams, V.: The global signal of the 11-year solar cycle in the stratosphere: Observations and models, J. Atmos. Sol.-Terr. Phy., 64, 203–210,, 2002. a

Langematz, U., Grenfell, J. L., Matthes, K., Mieth, P., Kunze, M., Steil, B., and Brühl, C.: Chemical effects in 11-year solar cycle simulations with the Freie Universität Berlin Climate Middle Atmosphere Model with online chemistry (FUB-CMAM-CHEM), Geophys. Res. Lett., 32, L13803,, 2005. a

Lary, D. J. and Pyle, J. A.: Diffuse radiation, twilight, and photochemistry – II, J. Atmos. Chem., 13, 393–406,, 1991. a

Lean, J., Rottman, G., Harder, J., and Kopp, G.: SORCE Contributions to New Understanding of Global Change and Solar Variability, Sol. Phys., 230, 27–53,, 2005. a, b, c

Marshall, G. J.: Half-century seasonal relationships between the Southern Annular Mode and Antarctic temperatures, Int. J. Climatol., 27, 373–383,, 2007. a

Masui, T., Matsumoto, K., Hijioka, Y., Kinoshita, T., Nozawa, T., Ishiwatari, S., Kato, E., Shukla, P. R., Yamagata, Y., and Kainuma, M.: An emission pathway for stabilization at 6 W m−2 radiative forcing, Climatic Change, 109, 59–76,, 2011. a

Maycock, A. C., Matthes, K., Tegtmeier, S., Thiéblemont, R., and Hood, L.: The representation of solar cycle signals in stratospheric ozone – Part 1: A comparison of recently updated satellite observations, Atmos. Chem. Phys., 16, 10021–10043,, 2016. a

Maycock, A. C., Matthes, K., Tegtmeier, S., Schmidt, H., Thiéblemont, R., Hood, L., Akiyoshi, H., Bekki, S., Deushi, M., Jöckel, P., Kirner, O., Kunze, M., Marchand, M., Marsh, D. R., Michou, M., Plummer, D., Revell, L. E., Rozanov, E., Stenke, A., Yamashita, Y., and Yoshida, K.: The representation of solar cycle signals in stratospheric ozone – Part 2: Analysis of global models, Atmos. Chem. Phys., 18, 11323–11343,, 2018. a

McGrath, M. P., Clemitshaw, K. C., Rowland, F. S., and Hehre, W. J.: Structures, relative stabilities, and vibrational spectra of isomers of chlorine oxide dimer (Cl2O2): the role of the chlorine oxide dimer in antarctic ozone depleting mechanisms, J. Phys. Chem., 94, 6126–6132, 1990. a

McLandress, C., Jonsson, A. I., Plummer, D. A., Reader, M. C., Scinocca, J. F., and Shepherd, T. G.: Separating the dynamical effects of climate change and ozone depletion. Part I: Southern hemisphere stratosphere, J. Climate, 23, 5002–5020,, 2010. a, b

McLandress, C., Shepherd, T. G., Scinocca, J. F., Plummer, D. A., Sigmond, M., Jonsson, A. I., and Reader, M. C.: Separating the dynamical effects of climate change and ozone depletion. Part II: Southern Hemisphere troposphere, J. Climate, 24, 1850–1868,, 2011. a

Meinshausen, M., Smith, S. J., Calvin, K., Daniel, J. S., Kainuma, M. L. T., Lamarque, J., Matsumoto, K., Montzka, S. A., Raper, S. C. B., Riahi, K., Thomson, A., Velders, G. J. M., and van Vuuren, D. P. P.: The RCP greenhouse gas concentrations and their extensions from 1765 to 2300, Climatic Change, 109, 213–241,, 2011. a

Met Office: Unified Model (UM), available at:, last access: 30 October 2018. a

Molina, M. J. and Rowland, F. S.: Stratospheric sink for chlorofluoromethanes: chlorine atomc-atalysed destruction of ozone, Nature, 249, 810–812,, 1974. a

Morgenstern, O., Braesicke, P., O'Connor, F. M., Bushell, A. C., Johnson, C. E., Osprey, S. M., and Pyle, J. A.: Evaluation of the new UKCA climate-composition model – Part 1: The stratosphere, Geosci. Model Dev., 2, 43–57,, 2009. a, b, c

Morgenstern, O., Giorgetta, M. A., Shibata, K., Eyring, V., Waugh, D. W., Shepherd, T. G., Akiyoshi, H., Austin, J., Baumgaertner, A. J. G., Bekki, S., Braesicke, P., Brühl, C., Chipperfield, M., Cugnet, D., Dameris, M., Dhomse, S., Frith, S. M., Garny, H., Gettelman, A., Hardiman, S. C., Hegglin, M. I., Jöckel, P., Kinnison, D. E., Lamarque, J. F., Mancini, E., Manzini, E., Marchand, M., Michou, M., Nakamura, T., Nielsen, J. E., Olivié, D., Pitari, G., Plummer, D. A., Rozanov, E., Scinocca, J. F., Smale, D., Teyssèdre, H., Toohey, M., Tian, W., and Yamashita, Y.: Review of the formulation of present-generation stratospheric chemistry-climate models and associated external forcings, J. Geophys. Res., 115, D00M02,, 2010. a

Morgenstern, O., Hegglin, M. I., Rozanov, E., O'Connor, F. M., Abraham, N. L., Akiyoshi, H., Archibald, A. T., Bekki, S., Butchart, N., Chipperfield, M. P., Deushi, M., Dhomse, S. S., Garcia, R. R., Hardiman, S. C., Horowitz, L. W., Jöckel, P., Josse, B., Kinnison, D., Lin, M., Mancini, E., Manyin, M. E., Marchand, M., Marécal, V., Michou, M., Oman, L. D., Pitari, G., Plummer, D. A., Revell, L. E., Saint-Martin, D., Schofield, R., Stenke, A., Stone, K., Sudo, K., Tanaka, T. Y., Tilmes, S., Yamashita, Y., Yoshida, K., and Zeng, G.: Review of the global models used within phase 1 of the Chemistry–Climate Model Initiative (CCMI), Geosci. Model Dev., 10, 639–671,, 2017. a

Neu, J. L., Prather, M. J., and Penner, J. E.: Global atmospheric chemistry: Integrating over fractional cloud cover, J. Geophys. Res., 112, D11306,, 2007. a

Newman, P. A., Daniel, J. S., Waugh, D. W., and Nash, E. R.: A new formulation of equivalent effective stratospheric chlorine (EESC), Atmos. Chem. Phys., 7, 4537–4552,, 2007. a

O'Connor, F. M., Johnson, C. E., Morgenstern, O., Abraham, N. L., Braesicke, P., Dalvi, M., Folberth, G. A., Sanderson, M. G., Telford, P. J., Voulgarakis, A., Young, P. J., Zeng, G., Collins, W. J., and Pyle, J. A.: Evaluation of the new UKCA climate-composition model – Part 2: The Troposphere, Geosci. Model Dev., 7, 41–91,, 2014. a

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., 108, 4407,, 2003. a

Reinsel, G. C., Weatherhead, E. C., Tiao, G. C., Miller, A. J., Nagatani, R. M., Wuebbles, D. J., and Flynn, L. E.: On detection of turnaround and recovery in trend for ozone, J. Geophys. Res., 107, ACH 1-1–ACH 1-12, 2002. a

Sen Gupta, A. and England, M. H.: Coupled Ocean-Atmosphere-Ice Response to Variations in the Southern Annular Mode, J. Climate, 19, 4457–4486, 2006. a

Shi, Q., Jayne, J. T., Kolb, C. E., and Worsnop, D. R.: Kinetic model for reaction of ClONO2 with H2O and HCl and HOCl with HCl in sulfuric acid solutions, J. Geophys. Res., 106, 24259–24274, 2001. a, b

Shindell, D. T., Schmidt, G. A., Miller, R. L., and Rind, D.: Northern Hemisphere winter climate response to greenhous gas, ozone, solar and volcanic forcing, J. Geophys. Res., 106, 7193–7210, 2001. a

Sinnhuber, B.-M., Sheode, N., Sinnhuber, M., Chipperfield, M. P., and Feng, W.: The contribution of anthropogenic bromine emissions to past stratospheric ozone trends: a modelling study, Atmos. Chem. Phys., 9, 2863–2871,, 2009. a, b

Solanki, S. K., Krivova, N. A., and Haigh, J. D.: Solar Irradiance Variability and Climate, Annu. Rev. Astron. Astr., 51, 311–351,, 2013. a

Solomon, S., Haskins, J., Ivy, D. J., and Min, F.: Fundamental differences between Arctic and Antarctic ozone depletion, P. Natl. Acad. Sci. USA, 111, 6220–6225,, 2014. a

Solomon, S., Kinnison, D., Bandoro, J., and Garcia, R.: Simulation of polar ozone depletion: An update, J. Geophys. Res.-Atmos., 120, 7958–7974,, 2015. a

Steele, H. M. and Hamill, P.: Effects of Temperature and Humidity on the Acid-Water Droplets in the Stratosphere, J. Aerosol Sci., 12, 517–528, 1981. a

Steinbrecht, W., Claude, H., and Winkler, P.: Enhanced upper stratospheric ozone: Sign of recovery or solar cycle effect?, J. Geophys. Res., 109, D14306,, 2004. a

Telford, P. J., Abraham, N. L., Archibald, A. T., Braesicke, P., Dalvi, M., Morgenstern, O., O'Connor, F. M., Richards, N. A. D., and Pyle, J. A.: Implementation of the Fast-JX Photolysis scheme (v6.4) into the UKCA component of the MetUM chemistry-climate model (v7.3), Geosci. Model Dev., 6, 161–177,, 2013. a

Thompson, D. W. J. and Solomon, S.: Interpretation of Recent Southern Hemisphere Climate Change, Science, 296, 895–899, 2002. a

Thompson, D. W. J., Solomon, S., Kushner, P. J., England, M. H., Grise, K. M., and Karoly, D. J.: Signatures of the Antarctic ozone hole in Southern Hemisphere surface climate change, Nat. Geosci., 4, 741–749,, 2011. a

Tourpali, K., Schuurmans, C. J. E., van Dorland, R., Steil, B., and Brühl, C.: Stratospheric and tropospheric response to enhanced solar UV radiation: A model study, Geophys. Res. Lett., 30, 1231,, 2003. a

Tourpali, K., Zerefos, C. S., Balis, D. S., and Bais, A. F.: The 11-year solar cycle in stratospheric ozone: Comparison between Umkehr and SBUVv8 and effects on surface erythemal irradiance, J. Geophys. Res., 112, D12306,, 2007. a

van Vuuren, D. P., Edmonds, J., Kainuma, M., Riahi, K., Thomson, A., Hibbard, K., Hurtt, G. C., Kram, T., Krey, V., Lamarque, J.-F., Masui, T., Meinshausen, M., Nakicenovic, N., Smith, S. J., and Rose, S. K.: The representative concentration pathways: an overview, Climatic Change, 109, 5–31,, 2011. a

Walters, D., Baran, A., Boutle, I., Brooks, M., Earnshaw, P., Edwards, J., Furtado, K., Hill, P., Lock, A., Manners, J., Morcrette, C., Mulcahy, J., Sanchez, C., Smith, C., Stratton, R., Tennant, W., Tomassini, L., Van Weverberg, K., Vosper, S., Willett, M., Browse, J., Bushell, A., Dalvi, M., Essery, R., Gedney, N., Hardiman, S., Johnson, B., Johnson, C., Jones, A., Mann, G., Milton, S., Rumbold, H., Sellar, A., Ujiie, M., Whitall, M., Williams, K., and Zerroukat, M.: The Met Office Unified Model Global Atmosphere 7.0/7.1 and JULES Global Land 7.0 configurations, Geosci. Model Dev. Discuss.,, in review, 2017. a

Wegner, T., Kinnison, D. E., Garcia, R. R., Solomon, S., and Journal, M.: Simulation of polar stratospheric clouds in the specified dynamics version of the whole atmosphere community climate model, J. Geophys. Res.-Atmos., 118, 4991–5002, 2013. a, b

Wild, O., Zhu, X., and Prather, J.: Fast-J: accurate simulation of the in- and below cloud photolysis in tropospheric chemical models, J. Atmos. Chem., 37, 245–282, 2000. a

Wilson, D. R. and Ballard, S. P.: A microphysically based precipitation scheme for the UK Meteorological Office Unified Model, Q. J. Roy. Meteor. Soc., 125, 1607–1636, 1999. a

WMO: Scientific Assessment of Ozone Depletion: 2010, Tech. Rep. 52, World Meteorological Organization, Geneva, Switzerland, 2011. a

Yang, X., Abraham, N. L., Archibald, A. T., Braesicke, P., Keeble, J., Telford, P. J., Warwick, N. J., and Pyle, J. A.: How sensitive is the recovery of stratospheric ozone to changes in concentrations of very short-lived bromocarbons?, Atmos. Chem. Phys., 14, 10431–10438,, 2014. a, b, c, d

Yin, J. H.: A consistent poleward shift of the storm tracks in simulations of 21st century climate, Geophys. Res. Lett., 32, L18701,, 2005. a

Zerefos, C. S., Tourpali, K., Bojkov, B. R., Balis, D. S., and Isaksen, I. S. A.: Solar activity-total column ozone relationships: Observations and model studies with heterogeneous chemistry, J. Geophys. Res., 102, 1561–1569, 1997. a

Zhang, R., Leu, M.-T., and Keyser, L. F.: Heterogeneous Reactions of ClONO2, HCl, and HOCl on Liquid Sulfuric Acid Surfaces, J. Phys. Chem., 98, 13583–13574,, 1994. a, b

Zhang, R., Leu, M.-T., and Keyser, L. F.: Hydrolysis of N2O5 and ClONO2 on the H2SO4/HNO3/H2O ternary solutions under stratospheric conditions, Geophys. Res. Lett., 22, 1493–1496, 1995. a

Publications Copernicus
Short summary
Two developments are made to the United Kingdom Chemistry and Aerosols (UKCA) model to improve simulation of stratospheric ozone. The first is the addition of a solar cycle. The influence on ozone from the solar cycle is found to be 1–2 %, which is consistent with other studies. The second is to the heterogeneous chemistry, the most significant change being the addition of reactions involving bromine species. This was shown to reduce ozone biases relative to observations in most regions.
Two developments are made to the United Kingdom Chemistry and Aerosols (UKCA) model to improve...