Implementation of the Fast-JX Photolysis scheme into the UKCA component of the MetUM chemistry climate model

Introduction Conclusions References


Introduction
In the recent past efficient comprehensive photolysis schemes have been developed that are fast enough to be run within global, three dimensional models, allowing a more interactive treatment of photolysis frequencies in composition modelling (Wild et al., 2000;Tie et al., 2003).These have been implemented into a range of global, three dimensional models (see for instance Liu et al., 2006;Voulgarakis et al., 2009) and have been demonstrated to provide an improved description of atmospheric composition, especially in regions where photolysis and variations in optical depth are important.
The UKCA module is a component of the UK Met Office's Unified Model (Me-tUM).A number of model configurations have been established and tested; one with stratospheric chemistry (Morgenstern et al., 2009), one with tropospheric chemistry (O'Connor et al., 2009;Telford et al., 2010) 1 and one with the GLOMAP-mode aerosol scheme (Mann et al., 2010).There are also whole atmosphere chemistry schemes that have been developed from the tropospheric and stratospheric chemistries (Archibald Figures Back Close Full  et al., 2012;Morgenstern et al., 2012).We describe the addition of the Fast-JX photolysis code into the MetUM framework and evaluate it against a range of observations.Although we make some comparisons in an idealised atmosphere based on the CCMVal model intercomparison (Chipperfield et al., 2010), the main evaluation here is carried out using a "nudged" tropospheric model simulation from 2004 through to 2008.The technique of nudging constrains the model to meteorological re-analyses, allowing the model to reproduce the observed weather (Telford et al., 2008).This allows us to make comparisons with an expanded range of data-sets, including those from campaigns and from satellites.Although most of our validation is performed on the tropospheric model, as this is where the interactive photolysis is more important because of the effects of clouds, we also include a short section on the performance of the stratospheric model.

Model description
The UKCA chemistry module is part of the Met Office's Unified Model (MetUM).We employ a version of the model based on the HadGEM3 configuration (Hewitt et al., 2011), but with the following configuration: a horizontal resolution of 3.75 • × 2.5 • in longitude and latitude.
-60 hybrid height levels in the vertical, from the surface up to a height of 84 km.
A time series of sea surface temperatures and sea ice coverage are prescribed from the HadISST dataset (Rayner et al., 2003).

Nudging
The technique of nudging is used to reproduce the atmospheric conditions over the period studied.The nudging is applied from 3 to 45 km with a relaxation time scale of Figures data (Dee et al., 2011) instead of the ERA-40 re-analyses.These are the latest reanalysis products from the ECMWF with the following changes from the ERA-40 reanalysis: higher horizontal resolution, 4D-var data assimilation, new humidity analyses, and improved model physics, and a different period of coverage (1979( -2012( compared to 1957( -2002 in ERA-40) in ERA-40).

Tropospheric chemistry
The tropospheric version of the model employs a medium sized chemistry scheme that simulates the O x , HO x and NO x chemical cycles and the oxidation of CO, ethane, propane, and isoprene (Telford et al., 2010).The Mainz isoprene mechanism (P öschl et al., 2000) is used to parameterise isoprene oxidation.In total the model has 168 chemical reactions and 56 chemical tracers.We also add the reaction between HO 2 and NO using the rates and yields of Butkovskaya et al. (2007) and hydrolysis of N 2 O 5 (Morgenstern et al., 2009).Concentrations of ozone and NO y are overwritten above 30 hPa.The upper boundary condition for NO y are taken from the Cambridge 2-D model (Law and Pyle, 1993a,b).The upper boundary conditions for O 3 are taken from the Rosenlof climatology (Dall'Amico et al., 2010).Methane is fixed to 1.76 ppmv throughout the atmosphere.Dry deposition and wet deposition are parameterised using the approaches of Giannakopoulos et al. (1999).
Seven chemical species (nitrogen oxide (NO), carbon monoxide (CO), formaldehyde (HCHO), ethane, propane, acetaldehyde and acetone ((CH 3 ) 2 CO)) are emitted in the manner of Zeng and Pyle (2003).Isoprene (C 5 H 8 ) emissions are included in a similar manner, though with a diurnal cycle as in Young (2007).The emissions are taken from a sum of anthropogenic and natural emissions, with anthropogenic and biomass burning emissions taken from Lamarque et al. (2010) using data for the year 2000.We also lump the emissions of ethene and ethyne with ethane, and propene with propane.
Aircraft emissions of NO are taken from Eyers et al. (2004).We also include biogenic emissions of CO, Me of lightning emissions distributed according to the parameterisation of Price and Rind (1994).Finally we include 5.6 Tg yr −1 of NO x from natural soil emissions distributed according to the empirical model of Yienger and Levy II (1995).

Stratospheric chemistry
The stratospheric version of the model has a comprehensive stratospheric chemistry, including chlorine and bromine chemistry, heterogeneous processes on polar stratospheric clouds (PSCs), and liquid sulphate aerosols as well as a simplified tropospheric chemistry (Morgenstern et al., 2009).More detailed descriptions of the stratospheric version of the chemistry can be found in Morgenstern et al. (2009).
There are a few differences between the stratospheric model used in this study and that in Morgenstern et al. (2009).The base version of the climate model has been updated to HadGEM3 from HadGEM1a.In addition N 2 O is now transported completely separately from the other odd nitrogen species and the reaction N 2 O + O( 1 D) has been updated with rates from a more recent JPL assessment (Sander et al., 2006).

Climatological photolysis schemes
In earlier versions of the model, photolysis frequencies were determined from frequencies calculated off-line under average conditions, tabulated and interpolated for use on-line.There were two schemes, one for use in the troposphere, and one for use in the stratosphere.The tropospheric configuration of the chemistry model exclusively used the tropospheric photolysis scheme.The stratospheric configuration of the chemistry model used the tropospheric photolysis scheme below 300 hPa, the stratospheric photolysis scheme above 200 hPa with a linear transition between the two schemes from 200 hPa and 300 hPa.The tropospheric photolysis scheme employed the two stream method of Hough (1988) with clouds included as a zonally averaged climatological distribution (Law et al., 1998).The photolysis frequencies depended only on latitude, local solar time, month of the year and pressure.The stratospheric scheme was based on the look up table approach of Lary and Pyle (1991), with some updated cross section measurements (Morgenstern et al., 2009).The frequencies in the stratospheric scheme responded to changes in the overhead ozone column.
Apart from the stratospheric ozone, these schemes did not respond to changes in optical depth, relying instead on average conditions for particular latitudes and altitudes.This precluded us from looking at a complete set of feedbacks caused by variations in the optical depth (from e.g.clouds, aerosols) onto atmospheric composition, limiting our ability to simulate the atmosphere.The addition of a more complete interactive photolysis scheme allows us to overcome this hurdle.

Fast-J interactive photolysis schemes
Recently photolysis schemes suitable for interactive use, including the Fast-J scheme, have been developed that are fast enough to be incorporated into global models.The original Fast-J scheme was developed for tropospheric photochemistry (Wild et al., 2000).It arranged light from wavelengths between 289 to 850 nm into seven discrete bins.A further development, Fast-J2, extended the scheme into the stratosphere (Bian and Prather, 2002) by adding 11 wavelength bins from 177-290 nm.To speed up the code the rayleigh scattering for the additional wavelengths was treated as pseudo absorption.In the absence of this scattering Fast-J2 fails to reproduce stratospheric photolysis frequencies at low solar zenith angles or at dusk, yielding problems in the high latitude winter stratosphere (Morgenstern et al., 2009) and has difficulties in describing the effects of high aerosol loadings.Fast-JX combines these two schemes, providing the full scattering calculation for all 18 wavelength bins (Neu et al., 2007).It also has several other improvements most notably a more efficient way of introducing extra Figures

Back Close
Full

Technical implementation
We take the Fast-JX code from the Fast-JX website2 , this will now be referred to as Prather et al. (2010).Given the optical depth from absorbing and scattering species this code determines how many photons of each wavelength are absorbed and scattered as light passes through the atmosphere.This is done using the "plane parallel assumption" which assumes that, for a grid box, the horizontal properties are constant and radiative transfer properties only depend on the vertical co-ordinate.Radiative fluxes between horizontally adjacent grid-boxes is also neglected.The calculation of the radiative properties thus divides into a series of columns.The path of the radiation is traced through this column, being scattered or absorbed according to the contents of the grid-box.The amount of scattering and absorption depends on the optical depth of the different scatterers in the grid-box.
The photolysis frequencies ("j " rates) for each reaction are determined from the wavelength bin resolved flux in each grid-box and the cross section of each species in each wavelength bin.These cross sections are evaluated from experimental measurements of the frequencies as described in Sect.3.2.Because Fast-JX still does not provide photolysis frequencies for wavelengths below 177 nm, which are important for some reactions in the upper stratosphere and mesosphere we evaluate the cross sections for these wavelengths using the original climatological scheme Lary and Pyle (1991)

Calculation of optical depth
The optical depth (τ) is determined as the sum of the optical depth of ice water clouds (τ w ), liquid water clouds (τ i ), aerosols, ozone and oxygen and rayleigh scattering.
The optical depth of liquid water clouds is calculated using the parameterisation of Slingo (1989), Here LWP is the liquid water path of the cloud, a i and b i are parameters and r eff is the effective radius of the water droplets.The parameters, which are taken to be −8.9 m 2 kg −1 and 1.67 × 10 3 m kg −1 , are updated from Edwards and Slingo (1996).The effective radius is assumed to be 6 µm over land and 12 µm over water.
The optical depth of ice water clouds is calculated using the non spherical parameterisation of Edwards et al. (2007), Here IWP is the ice water path of the cloud, c i , d i and e i are parameters updated to account for the different spectral bands and d eff is the effective diameter of the ice particles.The values of c i , d i and e i are −2.189m 2 kg −1 , 3.311 × 10 −3 m kg −1 and 3.611 × 10 −12 kg −1 , respectively.The effective diameter is taken to be 100 µm.Neu et al. (2007) noted the importance of accounting for the effects of overlapping cloud layers, describing several means of doing it.We account for this effect using the simple approach of Briegleb (1992) where the optical depth is modified by the factor f 3/2 , where f is the fraction of each grid box area covered by clouds.This has been demonstrated by Feng et al. (2004) to produce a computationally cheap representation of the effects of overlapping cloud layers being used extensively in global models (Liu et al., 2006;Voulgarakis et al., 2009) The only aerosol accounted for at present is sulphate (Bellouin et al., 2007).The effect of hygroscopic growth on particle size and optical depth is parameterised using the approach of Fitzgerald (1975) as used by Bellouin et al. (2007).
All optical depths are determined at 600 nm and then scaled for the different wavelength bins.For the purposes of determining the addition of extra levels in the presence of high optical depths, only the cloud optical depth is used.

Calculation of cross sections
The cross sections are based on v6.4 of the file available from the Fast-JX website (Prather et al., 2010).This data set is a combination from the JPL14 assessment (Sander et al., 2003) with the addition of IUPAC data for NO 2 and VOC photolysis frequencies and of other sundry additions such as the acetone photolysis of Blitz et al. (2004) and the infrared photolysis of HO 2 NO 2 of Jim énez et al. (2005).We have updated these frequencies to include data from JPL15 assessment (Sander et al., 2006), with the most notable changes including updates for N 2 O 5 and OClO photolysis.We summarise all reactions described, noting whether they are utilised by the stratospheric or tropospheric chemistry, in Table A1.

Methodology
To test the performance of the model we carry out a series of tests.The first of these is to investigate how the photolysis code performs under idealised conditions.To do this we employ the scenarios devised for the PHOTOCOMP photolysis study as part of the CCMVal chemistry climate model intercomparison (Chipperfield et al., 2010).
We then proceed to investigate how the MetUM model with the Fast-JX photolysis included performs.direct comparisons we run the model in "nudged mode" (Telford et al., 2008).This constrains the model to meteorological observations making comparisons with data over our chosen years more meaningful.
The first brief check we make is to ensure that nudging with the ERA-Interim dataset is as successful as nudging with the ERA-40 dataset.We do this by comparing biases, correlations and root mean squared differences between the nudged model and analyses as used by Telford et al. (2008).We then proceed to compare photolysis frequencies with observations from the INTEX-NA campaign, sampling the model at the time and location of the measurements.We make a quantitative assessment of the climatological and interactive photolysis schemes by comparing biases, correlations and root mean squared differences (RMSD) between the model and the observations.
We then study the effects these changes in photolysis frequencies have on atmospheric composition, focusing on ozone and carbon monoxide.Our principal comparisons are with observations from the TES satellite, comparing both the distributions of CO and O 3 and their correlation.We sample the model in a similar manner to Rodgers and Connor (2003), taking measurements at the appropriate time and location as the satellite and applying the instrument's averaging kernels to make comparisons more meaningful.To understand the changes in the chemical species distributions we also study the ozone budgets and changes in OH.

Idealised comparison
We first check the photolysis frequencies of NO, O the photolysis in different models.The three scenarios used clear sky conditions with a solar zenith angle of 15 • , 84 • and averaged over a day between 84 • and 96 • .
First we compare photolysis frequencies as a function of altitude in a series of idealised atmospheric profiles constructed for a multi model comparison study.The photolysis frequencies for NO, O 2 and O 3 and Cl 2 O  Given the general successful performance under idealised conditions we move on to studying the effects of the photolysis code within the wider chemistry climate model.

Validation of nudged tropospheric MetUM run
After demonstrating that the photolysis schemes works well under idealised conditions we test the performance of the scheme within a MetUM model run.Figures

Back Close
Full

Validation of nudging with ERA-Interim analyses
As this is the first time we have used ERA-Interim re-analysis data in the nudging we briefly re-evaluate the nudging performance.To do this we use the simple statistical tests of Telford et al. (2008), calculating the bias, correlation and root mean squared difference between the nudged model and the ERA-Interim data.These results are shown in Table 1.The performance of the nudging is broadly similar to that when using ERA-40 (compare to Table 1 in Telford et al., 2008), with the most notable difference being the reduction in the RMSE of potential temperature (θ), which probably reflects that improvements in the ERA-Interim analyses.From these results we conclude that the model nudged towards the ERA-Interim analyses can be used as a hindcast of the period 2004 to 2008.This is consistent with the results of Kipling et al. (2012) who show that nudging towards ERA-Interim analyses was able to improve modeled aerosol properties.

Comparison of cloud optical depth
The greatest values and variability in the tropospheric optical depth arise from clouds.
To demonstrate the model is able to provide a reasonable description of this cloud optical depth we compare our modelled distribution with satellite data.Because of its high temporal resolution, we chose to compare to the ISCCP (International Satellite Cloud Climatology Project) D2 data-sets (Rossow et al., 1996;Rossow and Schiffer, 1999), which are constructed from the output of several NOAA (the American National Oceanic Atmospheric Administration) satellites.The comparisons for December, January and February 2005 and June to August 2005 are shown in Fig. 2.
The model is able to capture many of the features of the observations.For instance the storm tracks in the Northern Atlantic and Pacific oceans are well reproduced as are the decks of stratus clouds off the coast of Peru.There are some discrepancies between the model and data, for instance the model underestimates the optical depth in the inter tropical convergence zone (ITCZ).Some of the other differences can be Introduction

Conclusions References
Tables Figures

Back Close
Full attributed to issues with the data, with errors introduced at high solar zenith angles, contamination from dust and an underestimate of the amount of optically thick cloud (Marchand et al., 2010).For instance the discrepancy between the model and observations over the Sahara is a result of factors such as dust in the measurements rather than a deficiency in clouds in the model.However some of the discrepancies can be attributed to the simplified calculation of cloud optical depth and the limited resolution of the model.From the general good agreement we conclude that the model is able to provide a reasonable description of the cloud optical depth, with the agreement being at least as good as that in other models (see e.g.Voulgarakis et al., 2009).

Comparison of photolysis frequencies with INTEX campaign data
Although the modelled photolysis performs well under an idealised scenario, to more fully evaluate its performance we compare the modelled photolysis frequencies, with the interactive and climatological photolysis schemes, to observational data.The observations we choose to use are from the INTEX-A flight campaign, conducted over North America and the Atlantic in the summer of 2004 (Singh et al., 2006).On these flights the actinic flux was measured in different wavelength regions using spectral radiometers and photolysis frequencies are obtained by multiplying by cross sections from Sander et al. (2003).This procedure to measure photolysis frequencies has been reported to be accurate to around 15-20 % (Shetter and M üller, 1999).The O 3 photolysis rate is shown to be significantly improved by the use of the interactive photolysis frequencies, with the low bias removed.The capture of the variability of the photolysis frequencies is also improved, reflecting the improved ability of the interactive photolysis frequencies to model the effects of variations in optical depth.The improvement in the NO 2 photolysis rate is not as dramatic, though it can be seen that the interactive photolysis reduces a smaller bias and captures the variability better.
Figure 3 only shows one flight from the entire campaign.Although we chose this at random, selecting the first flight in the data, it is possible that the conditions might favour one photolysis scheme in particular.Therefore we calculate the bias, correlation and, root mean squared difference (RMSD) between the modelled photolysis frequencies and the observed frequencies for all flights.The biases provide information on average frequencies.The correlations indicate the ability of the model to capture the variability in the photolysis.The RMSD combines information about average agreement and variability.The biases, correlations and RMSD are given in Table 2 for the two frequencies shown in Fig. 3, along with selected other frequencies.
The tabulated results are in accord with Fig. 3, with the interactive photolysis having lower biases and higher correlations with the data than the climatological photolysis.For several reactions, including O 3 photolysis, the discrepancy between the climatological photolysis and the data is considerably larger than the reported experimental uncertainty.For the interactive photolysis there are much smaller differences between the average model and observed values.Some of the reactions where there are differences between the model and the interactive photolysis and data, most notably HONO, are those where there have been updates to the reference cross sections.In the case of HONO the values (Sander et al., 2006) used in the interactive photolysis now incorporate additional measurements (Kenner and Stuhl, 1986;Stutz et al., 2000).This is relevant as the observed photolysis frequencies are not measured directly, but are a product of the measured radiation multiplied by cross sections from Sander et al.Introduction

Conclusions References
Tables Figures

Back Close
Full  2003).So some of the discrepancy can be attributed to an improved understanding of the experimental cross sections.Whilst some of the bias with the climatological scheme can be attributed to the use of old photolysis frequency measurements and a simple treatment of cloud optical depth, much can be attributed to an overestimation of the stratospheric ozone column.
By inspecting Fig. 3 we also see, as expected, that the global model, with either photolysis scheme, is not designed to describe small scale variations.However the use of the interactive photolysis, rather than the climatological scheme does produce notable improvements.The correlation between the interactive photolysis and data is always greater than 0.8 and only greater than 0.8 for the photolysis frequencies of one species, HCHO, with the climatological scheme.This is as expected as the photolysis frequencies are now sensitive to changes in cloud optical depth.
Although the Fast-JX code is designed to run using 18 wavelength bins it can be run using a subset of 8 or 12 of these, where the lower wavelength bins are excluded.We repeated the comparison between the model with the default Fast-JX model and with 8 and 12 bins.The largest differences are found in the ozone photolysis frequencies.The photolysis frequencies with 8 and 12 bins only differ slightly from each other, but at high altitudes are appreciably lower than the default scheme.The difference can be as large as 7 % with an average bias of 1 %.The correlation and RMSE are unchanged to two significant figures.However for studies focussed on the lower troposphere the reduced wavelength bin schemes are acceptable, though for studies that examine the upper troposphere and stratosphere the default setting, that of using all 18 wavelength bins, would be required.

Comparison of ozone and carbon monoxide to TES measurements
After demonstrating that the photolysis frequencies look reasonable we investigate their effects on the tracer distributions, comparing model simulations, with the climatological and interactive photolysis schemes, and measurements.The large spatial and Figures

Back Close
Full temporal scales of satellite observations make them ideal to be employed in this capacity.We compare to observations of ozone and carbon monoxide from the TES satellite.The Tropospheric Emission Spectrometer (TES) is an infrared Fourier transform spectrometer (Beer et al., 2001) onboard NASA's Aura satellite which was launched in 2004.In this study we utilise CO and ozone profile observations from the TES Global Survey mode.In this mode profiles have a nadir footprint of 5.3×8.3km and are spaced approximately 180 km apart along the orbital track with a 16 day repeat cycle.
In order to be able to make systematic comparisons between the model and the TES observations we have developed code that samples the model online at the times and locations of the satellite measurements.This differs slightly from the technique used in previous studies, such as those by Voulgarakis et al. (2011), who took relatively high frequency global (six hourly) averages and sampled these as close to the time and location of the TES measurements as possible.After interpolating the modelled values of O 3 and CO onto pressure levels used by the TES retrievals, we here we apply averaging kernels following the method of Rodgers and Connor (2003), a procedure which is shown to be vital by the study of Aghedo et al. (2011).We bin the O 3 and CO data from the model and the data into 4 • × 5 • bins and average from 800-400 hPa.
First we compare the ozone distributions between the model, with interactive and climatological photolysis schemes, and the TES observations (Nassar et al., 2008) (Fig. 4).For consistency with Voulgarakis et al. (2011) we compare for the periods July to August and January to February.With either photolysis schemes the model captures the main features such as the low ozone values over the Western tropical Pacific and the products of biomass burning from Southern Africa.The modelled ozone, with either photolysis schemes, is slightly higher (9 % averaged over both periods) than the observations.Validation of TES ozone profiles against ozonesondes and LIDAR observations have shown that TES itself has a mean high bias in this altitude region of 5-10 % (Nassar et al., 2008;Richards et al., 2008).This contrasts to the previous version of the model which had a slight low bias (see Fig. 2 in Voulgarakis et al., 2011) .This increase is believed to mainly arise from changes in the underlying climate model.To probe why the changes in photolysis have such a slight impact on the ozone we tabulate the budget of ozone in Table 3, where we compare the chemical production and loss of ozone, its source from the stratosphere (STE) and its sink to the surface by deposition.We also include the total ozone burden and methane lifetime.As well as the budget from the model with climatological and interactive photolysis we include a multi model average from Stevenson et al. (2006) and results from another version of the MetUM model that employs the interactive photolysis, albeit with simpler tropospheric chemistry (Morgenstern et al., 2012).
We can see that the change in photolysis scheme has a large effect on the chemical production and loss terms, but that the overall, global, effect is small as the changes balance each other.The interactive photolysis has production and loss terms more similar to the multi model average of Stevenson et al. (2006).The model of Morgenstern et al. (2012) has smaller chemical fluxes which can be accounted for by its simpler tropospheric chemistry.However, whilst using the same implementation of the interactive photolysis scheme, it produces a relatively low ozone burden and high methane lifetime, indicating that the low methane lifetime seen in the tropospheric chemistry scheme in Table 3 is not solely a result of the interactive photolysis.
Next we compare the carbon monoxide distributions between the model, with interactive and climatological photolysis schemes, and the TES observations (Luo et al., 2007;Lopez et al., 2008) (Fig. 5).Again we compare for the periods July to August and January to February.
The differences between the two versions of the model are now more dramatic.Again both versions capture the broad features with high values of carbon monoxide produced by biomass burning, in particular from Africa, and anthropogenic emissions from Asia.However whilst the model with climatological photolysis has only a small bias (2 %), when the interactive photolysis is employed there is a large negative bias (−24 %).Although there are small positive biases (less than 10 %) in the tropics and small negative biases in midlatitudes (less than 10 %) in the observations (Luo et al., 2007;Lopez et al., 2008) are unchanged between the two models then we ascribe these changes to the main sink of CO, its reaction with the hydroxyl radical, OH.Indeed if we look at the methane lifetime (Table 3), which is dominated by its reaction with OH, we see that this is greatly reduced by switching from the climatological to interactive photolysis scheme.We can look directly at the changes in the OH by plotting the OH field, weighted as in Lawrence et al. (2001), with the climatological and interactive photolysis schemes (Fig. 6).
We first note that the use of the interactive photolysis scheme increases the total OH burden, which explains the decreased CO concentrations and CH 4 lifetimes.We also note that the amount of OH is now significantly higher than that reported by Spivakovsky et al. (2000) leading us to believe that the amount of OH is excessive.To understand the increase we examined the production and loss terms of HO x (≡ OH+HO 2 ) and found that the increases in HO x are dominated by increases in its production via O( 1 D)+H 2 O.The water vapour is constant between the two simulations and in a recent study found to compare reasonably well with other models and observations (Russo et al., 2011).Therefore the increase must arise from increased production of O( 1 D), which is produced from the photolysis of ozone.This is greatly increased by using the interactive scheme as can be seen from Fig. 3. However we believe that the increased j (O( 1 D)) is realistic.
The apparently superior performance of the climatological photolysis scheme is a result of two conflicting factors, the high bias in ozone and the low bias in the ozone photolysis frequencies.Removing one of these factors, although an improvement in itself, thus worsens the performance of some aspects of the model.Indeed Morgenstern et al. ( 2012) see no such problems with high biases in OH despite using this implementation of the interactive photolysis scheme, albeit with a considerably different chemistry scheme.Results from a more complete whole atmospheric chemistry (Archibald et al., 2012) also indicate lower ozone concentrations produce lower OH concentrations, higher CO concentrations and a longer methane lifetime.
Finally, adopting the approach of Voulgarakis et al. (2011), we compare the O 3 -CO correlations using them to understand the models performance (Fig. 7).Like previous Introduction

Conclusions References
Tables Figures

Back Close
Full versions of the model (Voulgarakis et al., 2011) the O 3 -CO correlation is lower in the MetUM model than in the observations.This may reflect the use of climatological biomass burning emissions, which would tend to reduce correlations.On a global scale the two versions of the model perform similarly with "scores" of 0.3 for both photolysis schemes for both periods.The "scores" are obtained by correlating the tracer-tracer correlations in the model and observations.However there are different regions where the two photolysis schemes perform better, with the interactive photolysis producing better agreement near sources of biomass burning (S.America; Central Africa) and the climatological scheme performing more strongly in the extra-tropics.This worsening of performance of the interactive photolysis scheme can be linked to the high OH concentrations attenuating CO concentrations too rapidly with distance from sources.

Conclusions with regards to tropospheric run
We have evaluated the performance of the Fast-JX photolysis scheme in the tropospheric version of the MetUM model using a nudged run between 2004 and 2008.
After briefly checking that nudging to the ERA-Interim analyses is successful, and that the modelled cloud optical depths were realistic, we proceeded to demonstrate improvements to the modelled photolysis frequencies using comparisons to data from the INTEX-NA campaign.These results showed that the online photolysis model provides a significantly better description of the observed frequencies.We then investigated the effects that these improved photolysis frequencies have on the modelled chemical fields, comparing the changes to observations from the TES satellite.The improved photolysis frequencies, whilst producing some improvements to the chemical fields, produced a high bias in the global burden of OH, which itself affected CO concentrations and the CH 4 lifetime.This was as a result of the more realistic photolysis frequencies increasing the sensitivity of the model to high biases in ozone, the effects of which in the climatological photolysis scheme were masked by too small ozone photolysis frequencies.Introduction

Conclusions References
Tables Figures

Back Close
Full

Stratospheric chemistry
In addition to testing the model in the troposphere a further simple validation was performed in the stratosphere.In order to test the feedback from chemical changes onto the dynamics we ran the model unconstrained by nudging.The run we perform is similar to that used in Morgenstern et al. (2009), which itself follows the definition of the REF-B0 experiment of the CCMVal model intercomparison (Eyring et al., 2010).This uses forcings that represent perpetual 2000 conditions, although, as we are able to initialise the chemical fields from those in Morgenstern et al. (2009) we only perform a 10 yr run, discarding the first five years as "spin-up".
In Fig. 8 we examine the average annual cycle of the total ozone column.The values from the model are compared to a combination of data from a variety of sources including TOMS, SBUV, GOME and OMI, for the year 2000 (Bodeker et al., 2005).The model captures most of the features including the latitudinal gradient and timing and magnitude of the Northern Hemisphere spring maximum.
The previous version of the chemistry scheme had some discrepancies with observations around the tropopause, including high biases in temperature and ozone (Morgenstern et al., 2009).To assess how the updated model performs we compare ozone profiles with those from the SHADOZ network (Thompson et al., 2003a,b).We show results from one representative site, Fiji, in Fig. 9 displaying monthly mean ozone profiles as a function of pressure.The model is able to capture the seasonal cycle well, though it does still slightly overestimate UTLS ozone.Tropospheric ozone is slightly low, as should be expected, as this version of the model does not attempt to fully simulate tropospheric emissions and chemistry.

Discussion
Whilst we have demonstrated that the photolysis frequencies are improved by the use of the interactive photolysis scheme and that the stratospheric version of the model Introduction

Conclusions References
Tables Figures

Back Close
Full performs well with it there are still obvious issues with the tropospheric model.We have demonstrated that these arise from high biases in the ozone concentrations which, without the large low bias in photolysis frequencies seen in the climatological scheme, now produce a high bias in OH.
We have reduced the O 3 high bias by various measures, including the replacement of the original upper boundary condition for NO y as used by O'Connor et al. (2009) with one from the Cambridge-2D model, adding of the reaction between HO 2 and NO using the rates and yields of Butkovskaya et al. (2007) and adding the hydrolysis of N 2 O 5 (Morgenstern et al., 2009).Whilst these reduced the biases in ozone they did not eliminate them.Further work is being carried out at present to try and understand the sources of these biases, but they cannot obscure the improved description of observed photolysis rates by the interactive photolysis scheme as opposed to the climatological scheme.
Even with the resolution of the issues with biases in ozone there are still a potential ways of expanding the capability of the interactive photolysis scheme.One obvious omissions regards the limited use of aerosols in the model.This will change when the GLOMAP-mode scheme (Mann et al., 2010) is coupled to the photolysis.Although on a global scale aerosols are less important than clouds for the optical depth calculation there are specific episodes where this coupling will be useful.A further improvement, of more import for the stratosphere, will be the use of a variable solar constant, a driver of some of the variability in stratospheric ozone.Whilst our simple approach of accounting for overlapping cloud layers seems to work at present, increasing vertical resolution may require more sophisticated techniques to be employed (Neu et al., 2007).The new techniques we describe for validating the model will also be useful in assessing other aspects of the model performance.Introduction

Conclusions References
Tables Figures

Back Close
Full

Conclusions
We describe the implementation of the Fast-JX interactive photolysis scheme in the UK Met Office Unified Model.The interactive photolysis scheme is shown to improve the modelled photolysis frequencies from the previous, climatological scheme in comparison to observations and other models.The effect on atmospheric composition is less clear cut, with the interactive photolysis scheme producing too much OH.However this is believed to be a result of a previous cancellation of errors being removed, and not as a result of any inherent problem with the interactive photolysis scheme.As well as being a more realistic manner of modelling the photolysis the interactive photolysis will permit more detailed studies of the interactions between the atmosphere and its constituents.Figures

1 of
2 CO and C 5 H 8 and emissions of NO from natural soil emissions and lightning.The biogenic emissions are distributed according to Guenther et al. Discussion Paper | Discussion Paper | Discussion Paper | (1995) and total 40 Tg yr −1 of CO from the oceans and 45 Tg yr −1 of CO, 40 Tg yr −Me 2 CO and 570 Tg yr −1 of C 5 H 8 from the land.We include, on average, 5 Tg yr −1 Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | levels for very optically dense clouds.It is the implementation of this code into the MetUM model that we describe.
and add them to the Fast-JX reaction frequencies.For speed this extra calculation is only performed above 20 Pa as below this level the flux of these high energy photons is negligible.Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | To do this we perform a run with the tropospheric chemistry version of the model, from 2004 to 2008 and compare to different sets of observations including data from the INTEX campaign and the TES satellite.To enable us to make more Discussion Paper | Discussion Paper | Discussion Paper | 2 and O 3 3 under idealised conditions using the three assessed scenarios from the PHOTOCOMP assessment of the CCMVal model intercomparison (see Sect. 6.3.1 in Chipperfield et al., 2010).These employed idealised atmospheric conditions to allow systematic comparisons between Discussion Paper | Discussion Paper | Discussion Paper | 2 are shown in Fig.1, which can be compared to Fig.6.1 inChipperfield et al. (2010).The O 3 photolysis frequencies look similar to the PHOTOCOMP average, reducing from around 10 −2 s −1 at 0.1 hPa to between 10 −4 -10 −3 s −1 near the surface, depending on the solar zenith angle.Similarly the O 2 photolysis frequencies look similar to the PHOTOCOMP average, reducing from around 10 200 hPa depending on the solar zenith angle.The NO photolysis frequencies are the exception, being high biased throughout the atmosphere.At the top of the atmosphere the photolysis frequencies are around 10 −5 s −1 in contrast to the PHO-TOCOMP average, which is nearer to 10 −6 s −1 .This is consistent with the models that employ Fast-JX in the PHOTOCOMP study, which see good agreement with the multi model mean photolysis frequencies for O 2 and O 3 , but report a high bias in the case of NO.This discrepancy has been linked to factors such as a neglect of NO absorption in the calculation of the optical depth, leading to an overestimate of the modelled cross section.
Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | The use of the nudged version of the MetUM model allows us to make meaningful comparisons with the campaign data on a measurement by measurement basis.To do this we sample the model at the same time and location as the measurements were made, an approach similar to O'Connor et al. (2005) and employed by Kipling et al. (2012), interpolating the model to the pressure and co-ordinates where the measurements were made.As noted by O'Connor et al. (2005), the relatively low resolution that we use will prevent us from resolving small scale features in the observations.However we expect that we can produce the larger scale features relating to altitude, latitude and time of day.Figure3shows the O 3 and NO 2 photolysis frequencies from the first flight Introduction Discussion Paper | Discussion Paper | Discussion Paper | of the INTEX-NA campaign (1 July 2004) and in the MetUM model, with and without interactive photolysis frequencies.
Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | , this model bias is significant.As the emissions Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Dall'Amico, M., Gray, L., Rosenlof, K., Scaife, A., Shine, K., and Stott, P.: Stratospheric temperature trends: impact of ozone variability and the QBO, Clim.Dynam.
Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper |

Fig. 6 .Fig. 7 .Fig. 8 .Fig. 9 .
Fig. 6.OH distribution, weighted as in (Lawrence et al., 2001), as a function of height and latitude.The black numbers are the model values, the red numbers are from the Spivakovsky et al. (2000) climatology (obtained from Lamarque et al., 2012) and the numbers in brackets the standard deviation of the model.Top plot with climatological photolysis and bottom plot with interactive photolysis.

Table 1 .
Telford et al. (2008)nt of model performance in October 2005 with ERA-Interim nudging using the statistical assessments ofTelford et al. (2008)for potential temperature (θ) and zonal wind (u).This calculates the model mean on four representative model levels and determines the root mean squared error (RMSE) and correlations over time (TC) and space (SC) with respect to the ERA-Interim analyses on these same levels.− 0.13 m s −1 1.17 m s −1 0.97 1.00 35 16.19 − 0.22 m s −1 1.37 m s −1 0.97 1.00 Introduction

Table 2 .
Quantitative comparison between photolysis frequencies from INTEX data and modelled frequencies with interactive and climatological photolysis.