Free Access
Issue
A&A
Volume 531, July 2011
Article Number A6
Number of page(s) 20
Section The Sun
DOI https://doi.org/10.1051/0004-6361/201015843
Published online 30 May 2011

© ESO, 2011

1. Introduction

Radiative energy input from the Sun is the main external source of heat to the Earth’s coupled atmospheric/oceanic system (Hansen 2000; Haigh 2001). Variations in this external source obviously affect the climate, although there is still considerable uncertainty and debate on the level of this effect (Haigh 2001; Hansen et al. 2002, 2005; Solomon 2007; Camp & Tung 2007). One way to disentangle the influence of the Sun on the climate is to look for its influence on climate in pre-industrial times, when the influence of man-made greenhouse gases can be neglected and the natural variations can be isolated.

Direct measurements of atmospheric and oceanic parameters during the last century suggest the occurrence of climate shifts on several time scales (Trenberth 1990; Zhang & Song 2006; Solomon 2007). Holocene climate proxies also suggest the occurrence of millennial scale climate variations (Rodbell et al. 1999; Bond et al. 2001). For example, the intensity and frequency of the El Nino and Southern Oscillation (ENSO), which is the most important oscillation mode of the atmosphere/ocean coupled system, seems to have undergone significant climatic shifts in the Holocene (Rodbell et al. 1999), with strong ENSO events becoming established only about 5000 years ago. Furthermore, Holocene climate proxies also suggest substantial correlations between atmospheric and oceanic conditions and solar variations, but the physical mechanisms are not understood yet (Bond et al. 2001; Haigh 2001, 2007; Usoskin et al. 2005).

With the current generation of general circulation models it is now possible to carry out long-term numerical simulations of the evolution of the Earth’s climate system (Jungclaus et al. 2010; Schmidt 2010). Such studies require long time series of solar irradiance, ideally covering the whole Holocene. Although reconstructions of irradiance since the Maunder Minimum are common (e.g. Lean et al. 1995; Solanki & Fligge 1999; Fligge & Solanki 2000; Wang et al. 2005; Balmaceda 2007; Krivova et al. 2007, 2010a; Tapping et al. 2007; Crouch et al. 2008), on longer timescales few quantitative reconstructions are available due to the lack of sufficiently time-resolved solar data and appropriate techniques for reconstructing the irradiance when only data averaged over 10 years (i.e. roughly a solar cycle) are available (but see Steinhilber et al. 2009). Here we compute the total (and to a limited extend the spectral) irradiance for the Holocene employing the reconstructions of the open flux and sunspot number obtained from the cosmogenic isotope 14C. The motivation for this work came from a study carried out by Balmaceda (2007), who compared the annual and cycle averages of the TSI, total and open magnetic fluxes, as computed using the model described by Krivova et al. (2007), and the observed sunspot number. She found robust relationships between these quantities. Vieira et al. (2008) employed these relationships to make a preliminary computation of the total solar irradiance for the Holocene from the sunspot number and open magnetic flux estimated from 14C. More recently, Steinhilber et al. (2009) presented a reconstruction of the total solar irradiance for the Holocene based on the cosmogenic radionuclide 10 Be measured in ice cores. A shorter reconstruction has been presented by Delaygue & Bard (2010).

In the present work we derive the relationship between solar cycle averaged open magnetic flux and total solar irradiance for the physics-based model of Krivova et al. (2007, 2010a). This is then used to reconstruct the irradiance throughout the Holocene based on 14C data. The reconstructed values lying in the telescope era are compared with the cycle averaged values obtained by Krivova et al. (2010a). The present model is so far unique in that it is based on a consistently derived, physics-based relationship between TSI and the magnetic flux, which incidentally does not correspond to a direct, instantaneous relationship as is often incorrectly assumed. By following this approach the real physical processes are considered in all steps, but some parameters, which cannot be obtained directly from theory, are evaluated by fitting the model’s output to observational data.

We structure the paper as follows. In Sect. 2, we describe the physics-based model used to reconstruct the TSI and SSI, which is, in principle, identical to the models employed by Krivova et al. (2007, 2010a) to calculate the solar irradiance on time scales from days to centuries, but adapted to work with heavily time-averaged data. In Sect. 3, the reconstruction of the TSI and SSI for the Holocene is presented and sources of uncertainty are discussed. Finally, we present the conclusions in Sect. 4.

2. Model description

The primary source of information on solar activity in pre-telescopic times are concentrations of cosmogenic isotopes, such as 14C or 10Be, in natural archives. From these it is possible to deduce the Sun’s open magnetic flux and the sunspot number, as has been demonstrated in numerous publications (e.g., Usoskin et al. 2003, 2006a; Solanki et al. 2004; Steinhilber et al. 2010). In fact, when computing the sunspot number from the open flux, the magnetic flux in active regions, ephemeral regions and the total magnetic flux are computed. These are also the ingredients needed to estimate the total and spectral solar irradiance employing the model described by Krivova et al. (2007, 2010a). The main hurdle is that in general only the cycle-averaged open flux is known, so that the models originally developed for time-resolved data need to be modified in order to be applicable on these long time scales. For the magnetic flux the relevant relationships were most recently derived by Vieira & Solanki (2010). Here, we now take the next step and determine the Sun’s cycle averaged irradiance.

The aim of this section is to derive the relationship between the cycle-averaged TSI and a solar quantity as closely related to the observations (isotope concentrations) as possible. The open magnetic flux of the Sun is such a quantity. We stress that the open magnetic flux itself produces only a rather minor part of solar irradiance variations, but is closely linked to the quantities which produce larger fractions of the irradiance variations. In particular, the cycle averaged open magnetic flux closely follows the Sun’s total magnetic flux.

Before we begin, a brief note on terminology: SATIRE stands for Spectral and Total Irradiance Reconstructions. The basic philosophy of this model family is that all solar irradiance variations on time-scales longer than a day are due to the evolution of the solar surface magnetic field. Members of the SATIRE family of models are SATIRE-S, which models the irradiance during the satellite era with high precision (in its present version since 1974), and SATIRE-T, which computes irradiance during the telescope era (since 1610). Here we introduce the SATIRE-M model, where “M” stands for millennial time scales.

2.1. Model of solar irradiance variations

The Sun’s irradiance responds directly to the amount and distribution of magnetic flux on the solar disk, which gives rise to dark features, sunspots, and bright features, faculae and network. In order to quantitatively reproduce solar radiative flux it is sufficient to segment the solar disk, using magnetograms and brightness images, into dark (sunspot umbra and penumbra) and bright (faculae and network) magnetic components as well as the unaffected quiet Sun (Krivova et al. 2003; Wenzler et al. 2005, 2006). The computed irradiance depends both on the area covered by a given type of feature and the spatial location on the solar disk. In the period before magnetograms became available, a simplified version of the model had to be employed. In this model a homogeneous distribution of magnetic features over the solar surface is assumed, so that only the fractional area of the solar disk covered by each of them (global, or disk-averaged filling factor, α) and the specific flux, F, (i.e. the flux that would be emitted if the whole solar disk was covered by that type of feature) of each feature is required. The evolution of radiative flux at a given wavelengths, λ, can then be computed as follows (Krivova et al. 2007, 2010a), F(λ,t)=αq(t)Fq(λ)+αu(t)Fu(λ)+αp(t)Fp(λ)\begin{eqnarray} F(\lambda, t) & = & \alpha_{\rm q}(t) F_{\rm q}(\lambda) + \alpha_{\rm u}(t) F_{\rm u}(\lambda) + \alpha_{\rm p}(t) F_{\rm p}(\lambda) \nonumber \\ & & + \, \alpha_{\rm f}(t) F_{\rm f}(\lambda) + \alpha_{\rm n}(t) F_{\rm n}(\lambda), \label{IMeq01} \end{eqnarray}(1)where the Fq, Fu, Fp, Ff, and Fn are the time-independent radiative fluxes of the five component model: quiet Sun (q), sunspot umbra (u), sunspot penumbra (p), faculae (f), and network (n). Following Krivova et al. (2007) and others, the spectrum of each component that was calculated by Unruh et al. (1999) using the ATLAS9 code of Kurucz (1993). Here, we employ the same model to describe the brightness of facular and network components.

The fraction of the solar disk covered by each component is represented by αq, αu, αp, αf and αn, respectively. By taking this approach we are assuming a nearly homogeneous distribution of magnetic features on the solar disk. The filling factor of the quiet Sun can then be written as: αq(t)=1αu(t)αp(t)αf(t)αn(t).\begin{equation} \alpha_{\rm q}(t) = 1 - \alpha_{\rm u}(t) - \alpha_{\rm p}(t) - \alpha_{\rm f}(t) - \alpha_{\rm n}(t). \label{IMeq02} \end{equation}(2)Substituting Eqs. (2) in (1), we obtain F(λ,t)=(1αu(t)αp(t)αf(t)αn(t))Fq(λ)+αu(t)Fu(λ)\begin{eqnarray} F(\lambda, t) &= & (1 - \alpha_{\rm u}(t) - \alpha_{\rm p}(t) - \alpha_{\rm f}(t) - \alpha_{\rm n}(t)) F_{\rm q}(\lambda) + \alpha_{\rm u}(t) F_{\rm u}(\lambda) \nonumber \\ & & + \, \alpha_{\rm p}(t) F_{\rm p}(\lambda) + \alpha_{\rm f} F_{\rm f}(\lambda) + \alpha_{\rm n}(t) F_{\rm n}(\lambda). \label{IMeq03} \end{eqnarray}(3)This equation can be rewritten as F(λ,t)=Fq(λ)+αu(t)ΔFu(λ)+αp(t)ΔFp(λ)\begin{eqnarray} F(\lambda, t) &= & F_{\rm q}(\lambda) + \alpha_{\rm u}(t) \Delta F_{\rm u}(\lambda) + \alpha_{\rm p}(t) \Delta F_{\rm p}(\lambda) \nonumber \\ & &+ \, \alpha_{\rm f}(t) \Delta F_{\rm f}(\lambda) + \alpha_{\rm n}(t) \Delta F_{\rm n}(\lambda), \label{IMeq04} \end{eqnarray}(4)where ΔFu,p,f,n(λ)=Fu,p,f,n(λ)Fq(λ).\begin{equation} \Delta F_{\rm u,p,f,n}(\lambda) = F_{\rm u,p,f,n}(\lambda) - F_{\rm q}(\lambda). \label{IMeq05} \end{equation}(5)

2.2. Filling factors for sunspot umbrae and penumbrae

We compute the filling factors for the umbral and penumbral components by employing the approach described by Krivova et al. (2007). A fixed ratio of 0.2 between umbral and total spot area is assumed following Wenzler et al. (2005), which translates into αu/αs=αu/(αu+αp)=0.2.\begin{equation} \alpha_{\rm u} / \alpha_{\rm s} = \alpha_{\rm u} / (\alpha_{\rm u} + \alpha_{\rm p}) = 0.2. \label{IMeq08} \end{equation}(6)We compute the sunspot area, i.e. the fraction of the disk covered by all sunspots on the solar disk, by making use of a linear relationship to the sunspot number (R) (Fligge & Solanki 1997; Balmaceda et al. 2009; Hathaway 2010): αs=A1R+A2,\begin{equation} \alpha_{\rm s} = A_1 R + A_2, \label{IMeq09} \end{equation}(7)where A1 is the proportionality coefficient and A2 is the offset. The key step for extending the irradiance model to the whole Holocene is to compute the sunspot number as well as the magnetic flux in active regions (AR) and ephemeral regions (ER) from the open magnetic flux, since it is the first solar parameter obtained from the chain of models that convert cosmogenic isotopes data into solar parameters. This problem was first addressed by Usoskin et al. (2002) and later updated by Usoskin et al. (2003, 2007), Solanki et al. (2004), and Vieira & Solanki (2010). In the current work we employ the method described by Vieira & Solanki (2010) to compute the total solar surface magnetic flux and the sunspot number from the open flux. From the information recovered from cosmogenic isotopes, in particular the INTCAL98 14C record (Stuiver et al. 1998), we are limited to estimate 10-year averages of the sunspot area and the solar surface magnetic components. Following Vieira & Solanki (2010) we write for decadal values Rj=aRφopenj+bRφopenj+1,\begin{equation} \Rj =a_R \openj + b_R \openjj, \label{IMeq10} \end{equation}(8)where  ⟨ ... ⟩ j and  ⟨ ... ⟩ j + 1 are the jth and jth + 1 10-year averaged values (or data points) of a variable. Since 10-years is close to the average length of the solar cycle (approximately 11 years for the past 23 cycles), jth and jth + 1 also roughly refer to the jth and jth + 1 cycles. The coefficients aR and bR are given by aR=1cτ1\begin{equation} a_R = \frac{1}{c \tau_1} \label{IMeq11} \end{equation}(9)and bR=1cΔt,\begin{equation} b_R = \frac{1}{c \Delta t}, \label{IMeq12} \end{equation}(10)where Δt is the sampling interval, τ1 is given by 1τ1=1τopens1Δt,\begin{equation} \frac{1}{\tau_1} = \frac{1}{\tau_{\rm open}^{\rm s}} - \frac{1}{\Delta t}, \label{IMeq13} \end{equation}(11)and the constant c is given by c=[(1τacts+τopenrτopensτactr)τact+τephkXτephs]ϵactmax,21Rmax,21·\begin{equation} c = \left[ \left( \frac{1}{\tau_{\rm act}^{\rm s}} + \frac{\tau_{\rm open}^{\rm r}}{\tau_{\rm open}^{\rm s} \tau_{\rm act}^{\rm r}}\right) \tau_{\rm act} + \frac{\tau_{\rm eph}kX}{\tau_{\rm eph}^{\rm s}} \right] \frac{\epsilon_{\rm act}^{\rm max,21}}{R^{\rm max,21}}\cdot \label{IMeq14} \end{equation}(12)In Eqs. (8) − (12), τopenr\hbox{$\tau_{\rm open}^{\rm r}$} and τopens\hbox{$\tau_{\rm open}^{\rm s}$} are the decay times for the rapidly and slowly evolving components of the open flux. τacts\hbox{$\tau_{\rm act}^{\rm s}$} and τephs\hbox{$\tau_{\rm eph}^{\rm s}$} are the flux transfer times from active and ephemeral regions to the slowly evolving open magnetic flux, while τactr\hbox{$\tau_{\rm act}^{\rm r}$} is the flux transfer times from active regions to the rapidly evolving open magnetic flux. ϵactmax,21\hbox{$\epsilon_{\rm act}^{\rm max,21}$} and Rmax,21 are the maximum emergence rate of magnetic flux in active regions and the maximum value of the sunspot number observed in cycle 21, respectively. The parameter k is the ratio between sunspot number amplitude and the 10-year averaged value (Usoskin et al. 2007). X is the scaling factor for the flux emergence rate in ER relative to active regions. Here, we employ the set of parameter values employed by Krivova et al. (2010a) (see Table 1). The sampling interval, Δt, is 10 years. The effective decay time scales of active regions (AR), τact, and ephemeral regions (ER), τeph, are given by 1τact=1τact0+1τacts+1τactr\begin{equation} \frac{1}{\tau_{\rm act}} = \frac{1}{\tau_{\rm act}^0} + \frac{1}{\tau_{\rm act}^{\rm s}} + \frac{1}{\tau_{\rm act}^{\rm r}} \label{eq13} \end{equation}(13)and 1τeph=1τeph0+1τephs,\begin{equation} \frac{1}{\tau_{\rm eph}} = \frac{1}{\tau_{\rm eph}^0} + \frac{1}{\tau_{\rm eph}^{\rm s}}, \label{eq14} \end{equation}(14)where, τact0\hbox{$\tau_{\rm act}^0$} and τeph0\hbox{$\tau_{\rm eph}^0$} are the decay time scales of active and ephemeral regions, respectively.

Combining Eqs. (6) − (8), we obtain αuj=0.2(A1aRφopenj+A1bRφopenj+1+A2)=\begin{eqnarray} \alphauj & = &0.2 (A_1 a_R \openj + A_1 b_R \openjj + A_2) \nonumber \\ & = &a_{\rm u} \openj + b_{\rm u} \openjj + c_{\rm u} \label{IMeq15} \end{eqnarray}(15)and αpj=0.8(A1aRφopenj+A1bRφopenj+1+A2)=\begin{eqnarray} \alphapj &=& 0.8 (A_1 a_R \openj + A_1 b_R \openjj + A_2) \nonumber \\ &=& a_{\rm p} \openj + b_{\rm p} \openjj + c_{\rm p}. \label{IMeq16} \end{eqnarray}(16)The sunspots of a given cycle influence not only the open flux of that particular cycle, but also of the next cycle due to the slow decay of the open flux (Solanki et al. 2000). Consequently, it takes the knowledge of the open flux from 2 cycles to get the sunspot areas and the magnetic fluxes of active and ephemeral regions for one cycle, as indicated by Eqs. (15) and (16). It means that the part of the magnetic flux that emerges in active and ephemeral regions, and is dragged outward by the solar plasma and permeates the heliosphere, is implicitly represented in Eq. (8). The relative contribution between the jth and jth + 1 data points (aR/bR=(Δtτopens)/τopens\hbox{${a_R}/{b_R} = {(\Delta t - \tau_{\rm open}^{\rm s})}/{\tau_{\rm open}^{\rm s}}$}) depends on the decay time scale of the slow component of the open flux and on the sampling interval. For the set of parameters employed here (given in Table 1), the contribution of the jth data point relative to that of the jth + 1 data point is aRbR=aubu=apbp2.4,\begin{equation} \frac{a_R}{b_R} = \frac{a_{\rm u}}{b_{\rm u}} =\frac{a_{\rm p}}{b_{\rm p}} \approx 2.4, \end{equation}(17)which suggests that a direct, instantaneous relation between sunspot areas and the open flux is not appropriate.

Table 1

Parameters of the TSI model presented by Krivova et al. (2010a) and employed here.

2.3. Filling factors for faculae and network

We compute the 10-year averaged filling factors for faculae (⟨ αf ⟩ j) and the network (⟨ αn ⟩ j) from the solar surface magnetic flux employing the same scheme as in previous applications of the SATIRE models (Krivova et al. 2003, 2007, 2010a; Wenzler et al. 2004, 2005, 2006).  ⟨ αf ⟩ j and  ⟨ αn ⟩ j are proportional to the magnetic flux in faculae (⟨ φf ⟩ j) and network (⟨ φn ⟩ j), respectively, until a saturation limit φsat is reached. The saturation flux is a free parameter of the SATIRE models, which depends on the resolution and the noise level of the employed magnetic field maps. Observing that the 10-year averaged values  ⟨ φf ⟩ j and  ⟨ φn ⟩ j are always smaller than the saturation level, so that they never saturate, we can write αfj=φfjSsunBsat,f,\begin{equation} \alphafj = \frac{\facj}{S_{\rm sun} B_{\rm sat,f}}, \label{IMeq17} \end{equation}(18)and αnj=φnjSsunBsat,n,\begin{equation} \alphanj = \frac{\netj}{S_{\rm sun} B_{\rm sat,n}}, \label{IMeq18} \end{equation}(19)where Ssun is the Sun’s surface area. Bsat,f and Bsat,n are the saturation field strengths for faculae and network, respectively. Here we follow Krivova et al. (2007) and Krivova et al. (2010a) who incorporated a correction in Bsat,f that takes into account that the faculae mainly appear in the active belts.

The magnetic flux in active regions (⟨ φact ⟩ j) is the sum of the flux in sunspots (umbrae plus penumbrae) and faculae. Thus, the magnetic flux in faculae is given by φfj=φactjφujφpj.\begin{equation} \facj = \actj - \uj - \pj. \label{IMeq19} \end{equation}(20)Following Krivova et al. (2007), we compute φuj+φpj=αuj[Bz,u]Ssun+αpj[Bz,p]Ssun,\begin{equation} \uj + \pj= \alphauj [B_{z,{\rm u}}] S_{\rm sun}+ \alphapj [B_{z,{\rm p}}] S_{\rm sun}, \label{IMeq20} \end{equation}(21)where  [Bz,u]  and  [Bz,p]  are the area-averaged vertical components of the magnetic field strength in umbrae and penumbrae, respectively. Note that since  ⟨ αu ⟩ j and  ⟨ αp ⟩ j are indirectly obtained from  ⟨ R ⟩ j, it is assumed implicitely that averaged over 10 years there is no major difference between the number of sunspots on the visible part of the solar surface (solar disk) and on the hidden part.

From Eq. (25) of Vieira & Solanki (2010), the magnetic flux in active regions is approximately given by φactjτactϵactmax,21Rmax,21Rj.\begin{equation} \actj \approx \tau_{\rm act} \frac{\epsilon_{\rm act}^{\rm max,21}}{R^{\rm max,21}} \langle R \rangle _j. \label{IMeq21} \end{equation}(22)Combining Eqs. (15), (18) − (22), we obtain αfj=afφopenj+bfφopenj+1+cf.\begin{equation} \alphafj =a_{\rm f} \openj + b_{\rm f} \openjj + c_{\rm f}. \label{IMeq23} \end{equation}(23)The sum of ephemeral regions flux and of the open flux describes the evolution of the network: φnj=φephj+φopenj.\begin{equation} \netj = \ephj + \openj. \label{IMeq24} \end{equation}(24)From Eq. (29) of Vieira & Solanki (2010), the 10-year averaged magnetic flux in ephemeral regions is φephjτephϵactmax,21Rmax,21kXRj\begin{eqnarray} \ephj &\approx& \tau_{\rm eph} \frac{\epsilon_{\rm act}^ {\rm max,21}}{R^{\rm max,21}} k X \Rj \nonumber \\ &\approx& \tau_{\rm eph} \frac{\epsilon_{\rm act}^ {\rm max,21}}{R^{\rm max,21}} k X \left( a_R \openj + b_R \openjj \right). \label{IMeq25} \end{eqnarray}(25)The contribution from the next i.e. (j + 1)-st cycle is due to the extended period over which ephemeral regions emerge. From Eqs. (24) and (25), φnj=aφnφopenj+bφnφopenj+1,\begin{equation} \netj = a_{\phi_{\rm n}} \openj + b_{\phi_{\rm n}} \openjj, \label{IMeq26} \end{equation}(26)where aφn=τephϵactmax,21Rmax,21kXaR+1\begin{equation} a_{\phi_{\rm n}} = \tau_{\rm eph} \frac{\epsilon_{\rm act}^{\rm max,21}}{R^{\rm max,21}} k X a_R + 1 \label{IMeq27} \end{equation}(27)and bφn=τephϵactmax,21Rmax,21kXbR.\begin{equation} b_{\phi_{\rm n}} = \tau_{\rm eph} \frac{\epsilon_{\rm act}^{\rm max,21}}{R^{\rm max,21}} k X b_R. \label{IMeq28} \end{equation}(28)Substituting Eqs. (26) into (19), we obtain αnj=aφnφopenj+bφnφopenj+1SsunBsat,n=\begin{eqnarray} \alphanj &=& \frac{a_{\phi_{\rm n}} \openj + b_{\phi_{\rm n}} \openjj}{S_{\rm sun} B_{\rm sat,n}} \nonumber \\ &=& a_{\rm n} \openj + b_{\rm n} \openjj. \label{IMeq29} \end{eqnarray}(29)While the relative contributions between jth and jth + 1 data points of the open flux for the umbral, penumbral, and facular filling factors depend on the decay time of slowly evolving open flux and the sampling time, the relative contribution to the network is given by anbn=Δtτopens+Rmax,21cΔtϵactmax,21τephkX1.\begin{equation} \frac{a_{\rm n}}{b_{\rm n}}= \frac{\Delta t}{\tau_{\rm open}^{\rm s}} + \frac{R^{\rm max,21} c \Delta t } { \epsilon_{\rm act}^{\rm max,21} \tau_{\rm eph} k X} - 1. \label{IMeq30} \end{equation}(30)For the set of parameters employed here, this ratio is approximately 3.3. This indicates that although the contribution from the next cycle is lower than that from the cycle for which we want to know  ⟨ αn ⟩ j, it is still necessary to know the next state in order to evaluate the filling factor of the network accurately.

2.4. Estimate of the total solar irradiance

The solar irradiance at wavelength λ can be computed by substituting the expressions for the filling factors obtained for the components into Eq. (4) F(λ)j=aF(λ)φopenj+bF(λ)φopenj+1+Fq(λ),\begin{equation} \Fj = a_{F} (\lambda) \openj+ b_{F}(\lambda) \openjj + F_{\rm q}(\lambda), \label{IMeq31} \end{equation}(31)where aF(λ)=auΔFu(λ)+apΔFp(λ)+afΔFf(λ)+anΔFn(λ),\begin{equation} a_{F}(\lambda) = a_{\rm u} \Delta F_{\rm u}(\lambda) + a_{\rm p} \Delta F_{\rm p}(\lambda) + a_{\rm f} \Delta F_{\rm f}(\lambda) + a_{\rm n} \Delta F_{\rm n}(\lambda), \label{IMeq32} \end{equation}(32)and bF(λ)=buΔFu(λ)+bpΔFp(λ)+bfΔFf(λ)+bnΔFn(λ).\begin{equation} b_{F}(\lambda) = b_{\rm u} \Delta F_{\rm u}(\lambda) + b_{\rm p} \Delta F_{\rm p}(\lambda) + b_{\rm f} \Delta F_{\rm f}(\lambda) + b_{\rm n} \Delta F_{\rm n}(\lambda). \label{IMeq33} \end{equation}(33)In order to obtain the total solar irradiance one can either integrate  ⟨ F(λ) ⟩ j over all wavelengths or right from the beginning work with wavelength integrated radiative fluxes and flux differences Fqt=0Fq(λ)dλ,\begin{equation} F_{\rm q}^t = \int_{0}^{\infty} F_{\rm q}(\lambda) {\rm d}\lambda, \label{IMeq34} \end{equation}(34)and ΔFu,p,f,nt=0[Fu,p,f,n(λ)Fq(λ)]dλ.\begin{equation} \Delta F_{\rm u,p,f,n}^t = \int_{0}^{\infty} \left[ F_{\rm u,p,f,n}(\lambda) - F_{\rm q}(\lambda) \right] {\rm d}\lambda. \label{IMeq35} \end{equation}(35)Replacing the Fq and ΔFi in Eqs. (31) − (33) by Fqt\hbox{$F_{\rm q}^t$} and ΔFit\hbox{$\Delta F_{i}^t$}, i = u,p,f,n, one then obtains the total solar irradiance averaged over 10 years,  ⟨ Ft ⟩ j. We refer for clarity to  ⟨ Ft ⟩ j and Fitj\hbox{$\langle F_i^t \rangle _j$}, i = u,p,f,n, as TSI and TSIi.

Hence, the solar spectral and total irradiance can be computed from the jth and jth + 1 values of the open flux in a relatively straightforward manner.

2.5. Reconstruction of the open flux from cosmogenic data

The open flux in the past can be reconstructed from data of cosmogenic isotopes measured in terrestrial archives. Here we use data from 14C measured in tree-rings for the past millennia as compiled by the INTCAL04 collaboration (Reimer et al. 2004). We have applied a commonly used approach of inverting the cosmogenic isotope data into cosmic ray flux. This approach, proposed Usoskin et al. (2003, 2004), is now a standard method used to reconstruct of the solar or heliospheric parameters in the past (e.g. Solanki et al. 2004; Usoskin et al. 2006a, 2007; Vonmoos et al. 2006; Muscheler et al. 2007; Steinhilber et al. 2008). The method includes the following steps.

  • 1.

    First, the measured values of Δ14C are converted, using a carbon cycle model, into the global 14C production rate, Q14C. We used a standard multi-box carbon cycle model (Siegenthaler 1980) parameterized using the approach of Usoskin & Kromer (2005).

  • 2.

    Considering prescribed geomagnetic field data and a model of 14C production in the atmosphere (Castagnoli & Lal 1980), the production rate Q14C was converted into the cosmic ray spectrum quantified via the modulation potential Φ.

  • 3.

    Using a standard theory of cosmic ray transport in the Heliosphere one can relate the modulation potential Φ to the open magnetic flux assuming a constant 11-yr averaged solar wind speed (e.g. Usoskin et al. 2002).

Here we give a brief description of the method, while full details are provided elsewhere (Solanki et al. 2004; Usoskin et al. 2007).

The uncertainties of the reconstructed open flux σφopen have several sources. The major source of the systematic error is related to the uncertainties of the 14C production rate in the pre-Holocene period. Other sources of systematic uncertainties are related to the paleomagnetic models (considered below in detail), uncertainties of the carbon cycle model, uncertainties of the 14C production model, and uncertainties of the cosmic ray heliospheric transport model. Random errors are related to measurement errors in Δ14C and in geomagnetic data. All these uncertainties are described in detail by Solanki et al. (2004, Supplementary matarial). The total uncertainty in the open flux σφopen is of the order of 0.2  ×  1014 Wb for the last five millennia with relatively good geomagnetic data and rises up to 0.5  ×  1014 Wb in the early Holocene.

We note that computation of the open flux from 14C is not a part of the present work as it was performed using the same algorithms and codes as in Usoskin et al. (2007), but now applying different paleomagnetic models. We point out that the 14C-based open flux is consistent with the Ti data from meteorites (Usoskin et al. 2006b).

Table 2

Coefficients of the model to compute TSI for the Holocene (SATIRE-M).

thumbnail Fig. 1

a) Evolution of the solar open magnetic flux. The black and red lines present the estimates of the open flux from the geomagentic aa-index and from observations of the interplanetary magnetic field (Lockwood et al. 2009a,b), respectively. The blue line displays the modeled open flux from SATIRE-T. b) 44Ti activity in stony meteorites as a function of the time of fall (see Taricco et al. 2006, and references therein). Error bars refer to 1-sigma uncertainties. The red line corresponds to the modeled 44Ti based on the SATIRE-T open flux according to Usoskin et al. (2006b). The 44Ti activity is given in disintegrations per minute per kilogram (dpm kg-1) of Ni  +  Fe. The dashed red lines correspond to the estimated 1-sigma uncertanty of the model.

3. Results

3.1. Model setup and validation

The equations modeling the evolution of TSI during the Holocene have been derived in the previous section. In this section, we present some tests. We compute the coefficients for the solar spectral and total irradiance model for the Holocene according to the parameters of the SATIRE-T model presented by Krivova et al. (2010a). Table 2 lists the coefficients entering the equations given in Sect. 2 that are used to compute the solar irradiance from the open magnetic flux (which is deduced from 14C data). We stress that we have not introduced new free parameters when simplifying the model to work with 10-year averaged values. Thus, uncertainties in the parameters are a heritage from the SATIRE-T model described by Krivova et al. (2010a).

Krivova et al. (2010a) constrained the parameters of the model by comparing the model output with observations and reconstructions of five solar variables. While two data sets are magnetic flux variables, three data sets are related directly to the solar irradiance. The total solar surface magnetic flux was constrained by comparing the model output with estimates carried out at the Mt. Wilson Solar Observatory (MWO), National Solar Observatory Kitt Peak (KP NSO) and Wilcox Solar Observatory (WSO) over cycles 20 to 23 (Arge et al. 2002; Wang et al. 2005). The open flux was compared to the reconstruction based on the geomagnetic aa-index (Lockwood et al. 2009a,b) from 1904 to 2008. These two data sets provided constraints on the amplitude of the magnetic activity cycle and its long-term evolution. We stress that the evolution of the open flux as reconstructed based on the aa-index, which nearly doubled during the last century, is well reproduced by the magnetic field model (see Fig. 1a). In addition to the reconstruction based on the geomagnetic aa-index, the long-term evolution of solar magnetic activity can also be constrained by cosmogenic isotopes deposited in terrestrial and extra-terrestrial archives (Solanki et al. 2004; Usoskin et al. 2006b). While terrestrial archives can span the entire Holocene, the processes leading to the deposition of the cosmogenic isotopes are affected by the evolution of, e.g., the Earth’s magnetic field and of the highly coupled atmosphere-ocean system. Extraterrestrial archives, such as meteorites, are not affected by such processes and can be employed to test the long-term trend observed in the modeled open flux. In order to substantiate the analysis carried out by Krivova et al. (2010a) and Vieira & Solanki (2010), we test the reconstruction of the open flux by comparing the activity of cosmogenic isotope Ti in meteorites, which fell during the past 240 years, and the modeled activity based on the open flux computed in the framework of the above models. We carried out the analysis using the same method as described by Usoskin et al. (2006b) including evaluation of all known sources of uncertainties. Figure 1b presents a comparison between the Ti measurements in 19 meteorites, indicated by diamonds, and the modeled Ti activity at the time of the meteorite’s fall (see details in Usoskin et al. 2006b; Taricco et al. 2006). The Ti activity is given in disintegrations per minute per kilogram (dpm kg-1). The χ2 is about 7.6 for the 19 data points, i.e. about 0.4 per degree of freedom, that is comparable to the results previously discussed by Usoskin et al. (2006b). Thus, the open flux modeled employing the parameters of the SATIRE-T model follows the data reasonably well, which leads us to conclude that the parameterization used is consistent with the available longer term solar data.

thumbnail Fig. 2

a) Comparison between the observed (PMOD Composite; blue line) and the modeled (SATIRE-T; green line) TSI from Nov./1978 to Dec./2003. b) Difference between the observed and modeled TSI. The gray and black lines present the daily and 90-day averaged values, respectively. c) Distribution of the difference between daily values of the observed and modeled TSI (model error). The standard deviation and the median absolute deviation are indicated in the panel. The red curve represents the corresponding probability distribution function for the normal distribution with the mean and standard deviation of the model error. d) Cumulative distribution function (CDF) for the difference between observed and modeled TSI. e) Probability plot for normal distribution with the mean and standard deviation of the model error. f) Scatter plot of the observed versus modeled TSI. The red line indicates the set of expectation values of the model. The correlation coefficient and the coefficients of the best linear fitting (green line) are presented in the panel. The p-value for testing the hypothesis of no correlation is lower than 0.5, which indicates that the correlation is significant.

The model discussed by Krivova et al. (2010a) was further constrained by requiring the computed TSI variations to match the PMOD composite (Fröhlich 2006), which is obtained from space-based observations since 1978. Here, we employ (among other measurements) the difference between the SATIRE-T model output and the PMOD composite to estimate the uncertainty of the TSI model for the Holocene. We refer to this difference as the model error. Figure 2a presents a comparison between the observed (PMOD composite; blue line) and the modeled (SATIRE-T; green line) TSI. The plot covers the period from Nov./1978 to Dec./2003, which was employed by Krivova et al. (2010a) to constrain the SATIRE-T model. Note that as the absolute value of the TSI is not well determined we offset the model to match the PMOD composite level, although observations employing the SORCE/TIM instrument suggest that the absolute value during the minimum between cycles 23 − 24 is approximately 4 − 5 W/m2 lower than that measured by other current instruments (Kopp et al. 2005; Kopp & Lean 2011). Figure 2b presents the difference between the PMOD composite and the adjusted SATIRE-T model. The gray and black lines display daily values and 90-day means, respectively. Figure 2c presents the distribution of the model error. The standard deviation (σ) and the median absolute deviation (MAD) are indicated in the panel. The MAD value is defined as median(abs(X − median(X))), where X is the sample population. The red line represents the corresponding normal distribution for the estimated mean and standard deviation. While the empirical cumulative distribution function (CDF) of the model error is presented in Fig. 2d, the probability plot for a Normal Distribution is presented in Fig. 2e. In the probability plot, the marks (“x”) indicate the empirical probability versus the model error values. Note that the scale of the ordinate axis is not linear. The distance between the ordinate tick marks is proportional to the distance between the quantiles of a normal distribution. The dashed line connects the 25th and the 75th quantiles. The assumption of normality is reasonable if the data points are near the dashed line. However, this plot clearly indicates that the distribution is not normal. This is a consequence of the procedure employed to constrain the model, which reduced the weight of outliers to compute the χ2. Here, we use the MAD, which is a more robust estimator than the standard deviation and variance, as a measure of the dispersion of the model error distribution. In order to estimate the dispersion from the MAD value we multiply the MAD value by the scale parameter for the normal distribution (Mosteller & Tukey 1977). The value obtained is approximately 0.28 W/m2. Note that this uncertainty applies to daily data, whereas we are dealing with decadal averages. For these the difference between reconstruction and observations is very much smaller. Nonetheless, we retain 0.28 W/m2 as the model error, since it corresponds roughly to the uncertainty in the Maunder Minimum level (Krivova et al. 2007). Figure 2f presents the scatter plot of the observed (PMOD Composite) versus modeled TSI (SATIRE-T). The green line presents the best linear fit to the data set. The coefficients and their uncertainties are indicated in the figure as well as the correlation coefficient (Rc = 0.79).

We validate the SATIRE-M model during the telescope era by comparing its output with that of the SATIRE-T model. The two models are not independent, so that with these comparisons we mainly test the simplifications introduced in order to be able to adequately deal with the 10-year averaged data. To compare both models, first we compute 10-year averages from 1700 AD to 2000 AD of the open flux, total magnetic flux and TSI as well as of the contributions from faculae, network and sunspots (umbrae + penumbrae) to TSI from the SATIRE-T model. We then employ the 10-year averaged open flux from SATIRE-T as input to SATIRE-M. This allows us to compute the contribution of each atmospheric component to TSI employing SATIRE-M.

Figure 3 presents a comparison between the output of the models SATIRE-T and SATIRE-M. Each panel contains the following elements: the red lines indicate the set of expectation values for the model, while the blue lines are the best fit lines. The linear model coefficients and the correlation coefficient are also shown in each panel. The vertical error bars correspond to the dispersion of the SATIRE-T model error computed from the MAD value of the SATIRE-T error distribution and corrected for the normal distribution. The estimate of the SATIRE-M error is discussed in the next subsection. In addition to the dispersion of the TSI error distribution, the other constraints imposed by Krivova et al. (2010a) can also be employed to compute the dispersion of the model error distributions of the total magnetic flux and components of the TSI. Additionally to the constraint imposed by the PMOD composite, Krivova et al. (2010a) constrained the amplitude of the facular contribution to TSI by comparing the model output with the estimate computed by Wenzler et al. (2006). The dispersion of the difference between the models is about 0.27 W/m2 when corrected for the normal distribution. Although the sunspot deficit computed by Wenzler et al. (2006) was not employed to constrain the SATIRE-T model, it could be employed to estimate the uncertainty of the sunspot (umbral + penumbral) contribution to TSI. In this case, the dispersion is about 0.14 W/m2. In order to ensure that not only the total irradiance is reproduced correctly, we also compare the solar irradiance integrated over the wavelength range 220 − 240 nm (see Krivova et al. 2009, 2010b) with the SATIRE-T model output by Krivova et al. (2010a).

The SATIRE-M model reproduces quite accurately the evolution of the total magnetic flux as computed by the SATIRE-T model (Rc = 0.98; slope = 1.00  ±  0.02; Fig. 3), but it slightly overestimates the magnitude of the total magnetic flux (offset =  −1.88  ±  0.79 (× 1014 Wb)). The SATIRE-M model also reproduces well the total solar irradiance computed employing SATIRE-T (Rc = 0.98), but slightly overestimates low values and underestimates high values, which result in a slope = 1.09  ±  0.05). A similar behavior is also observed when comparing the facular and network brightness excess. The coefficient of correlation is approximately 0.98 while the slope is approximately 1.05. The sunspot (umbrae + penumbrae) brightness deficit is also well reproduced but the coefficient of correlation (Rc = 0.96) is lower than the one obtained for the comparison of the brightness excess from faculae and network. The slope of the linear regression is close to unity (slope = 0.98). The UV irradiance integrated from 220 to 240 nm follows the behavior of the facular and network contribution of the TSI, with a high coefficient of correlation (Rc = 0.99).

The correlation coefficients for all cases exceed 0.96 and the differences between the models display nearly Gaussian distributions. Thus, we conclude that the 10-year averages of the TSI (and UV irradiance) can be obtained from a linear combination of the jth and jth + 1 observations of the open magnetic flux with reasonable accuracy.

thumbnail Fig. 3

Comparison between 10-year averaged values of the parameters computed employing the SATIRE-M and SATIRE-T models from 1700 AD to 2000 AD. The following scatter plots are presented: a) total magnetic flux; b) total solar irradiance; c) contribution of the bright components (faculae + network) to the total solar irradiance; d) sunspot (umbra + penumbra) contribution to TSI; e) UV integrated in the range 220−240 nm. The red lines present the set of expectation values for the models and the blue lines are the best linear fits. The correlation coefficient and the coefficients of the linear regression are presented in the panels.

thumbnail Fig. 4

Reconstructions of the geomagnetic field intensity. The VADM reconstructions by Knudsen et al. (2008) and Genevey et al. (2008) are presented in panels a) and b). The VDM reconstructions by Genevey et al. (2008) and Korte & Constable (2005) are shown in panels c) and d), respectively. The dashed lines represent 1-σ statistical errors. For reference, the VADM reconstruction by Knudsen et al. (2008) is presented in all panels.

3.2. Reconstructions of the solar irradiance

As shown in Sect. 2, an estimate of the evolution of the TSI during the Holocene can be obtained from a linear combination of the jth and jth + 1 decadal values of the open flux, which can be computed from cosmogenic isotopes recorded in natural archives. As described in Sect. 2.5, the reconstructions of the open magnetic flux used here are based on the INTCAL04 Δ14C data set (Reimer et al. 2004) employing almost exactly the same algorithm as the one used by Solanki et al. (2004) and Usoskin et al. (2007, 2009). In addition to the uncertainties in computing the TSI from decadal values of the open flux employing SATIRE-T, there are also uncertainties related to the computation of the open flux from cosmogenic isotopes (see Sect. 2.5). The largest uncertainty is posed by the evolution of the geomagnetic field (Usoskin et al. 2006a, 2007; Snowball & Muscheler 2007). Information of the geomagnetic field prior to the era of direct measurements can be retrieved from archeological artifacts, lava flows, and sediments (see Donadini et al. 2010, and references therein). In order to reconstruct the evolution of the geomagnetic field it is necessary to determine at least one component of the magnetic field and the period that the information was stored in the natural archive. Both parameters introduce uncertainties to the determination of the Earth’s magnetic field, as discussed by Donadini et al. (2010). Paleomagnetic reconstructions of the Earth’s dipole moment can be based on different assumptions. Precise determination of the true dipole moment is hardly possible before the era of systematic geomagnetic measurements started in the 19th century. When the quality and spread of data is still good enough, a VDM (virtual dipole moment) can be estimated, assuming that the geomagnetic field can be approximated by a tilted geocentric dipole. When the spatial coverage of paleomagnetic samples becomes too patchy, only the VADM (virtual axial dipole moment) can be estimated by assuming that the geomagnetic axis is aligned with the geographical axis. Both VDM and VADM form a proxy for the true dipole moment (DM), but they include also non-dipole contributions assigned to the dipole strength. Since the VADM assumes the dipole to be aligned with the geographical axis, its values can be strongly distorted if the dipole axis is tilted. VADM and VDM estimates may differ significantly from the true DM, especially on short time scales due to the non-dipole influence in the VADMs caused by a strong regional bias in data distribution (Korte & Constable 2005). Generally VDM models give a better approximation to the true DM when the data coverage is sufficient but tend to underestimate it with spatially limited data sets. In contrast, VADM models tend to overestimate the true moment, especially for insufficient data coverage and tilted dipole axis (Donadini et al. 2010). Therefore, we apply the following approach to estimate the systematic uncertainty of the results due to paleomagnetic reconstructions: we consider the VDM to be a lower bound and VADM an upper bound to the true dipole moment. In order to be able to provide a single reconstruction of TSI over a certain period of time, we simply average the output from VDM and VADM based reconstructions of TSI. We do not claim that the average of the two groups is the most probable value, although the most probable value lies between these bounds. Table 3 presents the geomagnetic models employed here, references, and the time coverage. Note that we have employed the global reconstructions from Genevey et al. (2008). Figure 4 displays the reconstructions of the geomagnetic field intensity, including VADM reconstructions by Knudsen et al. (2008) and Genevey et al. (2008), and the VDM reconstructions by Genevey et al. (2008) and Korte & Constable (2005). Here, we consider this set of reconstructions to estimate the uncertainty of the TSI reconstruction caused by uncertainties in thegeomagnetic field.

Table 3

Paleogeomagnetic models.

thumbnail Fig. 5

a) Reconstructions of the total solar irradiance for the last 3000 years. All reconstructions are based on the INTCAL04 Δ14C data set, but employing different reconstructions of the geomagnetic dipole momentum. The VADM reconstructions by Knudsen et al. (2008) (KN08/blue line) and Genevey et al. (2008) (GN08-8k/thin red line) are presented. In addition, the VDM reconstructions by Genevey et al. (2008)(GN08/green line) and Korte & Constable (2005) (KC05/gray line) are also plotted. For reference, 10-year running means of the SATIRE-T reconstruction are displayed from 1700 to the present (K10/thick red line). b) Estimate of the error due to the uncertainties of the open flux reconstructions. The red dotted line presents the uncertainty of the TSI computed using SATIRE-T. c) Difference between the maximum and minimum values of the TSI computed employing different paleomagnetic models. d) Estimated total uncertainties of the reconstructed TSI employing SATIRE-M model.

3.2.1. Reconstructions for the last 3000 years

In Fig. 5, we present the reconstructions of the TSI for the last 3000 years, which is the period over which all reconstructions of the geomagnetic field overlap. Panel (a) presents the reconstructions of the TSI. The reconstructions based on the VADM models by Knudsen et al. (2008) and Genevey et al. (2008) are presented as blue and thin red lines, respectively. The reconstructions based on the VDM models by Korte & Constable (2005) and Genevey et al. (2008) are displayed as the gray and cyan lines, respectively. This color scheme applies for panels (a), (b) and (d). For reference, the SATIRE-T reconstruction is presented from 1700 to present (thick red line). The different VADM models give results very similar to each other over the full plotted period of time, whereas the VDM models give somewhat higher irradiances, in particular the VDM reconstruction based on Korte & Constable (2005), because of their systematically lower paleomagnetic dipole moment (see Fig. 4, Sect. 3.2 and Appendix A).

thumbnail Fig. 6

Reconstructions of the total solar irradiance for the Holocene. The figure has the same structure as Fig. 5, but for the period 9495 BC to the present. See the caption of Fig. 5 for more information. While the green line in panel c) indicates the number of reconstructions employed to compute the uncertainty of the TSI due to the uncertainty of the evolution of the geomagnetic dipole moment, the black line represents the long-term trend of the uncertainty of the evolution of the geomagnetic dipole moment (y =  −1.7  ×  10 − 5 ∗ t + 0.15 [W/m2]). The period for which only VADM reconstructions are available is indicated in panel d).

Uncertainties in the reconstructions are presented in Figs. 5b − d. Here, we take into account three major sources of uncertainties to the TSI computed from 10-year averaged values of open flux. The first one is the uncertainty in the TSI computed employing the SATIRE-T model (σTSI − SATIRE), which is represented by the red dotted lines in panels (b) and (c) and was discussed in Sect. 3.1. The second source of uncertainty lies in the computation of the open magnetic flux from cosmogenic isotopes (σφopen), which propagates into the TSI employing the procedure described in Sect. 2. Panel (b) presents the estimate of the uncertainty of the j data point due to the open flux (σTSIφopen2=(aTSIσφopen)2+(bTSIσφopen)2\hbox{$\sigma_{{\rm TSI}-\phi_{\rm open}}^2 = (a_{\rm TSI}\sigma_{\phi_{\rm open}})^2 + (b_{\rm TSI}\sigma_{\phi_{\rm open}})^2$}). The third source of uncertainty is the discrepancy between the reconstructions based on conceptually different paleomagnetic models (σTSI − GM). Panel (c) shows the difference between the maximum and minimum values of the jth TSI value computed employing different paleo-geomagnetic models (σTSI − GM). Panel (d) presents the total uncertainty that is given by the sum of the evaluated uncertainties (σtotal2=σTSISATIRE2+σTSIφopen2+σTSIGM2\hbox{$\sigma_{\rm total}^2 = \sigma_{\rm TSI-SATIRE}^2 + \sigma_{{\rm TSI}-\phi_{\rm open}}^2 + \sigma_{\rm TSI-GM}^2$}). For this period, the major sources of uncertainty are related to the irradiance model itself and to the conceptually different paleomagnetic models. The uncertainty introduced by the open flux is much lower. Moreover, the overall uncertainty during this period does not exceed 0.4 W/m2. Note that we have followed common practice and have assumed that there were no significant changes in ocean mixing during the Holocene (Stuiver 1991; Usoskin & Kromer 2005).

A summary of the statistical properties of the TSI reconstructions for the last 3000 years is plotted in Fig. A.1 of the Appendix. From this figure and statistical tests we found that the distributions do not correspond to normal distributions. Additionally, reconstructions from the same group (VADM or VDM) have the same median values. However, the reconstructions belonging to different groups disagree from each other (see Appendix A).

In the model, we assume that the network filling factor (αn) is proportional to the network magnetic flux (φopen + φeph). The proportionality coefficient (1/φsat,n) relates the physical properties of the network emission to the magnetic flux. However, it depends on the spatial resolution of the observations employed to compute it as well as on the estimate of the contribution of the bright components to the TSI. The SATIRE-T model suggests that the increase in the network emission since the Maunder Minimum is approximately 0.7 W/m2. This value is close to the increase of 0.9 W/m2 estimated by Fröhlich (2009) and employed by Steinhilber et al. (2009). On the other hand, the reconstruction of the bright component of the TSI by Wenzler et al. (2006) from magnetograms and continuum disk images indicates a contribution of approximately 0.2 W/m2 during the minima between cycles 21 − 22 and 22 − 23. This comparison stresses the need for higher resolution and more sensitive observations of the solar surface magnetic field, such as now being provided by the observations by the Helioseismic and Magnetic Imager (HMI) instrument on board of the Solar Dynamics Observatory (SDO). Such data are expected to improve the estimate of the contribution of the bright components of the solar atmosphere to the TSI.

thumbnail Fig. 7

a) Comparison between the TSI reconstruction by Steinhilber et al. (2009) (S09; VADM; green line) and two TSI reconstructions obtained in this paper based on VADM and VDM paleomagnetic reconstructions by Knudsen et al. (2008) (KN08; VADM; blue line) and Korte & Constable (2005) (KC05; VDM; gray line), respectively. For reference, 10-year running means of the SATIRE-T reconstruction are displayed from aprroximately 1640 to the present (K10, thick red line). b) Differences between 120-year running means of the models based on 14C and 10Be records. c) Scatter plot TSI-S09 versus TSI-KN08. d) Scatter plot TSI-S09 versus TSI-KC05.

thumbnail Fig. 8

Comparision between the TSI reconstructions based on 14C and 14Be from 500 AD to the present. Except for the considered time interval the structure of the figure is identical to panel a) of Fig. 7.

3.2.2. Reconstructions for the Holocene

Reconstructions of TSI for the whole Holocene are presented in Fig. 6. The structure and color scheme are the same of the Fig. 5. The first 4500 years are covered by the reconstructions based on the VADM values published by Knudsen et al. (2008) and partially by the one based on the VADM presented by Genevey et al. (2008). VADM and VDM reconstructions are jointly available from approximately 5000 BC to 1700 AD, which is the period covered by the Korte & Constable (2005) reconstruction.

We evaluate the uncertainties for these reconstructions in the same way as we described in the previous section, by considering three major sources. The uncertainties of each reconstruction due to the uncertainties of the computation of the open flux are shown in panel (b). It is noticeable that the uncertainties are higher in the early part of the time series. The uncertainty estimated for the KN08 reconstruction decreases after approximately 6000 BC to the level of the other reconstructions. In order to take into account the uncertainty related to the conceptually different paleomagnetic reconstructions, we consider separately the periods that VADM and VDM reconstructions are and are not jointly available. For the period that they are jointly available, i.e. after 5000 BC, we follow the method described in the previous section. For the period prior to jointly reconstructions, we linearly extrapolate the trends bound between 5000 BC and 1700 AD. The results are presented in panel (c). The dashed green line in this panel displays the number of reconstructions employed to compute the maximum and minimum values (scale on the right hand side). While the uncertainty due to the paleomagnetic models is on average at the same level as the uncertainty of the SATIRE-T model from 5000 BC to approximately 1000 BC, during the last 3000 years the uncertainty is lower. The total uncertainty for each reconstruction, which was computed as presented in Sect. 3.2.1, is shown in panel (d). The overall uncertainty is about 0.4 W/m2.

As can be gathered from Fig. A.2 the TSI distributions on this longer time scale are significant closer to Gaussians. The difference to the more skewed distributions found for the last 3000 years has to do with the irregular distributions of grand minima.

thumbnail Fig. 9

Comparison between the evolutions of the TSI after the grand minimum observed in approximately 765 BC and after the Maunder Minimum. The figure follows the format of Fig. 5a. The 10-year averaged SATIRE-T reconstruction (thick red line) is shifted in time by 2465 years in order to match the minimum value observed in 765 BC. The dashed blue line connects the beginning and the end of the SATIRE-T reconstruction and emphasizes the roughly linear rise in irradiance sind the end of the Maunder Minimum.

3.2.3. Comparison between TSI reconstructions based on 14C and 10Be records

In this section, we compare the TSI reconstruction based on 10Be records presented by Steinhilber et al. (2009) and two TSI reconstructions based on 14C computed in this paper. Figure 7a shows the TSI reconstructions by Steinhilber et al. (2009) (S09; VADM; green line) and the reconstructions based on the Knudsen et al. (2008) (KN08; VADM; blue line) and Korte & Constable (2005) (KC05; VDM; gray line) paleomagnetic models since approximately 9500 BC. For reference, 10-year running means of the SATIRE-T reconstruction are displayed from approximately 1640 AD to the present (K10, thick red line). The S09 reconstruction, which was obtained from 40-year running means of 10Be, was resampled to decadal values employing a linear interplolation. An offset of 1365.6 W/m2 was introduced to match the present value of the TSI (PMOD composite).

thumbnail Fig. 10

Evolution of the TSI during the period of high solar activity from 3000 BC to 1000 BC. The figure follows the format of Fig. 5a.

Table 4

SATIRE models employed to construct a representative TSI time series for the Holocene.

thumbnail Fig. 11

TSI weighted reconstruction since approximately 9500 BC. In order to provide a better visualization, the evolution since 1000 BC is displayed in panel b). The filled gray band represents region limited by the KN08-VADM and KC05-VDM reconstructions. For reference, the red lines represent the 10-year averaged reconstruction by Krivova et al. (2010a).

The TSI reconstructions based on different approaches are in the same variability range, while producing some discrepancies at mid to long-term time scales. Although the values of the three time series lie in the same interval and display similar trends prior to 1000 BC, the reconstructions based on 14C and 10Be diverge during the last 3000 years. Figure 7b offers a better visualization of the structural difference between the time series. In this figure, the differences between 120-year running means of the models based on 14C and 10Be records are presented. For the period analyzed, the standard deviation of the differences KN08-S09 and KC05-S09 are approximately 0.32 and 0.34 W/m2, respectively. Note that the differences between 14C and 10Be reconstructions are larger than the differences due to the uncertainty of the geomagnetic field. Panels (c) and (d) of Fig. 7 show the scatter plots of the S09 reconstruction versus the KN08 and KC05 reconstructions, respectively. The correlation coefficients are indicated in the panels. As a consequence of the structural differences of the time series, the correlation coefficients between the reconstructions based on 14C and 10Be records are relatively low.

During the last 3000 years the S09 reconstruction suggests a higher level of TSI compared to the KN08 and KC09 reconstructions. For a better visualization, we present in Fig. 8 the evolution of the time series from 500 AD to the present. Near the end of the Maunder Minimum, the difference is approximately 0.5 W/m2 if we require the TSI reconstruction by Steinhilber et al. (2009) to match the present value of the TSI. Since the two reconstructions are based on different approaches and cosmogenic isotopes, it is difficult to identify the major source of the discrepancy. Different long-term trends in the 14C and 10Be data certainly play an important role (Vonmoos et al. 2006; Usoskin et al. 2009), but methodological differences may also contribute. Note that the TSI reconstructed by S09 also shows a much flatter trend than that obtained by Krivova et al. (2010a) from telescopically measured group sunspot number (K10).

3.2.4. Discussion of examples of reconstructed irradiance trends

The reconstructions indicate that approximately 2450 years ago the TSI was at the same level as is observed at the present time. Figure 9 shows the TSI reconstructions from 1100 BC to 100 AD. According to the VDM and VADM reconstructions, the peak value ranges from approximately 1366.0 to 1366.2 W/m2, respectively. This peak was reached about 330 years after the preceding grand minimum, which was observed in 765 BC. The TSI values during this grand minima computed employing different paleomagnetic reconstructions ranges from 1364.6 to 1364.8 W/m2. This grand minimum was previously discussed by Eddy (1977b,a), Stuiver (1980), Stuiver & Braziunas (1989), Goslar (2003) and Usoskin et al. (2007). For reference, the decadally averaged SATIRE-T reconstruction (thick red line) is shifted 2455 years to match the minimum value observed in 765 BC. The dashed blue line connects the beginning and the end of the SATIRE-T reconstruction. The amplitude and time scale of the TSI evolution from the previous grand minimum, which occurred at approximately 765 BC, to the grand maximum is similar to the amplitude and time scale reconstructed from the end of Maunder Minimum to the present. However, there are also noticeable differences of the evolution of the two periods. For roughly 200 years after the 765 BC grand minimum the TSI is higher than for the same interval after the Maunder Minimum. This maximum around 430 BC was followed by the next grand minimum, which is seen as a large and very rapid decrease of the TSI. This extended minimum lasted for approximately 150 years. The minimum value of the TSI, which ranges from 1364.8 to 1365.0 W/m2, is reached in 345 BC. Note that all 3 plotted grand minima have the same TSI level within the uncertainty limits. The similarities between the evolution of the TSI for the ancient and the present periods cannot be used to predict the occurrence of a grand minimum in the near future, but shows that very rapid and deep drops in the irradiance can occur starting from a high level of TSI, such as the present.

The reconstructions also suggest that from about 3000 BC to 2000 BC the peaks in solar irradiance exhibited values at a similar level as observed at the present time. Figure 10 shows the reconstructions and uncertainties from 3000 BC to 1000 BC. As discussed in the previous subsection, the TSI reconstruction based on the VDM reconstruction (Korte & Constable 2005) displays higher values, which are for some periods slightly higher than the averaged value observed during the last three solar cycles. The structure of the signal during this period is clearly different from the episode described in the previous paragraph. With the exception of a period close to 2860 BC, pronounced grand minima are not observed, although large amplitude oscillations are present. Obviously, the TSI exhibited a rather different behavior during that period of time than in the last 400 years.

3.3. Construction of a TSI time series for investigations of the climate evolution during the Holocene

One of the main motivations for reconstructing the TSI over the Holocene is to produce a consistent record of the external heat input into the coupled atmosphere/ocean/land system that is needed to investigate the evolution of the Earth’s climate employing general circulation models (Schmidt 2010; Schmidt et al. 2010). However, in the previous sections we have produced multiple TSI time series whose validity and length depends on that of the particular geomagnetic field reconstruction underlying it. For use with climate models, however, it is desirable to construct an optimal time series extending through the whole Holocene. In this subsection, we describe a procedure for combining the SATIRE reconstructions during different periods of time. This is not completely trivial due to the different lengths of the time series we are averaging and the fact that they partly have different median values.

In order to produce a representative time series, we divide the Holocene into four periods according to the available data for each period. Table 4 summarizes the periods considered and the data available for each period. After the end of the Maunder Minimum we compute daily values, while prior to the end of the Maunder Minimum we compute 10-year averages. For the period for which both solar disk magnetograms and continuum images are available (period 1) we employ the SATIRE-S reconstruction. SATIRE-T reconstruction is used from the beginning of the Maunder Minimum (~1640 AD) to 1978 AD. Prior to 1640 AD reconstructions based on cosmogenic isotopes are used. As VDM paleomagnetic reconstructions are available since approximately 5000 BC, we have to treat separately periods 3 and 4. We emphasize that the reconstructions based on different proxies have different time resolutions.

For period 3, VADM and VDM reconstructions are available. The TSI reconstructions based on conceptually different paleomagnetic models differ significantly from each other. As discussed in Sect. 3.2, the paleomagnetic reconstructions based on the VADM and VDM models provide the upper and lower limits of the evolution of the geomagnetic field, respectively. Consequently, we need to combine them to produce a representative reconstruction for this period. As we have discussed previously, the median values of the reconstructions for the last 3000 years that belong to a given group are not significantly different. In this way, assuming that this result holds for the whole period starting in 5000 BC, we can compute the final TSI at time step j as an average of the two methods VADM and VDM taking into account that they represent lower (VDM) and upper (VADM) bounds to the true dipole moment. We do not claim that the average of the two groups is the true value, but since the true value lies between these bounds, if (the true value) is more likely to lie close to the average value than to either of the bounds. Here we employ the reconstructions based on the KC05 and KN08 dipole moments that are representative of the VDM and VADM groups, respectively.

For period 4, only VADM reconstructions are available. Here we employ the KN08 reconstruction as the reference. To be consistent with other periods, it is necessary to introduce an offset that can be estimated from the difference of the mean values of the KN08-VADM and KC05-VDM reconstructions for the period 3. In this way, we implicitly assume that the VDM paleomagnetic models would have lower values than the VADM paleomagnetic models in period 4 as they have had in more recent times as well. The numerical value of the offset is approximately 0.05 W/m2. In this way, discontinuities are not introduced in the interface between periods 3 and 4. Note that the upper bound limit in Fig. 11 is obtained by adding 2·offset to the VADM reconstruction. For more recent times the SATIRE-T and SATIRE-S based reconstructions are appended to the 14C based TSI reconstructions (see Table 4). The final TSI reconstruction is presented in Fig. 11. The periods considered and listed in Table 4 are indicated in the figure. To avoid overloading the figure, 10-year averages of the SATIRE-T model are plotted.

4. Concluding remarks

In order to understand the evolution of the Earth’s climate during the Holocene the variability of the internal and external drivers needs to be quantified. Among the external drivers, the solar irradiance plays a prominent role. Here, we have reconstructed the TSI over the Holocene. This is an extension of the SATIRE family of models to a longer time scale. It represents the first reconstruction of TSI during pre-telescopic times that is based on a physically consistent treatment, which builds on the assumption that the magnetic field is responsible for all TSI changes. This assumption implies that our model cannot handle changes in TSI that do not originate in surface magnetism (e.g. convective efficiency changes unrelated to surface magnetism or effects associated with r modes, i.e. oscillatory modes related with rotation, or changes associated with magnetic fields buried deep in the convection zone). On the other hand, it allows us to rigorously derive the proper relationship between irradiance and magnetic flux. Earlier reconstructions of TSI based on multi-decadal averages of cosmogenic data (e.g. Steinhilber et al. 2009) are based on linear extrapolations of an empirical relationship between open flux and TSI. However, our computations do not fully support a simple one-to-one relationship between TSI and solar open magnetic flux on decadal time-scale, since we find that magnetic flux from two cycles contributes to the TSI from one cycle. The reconstructions carried out here are based on this consistently derived relationship. Another major difference between the present and previous reconstructions of TSI over the Holocene lies in our careful assessment of the uncertainties in reconstructions of the geomagnetic field by considering the available geomagnetic field reconstructions that reaches back at least 3000 years. The various reconstructions of the geomagnetic field can lead to significant differences in TSI at earlier times.

The measured solar irradiance can be reproduced on the basis of the assumption that its variability is related to the evolution of the magnetic features lying on the solar surface (Krivova et al. 2003; Wenzler et al. 2006). However, the spatial and temporal resolution of records of magnetic features on the solar disk degrades progressively as we go back in time. For example, the solar irradiance can be reproduced with high accuracy from magnetograms and continuum images of the solar disk employing a physics-based model (SATIRE-S Krivova et al. 2003). However, these records are jointly available only for the last three solar cycles. To deal with this limitation, models of the evolution of the solar magnetic field can be employed to replace the direct observations of the evolution of the magnetic features on the solar disk. Most of these models are based on the records of sunspot number, which are available since approximately 1610 AD (e.g., SATIRE-T, Balmaceda 2007; Krivova et al. 2007, 2010a; and other models, Lean et al. 1995; Solanki & Fligge 1999; Wang et al. 2005; Tapping et al. 2007; Crouch et al. 2008), but contain no information on the spatial distribution of the magnetic features on the solar surface and suffer from loss of temporal resolution at earlier times. On longer time-scales, the evolution of the solar magnetic flux can be evaluated based on records of cosmogenic isotopes stored in natural archives, although generally at the cost of a further loss of temporal resolution. Here, we have shown that 10-year averages of the total and spectral irradiance can be obtained by employing the magnetic flux computed from cosmogenic isotopes as the input of the SATIRE model (SATIRE-M, where “M” stands for millennial time scales). Note that the core of the irradiance reconstruction method is the same for all variations of the SATIRE model. The solar atmosphere is divided into components: sunspots, umbrae and penumbrae, faculae, network and the quiet Sun. The evolution of the spectrum is described by the filling factors on different parts of the solar surface of each component and its time independent brightness. Only the evolution of the filling factors of the surface components is computed according to the available data.

The reconstructed 10-year averaged TSI displays variations over a range of approximately 1.5 W/m2. This is slightly larger than the 1.3 W/m2 between the Maunder Minimum and the current grand maximum found by Krivova et al. (2010a), mainly because an earlier Grand minimum (around 7500 BC) displays a lower TSI than the Maunder Minimum and the Grand Maximum at the beginning of the time series displays a larger TSI than the just ending Grand Maximum. We point out that this range is uncertain as there are significant uncertainties in the computation of the open flux at the earlier times. The major sources of uncertainty of the solar irradiance reconstructions during the Holocene come from the simplifications employed to derive the evolution of the magnetic flux emerging in AR and ER based on decadal values of the open flux and from the evolution of the Earth’s magnetic field. In order to evaluate the uncertainties due to the evolution of the Earth’s magnetic field, we compared an ensemble of reconstructions that are based on two conceptually different approaches to computing the evolution of the Earth’s magnetic dipole. For the period that all reconstructions are jointly available (i.e. during the last 3000 years) the two groups show significantly different median values, reflecting the fact that reconstructions of the virtual dipole moment may underestimate the true dipole and, consequently, overestimate the solar parameters. The opposite is expected for the virtual axial dipole moment reconstructions. In addition to these uncertainties, we recall that the reconstructions are based on the assumption that carbon cycle parameters have not changed significantly throughout the Holocene, which may not be exactly valid, particularly in the early Holocene, as discussed previously by Solanki et al. (2004), Ellison et al. (2006), Usoskin et al. (2007, 2009), and others. We stress that the TSI reconstructions presented in this paper are based on one individual record of 14C, and, as shown in Sect. 3.2.3, using 10Be records would probably lead to another long-term behaviour.

We have also produced a reconstruction that is a careful combination of the various TSI reconstructions we have carried out for the period prior to the Maunder Minimum, i.e. before 1640 AD. Between 1640 AD and 1974 AD it corresponds to the SATIRE-T reconstruction of Krivova et al. (2010a), while after that we employ a SATIRE-S reconstruction based on Kitt Peak (Wenzler et al. 2006) and MDI (following Krivova et al. 2003) magnetograms. This combined reconstruction, which is our best estimate of the TSI during the Holocene, is available in electronic form at the CDS. Alternatively, the irradianced data can also be downloaded from the website: http://www.mps.mpg.de/projects/sun-climate/data.html.

The present analysis is in agreement with earlier reconstructions of the solar activity during the Holocene regarding the occurrence of extended periods of low and high solar activity (e.g. Usoskin et al. 2007). These coincide with periods of low and high TSI. We noted some similarities between the evolution of the solar activity from the Maunder Minimum to the present and an episode observed approximately 2800 years ago. Both episodes displayed changes of the solar irradiance of about 1.3 W/m2 in an interval of approximately 300 years. The ancient episode ended rather abruptly in a grand minimum, which emphasizes that a very rapid drop in the irradiance can occur starting from a high level of TSI, such as the present one. However, we stress that we cannot employ these similarities to predict the occurrence of a grand minimum in the near future.

Online material

Appendix A: Statistical properties of the TSI reconstructions

Panels (a) to (d) of Fig. A.1 display the probability distributions of the 10-year averaged TSI values. Each panel shows the probability distribution functions (PDF) for the normal distribution computed employing the respective mean and standard deviation values that are also given in each panel. In order to verify if the reconstructed TSI values are normally distributed, we performed a test of the default null hypothesis that the values come from a normal distribution, against the hypothesis that they do not come from this distribution. The test statistic is T=max[abs(ecdf(TSI)ncdf(TSI))],\appendix \setcounter{section}{1} \begin{equation} T = \max \left[ {abs} \left (ecdf({\rm TSI}) -ncdf({\rm TSI}) \right) \right], \end{equation}(A.1)where ecdf(TSI) is the empirical cumulative distribution function (CDF) estimated from the TSI distribution and ncdf(TSI) is the normal CDF with mean and standard deviation equal to the mean and standard deviation of the TSI distribution. Recall that the value of the ecdf at a given value of the TSI is obtained from the PDF by integrating the latter from its lower end up to that particular TSI value. The T values are indicated in the respective panels. The null hypothesis is rejected at a given significance level when T is greater than a critical value. The test employed here uses a table of critical values computed using Monte-Carlo simulations, where the critical value for the significance level of 0.05 is 0.0547. We found that the null hypothesis of all reconstructions is rejected at this significance level, i.e., the distributions do not correspond to normal distributions. This is supported by a visual inspection which reveals the strong asymmetry in all the the empirical distributions. Besides the main peak located above 1365 W/m2 there is also a lower (and statistically less significant) secondary peak below 1365 W/m2, which corresponds to the relatively common grand minima in this period of time.

Figure A.1e compares the median values (red lines). Additionally, the upper and lower boundaries of the boxes display the respective quartile values. The dashed lines show the regions extending from the boundaries of each box to the most extreme values within 1.5 times the interquartile range. As expected, the VDM reconstructions display higher median values

than the VADM reconstructions. Furthermore, as also suggested by the empirical CDFs, the VDM distributions are shifted in relation to the VADM distributions, and also display a larger range of irradiance values.

We employ the Wilcoxon test to verify if two reconstructions come from identical distributions with equal medians, against the alternative that they do not have equal medians. Note that the distributions need not necessarily be normal. Table A.1 summarizes the results. The reconstructions from the same group (VADM or VDM) have the same median values. However, the reconstructions belonging to different groups disagree from each other.

The statistical properties of the TSI reconstructions for the whole Holocene based on the paleomagnetic reconstructions by Knudsen et al. (2008), Genevey et al. (2008), and Korte & Constable (2005) are presented in panels (a) − (c) of Fig. A.2, respectively. As in Fig. A.1, the red lines represent the probability distribution functions for the normal distribution. Note that each time series covers a different time span (see Table 3). Tests of the null hypothesis that the values come from a normal distribution, as discussed in Sect. 3.2.1, indicate that it is not possible to reject this hypothesis at the significance level of 0.05 for the reconstructions based on the Genevey et al. (2008) paleomagnetic reconstruction. A visual inspection shows that the reconstruction based on the paleomagnetic model by Korte & Constable (2005), that extends from 5000 BC to 1700, gives a slightly asymmetrical distribution. This behavior reflects the more frequent occurrence of grand minima in the later part of the time series, as can be gathered by comparing with the far more asymmetric distributions in Fig. A.1.

Table A.1

Wilcoxon rank sum test of the distribution medians.

thumbnail Fig. A.1

Comparison between the statistical properties of the TSI reconstructions for the period 1000 BC to 1700 AD. Panels a) to d) present the distributions of the TSI reconstructions. The mean, median and standard deviation values are indicated in the panels. The red lines represent the probability distributions functions for the normal distribution. Panel e) allows a visual comparison between the distributions. The median (red line), lower and upper quartile values of the distributions (lower and upper ends of the blue boxes) are plotted in each box. The black lines extend from the end of each box to the most extreme values within 1.5 times the interquartile range. Note that the interquartile range, which is also a measure of the statistical dispersion of the data, is equal to the difference between the third and first quartiles.

thumbnail Fig. A.2

Comparison between the statistical properties of the TSI reconstructions during the Holocene. The time coverage of each time series is presented in Table 3. Panels a) to d) present the distributions of the TSI reconstructions. The red line represents the normal distribution computed employing the sample mean and standard deviation values of as the overplotted sample (the values are indicated in each panel). Additionally, the test statistic (T) and the critical (Tcritical) values are listed.

Acknowledgments

We are grateful to C. Fröhlich and F. Steinhilber for providing the PMOD composite TSI and the TSI reconstruction based on 10Be records, respectively. We are also thankful to M. Lockwood for providing the open flux reconstructions based on the aa-index. This work was supported by the Deutsche Forschungs-gemeinschaft, DFG project number SO 711/1-2, by the WCU grant No. R31-10016 funded by the Korean Ministry of Education, Science and Technology, and by the European Commission’s Seventh Framework Programme (FP7/2007-2013) under the grant agreement No. 218816 (SOTERIA project).

References

  1. Arge, C. N., Hildner, E., Pizzo, V. J., & Harvey, J. W. 2002, J. Geophys. Res., 107, 1319 [NASA ADS] [CrossRef] [Google Scholar]
  2. Balmaceda, L. A. 2007, Ph.D. Thesis, Max Planck Inst. Sol. Syst. Res. [Google Scholar]
  3. Balmaceda, L. A., Solanki, S. K., Krivova, N. A., & Foster, S. 2009, J. Geophys. Res., 114, A07104 [NASA ADS] [CrossRef] [Google Scholar]
  4. Bond, G., Kromer, B., Beer, J., et al. 2001, Science, 294, 2130 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  5. Camp, C. D., & Tung, K. K. 2007, Geophys. Res. Lett., 34, L14703 [NASA ADS] [CrossRef] [Google Scholar]
  6. Castagnoli, G., & Lal, D. 1980, Radiocarbon, 22, 133 [Google Scholar]
  7. Crouch, A. D., Charbonneau, P., Beaubien, G., & Paquin-Ricard, D. 2008, ApJ, 677, 723 [Google Scholar]
  8. Delaygue, G., & Bard, E. 2010, Clim. Dyn., 1 [Google Scholar]
  9. Donadini, F., Korte, M., & Constable, C. 2010, Space Sci. Rev., 155, 219 [NASA ADS] [CrossRef] [Google Scholar]
  10. Eddy, J. 1977a, Climatic Change, 1, 173 [Google Scholar]
  11. Eddy, J. 1977b, Sci. Am., 236, 80 [NASA ADS] [CrossRef] [Google Scholar]
  12. Ellison, C. R. W., Chapman, M. R., & Hall, I. R. 2006, Science, 312, 1929 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  13. Fligge, M., & Solanki, S. K. 1997, Sol. Phys., 173, 427 [NASA ADS] [CrossRef] [Google Scholar]
  14. Fligge, M., & Solanki, S. K. 2000, Geophys. Res. Lett., 27, 2157 [NASA ADS] [CrossRef] [Google Scholar]
  15. Fröhlich, C. 2006, Space Sci. Rev., 125, 53 [Google Scholar]
  16. Fröhlich, C. 2009, A&A, 501, L27 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Genevey, A., Gallet, Y., Constable, C. G., Korte, M., & Hulot, G. 2008, Geochem. Geophys. Geosyst., 9, Q04038 [NASA ADS] [CrossRef] [Google Scholar]
  18. Goslar, T. 2003, in PAGES News, Past Global Changes, 11, 12 [Google Scholar]
  19. Haigh, J. D. 2001, Science, 294, 2109 [CrossRef] [PubMed] [Google Scholar]
  20. Haigh, J. D. 2007, Liv. Rev. Sol. Phys., URL cited on Aug./2010: http://www.livingreviews.org/lrsp-2007-2, 4 [Google Scholar]
  21. Hansen, J. E. 2000, Space Sci. Rev., 94, 349 [NASA ADS] [CrossRef] [Google Scholar]
  22. Hansen, J., Sato, M., Nazarenko, L., et al. 2002, J. Geophys. Res., 107, 4347 [CrossRef] [Google Scholar]
  23. Hansen, J., Sato, M., Ruedy, R., et al. 2005, J. Geophys. Res., 110, D18104 [NASA ADS] [CrossRef] [Google Scholar]
  24. Hathaway, D. H. 2010, Liv. Rev. Sol. Phys., URL cited on Aug./2010): http://www.livingreviews.org/lrsp-2010-1, 7 [Google Scholar]
  25. Jungclaus, J., Lorenz, S., Timmreck, C., et al. 2010, Clim. Past Discuss., 6, 1009 [Google Scholar]
  26. Knudsen, M. F., Riisager, P., Donadini, F., et al. 2008, Earth Planet. Sci. Lett., 272, 319 [NASA ADS] [CrossRef] [Google Scholar]
  27. Kopp, G., & Lean, J. L. 2011, Geophys. Res. Lett., 38, L01706 [Google Scholar]
  28. Kopp, G., Lawrence, G., & Rottman, G. 2005, Sol. Phys., 230, 129 [NASA ADS] [CrossRef] [Google Scholar]
  29. Korte, M., & Constable, C. 2005, Earth Planet. Sci. Lett., 236, 348 [CrossRef] [Google Scholar]
  30. Krivova, N. A., Solanki, S. K., Fligge, M., & Unruh, Y. C. 2003, A&A, 399, L1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  31. Krivova, N. A., Balmaceda, L., & Solanki, S. K. 2007, A&A, 467, 335 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  32. Krivova, N. A., Solanki, S. K., Wenzler, T., & Podlipnik, B. 2009, J. Geophys. Res., 114, D00I04 [CrossRef] [Google Scholar]
  33. Krivova, N., Vieira, L., & Solanki, S. 2010a, J. Geophys. Res., 115, A12112 [NASA ADS] [CrossRef] [Google Scholar]
  34. Krivova, N. A., Solanki, S. K., & Unruh, Y. C. 2010b, J. Atmos. Sol.-Terr. Phys., in press [Google Scholar]
  35. Kurucz, R. 1993, ATLAS9 Stellar Atmosphere Programs and 2 km s-1 grid, Kurucz CD-ROM No. 13, Cambridge, Mass., Smithsonian Astrophysical Observatory [Google Scholar]
  36. Lean, J., Beer, J., & Bradley, R. 1995, Geophys. Res. Lett., 22, 3195 [NASA ADS] [CrossRef] [Google Scholar]
  37. Lockwood, M., Owens, M. J., & Rouillard, A. P. 2009a, J. Geophys. Res., 114, A11103 [NASA ADS] [CrossRef] [Google Scholar]
  38. Lockwood, M., Owens, M. J., & Rouillard, A. P. 2009b, J. Geophys. Res., 114, A11104 [NASA ADS] [CrossRef] [Google Scholar]
  39. Mosteller, F., & Tukey, J. 1977, Data Analysis and Regression (Addison-Wesley) [Google Scholar]
  40. Muscheler, R., Joos, F., Beer, J., et al. 2007, Quat. Sci. Rev., 26, 82 [Google Scholar]
  41. Reimer, P., Baillie, M. G. L., Bard, E., et al. 2004, Radiocarbon, 46, 1029 [Google Scholar]
  42. Rodbell, D. T., Seltzer, G. O., Anderson, D. M., et al. 1999, Science, 283, 516 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  43. Schmidt, G. 2010, J. Quat. Sci., 25, 79 [CrossRef] [Google Scholar]
  44. Schmidt, G. A., Jungclaus, J. H., Ammann, C. M., et al. 2010, Geosci. Model Dev. Discuss., 3, 1549 [NASA ADS] [CrossRef] [Google Scholar]
  45. Siegenthaler, Ulrich, & Heimann, M. O. H. 1980, Radiocarbon, 22, 177 [Google Scholar]
  46. Snowball, I., & Muscheler, R. 2007, The Holocene, 17, 851 [CrossRef] [Google Scholar]
  47. Solanki, S. K., & Fligge, M. 1999, Geophys. Res. Lett., 26, 2465 [NASA ADS] [CrossRef] [Google Scholar]
  48. Solanki, S. K., Schüssler, M., & Fligge, M. 2000, Nature, 408, 445 [Google Scholar]
  49. Solanki, S. K., Usoskin, I. G., Kromer, B., Schüssler, M., & Beer, J. 2004, Nature, 431, 1084 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  50. Solomon, S. 2007, Climate Change 2007: The Physical Science Basis, Contribution of Working Group I to the 4th IPCC Report (Cambridge, USA: Cambridge University Press) [Google Scholar]
  51. Steinhilber, F., Abreu, J. A., & Beer, J. 2008, Astrophys. Space Sci. Trans., 4, 1 [NASA ADS] [CrossRef] [Google Scholar]
  52. Steinhilber, F., Beer, J., & Fröhlich, C. 2009, Geophys. Res. Lett., 36, L19704 [Google Scholar]
  53. Steinhilber, F., Abreu, J. A., Beer, J., & McCracken, K. G. 2010, J. Geophys. Res., 115, A01104 [NASA ADS] [Google Scholar]
  54. Stuiver, M. 1980, Nature, 286, 868 [NASA ADS] [CrossRef] [Google Scholar]
  55. Stuiver, M. 1991, Quat. Res., 35, 1 [Google Scholar]
  56. Stuiver, M., & Braziunas, T. F. 1989, Nature, 338, 405 [NASA ADS] [CrossRef] [Google Scholar]
  57. Stuiver, M., Reimer, P. J., Bard, E., et al. 1998, Radiocarbon, 40, 1041 [Google Scholar]
  58. Tapping, K., Boteler, D., Charbonneau, P., et al. 2007, Sol. Phys., 246, 309 [Google Scholar]
  59. Taricco, C., Bhandari, N., Cane, D., Colombetti, P., & Verma, N. 2006, J. Geophys. Res., 111, A08102 [Google Scholar]
  60. Trenberth, K. E. 1990, Bull. Am. Meteorol. Soc., 71, 988 [NASA ADS] [CrossRef] [Google Scholar]
  61. Unruh, Y. C., Solanki, S. K., & Fligge, M. 1999, A&A, 345, 635 [NASA ADS] [Google Scholar]
  62. Usoskin, I. G., Horiuchi, K., Solanki, S. K., Kovaltsov, G. A., & Bard, E. 2009, J. Geophys. Res., 114, A03112 [NASA ADS] [CrossRef] [Google Scholar]
  63. Usoskin, I. G. & Kromer, B. 2005, Radiocarbon, 47, 31 [Google Scholar]
  64. Usoskin, I. G., Mursula, K., Solanki, S. K., Schüssler, M., & Kovaltsov, G. A. 2002, J. Geophys. Res., 107, 1374 [NASA ADS] [CrossRef] [Google Scholar]
  65. Usoskin, I. G., Solanki, S. K., & Kovaltsov, G. A. 2007, A&A, 471, 301 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  66. Usoskin, I. G., Mursula, K., Solanki, S. K., Schüssler, M., & Alanko, K. 2004, A&A, 413, 745 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  67. Usoskin, I. G., Solanki, S. K., Schüssler, M., Mursula, K., & Alanko, K. 2003, Phys. Rev. Lett., 91, 211101 [Google Scholar]
  68. Usoskin, I. G., Schüssler, M., Solanki, S. K., & Mursula, K. 2005, J. Geophys. Res., 110, A10102 [NASA ADS] [CrossRef] [Google Scholar]
  69. Usoskin, I. G., Solanki, S. K., & Korte, M. 2006a, Geophys. Res. Lett., 33, L08103 [CrossRef] [Google Scholar]
  70. Usoskin, I. G., Solanki, S. K., Taricco, C., Bhandari, N., & Kovaltsov, G. A. 2006b, A&A, 457, L25 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  71. Vieira, L. E. A., & Solanki, S. K. 2010, A&A, 509, A100 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  72. Vieira, L. E. A., Krivova, N., Solanki, S. K., & Balmaceda, L. A. 2008, Eos Trans. AGU, 89, Abstract GP31C [Google Scholar]
  73. Vonmoos, M., Beer, J., & Muscheler, R. 2006, J. Geophys. Res., 111, A10105 [NASA ADS] [CrossRef] [Google Scholar]
  74. Wang, Y.-M., Lean, J. L., & Sheeley, N. R., Jr. 2005, ApJ, 625, 522 [Google Scholar]
  75. Wenzler, T., Solanki, S. K., Krivova, N. A., & Fluri, D. 2004, A&A, 427, 1031 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  76. Wenzler, T., Solanki, S. K., & Krivova, N. A. 2005, A&A, 432, 1057 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  77. Wenzler, T., Solanki, S. K., Krivova, N. A., & Fröhlich, C. 2006, A&A, 460, 583 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  78. Zhang, M. H., & Song, H. 2006, Geophys. Res. Lett., 33, l12701 [NASA ADS] [CrossRef] [Google Scholar]

All Tables

Table 1

Parameters of the TSI model presented by Krivova et al. (2010a) and employed here.

Table 2

Coefficients of the model to compute TSI for the Holocene (SATIRE-M).

Table 3

Paleogeomagnetic models.

Table 4

SATIRE models employed to construct a representative TSI time series for the Holocene.

Table A.1

Wilcoxon rank sum test of the distribution medians.

All Figures

thumbnail Fig. 1

a) Evolution of the solar open magnetic flux. The black and red lines present the estimates of the open flux from the geomagentic aa-index and from observations of the interplanetary magnetic field (Lockwood et al. 2009a,b), respectively. The blue line displays the modeled open flux from SATIRE-T. b) 44Ti activity in stony meteorites as a function of the time of fall (see Taricco et al. 2006, and references therein). Error bars refer to 1-sigma uncertainties. The red line corresponds to the modeled 44Ti based on the SATIRE-T open flux according to Usoskin et al. (2006b). The 44Ti activity is given in disintegrations per minute per kilogram (dpm kg-1) of Ni  +  Fe. The dashed red lines correspond to the estimated 1-sigma uncertanty of the model.

In the text
thumbnail Fig. 2

a) Comparison between the observed (PMOD Composite; blue line) and the modeled (SATIRE-T; green line) TSI from Nov./1978 to Dec./2003. b) Difference between the observed and modeled TSI. The gray and black lines present the daily and 90-day averaged values, respectively. c) Distribution of the difference between daily values of the observed and modeled TSI (model error). The standard deviation and the median absolute deviation are indicated in the panel. The red curve represents the corresponding probability distribution function for the normal distribution with the mean and standard deviation of the model error. d) Cumulative distribution function (CDF) for the difference between observed and modeled TSI. e) Probability plot for normal distribution with the mean and standard deviation of the model error. f) Scatter plot of the observed versus modeled TSI. The red line indicates the set of expectation values of the model. The correlation coefficient and the coefficients of the best linear fitting (green line) are presented in the panel. The p-value for testing the hypothesis of no correlation is lower than 0.5, which indicates that the correlation is significant.

In the text
thumbnail Fig. 3

Comparison between 10-year averaged values of the parameters computed employing the SATIRE-M and SATIRE-T models from 1700 AD to 2000 AD. The following scatter plots are presented: a) total magnetic flux; b) total solar irradiance; c) contribution of the bright components (faculae + network) to the total solar irradiance; d) sunspot (umbra + penumbra) contribution to TSI; e) UV integrated in the range 220−240 nm. The red lines present the set of expectation values for the models and the blue lines are the best linear fits. The correlation coefficient and the coefficients of the linear regression are presented in the panels.

In the text
thumbnail Fig. 4

Reconstructions of the geomagnetic field intensity. The VADM reconstructions by Knudsen et al. (2008) and Genevey et al. (2008) are presented in panels a) and b). The VDM reconstructions by Genevey et al. (2008) and Korte & Constable (2005) are shown in panels c) and d), respectively. The dashed lines represent 1-σ statistical errors. For reference, the VADM reconstruction by Knudsen et al. (2008) is presented in all panels.

In the text
thumbnail Fig. 5

a) Reconstructions of the total solar irradiance for the last 3000 years. All reconstructions are based on the INTCAL04 Δ14C data set, but employing different reconstructions of the geomagnetic dipole momentum. The VADM reconstructions by Knudsen et al. (2008) (KN08/blue line) and Genevey et al. (2008) (GN08-8k/thin red line) are presented. In addition, the VDM reconstructions by Genevey et al. (2008)(GN08/green line) and Korte & Constable (2005) (KC05/gray line) are also plotted. For reference, 10-year running means of the SATIRE-T reconstruction are displayed from 1700 to the present (K10/thick red line). b) Estimate of the error due to the uncertainties of the open flux reconstructions. The red dotted line presents the uncertainty of the TSI computed using SATIRE-T. c) Difference between the maximum and minimum values of the TSI computed employing different paleomagnetic models. d) Estimated total uncertainties of the reconstructed TSI employing SATIRE-M model.

In the text
thumbnail Fig. 6

Reconstructions of the total solar irradiance for the Holocene. The figure has the same structure as Fig. 5, but for the period 9495 BC to the present. See the caption of Fig. 5 for more information. While the green line in panel c) indicates the number of reconstructions employed to compute the uncertainty of the TSI due to the uncertainty of the evolution of the geomagnetic dipole moment, the black line represents the long-term trend of the uncertainty of the evolution of the geomagnetic dipole moment (y =  −1.7  ×  10 − 5 ∗ t + 0.15 [W/m2]). The period for which only VADM reconstructions are available is indicated in panel d).

In the text
thumbnail Fig. 7

a) Comparison between the TSI reconstruction by Steinhilber et al. (2009) (S09; VADM; green line) and two TSI reconstructions obtained in this paper based on VADM and VDM paleomagnetic reconstructions by Knudsen et al. (2008) (KN08; VADM; blue line) and Korte & Constable (2005) (KC05; VDM; gray line), respectively. For reference, 10-year running means of the SATIRE-T reconstruction are displayed from aprroximately 1640 to the present (K10, thick red line). b) Differences between 120-year running means of the models based on 14C and 10Be records. c) Scatter plot TSI-S09 versus TSI-KN08. d) Scatter plot TSI-S09 versus TSI-KC05.

In the text
thumbnail Fig. 8

Comparision between the TSI reconstructions based on 14C and 14Be from 500 AD to the present. Except for the considered time interval the structure of the figure is identical to panel a) of Fig. 7.

In the text
thumbnail Fig. 9

Comparison between the evolutions of the TSI after the grand minimum observed in approximately 765 BC and after the Maunder Minimum. The figure follows the format of Fig. 5a. The 10-year averaged SATIRE-T reconstruction (thick red line) is shifted in time by 2465 years in order to match the minimum value observed in 765 BC. The dashed blue line connects the beginning and the end of the SATIRE-T reconstruction and emphasizes the roughly linear rise in irradiance sind the end of the Maunder Minimum.

In the text
thumbnail Fig. 10

Evolution of the TSI during the period of high solar activity from 3000 BC to 1000 BC. The figure follows the format of Fig. 5a.

In the text
thumbnail Fig. 11

TSI weighted reconstruction since approximately 9500 BC. In order to provide a better visualization, the evolution since 1000 BC is displayed in panel b). The filled gray band represents region limited by the KN08-VADM and KC05-VDM reconstructions. For reference, the red lines represent the 10-year averaged reconstruction by Krivova et al. (2010a).

In the text
thumbnail Fig. A.1

Comparison between the statistical properties of the TSI reconstructions for the period 1000 BC to 1700 AD. Panels a) to d) present the distributions of the TSI reconstructions. The mean, median and standard deviation values are indicated in the panels. The red lines represent the probability distributions functions for the normal distribution. Panel e) allows a visual comparison between the distributions. The median (red line), lower and upper quartile values of the distributions (lower and upper ends of the blue boxes) are plotted in each box. The black lines extend from the end of each box to the most extreme values within 1.5 times the interquartile range. Note that the interquartile range, which is also a measure of the statistical dispersion of the data, is equal to the difference between the third and first quartiles.

In the text
thumbnail Fig. A.2

Comparison between the statistical properties of the TSI reconstructions during the Holocene. The time coverage of each time series is presented in Table 3. Panels a) to d) present the distributions of the TSI reconstructions. The red line represents the normal distribution computed employing the sample mean and standard deviation values of as the overplotted sample (the values are indicated in each panel). Additionally, the test statistic (T) and the critical (Tcritical) values are listed.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.