## A new SATIRE-S spectral solar irradiance reconstruction for solar cycles 21-23 and its implications for stratospheric ozone

## William T. Ball, ∗

Physics Department, Blackett Laboratory, Imperial College London, SW7 2AZ, UK

## Natalie A. Krivova

Max-Planck-Institut f¨r Sonnensystemforschung, 37077 G¨ttingen, Germany u o

## Yvonne C. Unruh

Physics Department, Blackett Laboratory, Imperial College London, SW7 2AZ, UK

Joanna D. Haigh

Physics Department, Blackett Laboratory, Imperial College London, SW7 2AZ, UK Grantham Institute, Imperial College London, SW7 2AZ, UK

## Sami K. Solanki

Max-Planck-Institut f¨r Sonnensystemforschung, 37077 G¨ttingen, Germany u o School of Space Research, Kyung Hee University, Yongin, Gyeonggi 446-701, Korea

## ABSTRACT

We present a revised and extended total and spectral solar irradiance (SSI) reconstruction, which includes a wavelength-dependent uncertainty estimate, spanning the last three solar cycles using the SATIRE-S model. The SSI reconstruction covers wavelengths between 115 and 160 000 nm and all dates between August 1974 and October 2009. This represents the first full-wavelength SATIRE-S reconstruction to cover the last three solar cycles without data gaps and with an uncertainty estimate. SATIRE-S is compared with the NRLSSI model and SORCE/SOLSTICE ultraviolet (UV) observations. SATIRE-S displays similar cycle behaviour to NRLSSI for wavelengths below 242 nm and almost twice the variability between 242 and 310 nm. During the decline of last solar cycle, between 2003 and 2008, SSI from SORCE/SOLSTICE version 12 and 10 typically displays more than three times the variability of SATIRE-S between 200 and 300 nm. All three datasets are used to model changes in stratospheric ozone within a 2D atmospheric model for a decline from high solar activity to solar minimum. The different flux changes result in different modelled ozone trends. Using NRLSSI leads to a decline in mesospheric ozone, while SATIRE-S and SORCE/SOLSTICE result in an increase. Recent publications have highlighted increases in mesospheric ozone when considering version 10 SORCE/SOLSTICE irradiances. The recalibrated SORCE/SOLSTICE version 12 irradiances result in a much smaller mesospheric ozone response than when using version 10 and now similar in magnitude to SATIRE-S. This shows that current knowledge of variations in spectral irradiance is not sufficient to warrant robust conclusions concerning the impact of solar variability on the atmosphere and climate.

## 1. Introduction

There is substantial evidence to suggest that changes in the solar irradiance influence variations in the temperature and circulation of the Earth's atmosphere over the 11-year solar cycle. Many of these results are based on correlations with the 10.7 cm solar flux (e.g. Labitzke and van Loon (1995); van Loon and Shea (1999)) or the wavelengthintegrated, or total, solar irradiance (TSI); see Haigh (2003) and references therein. While TSI is a good indicator of the total solar forcing on the climate, it cannot be used to understand the physical interaction between the solar radiation and the atmosphere since spectral solar irradiance (SSI) variability, and the altitude in the atmosphere at which it is absorbed, is highly wavelength-dependent (Meier 1991; Lean et al. 1997; Krivova et al. 2006).

There is a growing body of evidence to suggest that TSI, and as a consequence SSI, may vary on secular timescales exceeding the 11-year solar cycle. Fr¨hlich (2009), o

Lockwood et al. (2010) and Ball et al. (2012) provide some evidence that TSI may have been slightly lower in the re- cent minimum compared to the two prior to that, though the PhysikalischMeteorologisches Observatorium Davos (PMOD) composite of TSI observations (Fr¨hlich o 2006) and the modelled TSI by Ball et al. (2012) are consistent, within the error bars, with no change between the last three minima. Estimates of the increase in TSI since the 17th century vary widely (see Schmidt et al. (2012) and Solanki and Unruh (2013) and references therein), though most recent estimates lie in the range ∼ 1-1.5 Wm -2

(Wang et al. 2005, Krivova et al. 2007, Steinhilber et al. 2009, Krivova et al. 2010), a change similar to solar cycle variability.

A large proportion of the variability in TSI is due to very much larger relative variations at UV wavelengths compared to the longer visible and infra-red (IR) wavelengths. Wavelengths shorter than 400 nm account for less than 10% of the absolute value of TSI, but contribute 3060% to TSI variability, according to models (Lean et al. 1997; Krivova et al. 2006) and measurements by the SOLar Stellar Irradiance Comparison Experiment (SOLSTICE) (Rottman et al. 2001) and the Solar Ultraviolet Spectral Irradiance Monitor (SUSIM) (Floyd et al. 2003; Morrill et al. 2011) on the Upper Atmosphere Research Satellite (UARS), made prior to 2006, and the SCanning Imaging Absorption SpectroMeter for Atmospheric CHartographY

(SCHIAMACHY) on the ENVIronmental SATellite (ENVISAT) (Pagaran et al. 2009). The data from the Spectral Irradiance Monitor (SIM) (Harder et al. 2005) and SOLSTICE (Snow et al. 2005) instruments on board the SOlar Radiation and Climate Experiment (SORCE) indicate that this contribution might be as high as 180% (Harder et al. 2009), though it should be noted that SIM is currently undergoing a reanalysis. A value larger than 100% is possible because the SSI in the visible measured by SIM varies in antiphase to that in the UV.

The UV radiation influences many processes in the atmosphere. Of particular interest is the interaction between solar UV radiation and ozone, which is the largest contributor to heating in the stratosphere. Variation of solar UV radiation over secular timescales may have an effect on global temperature trends and the impact is important to quantify.

Haigh et al. (2010) and Merkel et al. (2011) both investigated the potential impact that the SSI changes observed by SORCE (Rottman 2005; Harder et al. 2009) could have on stratospheric ozone concentrations, respectively using a coupled chemistry climate 2D atmospheric model and the fully 3D general circulation Whole Atmosphere Community Climate Model (WACCM). Both studies obtained qualitatively similar results when using hybrid SORCE data from SOLSTICE and SIM, though the two studies adopted different wavelengths to transition between SOLSTICE and SIM. While the magnitude and the exact heights varied, both studies found that between 2004 and 2007, when solar UV output was declining, ozone concentrations increased above ∼ 45 km while they decreased below ∼ 40 km. It is interesting to note that trends in O 3 from the Microwave Limb Sounder (MLS) on the Aura satellite and the Sounding of the Atmosphere using Broadband Emission Radiometry (SABER) on Thermosphere Ionosphere Mesosphere Energetics Dynamics (TIMED) observations presented by Haigh et al. (2010) and Merkel et al. (2011), respectively, suggest that SORCE SSI may better capture solar variability than models due to the negative response in mesospheric ozone that is out-of-phase with solar irradiance changes. On the other hand, Austin et al. (2008) did not find a negative ozone response to cycle changes of the Sun using combined data from several satellites prior to 2004. Dhomse et al. (2013) suggest that the negative response in the lower mesosphere cannot be used to distinguish between SSI datasets due to the large uncertainties in the ozone observations.

The significant differences in SSI variability between SORCE data and different models, the latter partly relying on earlier observations (see Ermolli et al., 2013), indicate that there is still much uncertainty in our knowledge of how the Sun's irradiance varies spectrally. The larger UV irradiance variability, and an inverse solar-cycle trend in the visible measured by the SIM instrument on the SORCE satellite (Harder et al. 2005, 2009), may indicate that the solar cycle variability observed by previous missions, and the models that reproduce similar behaviour, may be incorrect. It may also indicate a change in the Sun during the recent cycle. However, recent studies suggest that incomplete accounting for instrument degradation may contribute to the SSI trends suggested by SIM data (Ball et al. 2011, Deland and Cebula 2012, Lean and DeLand 2012, Ermolli et al. 2013).

This paper presents an extended and recalibrated data set using the Spectral And Total Irradiance REconstruction (SATIRE-S) model (Fligge et al. 2000; Krivova et al. 2003; Wenzler et al. 2004; Ball et al. 2012), for wavelengths between 115 and 160 000 nm and on all days between August 1974 and October 2009, for use by the climate and atmospheric research communities.

SATIRE-S is the most detailed of the SATIRE family of models, where 'S' refers to the model designed for the satellite era (Krivova et al. 2011), and provides the most reliable reconstruction of TSI and SSI. It does, however, rely on the availability of magnetograms and continuum intensity images, which restricts its applicability to a comparatively short period of time. In the past, the TSI and SSI reconstructed with SATIRE-S have been further limited by the fact that magnetograms obtained from different instruments do not have the same spatial resolution, noise level or magnetic field calibration. It requires careful intercalibrations between various magnetographs (and imagers) to allow a homogeneous reconstruction of TSI. This

has been successfully done by Wenzler et al. (2004, 2006) for Kitt Peak Solar Observatory (KP) 512-channel magnetograph (512) (Livingston et al. 1976a,b) and spectromagnetograph (SPM) (Jones et al. 1992) instruments and by Ball et al. (2012) for 512, SPM and MDI instruments (Scherrer et al. 1995), the latter of which is onboard the Solar and Heliospheric Observatory (SoHO) spacecraft.

Here we compute the SSI using magnetograms from all three instruments, thus extending the SSI reconstructed by SATIRE-S to fully cover the last three solar cycles including, for the first time, the extended solar minimum in 2008. The reconstruction is compared with the NRLSSI model (Lean 2000; Lean et al. 2005) and data from the SOLSTICE instrument (McClintock et al. 2005), onboard the SORCE satellite. We then show how the different spectral irradiances of these datasets affect changes in stratospheric O using the atmospheric model based on Harwood and Pyle 3 (1975).

## 2. Modelling solar irradiance with SATIRE

The SATIRE-S model (Fligge et al. 2000; Krivova et al. 2003; Wenzler et al. 2006; Krivova et al. 2011) assumes that all irradiance variations are the result of changes in the surface photospheric magnetic flux. SATIRE-S identifies four solar surface components in magnetograms and continuum intensity images: the background quiet sun; the dark penumbral and umbral components of sunspots; and small-scale magnetic features, that appear predominantly bright, called faculae.

Daily irradiance spectra are produced by summing the intensities of the four components weighted according to their surface distribution. The component intensities (as functions of wavelength and limb angle) are calculated with the spectral synthesis program ATLAS9 (Kurucz 1993) assuming local thermodynamic equilibrium (LTE) conditions. We use time-independent model atmospheres

(Fligge et al. 2000; Krivova et al. 2003; Solanki and Unruh 2013) with effective temperatures of 5777 K, 5450 K and 4500 K for quiet Sun, penumbral and umbral intensities, respectively. For faculae, we use the FAL-P model atmosphere (Fontenla et al. 1993), as modified by Unruh et al. (1999). The wavelength grid of the daily spectra from SATIRE-S has a resolution of 1 nm below 290 nm, 2 nm from 290 to 1000 nm, 5 nm from 1000 to 1600 nm, 10 nm from 1600 to 3200 nm, 20 nm from 3200 to 6400 nm, 40 nm from 6400 to 10 000 nm and 20 000 nm for the remainder of the spectrum up to 160 µ m.

The model has one free parameter; this relates the magnetic flux registered in a magnetogram pixel to the fraction of the pixel filled by faculae. The free parameter is set to a fixed value for each observatory (i.e., for KP and SoHO) as outlined in the next section.

## a. Method to combine reconstructions

In order to maximise the length of the SSI timeseries, magnetograms and continuum intensity images are taken from three instruments: two at KP that are based on spectropolarimetry of the Fe i 868.8 nm line (Livingston et al. 1976a), the KP/512 (Livingston et al. 1976b) and KP/SPM instrument (Jones et al. 1992), and the SoHO/MDI instrument that uses the Ni i 676.8 nm line (Scherrer et al. 1995). The free parameter for each instrument is fixed by comparing the reconstructed TSI to either TSI observations or to a TSI reconstruction made using images from a different instrument. Broadly, three steps are involved in the intercalibration, which are outlined below (see also Ball et al., 2012). The uncertainties arising in this process are outlined in section 2c.

In step (i), we fix the free parameter for the MDI reconstruction by requiring a regression slope of unity between the reconstructed TSI and the SORCE Total Irradiance Monitor (TIM) TSI observations (Kopp and Lawrence 2005). In step (ii), we combine the KP and MDI magnetogram and continuum images using the KP/SPM and SoHO/MDI overlap period of 895 days between 1999 and 2003. This requires fixing the free parameter for SPM, so that the reconstructed TSI during the overlap period agrees with the TSI derived from the MDI images. However, while the resulting spectral irradiances are very well correlated for the overlap period ( r c &gt; 0 91 at all wavelengths), we see slightly dif-. ferent variability amplitudes in the two reconstructions at some wavelengths. The different instrument design, the use of different spectral lines with different magnetic sensitivities, as well as different telescope optics and detectors mean that there are non-linear, position-dependent differences in the KP and MDI instruments in response to magnetic flux that leads to the differing variability amplitude. The differences in amplitude of variability in the overlap period are typically 2% in the visible and near-IR and remain below 8% at all wavelengths. To avoid discontinuites in the SSI trends when changing between reconstructions based on MDI and KP magnetograms, we adjust the variability amplitudes of the KP reconstructions by rescaling them to those of the MDI reconstructions.

Step (iii) involves correcting for the change between the 512 and SPM instruments on KP. While the imaging quality for the KP/512 data is poorer, the two KP polarimeters show very similar flux registration so that the correction can be used to convert the KP/512 magnetogram signal to the KP/SPM level (see Wenzler et al., 2006, and Ball et al., 2012). Thus, the same filling factor can be used for both KP data sets. While the scaling factor introduces uncertainties regarding the long-term TSI behaviour, it does not affect its spectral distribution.

The ATLAS9 model intensities assume LTE conditions in the solar atmosphere; this can result in large errors

in the modelled irradiance variability in some wavelength regions, mainly below 270 nm and at the Mg I line at 285 nm. SATIRE-S does show, however, good agreement with SSI observations from the UARS satellite: the reconstructed SSI in the range 220-240 nm agrees well with UARS/SUSIM measurements (Krivova et al. 2006, 2009) and reasonably well with UARS/SOLSTICE measurements (Unruh et al. 2012).

To better reflect the spectral irradiance variability between 115 nm and 270 nm, we apply the empirical method outlined in Krivova et al. (2006). This method relies on the good agreement in the temporal variability of the 220 to 240 nm region as calculated by SATIRE-S and uses the scaling coefficients derived from spectral irradiance measurements, over the period 1997 to 2002, taken by the UARS/SUSIM instrument (Brueckner et al. 1993,

Floyd et al. 2003). Therefore, spectral regions in SATIRES below 220 nm and between 240 and 270 nm rely on SUSIM measurements and the close agreement in these regions is partly by design. In section c we show an example of this with the reconstructed Lyα irradiance, which is in agreement with the composite of Lyα measurements and proxy-models by Woods et al. (2000), and three integrated UV wavelength bands below 290 nm.

The TSI data set that is obtained by integrating our new SSI reconstruction is considered as an update of the TSI reconstruction presented in Ball et al. (2012). While both reconstructions are consistent within their uncertainty ranges, the updated TSI is now based upon the integral of the SSI that is self-consistent at every wavelength for the full reconstruction period. In Ball et al. (2012), the intercycle decline between 1996 and 2008 was estimated to be 0.20 +0 12 . -0 09 . Wm -2 , where the errors are one-sigma uncertainties. The new reconstruction revises this estimate down to 0.13 +0 07 . -0 10 . Wm -2 . Note that the reconstructions are calibrated using only the SORCE/TIM measurements and are thus independent of any TSI composite post-1990 and independent within the uncertainty range prior to this period (see Ball et al., 2012, for more details).

In a final step we adjust the absolute levels of SATIRE-S SSI so that the integrated SSI is in agreement with SORCE/ TIM at the solar minimum in December 2008. For this, the entire spectrum is multiplied by a factor of 1.0047. This small correction of 0.5% assures that the original variability as obtained directly from SATIRE-S, is not affected.

## b. Data gap filling

The new SATIRE-S SSI reconstruction now extends through the most recent and unusually long solar minimum period, whereas the previous version (Krivova et al. 2009) ended in 2007. For the period between 1974 December 10 and 2009 October 31, images are missing on ∼ 50% of dates, mostly within cycles 21 and 22; we fill these data gaps to provide fluxes on all dates over the entire period. To avoid any assumptions about the solar behaviour, we decompose each wavelength into short-term, or rotational, and longterm timeseries. The long-term timeseries is obtained by smoothing the original timeseries using a Gaussian window equivalent to a boxcar width of 135 days. We use this period, longer than the typical 81-days, to reduce the impact of short-term variability. The short-term timeseries, which captures rotational variability, is obtained by subtracting the long-term timeseries from the original.

Gaps in the long-term SATIRE-S time-series are filled by linear interpolation. Most gaps are short, with 90% of gaps being no longer than a solar rotation of 27 days in length, so the long-term trend is well approximated by a linear interpolation. Only five data gaps exceed two solar rotations, the longest of which is a 282 day period around the solar minimum of 1976.

Gaps in the detrended, rotational time series are filled using solar activity indices: the NOAA and LASP Mg II indices (Viereck et al. 2004; Snow et al. 2005), combined through linear regression; the Lyα composite by Woods et al. (2000); the Penticton F10.7 cm radio flux (data available through the National Geophysical Data Center at http://www.ngdc.noaa.gov/); the TSI from version d41 62 1003 of the PMOD composite (Fr¨hlich 2000) and o the sunspot area composite record by Balmaceda et al. (2009). Each index is indicative of the behaviour of some feature in the solar atmosphere, although it is not clear exactly how they relate to each wavelength of SSI (see Dudok de Wit et al., 2008, and supplementary material section b). Rotational variability at each wavelength is better approximated using multi-linear regression of two indices than using just one. We calculated the regression coefficient for every combination of two indices for each wavelength using dates when all indices and reconstructed SSI exist. Then, for each data gap, the available index-pair with the highest value of the coefficient of determination, r 2 c , at each wavelength is used to calculate the SSI in the missing gap.

Finally, the detrended and smoothed time series are added together to produce a spectral reconstruction that reflects the long-term variability of SATIRE-S while retaining rotational consistency (see supplementary material section b for examples). This procedure is expected to perform less well prior to 1978 because the TSI and the Mg II index, which generally have the highest combined correlation coefficients, are unavailable then. Also, this time period coincides with the longest data gaps in the reconstruction.

The change in SSI between maximum activity in solar cycle 23, in 2002, and the minimum in 2008 is plotted in Fig. 1. The blue curve indicates the final SATIRE-S reconstruction. This figure is described in greater detail and discussed in section 3a.

Figure 1 - The (top) percentage and (bottom) absolute change in flux between cycle 23 maximum and minimum, in terms of the maximum change in irradiance between three-month averages at February 2002 and December 2008 (see Fig. 2), for the SATIRE-S (blue) and NRLSSI (red) models; in the upper plot solid lines represent a decrease in flux while dotted lines represent an increase. The vertical dashed lines indicate 242 and 310 nm and the horizontal line marks no change in flux. (inset) The absolute change in flux in 10 nm bands between three month averages at May 1996 and February 2002 for SATIRE-S (blue), NRLSSI (red) and UARS/SUSIM (purple).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0365v1-figure-acfb8fb3-b029-42ea-a080-c07574fafcaf.jpg)

## c. Uncertainty estimate

An accurate error estimate for the modelled reconstruction is difficult to provide, since it depends partly on unknowns (such as the amount of magnetic flux missed by the magnetograms) and on uncertainties that cannot be precisely constrained within the scope of this paper (such as the accuracy of the model atmospheres employed, or the influence of neglecting non-LTE effects; see section c of the supplementary materials). Therefore, we attempt to provide a long-term SSI uncertainty range similar to the approach taken for TSI in Ball et al. (2012). This is an empirical approach that takes into account the uncertainties introduced in the calibration steps described in section 2a. Specifically, we account for step (i), the regression fitting between the TSI derived from MDI images and the SORCE/TIM measurements, step (ii), the regression fitting to combine the SATIRE-S reconstructions for MDI and KP data, and step (iii), the uncertainties in the correction factor for the KP/512 relative to the KP/SPM magnetograms.

For wavelengths below 270 nm we add, in quadrature, the uncertainty from the SATIRE-S reconstruction and the estimated relative uncertainty of the UARS/SUSIM measurements, which are estimated to be of the order of 5% below 142 nm and decreases to about 2% above 160 nm (Woods et al., 1996, and Linton Floyd, personal communication); we linearly interpolate the uncertainty between these wavelengths. These uncertainties are provided with the published reconstruction. We also flag a few wavelengths in the SATIRE-S spectrum, most notably the Mg I line at 285 nm, where our detailed comparisons with SORCE/SIM measurements on rotational time scales indicate that SATIRE-S overestimates the solar variability; see section c of the supplementary materials.

To summarise and illustrate the temporal behaviour of the uncertainties, we list the cycle amplitudes and their associated uncertainties in Tab. ?? for the TSI and for selected 'broadband' spectral irradiances. Column 4 lists the cycle amplitudes for cycle 23; these are based on reconstructions from SoHO/MDI images, i.e. these account for uncertainties in step (i) only. As shown in Ball et al. (2012), the error on the regression fitting for the SoHO/MDI free parameter is small and arises mainly from the longterm uncertainty in SORCE/

TIM. The uncertainty on the cycle amplitude (between sunspot maximum in March 2000 to the minimum in December 2008), is of the order of 100 ppm, comparable to the long-term uncertainty of SORCE/TIM. We note that the agreement between the TSI derived with SATIRE-S and the available composites is just as good as the agreement between the different composites (Ball et al. 2012).

Going back in time, the uncertainties increase, mainly due to the additional calibration steps (ii) and (iii). This is illustrated by the larger uncertainties for the amplitudes of cycles 21 and 22 (see columns 2 and 3 of Tab. ?? , respectively). As illustrated in Fig. 2 and Fig. 7 in Ball et al. (2012), the uncertainties are asymmetric and typically show slightly larger positive ranges. This is due to the different response of the magnetograms when connecting reconstructions from MDI and KP, as in step (ii), and the uncertainties in the correction factor from step (iii) that only lead to an increase in flux variability, not a decrease (see Ball et al., 2012).

While considerable progress has been made in determining the absolute value of the total solar irradiance (Kopp and Lean 2011), the absolute spectral solar irradiance is still poorly constrained and a number of different 'standard' absolute solar spectra are available (see Thuillier et al. (2003) for a discussion of this). For this reason, the uncertainties listed here and distributed with the reconstructions are for relative irradiances. Relative spectral irradiances in each band are much better constrained than the absolute accuracy, though the degradation of space instruments means that considerable uncertainties remain when going beyond rotational time scales (Unruh et al. 2012; Ermolli et al. 2013). If users of the SATIRE-S reconstruction wish to use a different absolute spectra as a basis upon which SATIRE-S variability is placed, we provide absolute spectra binned onto the SATIRES wavelength grid. We do this for the ATLAS 3 and Whole Heliosphere Interval solar reference spectra (Woods et al. 2009).

The SATIRE-S SSI data can be found at http://www. mps.mpg.de/projects/sun-climate/data.html.

## 3. Inter-comparison of SATIRE-S with other datasets

In this section, we compare the SATIRE-S model with the NRLSSI model (Lean 2000; Lean et al. 2005) and with SORCE/SOLSTICE observations. NRLSSI is the most widely used, empirically derived, model of SSI, so we perform a comparison with SATIRE-S here. The solar cycle spectral variability recorded by SORCE/SOLSTICE is larger than the variability seen by the UARS/SUSIM and UARS/SOLSTICE instruments. It is also larger than the variability inferred from the models. SORCE/SIM displays even larger variability than SORCE/SOLSTICE, for the overlapping range between 240-310 nm, and it is an interesting SSI dataset to consider. However, a comparison of SATIRE-S and an updated version of SORCE/SIM version 17 was made in Ball et al. (2011). The SORCE/SIM UV changes below 310 nm were up to five times larger than the changes seen in SATIRE-S. We note that Figure 11 of Ball et al. (2011) indicates that SORCE/SIM displays up to 10 times the variability of SATIRE-S at some wavelengths between 300 and 400 nm, though the integrated

Table 1 - Amplitude of flux variability over a solar cycle (SC) for selected wavelength bands and for total solar irradiance (TSI). The uncertainty of the amplitude is also given.

| λλ , nm   | SC21, mWm - 2                                                       | SC22, mWm - 2   | SC23, mWm - 2   |
|-----------|---------------------------------------------------------------------|-----------------|-----------------|
| 200-270   | 114 +25 - 8 509 +127 - 46 241 +141 - 72 138 +67 - 34 979 +394 - 183 | 114 +25 - 8     | 99 +4 - 4       |
| 270-400   | 114 +25 - 8 509 +127 - 46 241 +141 - 72 138 +67 - 34 979 +394 - 183 | 501 +135 - 47   | 445 +16 - 21    |
| 400-700   | 114 +25 - 8 509 +127 - 46 241 +141 - 72 138 +67 - 34 979 +394 - 183 | 190 +166 - 80   | 201 +37 - 51    |
| 700-1000  | 114 +25 - 8 509 +127 - 46 241 +141 - 72 138 +67 - 34 979 +394 - 183 | 116 +78 - 38    | 118 +18 - 24    |
| TSI       | 114 +25 - 8 509 +127 - 46 241 +141 - 72 138 +67 - 34 979 +394 - 183 | 876 +445 - 195  | 837 +85 - 113   |

Figure 2 - Smoothed (solid) wavelength-integrated time series of SATIRE-S (blue), NRLSSI (red) and the PMOD composite of TSI (black) between 1978 and 2010. PMOD and NRLSSI have been normalised to the minimum of SATIRE-S in 1986. The SATIRE-S uncertainty range is shown by light blue shading and the PMOD error bars are given at the cycle minima. Grey shading highlights the periods where gaps exceeding three solar rotations have been filled.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0365v1-figure-90b9a19a-88ee-478e-a525-166fb909294a.jpg)

variability of SORCE/SIM over this wavelength range was approximately 3.4 times larger than SATIRE-S. Ball et al. (2011) suggested that degradation of the instrument has not been properly accounted for and may be overestimating the long-term trends (see also Lean and DeLand (2012); Deland and Cebula (2012); Unruh et al. (2012); Ermolli et al. (2013)). SORCE/SIM is currently undergoing a reanalysis that may affect the cycle variability and its uncertainty estimates (Peter Pilewskie, personal communication). We, therefore, only consider SORCE/SOLSTICE in the following. We compare SATIRE-S and NRLSSI for wavelengths between 120 and 3000 nm and up to 310 nm when comparing with SORCE/SOLSTICE.

scribe the evolution of sunspots and faculae, respectively. For wavelengths &lt; 400 nm, spectral irradiances are computed from multiple regression analysis with UARS/ detrended, rotational data to avoid instrumental degradawavelengths above 400 nm, facular and sunspot contrasts, from the models by Solanki and Unruh (1998), are scaled

SOLSTICE observations. This analysis is performed on tion effects and, therefore, assumes that rotational variability scales with solar cycle changes in irradiance. For to agree with solar cycle TSI observations. NRLSSI's integrated flux is ∼ 4 Wm -2 higher than SORCE/TIM so we normalise NRLSSI to SORCE/TIM, as was done for SATIRE-S in section 2a.

## a. Comparison with the NRLSSI model

NRLSSI is an empirical model that uses the diskintegrated Mg II and photospheric sunspot indices to de-

It is worth briefly considering how the TSI derived by integrating the SSI differs between the two models. In Fig. 2, the smoothed, wavelength-integrated SATIRE-S (blue), NRLSSI (red) and PMOD composite of TSI (black) (Fr¨hlich 2003) are plotted between 1978 and 2009. NRLSSI o

and PMOD have been normalised to the absolute value of SATIRE-S averaged over three months centred at the minimum of 1986. The SATIRE-S uncertainty range is plotted with light blue shading and the cycle minima error bars from PMOD are plotted as black bars (Fr¨hlich 2009). Alo though other TSI composites exist, PMOD is now generally accepted as the most accurate TSI composite of observations, which is why we consider only it here; see Ball et al. (2012) for comparisons with all TSI composites on rotational and cyclical timescales. We find correlation coefficients between NRLSSI and SATIRE-S TSI with PMOD, of 0.92 and 0.96, respectively; detrending the timeseries, as described in section 2b, yields correlation coefficients of 0.87 and 0.96 for the rotational variability. These statistics suggests that SATIRE-S reproduces the PMOD composite of observations better than NRLSSI. However, there are periods when NRLSSI matches PMOD better, on yearly and longer timescales, than SATIRE-S. For example, the large difference between SATIRE-S and PMOD around 19911993 is due to the remaining uncertainties in the crosscalibration of the KP 512 and SPM magnetograms (see also Wenzler et al., 2006, and Ball et al., 2012).

The inter-cycle trends of the three datasets are subtly different. Whereas PMOD and SATIRE-S show a decline of ∼ 0.20 +0 16 . -0 26 . and 0.13 +0 07 . -0 10 . Wm -2 , respectively, between 1996 and 2008, NRLSSI exhibits no change over this period. NRLSSI's behaviour is most likely due to the use of the Mg II index as a proxy for long-term changes; this index does not exhibit any strong inter-cycle variation (Fr¨hlich o 2009). Note, however, that the uncertainty of the cycle minima in PMOD and SATIRE-S also encompasses the NRLSSI model estimate of no change, though the Mg II record is not entirely free of long-term uncertainty either (Snow et al., 2005; Marty Snow, personal communication). We note that while accurate TSI is important to act as a constraint for the SSI in both NRLSSI and SATIRE-S (and SORCE/SIM, which is much less well constrained, see Ball et al. (2011)), it does not ensure that the SSI is correct, in either case, as higher variability at some wavelengths can be compensated by lower variability at other wavelengths.

For the SSI comparison between NRLSSI and SATIRES, in Fig. 1 we consider the change in flux, ∆F, between two 81-day averaged periods centred on: 2002 February 1, the second and highest peak of cycle 23; and the cycle 23/24 minimum on 2008 December 15. This provides the largest range of change in cycle 23.

In the upper plot of Fig. 1 the percentage change between the maximum and minimum of solar cycle 23 is plotted on a logarithmic scale while the lower plot depicts the absolute change in flux on a linear y-axis. The spectral uncertainty in SATIRE-S is very small for cycle 23 and is virtually invisible on the plotted scales. We, therefore, plot the spectral variability and the uncertainties for each cycle in section a of the supplementary materials. The wavelength spacings in the two models are different. Consequently, NRLSSI has been interpolated onto the SATIRES wavelength grid (see section 2). The regions below 242 nm and between 242 and 310 nm are important in ozone production and destruction processes in the stratosphere (see section 4), so these have been highlighted with vertical dashed lines.

Below 242 nm, the two models agree well in the change of flux, with larger differences apparent only below ∼ 150 nm; these result from the use of different instruments to define the cycle variability, i.e. UARS/SUSIM for SATIRES (see section 2c) and the scaled rotational variability of UARS/SOLSTICE for NRLSSI (see above). Integrating over 120-242 nm, a region important for the photodissociation of O 2 and O , we find that NRLSSI shows the same 3 solar cycle change as SATIRE-S, though we note that the SATIRE-S Lymanα has ∼ 6% larger cycle amplitude than NRLSSI during the descending phase of cycle 23, but which is generally larger for NRLSSI in earlier solar cycles (see section c and Fig. 4). These differences in Lymanα cycle changes will impact on OH chemistry that indirectly affect ozone concentrations. Between 242 and 310 nm, important in O 3 photodissociation, the integrated change in flux, ∆F, is more than 50% larger in SATIRE-S than NRLSSI.

SATIRE-S shows up to three times larger cycle variability than NRLSSI for wavelengths between 300 and 400 nm. At longer wavelengths, between 400 and 1250 nm, i.e. in the visible and near-IR, NRLSSI generally displays larger cycle variability than SATIRE-S. At yet longer wavelengths, both NRLSSI and SATIRE-S display negative variability in the IR, though this occurs over a wider range and with higher variability in SATIRE-S than NRLSSI, the latter of which only shows negative variability between 1500 and 1850 nm. Where this transition to negative variability occurs is highly dependent on the facular model used (see supplementary materials section c and Unruh et al. (1999, 2000, 2008)). We note that negative variability in the IR region is also seen in SORCE/SIM in the IR above 970 nm, but this change is much larger than both NRLSSI and SATIRE-S. In Ball et al. (2011), the integrated IR region of 972-1630 nm shows ∼ 0.00 Wm -2 change in SATIRE-S between 2004 and 2008, while SORCE/SIM increases by 0.24 Wm -2 . For all wavelengths between 200 and 1600 nm, detailed comparisons between SORCE/SIM and the SATIRE-S and NRLSSI models are made by Ball et al. (2011) and Lean and DeLand (2012), respectively.

## b. Comparison with SORCE/SOLSTICE

SORCE/SOLSTICE observations cover the UV region 115-310 nm. In Fig. 3, we compare the modelled SSI with SORCE/SOLSTICE between two 81-day average periods centred on 2003 August 15 and 2008 December 2008. This period constitutes ∼ 60% of the full cycle variation in TSI

and allows for a direct comparison with SORCE/SOLSTICE observations that started on 2003 May 14, i.e. some time after the maximum of cycle 23. The left plot of Fig. 3 shows the absolute change in flux, while the right is the change in absolute flux relative to SATIRE-S. We plot version 10 of SORCE/SOLSTICE data in addition to the latest version 12 as we consider version 10 in the next section. SORCE/SOLSTICEdata have a reported long-term uncertainty of 0.5% per year which amounts to ∼ 2.7% for the period considered in Fig. 3 and is shown with grey shading for version 12 only; version 12 data have additional degradation corrections to prior versions that result in a change to the absolute level and relative change in spectral irradiance (Marty Snow, personal communication). When not otherwise specified, version 12 is referred to in the following.

From Fig. 3, SORCE/SOLSTICE cycle variability is in reasonable agreement with the models below 140 nm and near 165, 180 and 255 nm. At most other wavelengths below 290 nm, cycle variability is larger in SORCE/SOLSTICE than in the models, typically by at least a factor of two: it is twice that of SATIRE-S at ∼ 170 and 200 nm; between 200 and 250 nm SORCE/SOLSTICE increases from two up to five times the change in SATIRE-S; and between 260 and 290 nm it typically ranges between 1.7 and 3.5 times larger than SATIRE-S. Above 290 nm the long-term uncertainty of SORCE/SOLSTICE exceeds the flux change, ∆F(Marty Snow, personal communication), so results cannot be considered realistic or useful in comparison to the models.

## c. Timeseries comparison with NRLSSI and SORCE/SOLSTICE

We show examples of timeseries for three different wavelength bands: Lymanα at ∼ 121 nm, 176-242 nm and 242290 nm.

In Fig. 4, we compare the Lymanα composite by Woods et al. (2000) (black) with SATIRE-S (blue), NRLSSI (red) and SORCE/SOLSTICE version 12 (green) and 10 (yellow, dashed). In this figure, SATIRE-S and NRLSSI are normalised to the Lymanα composite for the solar minimum of 1986 (as was done for the TSI plot with PMOD in Fig. 2), while both versions of SORCE/SOLSTICE are normalised to the three month averaged period around December 2008. The grey shading is an estimated 10% uncertainty (Woods et al. 2000) on the cycle amplitude, or approximately 0.22 mWm -2 at the 1 σ level. The Lymanα composite agrees with SORCE/SOLSTICE v10 almost exactly during the declining phase of cycle 23 because that version was used in the Lymanα composite. Light blue shading is the uncertainty range of SATIRE-S at 121.5 nm and is the quadrature sum of uncertainty from the SATIRE-S reconstruction and an estimated 5% uncertainty from UARS/SUSIM for the correction applied using the method by Krivova et al. (2006).

The nominal values of NRLSSI and SATIRE-S are at the lower end of the Lymanα uncertainty range. Considering that there is uncertainty at the cycle minima and maxima, SATIRE-S and NRLSSI absolute values could be shifted up and remain within the Lymanα composite uncertainty in all three solar cycles. Therefore, both models agree with the Lymanα observations, within the uncertainty, at almost all times on annual to decadal time scales.

UARS/SUSIM monitored spectral irradiance between

115 and 410 nm for the period 1991 to 2005 and suggested lower solar cycle changes in UV SSI than SORCE/ SOLSTICE. In Fig. 5, we show timeseries for two integrated regions, between 176-242 nm and 242-290 nm (important in the production and destruction of stratospheric ozone, respectively), for SATIRE-S (blue), NRLSSI (red), SORCE/SOLSTICE version 10 (yellow) and version 12 (green) and also UARS/SUSIM (purple). The absolute fluxes of NRLSSI and UARS/SUSIM are shifted to SATIRES for the average of the period between 1997 and 2000 and SORCE/SOLSTICE are shifted to SATIRE-S by the average over the period between June and November 2003. Both wavelength bands show similar solar cycle behaviour for UARS/SUSIM and the two models; SORCE/SOLSTICE show much larger solar cycle trends. We note that the degradation correction issue discussed by Krivova et al. (2006) can be seen just prior to the solar minimum of 1996. This would lead to a change of the cycle maximum in 1991 relative to that of 2000.

However, the two bands in Fig. 5 are broad spectral bands: selecting different wavelength ranges might lead to different conclusions as to which model is better at reconstructing SSI solar cycle changes based on UARS/SUSIM measurements. The inset in the lower plot of Fig. 1 shows the solar cycle change in flux from three month averages, centred on the solar minimum in 1996 and solar maximum in 2002, in 10 nm bands between 120 and 400 nm. This plot illustrates the wavelength dependence of the solar cycle changes. All three datasets show good agreement below 250 nm. Above 250 nm, UARS/SUSIM and SATIRE-S typically show larger solar cycle changes than NRLSSI; SATIRE-S shows closer solar cycle changes to UARS/SUSIM than NRLSSI.

Understanding the effect that the Sun has on Earth's atmosphere and climate requires reliable estimates of SSI variability. The differences between SATIRE-S, NRLSSI and SORCE/SOLSTICE lead to different predicted impacts on the Earth's atmosphere, e.g. heating and photodissociation rates, as the following example considering changes in ozone concentration shows.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0365v1-figure-a8098c57-a498-4621-b179-0e34f224e4db.jpg)

Figure 3 - The SSI change between two 81-day periods centred on 2003 August 15 and 2008 December 15 (left) in absolute flux and (right) relative to SATIRE-S, on a logarithmic scale. Plotted are SORCE/SOLSTICE version 10 (orange) and 12 (green), SATIRE-S (blue) and NRLSSI (red). Dotted lines indicate an inverse trend in flux between 2003 and 2008; flux variations at wavelengths longer than 290 nm cannot be considered reliable as the long-term uncertainty exceeds the flux change, ∆F, at these wavelengths (see main text). The uncertainty range for SORCE/SOLSTICE version 12 is shown in grey shading for wavelengths below 290 nm. Wavelengths are 1 nm width bins below 290 nm and 2 nm bins above. Dashed lines indicate the wavelength 242 nm.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0365v1-figure-c9af1fde-f2ff-4d77-8f9d-8b0e11c6ea0a.jpg)

## 4. Modelled ∆ O from different input SSI 3

We employ the radiative-chemical transport atmospheric model based on that of Harwood and Pyle (1975) to compare the resultant change in ozone concentration, ∆O , 3 when using modelled and observational SSI data as a solar input. The atmospheric model has been used in many studies to investigate the dynamics and chemical interactions of the Earth's atmosphere (e.g. Bekki et al., 1996, Warwick et al., 2004, Haigh et al., 2010). It calculates the zonal mean temperature and winds and chemical-constituent concentration in a time-dependent 2D model with full radiative-chemical-dynamical coupling. The model considers a spherical Earth with seasons, but without land topography or oceans. The model is resolved into 19 latitudes and 29 pressure levels up to 95 km. In this study the only difference between model runs is the specification of the input SSI. Use of different input SSI affects photochemical reactions involved in ozone production and destruction which account for the largest contribution to heating in the stratosphere, through the absorption of solar flux below 310 nm. The spectral resolution of the input spectra for the atmospheric model decreases non-linearly from less than 1 nm at 116 nm to 5 nm at 300 nm; it remains at 5 nm resolution up to 650 nm, after which it has 10 nm resolution up to 730 nm.

The atmospheric model requires input SSI up to 730

nm. As SORCE/SOLSTICE data are not available above 310 nm, SATIRE-S fluxes are used. This is valid as there is little effect on ∆O 3 profiles from wavelengths above 310 nm. SORCE/SOLSTICE data show anomalous, large inverse trends in ∆F between 298 and 310 nm, indicated by the dotted lines in Fig. 3, but tests show that this narrow region also has very little impact on the magnitude and spatial distribution of ∆O 3 ; the results remain effectively unchanged if we use SATIRE-S at and above 298 nm, so we use the full wavelength range of SORCE/SOLSTICE in runs using SOLSTICE.

Tests show that the background spectrum (absolute flux) upon which the changes in SSI are superimposed makes very little difference to the O 3 results which depend more critically on the spectral shape of the change. This is because of the strong wavelength dependence of ozone production/destruction reactions. We therefore make no attempt to adjust fluxes to a common absolute level and we use the SSI data as published (see also section d of the supplementary material).

Using the same model that we use here, Haigh et al. (2010) investigated middle atmosphere ozone changes for SSI changes between 2004 and 2007 for NRLSSI and hybrid SORCEspectra with SIM data for λ &gt; 200 nm and an older version of SOLSTICE below. Merkel et al. (2011) presented similar work using the WACCM model, but switched between an older version of SORCE/SOLSTICE and SORCE/

Figure 4 - The Lymanα composite (black) by Woods et al. (2000) with an uncertainty of 10% of the cycle amplitude (grey shading) compared with SATIRE-S (blue) with the uncertainty range (light blue shading), the NRLSSI model (red) and SORCE/SOLSTICE version 12 (green) and 10 (yellow). SATIRE-S and NRLSSI are normalised to the composite at the solar minimum of 1986 and SORCE/SOLSTICE datasets are normalised to the composite at the solar minimum of December 2008.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0365v1-figure-c77d8455-ff52-4e6c-bb7c-cd24404ec921.jpg)

SIM datasets at 240 nm. The choice of wavelength at which to make the switch is somewhat arbitrary and makes a direct comparison with the runs we do here difficult. Also, the exact choice of dates over which SSI are averaged is different between Haigh et al. (2010) and Merkel et al. (2011). We do not try to reproduce their model setups exactly. Instead, we focus on making comparisons between the models and observations we present here. However, it is worth noting that the spatial distributions of ∆O 3 when using version 10 of SORCE/SOLSTICE show similar structure and magnitude to the results presented by Haigh et al. (2010) and Merkel et al. (2011). Both studies, when employing SORCE data between 2004 and 2007, found a negative response (i.e. an increasing ozone outof-phase with solar irradiance) in the mesosphere around 55 km, of 1.2% and 2% respectively. They also found a positive response in the middle stratosphere below 40 km, reaching ∼ 2%. Using NRLSSI in these studies resulted in a response in-phase with SSI changes at all altitudes between 30 and 60 km of up to ∼ 1%.

The input UV SSI for the atmospheric model are the 81-day average spectra as presented in Fig. 3. In Fig. 6 the change in ozone concentration produced by taking the difference between atmospheric model outputs that use SSI from 2003 and 2008 is shown for (a) NRLSSI, (b) SATIRES, (c) SORCE/SOLSTICE v10 and (d) SORCE/SOLSTICE v12. The period considered is during a decline in TSI and UV fluxes that is approximately 50% larger than for the 2004 to 2007 periods used by Haigh et al. (2010) and Merkel et al. (2011). In Fig. 6a the result for NRLSSI is qualitatively similar in spatial distribution to the NRLSSI result of Haigh et al. (2010), but values are approximately 50% larger, as expected for the larger UV change. Figure 6c, using SORCE/SOLSTICE version 10 also displays

Figure 5 - Timeseries of SATIRE-S (blue), NRLSSI (red), SORCE/SOLSTICE version 10 (yellow) and version 12 (green) and UARS/SUSIM (purple) are shown for the wavelength bands (top) 176-242 nm and (bottom) 242-290 nm between 1991 and 2009. UARS/SUSIM and NRLSSI are shifted to the mean of SATIRE-S for the period between 1997 and 2000 and SORCE/SOLSTICE v10 and v12 between June and November 2003.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0365v1-figure-ce146c5a-3cc8-4fe0-9697-b490efb2a4b5.jpg)

qualitatively similar ∆O 3 to those using SORCE presented by Haigh et al. (2010) and Merkel et al. (2011). The general spatial distribution of ∆O 3 in both versions of SORCE/ SOLSTICE presented in Figs. 6c and 6d are similar. However, there are two striking differences: the magnitude of ∆O 3 in the equatorial lower mesosphere is reduced by a factor of six at around 55 km, from -1.6% to -0.2%, for version 10 and 12, respectively; and the zero-line has also shifted up by ∼ 5 km in version 12 with an increased ∆O 3 maximum in the stratosphere at around 40 km. The change in ∆O 3 from version 10 to 12 reflects a general decrease of ∆F at wavelengths &gt; 242 nm and competing increases and decreases in ∆F at wavelengths &lt; 242 nm.

annum. So, even though there are these differences, in the ∆Fs of the two SORCE/SOLSTICE versions, they are within the stated uncertainty. Our work shows that knowledge of the detailed spectra to better than this accuracy is required to place any faith in derived atmospheric effects.

Despite showing different responses in the atmospheric model for these different datasets, version 10 and 12 of SORCE/SOLSTICE only differ at a small number of wavelengths by more than its stated uncertainty of 0.5% per

The run using SATIRE-S, in Fig. 6b, also displays a negative ∆O 3 above ∼ 50 km, though of lower magnitude than for SORCE/SOLSTICE. Relative to NRLSSI, the negative response in the mesosphere results mainly from the larger ∆F in SATIRE-S above 242 nm. The total concentration of O+O 3 is similar in SATIRE-S and NRLSSI due to similar ∆F for wavelengths &lt; 242 nm. UV radiation at wavelengths between 242 and 310 nm photodissociates ozone to produce O( D) and O . 1 2 This has two effects: it reduces the O /O ratio and O( D) reacts with H O to pro3 1 2 duce OH which catalytically destroys O . 3 Both processes tend to decrease the O 3 concentration at these altitudes.

Figure 6 - Contour plots of the change in stratospheric ozone, ∆O , between 2003 and 2008 using SSI input from (a) 3 the NRLSSI model, (b) the SATIRE-S model, (c) SORCE/SOLSTICE version 10 and (d) version 12. All results are for December 25. Contours are given in 0.1% levels in (a) and (b) and 0.2% in (c) and (d) with negative changes, i.e. increases with decreasing solar activity, are shown with dotted contours. Blue and purple reflect negative changes, while green, yellow and red show positive changes.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0365v1-figure-d23b68e2-6a5d-42f1-8817-00b8b8ecf6c7.jpg)

Although NRLSSI and SATIRE-S show similar cycle variability below 242 nm, the reduction in ozone concentration at all altitudes and negative response in the mesosphere of SATIRE-S, with respect to NRLSSI, is a result of the larger flux change in SATIRE-S at wavelengths longer than 242 nm. SORCE/SOLSTICE shows very different flux changes compared to SATIRE-S and NRLSSI, but the similar response in the mesosphere relative to SATIRE-S is caused by competing effects on ozone concentration at all wavelengths.

ties that span the range of differences shown by the model runs here. It is therefore not possible, with datasets and methods currently available, to make unequivocal statements about the 'true' ozone response to solar cycle changes. The absolute calibration is less important than the longterm trends in spectral irradiance, the accuracy of which is limited by the stability of the instrument sensitivity. Accurate cycle trends in spectral irradiance are crucial to have confidence in the solar effect on the Earth. At the moment this is lacking.

The modelling study by Swartz et al. (2012) shows similar responses in O 3 to the NRLSSI and SORCE spectra as Haigh et al. (2010) and Merkel et al. (2011). Analysis of observational ozone data, e.g. SBUV, SAGE and HALOE

from Soukharev and Hood (2006), Aura/MLS from

Haigh et al. (2010) and TIMED/SABER from Merkel et al. (2011), suggests a range of changes in ozone with uncertain-

## 5. Discussion and conclusions

In this paper, we present SATIRE-S spectral solar irradiances covering fully the last three solar cycles spanning 1974 to 2009. SATIRE-S is a semi-empirical model that assumes all irradiance changes are the result of changes in surface magnetic flux (Fligge et al. 2000; Krivova et al.

2003). Data are available for all dates over the period with an accompanying error estimate.

We compare our new SATIRE-S spectral irradiances with the NRLSSI model and SORCE/ SOLSTICE observations. NRLSSI and SATIRE-S are in good agreement between ∼ 150-242 nm. Between 242 and 400 nm, the solar cycle change in flux in SATIRE-S is generally ∼ 50% larger than NRLSSI. Both models display significantly lower cycle variability than SORCE/SOLSTICE at almost all wavelengths between 190 and 290 nm.

Both models provide consistent, long-term reconstructions of solar irradiance with different cycle variability. At the moment there is insufficient reliable observational spectral irradiance data to identify which model is more accurate (on solar cycle time scales). While NRLSSI uses solar indices, SATIRE-S provides a direct translation from observed solar images to irradiances. The results presented in this paper show that SATIRE-S reproduces the PMOD TSI composite values better than NRLSSI (though the uncertainties of SATIRE-S and PMOD both encompass NRLSSI). Above 250 nm, the SATIRE-S spectral response is closer to that shown by UARS/SUSIM than NRLSSI (see inset of Fig 1). We suggest that the SATIRE-S spectral irradiances should be considered in future climate studies, either on its own or in addition to NRLSSI.

We present an example of the physical implications that these different SSI datasets have within the stratosphere. We calculate the ozone concentration resulting from the use of the three datasets within a 2D atmospheric model. We find that: (a) the magnitude of mesospheric ∆O 3 response when using SORCE/SOLSTICE versions 10 and 12 is very different; (b) SATIRE-S mimics the small negative mesospheric solar cycle change that SORCE/SOLSTICE version 12 produces; and (c) NRLSSI displays positive changes at all heights.

Interestingly, a recent study by Wang et al. (2013) suggests that modelled solar cycle changes in OH concentration in the stratosphere agree better with observations when using SORCE than NRLSSI data as a model input. We note that Wang et al. (2013) used hybrid SORCE datasets with SIM data either above 240 or above 210 nm and SOLSTICE below these wavelengths; this is different to the spectra we use here or that used by Haigh et al. (2010) and Merkel et al. (2011). However, the different ozone changes in the mesosphere resulting from SSI inputs from two different models and two versions of the same observational dataset in this study, highlights the need for a better understanding of solar cycle changes at wavelengths important for both O , O 2 3 and other atmospheric photochemistry such as OH as in Wang et al. (2013). If the solar effect on ozone is to be isolated correctly, it is imperative that greater certainty is established in SSI variability.

This raises an important point that should be considered when making comparisons of observed with modelled ozone changes which use different SSI datasets. Data are usually updated to make an improvement on previous releases. Unfortunately, the fast pace at which data revisions are made makes comparisons with previous publications difficult. In the case presented here, the change in ∆O 3 between version 10 and 12 of SORCE/SOLSTICE leads to a result in the mesosphere that is so different that it is not yet possible to make robust conclusions about SSI or ozone based on the negative ozone mesospheric response. Many previous investigations have used different SORCE/SIM spectral data, and hybrids with SORCE/SOLSTICE, to investigate the atmospheric response. This makes direct comparisons with the results published in the literature difficult. We have focused the investigation to highlight what is different about SATIRE-S and to only consider the effect of one SORCE dataset in each run. A major result is that significant uncertainties remain in how SSI affects the atmosphere on time scales longer than the solar rotation.

The SATIRE-S SSI data are available for download at http://www.mps.mpg.de/projects/sun-climate/data.html. The SSI reconstruction ends in October 2009. Work is currently underway to extend the reconstruction using fulldisk images from the Solar Dynamics Observatory Helioseismic and Magnetic Imager (Schou et al. 2012) and will be made available for download in the future.

## Acknowledgments.

We thank Marty Snow, Linton Floyd, Tom Woods, Peter Pilewskie and Judith Lean for helpful discussions. We acknowledge the use of the Ottawa/Penticton 2800 MHz Solar Radio Flux and version d41 62 1003 of the PMOD dataset from PMOD/WRC, Davos, Switzerland and the unpublished data from the VIRGO Experiment on SoHO, a project of international cooperation between ESA and NASA. This work has been partly supported by the NERC SolCli consortium grant, the STFC grant ST/I001972/1, the WCU grant number R31-10016 funded by the Korean Ministry of Education, Science and Technology and FP7 SOLID.

## REFERENCES

Austin, J., et al., 2008: Coupled chemistry climate model simulations of the solar cycle in ozone and temperature. Journal of Geophysical Research (Atmospheres) , 113 (D12) , D11306, doi:10.1029/2007JD009391.

Ball, W. T., Y. C. Unruh, N. A. Krivova, S. Solanki, and J. W. Harder, 2011: Solar irradiance variability: A sixyear comparison between SORCE observations and the SATIRE model. Astronomy &amp; Astrophysics , 530 , A71, doi:10.1051/0004-6361/201016189.

Ball, W. T., Y. C. Unruh, N. A. Krivova, S. Solanki, T. Wenzler, D. J. Mortlock, and A. H. Jaffe, 2012: Reconstruction of total solar irradiance 1974-2009. Astronomy &amp; Astrophysics , 541 , A27, doi:10.1051/0004-6361/ 201118702, 1202.3554 .

Balmaceda, L. A., S. K. Solanki, N. A. Krivova, and S. Foster, 2009: A homogeneous database of sunspot areas covering more than 130 years. Journal of Geophysical Research (Space Physics) , 114 (A13) , A07104, doi: 10.1029/2009JA014299, 0906.0942 .

Bekki, S., J. A. Pyle, W. Zhong, R. Toumi, J. D. Haigh, and D. M. Pyle, 1996: The role of microphysical and chemical processes in prolonging the climate forcing of the Toba eruption. Geophysical Research Letters , 23 , 2669-2672, doi:10.1029/96GL02088.

Brueckner, G. E., K. L. Edlow, L. E. Floyd, IV, J. L. Lean, and M. E. Vanhoosier, 1993: The solar ultraviolet spectral irradiance monitor (SUSIM) experiment on board the Upper Atmosphere Research Satellite (UARS). Journal of Geophysical Research , 98 , 10695, doi:10.1029/ 93JD00410.

Deland, M. T. and R. P. Cebula, 2012: Solar UV variations during the decline of Cycle 23. Journal of Atmospheric and Solar-Terrestrial Physics , 77 , 225-234, doi:10.1016/ j.jastp.2012.01.007.

Dhomse, S. S., et al., 2013: Stratospheric O 3 changes during 2001-2010: the small role of solar flux variations in a CTM. ACPD , 13 , 12263-12286, doi:10.5194/ acpd-13-12263-2013.

Dudok de Wit, T., M. Kretzschmar, J. Aboudarham, P.-O. Amblard, F. Auch` ere, and J. Lilensten, 2008: Which solar EUV indices are best for reconstructing the solar EUV irradiance? Advances in Space Research , 42 , 903-911, doi:10.1016/j.asr.2007.04.019, arXiv:astro-ph/0702053 .

Ermolli, I., et al., 2013: Recent variability of the solar spectral irradiance and its impact on climate modelling. ACP , 13 , 3945-3977, doi:10.5194/acp-13-3945-2013, 1303.5577 .

Fligge, M., S. K. Solanki, and Y. C. Unruh, 2000: Modelling irradiance variations from the surface distribution of the solar magnetic field. Astronomy &amp; Astrophysics , 353 , 380-388.

Floyd, L. E., J. W. Cook, L. C. Herring, and P. C. Crane, 2003: SUSIM'S 11-year observational record of the solar UV irradiance. Advances in Space Research , 31 , 21112120, doi:10.1016/S0273-1177(03)00148-0.

Fontenla, J. M., E. H. Avrett, and R. Loeser, 1993: Energy balance in the solar transition region. III - Helium emission in hydrostatic, constant-abundance models with diffusion. The Astrophysical Journal , 406 , 319-345, doi: 10.1086/172443.

Fr¨hlich, C., 2000: o Observations of Irradiance Variations. Space Science Reviews , 94 , 15-24.

Fr¨hlich, C., 2003: Long-term behaviour of space radiomeo ters. Metrologia , 40 , 60, doi:10.1088/0026-1394/40/1/ 314.

Fr¨hlich, C., 2006: Solar Irradiance Variability Since 1978. o Revision of the PMOD Composite during Solar Cycle 21. Space Science Reviews , 125 , 53-65, doi:10.1007/ s11214-006-9046-5.

Fr¨hlich, C., 2009: o Evidence of a long-term trend in total solar irradiance. Astronomy &amp; Astrophysics , 501 , L27L30, doi:10.1051/0004-6361/200912318.

Haigh, J. D., 2003: The effects of solar variability on the Earth's climate. Royal Society of London Philosophical Transactions Series A , 361 , 95, doi:10.1098/rsta.2002. 1111.

Haigh, J. D., A. R. Winning, R. Toumi, and J. W. Harder, 2010: An influence of solar spectral variations on radiative forcing of climate. Nature , 467 , 696-699.

Harder, J., G. Lawrence, J. Fontenla, G. Rottman, and T. Woods, 2005: The Spectral Irradiance Monitor: Scientific Requirements, Instrument Design, and Operation Modes. Solar Physics , 230 , 141-167, doi:10.1007/ s11207-005-5007-5.

Harder, J. W., J. M. Fontenla, P. Pilewskie, E. C. Richard, and T. N. Woods, 2009: Trends in solar spectral irradiance variability in the visible and infrared. Geophysical Research Letters , 36 , 7801, doi:10.1029/2008GL036797.

Harwood, R. S. and J. A. Pyle, 1975: A two-dimensional mean circulation model for the atmosphere below 80km. Quarterly Journal of the Royal Meteorological Society , 101 , 723-747, doi:10.1002/qj.49710143003.

Jones, H. P., T. L. Duvall, Jr., J. W. Harvey, C. T. Mahaffey, J. D. Schwitters, and J. E. Simmons, 1992: The NASA/NSO spectromagnetograph. Solar Physics , 139 , 211-232, doi:10.1007/BF00159149.

Kopp, G. and G. Lawrence, 2005: The Total Irradiance Monitor (TIM): Instrument Design. Solar Physics , 230 , 91-109, doi:10.1007/s11207-005-7446-4.

Kopp, G. and J. L. Lean, 2011: A new, lower value of total solar irradiance: Evidence and climate significance. Geophysical Research Letters , 38 , L01706, doi:10.1029/ 2010GL045777.

Krivova, N. A., L. Balmaceda, and S. K. Solanki, 2007: Reconstruction of solar total irradiance since 1700 from the surface magnetic flux. Astronomy &amp; Astrophysics , 467 , 335-346, doi:10.1051/0004-6361:20066725.

Krivova, N. A., S. K. Solanki, M. Fligge, and Y. C. Unruh, 2003: Reconstruction of solar irradiance variations in cycle 23: Is solar surface magnetism the cause? Astronomy &amp; Astrophysics , 399 , L1-L4, doi: 10.1051/0004-6361:20030029.

Krivova, N. A., S. K. Solanki, and L. Floyd, 2006: Reconstruction of solar UV irradiance in cycle 23. Astronomy &amp; Astrophysics , 452 , 631-639, doi:10.1051/0004-6361: 20064809.

Krivova, N. A., S. K. Solanki, and Y. C. Unruh, 2011: Towards a long-term record of solar total and spectral irradiance. Journal of Atmospheric and Solar-Terrestrial Physics , 73 , 223-234, doi:10.1016/j.jastp.2009.11.013.

Krivova, N. A., S. K. Solanki, T. Wenzler, and B. Podlipnik, 2009: Reconstruction of solar UV irradiance since 1974. Journal of Geophysical Research , 114 (D13) , doi: 10.1029/2009JD012375, 0907.1500 .

Krivova, N. A., L. E. A. Vieira, and S. K. Solanki, 2010: Reconstruction of solar spectral irradiance since the Maunder minimum. Journal of Geophysical Research (Space Physics) , 115 (A14) , A12112, doi:10. 1029/2010JA015431.

Kurucz, R., 1993: ATLAS9 Stellar Atmosphere Programs and 2 km/s grid. ATLAS9 Stellar Atmosphere Programs and 2 km/s grid. Kurucz CD-ROM No. 13. Cambridge, Mass.: Smithsonian Astrophysical Observatory, 1993. , 13 .

Labitzke, K. and H. van Loon, 1995: Connection between the troposphere and stratosphere on a decadal scale. Tellus Series A , 47 , 275, doi:10.1034/j.1600-0870.1995. t01-1-00008.

Lean, J., 2000: Evolution of the Sun's spectral irradiance since the Maunder Minimum. Geophysical Research Letters , 27 , 2425-2428, doi:10.1029/2000GL000043.

Lean, J., G. Rottman, J. Harder, and G. Kopp, 2005: SORCE Contributions to New Understanding of Global Change and Solar Variability. Journal of Solar Physics , 230 , 27-53, doi:10.1007/s11207-005-1527-2.

Lean, J. L., G. J. Rottman, H. L. Kyle, T. N. Woods, J. R. Hickey, and L. C. Puga, 1997: Detection and parameterization of variations in solar mid- and near-ultraviolet radiation (200-400 nm). Journal of Geophysical Research , 102 , 29939-29956, doi:10.1029/97JD02092.

Lean, L., Judith and T. DeLand, Matthew, 2012: How Does the Sun's Spectrum Vary? . J. Climate , 25 , 25552560.

Livingston, W. C., J. Harvey, A. K. Pierce, D. Schrage, B. Gillespie, J. Simmons, and C. Slaughter, 1976a: Kitt Peak 60-cm vacuum telescope. Applied Optics , 15 , 3339, doi:10.1364/AO.15.000033.

Livingston, W. C., J. Harvey, C. Slaughter, and D. Trumbo, 1976b: Solar magnetograph employing integrated diode arrays. Applied Optics , 15 , 40-52, doi: 10.1364/AO.15.000040.

Lockwood, M., C. Bell, T. Woollings, R. G. Harrison, L. J. Gray, and J. D. Haigh, 2010: Top-down solar modulation of climate: evidence for centennial-scale change. Environmental Research Letters , 5 (3) , 034008, doi: 10.1088/1748-9326/5/3/034008.

McClintock, W. E., G. J. Rottman, and T. N. Woods, 2005: Solar-Stellar Irradiance Comparison Experiment II (Solstice II): Instrument Concept and Design. Solar Physics , 230 , 225-258, doi:10.1007/s11207-005-7432-x.

Meier, R. R., 1991: Ultraviolet spectroscopy and remote sensing of the upper atmosphere. Space Science Reviews , 58 , 1-185, doi:10.1007/BF01206000.

Merkel, A. W., J. W. Harder, D. R. Marsh, A. K. Smith, J. M. Fontenla, and T. N. Woods, 2011: The impact of solar spectral irradiance variability on middle atmospheric ozone. Geophysical Research Letters , 38 , L13802, doi:10.1029/2011GL047561.

Morrill, J. S., L. Floyd, and D. McMullin, 2011: The Solar Ultraviolet Spectrum Estimated Using the Mg II Index and Ca II K Disk Activity. Solar Physics , 269 , 253-267, doi:10.1007/s11207-011-9708-7.

Pagaran, J., M. Weber, and J. Burrows, 2009: Solar Variability from 240 to 1750 nm in Terms of Faculae Brightening and Sunspot Darkening from SCIAMACHY. The Astrophysical Journal , 700 , 1884-1895, doi:10.1088/0004-637X/700/2/1884.

Rottman, G., 2005: The SORCE Mission. Solar Physics , 230 , 7-25, doi:10.1007/s11207-005-8112-6.

Rottman, G., T. Woods, M. Snow, and G. Detoma, 2001: The solar cycle variation in ultraviolet irradiance. Advances in Space Research , 27 , 1927-1932, doi:10.1016/ S0273-1177(01)00272-1.

Scherrer, P. H., et al., 1995: The Solar Oscillations Investigation - Michelson Doppler Imager. Solar Physics , 162 , 129-188, doi:10.1007/BF00733429.

Schmidt, G. A., et al., 2012: Climate forcing reconstructions for use in PMIP simulations of the Last Millennium (v1.1). Geoscientific Model Development , 5 , 185191, doi:10.5194/gmd-5-185-2012.

Schou, J., et al., 2012: Design and Ground Calibration of the Helioseismic and Magnetic Imager (HMI) Instrument on the Solar Dynamics Observatory (SDO). Solar Physics , 275 , 229-259, doi:10.1007/s11207-011-9842-2.

Snow, M., W. E. McClintock, G. Rottman, and T. N. Woods, 2005: Solar Stellar Irradiance Comparison Experiment II (Solstice II): Examination of the Solar Stellar Comparison Technique. Solar Physics , 230 , 295-324, doi:10.1007/s11207-005-8763-3.

Solanki, S. K. and Y. C. Unruh, 1998: A model of the wavelength dependence of solar irradiance variations. Astronomy &amp; Astrophysics , 329 , 747-753.

Solanki, S. K. and Y. C. Unruh, 2013: Solar irradiance variability. Astronomische Nachrichten , 334 , 145, doi: 10.1002/asna.201211752, 1210.5911 .

Soukharev, B. E. and L. L. Hood, 2006: Solar cycle variation of stratospheric ozone: Multiple regression analysis of long-term satellite data sets and comparisons with models. Journal of Geophysical Research (Atmospheres) , 111 (D10) , D20314, doi:10.1029/2006JD007107.

Steinhilber, F., J. Beer, and C. Fr¨hlich, 2009: o Total solar irradiance during the Holocene. Geophysical Research Letters , 361 , L19704, doi:10.1029/2009GL040142.

Swartz, W. H., R. S. Stolarski, L. D. Oman, E. L. Fleming, and C. H. Jackman, 2012: Middle atmosphere response to different descriptions of the 11-yr solar cycle in spectral irradiance in a chemistry-climate model. Atmospheric Chemistry &amp; Physics , 12 , 5937-5948, doi: 10.5194/acp-12-5937-2012.

Thuillier, G., M. Hers´ e, D. Labs, T. Foujols, W. Peetermans, D. Gillotay, P. C. Simon, and H. Mandel, 2003: The Solar Spectral Irradiance from 200 to 2400 nm as Measured by the SOLSPEC Spectrometer from the Atlas and Eureca Missions. Solar Physics , 214 , 1-22, doi: 10.1023/A:1024048429145.

Unruh, Y. C., W. T. Ball, and N. A. Krivova, 2012: Solar Irradiance Models and Measurements: A Comparison in the 220-240 nm wavelength band. Surveys in Geophysics , 33 , 475-481, doi:10.1007/s10712-011-9166-7, 1111.2068 .

Unruh, Y. C., N. A. Krivova, S. K. Solanki, J. W. Harder, and G. Kopp, 2008: Spectral irradiance variations: comparison between observations and the SATIRE model on solar rotation time scales. Astronomy &amp; Astrophysics , 486 , 311-323, doi:10.1051/0004-6361:20078421, 0802.4178 .

Unruh, Y. C., S. K. Solanki, and M. Fligge, 1999: The spectral dependence of facular contrast and solar irradiance variations. Astronomy &amp; Astrophysics , 345 , 635-642.

Unruh, Y. C., S. K. Solanki, and M. Fligge, 2000: Modelling solar irradiance variations: Comparison with observations, including line-ratio variations. Space Science Reviews , 94 , 145-152.

van Loon, H. and D. J. Shea, 1999: A probable signal of the 11-year solar cycle in the troposphere of the northern hemisphere. Geophysical Research Letters , 26 , 28932896, doi:10.1029/1999GL900596.

Viereck, R. A., et al., 2004: A composite Mg II index spanning from 1978 to 2003. Space Weather , 2 , S10005, doi: 10.1029/2004SW000084.

Wang, S., et al., 2013: Midlatitude atmospheric oh response to the most recent 11-y solar cycle. PNAS , 110 , 1215-1220, doi:10.1073/pnas.1213389110.

Wang, Y.-M., J. L. Lean, and N. R. Sheeley, Jr., 2005: Modeling the Sun's magnetic field and irradiance since 1713. The Astrophysical Journal , 625 , 522-538, doi:10. 1086/429689.

Warwick, N. J., S. Bekki, E. G. Nisbet, and J. A. Pyle, 2004: Impact of a hydrogen economy on the stratosphere and troposphere studied in a 2-D model. Geophysical Research Letters , 31 , L05107, doi:10.1029/2003GL019224.

Wenzler, T., S. K. Solanki, N. A. Krivova, and D. M. Fluri, 2004: Comparison between KPVT/SPM and SoHO/MDI magnetograms with an application to solar irradiance reconstructions. Astronomy &amp; Astrophysics , 427 , 1031-1043, doi:10.1051/0004-6361:20041313.

Wenzler, T., S. K. Solanki, N. A. Krivova, and C. Fr¨hlich, o 2006: Reconstruction of solar irradiance variations in cycles 21-23 based on surface magnetic fields. Astronomy &amp; Astrophysics , 460 , 583-595, doi:10.1051/0004-6361: 20065752.

Woods, T. N., W. K. Tobiska, G. J. Rottman, and J. R. Worden, 2000: Improved solar Lyman α irradiance modeling from 1947 through 1999 based on UARS observations. Journal of Geophysical Research , 105 , 2719527216, doi:10.1029/2000JA000051.

Woods, T. N., et al., 1996: Validation of the UARS solar ultraviolet irradiances: Comparison with the ATLAS 1 and 2 measurements. Journal of Geophysical Research , 101 , 9541-9570, doi:10.1029/96JD00225.

Woods, T. N., et al., 2009: Solar Irradiance Reference Spectra (SIRS) for the 2008 Whole Heliosphere Interval (WHI). Geophysical Research Letters , 36 , L01101, doi:10.1029/2008GL036373.

-

## Supplementary Materials

## a. Cycle variability

The uncertainty of SATIRE-S is time-dependent and increases going backwards in time from 2008. This is due to the way the reconstructions from the different instruments are combined (see Ball et al., 2012, and section 2a of the main paper). We show the estimated cycle uncertainty, as a function of wavelength, for three solar cycles in Fig. S1. The bottom panel of Fig. S1 is similar to the bottom panel of Fig. 1 in the main paper, but is now between March 2000-December 2008 and has the 1 σ uncertainty estimate included as yellow shading. The uncertainty is very small in cycle 23 and is difficult to see. However, cycles 21 (top; December 1979-September 1986) and 22 (middle; July 1989-May 1996) show much larger uncertainty due to steps (ii) and (iii) of the construction method explained in section 2a. We note that even with the larger uncertainty in cycles 21 and 22, it does not generally encompass the cycle changes of NRLSSI for which there is no uncertainty estimate.

## b. Gap filling

We use five solar indices to fill gaps in the SATIRE-S time series. Here we detail the behaviour of these indices and show examples of the gap-filled time series in three spectral bands.

Each index is indicative of the behaviour of some feature (or combination of features) or source of radiative flux in the solar atmosphere, which in turn should correlate with activity and therefore variability in solar spectral irradiance. The variability of each index depends on the time scale being considered. We only use the detrended indices to avoid biasing the long-term SATIRE-S SSI behaviour to any index.

To account for the rotational variability and fill gaps in the detrended SSI time series, we compare each detrended index with each detrended SATIRE-S SSI wavelength. We do this for all wavelengths above 270 nm and for the integrated region between 220-240 nm. Below 270 nm, all time series are determined from the integrated 220-240 nm region so only this one time series needs to be filled to account for all wavelengths in this wavelength region (see section 2b of main paper). In Fig. S2, the coefficient of determination, r 2 c (the square of the correlation coefficient, r c ), shows how well each detrended index can account for the rotational variability at each wavelength. Colours are given in the legend. Also shown is the highest value r 2 c from the multi-linear regressions of all pairings of indices (black). The scale on the x-axis is linear between 115 and 2000 nm and logarithmic above. The constant value of r 2 c below 270 nm reflects the use of the 220-240 nm region, where the dots are located, to determine the wavelengths in this region.

Considering just the individual indices, we briefly describe the source of their variations (except for TSI), and why they agree with the SSI variability in various wavelength intervals. Below 300 nm, the Mg II index agrees best with SATIRE-S SSI wavelengths, accounting for around 70-80% of the rotational variability. The Mg II index is a measure of the ratio of the emission reversal in the line core of the Mg II line (at ∼ 280 nm) and the Mg II wings. The line core is formed higher up in the solar atmosphere than the wings and is associated with regions of strong magnetic flux, as found in faculae and plages. Faculae dominate spectral variability below 300 nm, so the agreement of the Mg II index with these wavelengths is to be expected.

Between 300 and 400 nm, there is a transition between the dominance of faculae and of sunspots on spectral variability. This is reflected in the behaviour of the indices over this spectral interval where correlation coefficients are highly variable and quite poor at some wavelengths. The lower rotational agreement is not so much of a concern for cycle-length changes in these wavelengths, which is more important in the analysis of the long-term impact of SSI on, e.g., stratospheric chemistry. By combining two proxies, which represent sunspot and facular variations, high agreement ( &gt; 80%) with SATIRE-S wavelengths in this region can be maintained, as shown with the black line in Fig. S2.

The regions between 400 and 1300 nm and 2300 and 4000 nm are best represented by TSI, which can represent more than 80% of the variability at almost all of these wavelengths. TSI is replaced by the sunspot area (SSA) and F10.7 cm radio flux between 1400 and 2200 nm as being the best representative single index of rotational variability; this region contains the opacity minimum and represents the deepest directly observable layers of the photosphere where sunspots and pores dominate and where faculae are expected to appear dark. Where the transition to dark faculae occurs is dependent on the facular model atmosphere (Fontenla et al. 1993; Unruh et al. 1999, 2008). As the TSI includes facular brightening, its correlation with the near-IR SSI decreases. We thus find that the SSI is best represented by the SSA and F10.7 cm radio flux. The F10.7 cm radio flux is formed from plasma trapped in coronal loops anchored to sunspots, so, although it is not formed in the photosphere, the F10.7 cm radio flux is highly correlated with the SSA on rotational time scales and hence the modelled irradiance at these wavelengths. Between 4000 and 10 000 nm none of the indices agree well with rotational variability in SATIRE-S SSI. As in the 300-400 nm region, there is a transition between a stronger influence from pores and sunspots below 4500 nm and faculae above; the contrast of faculae increases monotonically from here while penumbral and umbral contrasts remain

Figure S1 - As for the lower panel of Fig. 1 of the main paper, but with cycle amplitude uncertainty included as yellow shading for SATIRE-S. This is given for cycles 21 (top), 22 (middle) and 23 (bottom). Cycle maximum and minimum dates are taken from the monthly sunspot numbers. Cycle 23 in this plot is for the 81-day average fluxes between March 2000 and December 2008, whereas Fig. 1 is between February 2002 and December 2008, which had a higher mean irradiance than around March 2000.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0365v1-figure-52025897-1d52-4f24-9d27-5deec58f206b.jpg)

relatively constant. Once again, by combining two indices, good agreement with the SATIRE-S rotational variability is found.

In Fig. S3, we give three examples of the completed gapfilling process. The integrated region of 220-240 nm and the spectral regions centred at 501 (at 2 nm resolution) and 1605 nm (10 nm) are plotted with the original SATIRE-S time series (black) and the additional daily gap-filled data (blue). The time series is consistent on long- and shorttime scales.

ture of the solar atmosphere and are constrained using spatially resolved observations. Some of the differences between NRLSSI and SATIRE-S at wavelengths above 400 nm are due to the use of model atmospheres with slightly different atmospheric profiles. For example, the extent of the near-IR variability, and even its sign (correlated or anti-correlated with TSI), depends very sensitively on the adopted facular model. We note that the intensities are derived from one-dimensional plane parallel atmospheres and thus do not take the geometry of magnetic features into account.

## c. Uncertainties in model spectrum

The uncertainties we consider in the spectral reconstruction do not include an uncertainty in the model atmospheres themselves, which are estimates of the struc-

Another consideration is the assumption of local thermodynamic equilibrium (LTE) conditions in the solar atmosphere when employing the ATLAS9 code to calculate model intensities. This leads to errors in the modelled ir-

Figure S2 - The coefficient of determination, r 2 c , between the detrended, rotational variability of each wavelength in SATIRE-S and five solar indices: the Mg II index (blue), the PMOD TSI composite (red), the Lymanα composite (purple), the F10.7 cm radio flux (yellow) and the sunspot area (green). Also shown, in black, is the highest value of r 2 c from the multi-linear regression of all pairs of the indices.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0365v1-figure-8be4c700-1f38-43fc-8c5e-82568f3f48df.jpg)

radiance variability in some wavelength regions, mainly in the UV. In the data release, we flag the wavelength regions where the uncertainty due to the use of LTE is large, such as, e.g., Mg I line at 285 nm. In Fig. 8a of (Unruh et al. 2008), the variability of SATIRE-S and SORCE/SIM wavelengths between 220 and 290 nm is shown. Above 270 nm, only the Mg I line in SATIRE-S has much larger rotational variability than SORCE/SIM (more than four times); other lines showing much larger rotational variability in the UV below 270 nm are corrected using the method described in (Krivova et al. 2006) and section 2a of this paper. For wavelengths longer than 285 nm, we flag those that show rotational variability exceeding SORCE/SIM or the instrument noise by 25%. This includes the Ca II H and K lines and the wavelength region around 385 nm (see Fig. 8b of Unruh et al., 2008).

## d. Effect of absolute flux on ozone and temperature

We show here that the choice of absolute flux has very little impact on both the cycle change in ozone and temperature in the stratosphere and lower mesosphere. We do this to support the claim made in section 4 of the main paper that it is the relative change between cycle maximum and minimum that is more important when considering cycle changes in ozone. Figure S4 gives the mean absolute fluxes of SATIRE-S (blue), NRLSSI (red) and SORCE/SOLSTICE versions 12 (green) and 10 (yellow) for the three month period centred on December 2008. The left y-axis scale corresponds to wavelengths below 210 nm and the right scale is for wavelengths longer than this. We see that all four datasets show a difference in absolute flux, but the largest differences are above 290 nm, where relative cycle variability is lowest in the models. SATIRE-S has the largest absolute flux at these wavelengths, while SORCE/SOLSTICE (at 290-300 nm) and NRLSSI (at 300310 nm) show the lowest.

Figure S5 gives the change in temperature, in degrees Kelvin, between 2003 and 2008, corresponding to the same runs that produced Fig. 6 in the main paper. All plots in Fig. S5 show an increasing change in temperature with higher altitude, though NRLSSI and SATIRE-S both show about half the change in temperature compared to SORCE/ SOLSTICE. SORCE/SOLSTICE version 10 shows a slightly higher change in temperature than version 12 at all altitudes.

In Fig. S6, we show that the choice of absolute flux has little impact on the change in temperature and ozone concentration. We demonstrate this for three cases: (i) by putting the ∆F between 2003 and 2008 from NRLSSI on the absolute flux of SATIRE-S with the resulting ozone and temperature response given in Fig. S6a and b, respectively; (ii) by using the ∆F from SORCE/SOLSTICE v12 with the absolute flux of SATIRE-S in Fig. S6c and d; and (iii) the same as (ii), but with the absolute flux of SORCE/SOLSTICE version 10 in Fig. S6e and f. Comparing Fig. 6a of the main paper with Fig. S6a (for ozone concentration) and Fig. S5a with Fig. S6b (for temperature), it can be seen that the change resulting from using either the absolute flux of NRLSSI or SATIRE-S does make a small change to the contours, but these changes are much less than the difference between contour levels. Minor differences are more easily detectable for the examples using SORCE/SOLSTICE v12 changes with the ab-

Figure S3 - Time series for three wavelength regions from SATIRE-S showing both the reconstruction direct from the model (black) and the gaps in the data that have been filled with proxy data (blue). The top plot shows the integrated 220-240 nm flux, the middle plot is for the 2 nm band at 501 nm and the bottom plot is the 10 nm band at 1605

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0365v1-figure-c25c8fac-5927-4369-958e-0b94dfedbcf6.jpg)

solute fluxes of SORCE/SOLSTICE v12 (Figs. 6d of the main paper and S5d), SATIRE-S (Figs. S6c and d) and SORCE/SOLSTICE v10 (Figs. S6e and f). However, once again, these differences are very small and only stand out at the winter pole at 0.3 hPa. It is, therefore, the case that the relative variations in the solar spectral irradiance have a much more significant effect on changes in ozone concentration and temperature in the stratosphere than the choice of absolute flux, at least for the range of absolute fluxes from datasets considered in this study.

Figure S4 - The absolute flux for the 81-day period centred on the solar minimum of December 2008 for SATIRE-S (blue), NRLSSI (red), SORCE/SOLSTICE v12 (yellow) and SORCE/SOLSTICE v10 (green).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0365v1-figure-9a99afd6-d1e2-4539-8d39-8827670f610c.jpg)

Figure S5 - The same as for Fig. 6 of the main paper, but showing the change in stratospheric temperature. Contours are given in 0.1 K levels.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0365v1-figure-296bec46-f270-4f52-9479-6a3c85a0e00e.jpg)

Figure S6 - As for Fig. 6 of the main paper and Fig. S5, contour plots of the change in stratospheric ozone (left panels) and the corresponding temperature change (right panels). The top row shows ozone and temperature changes when employing NRLSSI in the HP model with the absolute flux level of SATIRE-S. The middle row is using SORCE/SOLSTICE v12 change in flux with the absolute flux level of SATIRE-S. The bottom row is using SORCE/SOLSTICE v12 change in flux with the absolute flux of SORCE/SOLSTICE v10.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0365v1-figure-cbee373a-02b7-4cd2-9404-b35cd3a9dada.jpg)
William T. Ball Natalie A. Krivova Yvonne C. Unruh Joanna D. Haigh Sami K. Solanki
A new SATIRE-S spectral solar irradiance reconstruction for solar cycles 21--23 and its implications for stratospheric ozone
We present a revised and extended total and spectral solar irradiance (SSI)
reconstruction, which includes a wavelength-dependent uncertainty estimate,
spanning the last three solar cycles using the SATIRE-S model. The SSI
reconstruction covers wavelengths between 115 and 160,000 nm and all dates
between August 1974 and October 2009. This represents the first full-wavelength
SATIRE-S reconstruction to cover the last three solar cycles without data gaps
and with an uncertainty estimate. SATIRE-S is compared with the NRLSSI model
and SORCE/SOLSTICE ultraviolet (UV) observations. SATIRE-S displays similar
cycle behaviour to NRLSSI for wavelengths below 242 nm and almost twice the
variability between 242 and 310 nm. During the decline of last solar cycle,
between 2003 and 2008, SSI from SORCE/SOLSTICE version 12 and 10 typically
displays more than three times the variability of SATIRE-S between 200 and 300
nm. All three datasets are used to model changes in stratospheric ozone within
a 2D atmospheric model for a decline from high solar activity to solar minimum.
The different flux changes result in different modelled ozone trends. Using
NRLSSI leads to a decline in mesospheric ozone, while SATIRE-S and
SORCE/SOLSTICE result in an increase. Recent publications have highlighted
increases in mesospheric ozone when considering version 10 SORCE/SOLSTICE
irradiances. The recalibrated SORCE/SOLSTICE version 12 irradiances result in a
much smaller mesospheric ozone response than when using version 10 and now
similar in magnitude to SATIRE-S. This shows that current knowledge of
variations in spectral irradiance is not sufficient to warrant robust
conclusions concerning the impact of solar variability on the atmosphere and
climate.