Scholarly article on topic 'A History of Solar Activity over Millennia'

A History of Solar Activity over Millennia Academic research paper on "Earth and related environmental sciences"

Share paper
Academic journal
Living Rev. Solar Phys.

Academic research paper on topic "A History of Solar Activity over Millennia"

Living Rev. Solar Phys., 10, (2013), 1 doi:10.12942/lrsp-2013-1

LIVING ^jp REVIEWS in solar physics

A History of Solar Activity over Millennia

Ilya G. Usoskin

Sodankyla Geophysical Observatory (Oulu unit) FIN-90014 University of Oulu, Finland email:

Accepted: 7 March 2013 Published: 21 March 2013

(Update of lrsp-2008-3)


Presented here is a review of present knowledge of the long-term behavior of solar activity on a multi-millennial timescale, as reconstructed using the indirect proxy method. The concept of solar activity is discussed along with an overview of the special indices used to quantify different aspects of variable solar activity, with special emphasis upon sunspot number.

Over long timescales, quantitative information about past solar activity can only be obtained using a method based upon indirect proxies, such as the cosmogenic isotopes 14C and 10Be in natural stratified archives (e.g., tree rings or ice cores). We give an historical overview of the development of the proxy-based method for past solar-activity reconstruction over millennia, as well as a description of the modern state. Special attention is paid to the verification and cross-calibration of reconstructions. It is argued that this method of cosmogenic isotopes makes a solid basis for studies of solar variability in the past on a long timescale (centuries to millennia) during the Holocene.

A separate section is devoted to reconstructions of strong solar energetic-particle (SEP) events in the past, that suggest that the present-day average SEP flux is broadly consistent with estimates on longer timescales, and that the occurrence of extra-strong events is unlikely.

Finally, the main features of the long-term evolution of solar magnetic activity, including the statistics of grand minima and maxima occurrence, are summarized and their possible implications, especially for solar/stellar dynamo theory, are discussed.

Keywords: solar activity, paleo-astrophysics, cosmogenic isotopes, solar-terrestrial relations, solar physics, long-term reconstructions, solar dynamo

This review is licensed under a Creative Commons Attribution-Non-Commercial 3.0 Germany License.

Imprint / Terms of Use

Living Reviews in Solar Physics is a peer reviewed open access journal published by the Max Planck Institute for Solar System Research, Max-Planck-Str. 2, 37191 Katlenburg-Lindau, Germany. ISSN 1614-4961.

This review is licensed under a Creative Commons Attribution-Non-Commercial 3.0 Germany License: Figures that have been previously published elsewhere may not be reproduced without consent of the original copyright holders.

Because a Living Reviews article can evolve over time, we recommend to cite the article as follows:

Ilya G. Usoskin, "A History of Solar Activity over Millennia", Living Rev. Solar Phys., 10, (2013), 1. URL (accessed <date>):

The date given as <date> then uniquely identifies the version of the article you are referring to.

Article Revisions

Living Reviews supports two ways of keeping its articles up-to-date:

Fast-track revision A fast-track revision provides the author with the opportunity to add short notices of current research results, trends and developments, or important publications to the article. A fast-track revision is refereed by the responsible subject editor. If an article has undergone a fast-track revision, a summary of changes will be listed here.

Major update A major update will include substantial changes and additions and is subject to full external refereeing. It is published with a new publication number.

For detailed documentation of an article's evolution, please refer to the history document of the

article's online version at

21 March 2013: The review has been thoroughly revised and updated. Added Sections 3.8 and 5.1 and 8 new figures (3 were removed). 55 new references have been included (4 were removed).


1 Introduction 7

2 Solar Activity: Concept and Observations 9

2.1 The concept of solar activity ............................................................9

2.2 Indices of solar activity ..................................................................9

2.2.1 Direct solar indices................................................................9

2.2.2 Indirect indices....................................................................13

2.3 Solar activity observations in the pre-telescopic epoch ................................14

2.3.1 Instrumental observations: Camera obscura....................................14

2.3.2 Naked-eye observations ..........................................................14

2.3.3 Mathematical/statistical extrapolations........................................15

2.4 The solar cycle and its variations........................................................16

2.4.1 Quasi-periodicities................................................................16

2.4.2 Randomness vs. regularity........................................................17

2.4.3 A note on solar activity predictions..............................................18

2.5 Summary..................................................................................19

3 The Proxy Method of Past Solar-Activity Reconstruction 21

3.1 The physical basis of the method ........................................................21

3.1.1 Heliospheric modulation of cosmic rays..........................................21

3.1.2 Geomagnetic shielding............................................................23

3.1.3 Cosmic-ray-induced atmospheric cascade ......................................24

3.1.4 Transport and deposition........................................................26

3.2 Radioisotope 14C..........................................................................26

3.2.1 Measurements ....................................................................27

3.2.2 Production........................................................................28

3.2.3 Transport and deposition ............................ 31

3.2.4 The Suess effect and nuclear bomb tests........................................32

3.2.5 The effect of the geomagnetic field..............................................33

3.3 Cosmogenic isotope 10Be..................................................................34

3.3.1 Measurements ....................................................................34

3.3.2 Production........................................................................34

3.3.3 Atmospheric transport .............................. 35

3.3.4 Effect of the geomagnetic field ..................................................36

3.4 Other potential proxy....................................................................38

3.5 Towards a quantitative physical model..................................................38

3.5.1 Regression models................................................................39

3.5.2 Reconstruction of heliospheric parameters......................................40

3.5.3 A link to sunspot numbers ......................................................42

3.6 Solar activity reconstructions .............................. 43

3.7 Verification of reconstructions............................................................44

3.7.1 Comparison with direct data....................................................44

3.7.2 Meteorites and lunar rocks: A direct probe of the galactic cosmic-ray flux . 45

3.7.3 Comparison between isotopes....................................................47

3.8 Composite reconstruction ................................ 47

3.9 Summary ......................................... 49

4 Variability of Solar Activity Over Millennia 50

4.1 Quasi-periodicities and characteristic times ............................................50

4.2 Grand minima of solar activity..........................................................52

4.2.1 The Maunder minimum..........................................................52

4.2.2 Grand minima on a multi-millennial timescale..................................53

4.3 Grand maxima of solar activity..........................................................56

4.3.1 The modern episode of active sun................................................56

4.3.2 Grand maxima on a multi-millennial timescale ................................56

4.4 Related implications......................................................................57

4.4.1 Theoretical constrains............................................................58

4.4.2 Solar-terrestrial relations ........................................................58

4.4.3 Other issues........................................................................59

4.5 Summary..................................................................................59

5 Solar Energetic Particles in the Past 61

5.1 Cosmogenic isotopes......................................................................62

5.2 Lunar and meteoritic rocks ..............................................................65

5.3 Nitrates in polar ice ......................................................................66

5.4 Summary ......................................... 67

6 Conclusions 68

7 Acknowledgements 70 References 71

List of Tables

1 Approximate dates (in -BC/AD) of grand minima in reconstructed solar activity. . 54

2 Approximate dates (in -BC/AD) of grand maxima in the SN-L series................57

3 A list of candidates for extreme SEP events found in different cosmogenic isotope records throughout the Holocene..........................................................63

4 Estimates of 4^ omni-directional integral (above 30 MeV) flux........................66

1 Introduction

The concept of the perfectness and constancy of the sun, postulated by Aristotle, was a strong belief for centuries and an official doctrine of Christian and Muslim countries. However, as people had noticed even before the time of Aristotle, some slight transient changes of the sun can be observed even with the naked eye. Although scientists knew about the existence of "imperfect" spots on the sun since the early 17th century, it was only in the 19th century that the scientific community recognized that solar activity varies in the course of an 11-year solar cycle. Solar variability was later found to have many different manifestations, including the fact that the "solar constant", or the total solar irradiance, TSI, (the amount of total incoming solar electromagnetic radiation in all wavelengths per unit area at the top of the atmosphere) is not a constant. The sun appears much more complicated and active than a static hot plasma ball, with a great variety of nonstationary active processes going beyond the adiabatic equilibrium foreseen in the basic theory of sun-as-star. Such transient nonstationary (often eruptive) processes can be broadly regarded as solar activity, in contrast to the so-called "quiet" sun. Solar activity includes active transient and long-lived phenomena on the solar surface, such as spectacular solar flares, sunspots, prominences, coronal mass ejections (CMEs), etc.

The very fact of the existence of solar activity poses an enigma for solar physics, leading to the development of sophisticated models of an upper layer known as the convection zone and the solar corona. The sun is the only star, which can be studied in great detail and thus can be considered as a proxy for cool stars. Quite a number of dedicated ground-based and space-borne experiments are being carried out to learn more about solar variability. The use of the sun as a paradigm for cool stars leads to a better understanding of the processes driving the broader population of cool sun-like stars. Therefore, studying and modelling solar activity can increase the level of our understanding of nature.

On the other hand, the study of variable solar activity is not of purely academic interest, as it directly affects the terrestrial environment. Although changes in the sun are barely visible without the aid of precise scientific instruments, these changes have great impact on many aspects of our lives. In particular, the heliosphere (a spatial region of about 200-300 astronomical units across) is mainly controlled by the solar magnetic field. This leads to the modulation of galactic cosmic rays (GCRs) by the solar magnetic activity. Additionally, eruptive and transient phenomena in the sun/corona and in the interplanetary medium can lead to sporadic acceleration of energetic particles with greatly enhanced flux. Such processes can modify the radiation environment on Earth and need to be taken into account for planning and maintaining space missions and even transpolar jet flights. Solar activity can cause, through coupling of solar wind and the Earth's magnetosphere, strong geomagnetic storms in the magnetosphere and ionosphere, which may disturb radio-wave propagation and navigation-system stability, or induce dangerous spurious currents in long pipes or power lines. Another important aspect is the link between solar-activity variations and the Earth's climate (see, e.g., reviews by Haigh, 2007; Gray et al., 2010).

It is important to study solar variability on different timescales. The primary basis for such studies is observational (or reconstructed) data. The sun's activity is systematically explored in different ways (solar, heliospheric, interplanetary, magnetospheric, terrestrial), including ground-based and space-borne experiments and dedicated missions during the last few decades, thus covering 3-4 solar cycles. However, it should be noted that the modern epoch is characterized by unusually-high solar activity dominated by an 11-year cyclicity, and it is not straightforward to extrapolate present knowledge (especially empirical and semi-empirical relationships and models) to a longer timescale. The current cycle 24 indicates the return to the normal moderate level of solar activity, as manifested, e.g., via the extended and weak solar minimum in 2008-2009 and weak solar and heliospheric parameters, which are unusual for the space era but may be quite typical for the normal activity (see, e.g., Gibson et al., 2011). Thus, we may experience, in the near

future, the interplanetary conditions quite different with respect to those we got used to during the last decades.

Therefore, the behavior of solar activity in the past, before the era of direct measurements, is of great importance for a variety of reasons. For example, it allows an improved knowledge of the statistical behavior of the solar-dynamo process, which generates the cyclically-varying solar-magnetic field, making it possible to estimate the fractions of time the sun spends in states of very-low activity, what are called grand minima. Such studies require a long time series of solar-activity data. The longest direct series of solar activity is the 400-year-long sunspot-number series, which depicts the dramatic contrast between the (almost spotless) Maunder minimum and the modern period of very high activity. Thanks to the recent development of precise technologies, including accelerator mass spectrometry, solar activity can be reconstructed over multiple millennia from concentrations of cosmogenic isotopes 14C and 10Be in terrestrial archives. This allows one to study the temporal evolution of solar magnetic activity, and thus of the solar dynamo, on much longer timescales than are available from direct measurements.

This paper gives an overview of the present status of our knowledge of long-term solar activity, covering the period of Holocene (the last 11 millennia). A description of the concept of solar activity and a discussion of observational methods and indices are presented in Section 2. The proxy method of solar-activity reconstruction is described in some detail in Section 3. Section 4 gives an overview of what is known about past solar activity. The long-term averaged flux of solar energetic particles is discussed in Section 5. Finally, conclusions are summarized in Section 6.

2 Solar Activity: Concept and Observations

2.1 The concept of solar activity

The sun is known to be far from a static state, the so-called "quiet" sun described by simple stellar-evolution theories, but instead goes through various nonstationary active processes. Such nonstationary and nonequilibrium (often eruptive) processes can be broadly regarded as solar activity. Whereas the concept of solar activity is quite a common term nowadays, it is neither straightforwardly interpreted nor unambiguously defined. For instance, solar-surface magnetic variability, eruption phenomena, coronal activity, radiation of the sun as a star or even interplanetary transients and geomagnetic disturbances can be related to the concept of solar activity. A variety of indices quantifying solar activity have been proposed in order to represent different observables and caused effects. Most of the indices are highly correlated to each other due to the dominant 11-year cycle, but may differ in fine details and/or long-term trends. In addition to the solar indices, indirect proxy data is often used to quantify solar activity via its presumably known effect on the magnetosphere or heliosphere. The indices of solar activity that are often used for long-term studies are reviewed below.

2.2 Indices of solar activity

Solar (as well as other) indices can be divided into physical and synthetic according to the way they are obtained/calculated. Physical indices quantify the directly-measurable values of a real physical observable, such as, e.g., the radioflux, and thus have clear physical meaning as they quantify physical features of different aspects of solar activity and their effects. Synthetic indices (the most common being sunspot number) are calculated (or synthesized) using a special algorithm from observed (often not measurable in physical units) data or phenomena. Additionally, solar activity indices can be either direct (i.e., directly relating to the sun) or indirect (relating to indirect effects caused by solar activity), as discussed in subsequent Sections 2.2.1 and 2.2.2.

2.2.1 Direct solar indices

The most commonly used index of solar activity is based on sunspot number. Sunspots are dark areas on the solar disc (of size up to tens of thousands of km, lifetime up to half-a-year), characterized by a strong magnetic field, which leads to a lower temperature (about 4000 K compared to 5800 K in the photosphere) and observed as darkening.

Sunspot number is a synthetic, rather than a physical, index, but it has still become quite a useful parameter in quantifying the level of solar activity. This index presents the weighted number of individual sunspots and/or sunspot groups, calculated in a prescribed manner from simple visual solar observations. The use of the sunspot number makes it possible to combine together thousands and thousands of regular and fragmentary solar observations made by earlier professional and amateur astronomers. The technique, initially developed by Rudolf Wolf, yielded the longest series of directly and regularly-observed scientific quantities. Therefore, it is common to quantify solar magnetic activity via sunspot numbers. For details see the review on sunspot numbers and solar cycles (Hathaway and Wilson, 2004; Hathaway, 2010).

Wolf sunspot number (WSN) series

The concept of the sunspot number was developed by Rudolf Wolf of the Zürich observatory in the middle of the 19th century. The sunspot series, initiated by him, is called the Zürich or Wolf sunspot number (WSN) series. The relative sunspot number Rz is defined as

Rz = k (10 G + N), (1)

where G is the number of sunspot groups, N is the number of individual sunspots in all groups visible on the solar disc and k denotes the individual correction factor, which compensates for differences in observational techniques and instruments used by different observers, and is used to normalize different observations to each other.

Figure 1: Sunspot numbers since 1610. a) Monthly (since 1749) and yearly (1700-1749) Wolf sunspot number series. b) Monthly group sunspot number series. The grey line presents the 11-year running mean after the Maunder minimum. Standard (Zurich) cycle numbering as well as the Maunder minimum (MM) and Dalton minimum (DM) are shown in the lower panel.

The value of Rz (see Figure 1a) is calculated for each day using only one observation made by the "primary" observer (judged as the most reliable observer during a given time) for the day. The primary observers were Staudacher (1749-1787), Flaugergues (1788-1825), Schwabe (18261847), Wolf (1848-1893), Wolfer (1893-1928), Brunner (1929-1944), Waldmeier (1945-1980) and Koeckelenbergh (since 1980). If observations by the primary observer are not available for a certain day, the secondary, tertiary, etc. observers are used (see the hierarchy of observers in Waldmeier, 1961). The use of only one observer for each day aims to make Rz a homogeneous time series. As a drawback, such an approach ignores all other observations available for the day, which constitute a large fraction of the existing information. Moreover, possible errors of the primary observer cannot be caught or estimated. The observational uncertainties in the monthly Rz can be up to 25% (e.g., Vitinsky et al., 1986). The WSN series is based on observations performed at the Zurich Observatory during 1849-1981 using almost the same technique. This part of the series is fairly stable and homogeneous although an offset due to the change of the weighting procedure might have been introduced in 1945-1946 (Svalgaard, 2012). However, prior to that there have been many gaps in the data that were interpolated. If no sunspot observations are available for some period, the data gap is filled, without note in the final WSN series, using an interpolation between the available data and by employing some proxy data. In addition, earlier parts of the sunspot series were "corrected" by Wolf using geomagnetic observation (see details in Svalgaard, 2012), which makes the series less homogeneous. Therefore, the WSN series is a combination of

direct observations and interpolations for the period before 1849, leading to possible errors and inhomogeneities as discussed, e.g, by Vitinsky et al. (1986); Wilson (1998); Letfus (1999); Svalgaard (2012). The quality of the Wolf series before 1749 is rather poor and hardly reliable (Hoyt et al., 1994; Hoyt and Schatten, 1998; Hathaway and Wilson, 2004).

Note that the sun has been routinely photographed since 1876 so that full information on daily sunspot activity is available (the Greenwich series) with observational uncertainties being negligible for the last 140 years.

The routine production of the WSN series was terminated in Zurich in 1982. Since then, the sunspot number series is routinely updated as the International sunspot number Ri, provided by the Solar Influences Data Analysis Center in Belgium (Clette et al., 2007). The international sunspot number series is computed using the same definition (Equation 1) as WSN but it has a significant distinction from the WSN; it is based not on a single primary solar observation for each day but instead uses a weighted average of more than 20 approved observers.

In addition to the standard sunspot number Ri, there is also a series of hemispheric sunspot numbers RN and RS, which account for spots only in the northern and southern solar hemispheres, respectively (note that Ri = RN + RS). These series are used to study the N-S asymmetry of solar activity (Temmer et al., 2002).

Group sunspot number (GSN) series

Since the WSN series is of lower quality before the 1850s and is hardly reliable before 1750, there was a need to re-evaluate early sunspot data. This tremendous work has been done by Hoyt and Schatten (1996, 1998), who performed an extensive archive search and nearly doubled the amount of original information compared to the Wolf series. They have produced a new series of sunspot activity called the group sunspot numbers (GSN - see Figure 1b), including all available archival records. The daily group sunspot number Rg is defined as follows:

r9=12T £H , (2)

where Gj is the number of sunspot groups recorded by the j-th observer, k' is the observer's individual correction factor, n is the number of observers for the particular day, and 12.08 is a normalization number scaling Rg to Rz values for the period of 1874-1976. Rg is more robust than Rz since it is based on more easily identified sunspot groups and does not include the number of individual spots. The GSN series includes not only one "primary" observation, but all available observations, and covers the period since 1610, being, thus, 140 years longer than the original WSN series. It is particularly interesting that the period of the Maunder minimum (1645-1715) was surprisingly well covered with daily observations (Ribes and Nesme-Ribes, 1993; Hoyt and Schatten, 1996) allowing for a detailed analysis of sunspot activity during this grand minimum (see also Section 4.2). Systematic uncertainties of the Rg values are estimated to be about 10% before 1640, less than 5% from 1640-1728 and from 1800-1849, 15-20% from 1728-1799, and about 1% since 1849 (Hoyt and Schatten, 1998). The GSN series is more reliable and homogeneous than the WSN series before 1849. The two series are nearly identical after the 1870s (Hoyt and Schatten, 1998; Letfus, 1999; Hathaway and Wilson, 2004). However, the GSN series still contains some lacunas, uncertainties and possible inhomogeneities (see, e.g., Letfus, 2000; Usoskin et al., 2003a; Vaquero et al., 2012).

The search for other lost or missing records of past solar instrumental observations has not ended even since the extensive work by Hoyt and Schatten. Archival searches still give new interesting findings of forgotten sunspot observations, often outside major observatories - see a detailed review book by Vaquero and Vazquez (2009) and original papers by Casas et al. (2006); Vaquero et al. (2005, 2007); Arlt (2008, 2009). Interestingly, not only sunspot counts but also regular drawings, forgotten for centuries, are being restored nowadays in dusty archives. A very interesting work has

been done by Rainer Arlt (Arlt, 2008, 2009; Arlt and Abdolvand, 2011; Arlt, 2013) on recovering, digitizing, and analyzing regular drawings by S.H. Schwabe of 1825-1867 and J.C. Staudacher of 1749-1796. This work led to the extension of the Maunder butterfly diagram for several solar cycles backwards (Arlt, 2009; Usoskin et al., 2009c; Arlt and Abdolvand, 2011; Arlt, 2013) - see a newly built diagram for solar cycles Nos. 7-10 shown in Figure 2. In particular, this data confirms that GSN series is more homogenous before 1874 that WSN. A recent finding of the lost data by G. Marcgraf and correcting some earlier uncertain data for the period 1636-1642 by Vaquero et al. (2011) made it possible to revise the pattern of the beginning of the Maunder minimum.

1826 1828 1830 1832 1834 1836 1838 1840 1842 1844 1846 1848 1850 1852 1854 1856 1858 1860 1862 1864 1866

Figure 2: Maunder butterfly diagram of sunspot occurrence reconstructed by Arlt (2013) for 1825-1867 using recovered drawing of S.H. Schwabe.

Other indices

An example of a synthetic index of solar activity is the flare index, representing solar flare activity (e.g., Ozgiic et al., 2003; Kleczek, 1952). The flare index quantifies daily flare activity in the following manner; it is computed as a product of the flare's relative importance I in the Ha-range and duration t, Q = It, thus being a rough measure of the total energy emitted by the flare. The daily flare index is produced by Bogazici University (Ozguc et al., 2003) and is available since 1936.

A traditional physical index of solar activity is related to the radioflux of the sun in the wavelength range of 10.7 cm and is called the F10.7 index (e.g., Tapping and Charrois, 1994). This index represents the flux (in solar flux units, 1 sfu = 10-22 Wm-2 Hz-1) of solar radio emission at a centimetric wavelength. There are at least two sources of 10.7 cm flux - free-free emission from hot coronal plasma and gyromagnetic emission from active regions (Tapping, 1987). It is a good quantitative measure of the level of solar activity, which is not directly related to sunspots. Close correlation between the F10.7 index and sunspot number indicates that the latter is a good index of general solar activity, including coronal activity. The solar F10.7 cm record has been measured continuously since 1947.

Another physical index is the coronal index (e.g., Rybansky et al., 2005), which is a measure of the irradiance of the sun as a star in the coronal green line. Computation of the coronal index is based on observations of green corona intensities (Fe XIV emission line at 530.3 nm wavelength) from coronal stations all over the world, the data being transformed to the Lomnicky Stit photometric scale. This index is considered a basic optical index of solar activity. A synthesized

homogeneous database of the Fe XIV 530.3 nm coronal-emission line intensities has existed since 1943 and covers seven solar cycles.

Often sunspot area is considered as a physical index representing solar activity (e.g., Baranyi et al., 2001; Balmaceda et al., 2005). This index gives the total area of visible spots on the solar disc in units of millionths of the sun's visible hemisphere, corrected for apparent distortion due to the curvature of the solar surface. The area of individual groups may vary between tens of millionths (for small groups) up to several thousands of millionths for huge groups. This index has a physical meaning related to the solar magnetic flux emerging at sunspots. Sunspot areas are available since 1874 in the Greenwich series obtained from daily photographic images of the sun. In addition, some fragmentary data of sunspot areas, obtained from solar drawings, are available for earlier periods (Vaquero et al., 2004; Arlt, 2008).

An important quantity is solar irradiance, total and spectral (Frohlich, 2012). Irradiance variations are physically related to solar magnetic variability (e.g., Solanki et al., 2000), and are often considered manifestations of solar activity, which is of primary importance for solar-terrestrial relations.

Other physical indices include spectral sun-as-star observations, such as the Ca II-K index (e.g., Donnelly et al., 1994; Foukal, 1996), the space-based Mg II core-to-wing ratio as an index of solar UVI (e.g., Donnelly et al., 1994; Viereck and Puga, 1999; Snow et al., 2005) and many others.

All the above indices are closely correlated to sunspot numbers on the solar-cycle scale, but may depict quite different behavior on short or long timescales.

2.2.2 Indirect indices

Sometimes quantitative measures of solar-variability effects are also considered as indices of solar activity. These are related not to solar activity per se, but rather to its effect on different environments. Accordingly, such indices are called indirect, and can be roughly divided into terrestrial/geomagnetic and heliospheric/interplanetary.

Geomagnetic indices quantify different effects of geomagnetic activity ultimately caused by solar variability, mostly by variations of solar-wind properties and the interplanetary magnetic field. For example, the aa-index, which provides a global index of magnetic activity relative to a quiet-day curve for a pair of antipodal magnetic observatories (in England and Australia), is available from 1868 (Mayaud, 1972). An extension of the geomagnetic series is available from the 1840s using the Helsinki Ak(H) index (Nevanlinna, 2004a,b). Although the homogeneity of the geomagnetic series is compromised (e.g., Lukianova et al., 2009; Love, 2011), it still remains an important indirect index of solar activity. A review of the geomagnetic effects of solar activity can be found, e.g., in Pulkkinen (2007). It is noteworthy that geomagnetic indices, in particular low-latitude aurorae (Silverman, 2006), are associated with coronal/interplanetary activity (high-speed solarwind streams, interplanetary transients, etc.) that may not be directly related to the sunspot-cycle phase and amplitude, and therefore serve only as an approximate index of solar activity. One of the earliest instrumental geomagnetic indices is related to the daily magnetic declination range, the range of diurnal variation of magnetic needle readings at a fixed location, and is available from the 1780s (Nevanlinna, 1995). However, this data exists as several fragmentary sets, which are difficult to combine into a homogeneous data series.

Heliospheric indices are related to features of the solar wind or the interplanetary magnetic field measured (or estimated) in the interplanetary space. For example, the time evolution of the total (or open) solar magnetic flux is extensively debated (e.g., Lockwood et al., 1999; Wang et al., 2005; Krivova et al., 2007).

A special case of heliospheric indices is related to the galactic cosmic-ray intensity recorded in natural terrestrial archives. Since this indirect proxy is based on data recorded naturally throughout the ages and revealed now, it makes possible the reconstruction of solar activity changes on

long timescales, as discussed in Section 3.

2.3 Solar activity observations in the pre-telescopic epoch

Instrumental solar data is based on regular observation (drawings or counting of spots) of the sun using optical instruments, e.g., the telescope used by Galileo in the early 17th century. These observations have mostly been made by professional astronomers whose qualifications and scientific thoroughness were doubtless. They form the basis of the Group sunspot-number series (Hoyt and Schatten, 1998), which can be more-or-less reliably extended back to 1610 (see discussion in Section 2.2.1). However, some fragmentary records of qualitative solar and geomagnetic observations exist even for earlier times, as discussed below (Sections 2.3.1-2.3.2).

2.3.1 Instrumental observations: Camera obscura

The invention of the telescope revolutionized astronomy. However, another solar astronomical instrument, the camera obscura, also made it possible to provide relatively good solar images and was still in use until the late 18th century. Camera obscuras were known from early times, and they have been used in major cathedrals to define the sun's position (see the review by Vaquero, 2007; Vaquero and Vázquez, 2009). The earliest known drawing of the solar disc was made by Frisius, who observed the solar eclipse in 1544 using a camera obscura. That observation was performed during the Sporer minimum and no spots were observed on the sun. The first known observation of a sunspot using a camera obscura was done by Kepler in May 1607, who erroneously ascribed the spot on the sun to a transit of Mercury. Although such observations were sparse and related to other phenomena (solar eclipses or transits of planets), there were also regular solar observations by camera obscura. For example, about 300 pages of logs of solar observations made in the cathedral of San Petronio in Bologna from 1655-1736 were published by Eustachio Manfredi in 1736 (see the full story in Vaquero, 2007).

Therefore, observations and drawings made using camera obscura can be regarded as instrumental observations.

2.3.2 Naked-eye observations

Even before regular professional observations performed with the aid of specially-developed instruments (what we now regard as scientific observations) people were interested in unusual phenomena. Several historical records exist based on naked-eye observations of transient phenomena on the sun or in the sky.

From even before the telescopic era, a large amount of evidence of spots being observed on the solar disc can be traced back as far as to the middle of the 4th century BC (Theophrastus of Athens). The earliest known drawing of sunspots is dated to December 8, 1128 AD as published in "The Chronicle of John of Worcester" (Willis and Stephenson, 2001). However, such evidence from occidental and Moslem sources is scarce and mostly related to observations of transits of inner planets over the sun's disc, probably because of the dominance of the dogma on the perfectness of the sun's body, which dates back to Aristotle's doctrine (Bray and Loughhead, 1964). Oriental sources are much richer for naked-eye sunspot records, but that data is also fragmentary and irregular (see, e.g., Clark and Stephenson, 1978; Wittmann and Xu, 1987; Yau and Stephenson, 1988). Spots on the sun are mentioned in official Chinese and Korean chronicles from 165 BC to 1918 AD. While these chronicles are fairly reliable, the data is not straightforward to interpret since it can be influenced by meteorological phenomena, e.g., dust loading in the atmosphere due to dust storms (Willis et al., 1980) or volcanic eruptions (Scuderi, 1990) can facilitate sunspots observations. Direct comparison of Oriental naked-eye sunspot observations and European telescopic data shows that naked-eye observations can serve only as a qualitative indicator of sunspot activity,

but can hardly be quantitatively interpreted (see, e.g., Willis et al., 1996, and references therein). Moreover, as a modern experiment of naked-eye observations (Mossman, 1989) shows, Oriental chronicles contain only a tiny (1/2oo - Viooo) fraction of the number of sunspots potentially visible with the naked eye (Eddy et al., 1989). This indicates that records of sunspot observations in the official chronicles were highly irregular (Eddy, 1983) and probably dependent on dominating traditions during specific historical periods (Clark and Stephenson, 1978). Although naked-eye observations tend to qualitatively follow the general trend in solar activity according to a posteriori information (e.g., Vaquero et al., 2002), extraction of any independent quantitative information from these records seems impossible.

Visual observations of aurorae borealis at middle latitudes form another proxy for solar activity (e.g., Siscoe, 1980; Schove, 1983; Krivsky, 1984; Silverman, 1992; Schroder, 1992; Lee et al, 2004; Basurah, 2004; Vázquez and Vaquero, 2010). Fragmentary records of aurorae can be found in both occidental and oriental sources since antiquity. The first known dated notation of an aurora is from March 12, 567 BC from Babylon (Stephenson et al., 2004). Aurorae may appear at middle latitudes as a result of enhanced geomagnetic activity due to transient interplanetary phenomena. Although auroral activity reflects coronal and interplanetary features rather than magnetic fields on the solar surface, there is a strong correlation between long-term variations of sunspot numbers and the frequency of aurora occurrences. Because of the phenomenon's short duration and low brightness, the probability of seeing aurora is severely affected by other factors such as the weather (sky overcast, heat lightnings), the Moon's phase, season, etc. The fact that these observations were not systematic in early times (before the beginning of the 18th century) makes it difficult to produce a homogeneous data set. Moreover, the geomagnetic latitude of the same geographical location may change quite dramatically over centuries, due to the migration of the geomagnetic axis, which also affects the probability of watching aurorae (Siscoe and Verosub, 1983; Oguti and Egeland, 1995). For example, the geomagnetic latitude of Seoul (37.5° N 127° E), which is currently less than 30°, was about 40° a millennium ago (Kovaltsov and Usoskin, 2007). This dramatic change alone can explain the enhanced frequency of aurorae observations recorded in oriental chronicles.

2.3.3 Mathematical/statistical extrapolations

Due to the lack of reliable information regarding solar activity in the pre-instrumental era, it seems natural to try to extend the sunspot series back in time, before 1610 AD, by means of extrapolating its statistical properties. Indeed, numerous attempts of this kind have been made even recently (e.g., Nagovitsyn, 1997; de Meyer, 1998; Rigozo et al., 2001). Such models aim to find the main feature of the actually-observed sunspot series, e.g., a modulated carrier frequency or a multi-harmonic representation, which is then extrapolated backwards in time. The main disadvantage of this approach is that it is not a reconstruction based upon measured or observed quantities, but rather a "post-diction" based on extrapolation. This method is often used for short-term predictions, but it can hardly be used for the reliable long-term reconstruction of solar activity. In particular, it assumes that the sunspot time series is stationary, i.e., a limited time realization contains full information on its future and past. Clearly such models cannot include periods exceeding the time span of observations upon which the extrapolation is based. Hence, the pre- or post-diction becomes increasingly unreliable with growing extrapolation time and its accuracy is hard to estimate.

Sometimes a combination of the above approaches is used, i.e., a fit of the mathematical model to indirect qualitative proxy data. In such models a mathematical extrapolation of the sunspot series is slightly tuned and fitted to some proxy data for earlier times. For example, Schove (1955, 1979) fitted the slightly variable but phase-locked carrier frequency (about 11 years) to fragmentary data from naked-eye sunspot observations and auroral sightings. The phase locking is achieved by assuming exactly nine solar cycles per calendar century. This series, known as Schove

series, reflects qualitative long-term variations of the solar activity, including some grand minima, but cannot pretend to be a quantitative representation in solar activity level. The Schove series played an important historical role in the 1960s. In particular, a comparison of the A14C data with this series succeeded in convincing the scientific community that secular variations of 14C in tree rings have solar and not climatic origins (Stuiver, 1961). This formed a cornerstone of the precise method of solar-activity reconstruction, which uses cosmogenic isotopes from terrestrial archives. However, attempts to reconstruct the phase and amplitude of the 11-year cycle, using this method, were unsuccessful. For example, Schove (1955) made predictions of forthcoming solar cycles up to 2005, which failed. We note that all these works are not able to reproduce, for example, the Maunder minimum (which cannot be represented as a result of the superposition of different harmonic oscillations), yielding too high sunspot activity compared to that observed. From the modern point of view, the Schove series can be regarded as archaic, but it is still in use in some studies.

2.4 The solar cycle and its variations 2.4.1 Quasi-periodicities

The main feature of solar activity is its pronounced quasi-periodicity with a period of about 11 years, known as the Schwabe cycle. However, the cycle varies in both amplitude and duration. The first observation of a possible regular variability in sunspot numbers was made by the Danish astronomer Christian Horrebow in the 1770s on the basis of his sunspot observations from 1761 — 1769 (see details in Gleissberg, 1952; Vitinsky, 1965), but the results were forgotten. It took over 70 years before the amateur astronomer Schwabe announced in 1844 that sunspot activity varies cyclically with a period of about 10 years. This cycle, called the 11-year or Schwabe cycle, is the most prominent variability in the sunspot-number series. It is recognized now as a fundamental feature of solar activity originating from the solar-dynamo process. This 11-year cyclicity is prominent in many other parameters including solar, heliospheric, geomagnetic, space weather, climate and others. The background for the 11-year Schwabe cycle is the 22-year Hale magnetic polarity cycle. Hale found that the polarity of sunspot magnetic fields changes in both hemispheres when a new 11-year cycle starts (Hale et al., 1919). This relates to the reversal of the global magnetic field of the sun with the period of 22 years. It is often considered that the 11-year Schwabe cycle is the modulo of the sign-alternating Hale cycle (e.g., Sonett, 1983; Bracewell, 1986; Kurths and Ruz-maikin, 1990; de Meyer, 1998; Mininni et al., 2001), but this is only a mathematical representation. A detailed review of solar cyclic variability can be found in (Hathaway, 2010).

Sometimes the regular time evolution of solar activity is broken up by periods of greatly depressed activity called grand minima. The last grand minimum (and the only one covered by direct solar observations) was the famous Maunder minimum from 1645-1715 (Eddy, 1976, 1983). Other grand minima in the past, known from cosmogenic isotope data, include, e.g., the Sporer minimum around 1450-1550 and the Wolf minimum around the 14th century (see the detailed discussion in Section 4.2). Sometimes the Dalton minimum (ca. 1790-1820) is also considered to be a grand minimum. However, sunspot activity was not completely suppressed and still showed Schwabe cyclicity during the Dalton minimum. As suggested by Schiissler et al. (1997), this can be a separate, intermediate state of the dynamo between the grand minimum and normal activity, or an unsuccessful attempt of the sun to switch to the grand minimum state (Frick et al., 1997; Sokoloff, 2004). This is observed as the phase catastrophe of solar-activity evolution (e.g., Vitinsky et al., 1986; Kremliovsky, 1994). A peculiarity in the phase evolution of sunspot activity around 1800 was also noted by Sonett (1983), who ascribed it to a possible error in Wolf sunspot data and by Wilson (1988a), who reported on a possible misplacement of sunspot minima for cycles 4-6 in the WSN series. It has been also suggested that the phase catastrophe can be related to a tiny cycle, which might have been lost at the end of the 18th century because of very sparse

observations (Usoskin et al., 2001b, 2002a, 2003b; Zolotova and Ponyavin, 2007). We note that a new independent evidence proving the existence of the lost cycle has been found recently in the reconstructed sunspot butterfly diagram for that period (Usoskin et al., 2009c).

The long-term change (trend) in the Schwabe cycle amplitude is known as the secular Gleissberg cycle (Gleissberg, 1939) with the mean period of about 90 years. However, the Gleissberg cycle is not a cycle in the strict periodic sense but rather a modulation of the cycle envelope with a varying timescale of 60-120 years (e.g., Gleissberg, 1971; Kuklin, 1976; Ogurtsov et al., 2002).

Longer (super-secular) cycles cannot be studied using direct solar observations, but only indica-tively by means of indirect proxies such as cosmogenic isotopes discussed in Section 3. Analysis of the proxy data also yields the Gleissberg secular cycle (Feynman and Gabriel, 1990; Peristykh and Damon, 2003), but the question of its phase locking and persistency/intermittency still remains open. Several longer cycles have been found in the cosmogenic isotope data. A cycle with a period of 205-210 years, called the de Vries or Suess cycle in different sources, is a prominent feature, observed in various cosmogenic data (e.g., Suess, 1980; Sonett and Finney, 1990; Zhentao, 1990; Usoskin et al., 2004). Sometimes variations with a characteristic time of 600-700 years or 1000-1200 years are discussed (e.g., Vitinsky et al., 1986; Sonett and Finney, 1990; Vasiliev and Dergachev, 2002; Steinhilber et al., 2012; Abreu et al., 2012), but they are intermittent and can hardly be regarded as a typical feature of solar activity. A 2000 -2400-year cycle is also noticeable in radiocarbon data series (see, e.g., Vitinsky et al., 1986; Damon and Sonett, 1991; Vasiliev and Dergachev, 2002). However, the non-solar origin of these super-secular cycles (e.g., geomagnetic or climatic variability) cannot be excluded.

2.4.2 Randomness vs. regularity

The short-term (days - months) variability of sunspot numbers is greater than the observational uncertainties indicating the presence of random fluctuations (noise). As typical for most real signals, this noise is not uniform (white), but rather red or correlated noise (e.g., Ostryakov and Usoskin, 1990; Oliver and Ballester, 1996; Frick et al., 1997), namely, its variance depends on the level of the signal. While the existence of regularity and randomness in sunspot series is apparent, their relationship is not clear (e.g., Wilson, 1994) - are they mutually independent or intrinsically tied together? Moreover, the question of whether randomness in sunspot data is due to chaotic or stochastic processes is still open.

Earlier it was common to describe sunspot activity as a multi-harmonic process with several basic harmonics (e.g., Vitinsky, 1965; Sonett, 1983; Vitinsky et al., 1986) with an addition of random noise, which plays no role in the solar-cycle evolution. However, it has been shown (e.g., Rozelot, 1994; Weiss and Tobias, 2000; Charbonneau, 2001; Mininni et al., 2002) that such an oversimplified approach depends on the chosen reference time interval and does not adequately describe the long-term evolution of solar activity. A multi-harmonic representation is based on an assumption of the stationarity of the benchmark series, but this assumption is broadly invalid for solar activity (e.g., Kremliovsky, 1994; Sello, 2000; Polygiannakis et al., 2003). Moreover, a multi-harmonic representation cannot, for an apparent reason, be extrapolated to a timescale larger than that covered by the benchmark series. The fact that purely mathematical/statistical models cannot give good predictions of solar activity (as will be discussed later) implies that the nature of the solar cycle is not a multi-periodic or other purely deterministic process, but random (chaotic or stochastic) processes play an essential role in sunspot cycle formation (e.g., Moss et al., 2008; Kapyla et al., 2012). An old idea of the possible planetary influence on the dynamo has received a new pulse recently with some unspecified torque effect on the assumed quasi-rigid non-axisymmetric tahocline (Abreu et al., 2012). If confirmed this idea would imply a significant multi-harmonic driver of the solar activity, but the question is still open. Different numeric tests, such as an analysis of the Lyapunov exponents (Ostriakov and Usoskin, 1990; Mundt et al., 1991;

Kremliovsky, 1995; Sello, 2000), Kolmogorov entropy (Carbonell et al., 1994; Sello, 2000) and Hurst exponent (Ruzmaikin et al., 1994; Oliver and Ballester, 1998), confirm the chaotic/stochastic nature of the solar-activity time evolution (see, e.g., the recent review by Panchev and Tsekov, 2007).

It was suggested quite a while ago that the variability of the solar cycle may be a temporal realization of a low-dimensional chaotic system (e.g., Ruzmaikin, 1981). This concept became popular in the early 1990s, when many authors considered solar activity as an example of low-dimensional deterministic chaos, described by the strange attractor (e.g., Kurths and Ruzmaikin, 1990; Ostri-akov and Usoskin, 1990; Morfill et al., 1991; Mundt et al., 1991; Rozelot, 1995; Salakhutdinova, 1999; Serre and Nesme-Ribes, 2000; Hanslmeier et al., 2013). Such a process naturally contains randomness, which is an intrinsic feature of the system rather than an independent additive or multiplicative noise. However, although this approach easily produces features seemingly similar to those of solar activity, quantitative parameters of the low-dimensional attractor have varied greatly as obtained by different authors. Later it was realized that the analyzed data set was too short (Carbonell et al., 1993, 1994), and the results were strongly dependent on the choice of filtering methods (Price et al., 1992). Developing this approach, Mininni et al. (2000, 2001) suggest that one consider sunspot activity as an example of a 2D Van der Pol relaxation oscillator with an intrinsic stochastic component.

Such phenomenological or basic principles models, while succeeding in reproducing (to some extent) the observed features of solar-activity variability, do not provide insight into the nature of regular and random components of solar variability. In this sense efforts to understand the nature of randomness in sunspot activity in the framework of dynamo theory are more advanced. Corresponding theoretical dynamo models have been developed (see reviews by Ossendrijver, 2003; Charbonneau, 2010), which include stochastic processes (e.g., Weiss et al., 1984; Feynman and Gabriel, 1990; Schmalz and Stix, 1991; Moss et al., 1992; Hoyng, 1993; Brooke and Moss, 1994; Lawrence et al., 1995; Schmitt et al., 1996; Charbonneau and Dikpati, 2000; Brandenburg and Sokoloff, 2002). For example, Feynman and Gabriel (1990) suggest that the transition from a regular to a chaotic dynamo passes through bifurcation. Charbonneau and Dikpati (2000) studied stochastic fluctuations in a Babcock-Leighton dynamo model and succeeded in the qualitative reproduction of the anti-correlation between cycle amplitude and length (Waldmeier rule). Their model also predicts a phase-lock of the Schwabe cycle, i.e., that the 11-year cycle is an internal "clock" of the sun. Most often the idea of fluctuations is related to the a-effect, which is the result of the electromotive force averaged over turbulent vortices, and thus can contain a fluctuating contribution (e.g., Hoyng, 1993; Ossendrijver et al., 1996; Brandenburg and Spiegel, 2008; Moss et al., 2008). Note that a significant fluctuating component (with the amplitude more than 100% of the regular component) is essential in all these model.

2.4.3 A note on solar activity predictions

Randomness (see Section 2.4.2) in the SN series is directly related to the predictability of solar activity. Forecasting solar activity has been a subject of intense study for many years (e.g., Yule, 1927; Newton, 1928; Gleissberg, 1948; Vitinsky, 1965) and has greatly intensified recently with a hundred of journal articles being published to predict the solar cycle No. 24 maximum (see, e.g., the review by Pesnell, 2012), following the boost of space-technology development and increasing debates on solar-terrestrial relations. In fact, the situation has not been improved since the previous cycle, No. 23. The predictions for the peak sunspot number of solar cycle No. 24 range by a factor of 5, between 40 and 200, reflecting the lack of a reliable consensus method (Tobias et al., 2006). Detailed review of the solar activity prediction methods and results have been recently provided by (Hathaway, 2009; Petrovay, 2010; Pesnell, 2012).

A detailed classification of the prediction methods is given by Pesnell (2012) who separates climatology, precursor, theoretical (dynamo model), spectral, neural network, and stock market

prediction methods. All prediction methods can be generically divided into precursor and statistical (including the majority of the above classifications) techniques or their combinations (Hathaway et al., 1999). The fact that the prediction of solar cycle is not improved with adding more data (the new solar cycle) suggests that such methods are not able to give reliable prognoses.

The precursor methods are usually based on phenomenological, but sometimes physical, links between the poloidal solar-magnetic field, estimated, e.g., from geomagnetic activity in the declining phase of the preceding cycle or in the minimum time (e.g., Hathaway, 2009), with the toroidal field responsible for sunspot formation. These methods usually yield better short-term predictions of a forthcoming cycle maximum than the statistical methods, but cannot be applied to timescales longer than one solar cycle.

Statistical methods, including a low-dimensional solar-attractor representation (Kurths and Ruzmaikin, 1990), are based solely on the statistical properties of sunspot activity and may give a reasonable result for short-term forecasting, but yield very poor results for long-term predictions (see reviews by, e.g., Conway, 1998; Hathaway et al., 1999; Li et al., 2001; Usoskin and Mursula, 2003; Kane, 2007) because of chaotic/stochastic behavior (see Section 2.4.2).

A new method based on sophisticated dynamo numerical simulations emerges (e.g., Dikpati and Gilman, 2006; Dikpati et al., 2008; Choudhuri et al., 2007; Jiang et al., 2007), but the results are contradictory with each other. Prospectives of this approach are also not clear because of the stochastic component, which drives the dynamo out of the deterministic regime, and uncertainties in the input parameters (Tobias et al., 2006; Bushby and Tobias, 2007; Karak and Nandy, 2012).

Some models, mostly based on precursor method, succeed in reasonable predictions of a forthcoming solar cycle (i.e., several years ahead), but they do not pretend to extend further in time. On the other hand, many claims of the solar activity forecast for 40-50 years ahead and even beyond have been made recently, often without sensible argumentation. However, so far there is no evidence of any method giving a reasonable prediction of solar activity beyond the solar-cycle scale (see, e.g., Section 2.3.3), probably because of the intrinsic limit of solar-activity predictability due to its stochastic/chaotic nature (Kremliovsky, 1995; Tobias et al., 2006). Accordingly, such attempts can be regarded as speculative, unless they are verified by the actual behavior of solar activity. Note that even an exact prediction of the amplitude of one solar cycle can be just a random coincidence and cannot serve as a proof of the method's veracity. Only a sequence of successful predictions can form a basis for confidence, which requires several decades.

Note that several "predictions" of the general decline of the coming solar activity have been made recently (Solanki et al., 2004; Abreu et al., 2008; Lockwood et al., 2011), however, these are not really true predictions but rather the acknowledge of the fact that the Modern Grand maximum (Usoskin et al., 2003c; Solanki et al., 2004) must cease. Similar caution can be made about predictions of a Grand minimum (e.g., Lockwood et al., 2011; Miyahara et al., 2010) -a grand minimum should appear soon or later, but presently we are hardly able to predict its occurrence.

2.5 Summary

In this section, the concept of solar activity and quantifying indices is discussed, as well as the main features of solar-activity temporal behavior.

The concept of solar activity is quite broad and covers non-stationary and non-equilibrium (often eruptive) processes, in contrast to the "quiet" sun concept, and their effects upon the terrestrial and heliospheric environment. Many indices are used to quantify different aspects of variable solar activity. Quantitative indices include direct (i.e., related directly to solar variability) and indirect (i.e., related to terrestrial and interplanetary effects caused by solar activity), they can be physical or synthetic. While all indices depict the dominant 11-year cyclic variability, their relationships on other timescales (short scale or long-term trends) may vary to a great extent.

The most common and the longest available index of solar activity is the sunspot number, which is a synthetic index and is very useful for the quantitative representation of overall solar activity outside the grand minimum. During the grand Maunder minimum, however, it may give only a clue about solar activity whose level may drop below the sunspot formation threshold. The sunspot number series is available for the period from 1610 AD, after the invention of the telescope, and covers, in particular, the Maunder minimum in the late 17th century. Fragmentary non-instrumental observations of the sun before 1610, while giving a possible hint of relative changes in solar activity, cannot be interpreted in a quantitative manner.

Solar activity in all its manifestations is dominated by the 11-year Schwabe cycle, which has, in fact, a variable length of 9-14 years for individual cycles. The amplitude of the Schwabe cycle varies greatly - from the almost spotless Maunder minimum to the very high cycle 19, possibly in relation to the Gleissberg or secular cycle. Longer super-secular characteristic times can also be found in various proxies of solar activity, as discussed in Section 4.

Solar activity contains essential chaotic/stochastic components, that lead to irregular variations and make the prediction of solar activity for a timescale exceeding one solar cycle impossible.

3 The Proxy Method of Past Solar-Activity Reconstruction

In addition to direct solar observations, described in Section 2.2.1, there are also indirect solar proxies, which are used to study solar activity in the pre-telescopic era. Unfortunately, we do not have any reliable data that could give a direct index of solar variability before the beginning of the sunspot-number series. Therefore, one must use indirect proxies, i.e., quantitative parameters, which can be measured nowadays but represent different effects of solar magnetic activity in the past. It is common to use, for this purpose, signatures of terrestrial indirect effects induced by variable solar-magnetic activity, that is stored in natural archives. Such traceable signatures can be related to nuclear (used in the cosmogenic-isotope method) or chemical (used, e.g., in the nitrate method) effects caused by cosmic rays (CRs) in the Earth's atmosphere, lunar rocks or meteorites.

The most common proxy of solar activity is formed by the data on cosmogenic radionuclides (e.g., 10Be and 14C), which are produced by cosmic rays in the Earth's atmosphere (e.g, Stuiver and Quay, 1980; Beer et al., 1990; Bard et al., 1997; Beer, 2000). Other cosmogenic nuclides, which are used in geological and paleomagnetic dating, are less suitable for studies of solar activity (see e.g., Beer, 2000; Beer et al., 2012). Cosmic rays are the main source of cosmogenic nuclides in the atmosphere (excluding anthropogenic factors during the last decades) with the maximum production being in the upper troposphere/stratosphere. After a complicated transport in the atmosphere, the cosmogenic isotopes are stored in natural archives such as polar ice, trees, marine sediments, etc. This process is also affected by changes in the geomagnetic field and climate. Cosmic rays experience heliospheric modulation due to solar wind and the frozen-in solar magnetic field. The intensity of modulation depends on solar activity and, therefore, cosmic-ray flux and the ensuing cosmogenic isotope intensity depends inversely on solar activity. An important advantage of the cosmogenic data is that primary archiving is done naturally in a similar manner throughout the ages, and these archives are measured nowadays in laboratories using modern techniques. If necessary, all measurements can be repeated and improved, as has been done for some radiocarbon samples. In contrast to fixed historical archival data (such as sunspot or auroral observations) this approach makes it possible to obtain homogeneous data sets of stable quality and to improve the quality of data with the invention of new methods (such as accelerator mass spectrometry). Cosmogenic isotope data is the main regular indicator of solar activity on the very long-term scale but it cannot resolve the details of individual solar cycles. The redistribution of nuclides in terrestrial reservoirs and archiving may be affected by local and global climate/circulation processes, which are, to a large extent, unknown for the past. However, a combined study of different nuclides data, whose responses to terrestrial effects are very different, may allow for disentangling external and terrestrial signals.

3.1 The physical basis of the method 3.1.1 Heliospheric modulation of cosmic rays

The flux of cosmic rays (highly energetic fully ionized nuclei) is considered roughly constant (at least at the time scales relevant for the present study) in the vicinity of the Solar system. However, before reaching the vicinity of Earth, galactic cosmic rays experience complicated transport in the heliosphere that leads to modulation of their flux. Heliospheric transport of GCR is described by Parker's theory (Parker, 1965; Toptygin, 1985) and includes four basic processes: the diffusion of particles due to their scattering on magnetic inhomogeneities, the convection of particles by out-blowing solar wind, adiabatic energy losses in expanding solar wind, drifts of particles in the magnetic field, including the gradient-curvature drift in the regular heliospheric magnetic field, and the drift along the heliospheric current sheet, which is a thin magnetic interface between the two heliomagnetic hemispheres. Because of variable solar-magnetic activity, CR flux in the vicinity of Earth is strongly modulated (see Figure 3). The most prominent feature in CR modulation is

the 11-year cycle, which is in inverse relation to solar activity. The 11-year cycle in CR is delayed (from a month up to two years) with respect to the sunspots (Usoskin et al., 1998). The time profile of cosmic-ray flux as measured by a neutron monitor (NM) is shown in Figure 3 (panel b) together with the sunspot numbers (panel a). Besides the inverse relation between them, some other features can also be noted. A 22-year cyclicity manifests itself in cosmic-ray modulation through the alteration of sharp and flat maxima in cosmic-ray data, originated from the charge-dependent drift mechanism. One may also note short-term fluctuations, which are not directly related to sunspot numbers but are driven by interplanetary transients caused by solar eruptive events, e.g., flares or CMEs. An interesting feature is the increase of CR flux in 2009, when it was the highest ever recorded by NMs (Moraal and Stoker, 2010), as caused by the favorable heliospheric conditions (unusually weak heliospheric magnetic field and the flat heliospheric current sheet) (McDonald et al., 2010). For the previous 50 years of high and roughly-stable solar activity, no trends have been observed in CR data; however, as will be discussed later, the overall level of CR has changed significantly on the centurial-millennial timescales.

Figure 3: Cyclic variations since 1951. Panel a: Time profiles of sunspot numbers (http://sidc.oma. be/sunspot-data/); Panel b: Cosmic-ray flux as the count rate of a polar neutron monitor (Oulu NM, Climax NM data used before 1964), 100% NM count rate corresponds to May 1965.

Full solution of the CR transport problems is a complicated task and requires sophisticated 3D time-dependent self-consistent modelling. However, the problem can be essentially simplified for applications at a long-timescale. An assumption on the azimuthal symmetry (requires times longer that the solar-rotation period) and quasi-steady changes reduces it to a 2D quasi-steady problem. Further assumption of the spherical symmetry of the heliosphere reduces the problem to a 1D case. This approximation can be used only for rough estimates, since it neglects the drift effect, but it is useful for long-term studies, when the heliospheric parameters cannot be evaluated independently. Further, but still reasonable, assumptions (constant solar-wind speed, roughly power-law CR energy spectrum, slow spatial changes of the CR density) lead to the force-field approximation (Gleeson and Axford, 1968), which can be solved analytically. The differential

intensity Ji of the cosmic-ray nuclei of type i with kinetic energy T at 1 AU is given in this case as

MT,4) = JLis,i(T + ^+ ^ ^), (3)

where = (Zje/Ai)4 for a cosmic nuclei of j-th type (charge and mass numbers are Zi and Ai), T and 4 are expressed in MeV/nucleon and in MV, respectively, Tr = 938 MeV. T is the CR particle's kinetic energy, and 4 is the modulation potential. The local interstellar spectrum (LIS) JLIS forms the boundary condition for the heliospheric transport problem. Since LIS is not measured directly, i.e., outside the heliosphere, it is not well known in the energy range affected by CR modulation (below 100 GeV). Presently-used approximations for LIS (e.g., Garcia-Munoz et al., 1975; Burger et al., 2000; Webber and Higbie, 2003, 2009) agree with each other for energies above 20 GeV but may contain uncertainties of up to a factor of 1.5 around 1 GeV. These uncertainties in the boundary conditions make the results of the modulation theory slightly model-dependent (see discussion in Usoskin et al., 2005; Herbst et al., 2010) and require the LIS model to be explicitly cited. This approach gives results, which are at least dimensionally consistent with the full theory and can be used for long-term studies1 (Usoskin et al., 2002b; Caballero-Lopez and Moraal, 2004). Differential CR intensity is described by the only time-variable parameter, called the modulation potential 4, which is mathematically interpreted as the averaged rigidity (i.e., the particle's momentum per unit of charge) loss of a CR particle in the heliosphere. However, it is only a formal spectral index whose physical interpretation is not straightforward, especially on short timescales and during active periods of the sun (Caballero-Lopez and Moraal, 2004). Despite its cloudy physical meaning, this force-field approach provides a very useful and simple single-parametric approximation for the differential spectrum of GCR, since the spectrum of different GCR species directly measured near the Earth can be perfectly fitted by Equation 3 using only the parameter 4 in a wide range of solar activity levels (Usoskin et al., 2011). Therefore, changes in the whole energy spectrum (in the energy range from 100 MeV/nucleon to 100 GeV/nucleon) of cosmic rays due to the solar modulation can be described by this single number within the framework of the adopted LIS. The concept of modulation potential is a key concept for the method of solar-activity reconstruction by cosmogenic isotope proxy as it makes it possible to parameterize the GCR with one single parameter.

3.1.2 Geomagnetic shielding

Cosmic rays are charged particles and therefore are affected by the Earth's magnetic field. Thus the geomagnetic field puts an additional shielding on the incoming flux of cosmic rays. The shielding effect of the geomagnetic field is usually expressed in terms of the cutoff rigidity Pc, which is the minimum rigidity a vertically incident CR particle must posses (on average) in order to reach the ground at a given location and time (Cooke et al., 1991). Neglecting such effects as the East-West asymmetry, which is roughly averaged out for the isotropic particle flux, or nondipole magnetic momenta, which decay rapidly with distance, one can come to a simple approximation, called the Stormer's equation, that describes the vertical geomagnetic cutoff rigidity Pc:

Pc « 1.9 M (Ro/Rf cos4 XG [GV], (4)

where M is the geomagnetic dipole moment (in 1025 G cm3), Ro is the Earth's mean radius, R is the distance from the given location to the dipole center, and Xq is the geomagnetic latitude. The cutoff concept works like a Heaviside step-function so that all cosmic rays whose rigidity is below the cutoff are not allowed to enter the atmosphere while all particles with higher rigidity

1 Note that the famous work by Castagnoli and Lal (1980) contains an inconsistency in the force-field formula —

see details in Usoskin et al. (2005),

can penetrate. This approximation provides a good compromise between simplicity and reality, especially when using the eccentric dipole description of the geomagnetic field (Fraser-Smith, 1987). The eccentric dipole has the same dipole moment and orientation as the centered dipole, but the dipole's center and consequently the poles, defined as crossings of the axis with the surface, are shifted with respect to geographical ones.

The shielding effect is the strongest at the geomagnetic equator, where the present-day value of Pc may reach up to 17 GV in the region of India. There is almost no cutoff in the geomagnetic polar regions (Ag > 60°). However, even in the latter case the atmospheric cutoff becomes important, i.e., particles must have rigidity above 0.5 GV in order to initiate the atmospheric cascade which can reach ground (see Section 3.1.3).

The geomagnetic field is seemingly stable on the short-term scale, but it changes essentially on centurial-to-millennial timescales (e.g., Korte and Constable, 2006). Such past changes can be evaluated based on measurements of the residual magnetization of independently-dated samples. These can be paleo- (i.e., natural stratified archives such as lake or marine sediments or volcanic lava) or archaeological (e.g., clay bricks that preserve magnetization upon baking) samples. Most paleo-magnetic data preserve not only the magnetic field intensity but also the direction of the local field, while archeo-magnetic samples provide information on the intensity only. Using a large database of such samples, it is possible to reconstruct (under reasonable assumptions) the large-scale magnetic field of the Earth. Data available provides good global coverage for the last 3 millennia, allowing for a reliable paleomagnetic reconstruction of the true dipole moment (DM) or virtual dipole moment2 (VDM) and its orientation (the ArcheoInt collaboration - Genevey et al., 2008). Less precise, but still reliable reconstructions of the DM and its orientation are possible for the last seven millennia (the CALS7K.2 model by Korte and Constable, 2005), however they may slightly underestimate the dipole moment, especially in the earlier part of the period (Korte and Constable, 2008). Directional paleomagnetic reconstruction are less reliable on a longer timescale, because of the spatial sparseness of the paleo/archeo-magnetic samples in the earlier part of the Holocene (Korte et al., 2011). Some paleomagnetic reconstructions are shown in Figure 4. All paleomagnetic models depict a similar long-term trend - an enhanced intensity during the period between 1500 BC and 500 AD and a significantly lower field before that.

Changes in the dipole moment M inversely modulate the flux of CR at Earth, with strong effects in tropical regions and globally. The migration of the geomagnetic axis, which changes the geomagnetic latitude Xq of a given geographical location is also important; while not affecting the global flux of CR, it can dramatically change the CR effect regionally, especially at middle and high latitude. These changes affect the flux of CR impinging on the Earth's atmosphere both locally and globally and must be taken into account when reconstructing solar activity from terrestrial proxy data (Usoskin et al., 2008, 2010). Accounting for these effects is quite straightforward provided the geomagnetic changes in the past are known independently, e.g., from archeo and paleo-magnetic studies (Donadini et al., 2010). However, because of progressively increasing uncertainties of paleomagnetic reconstructions back in time, it presently forms the main difficulty for the proxy method on the long-term scale (Snowball and Muscheler, 2007), especially in the early part of the Holocene. On the other hand, the geomagnetic field variations are relatively well known for the last few millennia (Genevey et al., 2008; Korte and Constable, 2008).

3.1.3 Cosmic-ray—induced atmospheric cascade

When an energetic CR particle enters the atmosphere, it first moves straight in the upper layers, suffering mostly from ionization energy losses that lead to the ionization of the ambient rarefied air and gradual deceleration of the particles. However, after traversing some amount of matter (the

2 The concept of VDM assumes that the geomagnetic dipole is centered at the planet's center and its axis is aligned with the true magnetic axis.

Figure 4: Geomagnetic field intensity over millennia: VADM compilation by Yang et al. (2000, - Y00 curve with 1a statistical errors of the sample distribution); dipole moment according to Hongre et al. (1998,

- HBI red dots) since 800 AD, dipole moment according to CALS7K.2 model (Korte and Constable, 2005,

- K05 magenta curve with 1a error band) as well as a recent ArcheoInt compilations of VADM/VDM (Genevey et al., 2008, - A08 azure diamonds).

nuclear interaction mean-free path is on the order of 100 g/cm2 for a proton in the air) the CR particle may collide with a nucleus in the atmosphere, producing a number of secondaries. These secondaries have their own fate in the atmosphere, in particular they may suffer further collisions and interactions forming an atmospheric cascade (e.g., Dorman, 2004). Because of the thickness of the Earth's atmosphere (1033 g/cm2 at sea level) the number of subsequent interactions can be large, leading to a fully-developed cascade (also called an air shower) consisting of secondary rather than primary particles. A schematic view of the atmospheric cascade is shown in Figure 5. Three main components can be separated in the cascade:

• The "hadronic" nucleonic component is formed by the products of nuclear collisions of primary cosmic rays and their secondaries with the atmospheric nuclei, and consists mostly of superthermal protons and neutrons.

• The "soft" or electromagnetic component consists of electrons, positrons and photons.

• The "hard" or muon component consists mostly of muons; pions are short lived and decay almost immediately upon production, feeding muons and the "soft" component.

The development of the cascade depends mostly on the amount of matter traversed and is usually linked to residual atmospheric depth, which is very close to the static barometric pressure, rather than to the actual altitude, that may vary depending on the exact atmospheric density profile.

Cosmogenic isotopes are a by-product of the hadronic branch of the cascade (details are given below). Accordingly, in order to evaluate cosmic-ray flux from the cosmogenic isotope data, one needs to know the physics of cascade development. Several models have been developed for this cascade, in particular its hadronic branch with emphasis on the generation of cosmogenic isotope production. The first models were simplified quasi-analytical (e.g., Lingenfelter, 1963; O'Brien and Burke, 1973) or semi-empirical models (e.g., Castagnoli and Lal, 1980). With the fast advance of computing facilities it became possible to exploit the best numerical method suitable for such problems - Monte-Carlo (e.g., Masarik and Beer, 1999, 2009; Webber and Higbie, 2003; Webber et al., 2007; Usoskin and Kovaltsov, 2008a; Kovaltsov and Usoskin, 2010; Kovaltsov et al., 2012;

Figure 5: Schematic view of an atmospheric cascade caused by energetic cosmic rays in the atmosphere. Left-to-right are denoted, respectively, the soft, muon and hadronic components of the cascade. Symbols "N, p, n, y, -k, e-, e+, and 7" denote nuclei, protons, neutrons, muons, pions, electrons, positrons, and photons, respectively. Stars denote nuclear collisions, ovals - decay processes. This sketch does not represent the full development of the cascade and serves solely as an illustration for the processes discussed in the text. Image reproduced by permission from Usoskin (2011), copyright by SAIt.

Argento et al., 2013). The fact that models, based on different independent Monte-Carlo packages, namely, a general GEANT tool and a specific CORSIKA code, yield similar results provides additional verification of the approach.

3.1.4 Transport and deposition

A scheme for the transport and redistribution of the two most useful cosmogenic isotopes, 14 C and 10Be, is shown in Figure 6. After a more-or-less similar production, the two isotopes follow quite different fates, as discussed in detail in Sections 3.2.3 and 3.3.3. Therefore, expected terrestrial effects are quite different for the isotopes and comparing them with each other can help in disentangling solar and climatic effects (see Section 3.7.3). A reader can find great detail also in a book by Beer et al. (2012).

3.2 Radioisotope 14C

The most commonly used cosmogenic isotope is radiocarbon 14 C. This radionuclide is an unstable isotope of carbon with a half-life (T1/2) of about 5730 years. Since the radiocarbon method is extensively used in other science disciplines where accurate dating is a key issue (e.g., archeology, paleoclimatology, quaternary geology), it was developed primarily for this task. The solar-activity-reconstruction method, based on radiocarbon data, was initially developed as a by-product of the dating techniques used in archeology and Quaternary geology, in an effort to improve the quality of the dating by means of better information on the 14C variable source function. The present-day radiocarbon calibration curve, based on a dendrochronological scale, uninterruptedly covers

Cosmic Rays

Figure 6: Schematic representation of 14C (left) and 10Be (right) production chains. The flux of cosmic rays impinging on the Earth is affected by both heliospheric modulation and geomagnetic field changes. The climate may affect the redistribution of the isotopes between different reservoirs. Dashed line denotes a possible influence of solar activity on climate.

the whole Holocene (and extending to 50,000 BP - Reimer et al., 2009) and provides a solid quantitative basis for studying solar activity variations on the multi-millennial time scale.

3.2.1 Measurements

Radiocarbon is usually measured in tree rings, which allows an absolute dating of the samples by means of dendrochronology. Using a complicated technique, the 14C activity3 A is measured in an independently dated sample, which is then corrected for age as

A* = A ■ exp

/ 0.693A

V T1/2 )

where t and T1/2 are the age of the sample and the half-life of the isotope, respectively. Then the relative deviation from the standard activity Aa of oxalic acid (the National Bureau of Standards) is calculated:

S14C = ^~ ■ 1000 . (6)

After correction for the carbon isotope fractionating (account for the 13C isotope) of the sample, the radiocarbon value of A14C is calculated (see details in Stuiver and Pollach, 1977).

A14C = ¿14C - (2 ■ ¿13C + 50) ■ (1 + ¿14C/1000),

where £13C is the per mille deviation of the 13C content in the sample from that in the standard belemnite sample calculated similarly to Equation 6. The value of A14C (measured in per mille %o) is further used as the index of radiocarbon relative activity. The series of A14C for the Holocene is presented in Figure 7A as published by the INTCAL04 collaboration of 21 dating laboratories

3 Isotope's activity quantifies (a) in the radiometric 14 technique its decay rate, and is usually given in terms of disintegrations per minute per gram of carbon, and (b) in the AMS technique, the 14C/12C ratio, all normalized to the standard.

as a result of systematic precise measurements of dated samples from around the world (Reimer et al, 2004). The most recent INTCAL09 dataset is available at IntCal09.htm.

Figure 7: Radiocarbon series for the Holocene. Upper panel: Measured content of A14C in tree rings by INTCAL-98/04 collaboration (Stuiver et al, 1998; Reimer et al, 2004). The long-term trend is caused by the geomagnetic field variations and the slow response of the oceans. Lower panel: Production rate of 14C in the atmosphere, reconstructed from the measured A14C (Usoskin and Kromer, 2005).

A potentially interesting approach has been made by Lai et al. (2005), who measured the amount of 14C directly produced by CR in polar ice. Although this method is free of the carbon-cycle influence, the first results, while being in general agreement with other methods, are not precise.

3.2.2 Production

The main source of radioisotope 14C (except anthropogenic sources during the last decades) is cosmic rays in the atmosphere. It is produced as a result of the capture of a thermal neutron by atmospheric nitrogen

14N + n 14C + p. (8)

Neutrons are always present in the atmosphere as a product of the cosmic-ray-induced cascade (see Section 3.1.3) but their flux varies in time along with the modulation of cosmic-ray flux. This provides continuous source of the isotope in the atmosphere, while the sinks are isotope decay and transport into other reservoirs as described below (the carbon cycle).

The connection between the cosmogenic-isotope-production rate, Q, at a given location (quantified via the geomagnetic latitude Ag ) and the cosmic-ray flux is given by


Q= S{P,<f>) Y(P) dP, (9)


where Pc is the local cosmic-ray-rigidity cutoff (see Section 3.1.2), S(P, 4>) is the differential energy spectrum of CR (see Section 3.1.1) and Y(P) is the differential yield function of cosmogenic isotope production, calculated using a Monte-Carlo simulation of the cosmic-ray-induced atmospheric cascade (Kovaltsov et al., 2012). Because of the global nature of the carbon cycle and its long attenuation time, the radiocarbon is globally mixed before the final deposition, and Equation 9 should be integrated over the globe. The yield function Y(P) of the 14C production is shown in Figure 8A together with those for 10Be (see Section 3.3.2) and for a ground-based neutron monitor (NM), which is the main instrument for studying cosmic-ray variability during the modern epoch. One can see that the yield function increases with the energy of CR. On the other hand, the energy spectrum of CR decreases with energy. Accordingly, the differential production rate (i.e., the product of the yield function and the spectrum, F = Y • S - the integrand of Equation 9), shown in Figure 8B, is more informative. The differential production rate reflects the sensitivity to cosmic rays, and the total production rate is simply an integral of F over energy above the geomagnetic threshold.

Energy (GeV)

Figure 8: Differential production rate for cosmogenic isotopes and ground-based neutron monitors as a function of cosmic-ray energy. Panel A: Yield functions of the globally averaged and polar 10Be production (Webber and Higbie, 2003), global 14C production (Kovaltsov et al., 2012), polar neutron monitor (Clem and Dorman, 2000) as well as the energy spectrum of galactic cosmic protons for medium modulation (0 = 550 MV). Panel B: The differential production rate for global and polar 10Be production, global 14C production, and polar neutron monitor.

Thanks to the development of atmospheric cascade models (Section 3.1.3), there are numerical

models that allow one to compute the radiocarbon production rate as a function of the modulation potential 4> and the geomagnetic dipole moment M. The overall production of 14 C is shown in Figure 9.

The production rate of radiocarbon, Qi4C, can vary as affected by different factors (see, e.g., Damon and Sonett, 1991):

• Variations of the cosmic-ray flux on a geological timescale due to the changing galactic background (e.g., a nearby supernova explosion or crossing the dense galactic arm).

• Secular-to-millennial variations are caused by the slowly-changing geomagnetic field. This is an important component of the variability, which needs to be independently evaluated from paleo and archeo-magnetic studies.

• Modulation of cosmic rays in the heliosphere by solar magnetic activity. This variation is the primary aim of the present method.

• Short-term variability of CR on a daily scale (suppression due to interplanetary transients or enhancement due to solar energetic-particle events) can be hardly resolved in radiocarbon data.

Figure 9: Globally-averaged production rate of 14C as a function of the modulation potential <f> and geomagnetic dipole moment M, computed using the yield function by Kovaltsov et al. (2012), LIS by Burger et al. (2000) and cosmic-ray-modulation model by Usoskin et al. (2005). Another often used model (Masarik and Beer, 2009) yields a similar result.

Therefore, the production rate of 14C in the atmosphere can be modelled for a given time (namely, the modulation potential and geomagnetic dipole moment) and location. The global production rate Q is then obtained as a result of global averaging.

There is still a small unresolved discrepancy in the absolute value of the modeled 14 C production rate. Different models yield the global-average production rate of 1.7-2.3 atoms cm-2 s_1 (see,

e.g., O'Brien, 1979; Masarik and Beer, 1999; Goslar, 2001; Usoskin et al., 2006b; Kovaltsov et al., 2012, and references therein), which is consistent with a direct estimate of the radiocarbon reservoir, based on analyses of the specific 14C activity on the ground, 1.76-2.0 (Suess, 1965; Damon et al., 1978; O'Brien, 1979). On the other hand, the steady-state production calculated from the

14C inventory in the carbon-cycle model (see Section 3.2.3) typically yields 1.6-1.7 atoms cm-2 s-1 for the pre-industrial period (e.g., Goslar, 2001, and references therein). Thus, results obtained from the carbon cycle models and production models agree only marginally in the absolute values, and this needs further detailed studies. The situation has been improved recently, when the newest numerical model (Kovaltsov et al., 2012) yields the global average production of 1.64 and 1.88 atoms cm-2 s-1 for the modern and pre-industrial periods, respectively. In 14C-based reconstructions, the pre-industrial steady-state production is commonly used.

3.2.3 Transport and deposition

Upon production cosmogenic radiocarbon gets quickly oxidized to carbon dioxide CO2 and takes part in the regular carbon cycle of interrelated systems: atmosphere-biosphere-ocean (Figure 6). Because of the long residence time, radiocarbon becomes globally mixed in the atmosphere and involved in an exchange with the ocean. It is common to distinguish between an upper layer of the ocean, which can directly exchange CO2 with the air and deeper layers. The measured A14C comes from the biosphere (trees), which receives radiocarbon from the atmosphere. Therefore, the processes involved in the carbon cycle are quite complicated. The carbon cycle is usually described using a box model (Oeschger et al., 1974; Siegenthaler et al., 1980), where it is represented by fluxes between different carbon reservoirs and mixing within the ocean reservoir(s), as shown in Figure 10. Production and radioactive decay are also included in box models. Free parameters in a typical box model are the 14C production rate Q, the air-sea exchange rate (expressed as turnover rate k), and the vertical-eddy-diffusion coefficient K, which quantifies ocean ventilation. Starting from the original representation (Oeschger et al., 1974), a variety of box models have been developed, which take into account subdivisions of the ocean reservoir and direct exchange between the deep ocean and the atmosphere at high latitudes. More complex models, including a diffusive approach, are able to simulate more realistic scenarios, but they require knowledge of a large number of model parameters. These parameters can be evaluated for the present time using the bomb test - studying the transport and distribution of the radiocarbon produced during the atmospheric nuclear tests. However, for long-term studies, only the production rate is considered variable, while the gas-exchange rate and ocean mixing are kept constant. Under such assumptions, there is no sense in subdividing reservoirs or processes, and a simple carbon box model is sufficient.

Using the carbon cycle model and assuming that all its parameters are constant in time, one can evaluate the production rate Q from the measured A14C data. This assumption is well validated for the the Holocene (Damon et al., 1978; Stuiver et al., 1991) as there is no evidence of considerable oceanic change or other natural variability of the carbon cycle (Gerber et al., 2002), and accordingly all variations of A14C predominantly reflect the production rate. This is supported by the strong similarity of the fluctuations of the 10Be data in polar ice cores (Section 3.3) compared to 14C, despite their completely different geochemical fate (Bard et al., 1997; Steinhilber et al., 2012). However, the changes in the carbon cycle during the last glaciation and deglaciation were dramatic, especially regarding ocean ventilation; this and the lack of independent information about the carbon cycle parameters, make it hardly possible to qualitatively estimate solar activity from 14 C before the Holocene.

First attempts to extract information on production-rate variations from measured A14C were based on simple frequency separations of the signals. All slow changes were ascribed to climatic and geomagnetic variations, while short-term fluctuations were believed to be of solar origin. This was done by removing the long-term trend from the A14C series and claiming the residual as being

Figure 10: A 12-box model of the carbon cycle (Broeker and Peng, 1986; Siegenthaler et al., 1980). The number on each individual box is the steady-state A14C of this particular reservoir expressed in per mil. Image reproduced by permission from Bard et al. (1997), copyright by Elsevier.

a series of solar variability (e.g., Peristykh and Damon, 2003). This oversimplified approach was natural at earlier times, before the development of carbon cycle models, but later it was replaced by the inversion of the carbon cycle (i.e., the reconstruction of the production rate from the measured 14C concentration). Although mathematically this problem can be solved correctly as a system of linear differential equations, the presence of fluctuating noise with large magnitude makes it not straightforward, since the time derivative cannot be reliably identified leading thus to possible amplification of the high-frequency noise in A14C data. One traditional approach (e.g., Stuiver and Quay, 1980) is based on an iterative procedure, first assuming a constant production rate, and then fitting the calculated A14C variations to the actual measurements using a feedback scheme. A concurrent approach based on the presentation of the carbon cycle as a Fourier filter (Usoskin and Kromer, 2005) produces similar results. Roughly speaking, the carbon cycle acts as an attenuating and delaying filter for the 14C signal (see Figure 11). The higher the frequency is, the greater the signal is attenuated. In particular, the large 11-year solar cycle expected in the 14C is attenuated by a factor of hundred in the measured A14C data, making it hardly detectable. Because of the slow oceanic response, the 14C data is also delayed with respect to the production signal. The production rate QuC for the Holocene is shown in Figure 7 and depicts both short-term fluctuations as well as slower variations, mostly due to geomagnetic field changes (see Section 3.2.5).

3.2.4 The Suess effect and nuclear bomb tests

Unfortunately, cosmogenic 14 C data cannot be easily used for the last century, primarily because of the extensive burning of fossil fuels. Since fossil fuels do not contain 14C, the produced CO2 dilutes the atmospheric 14CO2 concentration with respect to the pre-industrial epoch. Therefore, the measured A14C cannot be straightforwardly translated into the production rate Q after the late 19th century, and a special correction for fossil fuel burning is needed. This effect, known as the Suess effect (e.g., Suess, 1955), can be up to —25%o in A14C in 1950 (Tans et al., 1979), which is an order of magnitude larger than the amplitude of the 11-year cycle of a few %o. Moreover, while the cosmogenic production of 14 C is roughly homogeneous over the globe and time, the use of fossil fuels is highly nonuniform (e.g., de Jong and Mook, 1982) both spatially (de-

1E-4 1E-3 1E-2

Frequency (year- )

<5 30 A £ k

Frequency (year- )

Figure 11: The frequency characteristics of the carbon cycle: attenuation (left-hand panel) and phase shift (right-hand panel) as a function of the frequency of the 14C production signal. Lines stand for a classical Oeschger-Siegenthaler box model (Siegenthaler et al., 1980), and open circles for a sophisticated PANDORA model (Bard et al., 1997).

veloped countries, in the northern hemisphere) and temporarily (World Wars, Great Depression, industrialization, etc.). This makes it very difficult to perform an absolute normalization of the radiocarbon production to the direct measurements. Sophisticated numerical models (e.g., Sabine et al., 2004; Mikaloff Fletcher et al., 2006) aim to account for the Suess effect and make good progress. However, the results obtained indicate that the determination of the Suess effect does not yet reach the accuracy required for the precise modelling and reconstruction of the 14 C production for the industrial epoch. As noted by Matsumoto et al. (2004), "... not all is well with the current generation of ocean carbon cycle models. At the same time, this highlights the danger in simply using the available models to represent state-of-the-art modeling without considering the credibility of each model." Note that the atmospheric concentration of another carbon isotope 13C is partly affected by land use, which has also been modified during the last century.

Another anthropogenic activity greatly disturbing the natural variability of 14C is related to the atmospheric nuclear bomb tests actively performed in the 1960s. For example, the radiocarbon concentration nearly doubled in the early 1960s in the northern hemisphere after nuclear tests performed by the USSR and the USA in 1961 (Damon et al., 1978). On one hand, such sources of momentary spot injections of radioactive tracers (including 14C) provide a good opportunity to verify and calibrate the exchange parameters for different carbon-cycle reservoirs and circulation models (e.g., Bard et al., 1987; Sweeney et al., 2007). Thus, the present-day carbon cycle is more-or-less known. On the other hand, the extensive additional production of isotopes during nuclear tests makes it hardly possible to use the 14C as a proxy for solar activity after the 1950s (Joos, 1994).

These anthropogenic effects do not allow one to make a straightforward link between pre-industrial data and direct experiments performed during more recent decades. Therefore, the question of the absolute normalization of 14C model is still open (see, e.g., the discussion in Solanki et al., 2004, 2005; Muscheler et al., 2005).

3.2.5 The effect of the geomagnetic field

As discussed in Section 3.1.2, knowledge of geomagnetic shielding is an important aspect of the cos-mogenic isotope method. Since radiocarbon is globally mixed in the atmosphere before deposition, its production is affected by changes in the geomagnetic dipole moment M, while magnetic-axis migration plays hardly any role in 14C data.

The crucial role of paleomagnetic reconstructions has long been known (e.g., Elsasser et al.,

1956; Kigoshi and Hasegawa, 1966). Many earlier corrections for possible geomagnetic-field changes were performed by detrending the measured A14C abundance or production rate Q (Stuiver and Quay, 1980; Voss et al., 1996; Peristykh and Damon, 2003), under the assumption that geomagnetic and solar signals can be disentangled from the production in the frequency domain. Accordingly, the temporal series of either measured A14C or its production rate Q is decomposed into the slow changing trend and faster oscillations. The trend is supposed to be entirely due to geomagnetic changes, while the oscillations are ascribed to solar variability. Such a method, however, obliterates all information on possible long-term variations of solar activity. Simplified empirical correction factors were also often used (e.g., Stuiver and Quay, 1980; Stuiver et al., 1991). The modern approach is based on a physics-based model (e.g., Solanki et al., 2004; Vonmoos et al., 2006) and allows the quantitative reconstruction of solar activity, explicitly using independent reconstructions of the geomagnetic field. In this case the major source of errors in solar activity reconstructions is related to uncertainties in the paleomagnetic data (Snowball and Muscheler, 2007). These errors are insignificant for the last millennium (Usoskin et al., 2006a), but become increasingly important for earlier times.

3.3 Cosmogenic isotope 10Be

3.3.1 Measurements

The cosmogenic isotope 10Be is useful for long-term studies of solar activity because of its long halflife of around 1.5 x 106 years. Its concentration is usually measured in stratified ice cores allowing for independent dating. Because of its long life, the beryllium concentration is difficult to measure by the decay rate (Beer, 2000). Accordingly, the 10Be/9Be ratio needs to be precisely measured at an accuracy better than 10-13. This can be done using AMS (Accelerator Mass Spectrometry) technique, which makes the measurements complicated and expensive. Correction for the decay is straightforward and does not include isotope fractionating. From the measured samples, first the 10Be concentration is defined, usually in units of 104 atoms/g. Sometimes, a correction for the snow precipitation amount is considered leading to the observable 10Be flux, which is the number of atoms, precipitating to the surface per cm2 per second.

There exist different 10Be series suitable for studies of long-term solar activity, coming from ice cores in Greenland and Antarctica. They have been obtained from different cores with different resolutions, and include data from Milcent, Greenland (Beer et al., 1983), Camp Century, Greenland (Beer et al., 1988), Dye 3, Greenland (Beer et al., 1990), Dome Concordia and South Pole, Antarctica (Raisbeck et al., 1990), GRIP, Greenland (Yiou et al., 1997), GISP2, Greenland (Finkel and Nishiizumi, 1997), Dome Fuji, Antarctica (Horiuchi et al., 2007, 2008), Dronning Maud Land, Antarctica (Ruth et al., 2007), etc. We note that data on 10Be in other archives, e.g., lake sediments, is usually more complicated to interpret because of the potential influence of the climate (Horiuchi et al., 1999; Belmaker et al., 2008).

Details of the 10Be series and their comparison with each other can be found in Beer (2000), Muscheler et al. (2007), and Beer et al. (2012).

3.3.2 Production

The isotope 10Be is produced as a result of spallation of atmospheric nitrogen and oxygen (carbon is less abundant by far in the atmosphere and makes a negligible contribution) by the nucleonic component of the cosmic-ray-induced atmospheric cascade (Section 3.1.3).

A small contribution may also exist from photo-nuclear reactions (Bezuglov et al., 2012). The cross section (a few mb) of the spallation reactions is almost independent of the energy of impacting particles and has a threshold of about 15 MeV. Thus, the production of 10Be is defined mostly by the multiplicity of the nucleonic component, which increases with the energy of primary cosmic rays

(see Figure 8). Maximum production occurs at an altitude of 10-15 km due to a balance between the total energy of the cascade (which increases with altitude) and the number of secondaries (decreasing with altitude). Most of the global 10Be is produced in the stratosphere (55-70%) and the rest in the troposphere (Lal and Peters, 1967; Masarik and Beer, 1999, 2009; Usoskin and Kovaltsov, 2008a; Kovaltsov and Usoskin, 2010).

Computation of 10Be isotope production is straightforward, provided a model of the atmospheric cascade is available. The first consistent model was developed by D. Lal et al. (Bhandari et al., 1966; Lal and Peters, 1967; Lal and Suess, 1968), using an empirical approach based on fitting simplified model calculations to measurements of the isotope concentrations and "star" (inelastic nuclear collisions) formations in the atmosphere. Next was an analytical model by O'Brien (1979), who solved the problem of the GCR-induced cascade in the atmosphere using an analytical stationary approximation in the form of the Boltzmann equation. Those models were based on calculating the rate of inelastic collisions or "stars" and then applying the mean spallation yield per "star". A new step in the modelling of isotope production was made by Masarik and Beer (1999), who performed a full Monte-Carlo simulation of a GCR-initiated cascade in the atmosphere and used cross sections of spallation reactions directly instead of the average "star" efficiency. Modern models (Webber and Higbie, 2003; Webber et al., 2007; Usoskin and Kovaltsov, 2008a; Kovaltsov and Usoskin, 2010) are based on a full Monte-Carlo simulation of the atmospheric cascade, using improved cross sections. The global production rate of 10Be is about 0.02-0.03 atoms cm-2 s-1 (Masarik and Beer, 1999; Webber et al., 2007; Kovaltsov and Usoskin, 2010), which is lower than that for 14C by two orders of magnitude (about 2 atoms cm-2 s-1; see Section 3.2.2). The yield function of 10Be production is shown in Figure 8A and the differential production rate in Figure 8B. One can see that the peak of 10Be sensitivity, especially in polar regions, is shifted towards lower energies (below 1 GeV) compared with a neutron monitor. This implies that the 10Be isotope is relatively more sensitive to less energetic CR and is, therefore, more affected by solar energetic particles (Usoskin et al., 2006b). Comparison of model computations with direct beryllium production experiments (Usoskin and Kovaltsov, 2008a; Kovaltsov and Usoskin, 2010), and also the results of modelling of the short-living 7Be isotope (Usoskin et al., 2009a) suggest that some numerical models (Masarik and Beer, 1999; Webber and Higbie, 2003; Webber et al., 2007) tend to underestimate the production.

Although the production of 10Be can be more or less precisely modelled, a simple normalization "surface", similar to that shown in Figure 9 for 14C, is not easy to produce because of partial mixing in the atmosphere (see Section 3.3.3). Simplified models, assuming either only global (e.g., Beer, 2000) or polar production (Bard et al., 1997; Usoskin et al., 2004), have been used until recently. However, it has been recognized that a more realistic model of the limited atmospheric mixing should be used. Without detailed knowledge of 10Be transport in the atmosphere, it is impossible to relate the quantitatively-measured concentration to the production (as done for 14C using the carbon cycle), and one has to assume that the measured abundance is proportional (with an unknown coefficient) to the production rate in a specific geographical region (see Section 3.3.3).

3.3.3 Atmospheric transport

After production, the 10Be isotope has a seemingly simple (Figure 6) but difficult-to-account-for fate in the atmosphere. Its atmospheric residence time depends on scavenging, stratosphere-troposphere exchange and inter-tropospheric mixing (e.g., McHargue and Damon, 1991). Soon after production, the isotope becomes attached to atmospheric aerosols and follows their fate. In addition, it may be removed from the lower troposphere by wet deposition (rain and snow). The mean residence time of the aerosol-bound radionuclide in the atmosphere is quite different for the troposphere, being a few weeks, and stratosphere, where it is one to two years (Raisbeck et al., 1981). Accordingly, 10Be produced in the troposphere is deposited mostly locally, i.e., in the

polar regions, while stratospheric 10Be can be partly or totally mixed. In addition, because of the seasonal (usually Spring) intrusion of stratospheric air into the troposphere at mid-latitudes, there is an additional contribution of stratospheric 10Be. Therefore, the measured 10Be concentration (or flux) in polar ice is modulated not only by production but also by climate/precipitation effects (e.g., Steig et al., 1996; Bard et al., 1997). This led Lal (1987) to the extreme conclusion that variations of polar 10Be reflect a meteorological, rather than solar, signal. However, comparison between Greenland and Antarctic 10Be series and between 10Be and 14C data (e.g., Bard et al., 1997; Horiuchi et al., 2008; Beer et al., 2012) suggests that the beryllium data mostly depicts production variations (i.e., solar signal) on top of which some meteorological effects can be superposed (see also Section 3.7.3).

Since both assumptions of the global and purely-local polar production of 10Be archived in polar ice are over-simplified, several attempts have been made to overcome this problem. For instance, McCracken (2004) proposed several simple mathematical models of partial atmospheric mixing (without division in the troposphere and stratosphere) and compared them with observed data. From this semi-empirical approach McCracken concluded that M2 (full mixing above 60° latitude and a limited mixing between 40° and 60° latitude) is a reasonable model for Antarctica. Vonmoos et al. (2006) assumed that the production of 10Be recorded in Greenland is related to the entire hemisphere in the stratosphere (i.e, global stratospheric mixing) but is limited to latitudes above 40° latitude in the troposphere (partial tropospheric mixing). This approach uses either semi-empirical or indirect arguments in choosing the unknown degree of mixing.

Recent efforts in employing modern atmospheric 3D circulation models for simulations of 10Be transport and deposition, including realistic air-mass transport and dry-vs-wet deposition (Field et al., 2006; Heikkila et al., 2008, 2009), look more promising. An example of 10Be deposition computed on the world grid using the NASA GISS model (Field et al., 2006) is shown in Figure 12. Precision of the models allows one to distinguish local effects, e.g., for Greenland (Heikkila et al., 2008). A simulation performed by combining a detailed 10Be-production model with an air-dynamics model can result in an absolute model relating production and deposition of the radionuclide. We may expect this breakthrough to occur in the near future. The validity and usefulness of this approach has been recently demonstrated by Usoskin et al. (2009a), who directly modeled production (using the CRAC model - Usoskin and Kovaltsov, 2008a) and transport (using the GISS ModelE - Koch et al., 2006) of a short-living beryllium isotope 7Be and showed that such a combined model is able to correctly reproduce both the absolute level and temporal variations of the 7Be concentration measured in near ground air around the globe. Keeping in mind the similarity between production and transport of the two beryllium isotopes, 7Be and 10Be, this serves as support for the advanced modelling of 10Be transport. A similar general agreement between measured and modelled seasonal variability has been recently found for 10Be in an Antarctic ice core (Pedro et al., 2011).

3.3.4 Effect of the geomagnetic field

In order to properly account for geomagnetic changes (Section 3.1.2), one needs to know the effective region in which the radionuclide is produced before being stored in the archive analyzed. For instance, if the concentration of 10Be measured in polar ice reflects mainly the isotope's production in the polar atmosphere (as, e.g., assumed by Usoskin et al., 2003c), no strong geomagnetic signal is expected to be observed, since the geographical poles are mostly related to high geomagnetic latitudes. On the other hand, assuming global mixing of atmospheric 10Be before deposition in polar ice (e.g., Masarik and Beer, 1999), one expects that only changes in the geomagnetic dipole moment affect will the signal. However, because of partial mixing, which can be different in the stratosphere and troposphere, taking into account migration and displacement of the geomagnetic dipole axis may be essential for a reliable reconstruction of solar variability from 10Be data

(a) Annual mean wet 10Be deposition (10 27 kg/m2/s)

660 1320 1980 2640 3300 3960 4620 5280 5940

(b) Annual mean dry 10Be deposition (10~27 kg/m2/s)

660 1320 1980 2640 3300 3960 4620 5280 5940

Figure 12: Wet (panel a) and dry (panel b) deposition of 10Be, computed using the NASA GISS model (Field et al., 2006) for a fixed sea-surface temperature.

(McCracken, 2004). Therefore, only a full combination of the transport and production models, the latter explicitly including geomagnetic effects estimated from paleomagnetic reconstructions, can adequately account for geomagnetic changes and separate the solar signal. These will form the next generation of physics-based models for the cosmogenic-isotope proxy method. We note that paleomagnetic data should ideally not only provide the dipole moment (VADM or VDM) but should also provide estimates of the geomagnetic axis attitude and displacement of the dipole center (Korte et al., 2011).

3.4 Other potential proxy

An interesting new potential proxy for solar activity (or cosmic ray) variability on the long-term centennial-to-millennial time scale has been proposed recently by Traversi et al. (2012). This is the nitrate content in a polar ice core Talos Dome in Antarctica, which has a favorite location in the sense of snow accumulation and conservation of such volatile specie as nitrate. Nitrate-related species are partly produced in the stratosphere/troposphere as a result of the ionization of the atmospheric air by cosmic rays and, partly, via terrestrial sources (e.g., lightnings) and are subject to air transport (Rozanov et al., 2012). As shown by Traversi et al. (2012), the nitrate concentration/flux measured in the Talos Dome ice core for the Holocene period agrees well with the cosmogenic data of 14C in tree rings and 10Be in both Antarctic and Greenland ice cores, on the time scales from centennia to millennia. Due to the large errors of the ice core dating, 200-300 years (ScMipbach et al., 2011), shorter time scales cannot be considered. The level of the nitrate variability is generally consistent with that predicted by theoretical models assuming its production by GCR in the atmosphere (Semeniuk et al., 2011; Rozanov et al., 2012). Thus, the nitrate in an ice core provides a potential to become a new proxy of long-term solar activity, with independent atmospheric fate, which would strengthen the robustness of the reconstructions. However, an independent confirmation of the result and a more detailed model are needed before it can serve as a new quantitative proxy. Note that the mechanism of the nitrate production and transport is not related to the possible nitrate peaks claimed to be caused by strong solar energetic-particle events (see Section 5.3).

3.5 Towards a quantitative physical model

Several methods have been developed historically to convert measured cosmogenic-isotope data into a solar activity index, ranging from very simple regressions to physics-based models. A new step in long-term solar-activity reconstruction has been made recently, which is the development of the proxy method in which physics-based models are used, instead of a phenomenological regression, to link SN with cosmogenic-isotope production (Usoskin et al., 2003c, 2007; Solanki et al., 2004; Vonmoos et al., 2006; Muscheler et al., 2007; Steinhilber et al., 2008). Due to recent theoretical developments, it is now possible to construct a chain of physical models to model the entire relationship between solar activity and cosmogenic data.

The physics-based reconstruction of solar activity (in terms of sunspot numbers) from cosmo-genic proxy data includes several steps:

• Computation of the isotope's production rate in the atmosphere from the measured concentration in the archive (Sections 3.2.2 and 3.3.2);

• Computation, considering independently-known secular geomagnetic changes (see Section 3.2.5) and a model of the CR-induced atmospheric cascade, of the GCR spectrum parameter quantified via the modulation potential ^ (Section 3.5.2), some reconstructions being terminated at this point;

• Computation of a heliospheric index, whether of the open solar magnetic flux or of the average HMF intensity at the Earth's orbit (Section 3.5.2)

• Computation of a solar index (sunspot number series), corresponding to the above-derived heliospheric parameter (Section 3.5.3).

Presently, all these steps can be completed using appropriate models. Some models stop after computations of the modulation potential as its translation into the solar index may include additional uncertainties. Although the uncertainties of the models may be considerable, the models allow a full basic quantitative reconstruction of solar activity in the past. However, much needs to be done, both theoretically and experimentally, to obtain an improved reconstruction.

3.5.1 Regression models

Mathematical regression is the most apparent and often used (even recently) method of solar-activity reconstruction from proxy data (see, e.g., Stuiver and Quay, 1980; Ogurtsov, 2004). The reconstruction of solar activity is performed in two consecutive steps. First, a phenomenological regression (either linear or nonlinear) is built between a proxy data set and a direct solar-activity index for the available "training" period (e.g., since 1750 for WSN or since 1610 for GSN). Then this regression is extrapolated backwards to evaluate SN from the proxy data. The main shortcoming of the regression method is that it depends on the time resolution and choice of the "training" period. The former is illustrated by Figure 13, which shows the scatter plot of the 10Be concentration vs. GSN for the annual and 11-year smoothed data. One can see that the slope of the 10Be-vs-GSN

Be (104 at/g)

X X X O 11-y — —Regression, annual Regression, 11y b)

1 1 N^ O

1 10Be (104 at/g) 15

Figure 13: Scatter plot of smoothed group sunspot numbers vs. (2-year delayed) 10Be concentration. a) Annual (connected small dots) and 11-year averaged (big open dots) values. b) Best-fit linear regressions between the annual (dashed line) and 11-year averaged values (solid line). The dots are the same as in panel (a). (After Usoskin and Kovaltsov, 2004).

relation (about -500 g/atom) within individual cycles is significantly different from the slope of the long-term relation (about -100 g/atom), i.e., individual cycles do not lie on the line of the 11-year averaged cycles. Moreover, the slope of the regression for individual 11-year cycles varies essentially depending on the solar activity level. Therefore, a formal regression built using the annual data for 1610-1985 yields a much stronger GSN-vs-10Be dependence than for the cycle-averaged data (see Figure 13b), leading to a potentially-erroneous evaluation of the sunspot number from the 10 Be proxy data.

It is equally dangerous to evaluate other solar/heliospheric/terrestrial indices from sunspot numbers, by extrapolating an empirical relation obtained for the last few decades back in time.

This is because the last decades (after the 1950s), which are well covered by direct observations of solar, terrestrial and heliospheric parameters, correspond to a very high level of solar activity. After a steep rise in activity level between the late 19th and mid 20th centuries, the activity remained at a roughly constant high level, being totally dominated by the 11-year cycle without a long-term trends. Accordingly, all empirical relations built based on data for this period are focused on the 11-year variability and can overlook possible long-term trends (Mursula et al., 2003). This may affect all regression-based reconstructions, whose results cannot be independently (directly or indirectly) tested. In particular, this may be related to solar irradiance reconstructions, which are often based on regression-like models, built and verified using data for the last three solar cycles, when there was no strong trend in solar activity.

£ -5 Í -10 =S -15 g -20

1650 1700 1750 1800 1850 1900 1950 years

Figure 14: An unsuccessful attempt of the reconstruction of cosmic-ray intensity in the past using a regression with sunspot numbers. Dots represent the observed cosmic-ray intensity since 1951. Note the absence of a long-term trend. Image reproduced by permission from Belov et al. (2006), copyright by Elsevier.

As an example let us consider an attempt (Belov et al., 2006) to reconstruct cosmic-ray intensity since 1610 from sunspot numbers using a (nonlinear) regression. The regression between the count rate of a neutron monitor and sunspot numbers, established for the last 30 years, yields an agreement at a 95% confidence level for the period 1976-2003 (see Figure 3). Based on that, Belov et al. (2006) extrapolated the regression back in time to produce a reconstruction of cosmic-ray intensity (quantified in NM count rate) to 1560 (see Figure 14). One can see that there is no notable long-term trend in the reconstruction, and the fact that all CR maxima essentially lie at the same level, from the Maunder minimum to modern times, is noteworthy. It would be difficult to dispute such a result if there was no direct test for CR levels in the past. Independent reconstructions based on cosmogenic isotopes or theoretical considerations (e.g., Usoskin et al., 2002b; Scherer et al., 2004; Scherer and Fichtner, 2004) provide clear evidence that cosmic-ray intensity was essentially higher during the Maunder minimum than nowadays. This example shows how easy it is to overlook an essential feature in a reconstruction based on a regression extrapolated far beyond the period it is based on. Fortunately, for this particular case we do have independent information that can prevent us from making big errors. In many other cases, however, such information does not exist (e.g., for total or spectral solar irradiance), and those who make such unverifiable reconstructions should be careful about the validity of their models beyond the range of the established relations.

3.5.2 Reconstruction of heliospheric parameters

The modulation potential ^ (see Section 3.1.1) is directly related to cosmogenic isotope production in the atmosphere. It is a parameter describing the spectrum of galactic cosmic rays (see the definition and full description of this index in Usoskin et al., 2005) and is sometimes used as a

stand-alone index of solar (or, actually, heliospheric) activity. We note that, provided the isotope production rate Q is estimated and geomagnetic changes can be properly accounted for, it is straightforward to obtain a time series of the modulation potential, using, e.g., the relation shown in Figure 9. Several reconstructions of modulation potential for the last few centuries are shown in Figure 15. While being quite consistent in the relative changes, they differ in the absolute level and fine details. Reconstructions of solar activity often end at this point, representing solar activity by the modulation potential, as some authors (e.g., Beer et. al, 2003; Vonmoos et. al., 2006; Muscheler et al., 2007) believe that further steps (see Section 3.5.3) may introduce additional uncertainties. However, since 4> is a heliospheric, rather than solar, index, the same uncertainties remain when using it as an index of solar activity. Moreover, the modulation potential is a model-dependent quantity (see discussion in Section 3.1.1) and therefore does not provide an unambiguous measure of heliospheric activity. In addition, the modulation potential is not a physical index but rather a formal fitting parameter to describe the GCR spectrum near Earth (Usoskin et. al., 2005) and, thus, is not a universal solar-activity index.

Figure 15: Several reconstructions of the decade-averaged modulation potential for the last few centuries: from sunspot numbers (SN(U02) - Usoskin et al., 2002b), from 14C data (14C(S04), 14C(M05), 14C(M07) - Solanki et al., 2004; Muscheler et al., 2005, 2007, respectively), from Antarctic 10Be data. (10Be(U03), 10Be(MC04) - Usoskin et al., 2003c; McCracken et al., 2004, respectively). The thick black NM curve is based on direct cosmic-ray measurements by neutron monitors since 1951 (Usoskin et al, 2011) and ionization chambers since 1936 (McCracken and Beer, 2007).

Modulation of GCR in the heliosphere (see Section 3.1.1) is mostly defined by the turbulent heliospheric magnetic field (HMF), which ultimately originates from the sun and is thus related to solar activity. It has been shown, using a theoretical model of the heliospheric transport of cosmic

rays (e.g., Usoskin et al., 2002b), that on the long-term scale (beyond the 11-year solar cycle) the modulation potential ^ is closely related to the open solar magnetic flux Fo, which is a physical quantity describing the solar magnetic variability (e.g., Solanki et al., 2000; Krivova et al., 2007).

Sometimes, instead of the open magnetic flux, the mean HMF intensity at Earth orbit, B, is used as a heliospheric index (Caballero-Lopez and Moraal, 2004; McCracken, 2007; Steinhilber et al., 2010). Note that B is linearly related to Fo assuming constant solar-wind speed, which is valid on long-term scales. An example of HMF reconstruction for the last 600 years is shown in Figure 16. In addition, the count rate of a "pseudo" neutron monitor (i.e., a count rate of a neutron monitor if it was operated in the past) is considered as a solar/heliospheric index (e.g., Beer, 2000; McCracken and Beer, 2007).

Figure 16: An example of reconstruction of the heliospheric magnetic field at Earth orbit for the last 600 years from 10Be data. Image reproduced by permission from McCracken (2007), copyright by AGU.

3.5.3 A link to sunspot numbers

The open solar magnetic flux Fo described above is related to the solar surface magnetic phenomena such as sunspots or faculae. Modern physics-based models allow one to calculate the open solar magnetic flux from data of solar observation, in particular sunspots (Solanki et al., 2000, 2002; Krivova et al., 2007). Besides the solar active regions, the model includes ephemeral regions. Although this model is based on physical principals, it contains some unknowns like the decay time of the open flux, which cannot be measured or theoretically calculated and has to be found by means of fitting the model to data. This free parameter has been determined by requiring the model output to reproduce the best available data sets for the last 30 years with the help of a genetic algorithm. Inversion of the model, i.e., the computation of sunspot numbers for given Fo values is formally a straightforward solution of a system of linear differential equations, however, the presence of noise in the real data makes it only possible in a numerical-statistical way (see, e.g., Usoskin et al., 2004, 2007). By inverting this model one can compute the sunspot-number series corresponding to the reconstructed open flux, thus forging the final link in a chain quantitatively

connecting solar activity to the measured cosmogenic isotope abundance. A sunspot-number series reconstructed for the Holocene using 14C isotope data is shown in Figure 17. While the definition of the grand minima (Section 4.2) is virtually insensitive to the uncertainties of paleomagnetic data, the definition of grand maxima depends on the paleomagnetic model used (Usoskin et al., 2007). Since the Y00 paleomagnetic model forms an upper bound for the true geomagnetic strength (Section 3.1.2), the corresponding solar-activity reconstructions may underestimate the solar-activity level. Accordingly, the grand maxima defined using the Y00 model are robust and can be regarded as "maximum maximorum" (see Section 4.3).

0-8000 -6000 -4000 -2000 0 2000

Years (-BC/AD)

Figure 17: Long-term sunspot-number reconstruction from 14C data (after Usoskin et al., 2007). All data are decade averages. Solid (denoted as 'Y00') and grey ('K05') curves are based on the paleo-geomagnetic reconstructions of Yang et al. (2000) and Korte and Constable (2005), respectively. Observed group sunspot numbers (Hoyt and Schatten, 1998) are shown after 1610.

As very important for the climate research, the variations of the total solar irradiance (TSI) are sometimes reconstructed from the solar proxy data (Steinhilber et al., 2009; Vieira et al., 2011). However, the absolute range of the TSI variability on the centennial-millennial time scales still remains unknown (Schmidt et al., 2012).

3.6 Solar activity reconstructions

Detailed computational models of cosmogenic isotope production in the atmosphere (e.g., Masarik and Beer, 1999) have opened up a new possibility for long-term solar-activity reconstruction (e.g., Beer, 2000). The first quantitative reconstructions of solar activity from cosmogenic proxy appeared in the early 2000s based on 10Be deposited in polar ice (Beer et al., 2003; Usoskin et al., 2003c).

Beer et al. (2003) reconstructed the modulation potential on a multi-millennial timescale using the model computations by Masarik and Beer (1999) and the 10Be data from the GISP2 core in Greenland. This result has been extended, even including the 14C data set, and presently covers the whole Holocene (Vonmoos et al., 2006; Steinhilber et al., 2008, 2010). Usoskin et al. (2003c) presented the reconstruction of sunspot activity over the last millennium, based on 10Be data from both Greenland and Antarctica, using a physics-based model described in detail in Usoskin et al. (2004). This result reproduces the four known grand minima of solar activity - Maunder, Sporer, Wolf and Oort minima (see Section 4.2). Later Solanki et al. (2004) reconstructed 10-year-averaged sunspot numbers from the 14 C content in tree rings throughout the Holocene and estimated its uncertainties (see Figure 18). This result was disputed by Muscheler et al. (2005), whose concurrent model, however, rested on an erroneous normalization, as argued in Solanki et al.

(2005). The reconstruction of Solanki et al. (2004) has been recently updated by Usoskin et al. (2006a), using a newer paleomagnetic reconstruction by Korte and Constable (2005), and was later slightly revised (Usoskin et al., 2007), considering an updated model of the solar open magnetic flux by Krivova et al. (2007). Reconstruction of the HMF from 10Be data has been performed by Caballero-Lopez and Moraal (2004), using a model of CR modulation in the heliosphere and a

10Be production model by Webber and Higbie (2003). Recently it was revised (McCracken, 2007) to present a detailed reconstruction of HMF intensity since 1428. The most recent reconstruction of the heliospheric modulation potential was done by Steinhilber et al. (2012) using the combined principal component analysis of several data sets.

The obtained results are discussed in Section 4.

Figure 18: 10-year averaged sunspot numbers: Actual group sunspot numbers (thick grey line) and the reconstructions based on 10Be (thin curve, Usoskin et al., 2003c) and on 14C (thick curve with error bars, Solanki et al., 2004). The horizontal dotted line depicts the high activity threshold.

3.7 Verification of reconstructions

Because of the diversity of the methods and results of solar-activity reconstruction, it is vitally important to verify them. Even though a full verification is not possible, there are different means of indirect or partial verification, as discussed below. Several solar-activity reconstructions on the millennium timescale, which differ from each other to some degree and are based on terrestrial cosmogenic isotope data, have been published recently by various groups. Also, they may suffer from systematic effects. Therefore, there is a need for an independent method to verify/calibrate these results in order to provide a reliable quantitative estimate of the level of solar activity in the past, prior to the era of direct observations.

3.7.1 Comparison with direct data

The most direct verification of solar-activity reconstruction is a comparison with the actual GSN sunspot data for the last few centuries. However, regression-based models (see Section 3.5.1) cannot be tested in this way, since it would require a long set of independent direct data outside the "training" interval. It is usual to include all available data into the "training" period to increase the statistics of the regression, which rules out the possibility of testing the model. On the other hand, such a comparison to the actual GSN since 1610 can be regarded as a direct test for a physics-based model since it does not include phenomenological links over the same

time interval. The period of the last four centuries is pretty good for testing purposes since it includes the whole range of solar activity levels from the nearly spotless Maunder minimum to the modern period of a very active sun. As an example, a comparison between the observed GSN and the 14C-based (Solanki et al., 2004) and 10Be-based (Usoskin et al., 2003c) reconstructions is shown in Figure 18. The agreement between the actual and reconstructed sunspot numbers is quite good, the correlation coefficient for the 14C-based series is r = 0.93 with the RMS deviation between the two series being six for the period of 1610-1900 (Solanki et al., 2004). We want to stress that this reconstruction is fully physics based and does not include any fitting to the whole GSN data series; thus this comparison verifies the approach in both absolute level and relative variations. The agreement between GSN and 10Be-based reconstructions (Usoskin et al., 2003c) is also good (r = 0.78, RMS = 10 for 1700-1985). In this case, however, the comparison can only test the relative variation because of the unknown proportionality coefficient between the measured concentration of 10Be and the production rate (Section 3.3.3), which is fitted to match the overall level of the reconstructed solar activity. One can see that the reconstructed sunspot series generally follows the real GSN series, depicting the same main features, namely, the Maunder minimum, the tiny Dalton minimum, a slight decrease of activity around 1900 (sometimes called the modern minimum) as well as a steep rise in the first half of the 20th century. This validates the reliability of the physics-based reconstruction of sunspot numbers. Note, however, that individual 11-year cycles are poorly resolved in these reconstructions.

Models focused on the reconstruction of heliospheric parameters (HMF or the modulation potential cannot be verified in this manner since no heliospheric data exists before the middle of the 20th century. Comparison to direct cosmic-ray data after the 1950s (or, with caveats, after the 1930s - McCracken and Beer, 2007) is less conclusive, since the latter are of shorter length and correspond to a period of high solar activity, leading to larger uncertainties during grand minima.

It is important that some (semi)empirical relations forming the basis for the proxy method are established for the recent decades of high solar activity. The end of the Modern grand maximum of activity and the current low level of activity, characterized by the highest ever observed cosmic ray flux as recorded by ground-based neutron monitors, the very low level of the HMF and geomagnetic activity, should help to verify the connections between solar activity, cosmic ray fluxes, geomagnetic activity, the heliospheric magnetic field, and open field. Since some of these connections are somewhat controversial, these extreme conditions should help to quantify them better.

3.7.2 Meteorites and lunar rocks: A direct probe of the galactic cosmic-ray flux

Another more-or-less direct test of solar/heliospheric activity in the past comes from cosmogenic isotopes measured in lunar rock or meteorites. Cosmogenic isotopes, produced in meteoritic or lunar rocks during their exposure to CR in interplanetary space, provide a direct measure of cosmic-ray flux. Uncertainties due to imprecisely known terrestrial processes, including the geomagnetic shielding and redistribution process, are naturally avoided in this case, since the nuclides are directly produced by cosmic rays in the body of the rock, where they remain until they are measured, without any transport or redistribution. The activity of a cosmogenic isotope in meteorite/lunar rock represents an integral of the balance between the isotope's production and decay, thus representing the time-integrated CR flux over a period determined by the mean life of the radioisotope. The results of different analyses of measurements of cosmogenic isotopes in meteoritic and lunar rocks show that the average GCR flux remained roughly constant - within 10% over the last million years and within a factor of 1.5 for longer periods of up to 109 years (e.g., Vogt et al., 1990; Grieder, 2001).

By means of measuring the abundance of relatively short-lived cosmogenic isotopes in meteorites, which fell through the ages, one can evaluate the variability of the CR flux, since the production of cosmogenic isotopes ceases after the fall of the meteorite. A nearly ideal isotope

for studying centurial-scale variability is 44Ti with a half-life of 59.2 ± 0.6 yr (a lifetime of about 85 years). The isotope is produced in nuclear interactions of energetic CR with nuclei of iron and nickel in the body of a meteorite (Bonino et al., 1995; Taricco et al., 2006). Because of its mean life, 44Ti is relatively insensitive to variations in cosmic-ray flux on decade (11-year Schwabe cycle) or shorter timescales, but is very sensitive to the level of CR flux and its variations on a centurial scale. Using a full model of 44Ti production in a stony meteorite (Michel and Neumann, 1998) and data on the measured activity of cosmogenic isotope 44 Ti in meteorites, which fell during the past 235 years (Taricco et al., 2006), Usoskin et al. (2006c) tested, in a straightforward manner, several recent reconstructions of heliospheric activity after the Maunder minimum. First, the expected 44Ti activity has been calculated from the reconstructed series of the modulation potential, and then compared with the results of actual measurements (see Figure 19). It has been shown that 44Ti data can distinguish between various reconstructions of past solar activity, allowing unrealistic models to be ruled out. Since the life-time of the 44Ti is much longer than the 11-year cycle, this method does not allow for the reconstruction of solar/heliospheric activity, but it serves as a direct way to test existing reconstructions independently. Most of the reconstructions appear consistent with the measured 44Ti activity in meteorites, including the last decades, thus validating their veracity. The only apparently-inconsistent model is the one by Muscheler et al. (2005), which is based on erroneous normalization (as discussed in Solanki et al., 2005). In particular, the 44Ti data confirms significant secular variations of the solar magnetic flux during the last century (cf. Lockwood et al., 1999; Solanki et al., 2000; Wang et al., 2005).

Figure 19: Immediate 44Ti activity in stony meteorites as a function of time of fall. Dots with error bars correspond to measured values (Taricco et al., 2006). Curves correspond to the theoretically expected 44Ti activity, computed using the method of Usoskin et al. (2006c) and different reconstructions of 0 shown in Figure 15.

3.7.3 Comparison between isotopes

As an indirect test of the solar-activity reconstruction, one can compare different isotopes. The idea behind this test is that two isotopes, 14C and 10Be, have essentially different terrestrial fates, so that only the production signal, namely, solar modulation of cosmic rays, can be regarded as common in the two series. Processes of transport/deposition are completely different (moreover, the 14C series is obtained as an average of the world-wide-distributed samples). The effect of changing geomagnetic fields is also different (although not completely) for the two isotopes, since radiocarbon is globally mixed, while 10Be is only partly mixed before being stored in an archive. Even comparison between data of the same 10Be isotope, but measured in far-spaced ice cores (e.g., Greenland and Antarctica), may help in separating climatic and extraterrestrial factors, since meteorology in the two opposite polar areas is quite different.

The first thorough consistent comparison between 10Be and 14C records for the last millennium was performed by Bard et al. (1997). They assumed that the measured 10Be concentration in Antarctica is directly related to CR variations. Accordingly, 14C production was considered as proportional to 10Be data. Then, applying a 12-box carbon-cycle model, Bard et al. (1997) computed the expected A14C synthetic record. Finally, these 10Be-based A14C variations were compared with the actual measurements of A14C in tree rings, which depicted a close agreement in the profile of temporal variation (coefficient of linear correlation r = 0.81 with exact phasing). Despite some fine discrepancies, which can indicate periods of climatic influence in either (or both) of the series, that result has clearly proven the dominance of solar modulation of cosmogenic nu-clide production variations during the last millennium. This conclusion has been confirmed (e.g., Usoskin et al., 2003c; Muscheler et al., 2007) in the sense that quantitative solar-activity reconstructions, based on 10Be and 14C data series for the last millennium, yield very similar results, which differ only in small details. However, a longer comparison over the entire Holocene timescale suggests that, while centennial variations of solar activity reconstructed from the two isotopes are very close to each other, there might be a discrepancy in the very long-term trend (Vonmoos et al., 2006; Muscheler et al., 2007), whose nature is not clear (climate changes, geomagnetic effects or model uncertainties).

Recently, Usoskin et al. (2009b) studied the dominance of the solar signal in different cosmogenic isotope data on different time scales. They compared the expected 10Be variations computed from 14C-based reconstruction of cosmic ray intensity with the actually measured 10Be abundance at the sites and found that: (1) There is good agreement between the 14C and 10Be data sets, on different timescales and at different locations, confirming the existence of a common solar signal in both isotope data; (2) The 10Be data are driven by the solar signal on timescales from about centennial to millennial time scales; (3) The synchronization is lost on short (< 100 years) timescales, either due to local climate or chronological uncertainties (Delaygue and Bard, 2011) but the solar signal becomes important even at short scales during periods of Grand minima of solar activity, (4) There is an indication of a possible systematic uncertainty in the early Holocene (cf., Vonmoos et al., 2006), likely due to a not-perfectly-stable thermohaline circulation. Overall, both 14C- and 10Be-based records are consistent with each other over a wide range of timescales and time intervals.

Thus, comparison of the results obtained from different sources implies that the variations of cosmogenic nuclides on the long-term scale (centuries to millennia) during the Holocene are primarily defined by the solar modulation of CR.

3.8 Composite reconstruction

Most of the earlier solar activity reconstructions are based on single proxy records, either 14 C or 10Be. Although they are dominated by the same production signal, viz. solar activity, (see Section 3.7.3), they still contain essential fractions of noise.

A promising first step in the direction of extracting the common solar signal from different proxy records was made recently by Steinhilber et al. (2012). They combined, in a composite reconstruction, different 10Be ice core records from Greenland and Antarctica with the global 14C tree ring record. The composite was made in a mathematical way, using the principal component analysis as a numerical tool. This analysis formally finds the common variability in different series, that is assumed to be the solar signal. The resultant reconstruction is shown in Figure 20 and is the most consistent reconstruction of the cosmic ray (and thus solar magnetic activity) variability over the Holocene, available up-to-date. However, since the used mathematical tool can only work with the relative variability, the reconstruction also yields the relative values rather than absolute values, and it is not available in the terms of sunspot numbers.

9000 6000 7000 6000 5000 4000 3000 2000 1000 0

Time (years BP)

Figure 20: Reconstruction of the solar modulation potential over the Holocene based on a composite record (Section 3.8), along with la uncertainties. Time is given in year BP. Image reproduced by permission from Steinhilber et al. (2012).

A full physics-based multi-proxy composite reconstruction of the solar activity on the millennial time scale is still pending.

3.9 Summary

In this section, a proxy method of past-solar-activity reconstruction is described in detail.

This method is based on the use of indirect proxies of solar activity, i.e., quantitative parameters, which can be measured now, but represent signatures, stored in natural archives, of the different effects of solar magnetic activity in the past. Such traceable signatures can be related to nuclear or chemical effects caused by cosmic rays in the Earth's atmosphere, lunar rocks or meteorites. This approach allows one to obtain homogeneous data sets with stable quality and to improve the quality of data when new measurement techniques become available. It provides the only possible regular indicator of solar activity on a very long-term scale.

The most common proxy of solar activity is formed by data of the cosmogenic radionuclides, 10Be and 14C, produced by cosmic rays in the Earth's atmosphere. After a complicated transport in the atmosphere, these cosmogenic isotopes are stored in natural archives such as polar ice, trees, marine sediments, from where they can now be measured. This process is also affected by changes in the geomagnetic field and the climate.

Radioisotope 14C, measured in independently-dated tree rings, forms a very useful proxy for long-term solar-activity variability. It participates in the complicated carbon cycle, which smoothes out spatial and short-term variability of isotope production. For the Holocene period, with its stable climate, it provides a useful tool for studying solar activity in the past. Existing models allow the quantitative conversion between the measured relative abundance of 14 C and the production rate in the atmosphere. The use of radiocarbon for earlier periods, the glacial and deglaciation epochs, is limited by severe climate and ocean ventilation changes. Radiocarbon data cannot be used after the end of the 19th century because of the Suess effect and atmospheric nuclear tests.

Another solar activity proxy is the cosmogenic 10Be isotope measured in stratified polar ice cores. Atmospheric transport of 10Be is relatively straightforward, but its details are as of yet unresolved, leading to the lack of a reliable quantitative model relating the measured isotope concentration in ice to the atmospheric production. Presently, it is common to assume that the production rate is proportional, with an unknown coefficient, to the measured concentration. However, a newly-developed generation of models, which include 3D atmospheric-circulation models, will hopefully solve this problem soon.

Recently a new proxy, nitrate concentration measured in an Antarctic ice core, has been proposed for long-term solar activity reconstructions, but it still needs verification and model support.

Modern physics-based models make it possible to build a chain, which quantitatively connects isotope production rate and sunspot activity, including subsequently the GCR flux quantified via the modulation potential, the heliospheric index, quantified via the open solar magnetic flux or the average HMF intensity at the Earth's orbit, and finally the sunspot-number series. Presently, all these steps can be made using appropriate models allowing for a full basic quantitative reconstruction of solar activity in the past. The main uncertainties in the solar-activity reconstruction arise from paleo-magnetic models and the overall normalization.

An independent verification of the reconstructions, including direct comparison with sunspot numbers, cosmogenic isotopes in meteorites and the comparison of different models with each other, confirms their veracity in both relative variations and absolute level. It also implies that the variations in cosmogenic nuclides on the long-term scale (centuries to millennia) during the Holocene are primarily defined by the solar modulation of CR.

4 Variability of Solar Activity Over Millennia

Several reconstructions of solar activity on multi-millennial timescales have been performed recently using physics-based models (see Section 3) from measurements of 14C in tree rings and 10Be in polar ice. The validity of these models for the last few centuries was discussed in Section 3.7. In this section we discuss the temporal variability of thus-reconstructed solar activity on a longer scale.

Here we consider the 14C-based decade reconstruction of sunspot numbers (shown in Figure 21). It is identical to that shown in Figure 17, but includes also a Gleissberg 1-2-2-2-1 filter in order to suppress noise and short-term fluctuations. This series forms the basis for the forthcoming analysis, while differences related to the use of other reconstructions are discussed.

Years (-BC/AD)

Figure 21: Sunspot activity (over decades, smoothed with a 12221 filter) throughout the Holocene, reconstructed from 14C by Usoskin et al. (2007) using geomagnetic data by Yang et al. (2000). Blue and red areas denote grand minima and maxima, respectively.

4.1 Quasi-periodicities and characteristic times

In order to discuss spectral features of long-term solar-activity dynamics, we show in Figure 22 a wavelet spectral decomposition of the sunspot number reconstruction throughout the Holocene shown in Figure 21. The left-hand panels show the conventional wavelet decomposition in the time-frequency domain, while the right-hand panels depict the global spectrum, namely, an integral over the time domain, which is comparable to a Fourier spectrum. The peak in the global spectrum at about an 80-year period corresponds to the Gleissberg periodicity, known from a simple Fourier analysis of the A14C series (Peristykh and Damon, 2003). The peak at an approximately 150 year period does not correspond to a persistent periodicity, but is formed by a few time intervals (mostly 6000-4000 BC) and can be related to another "branch" of the secular cycle, according to

Ogurtsov et al. (2002). The de Vries/Suess cycle, with a period of about 210 years, is prominent in the global spectrum, but it is intermittent and tends to become strong with around 2400 clustering time (Usoskin and Kovaltsov, 2004). Another variation with a period of around 350 years can be observed after 6000 BC (cf. Steinhilber et al., 2012). Variations with a characteristic time of 600-700 years are intermittent and can be hardly regarded as a typical feature of solar activity. There is also a weak millennial periodicity (Steinhilber et al., 2012). Of special interest is the 2000-2400 year Hallstatt cycle (see, e.g., Vitinsky et al., 1986; Damon and Sonett, 1991; Vasiliev and Dergachev, 2002), which is relatively stable and mostly manifests itself as a modulation of long-term solar activity, leading to the clustering of grand minima (Usoskin et al., 2007).

On the other hand, an analysis of the occurrence of grand minima (see Section 4.2) shows no clear periodicity except for a marginal 2400 year clustering, implying that the occurrence of grand minima and maxima is not a result of long-term cyclic variability but is defined by stochastic/chaotic processes.

5000 4000

500 400

300 200


-6000 Time (year)

1000 Power

-10000 -8000 -6000 -4000 Time (year)

Figure 22: Wavelet (Morlet basis) spectrum of the sunspot-number reconstruction shown in Figure 21. Left and right-hand panels depict 2D and global wavelet spectra, respectively. Upper and lower panels correspond to period ranges of 500-5000 years and 80-500 years, respectively. Dark/light shading denotes high/low power. Hatched areas depict the cone of influence where the result is not fully reliable because of the proximity of the edges of the time series.

4.2 Grand minima of solar activity

A very particular type of solar activity is the grand minimum, when solar activity is greatly reduced. The most famous is the Maunder minimum in the late 17th century, which is discussed below in some detail (for a detailed review see the book by Soon and Yaskell, 2003). Grand minima are believed to correspond to a special state of the dynamo (Sokoloff, 2004; Miyahara et al., 2006b), and its very existence poses a challenge for the solar-dynamo theory. It is noteworthy that dynamo models do not agree on how often such episodes occur in the sun's history and whether their appearance is regular or random. For example, the commonly used mean-field dynamo yields a fairly-regular 11-year cycle (Charbonneau, 2010), while dynamo models including a stochastic driver predict the intermittency of solar magnetic activity (Choudhuri, 1992; ScMissler et al., 1994; Schmitt et al., 1996; Ossendrijver, 2000; Weiss and Tobias, 2000; Mininni et al., 2001; Charbonneau, 2001). Most of the models predict purely random occurrence of the grand minima, without any intrinsic long-term memory (Moss et al., 2008). Although cosmogenic isotope data suggest the possible existence of such memory (Usoskin et al., 2007), statistics is not sufficient to distinguish between the two cases (Usoskin et al., 2009d).

4.2.1 The Maunder minimum

The Maunder minimum is a representative of grand minima in solar activity (e.g., Eddy, 1976), when sunspots have almost completely vanished from the solar surface, while the solar wind kept blowing, although at a reduced pace (Cliver et al., 1998; Usoskin et al., 2001a). There is some uncertainty in the definition of its duration; the "formal" duration is 1645-1715 (Eddy, 1976), while its deep phase with the absence of apparent sunspot cyclic activity is often considered as 1645-1700, with the low, but very clear, solar cycle of 1700-1712 being ascribed to a recovery or transition phase (Usoskin et al., 2000). The Maunder minimum was amazingly well covered (more than 95% of days) by direct sunspot observations (Hoyt and Schatten, 1996), especially in its late phase (Ribes and Nesme-Ribes, 1993). On the other hand, sunspots appeared rarely (during ~ 2% of the days) and seemingly sporadically, without an indication of the 11-year cycle (Usoskin and Mursula, 2003). This makes it almost impossible to apply standard methods of time-series analysis to sunspot data during the Maunder minimum (e.g., Frick et al., 1997)). Therefore, special methods such as the distribution of spotless days vs. days with sunspots (e.g., Harvey and White, 1999; Kovaltsov et al., 2004) or an analysis of sparsely-occurring events (Usoskin et al., 2000) should be applied in this case. Using these methods, Usoskin et al. (2001a) have shown that sunspot occurrence during the Maunder minimum was gathered into two large clusters (1652-1662 and 1672-1689), with the mass centers of these clusters being in 1658 and 1679-1680. Together with the sunspot maxima before (1640) and after (1705) the deep Maunder minimum, this implies a dominant 22-year periodicity in sunspot activity throughout the Maunder minimum (Mursula et al., 2001), with a subdominant 11-year cycle emerging towards the end of the Maunder minimum (Ribes and Nesme-Ribes, 1993; Mendoza, 1997; Usoskin et al., 2000) and becoming dominant again after 1700. Similar behavior of a dominant 22-year cycle and a weak subdominant Schwabe cycle during the Maunder minimum has been found in other indirect solar proxy data: auroral occurrence (Krivsky and Pejml, 1988; Schlamminger, 1990; Silverman, 1992) and 14C data (Stuiver and Braziunas, 1993; Kocharov et al., 1995; Peristykh and Damon, 1998; Miyahara et al., 2006b). This is in general agreement with the concept of "immersion" of 11-year cycles during the Maunder minimum (Vitinsky et al., 1986, and references therein). This concept means that full cycles cannot be resolved and sunspot activity only appears as pulses around cycle-maximum times.

An analysis of 10Be data (Beer et al., 1998) implied that the 11-year cycle was weak but fairly regular during the Maunder minimum, but its phase was inverted (Usoskin et al., 2001a). A recent theoretical study (Owens et al., 2012; Wang and Sheeley Jr, 2013) confirms that such a phase change between cosmic rays and solar activity can appear for very weak cycles.

The time behavior of sunspot activity during the Maunder minimum yielded the following general scenario (Vitinsky et al., 1986; Ribes and Nesme-Ribes, 1993; Sokoloff and Nesme-Ribes, 1994; Usoskin et al., 2000, 2001a; Miyahara et al., 2006b). Transition from the normal high activity to the deep minimum did not have any apparent precursor. On the other hand, newly recovered data suggest that the start of the Maunder minimum might had been not very sudden but via a regular cycle of reduced height (Vaquero et al., 2011). A 22-year cycle was dominant in sunspot occurrence during the deep minimum (1645-1700), with the subdominant 11-year cycle, which became visible only in the late phase of the Maunder minimum. There is an indication that the length of solar cycle may slightly extend during and already slightly before a grand minimum (Miyahara et al., 2004; Nagaya et al., 2012). The 11-year Schwabe cycle started dominating solar activity after 1700. Recovery of sunspot activity from the deep minimum to normal activity was gradual, passing through a period of nearly-linear amplification of the 11-year cycle. It is interesting to note that such a qualitative evolution of a grand minimum is consistent with predictions of the stochastically-forced return map (Charbonneau, 2001).

Although the Maunder minimum is the only one with available direct sunspot observations, its predecessor, the Sporer minimum from 1450-1550, is covered by precise bi-annual measurements of 14C (Miyahara et al., 2006a). An analysis of this data (Miyahara et al., 2006a,b) reveals a similar pattern with the dominant 22-year cycle and suppressed 11-year cycle, thus supporting the idea that the above general scenario may be typical for a grand minimum. A similar pattern has been recently also for an un-named grand minima in the 4-th century BC (Nagaya et al., 2012).

A very important feature of sunspot activity during the Maunder minimum was its strong north-south asymmetry, as sunspots were only observed in the southern solar hemisphere during the end of the Maunder minimum (Ribes and Nesme-Ribes, 1993; Sokoloff and Nesme-Ribes, 1994). This observational fact has led to intensive theoretical efforts to explain a significant asymmetry of the sun's surface magnetic field in the framework of the dynamo concept (see the review by Sokoloff, 2004, and references therein). Note that a recent discovery (Arlt, 2008, 2009) of the Staudacher's original drawings of sunspots in late 18th century shows that similarly asymmetric sunspot occurrence existed also in the beginning of the Dalton minimum in 1790s (Usoskin et al., 2009c). However, the northern hemisphere dominated at that period contrary to the situation during the Maunder minimum.

4.2.2 Grand minima on a multi-millennial timescale

The presence of grand minima in solar activity on the long-term scale has been mentioned numerously (e.g., Eddy, 1977a; Solanki et al., 2004), using the radioisotope 14C data in tree rings. For example, Eddy (1977b) identified major excursions in the detrended 14C record as grand minima and maxima of solar activity and presented a list of six grand minima and five grand maxima for the last 5000 years (see Table 1). Stuiver and Braziunas (1989) and Stuiver et al. (1991) also studied grand minima as systematic excesses of the high-pass filtered 14 C data and suggested that the minima are generally of two distinct types: short minima of duration 50-80 years (called Maunder-type) and longer minima collectively called Sporer-like minima. Using the same method of identifying grand minima as significant peaks in high-pass filtered A14C series, Voss et al. (1996) provided a list of 29 such events for the past 8000 years. A similar analysis of bumps in the 14C production rate was presented recently by Goslar (2003). However, such studies retained a qualitative element, since they are based on high-pass-filtered 14C data and thus implicitly assume that 14C variability can be divided into short-term solar variations and long-term changes attributed solely to the slowly-changing geomagnetic field. This method ignores any possible long-term changes in solar activity on timescales longer than 500 years (Voss et al., 1996). The modern approach, based on physics-based modelling (Section 3), allows for the quantitative reconstruction of the solar activity level in the past, and thus, for a more realistic definition of the periods of grand minima or maxima.

Table 1: Approximate dates (in -BC/AD) of grand minima in reconstructed solar activity.

No. center duration comment

1 1680 80 Maunder

2 1470 160 Sporer

3 1305 70 Wolf

4 1040 60 a, d)

5 685 70 b, d)

6 -360 60 a, b, c, d)

7 -765 90 a, b, c, d)

8 -1390 40 b, d)

9 -2860 60 a, c, d)

10 -3335 70 a, b, c, d)

11 -3500 40 a, b, c, d)

12 -3625 50 a, b, d)

13 -3940 60 a, c, d)

14 -4225 30 c, d)

15 -4325 50 a, c, d)

16 -5260 140 a, b, d)

17 -5460 60 c, d)

18 -5620 40 d)

19 -5710 20 c, d)

20 -5985 30 a, c, d)

21 -6215 30 c, d, e)

22 -6400 80 a, c, d, e)

23 -7035 50 a, c, d)

24 -7305 30 c, d)

25 -7515 150 a, c, d)

26 -8215 110 d)

27 -9165 150 d)

a) According to Stuiver and Quay (1980); Stuiver and Braziunas (1989).

b) According to Eddy (1977a,b).

c) According to Goslar (2003).

d) According to Usoskin et al. (2007).

e) Exact duration is uncertain.

A list of 27 grand minima, identified in the quantitative solar-activity reconstruction of the last 11,000 years, shown in Figure 21, is presented in Table 1 (after Usoskin et al., 2007). The cumulative duration of the grand minima is about 1900 years, indicating that the sun in its present evolutionary stage spends ~ 1/6 (17%) of its time in a quiet state, corresponding to grand minima. Note that the definition of grand minima is quite robust.

The question of whether the occurrence of grand minima in solar activity is a regular or chaotic process is important for understanding the action of the solar-dynamo machine. Even a simple deterministic numerical dynamo model can produce events comparable with grand minima (Brandenburg et al., 1989). Such models can also simulate a sequence of grand minima occurrences, which are irregular and seemingly chaotic (e.g., Jennings and Weiss, 1991; Tobias et al., 1995; Covas et al., 1998). The presence of long-term dynamics in the dynamo process is often explained in terms of the a-effect, which, being a result of the electromotive force averaged over turbulent vortices, can contain a fluctuating part (e.g., Hoyng, 1993; Ossendrijver et al., 1996) leading to

irregularly occurring grand minima (e.g., Brandenburg and Spiegel, 2008). The present dynamo models can reproduce almost all the observed features of the solar cycle under ad hoc assumptions (e.g., Pipin et al., 2012), although it is still unclear what leads to the observed variability. Most of these models predict that the occurrence of grand minima is a purely random "memoryless" Poisson-like process, with the probability of a grand minimum occurring being constant at any given time. This unambiguously leads to the exponential shape of the waiting-time distribution (waiting time is the time interval between subsequent events) for grand minima.

Usoskin et al. (2007) performed a statistical analysis of grand minima occurrence time (Table 1) and concluded that their occurrence is not a result of long-term cyclic variations, but is defined by stochastic/chaotic processes. Moreover, waiting-time distribution deviates from the exponential law. This implies that the event occurrence is still random, but the probability is nonuniform in time and depends on the previous history. In the time series it is observed as a tendency of the events to cluster together with a relatively-short waiting time, while the clusters are separated by long event-free intervals (cf. Section 4.1). Such behavior can be interpreted in different ways, e.g., self-organized criticality or processes related to accumulation and release of energy. This poses a strong observational constraint on theoretical models aiming to explain the long-term evolution of solar activity (Section 4.4.1). However, as discussed by Moss et al. (2008) and Usoskin et al. (2009d), the observed feature can be an artefact of the small statistics (only 27 grand minima are identified during the Holocene), making this result only indicative and waiting for a more detailed investigation.

A histogram of the duration of grand minima from Table 1 is shown in Figure 23. The mean duration is 70 year but the distribution is bimodal. The minima tend to be either of a short (30-90 years) duration similar to the Maunder minimum, or rather long (> 100 years), similar to the Sporer minimum, in agreement with earlier conclusions (Stuiver and Braziunas, 1989). This suggests that grand minima correspond to a special state of the dynamo. Once falling into a grand minimum as a result of a stochastic/chaotic, but non-Poisson process, the dynamo is "trapped" in this state and its behavior is driven by deterministic intrinsic features.


Duration (Years)

Figure 23: Histogram of the duration of grand minima from Table 1.

4.3 Grand maxima of solar activity

4.3.1 The modern episode of active sun

In the last decades we were living in a period of a very active sun with a level of activity that is unprecedentedly high for the last few centuries covered by direct solar observation. The sunspot number was growing rapidly between 1900 and 1940, with more than a doubling average group sunspot number, and has remained at that high level until recently (see Figure 1). Note that growth comes mostly from raising the cycle maximum amplitude, while sunspot activity always returns to a very low level around solar cycle minima. While the average group sunspot number for the period 1750-1900 was 35 ± 9 (39 ± 6, if the Dalton minimum in 1797-1828 is not counted), it stands high at the level of 75 ± 3 for 1950-2000. Therefore, the modern active sun episode, which started in the 1940s, can be regarded as the modern grand maximum of solar activity, as opposed to a grand minimum (Wilson, 1988b). As first shown by Usoskin et al. (2003c) and Solanki et al. (2004), such high activity episodes occur quite seldom.

However, as we can securely say now, after the very weak solar minimum in 2008-2009 (e.g., Gibson et al., 2011), solar activity returns to its normal moderate level, or perhaps even to a low-activity stage, comparable to the Dalton minimum in the turn of 18-19th centuries (e.g., Lockwood et al., 2011). Thus, the high activity episode known as the Modern grand maximum is over.

Is such high solar activity typical or is it something extraordinary? While it is broadly agreed that the modern active sun episode is a special phenomenon, the question of how (a)typical such upward bumps are from "normal" activity is a topic of hot debate.

4.3.2 Grand maxima on a multi-millennial timescale

The question of how often grand maxima occur and how strong they are, cannot be studied using the 400-year-long series of direct observations. An increase in solar activity around 1200 AD, also related to the Medieval temperature optimum, is sometimes qualitatively regarded as a grand maximum (Wilson, 1988b; de Meyer, 1998), but its magnitude is lower than the modern maximum (Usoskin et al., 2003c). Accordingly, it was not included in a list of grand maxima by Eddy (1977a,b).

A quantitative analysis is only possible using proxy data, especially cosmogenic isotope records. Using a physics-based analysis of solar-activity series reconstructed from 10Be data from polar (Greenland and Antarctica) archives, Usoskin et al. (2003c, 2004) stated that the modern maximum is unique in the last millennium. Then, using a similar analysis of the 14C calibrated series, Solanki et al. (2004) found that the modern activity burst is not unique, but a very rare event, with the previous burst occurring about 8 millennia ago. An update (Usoskin et al., 2006a) of this result, using a more precise paleo-magnetic reconstruction by Korte and Constable (2005) since 5000 BC, suggests that an increase of solar activity comparable with the modern episode might have taken place around 2000 BC, i.e., around 4 millennia ago. This result is confirmed by the most recent composite reconstruction by Steinhilber et al. (2012). The result by Solanki et al. (2004) has been disputed by Muscheler et al. (2005) who claimed that equally high (or even higher) solar-activity bursts occurred several times during the last millennium, circa 1200 AD, 1600 AD and at the end of the 19th century. We note that the latter claimed peak (ca. 1860) is not confirmed by direct solar or geomagnetic data. However, as argued by Solanki et al. (2005), the level of solar activity reconstructed by Muscheler et al. (2005) was overestimated because of an erroneous normalization to the data of ground-based ionization chambers (see also McCracken and Beer, 2007). This indicates that the definition of grand maxima is less robust than grand minima and is sensitive to other parameters such as geomagnetic field data or overall normalization.

Keeping possible uncertainties in mind, let us consider a list of the largest grand maxima

Table 2: Approximate dates (in -BC/AD) of grand maxima in the SN-L series (after Usoskin et al., 2007).

No. center duration

11 1960 80

2 -445 40

3 -1790 20

4 -2070 40

5 -2240 20

6 -2520 20

7 -3145 30

8 -6125 20

9 -6530 20

10 -6740 100

11 -6865 50

12 -7215 30

13 -7660 80

14 -7780 20

15 -7850 20

16 -8030 50

17 -8350 70

18 -8915 190

19 -9375 130

f Center and duration of the modern maximum are preliminary since it is still ongoing.

(the 50 year smoothed sunspot number stably exceeding 50), identified for the last 11,400 years using 14C data, as shown in Table 2 (after Usoskin et al., 2007). A total of 19 grand maxima have been identified with a total duration of around 1030 years, suggesting that the sun spends around 10% of its time in an active state. A statistical analysis of grand-maxima-occurrence time suggests that they do not follow long-term cyclic variations, but like grand minima, are defined by stochastic/chaotic processes. The distribution of the waiting time between consecutive grand maxima is not unambiguously clear, but also hints at a deviation from exponential law. The duration of grand maxima has a smooth distribution, which nearly exponentially decreases towards longer intervals. Most of the reconstructed grand maxima (about 75%) were not longer than 50 years, and only four grand minima (including the modern one) have been longer than 70 years (cf. Barnard et al., 2011). This suggest that the probability of the modern active-sun episode continuing is low4 (cf. Solanki et al., 2004; Abreu et al., 2008).

4.4 Related implications

Reconstructions of long-term solar activity have different implications in related areas of science. The results, discussed in this overview, can be used in such diverse research disciplines as theoretical astrophysics, solar-terrestrial studies, paleo-climatology, and even archeology and geology. We will not discuss all possible implications of long-term solar activity in great detail but only briefly mention them here.

4 This is not a prediction of future solar activity, but only a statistical estimate.

4.4.1 Theoretical constrains

The basic principles of the occurrence of the 11-year Schwabe cycle are more-or-less understood in terms of the solar dynamo, which acts, in its classical form (e.g., Parker, 1955), as follows (see detail in Charbonneau, 2010). Differential rotation Q produces a toroidal magnetic field from a poloidal one, while the "a-effect", associated with the helicity of the velocity field or Joy's Law tilt of active regions, produces a poloidal magnetic field from a toroidal one. This classical model results in a periodic process in the form of propagation of a toroidal field pattern in the latitudinal direction (the "butterfly diagram"). As evident from observation, the solar cycle is far from being a strictly periodic phenomenon, with essential variations in the cycle length and especially in the amplitude, varying dramatically between nearly spotless grand minima and very large values during grand maxima. The mere fact of such great variability, known from sunspot data, forced solar physicists to develop dynamo models further. Simple deterministic numerical dynamo models, developed on the basis of Parker's migratory dynamo, can simulate events, which are seemingly comparable with grand minima/maxima occurrence (e.g., Brandenburg et al., 1989). However, since variations in the solar-activity level, as deduced from cosmogenic isotopes, appear essentially nonperiodic and irregular, appropriate models have been developed to reproduce irregularly-occurring grand minima (e.g., Jennings and Weiss, 1991; Tobias et al., 1995; Covas et al., 1998). Models, including an ad hoc stochastic driver (Choudhuri, 1992; Schmitt et al., 1996; Ossendrijver, 2000; Weiss and Tobias, 2000; Mininni et al., 2001; Charbonneau, 2001; Charbonneau et al., 2004), are able to reproduce the great variability and intermittency found in the solar cycle (see the review by Charbonneau, 2010). A recent statistical result of grand minima occurrence (Usoskin et al., 2007, Section 4.3.2) shows disagreement between observational data, depicting a degree of self-organization or "memory", and the above dynamo model, which predicts a pure Poisson occurrence rate for grand minima (see Section 4.2). This poses a new constraint on the dynamo theory, responsible for long-term solar-activity variations (Sokoloff, 2004; Moss et al., 2008).

In general, the following additional constraints can be posed on dynamo models aiming to describe the long-term (during the past 11,000 years) evolution of solar magnetic activity.

• The sun spends about 3/4 of its time at moderate magnetic-activity levels, about 1/6 of its time in a grand minimum and about 1/5 -1/10 in a grand maximum. Recent solar activity corresponds to a grand maximum, which has ceased after solar cycle 23.

• Occurrence of grand minima and maxima is not a result of long-term cyclic variations but is defined by stochastic/chaotic processes.

• Observed statistics of the occurrence of grand minima and maxima display deviation from a "memory-less" Poisson-like process, but tend to either cluster events together or produce long event-free periods. This can be interpreted in different ways, such as self-organized criticality (e.g., de Carvalho and Prado, 2000), a time-dependent Poisson process (e.g., Wheatland, 2003), or some memory in the driving process (e.g., Mega et al., 2003).

• Grand minima tend to be of two different types: short minima of Maunder type and long minima of Sporer type. This suggests that a grand minimum is a special state of the dynamo.

• Duration of grand maxima resemble a random Possion-like process, in contrast to grand minima.

4.4.2 Solar-terrestrial relations

The sun ultimately defines the climate on Earth supplying it with energy via radiation received by the terrestrial system, but the role of solar variability in climate variations is far from being clear. Solar variability can affect the Earth's environment and climate in different ways (see, e.g., reviews

by Haigh, 2007; Gray et al., 2010). Variability of total solar irradiance (TSI) measured during recent decades is known to be too small to explain observed climate variations (e.g., Foukal et al., 2006; Frohlich, 2006). On the other hand, there are other ways solar variability may affect the climate, e.g., an unknown long-term trend in TSI (Solanki and Krivova, 2004; Wang et al., 2005) or a terrestrial amplifier of spectral irradiance variations (Shindell et al., 1999; Haigh et al., 2010). Uncertainties in the TSI/SSI reconstructions remain large (Shapiro et al., 2011; Schmidt et al., 2012), making it difficult to assess climate models on the long-term scale. Alternatively, an indirect mechanism also driven by solar activity, such as ionization of the atmosphere by CR (Usoskin and Kovaltsov, 2006) or the global terrestrial current system (Tinsley and Zhou, 2006) can modify atmospheric properties, in particular cloud cover (Ney, 1959; Svensmark, 1998; Usoskin and Kovaltsov, 2008b). Even a small change in cloud cover modifies the transparency/absorption/reflectance of the atmosphere and affects the amount of absorbed solar radiation, even without changes in the solar irradiance. However, the direct role of this effect is estimated to be small (Usoskin et al., 2008; Gray et al., 2010).

Accordingly, improved knowledge of the solar driver's variability may help in disentangling various effects in the very complicated system that is the terrestrial climate (e.g., de Jager, 2005; Versteegh, 2005; Gray et al., 2010). It is of particular importance to know the driving forces in the pre-industrial era, when all climate changes were natural. Knowledge of the natural variability can lead to an improved understanding of anthropogenic effects upon the Earth's climate.

Studies of the long-term solar-terrestrial relations are mostly phenomenological, lacking a clear quantitative physical mechanism. Even phenomenological and empirical studies suffer from large uncertainties, related to the quantitative interpretation of proxy data, temporal and spatial resolution (Versteegh, 2005). Therefore, more precise knowledge of past solar activity, especially since it is accompanied by continuous efforts of the paleo-climatic community on improving climatic data sets, is crucial for improved understanding of the natural (including solar) variability of the terrestrial environment.

4.4.3 Other issues

The proxy method of solar-activity reconstruction, based on cosmogenic isotopes, was developed from the radiocarbon dating method, when it was recognized that the production rate of 14 C is not constant and may vary in time due to solar variability and geomagnetic field changes. Neglect of these effects can lead to inaccurate radiocarbon (or more generally, cosmogenic nuclide) dating, which is a key for, e.g., archeology and Quaternary geology. Thus, knowledge of past solar activity and geomagnetic changes allows for the improvement of the quality of calibration curves, such as the IntCal (Stuiver et al., 1998; Reimer et al., 2004, 2009) for radiocarbon, eventually leading to more precise dating.

Long-term variations in the geomagnetic field are often evaluated using cosmogenic isotope data. Knowledge of source variability due to solar modulation is important for better results.

4.5 Summary

In this section, solar activity on a longer scale is discussed, based on recent reconstructions.

According to these reconstructions, the sun has spent about 70% of its time during the Holocene, which is ongoing, in a normal state characterized by medium solar activity. About 15-20% of the time the sun has experienced a grand minimum, while 10-15% of the time has been taken up by periods of very high activity.

One of the main features of long-term solar activity is its irregular behavior, which cannot be described by a combination of quasi-periodic processes as it includes an essentially random component.

Grand minima, whose typical representative is the Maunder minimum of the late 17th century, are typical solar phenomena. A total of 27 grand minima have been identified in reconstructions of the Holocene period. Their occurrence suggests that they appear not periodically, but rather as the result of a chaotic process within clusters separated by 2000-2500 years. Grand minima tend to be of two distinct types: short (Maunder-like) and longer (Sporer-like). The appearance of grand minima can be reproduced by modern stochastic-driven dynamo models to some extent, but some problems still remain to be resolved.

The modern level of solar activity (after the 1940s) was very high, corresponding to a grand maximum, which are typical but rare and irregularly-spaced events in solar behavior. However, this grand maximum has ceased after solar cycle 23. The duration of grand maxima resembles a random Possion-like process, in contrast to grand minima.

These observational features of the long-term behavior of solar activity have important implications, especially for the development of theoretical solar-dynamo models and for solar-terrestrial studies.

5 Solar Energetic Particles in the Past

In addition to galactic cosmic rays, which are always present in the Earth's vicinity, sometimes sporadic solar energetic-particle (SEP) events with a greatly enhanced flux of less energetic particles in the interplanetary medium also occur (e.g., Klecker et al., 2006). Strong SEP events mostly originate from CME-related shocks propagating in the solar corona and interplanetary medium, that lead to effective bulk acceleration of charged particles (e.g., Cane and Lario, 2006). Although these particles are significantly less energetic than GCRs, they can occasionally be accelerated to an energy reaching up to several GeV, which is enough to initiate the atmospheric cascade. Peak intensity of SEP flux can be very high, up to 104 particles (with energy > 30 MeV) per cm2 per second. In fact, the long-term average flux (or fluence) of SEP is mostly defined by rare major events, which occur a few times per solar cycle, with only minor contributions from a large number of weak events (Shea and Smart, 1990, 2002). As an example, energy spectra of GCR and SEP are shown in Figure 24 for the day of January 20, 2005, when an extreme SEP event took place. Such SEPs dominate the low-energy section of cosmic rays (below hundreds of MeV of a particle's kinetic energy), which is crucial for the radiation environment, and play an important role in solar-terrestrial relations. For many reasons it is important to know the variations of SEPs on long-term scales.

10 1091 ioS

w 10o!

E 10% o

J 103,

0.01 0.1 1

Energy (GeV)

Figure 24: Daily fluence of solar energetic particles (dashed curve - Tylka and Dietrich, 2009) and galactic cosmic rays (solid curve) for the day of January 20, 2005. Open circles represent space-borne measurements (Mewaldt, 2006; Mewaldt et al., 2012).

It is not straightforward to evaluate the average SEP flux even for the modern instrumental epoch of direct space-borne measurements (e.g., Mewaldt et al., 2007). For example, estimates for the average flux of SEPs with an energy above 30 MeV (called f 30 henceforth) for individual cycles may vary by an order of magnitude, from 10 cm-2 s-1 for cycle 21 up to 70 cm-2 s-1 for cycle 19 (Reedy, 2012). Moreover, estimates of the SEP flux were quite uncertain during the earlier years of space-borne measurements because of two effects, which are hard to account for (e.g., Reeves et al., 1992; Tylka et al., 1997). One is related to the very high flux intensities of SEPs during the peak phase of events, when a detector can be saturated because of the dead-time effect (the maximum trigger rate of the detector is exceeded). The other is related to events with high energy solar particles, which can penetrate into the detector through the walls of the collimator or the

detector, leading to an enhanced effective acceptance cone with respect to the "expected" one. Since the SEP fluence is defined by major events, these effects may lead to an underestimate of the average flux of SEPs. The modern generation of detectors are better suited for measuring high fluxes. The average f 30 flux for the last five solar cycles (1954-2006) is estimated at about 35 cm-2 s-1 (Smart and Shea, 2002; Shea et al., 2006).

5.1 Cosmogenic isotopes

The development of the method of cosmogenic isotopes makes it possible to estimate occurrence of extreme SEP events in the past. Some earlier attempts were inconclusive. For example, Usoskin et al. (2006b) found that a typical strong SEP event leaves no distinguishable signature in 14C but may be observed from ice core 10Be records. However, the question of the possible rare occurrence of extreme SEP events on the millennial time scale is important not only from the theoretical point of view, but also for assessment of radiation risks for space-borne missions, especially manned ones. What can be the strongest SEP event originated from the sun, how often they can occur? These questions need to be answered. Several attempts have been made to evaluate that from the cosmogenic isotope data (Lingenfelter and Hudson, 1980; Usoskin et al., 2006b; Webber et al., 2007), but the result was grossly uncertain (Hudson, 2010; Schrijver et al., 2012), mostly because of the large model uncertainties of the radionuclide production.

A new step forward has been done recently by Usoskin and Kovaltsov (2012), who analyzed two 14C and five 10Be records over the last millennia and searched for possible signatures of extreme SEP events.

Years AD

Figure 25: Time profiles of the measured A14C content in Japanese cedar (M12 - Miyake et al., 2012) and German oak (ETH ZUrich & Mannheim AMS - Usoskin et al., 2013) trees for the period around 775 AD. Smooth black and grey lines depict a family of best fit A14C profiles, calculated using a family of realistic carbon cycle models for an instantaneous injection of 14C into the stratosphere (Usoskin and Kovaltsov, 2012). Image after Usoskin et al. (2013).

While the response of 10Be to an SEP event is simply a 1 -2-yr long peak, because of the simple atmospheric transport/deposition (see Section 3.3.3), the response of 14C has a typical shape shown in Figure 25 - with a sharp peak and exponential decay of the length of several decades, due to the carbon cycle (see Section 3.2.3). Usoskin and Kovaltsov (2012) checked all the available cosmogenic isotope data through the entire Holocene looking for a potential SEP signatures, and came up with a list of candidates of extreme SEP events and assessments of their strength (Table 3).

Table 3: A list of candidates for extreme SEP events found in different cosmogenic isotope records throughout the Holocene: approximate year, dataset used (Dye3 - McCracken et al. (2004); NGRIP -Berggren et al. (2009); IntCal09 - Reimer et al. (2009); GRIP - Yiou et al. (1997); Dome Fuji - Horiuchi et al. (2008); South Pole - Raisbeck et al. (1990); M12 - Miyake et al. (2012)), and the F30 fluence [cm-2]. Table after Usoskin and Kovaltsov (2012).

SPE year Series F 3 0

1460-1462 AD NGRIP(1460) 1.5 x 1010

Dye3 (1462) 9.7 x 109

1505 AD Dye3 1.3 x 1010

1719 AD NGRIP 1x 1010

1810 AD NGRIP 1x 1010

8910 BC IntCal09 2.0 x 1010

8155 BC IntCal09 1.3 x 1010

8085 BC IntCal09 1.5 x 1010

7930 BC IntCal09 1.3 x 1010

7570 BC IntCal09 2.0 x 1010

7455 BC IntCal09 1.5 x 1010

6940 BC IntCal09 1.1 x 1010

6585 BC IntCal09 1.7 x 1010

5835 BC IntCal09 1.5 x 1010

5165 BC GRIP 2.4 x 1010

4680 BC IntCal09 1.6 x 1010

3260 BC IntCal09 2.4 x 1010

2615 BC IntCal09 1.2 x 1010

2225 BC IntCal09 1.2 x 1010

1485 BC IntCal09 2.0 x 1010

95 AD GRIP 2.6 x 1010

265 AD IntCal09 2.0 x 1010

785 AD IntCal09 2.4 x 1010

Dome Fuji 5.3 x 1010 f

M12 4 x 1010 f

1455 AD South Pole 7.0 x 1010 f

f Upper bound.

The list includes 23 candidates for extreme SEP events with the fluence F30 exceeding 1010 cm-2, viz. the greatest fluence observed for the space era in 1960 (Shea and Smart, 1990). Note that only two of these candidates appear in more than one series - the events of ca. 1460 AD and ca. 780 AD. The former had signatures in two annual 10Be series, NGRIP and Dye3. The later was observed in two 14C series, biennial M12 and 5-yr IntCal09, and in quasi-decadal Dome Fuji 10Be series. The quasi-decadal South Pole 10Be series does not show an increase ca. 780 AD placing an upper limit on the strength of the event.

We note that the event of ca. 775 AD was analyzed using biennial 14C data by Miyake et al. (2012), who suggested that the event was probably caused by 7-rays from an unknown nearby supernova. This event is confirmed by annual 14C data from a German oak tree (Usoskin et al., 2013). However, because of the use of an inappropriate carbon cycle model, Miyake et al. (2012) grossly (by a factor of 5) overestimated the corresponding 14C production, leading to the need of a supernova. Moreover, this leads to a strong disagreement between 14C and 10Be data sets, since this event is not observed in the South Pole record and is not exceptionally strong in the Dome

Fuji record. However, if an appropriate model of the carbon cycle is used, the production of 14C appears in a reasonable range, being consistent with 10Be (Usoskin et al, 2013). Therefore, there is no need to involve such an exotic object as a nearby supernova whose remnants are unknown for us - the event of ca. 775-780 AD can be consistently explained by a extreme but not exceptional SEP event.

The integral probability distribution of the occurrence of strong SEP events, as revealed from the cosmogenic isotope data, is shown in Figure 26.

cu o.oi -

1E-4 -

1E-5 -

1 I y X t

V / Y ■A I * O / •

Space era Isotope (annual) Isotope (decadal)

Annual fluence F30 [10 cm" ]

Figure 26: Cumulative probability (with the 90% confidence interval) of occurrence of a SEP event with fluence (> 30 MeV) exceeding the given value F30, as assessed from the data, for the space era 1956-2008 (triangles), cosmogenic isotope annual data, (stars), and cosmogenic isotope decadal data, (circles). Gray dotted curve depicts the best-fit exponent. Image reproduced by permission from Usoskin and Kovaltsov (2012), copyright by AAS.

One can see that the break in the distribution marginally hinted in the directly observed SEP events at around F30 = (5-7) xlO9 cm-2 (nonproportiona.lly fewer strong events observed) is confidently confirmed by the cosmogenic isotope data.. In particular, no event with F30 > 2 x 1010 citT2 was found over the last 600 years using annually resolved 10Be data.. It is noteworthy that the idea, of an possible extreme Ca.rrington SPE of 1859 AD (McCra.cken et al, 2001) is discarded (see also Wolff et al, 2012). On the longer time scale of 11 millennia., no event with F-^o > 5 x 1010 citT2 has been found. This gives a. new strict observational constraint on the occurrence probability of extreme SPEs.

According to Usoskin and Kovaltsov (2012) practical limits can be set as F30 ~ 1, 2-3 and 5 x 1010 citT2 (10, 20-30 and 50 times greater than the SEP event of February 23, 1956), for the occurrence probability of 10~2, 10~3, and 10™4 yr_1, respectively. The mean SEP flux is found as « 40 (cm2 s)_1 in agreement with estimates from the lunar rocks. On average, extreme SPEs contribute about 10% to the total SEP fluence.

5.2 Lunar and meteoritic rocks

Since energy spectra of SEP and GCR are dramatically different, one may think of a natural spectrometer to separate their effects and thus evaluate their fluxes independently. A spectrometer that is able to separate cosmic rays is lunar (or meteoritic) rocks.

depth (g/cm2)

Figure 27: Measured (dots) and calculated (curves) 14C activity in a lunar sample 68815 (Jull et al., 1998). The big diamond implies contamination of a thin surface layer by 14C implanted from solar wind. The dotted curve represents the expected production due to GCR, while the solid curve is the best fit SEP+GCR model production.

Figure 27 depicts an example of 14C measured in a lunar sample (Jull et al., 1998). The dotted line shows the expected production of radiocarbon by GCR. The production increases with depth due to the development of a nucleonic cascade in the matter, initiated by energetic GCR particles, similar to the atmospheric cascade. Less energetic particles of solar origin produce the isotope only in upper layers of the rock, since their low energy does not allow them to initiate a cascade. On the other hand, thanks to their high flux in the lower energy range, the production of 14C in the upper layers is much higher than that from GCR. Thus, by first measuring the isotope activity in deep layers one can evaluate the average GCR flux, and then the measured excess in the upper level yields an estimate for the SEP flux in both integral intensity and spectral shape. The result is based on model computations and therefore is slightly model dependent but makes it possible to give a robust estimate of the GCR and SEP in the past.

A disadvantage of this approach is that lunar samples are not stratified and do not allow for temporal separation. The measured isotope activity is a balance between production and decay and, therefore, represents the production (and the ensuing flux) integrated over the life-time of the isotope before the sample has been measured. However, using different isotopes with different life times, one can evaluate the cosmic-ray flux integrated over different timescales.

Estimates of the average SEP flux f 30 on different timescales, as obtained from various isotopes measured in lunar samples, are collected in Table 4. Based on isotopes with different life-times (see Table 4) one can evaluate the average flux of SEP on different time scale (see Figure 28). The average f 30 flux for the last five solar cycles (1954-2008) is consistent with the average flux estimated in the past for longer timescales from 103 to 107 years (cf. Reedy, 2002, 2012).

However, this method is not able to provide an estimate of the occurrence rate of extreme SEP events. If one assumes that the entire average SEP flux is produced within one extreme event occurring at half of the isotope's life-time ago (Reedy, 1996), an upper limit for the occurrence of extreme SEP events can be placed. This is an unrealistically extreme assumption, which may

Table 4: Estimates of 4k omni-directional integral (above 30 MeV) flux, f 30 in [cm2 s] 1, of solar energetic particles, obtained from different sources.

Timescale Method Source Reference f 30 (cm 2 s x)

1954-2008 measurements space-borne Reedy (2012) 35

104 yr 14C lunar rock Jull et al. (1998) 42

105 yr 41 Ca lunar rock Fink et al. (1998) 56

5 x 105 yr 36Cl lunar rock Nishiizumi et al. (2009) 46

106 yr 26 Al lunar rock Kohl et al. (1978) 25

106 yr 26 Al lunar rock Grismore et al. (2001) 55

106 yr 10Be, 26Al lunar rock Michel et al. (1996) 24

106 yr 10Be, 26Al lunar rock Fink et al. (1998) 32

106 yr 10Be, 26Al lunar rock Nishiizumi et al. (2009) 24

2 x 106 yr 10Be, 26Al lunar rock Nishizumi et al. (1997) - 35

5 x 106 yr 53Mn lunar rock Kohl et al. (1978) 25

2 x 106 yr 21Ne, 22Ne, 38Ar lunar rock Rao et al. (1994) 22

lead to an overestimate by many orders of magnitude, but it sets the very conservative upper limit which cannot be exceeded.

5.3 Nitrates in polar ice

It has been discussed until recently that another quantitative index of strong SEP events (with F30 x 109 cm-2) might be related to nitrate (NO-) records measured in polar ice cores. The concentration of nitrates has been measured in polar ice from both the Southern (South Pole, e.g., Dreschhoff and Zeller, 1990) and Northern (Greenland, e.g., Zeller and Dreschhoff, 1995; Dreschhoff and Zeller, 1998) polar caps, depicting pronounced spikes associated with strong SEP events (McCracken et al., 2001). As a result of the analysis a list of large SEP events since 1560 and their fluences have been published (see Table 1 in McCracken et al., 2001) and widely used.

However, as shown by several independent recent studies (Wolff et al., 2012; Usoskin and Kovaltsov, 2012) on the example of the Carrington event (September 1859), the nitrate spikes are not related to SEP events. According to McCracken et al. (2001), the nitrate spike and the associated SEP event was the strongest in the entire record (F30 « 2 x 1010 cm-2). Wolff et al. (2012) have measured, with high resolution, nitrate content in 14 ice cores from Antarctic and Greenland for a few decades around 1859. Only one Greenland series depicts a spike which can be associated with the event, all other series have no signatures. Moreover, all similar spikes found in Greenland datasets are accompanied by chemical tracers (ammonium, formate, black carbon, etc.) clearly pointing to the anthropogenic source of nitrates - biomass burning plumes. No significant spikes have been found in the Antarctic records. Wolff et al. (2012) concluded that "Nitrate spikes cannot be used to derive the statistics of SEPs." Another confirmation of this conclusion was made by Usoskin and Kovaltsov (2012), who calculated, from the F30 fluence proposed by McCracken et al. (2001) for the Carrington event, the 10Be production. If the Carrington SEP event was so strong, it would have necessarily left its clear signature in the annually resolved 10Be record, which however contradicts to the real data from NGRIP and Dye3 ice cores.

Thus, the nitrate record in polar ice cannot serve as an index of SEP events. On the other hand, it may be used to study long-term variability of GCR (see Section 3.4).

5.4 Summary

In this section, estimates of the averaged long-term flux of SEPs are discussed.

Measurements of cosmogenic isotopes with different life times in lunar and meteoritic rocks allow one to make rough estimates of the SEP flux over different timescales. The directly space-borne-measured SEP flux for past decades is broadly consistent with estimates on longer timescales - up to millions of years. The same measurements can provide a very conservative upper estimate for the occurrence rate of extreme SEP events. Terrestrial cosmogenic isotope data in dated archives (tree trunks, ice cores) give a possibility to assess the occurrence rate of strong SEP events on the time scales up to ten of millennia. Measurements of nitrates in polar ice have been shown to be an invalid index of strong SEP events in the past.

Different estimates of the extreme (quantified as the fluence of SEP with energy above 10 MeV) SEP event occurrence probability are summarized in Figure 28).

An analysis of various kinds of data suggests that the distribution of the intensity of SEP events has a break, and the occurrence of extra-strong events (with the F30 fluence exceeding 5 x 1010 cm-2) is unlikely on the multi-millennial time scale.

104 106 108 1010 1012 1014 1016 Fluence (cm-2, above 10 MeV)

Figure 28: Cumulative frequency distribution of SEP events with fluences greater than F10 (for particles with energies above 10 MeV). Red histogram: satellite-based direct observations; Blue diamonds: conservative upper limits derived from lunar isotopes (see Section 5.2); Blue dashed line: upper limit based on 14C record (Hudson, 2010); Image reproduced by permission from Schrijver et al. (2012), copyright by AGU.

6 Conclusions

In this review the present knowledge of long-term solar activity on a multi-millennial timescale, as reconstructed using the indirect proxy method, is discussed.

Although the concept of solar activity is intuitively understandable as a deviation from the "quiet" sun concept, there is no clear definition for it, and different indices have been proposed to quantify different aspects of variable solar activity. One of the most common and practical indices is sunspot number, which forms the longest available series of direct scientific observations. While all other indices have a high correlation with sunspot numbers, dominated by the 11-year cycle, the relationship between them at other timescales (short- and long-term trends) may vary to a great extent.

On longer timescales, quantitative information of past solar activity can only be obtained using the method based upon indirect proxy, i.e., quantitative parameters, which can be measured nowadays but represent the signatures, stored in natural archives, of the different effects of solar magnetic activity in the past. Such traceable signatures can be related to nuclear or chemical effects caused by cosmic rays in the Earth's atmosphere, lunar rocks or meteorites. The most common proxy of solar activity is formed by data from the cosmogenic radionuclides, 10Be and 14C, produced by cosmic rays in the Earth's atmosphere and stored in independently-dated stratified natural archives, such as tree rings or ice cores. Using a recently-developed physics-based model it is now possible to reconstruct the temporal behavior of solar activity in the past, over many millennia. The most robust results can be obtained for the Holocene epoch, which started more than 11,000 years ago, whose stable climate minimizes possible uncertainties in the reconstruction. An indirect verification of long-term solar-activity reconstructions supports their veracity and confirms that variations of cosmogenic nuclides on the long-term scale (centuries to millennia) during the Holocene make a solid basis for studies of solar variability in the past. However, such reconstructions may still contain systematic uncertainties related to unknown changes in the geomagnetic field or climate of the past, especially in the early part of the Holocene.

Measurements of the concentration of different cosmogenic isotopes in lunar and meteoritic rocks make it possible to estimate the SEP flux on different timescales. Directly space-borne-measured SEP flux for recent decades is broadly consistent with estimates on longer timescales - up to millions of years. The occurrence of extra-strong events, with the fluence of SEP (with energy greater than 30 MeV) exceeding 5 x 1010 cm-2 is unlikely on the multimillenial time scale.

In general, the following main features are observed in the long-term evolution of solar magnetic activity.

• Solar activity is dominated by the 11-year Schwabe cycle on an interannual timescale. Some additional longer characteristic times can be found, including the Gleissberg secular cycle, de Vries/Suess cycle, and a quasi-cycle of 2000-2400 years. However, all these longer cycles are intermittent and cannot be regarded as strict phase-locked periodicities.

• One of the main features of long-term solar activity is that it contains an essential chaotic/ stochastic component, which leads to irregular variations and makes solar-activity predictions impossible for a scale exceeding one solar cycle.

• The sun spends about 70% of its time at moderate magnetic activity levels, about 15-20% of its time in a grand minimum and about 10-15% in a grand maximum.

• Grand minima are a typical but rare phenomena in solar behavior. Their occurrence appears not periodically, but rather as the result of a chaotic process within clusters separated by 2000-2500 years. Grand minima tend to be of two distinct types: short (Maunder-like) and longer (Sporer-like).

• The recent level of solar activity (after the 1940s) was very high, corresponding to a prolonged grand maximum, but it is ceasing now to the normal moderate level. Grand maxima are also rare and irregularly occurring events, though the exact rate of their occurrence is still a subject of debates.

These observational features of the long-term behavior of solar activity have important implications, especially for the development of theoretical solar-dynamo models and for solar-terrestrial studies.

7 Acknowledgements

I thank the editorial board of the Living Reviews in Solar Physics for the invitation to prepare this review and organizational aid. This work would be impossible without encouraging discussions, as well as direct and indirect help, from my colleagues Gennady Kovaltsov and Sami Solanki, whose invaluable support is acknowledged with my greatest gratitude. I have also greatly benefited from lively and open discussions and debates with Rainer Arlt, Edouard Bard, Jiirg Beer, Narendra Bhandari, late Giuliana Cini Castagnoli, Paul Charbonneau, Ed Cliver, Vincent Cour-tillot, Cornelis de Jager, Gisella Dreschhoff, Erwin Fluckiger, Peter Foukal, Yves Gallet, Agnes Genevey, David Hathaway, Monika Korte, Natalia Krivova, Bernd Kromer, Devendra Lal, Michael Lockwood, Ken McCracken, David Moss, Kalevi Mursula, Raimund Muscheler, Hiroko Miyahara, Heikki Nevanlinna, Keran O'Brien, Valery Ostryakov, Alexander Ruzmaikin, Dieter Schmitt, Manfred Schiissler, Margaret Shea, Don Smart, Ian Snowball, Dmitry Sokoloff, Willie Soon, Leif Sval-gaard, Carla Taricco, Rita Traversi, Jose Vaquero, and many others. I thank Evguenia Usoskina for help with editing this text. I am happy to acknowledge those colleagues who work in this field and make it living and vibrant.


Abreu, J.A., Beer, J., Steinhilber, F., Tobias, S.M. and Weiss, N.O., 2008, "For how long will the current grand maximum of solar activity persist?", Geophys. Res. Lett., 35, L20109. [DOI], [ADS] (Cited on pages 19 and 57.)

Abreu, J.A., Beer, J., Ferriz-Mas, A., McCracken, K.G. and Steinhilber, F., 2012, "Is there a planetary influence on solar activity?", Astron. Astrophys., 548, A88. [DOI], [ADS] (Cited on page 17.)

Argento, D.C., Reedy, R.C. and Stone, J.O., 2013, "Modeling the earth's cosmic radiation", Nucl. Instrum. Methods B, 294, 464-469. [DOI], [ADS] (Cited on page 26.)

Arlt, R., 2008, "Digitization of Sunspot Drawings by Staudacher in 1749-1796", Solar Phys., 247, 399410. [DOI], [ADS] (Cited on pages 11, 12, 13, and 53.)

Arlt, R., 2009, "The Butterfly Diagram in the Eighteenth Century", Solar Phys., 255, 143-153. [DOI], [ADS], [arXiv:0812.2233] (Cited on pages 11, 12, and 53.)

Arlt, R. and Abdolvand, A., 2011, "First solar butterfly diagram from Schwabe's observations in 18251867", in Physics of Sun and Star Spots, Proceedings of IAU Symposium 273 held in Ventura, California, USA, August 22-26, 2010, (Eds.) Choudhary, D.P., Strassmeier, K.G., IAU Symposium, S273, pp. 286289, Cambridge University Press, Cambridge; New York. [DOI], [ADs], [arXiv:1010.3131 [astro-ph.SR]| (Cited on page 12.)

Arlt, R. et al., 2013, Mon. Not. R. Astron. Soc, submitted (Cited on page 12.)

Balmaceda, L.A., Solanki, S.K. and Krivova, N.A., 2005, "A cross-calibrated sunspot areas time series since 1874", Mem. Soc. Astron. Ital, 76, 929-932. [ADS] (Cited on page 13.)

Baranyi, T., Gyori, L., Ludmany, A. and Coffey, H.E., 2001, "Comparison of sunspot area data bases", Mon. Not. R. Astron. Soc., 323, 223-230. [DOI], [ADS] (Cited on page 13.)

Bard, E., Arnold, M., Maurice, P. and Duplessy, J.-C., 1987, "Measurements of bomb radiocarbon in the ocean by means of accelerator mass spectrometry: Technical aspects", Nucl. Instrum. Methods B, 29, 297-301. [DOI], [ADS] (Cited on page 33.)

Bard, E., Raisbeck, G.M., Yiou, F. and Jouzel, J., 1997, "Solar modulation of cosmogenic nuclide production over the last millennium: comparison between 14C and 10Be records", Earth Planet. Sci. Lett., 150, 453-462. [DOI], [ADS] (Cited on pages 21, 31, 32, 33, 35, 36, and 47.)

Barnard, L., Lockwood, M., Hapgood, M.A., Owens, M.J., Davis, C.J. and Steinhilber, F., 2011, "Predicting space climate change", Geophys. Res. Lett., 38, L16103. [DOI], [ADS] (Cited on page 57.)

Basurah, H.M., 2004, "Auroral Evidence for Early High Solar Activities", Solar Phys., 225, 209-212. [DOI], [ADS] (Cited on page 15.)

Beer, J., 2000, "Long-term indirect indices of solar variability", Space Sci. Rev., 94, 53-66. [ADS] (Cited on pages 21, 34, 35, 42, and 43.)

Beer, J., Siegenthaler, U., Oeschger, H., Andrée, M., Bonani, G., Suter, M., Wolfli, W., Finkel, R.C. and Langway, C.C., 1983, "Temporal 10Be Variations", in 18th International Cosmic Ray Conference, Bangalore, India, 22 August-3 September, 1983, (Ed.) Durgaprasad, D. et al., 9, pp. 317-320, TIFR, Bombay. [ADS] (Cited on page 34.)

Beer, J., Siegenthaler, U., Oeschger, H., Bonani, G. and Finkel, R.C., 1988, "Information on past solar activity and geomagnetism from Be-10 in the Camp Century ice core", Nature, 331, 675-679. [DOI], [ADS] (Cited on page 34.)

Beer, J., Blinov, A., Bonani, G., Finkel, R.C., Hofmann, H.J., Lehmann, B., Oeschger, H., Sigg, A., Schwander, J., Staffelbach, T., Stauffer, B., Suter, M. and Wotfli, W., 1990, "Use of 10Be in polar ice to trace the 11-year cycle of solar activity", Nature, 347, 164-166. [DOI], [ADS] (Cited on pages 21 and 34.)

Beer, J., Tobias, S. and Weiss, N., 1998, "An Active Sun Throughout the Maunder Minimum", Solar Phys, 181, 237-249. [ADS] (Cited on page 52.)

Beer, J., Vonmoos, M.V., Muscheler, R., McCracken, K.G. and Mende, W., 2003, "Heliospheric Modulation over the past 10,000 Years as derived from Cosmogenic Nuclides", in 28th International Cosmic Ray Conference, July 31-August 7, 2003, Trukuba, Japan, (Eds.) Kajita, T., Asaoka, Y., Kawachi, A., Matsubara, Y., Sasaki, M., 7, pp. 4147-4150, Universal Academy Press, Tokyo, Japan. [ADS] (Cited on pages 41 and 43.)

Beer, J., McCracken, K.G. and von Steiger, R., 2012, Cosmogenic Radionuclides: Theory and Applications in the Terrestrial and Space Environments, Physics of Earth and Space Environments, Springer, Berlin; Heidelberg (Cited on pages 21, 26, 34, and 36.)

Belmaker, R., Lazar, B., Tepelyakov, N., Stein, M. and Beer, J., 2008, "10Be in Lake Lisan sediments - A proxy for production or climate?", Earth Planet. Sci. Lett., 269, 447-456. [DOI] (Cited on page 34.)

Belov, A.V., Gushchina, R.T., Obridko, V.N., Shelting, B.D. and Yanke, V.G., 2006, "Long-term variations of galactic cosmic rays in the past and future from observations of various solar activity characteristics", J. Atmos. Terr. Phys., 68, 1161-1166. [DOI], [ADS] (Cited on page 40.)

Berggren, A.-M., Beer, J., Possnert, G., Aldahan, A., Kubik, P., Christl, M., Johnsen, S.J., Abreu, J. and Vinther, B.M., 2009, "A 600-year annual 10Be record from the NGRIP ice core, Greenland", Geophys. Res. Lett., 36, L11801. [DOI], [ADS] (Cited on page 63.)

Bezuglov, M.V., Malyshevsky, V.S., Malykhina, T.V., Torgovkin, A.V., Fomin, G.V. and Shramenko, B.I., 2012, "Photonuclear channel of 7Be production in the Earth's atmosphere", Phys. Atom. Nucl., 75, 393-397. [DOI], [ADS] (Cited on page 34.)

Bhandari, P.A., Lal, D. and Rama, 1966, "Stratospheric circulation studies based on natural and artificail radioactive tracer elements", Tellus, 18, 391-406 (Cited on page 35.)

Bonino, G., Castagnoli, G.C., Bhandari, N. and Taricco, C., 1995, "Behavior of the Heliosphere over prolonged solar quiet periods by 44Ti measurements in meteorites", Science, 270, 1648-1650. [DOI], [ADS] (Cited on page 46.)

Bracewell, R.N., 1986, "Simulating the sunspot cycle", Nature, 323, 516-519. [DOI], [ADS] (Cited on page 16.)

Brandenburg, A. and Sokoloff, D., 2002, "Local and Nonlocal Magnetic Diffusion and Alpha-Effect Tensors in Shear Flow Turbulence", Geophys. Astrophys. Fluid Dyn., 96, 319-344. [DOI], [ADS], [arXiv:astro-ph/0111568] (Cited on page 18.)

Brandenburg, A. and Spiegel, E.A., 2008, "Modeling a Maunder minimum", Astron. Nachr., 329, 351-358. [ADS], [arXiv:0801.2156] (Cited on pages 18 and 55.)

Brandenburg, A., Krause, F., Meinel, R., Moss, D. and Tuominen, I., 1989, "The stability of nonlinear dynamos and the limited role of kinematic growth rates", Astron. Astrophys., 213, 411-422. [ADS] (Cited on pages 54 and 58.)

Bray, R.J. and Loughhead, R.E., 1964, Sunspots, International Astrophysics Series, 7, Chapman & Hall, London (Cited on page 14.)

Broeker, W.S. and Peng, T.-H., 1986, "Carbon cycle: 1985 Glacial to interglacial changes in the operation of the global carbon cycle", Radiocarbon, 28, 309-327 (Cited on page 32.)

Brooke, J.M. and Moss, D., 1994, "Non-linear dynamos in torus geometry: transition to chaos", Mon. Not. R. Astron. Soc., 266, 733-739. [ADS] (Cited on page 18.)

Burger, R.A., Potgieter, M.S. and Heber, B., 2000, "Rigidity dependence of cosmic ray proton latitudinal gradients measured by the Ulysses spacecraft: Implications for the diffusion tensor", J. Geophys. Res., 105, 27,447-27,456. [DOI], [ADS] (Cited on pages 23 and 30.)

Bushby, P.J. and Tobias, S.M., 2007, "On Predicting the Solar Cycle Using Mean-Field Models", Astrophys. J., 661, 1289-1296. [DOI], [ADS], [arXiv:0704.2345 [astro-ph]] (Cited on page 19.)

Caballero-Lopez, R.A. and Moraal, H., 2004, "Limitations of the force field equation to describe cosmic ray modulation", J. Geophys. Res., 109, A01101. [DOI], [ADS] (Cited on pages 23, 42, and 44.)

Cane, H.V. and Lario, D., 2006, "An Introduction to CMEs and Energetic Particles", Space Sci. Rev., 123, 45-56. [DOI], [ADS] (Cited on page 61.)

Carbonell, M., Oliver, R. and Ballester, J.L., 1993, "On the asymmetry of solar activity", Astron. Astrophys., 274, 497-504. [ADS] (Cited on page 18.)

Carbonell, M., Oliver, R. and Ballester, J.L., 1994, "A search for chaotic behaviour in solar activity", Astron. Astrophys., 290, 983-994. [ADS] (Cited on page 18.)

Casas, R., Vaquero, J.M. and Vazquez, M., 2006, "Solar Rotation in the 17th century", Solar Phys., 234, 379-392. [DOI], [ADS] (Cited on page 11.)

Castagnoli, G. and Lal, D., 1980, "Solar modulation effects in terrestrial production of carbon-14", Radiocarbon, 22, 133-158 (Cited on pages 23 and 25.)

Charbonneau, P., 2001, "Multiperiodicity, Chaos, and Intermittency in a Reduced Model of the Solar Cycle", Solar Phys., 199, 385-404. [ADS] (Cited on pages 17, 52, 53, and 58.)

Charbonneau, P., 2010, "Dynamo Models of the Solar Cycle", Living Rev. Solar Phys., 7, lrsp-2010-3. [ADS]. URL (accessed 25 August 2012): (Cited on pages 18, 52, and 58.)

Charbonneau, P. and Dikpati, M., 2000, "Stochastic Fluctuations in a Babcock-Leighton Model of the Solar Cycle", Astrophys. J., 543, 1027-1043. [DOI], [ADS] (Cited on page 18.)

Charbonneau, P., Blais-Laurier, G. and St-Jean, C., 2004, "Intermittency and Phase Persistence in a Babcock-Leighton Model of the Solar Cycle", Astrophys. J. Lett., 616, L183-L186. [DOI], [ADS] (Cited on page 58.)

Choudhuri, A.R., 1992, "Stochastic fluctuations of the solar dynamo", Astron. Astrophys., 253, 277-285. [ADS] (Cited on pages 52 and 58.)

Choudhuri, A.R., Chatterjee, P. and Jiang, J., 2007, "Predicting Solar Cycle 24 With a Solar Dynamo Model", Phys. Rev. Lett., 98, 131103. [DOI], [ADS], [arXiv:astro-ph/0701527] (Cited on page 19.)

Clark, D.H. and Stephenson, F.R., 1978, "An Interpretation of the Pre-Telescopic Sunspot Records from the Orient", Quart. J. R. Astron. Soc., 19, 387-410. [ADS] (Cited on pages 14 and 15.)

Clem, J.M. and Dorman, L.I., 2000, "Neutron Monitor Response Functions", Space Sci. Rev., 93, 335-359. [DOI], [ADS] (Cited on page 29.)

Clette, F., Berghmans, D., Vanlommel, P., van der Linden, R.A.M., Koeckelenbergh, A. and Wauters, L., 2007, "From the Wolf number to the International Sunspot Index: 25 years of SIDC", Adv. Space Res., 40, 919-928. [DOI], [ADS] (Cited on page 11.)

Cliver, E.W., Boriakoff, V. and Bounar, K.H., 1998, "Geomagnetic activity and the solar wind during the Maunder Minimum", Geophys. Res. Lett., 25, 897-900. [DOI], [ADS] (Cited on page 52.)

Conway, A.J., 1998, "Time series, neural networks and the future of the Sun", New Astron. Rev., 42, 343-394. [DOI], [ADS] (Cited on page 19.)

Cooke, D.J., Humble, J.E., Shea, M.A., Smart, D.F., Lund, N., Rasmussen, I.L., Byrnak, B., Goret, P. and Petrou, N., 1991, "On cosmic-ray cut-off terminology", Nuovo Cimento C, 14, 213-234. [DOI], [ADS] (Cited on page 23.)

Covas, E., Tavakol, R., Tworkowski, A. and Brandenburg, A., 1998, "Axisymmetric mean field dynamos with dynamic and algebraic a-quenching", Astron. Astrophys., 329, 350-360. [ADS] (Cited on pages 54 and 58.)

Damon, P.E. and Sonett, C.P., 1991, "Solar and terrestrial components of the atmospheric 14C variation spectrum", in The Sun in Time, (Eds.) Sonett, C.P., Giampapa, M.S., Matthews, M.S., pp. 360-388, University of Arizona Press, Tucson. [ADS] (Cited on pages 17, 30, and 51.)

Damon, P.E., Lerman, J.C. and Long, A., 1978, "Temporal Fluctuations of Atmospheric 14C: Causal Factors and Implications", Annu. Rev. Earth Planet. Sci., 6, 457-494. [DOI], [ADS] (Cited on pages 31 and 33.)

de Carvalho, J.X. and Prado, C.P.C., 2000, "Self-Organized Criticality in the Olami-Feder-Christensen Model", Phys. Rev. Lett., 84, 4006-4009. [DOI], [ADS] (Cited on page 58.)

de Jager, C., 2005, "Solar Forcing of Climate. 1: Solar Variability", Space Sci. Rev., 120, 197-241. [DOI], [ADS] (Cited on page 59.)

de Jong, A.F.M. and Mook, W.G., 1982, "An anomalous Suess effect above Europe", Nature, 298, 641-644. [DOI], [ADS] (Cited on page 32.)

de Meyer, F., 1998, "Modulation of the Solar Magnetic Cycle", Solar Phys., 181, 201-219. [DOI], [ADS] (Cited on pages 15, 16, and 56.)

Delaygue, G. and Bard, E., 2011, "An Antarctic view of Beryllium-10 and solar activity for the past millennium", Climate Dyn, 36, 2201-2218. [DOI] (Cited on page 47.)

Dikpati, M. and Gilman, P.A., 2006, "Simulating and Predicting Solar Cycles Using a Flux-Transport Dynamo", Astrophys. J., 649, 498-514. [DOI], [ADS] (Cited on page 19.)

Dikpati, M., de Toma, G. and Gilman, P.A., 2008, "Polar Flux, Cross-equatorial Flux, and Dynamogenerated Tachocline Toroidal Flux as Predictors of Solar Cycles", Astrophys. J., 675, 920-930. [DOI], [ADS] (Cited on page 19.)

Donadini, F., Korte, M. and Constable, C., 2010, "Millennial Variations of the Geomagnetic Field: from Data Recovery to Field Reconstruction", Space Sci. Rev., 155, 219-246. [DOI], [ADS] (Cited on page 24.)

Donnelly, R.F., White, O.R. and Livingston, W.C., 1994, "The solar Ca II K index and the Mg II core-to-wing ratio", Solar Phys., 152, 69-76. [DOI], [ADS] (Cited on page 13.)

Dorman, L.I., 2004, Cosmic Rays in the Earth's Atmosphere and Underground, Kluwer Academic Publishers, Dordrecht. [Google Books] (Cited on page 25.)

Dreschhoff, G.A.M. and Zeller, E.J., 1990, "Evidence of individual solar proton events in Antarctic snow", Solar Phys., 127, 333-346. [DOI], [ADS] (Cited on page 66.)

Dreschhoff, G.A.M. and Zeller, E.J., 1998, "Ultra-High Resolution Nitrate in Polar Ice as Indicator of Past Solar Activity", Solar Phys., 177, 365-374. [DOI], [ADS] (Cited on page 66.)

Eddy, J.A., 1976, "The Maunder Minimum", Science, 192, 1189-1202. [DOI], [ADS] (Cited on pages 16 and 52.)

Eddy, J.A., 1977a, "The case of the missing sunspots", Scientific American, 236, 80-88. [ADS] (Cited on pages 53, 54, and 56.)

Eddy, J.A., 1977b, "Climate and the changing sun", Climatic Change, 1, 173-190. [DOI] (Cited on pages 53, 54, and 56.)

Eddy, J.A., 1983, "The Maunder Minimum: A reappraisal", Solar Phys., 89, 195-207. [DOI], [ADS] (Cited on pages 15 and 16.)

Eddy, J.A., Stephenson, F.R. and Yau, K.K.C., 1989, "On pre-telescopic sunspot records", Quart. J. R. Astron. Soc., 30, 65-73. [ADS] (Cited on page 15.)

Elsasser, W., Nay, E.P. and Winkler, J.R., 1956, "Cosmic-ray intensity and geomagnetism", Nature, 178, 1226-1227. [DOI], [ADS] (Cited on page 33.)

Feynman, J. and Gabriel, S.B., 1990, "Period and phase of the 88-year solar cycle and the Maunder minimum: Evidence for a chaotic Sun", Solar Phys., 127, 393-403. [DOI], [ADS] (Cited on pages 17 and 18.)

Field, C.V., Schmidt, G.A., Koch, D. and Salyk, C., 2006, "Modeling production and climate-related impacts on 10Be concentration in ice cores", J. Geophys. Res., 111, D15107. [DOI], [ADS] (Cited on pages 36 and 37.)

Fink, D., Klein, J., Middleton, R., Vogt, S., Herzog, G.F. and Reedy, R.C., 1998, "41Ca, 26Al, and 10Be in lunar basalt 74275 and 10Be in the double drive tube 74002/74001", Geochim. Cosmochim. Acta,, 62, 2389-2402. [DOI], [ADS] (Cited on page 66.)

Finkel, R.C. and Nishiizumi, K., 1997, "Beryllium 10 concentrations in the Greenland Ice Sheet Project 2 ice core from 3-40 ka", J. Geophys. Res., 102, 26,699-26,706. [DOI], [ADS] (Cited on page 34.)

Foukal, P., 1996, "The behavior of solar magnetic plages measured from Mt. Wilson observations between 1915-1984", Geophys. Res. Lett., 23, 2169-2172. [DOI], [ADS] (Cited on page 13.)

Foukal, P., Fröhlich, C., Spruit, H. and Wigley, T.M.L., 2006, "Variations in solar luminosity and their effect on the Earth's climate", Nature, 443, 161-166. [DOI], [ADS] (Cited on page 59.)

Fraser-Smith, A.C., 1987, "Centered and eccentric geomagnetic dipoles and their poles, 1600-1985", Rev. Geophys., 25, 1-16. [DOI], [ADS] (Cited on page 24.)

Frick, P., Galyagin, D., Hoyt, D.V., Nesme-Ribes, E., Schatten, K.H., Sokoloff, D. and Zakharov, V., 1997, "Wavelet analysis of solar activity recorded by sunspot groups", Astron. Astrophys., 328, 670-681. [ADS] (Cited on pages 16, 17, and 52.)

Fröhlich, C., 2006, "Solar Irradiance Variability since 1978", Space Sci. Rev., 125, 53-65. [ADS] (Cited on page 59.)

Fröhlich, C., 2012, "Total Solar Irradiance Observations", Surv. Geophys., 33, 453-473. [DOI], [ADS] (Cited on page 13.)

Garcia-Munoz, M., Mason, G.M. and Simpson, J.A., 1975, "The anomalous 4He component in the cosmic-ray spectrum at <50 MeV per nucleon during 1972-1974", Astrophys. J., 202, 265-275. [DOI], [ADS] (Cited on page 23.)

Genevey, A., Gallet, Y., Constable, C.G., Korte, M. and Hulot, G., 2008, "ArcheoInt: An upgraded compilation of geomagnetic field intensity data for the past ten millennia", Geochem. Geophys. Geosyst., 9, Q04038. [DOI] (Cited on pages 24 and 25.)

Gerber, S., Joos, F., Brögger, P., Stocker, T.F., Mann, M.E., Sitch, S. and Scholze, M., 2002, "Constraining temperature variations over the last millennium by comparing simulated and observed atmospheric CO2", Climate Dyn, 20, 281-299. [ADS] (Cited on page 31.)

Gibson, S.E., de Toma, G., Emery, B., Riley, P., Zhao, L., Elsworth, Y., Leamon, R.J., Lei, J., McIntosh, S., Mewaldt, R.A., Thompson, B.J. and Webb, D., 2011, "The Whole Heliosphere Interval in the Context of a Long and Structured Solar Minimum: An Overview from Sun to Earth", Solar Phys., 274, 5-27. [DOI], [ADS] (Cited on pages 7 and 56.)

Gleeson, L.J. and Axford, W.I., 1968, "Solar Modulation of Galactic Cosmic Rays", Astrophys. J., 154, 1011-1026. [DOI], [ADS] (Cited on page 22.)

Gleissberg, W., 1939, "A long-periodic fluctuation of the sun-spot numbers", Observatory, 62, 158-159. [ADS] (Cited on page 17.)

Gleissberg, W., 1948, "A preliminary forecast of solar activity", Popular Astron., 56, 399. [ADS] (Cited on page 18.)

Gleissberg, W., 1952, Die Häufigkeit der Sonnenflecken, Akademie-Verlag, Berlin (Cited on page 16.)

Gleissberg, W., 1971, "The Probable Behaviour of Sunspot Cycle 21", Solar Phys., 21, 240-245. [DOI], [ADS] (Cited on page 17.)

Goslar, T., 2001, "Absolute production of radiocarbon and the long-term trend of atmospheric radiocarbon", Radiocarbon, 43, 743-749 (Cited on page 31.)

Goslar, T., 2003, "14C as an Indicator of Solar Variability", PAGES News, 11(2/3), 12-14. URL (accessed

26 February 2008): (Cited on pages 53 and 54.)

Gray, L.J., Beer, J., Geller, M., Haigh, J.D., Lockwood, M., Matthes, K., Cubasch, U., Fleitmann, D., Harrison, G., Hood, L., Luterbacher, J., Meehl, G.A., Shindell, D., van Geel, B. and White, W., 2010, "Solar Influences on Climate", Rev. Geophys, 48, RG4001. [DOI], [ADS] (Cited on pages 7 and 59.)

Grieder, P.K.F., 2001, Cosmic Rays at Earth: Researcher's Reference Manual and Data Book, Elsevier Science, Amsterdam; New York (Cited on page 45.)

Grismore, R., Llewellyn, R.A., Brown, M.D., Dowson, S.T. and Cumblidge, K., 2001, "Measurements of the concentrations of 26Al in lunar rocks 15555 and 60025", Earth Planet. Sci. Lett., 187, 163-171. [DOI], [ADS] (Cited on page 66.)

Haigh, J.D., 2007, "The Sun and the Earth's Climate", Living Rev. Solar Phys., 4, lrsp-2007-2. [ADS]. URL (accessed 25 February 2008): (Cited on pages 7 and 59.)

Haigh, J.D., Winning, A.R., Toumi, R. and Harder, J.W., 2010, "An influence of solar spectral variations on radiative forcing of climate", Nature, 467, 696-699. [DOI], [ADS] (Cited on page 59.)

Hale, G.E., Ellerman, F., Nicholson, S.B. and Joy, A.H., 1919, "The Magnetic Polarity of Sun-Spots", Astrophys. J, 49, 153-178. [DOI], [ADS] (Cited on page 16.)

Hanslmeier, A., Brajsa, R., Calogovic, J., Vrsnak, B., RuZdjak, D., Steinhilber, F., MacLeod, C.L., Ivezic, Z. and Skokic, I., 2013, "The chaotic solar cycle. II. Analysis of cosmogenic 10Be data", Astron. Astrophys., 550, A6. [DOI], [ADS] (Cited on page 18.)

Harvey, K.L. and White, O.R., 1999, "What is solar cycle minimum?", J. Geophys. Res., 104(A9), 19,75919,764. [DOI], [ADS] (Cited on page 52.)

Hathaway, D.H., 2009, "Solar Cycle Forecasting", Space Sci. Rev., 144, 401-412. [DOI], [ADS] (Cited on pages 18 and 19.)

Hathaway, D.H., 2010, "The Solar Cycle", Living Rev. Solar Phys., 7, lrsp-2010-1. [ADS]. URL (accessed 25 March 2010): (Cited on pages 9 and 16.)

Hathaway, D.H. and Wilson, R.M., 2004, "What the Sunspot Record Tells Us About Space Climate", Solar Phys., 224, 5-19. [DOI] (Cited on pages 9 and 11.)

Hathaway, D.H., Wilson, R.M. and Reichmann, E.J., 1999, "A Synthesis of Solar Cycle Prediction Techniques", J. Geophys. Res., 104, 22,375-22,388. [DOI], [ADS] (Cited on page 19.)

Heikkila, U., Beer, J. and Feichter, J., 2008, "Modeling cosmogenic radionuclides 10Be and 7Be during the Maunder Minimum using the ECHAM5-HAM General Circulation Model", Atmos. Chem. Phys., 8, 2797-2809. [DOI], [ADS]. URL (accessed 11 September 2012): (Cited on page 36.)

Heikkila, U., Beer, J. and Feichter, J., 2009, "Meridional transport and deposition of atmospheric 10Be", Atmos. Chem. Phys., 9, 515-527. [DOI], [ADS]. URL (accessed 11 September 2012): (Cited on page 36.)

Herbst, K., Kopp, A., Heber, B., Steinhilber, F., Fichtner, H., Scherer, K. and Matthia, D., 2010, "On the importance of the local interstellar spectrum for the solar modulation parameter", J. Geophys. Res., 115, D00I20. [DOI], [ADS] (Cited on page 23.)

Hongre, L., Hulot, G. and Khokhlov, A., 1998, "An analysis of the geomagnetic field over the past 2000 years", Phys. Earth Planet Inter., 106, 311-335. [DOI], [ADS] (Cited on page 25.)

Horiuchi, K., Minoura, K., Kobayashi, K., Nakamura, T., Hatori, S., Matsuzaki, H. and Kawai, T., 1999, "Last-glacial to post-glacial 10Be fluctuations in a sediment core from the Academician Ridge, Lake Baikal", Geophys. Res. Lett., 26, 1047-1050. [DOI], [ADS] (Cited on page 34.)

Horiuchi, K., Ohta, A., Uchida, T., Matsuzaki, H., Shibata, Y. and Motoyama, H., 2007, "Concentration of 10Be in an ice core from the Dome Fuji station, Eastern Antarctica: Preliminary results from 1500 to 1810 yr AD", Nucl. Instrum. Methods B, 259, 584-587. [DOI], [ADS] (Cited on page 34.)

Horiuchi, K., Uchida, T., Sakamoto, Y., Ohta, A., Matsuzaki, H., Shibata, Y. and Motoyama, H., 2008, "Ice core record of 10Be over the past millennium from Dome Fuji, Antarctica: A new proxy record of past solar activity and a powerful tool for stratigraphic dating", Quat. Geochronology, 3, 253-261. [DOI] (Cited on pages 34, 36, and 63.)

Hoyng, P., 1993, "Helicity fluctuations in mean field theory: an explanation for the variability of the solar cycle?", Astron. Astrophys., 272, 321-339. [ADS] (Cited on pages 18 and 54.)

Hoyt, D.V. and Schatten, K., 1998, "Group Sunspot Numbers: A New Solar Activity Reconstruction", Solar Phys., 179, 189-219. [ADS] (Cited on pages 11, 14, and 43.)

Hoyt, D.V. and Schatten, K.H., 1996, "How Well Was the Sun Observed during the Maunder Minimum?", Solar Phys., 165, 181-192. [DOI], [ADS] (Cited on pages 11 and 52.)

Hoyt, D.V., Schatten, K.H. and Nesme-Ribes, E., 1994, "The one hundredth year of Rudolf Wolf's death: Do we have the correct reconstruction of solar activity?", Geophys. Res. Lett., 21, 2067-2070. [DOI], [ADS] (Cited on page 11.)

Hudson, H. S., 2010, "Solar flares add up", Nature Phys., 6, 637-638. [DOI], [ADS] (Cited on pages 62 and 67.)

Jennings, R.L. and Weiss, N.O., 1991, "Symmetry breaking in stellar dynamos", Mon. Not. R. Astron. Soc., 252, 249-260. [ADS] (Cited on pages 54 and 58.)

Jiang, J., Chatterjee, P. and Choudhuri, A.R., 2007, "Solar activity forecast with a dynamo model", Mon. Not. R. Astron. Soc., 381, 1527-1542. [DOI], [ADS], [arXiv:0707.2258 [astro-ph]] (Cited on page 19.)

Joos, F., 1994, "Imbalance in the budget", Nature, 370, 181-182. [DOI] (Cited on page 33.)

Jull, A.J.T., Cloudt, S., Donahue, D.J., Sisterson, J.M., Reedy, R.C. and Masarik, J., 1998, "14C depth profiles in Apollo 15 and 17 cores and lunar rock 68815", Geochim. Cosmochim. Acta, 62, 3025-3036. [DOI], [ADS] (Cited on pages 65 and 66.)

Kane, R.P., 2007, "A Preliminary Estimate of the Size of the Coming Solar Cycle 24, based on Ohl's Precursor Method", Solar Phys., 243, 205-217. [DOI], [ADS] (Cited on page 19.)

Käpylä, P.J., Mantere, M.J. and Brandenburg, A., 2012, "Cyclic Magnetic Activity due to Turbulent Convection in Spherical Wedge Geometry", Astrophys. J. Lett., 755, L22. [DOI], [ADS] (Cited on page 17.)

Karak, B.B. and Nandy, D., 2012, "Turbulent Pumping of Magnetic Flux Reduces Solar Cycle Memory and thus Impacts Predictability of the Sun's Activity", Astrophys. J. Lett., 761, L13. [DOI], [ADS], [arXiv:1206.2106 [astro-ph.SR]] (Cited on page 19.)

Kigoshi, K. and Hasegawa, H., 1966, "Secular variation of atmospheric radiocarbon concentration and its dependence on geomagnetism", J. Geophys. Res., 71, 1065-1071. [ADS] (Cited on page 34.)

Klecker, B., Kunow, H., Cane, H.V., Dalla, S., Heber, B., Kecskemety, K., Klein, K.-L., Kota, J., Kucharek, H., Lario, D., Lee, M.A., Popecki, M.A., Posner, A., Rodriguez-Pacheco, J., Sanderson, T., Simnett, G.M. and Roelof, E.C., 2006, "Energetic Particle Observations", Space Sci. Rev., 123, 217-250. [DOI], [ADS] (Cited on page 61.)

Kleczek, J., 1952, "Solar flare index calculations", Publ. Centr. Astron. Inst. Czechoslovakia,, 22 (Cited on page 12.)

Koch, D., Schmidt, G.A. and Field, C.V., 2006, "Sulfur, sea salt, and radionuclide aerosols in GISS ModelE", J. Geophys. Res., 111(D10), D06206. [DOI], [ADS] (Cited on page 36.)

Kocharov, G.E., Ostryakov, V.M., Peristykh, A.N. and Vasil'ev, V.A., 1995, "Radiocarbon Content Variations and Maunder Minimum of Solar Activity", Solar Phys., 159, 381-391. [DOI], [ADS] (Cited on page 52.)

Kohl, C.P., Murrell, M.T., Russ III, G.P. and Arnold, J.R., 1978, "Evidence for the constancy of the solar cosmic ray flux over the past ten million years: 53Mn and 26Al measurements", in Lunar and Planetary Science IX, Proceedings of the conference, Houston, TX, March 13-17, 1978, Geochim. Cosmochim. Acta Suppl., 10, pp. 2299-2310, Pergamon Press, New York. [ADS] (Cited on page 66.)

Korte, M. and Constable, C.G., 2005, "The geomagnetic dipole moment over the last 7000 years - new results from a global model", Earth Planet. Sci. Lett., 236, 348-358. [DOI], [ADS] (Cited on pages 24,

25, 43, 44, and 56.)

Korte, M. and Constable, C.G., 2006, "Centennial to millennial geomagnetic secular variation", Geophys. J. Int., 167, 43-52. [DOI], [ADS] (Cited on page 24.)

Korte, M. and Constable, C.G., 2008, "Spatial and temporal resolution of millennial scale geomagnetic field models", Adv. Space Res., 41, 57-69. [DOI], [ADS] (Cited on page 24.)

Korte, M., Constable, C., Donadini, F. and Holme, R., 2011, "Reconstructing the Holocene geomagnetic field", Earth Planet. Sci. Lett., 312, 497-505. [DOI], [ADS] (Cited on pages 24 and 38.)

Kovaltsov, G.A. and Usoskin, I.G., 2007, "Regional cosmic ray induced ionization and geomagnetic field changes", Adv. Geosci, 13, 31-35. [ADS]. URL (accessed 14 October 2008): (Cited on page 15.)

Kovaltsov, G.A. and Usoskin, I.G., 2010, "A new 3D numerical model of cosmogenic nuclide 10 Be production in the atmosphere", Earth Planet. Sci. Lett., 291, 182-188. [DOI], [ADS] (Cited on pages 25 and 35.)

Kovaltsov, G.A., Usoskin, I.G. and Mursula, K., 2004, "An Upper Limit on Sunspot Activity During the Maunder Minimum", Solar Phys., 224, 95-101. [DOI], [ADS] (Cited on page 52.)

Kovaltsov, G.A., Mishev, A. and Usoskin, I.G., 2012, "A new model of cosmogenic production of radiocarbon 14C in the atmosphere", Earth Planet. Sci. Lett., 337, 114-120. [DOI], [ADS], [arXiv:1206.6974 []] (Cited on pages 25, 29, 30, and 31.)

Kremliovsky, M.N., 1994, "Can we understand time scales of solar activity?", Solar Phys., 151, 351-370. [DOI], [ADS] (Cited on pages 16 and 17.)

Kremliovsky, M.N., 1995, "Limits of Predictability of Solar Activity", Solar Phys., 159, 371-380. [DOI], [ADS] (Cited on pages 18 and 19.)

Krivova, N.A., Balmaceda, L. and Solanki, S.K., 2007, "Reconstruction of solar total irradiance since 1700 from the surface magnetic flux", Astron. Astrophys., 467, 335-346. [DOI], [ADS] (Cited on pages 13, 42, and 44.)

Krivsky, L., 1984, "Long-term fluctuations of solar activity during the last thousand years", Solar Phys., 93, 189-194. [ADS] (Cited on page 15.)

Krivsky, L. and Pejml, K., 1988, "Solar Activity, Aurorae and Climate in Central Europe in the Last 1000 Years", Publ. Astron. Inst. Czechoslovak Acad. Sci., 75, 32 (Cited on page 52.)

Kuklin, G.V., 1976, "Cyclical and Secular Variations of Solar Activity", in Basic Mechanisms of Solar Activity, Symposium no. 71 held in Prague, Czechoslovakia, 25-29 August 1975, (Eds.) Bumba, V., Kleczek, J., IAU Symposium, 71, pp. 147-148, D. Reidel, Dordrecht; Boston. [ADS] (Cited on page 17.)

Kurths, J. and Ruzmaikin, A.A., 1990, "On forecasting the sunspot numbers", Solar Phys., 126, 407-410. [DOI], [ADS] (Cited on pages 16, 18, and 19.)

Lal, D., 1987, "10Be in polar ice: data reflect changes in cosmic ray flux or polar meteorology", Geophys. Res. Lett., 14, 785-788. [DOI], [ADS] (Cited on page 36.)

Lal, D. and Peters, B., 1967, "Cosmic Ray Produced Radioactivity on the Earth", in Kosmische Strahlung II / Cosmic Rays II, (Ed.) Sittle, K., Handbuch der Physik / Encyclopedia of Physics, 9/46/2, pp. 551-612, Springer, Berlin. [DOI], [ADS] (Cited on page 35.)

Lal, D. and Suess, H.E., 1968, "The radioactivity of the atmosphere and hydrosphere", Annu. Rev. Nucl. Sci., 18, 407-434. [DOI] (Cited on page 35.)

Lal, D., Jull, A.J.T., Pollard, D. and Vacher, L., 2005, "Evidence for large century time-scale changes in solar activity in the past 32 Kyr, based on in-situ cosmogenic 14C in ice at Summit, Greenland", Earth Planet. Sci. Lett., 234, 335-349. [DOI], [ADS] (Cited on page 28.)

Lawrence, J.K., Cadavid, A.C. and Ruzmaikin, A.A., 1995, "Turbulent and Chaotic Dynamics Underlying Solar Magnetic Variability", Astrophys. J., 455, 366. [DOI], [ADS] (Cited on page 18.)

Lee, E.H., Ahn, Y.S., Yang, H.J. and Chen, K.Y., 2004, "The Sunspot and Auroral Activity Cycle Derived from Korean Historical Records of the 11th 18th Century", Solar Phys., 224, 373-386. [DOI], [ADS] (Cited on page 15.)

Letfus, V., 1999, "Daily relative sunspot numbers 1749-1848: reconstruction of missing observations", Solar Phys., 184, 201-211. [DOI], [ADS] (Cited on page 11.)

Letfus, V., 2000, "Relative sunspot numbers in the first half of eighteenth century", Solar Phys., 194, 175-184. [ADS] (Cited on page 11.)

Li, K.J., Yun, H.S. and Gu, X.M., 2001, "Latitude Migration of Sunspot Groups", Astron. J., 122, 21152117. [DOI], [ADS] (Cited on page 19.)

Lingenfelter, R.E., 1963, "Production of Carbon 14 by Cosmic-Ray Neutrons", Rev. Geophys. Space Phys., 1, 35-55. [DOI], [ADS] (Cited on page 25.)

Lingenfelter, R.E. and Hudson, H.S., 1980, "Solar particle fluxes and the ancient sun", in The Ancient Sun: Fossil Record in the Earth, Moon and Meteorites, Proceedings of the conference held at Boulder, Colorado, October 16-19, 1979, (Eds.) Pepin, R.O., Eddy, J.A., Merrill, R.B., pp. 69-79, Pergamon Press, New York; Oxford. [ADS] (Cited on page 62.)

Lockwood, M., Stamper, R. and Wild, M.N., 1999, "A doubling of the Sun's coronal magnetic field during the past 100 years", Nature, 399, 437-439. [DOI], [ADS] (Cited on pages 13 and 46.)

Lockwood, M., Owens, M.J., Barnard, L., Davis, C.J. and Steinhilber, F., 2011, "The persistence of solar activity indicators and the descent of the Sun into Maunder Minimum conditions", Geophys. Res. Lett., 38, L22105. [DOI], [ADS] (Cited on pages 19 and 56.)

Love, J.J., 2011, "Secular trends in storm-level geomagnetic activity", Ann. Geophys., 29, 251-262. [DOI], [ADS] (Cited on page 13.)

Lukianova, R., Alekseev, G. and Mursula, K., 2009, "Effects of station relocation in the aa index", J. Geophys. Res., 114, A02105. [DOI], [ADS] (Cited on page 13.)

Masarik, J. and Beer, J., 1999, "Simulation of particle fluxes and cosmogenic nuclide production in the Earth's atmosphere", J. Geophys. Res., 104(D10), 12,099-12,111. [DOI], [ADS] (Cited on pages 25, 31, 35, 36, and 43.)

Masarik, J. and Beer, J., 2009, "An updated simulation of particle fluxes and cosmogenic nuclide production in the Earth's atmosphere", J. Geophys. Res., 114(D13), D11103. [DOI], [ADS] (Cited on pages 25,

30, and 35.)

Matsumoto, K., Sarmiento, J.L., Key, R.M., Aumont, O., Bullister, J.L., Caldeira, K., Campin, J.-M., Doney, S.C., Drange, H., Dutay, J.-C., Follows, M., Gao, Y., Gnanadesikan, A., Gruber, N., Ishida, A., Joos, F., Lindsay, K., Maier-Reimer, E., Marshall, J.C., Matear, R.J., Monfray, P., Mouchet, A., Najjar, R., Plattner, G.-K., Schlitzer, R., Slater, R., Swathi, P.S., Totterdell, I.J., Weirig, M.-F., Yamanaka, Y., Yool, A. and Orr, J.C., 2004, "Evaluation of ocean carbon cycle models with data-based metrics", Geophys. Res. Lett., 31, L07303. [DOI], [ADS] (Cited on page 33.)

Mayaud, P.-N., 1972, "The aa indices: A 100-year series characterizing the magnetic activity", J. Geophys. Res., 77, 6870-6874. [DOI], [ADS] (Cited on page 13.)

McCracken, K.G., 2004, "Geomagnetic and atmospheric effects upon the cosmogenic 10Be observed in polar ice", J. Geophys. Res., 109(A18), A04101. [DOI], [ADS] (Cited on pages 36 and 38.)

McCracken, K.G., 2007, "Heliomagnetic field near Earth, 1428-2005", J. Geophys. Res., 112(A11), A09106. [DOI], [ADS] (Cited on pages 42 and 44.)

McCracken, K.G. and Beer, J., 2007, "Long-term changes in the cosmic ray intensity at Earth, 1428-2005", J. Geophys. Res., 112(A11), A10101. [DOI], [ADS] (Cited on pages 41, 42, 45, and 56.)

McCracken, K.G., Dreschhoff, G.A.M., Zeller, E.J., Smart, D.F. and Shea, M.A., 2001, "Solar cosmic ray events for the period 1561-1994: 1. Identification in polar ice, 1561-1950", J. Geophys. Res., 106, 21,585-21,598. [DOI], [ADS] (Cited on pages 64 and 66.)

McCracken, K.G., McDonald, F.B., Beer, J., Raisbeck, G.M. and Yiou, F., 2004, "A phenomenological study of the long-term cosmic ray modulation, 850-1958 AD", J. Geophys. Res., 109(A18), A12103. [DOI], [ADS] (Cited on pages 41 and 63.)

McDonald, F.B., Webber, W.R. and Reames, D.V., 2010, "Unusual time histories of galactic and anomalous cosmic rays at 1 AU over the deep solar minimum of cycle 23/24", Geophys. Res. Lett., 37, L18101. [DOI], [ADS] (Cited on page 22.)

McHargue, L.R. and Damon, P.E., 1991, "The Global Beryllium 10 Cycle", Rev. Geophys, 29, 141-158. [DOI], [ADS] (Cited on page 35.)

Mega, M.S., Allegrini, P., Grigolini, P., Latora, V., Palatella, L., Rapisarda, A. and Vinciguerra, S., 2003, "Power-Law Time Distribution of Large Earthquakes", Phys. Rev. Lett., 90, 188501. [DOI], [ADS], [arXiv:cond-mat/0212529] (Cited on page 58.)

Mendoza, B., 1997, "Geomagnetic activity and wind velocity during the Maunder minimum", Ann. Geophys., 15, 397-402. [DOI], [ADS]. URL (accessed 10 October 2008): (Cited on page 52.)

Mewaldt, R.A., 2006, "Solar energetic particle composition, energy spectra, and space weather", Space Sci. Rev, 124, 303-316. [DOI], [ADS] (Cited on page 61.)

Mewaldt, R.A., Cohen, C.M.S., Mason, G.M., Haggerty, D.K. and Desai, M.I., 2007, "Long-Term Fluences of Solar Energetic Particles from H to Fe", Space Sci. Rev., 130, 323-328. [DOI], [ADS] (Cited on page 61.)

Mewaldt, R.A., Looper, M.D., Cohen, C.M.S., Haggerty, D.K., Labrador, A.W., Leske, R.A., Mason, G.M., Mazur, J.E. and von Rosenvinge, T.T., 2012, "Energy Spectra, Composition, and Other Properties of Ground-Level Events During Solar Cycle 23", Space Sci. Rev., 171, 97-120. [DOI], [ADS] (Cited on page 61.)

Michel, R. and Neumann, S., 1998, "Interpretation of cosmogenic nuclides in meteorites on the basis of accelerator experiments and physical model calculations", Proc. Indian Acad. Sci. (Earth Planet. Sci.), 107, 441-457. [ADS] (Cited on page 46.)

Michel, R., Leya, I. and Borges, L., 1996, "Production of cosmogenic nuclides in meteoroids: accelerator experiments and model calculations to decipher the cosmic ray record in extraterrestrial matter", Nucl. Instrum.. Methods B, 113, 434-444. [DOI], [ADS] (Cited on page 66.)

Mikaloff Fletcher, S.E., Gruber, N., Jacobson, A.R., Doney, S.C., Dutkiewicz, S., Gerber, M., Follows, M., Joos, F., Lindsay, K., Menemenlis, D., Mouchet, A., Muller, S.A. and Sarmiento, J.L., 2006, "Inverse estimates of anthropogenic CO2 uptake, transport, and storage by the ocean", Global Biogeochem. Cycles, 20, GB2002. [DOI], [ADS] (Cited on page 33.)

Mininni, P.D., Gomez, D.O. and Mindlin, G.B., 2000, "Stochastic Relaxation Oscillator Model for the Solar Cycle", Phys. Rev. Lett., 85, 5476-5479. [DOI], [ADS] (Cited on page 18.)

Mininni, P.D., Gomez, D.O. and Mindlin, G.B., 2001, "Simple Model of a Stochastically Excited Solar Dynamo", Solar Phys., 201, 203-223. [DOI], [ADS] (Cited on pages 16, 18, 52, and 58.)

Mininni, P.D., Gomez, D.O. and Mindlin, G.B., 2002, "Biorthogonal Decomposition Techniques Unveil the Nature of the Irregularities Observed in the Solar Cycle", Phys. Rev. Lett., 89(6), 061101. [ADS] (Cited on page 17.)

Miyahara, H., Masuda, K., Muraki, Y., Furuzawa, H., Menjo, H. and Nakamura, T., 2004, "Cyclicity of Solar Activity During the Maunder Minimum Deduced from Radiocarbon Content", Solar Phys., 224, 317-322. [DOI], [ADS] (Cited on page 53.)

Miyahara, H., Masuda, K., Muraki, Y., Kitagawa, H. and Nakamura, T., 2006a, "Variation of solar cyclicity during the Spoerer Minimum", J. Geophys. Res., 111(A10), A03103. [DOI], [ADS] (Cited on page 53.)

Miyahara, H., Sokoloff, D. and Usoskin, I.G., 2006b, "The Solar Cycle at the Maunder Minimum Epoch", in Advances in Geosciences, Vol. 2: Solar Terrestrial (ST), (Eds.) Ip, W.-H., Duldig, M., pp. 1-20, World Scientific, Singapore; Hackensack, NJ. [ADS] (Cited on pages 52 and 53.)

Miyahara, H., Kitazawa, K., Nagaya, K., Yokoyama, Y., Matsuzaki, H., Masuda, K., Nakamura, T. and Muraki, Y., 2010, "Is the Sun heading for another Maunder Minimum? Precursors of the grand solar minima", J. Cosmol., 8, 1970-1982. [ADS] (Cited on page 19.)

Miyake, F., Nagaya, K., Masuda, K. and Nakamura, T., 2012, "A signature of cosmic-ray increase in ad 774-775 from tree rings in Japan", Nature, 486, 240-242. [DOI], [ADS] (Cited on pages 62 and 63.)

Moraal, H. and Stoker, P.H., 2010, "Long-term neutron monitor observations and the 2009 cosmic ray maximum", J. Geophys. Res., 115, A12109. [DOI], [ADS] (Cited on page 22.)

Morfill, G.E., Scheingraber, H., Voges, W. and Sonett, C.P., 1991, "Sunspot number variations - Stochastic or chaotic", in The Sun in Time, (Eds.) Sonett, C.P., Giampapa, M.S., Matthews, M.S., pp. 30-58, University of Arizona Press, Tucson. [ADS] (Cited on page 18.)

Moss, D., Brandenburg, A., Tavakol, R. and Tuominen, I., 1992, "Stochastic effects in mean-field dynamos", Astron. Astrophys., 265, 843-849. [ADS] (Cited on page 18.)

Moss, D., Sokoloff, D., Usoskin, I.G. and Tutubalin, V., 2008, "Solar Grand Minima and Random Fluctuations in Dynamo Parameters", Solar Phys., 250, 221-234. [DOI], [ADS], [arXiv:0806.3331] (Cited on pages 17, 18, 52, 55, and 58.)

Mossman, J.E., 1989, "A comprehensive search for sunspots without the aid of a telescope, 1981-1982", Quart. J. R. Astron. Soc., 30, 59-64. [ADS] (Cited on page 15.)

Mundt, M.D., Maguire II, W.B. and Chase, R.R.P., 1991, "Chaos in the Sunspot Cycle: Analysis and Prediction", J. Geophys. Res., 96, 1705-1716. [DOI], [ADS] (Cited on pages 17 and 18.)

Mursula, K., Usoskin, I.G. and Kovaltsov, G.A., 2001, "Persistent 22-year cycle in sunspot activity: Evidence for a relic solar magnetic field", Solar Phys., 198, 51-56. [DOI], [ADS] (Cited on page 52.)

Mursula, K., Usoskin, I.G. and Kovaltsov, G.A., 2003, "Reconstructing the long-term cosmic ray intensity: linear relations do not work", Ann. Geophys., 21, 863-867. [ADS]. URL (accessed 10 October 2008): (Cited on page 40.)

Muscheler, R., Joos, F., Muller, S.A. and Snowball, I., 2005, "Climate: How unusual is today's solar activity?", Nature, 436, E3-E4. [DOI], [ADS] (Cited on pages 33, 41, 43, 46, and 56.)

Muscheler, R., Joos, F., Beer, J., Muller, S.A., Vonmoos, M. and Snowball, I., 2007, "Solar activity during the last 1000 yr inferred from radionuclide records", Quat. Sci. Rev., 26, 82-97. [DOI], [ADS] (Cited on pages 34, 38, 41, and 47.)

Nagaya, K., Kitazawa, K., Miyake, F., Masuda, K., Muraki, Y., Nakamura, T., Miyahara, H. and Mat-suzaki, H., 2012, "Variation of the Schwabe Cycle Length During the Grand Solar Minimum in the 4th Century BC Deduced from Radiocarbon Content in Tree Rings", Solar Phys., 280, 223-236. [DOI], [ADS] (Cited on page 53.)

Nagovitsyn, Y.A., 1997, "A nonlinear mathematical model for the solar cyclicity and prospects for reconstructing the solar activity in the past", Astron. Lett., 23, 742-748. [ADS] (Cited on page 15.)

Nevanlinna, H., 1995, "Auroral observations in Finland - visual sightings during the 18th and 19th centuries", J. Geomag. Geoelectr., 47, 953-960 (Cited on page 13.)

Nevanlinna, H., 2004a, "Historical Space Climate Data from Finland: Compilation and Analysis", Solar Phys., 224, 395-405. [DOI], [ADS] (Cited on page 13.)

Nevanlinna, H., 2004b, "Results of the Helsinki magnetic observatory 1844-1912", Ann. Geophys., 22, 1691-1704. [ADS]. URL (accessed 14 October 2008): (Cited on page 13.)

Newton, H.W., 1928, "The Sun's cycle of activity", Quart. J. R. Meteorol. Soc., 54(227), 161-174 (Cited on page 18.)

Ney, E.P., 1959, "Cosmic Radiation and the Weather", Nature, 183, 451-452. [DOI] (Cited on page 59.)

Nishiizumi, K., Arnold, J.R., Kohl, C.P., Caffee, M.W., Masarik, J. and Reedy, R.C., 2009, "Solar cosmic ray records in lunar rock 64455", Geochim. Cosmochim. Acta, 73, 2163-2176. [ADS] (Cited on page 66.)

Nishizumi, K., Caffee, M.W. and Arnold, J.R., 1997, "10Be from the active Sun", in Lunar and Planetary Science XXVIII, Proceedings of the conference, March 17-21, 1997, 1027, Lunar and Planetary Institute, Houston. [ADS] (Cited on page 66.)

O'Brien, K., 1979, "Secular variations in the production of cosmogenic isotopes in the earth's atmosphere", J. Geophys. Res., 84, 423-431. [DOI], [ADS] (Cited on pages 31 and 35.)

O'Brien, K. and Burke, G.D.P., 1973, "Calculated cosmic ray neutron monitor response to solar modulation of galactic cosmic rays", J. Geophys. Res., 78, 3013-3019. [DOI], [ADS] (Cited on page 25.)

Oeschger, H., Siegenthaler, U., Schotterer, U. and Gugelmann, A., 1974, "A box diffusion model to study the carbon dioxide exchange in nature", Tellus, 27, 168-192 (Cited on page 31.)

Ogurtsov, M.G., 2004, "New Evidence for Long-Term Persistence in the Sun's Activity", Solar Phys., 220, 93-105. [DOI], [ADS] (Cited on page 39.)

Ogurtsov, M.G., Nagovitsyn, Y.A., Kocharov, G.E. and Jungner, H., 2002, "Long-period cycles of the Sun's activity recorded in direct solar data and proxies", Solar Phys., 211, 371-394. [DOI], [ADS] (Cited on pages 17 and 51.)

Oguti, T. and Egeland, A., 1995, "Auroral occurrence in Norwegian archives", J. Geomag. Geoelectr., 47, 353-359 (Cited on page 15.)

Oliver, R. and Ballester, J.L., 1996, "Rescaled Range Analysis of the Asymmetry of Solar Activity", Solar Phys., 169, 215-224. [DOi], [ADS] (Cited on page 17.)

Oliver, R. and Ballester, J.L., 1998, "Is there memory in solar activity?", Phys. Rev. E, 58, 5650-5654. [DOI], [ADS] (Cited on page 18.)

Ossendrijver, M., 2003, "The solar dynamo", Astron. Astrophys. Rev., 11, 287-367. [DOI], [ADS] (Cited on page 18.)

Ossendrijver, M.A.J.H., 2000, "The dynamo effect of magnetic flux tubes", Astron. Astrophys., 359, 12051210. [ADS] (Cited on pages 52 and 58.)

Ossendrijver, M.A.J.H., Hoyng, P. and Schmitt, D., 1996, "Stochastic excitation and memory of the solar dynamo", Astron. Astrophys., 313, 938-948. [ADS] (Cited on pages 18 and 54.)

Ostriakov, V.M. and Usoskin, I.G., 1990, "On the dimension of solar attractor", Solar Phys., 127, 405-412. [DOI], [ADS] (Cited on pages 17 and 18.)

Ostryakov, V.M. and Usoskin, I.G., 1990, "Correlation dimensions of structured signals", Sov. Tech. Phys. Lett., 16, 658-659 (Cited on page 17.)

Owens, M.J., Usoskin, I.G. and Lockwood, M., 2012, "Heliospheric modulation of galactic cosmic rays during grand solar minima: Past and future variations", Geophys. Res. Lett., 39, L19102. [DOI], [ADS] (Cited on page 52.)

Ozguc, A., Atac, T. and Rybak, J., 2003, "Temporal variability of the flare index (1966-2001)", Solar Phys., 214, 375-396. [DOI], [ADS] (Cited on page 12.)

Panchev, S. and Tsekov, M., 2007, "Empirical evidences of persistence and dynamical chaos in solar terrestrial phenomena", J. Atmos. Sol.-Terr. Phys., 69, 2391-2404. [DOI], [ADS] (Cited on page 18.)

Parker, E.N., 1955, "Hydromagnetic Dynamo Models", Astrophys. J., 122, 293-314. [DOI], [ADS] (Cited on page 58.)

Parker, E.N., 1965, "The passage of energetic charged particles through interplanetary space", Planet. Space Sci., 13, 9-49. [DOI], [ADS] (Cited on page 21.)

Pedro, J.B., Heikkila, U.E., Klekociuk, A., Smith, A.M., van Ommen, T.D. and Curran, M.A.J., 2011, "Beryllium-10 transport to Antarctica: Results from seasonally resolved observations and modeling", J. Geophys. Res., 116, D23120. [DOI], [ADS] (Cited on page 36.)

Peristykh, A.N. and Damon, P.E., 1998, "Modulation of atmospheric 14C concentration by the solar wind and irradiance components of the Hale and Schwabe solar cycles", Solar Phys., 177, 343-343. [DOI], [ADS] (Cited on page 52.)

Peristykh, A.N. and Damon, P.E., 2003, "Persistence of the Gleissberg 88-year solar cycle over the last ^12,000 years: Evidence from cosmogenic isotopes", J. Geophys. Res., 108, 1003. [DOI], [ADS] (Cited on pages 17, 32, 34, and 50.)

Pesnell, W.D., 2012, "Solar Cycle Predictions (Invited Review)", Solar Phys., 281, 507-532. [DOI], [ADS] (Cited on page 18.)

Petrovay, K., 2010, "Solar Cycle Prediction", Living Rev. Solar Phys., 7, 6. [ADS], [arXiv:1012.5513 [astro-ph.SR]] (Cited on page 18.)

Pipin, V.V., Sokoloff, D.D. and Usoskin, I.G., 2012, "Variations of the solar cycle profile in a solar dynamo with fluctuating dynamo governing parameters", Astron. Astrophys., 542, A26. [DOI], [ADS], [arXiv:1112.6218 [astro-ph.SR]] (Cited on page 55.)

Polygiannakis, J., Preka-Papadema, P. and Moussas, X., 2003, "On signal-noise decomposition of time-series using the continuous wavelet transform: application to sunspot index", Mon. Not. R. Astron. Soc, 343, 725-734. [DOI], [ADS], [arXiv:physics/0301030] (Cited on page 17.)

Price, C.P., Prichard, D. and Hogenson, E.A., 1992, "Do the sunspot numbers form a 'chaotic' set?", J. Geophys. Res., 97, 19,113-19,120. [DOI], [ADS] (Cited on page 18.)

Pulkkinen, T., 2007, "Space Weather: Terrestrial Perspective", Living Rev. Solar Phys., 4, lrsp-2007-1. URL (accessed 25 February 2008): (Cited on page 13.)

Raisbeck, G.M., Yiou, F., Fruneau, M., Loiseaux, J.M., Lieuvin, M. and Ravel, J.C., 1981, "Cosmogenic

10Be/7Be as a probe of atmospheric transport processes", Geophys. Res. Lett., 8, 1015-1018. [DOI], [ADS] (Cited on page 35.)

Raisbeck, G.M., Yiou, F., Jouzel, J. and Petit, J.R., 1990, "10Be and 5 2H in Polar Ice Cores as a Probe of the Solar Variability's Influence on Climate", Philos. Trans. R. Soc. London, Ser. A, 330, 463-469. [ADS] (Cited on pages 34 and 63.)

Rao, M.N., Garrison, D.H., Bogard, D.D. and Reedy, R.C., 1994, "Determination of the flux and energy distribution of energetic solar protons in the past 2 Myr using lunar rock 68815", Geochim. Cosmochim.. Acta, 58, 4231-4245. [DOI], [ADS] (Cited on page 66.)

Reedy, R.C., 1996, "Constraints on Solar Particle Events from Comparisons of Recent Events and Million-Year Averages", in Solar Drivers of the Interplanetary and Terrestrial Disturbances, Proceedings of the 16th International Workshop, National Solar Observatory/Sacramento Peak, Sunspot, New Mexico, USA, 16-20 October 1995, (Eds.) Balasubramaniam, K.S., Keil, S.L., Smartt, R.N., ASP Conference Series, 95, pp. 429-436, Astronomical Society of the Pacific, San Francisco. [ADS] (Cited on page 65.)

Reedy, R.C., 2002, "Recent Solar Energetic Particles: Updates and Trends", in 33rd Lunar and Planetary Science Conference, League City, TX, March 11-15, 2002, 1938, Lunar and Planetary Institute, Houston. [ADS]. URL (accessed 5 March 2013): (Cited on page 65.)

Reedy, R.C., 2012, "Update on Solar-Proton Fluxes During the Last Five Solar Activity Cycles", in 43rd Lunar and Planetary Science Conference, The Woodlands, TX, March 19-23, 2012, 43, 1285, Lunar and Planetary Institute, Houston. [ADS]. URL (accessed 5 March 2013): (Cited on pages 61, 65, and 66.)

Reeves, G.D., Cayton, T.E., Gary, S.P. and Belian, R.D., 1992, "The great solar energetic particle events of 1989 observed from geosynchronous orbit", J. Geophys. Res., 97, 6219-6226. [DOI], [ADS] (Cited on page 61.)

Reimer, P.J., Baillie, M.G.L., Bard, E., Bayliss, A., Beck, J.W., Bertrand, C.J.H., Blackwell, P.G., Buck, C.E., Burr, G.S., Cutler, K.B., Damon, P.E., Edwards, R.L., Fairbanks, R.G., Friedrich, M., Guilderson, T.P., Hogg, A.G., Hughen, K.A., Kromer, B., McCormac, G., Manning, S., Ramsey, C.B., Reimer, R.W., Remmele, S., Southon, J.R., Stuiver, M., Talamo, S., Taylor, F.W., van der Plicht, J. and Weyhenmeyer, C.E., 2004, "IntCal04 terrestrial radiocarbon age calibration, 0-26 cal kyr BP", Radiocarbon, 46(3), 1029-1058. URL (accessed 11 September 2012): (Cited on pages 28 and 59.)

Reimer, P.J., Baillie, M.G.L., Bard, E., Bayliss, A., Beck, J.W., Blackwell, P.G., Ramsey, C. Bronk, Buck, C.E., Burr, G.S., Edwards, R.L., Friedrich, M., Grootes, P.M., Guilderson, T.P., Hajdas, I., Heaton, T.J., Hogg, A.G., Hughen, K.A., Kaiser, K.F., Kromer, B., McCormac, F.G., Manning, S.W., Reimer, R.W., Richards, D.A., Southon, J.R., Talamo, S., Turney, C.S.M., van der Plicht, J. and Weyhenmeye, C.E., 2009, "IntCal09 and Marine09 Radiocarbon Age Calibration Curves, 0-50,000 Years cal BP", Radiocarbon, 51(4), 1111-1150. URL (accessed 11 September 2012): (Cited on

pages 27, 59, and 63.)

Ribes, J.C. and Nesme-Ribes, E., 1993, "The solar sunspot cycle in the Maunder minimum AD1645 to AD1715", Astron. Astrophys, 276, 549-563. [ADS] (Cited on pages 11, 52, and 53.)

Rigozo, N.R., Echer, E., Vieira, L.E.A. and Nordemann, D.J.R., 2001, "Reconstruction of Wolf Sunspot Numbers on the Basis of Spectral Characteristics and Estimates of Associated Radio Flux and Solar Wind Parameters for the Last Millennium", Solar Phys., 203, 179-191. [DOI], [ADS] (Cited on page 15.)

Rozanov, E., Calisto, M., Egorova, T., Peter, T. and Schmutz, W., 2012, "Influence of the Precipitating Energetic Particles on Atmospheric Chemistry and Climate", Surv. Geophys., 33, 483-501. [DOI], [ADS] (Cited on page 38.)

Rozelot, J.P., 1994, "On the stability of the 11-year solar cycle period (and a few others)", Solar Phys., 149, 149-154. [DOI], [ADS] (Cited on page 17.)

Rozelot, J.P., 1995, "On the chaotic behaviour of the solar activity", Astron. Astrophys., 297, L45-L48. [ADS] (Cited on page 18.)

Ruth, U., Barnola, J.-M., Beer, J., Bigler, M., Blunier, T., Castellano, E., Fischer, H., Fundel, F., Huy-brechts, P., Kaufmann, P., Kipfstuhl, S., Lambrecht, A., Morganti, A., Oerter, H., Parrenin, F., Rybak, O., Severi, M., Udisti, R., Wilhelms, F. and Wolff, E.W., 2007, "'EDML1': a chronology for the EPICA deep ice core from Dronning Maud Land, Antarctica, over the last 150 000 years", Clim. Past, 3, 475484. [DOI], [ADS]. URL (accessed 11 September 2012): (Cited on page 34.)

Ruzmaikin, A., Feynman, J. and Robinson, P., 1994, "Long-term persistence of solar activity", Solar Phys., 149, 395-403. [DOI], [ADS] (Cited on page 18.)

Ruzmaikin, A.A., 1981, "The solar cycle as a strange attractor", Comments Astrophys., 9, 85-93. [ADS] (Cited on page 18.)

Rybansky, M., Rusin, V., Minarovjech, M., Klocok, L. and Cliver, E.W., 2005, "Reexamination of the coronal index of solar activity", J. Geophys. Res., 110(A9), 8106. [DOI], [ADS] (Cited on page 12.)

Sabine, C.L., Feely, R.A., Gruber, N., Key, R.M., Lee, K., Bullister, J.L., Wanninkhof, R., Wong, C.S., Wallace, D.W.R., Tilbrook, B., Millero, F.J., Peng, T.-H., Kozyr, A., Ono, T. and Rios, A.F., 2004, "The Oceanic Sink for Anthropogenic CO2", Science, 305, 367-371. [DOI], [ADS] (Cited on page 33.)

Salakhutdinova, I.I., 1999, "Identifying the quasi-regular and stochastic components of solar cyclicity and their properties", Solar Phys, 188, 377-396. [DOI], [ADS] (Cited on page 18.)

Scherer, K. and Fichtner, H., 2004, "Constraints on the heliospheric magnetic field variation during the Maunder Minimum from cosmic ray modulation modelling", Astron. Astrophys., 413, L11-L14. [DOI], [ADS] (Cited on page 40.)

Scherer, K., Fahr, H.-J., Fichtner, H. and Heber, B., 2004, "Long-Term Modulation of Cosmic Rays in the Heliosphere and its Influence at Earth", Solar Phys., 224, 305-316. [DOI], [ADS] (Cited on page 40.)

Schlamminger, L., 1990, "Aurora borealis during the Maunder minimum", Mon. Not. R. Astron. Soc., 247, 67-69. [ADS] (Cited on page 52.)

Schmalz, S. and Stix, M., 1991, "An aO dynamo with order and chaos", Astron. Astrophys., 245, 654-661. [ADS] (Cited on page 18.)

Schmidt, G.A., Jungclaus, J.H., Ammann, C.M., Bard, E., Braconnot, P., Crowley, T.J., Delaygue, G., Joos, F., Krivova, N.A., Muscheler, R., Otto-Bliesner, B.L., Pongratz, J., Shindell, D.T., Solanki, S.K., Steinhilber, F. and Vieira, L.E.A., 2012, "Climate forcing reconstructions for use in PMIP simulations of the Last Millennium (v1.1)", Geosci. Model Dev., 5, 185-191. [DOI], [ADS]. URL (accessed 11 September 2012): (Cited on pages 43 and 59.)

Schmitt, D., Schussler, M. and Ferriz-Mas, A., 1996, "Intermittent solar activity by an on-off dynamo", Astron. Astrophys., 311, L1-L4. [ADS] (Cited on pages 18, 52, and 58.)

Schove, D.J., 1955, "The Sunspot Cycle, 649 B.C. to A.D. 2000", J. Geophys. Res., 60, 127-146. [DOI], [ADS] (Cited on pages 15 and 16.)

Schove, D.J., 1979, "Sunspot Turning-Points and Aurorae Since A.D. 1510", Solar Phys., 63, 423-432. [DOI], [ADS] (Cited on page 15.)

Schove, D.J., 1983, "Sunspot, auroral, radiocarbon and climatic fluctuations since 7000 BC", Ann. Geophys., 1, 391-396. [ADS] (Cited on page 15.)

Schrijver, C.J., Beer, J., Baltensperger, U., Cliver, E.W., Gudel, M, Hudson, H.S., McCracken, K.G., Osten, R.A., Peter, T., Soderblom, D.R., Usoskin, I.G. and Wolff, E.W., 2012, "Estimating the frequency of extremely energetic solar events, based on solar, stellar, lunar, and terrestrial records", J. Geophys. Res., 117, A08103. [DOI], [ADS] (Cited on pages 62 and 67.)

Schröder, W., 1992, "On the existence of the 11-year cycle in solar and auroral activity before and after the so-called Maunder minimum", J. Geomag. Geoelectr., 44, 119-128 (Cited on page 15.)

Schupbach, S., Federer, U., Bigler, M., Fischer, H. and Stocker, T.F., 2011, "A refined TALDICE-1a age scale from 55 to 112 ka before present for the Talos Dome ice core based on high-resolution methane measurements", Glim. Past, 7, 1001-1009. [DOI], [ADS]. URL (accessed 11 September 2012): (Cited on page 38.)

Schussler, M., Caligari, P., Ferriz-Mas, A. and Moreno-Insertis, F., 1994, "Instability and eruption of magnetic flux tubes in the solar convection zone", Astron. Astrophys., 281, L69-L72. [ADS] (Cited on page 52.)

Schussler, M., Schmitt, D. and Ferriz-Mas, A., 1997, "Long-term Variation of Solar Activity by a Dynamo Based on Magnetic Flux Tubes", in 1st Advances in Solar Physics Euroconference: Advances in the Physics of Sunspots, (Eds.) Schmieder, B., del Toro Iniesta, J.C., Vazquez, M., ASP Conference Series, 118, pp. 39-44, Astronomical Society of the Pacific, San Francisco. [ADS] (Cited on page 16.)

Scuderi, L.A., 1990, "Oriental sunspot observations and volcanism", Quart. J. R. Astron. Soc., 31, 109120. [ADS] (Cited on page 14.)

Sello, S., 2000, "Wavelet entropy as a measure of solar cycle complexity", Astron. Astrophys., 363, 311-315. [ADS], [arXiv:astro-ph/0005334] (Cited on pages 17 and 18.)

Semeniuk, K., Fomichev, V.I., McConnell, J.C., Fu, C., Melo, S.M.L. and Usoskin, I.G., 2011, "Middle atmosphere response to the solar cycle in irradiance and ionizing particle precipitation", Atmos. Chem. Phys., 11, 5045-5077. [DOI], [ADS]. URL (accessed 11 September 2012): (Cited on page 38.)

Serre, T. and Nesme-Ribes, E., 2000, "Nonlinear analysis of solar cycles", Astron. Astrophys., 360, 319330. [ADS] (Cited on page 18.)

Shapiro, A.I., Schmutz, W., Rozanov, E., Schoell, M., Haberreiter, M., Shapiro, A.V. and Nyeki, S., 2011, "A new approach to the long-term reconstruction of the solar irradiance leads to large historical solar forcing", Astron. Astrophys., 529, A67. [DOI], [ADS], [arXiv:1102.4763 [astro-ph.SR]] (Cited on page 59.)

Shea, M.A. and Smart, D.F., 1990, "A summary of major solar proton events", Solar Phys., 127, 297-320. [DOI], [ADS] (Cited on pages 61 and 63.)

Shea, M.A. and Smart, D.F., 2002, "Solar proton event patterns: the rising portion of five solar cycles", Adv. Space Res., 29, 325-330. [DOI], [ADS] (Cited on page 61.)

Shea, M.A., Smart, D.F., McCracken, K.G., Dreschhoff, G.A.M. and Spence, H.E., 2006, "Solar proton events for 450 years: The Carrington event in perspective", Adv. Space Res., 38, 232-238. [DOI], [ADS] (Cited on page 62.)

Shindell, D., Rind, D., Balachandran, N., Lean, J.L. and Lonergan, P., 1999, "Solar Cycle Variability, Ozone, and Climate", Science, 284, 305-308. [DOI], [ADS] (Cited on page 59.)

Siegenthaler, U., Heimann, M. and Oeschger, H., 1980, "14C variations caused by changes in the global carbon cycle", Radiocarbon, 22, 177-191 (Cited on pages 31, 32, and 33.)

Silverman, S.M., 1992, "Secular variation of the aurora for the past 500 years", Rev. Geophys, 30, 333-351. [DOI] (Cited on pages 15 and 52.)

Silverman, S.M., 2006, "Comparison of the aurora of September 1/2, 1859 with other great auroras", Adv. Space Res., 38, 136-144. [DOI] (Cited on page 13.)

Siscoe, G.L., 1980, "Evidence in the auroral record for secular solar variability", Rev. Geophys. Space Phys., 18, 647-658. [DOI], [ADS] (Cited on page 15.)

Siscoe, G.L. and Verosub, K.L., 1983, "High medieval auroral incidence over China and Japan: implications for the medieval site of the geomagnetic pole", Geophys. Res. Lett., 10, 345-348. [DOI], [ADS] (Cited on page 15.)

Smart, D.F. and Shea, M.A., 2002, "A review of solar proton events during the 22nd solar cycle", Adv. Space Res., 30, 1033-1044. [DOI], [ADS] (Cited on page 62.)

Snow, M., McClintock, W.E., Woods, T.N., White, O.R., Harder, J.W. and Rottman, G., 2005, "The Mg II Index from SORCE", Solar Phys., 230, 325-344. [DOI], [ADS] (Cited on page 13.)

Snowball, I. and Muscheler, R., 2007, "Palaeomagnetic intensity data: an Achilles heel of solar activity reconstructions", Holocene, 17, 851-859. [DOI] (Cited on pages 24 and 34.)

Sokoloff, D., 2004, "The Maunder Minimum and the Solar Dynamo", Solar Phys., 224, 145-152. [DOI], [ADS] (Cited on pages 16, 52, 53, and 58.)

Sokoloff, D. and Nesme-Ribes, E., 1994, "The Maunder minimum: A mixed-parity dynamo mode?", Astron. Astrophys, 288, 293-298. [ADS] (Cited on page 53.)

Solanki, S.K. and Krivova, N.A., 2004, "Solar Irradiance Variations: From Current Measurements to Long-Term Estimates", Solar Phys., 224, 197-208. [DOI], [ADS] (Cited on page 59.)

Solanki, S.K., Schussler, M. and Fligge, M., 2000, "Evolution of the Sun's large-scale magnetic field since the Maunder minimum", Nature, 408, 445-447. [DOI], [ADS] (Cited on pages 13, 42, and 46.)

Solanki, S.K., Schussler, M. and Fligge, M., 2002, "Secular variation of the Sun's magnetic flux", Astron. Astrophys., 383, 706-712. [DOI], [ADS] (Cited on page 42.)

Solanki, S.K., Usoskin, I.G., Kromer, B., Schussler, M. and Beer, J., 2004, "Unusual activity of the Sun during recent decades compared to the previous 11,000 years", Nature, 431, 1084-1087. [DOI], [ADS] (Cited on pages 19, 33, 34, 38, 41, 43, 44, 45, 53, 56, and 57.)

Solanki, S.K., Usoskin, I.G., Kromer, B., Schussler, M. and Beer, J., 2005, "Climate: How unusual is today's solar activity? Reply", Nature, 436, E4-E5. [DOI], [ADS] (Cited on pages 33, 43, 46, and 56.)

Sonett, C.P., 1983, "The great solar anomaly ca. 1780-1800: An error in compiling the record?", J. Geophys. Res., 88, 3225-3228. [DOI], [ADS] (Cited on pages 16 and 17.)

Sonett, C.P. and Finney, S.A., 1990, "The Spectrum of Radiocarbon", Philos. Trans. R. Soc. London, Ser. A, 330, 413-425. [ADS] (Cited on page 17.)

Soon, W.W.-H. and Yaskell, S.H., 2003, The Ma,und,er Minimum and the Variable Sun-Earth Connection, World Scientific, Singapore; River Edge, NJ. [Google Books] (Cited on page 52.)

Steig, E.J., Polissar, P.J., Stuiver, M., Grootes, P.M. and Finkel, R.C., 1996, "Large amplitude solar modulation cycles of 10Be in Antarctica: Implications for atmospheric mixing processes and interpretation of the ice core record", Geophys. Res. Lett., 23, 523-526. [DOI], [ADS] (Cited on page 36.)

Steinhilber, F., Abreu, J.A. and Beer, J., 2008, "Solar modulation during the Holocene", Astrophys. Space Sci. Trans., 4, 1-6. [ADS] (Cited on pages 38 and 43.)

Steinhilber, F., Beer, J. and Frohlich, C., 2009, "Total solar irradiance during the Holocene", Geophys. Res. Lett., 36, L19704. [DOI], [ADS] (Cited on page 43.)

Steinhilber, F., Abreu, J.A., Beer, J. and McCracken, K.G., 2010, "Interplanetary magnetic field during the past 9300 years inferred from cosmogenic radionuclides", J. Geophys. Res., 115, A01104. [DOI], [ADS] (Cited on pages 42 and 43.)

Steinhilber, F., Abreu, J.A., Beer, J., Brunner, I., Christl, M., Fischer, H., Heikkila, U., Kubik, P.W., Mann, M., McCracken, K.G., Miller, H., Miyahara, H., Oerter, H. and Wilhelms, F., 2012, "9,400 years of cosmic radiation and solar activity from ice cores and tree rings", Proc. Natl. Acad. Sci. USA, 109(16), 5967-5971. [DOI], [ADS] (Cited on pages 17, 31, 44, 48, 51, and 56.)

Stephenson, F.R., Willis, D.M. and Hallinan, T.J., 2004, "Aurorae: The earliest datable observation of the aurora borealis", Astron. Geophys., 45, 6.15-6.17. [DOI], [ADS] (Cited on page 15.)

Stuiver, M., 1961, "Variations in Radiocarbon Concentration and Sunspot Activity", J. Geophys. Res., 66, 273-276. [DOI], [ADS] (Cited on page 16.)

Stuiver, M. and Braziunas, T.F., 1989, "Atmospheric 14C and century-scale solar oscillations", Nature, 338, 405-408. [DOI], [ADS] (Cited on pages 53, 54, and 55.)

Stuiver, M. and Braziunas, T.F., 1993, "Sun, ocean, climate and atmospheric 14CO2: an evaluation of casual and spectral relationships", Holocene, 3(4), 289-305. [DOI] (Cited on page 52.)

Stuiver, M. and Pollach, H., 1977, "Discussion: Reporting of 14C data", Radiocarbon, 19, 355-363 (Cited on page 27.)

Stuiver, M. and Quay, P.D., 1980, "Changes in Atmospheric Carbon-14 Attributed to a Variable Sun", Science, 207, 11-19. [DOI], [ADS] (Cited on pages 21, 32, 34, 39, and 54.)

Stuiver, M., Braziunas, T.F., Becker, B. and Kromer, B., 1991, "Climatic, Solar, Oceanic, and Geomagnetic Influences on Late-Glacial and Holocene Atmospheric 14C/12C Change", Quat. Res., 35, 1-24. [DOI] (Cited on pages 31, 34, and 53.)

Stuiver, M., Reimer, P.J., Bard, E., Burr, G.S., Hughen, K.A., Kromer, B., McCormac, G., van der Plicht, J. and Spurk, M., 1998, "INTCAL98 radiocarbon age calibration, 24,000-0 cal BP", Radiocarbon, 40(3), 1041-1083. URL (accessed 11 September 2012): (Cited on

pages 28 and 59.)

Suess, H.E., 1955, "Radiocarbon concentration in modern wood", Science, 122, 415-417. [DOI] (Cited on page 32.)

Suess, H.E., 1965, "Secular Variations of the Cosmic-Ray-Produced Carbon 14 in the Atmosphere and Their Interpretations", J. Geophys. Res., 70, 5937. [DOI], [ADS] (Cited on page 31.)

Suess, H.E., 1980, "The radiocarbon record in tree rings of the last 8000 years", Radiocarbon, 22, 200-209 (Cited on page 17.)

Svalgaard, L., 2012, "How well do we know the sunspot number?", Proc. IAU, 286, 27-33. [DOI], [ADS] (Cited on pages 10 and 11.)

Svensmark, H., 1998, "Influence of Cosmic Rays on Earth's Climate", Phys. Rev. Lett., 81, 5027-5030. [DOI], [ADS] (Cited on page 59.)

Sweeney, C., Gloor, E., Jacobson, A.R., Key, R.M., McKinley, G., Sarmiento, J.L. and Wanninkhof, R., 2007, "Constraining global air-sea gas exchange for CO2 with recent bomb 14C measurements", Global Biogeochem. Cycles, 21, GB2015. [DOI], [ADS] (Cited on page 33.)

Tans, P.P., de Jong, A.F.M. and Mook, W.G., 1979, "Natural atmospheric 14C variation and the Suess effect", Nature, 280, 826-828. [DOI] (Cited on page 32.)

Tapping, K.F., 1987, "Recent solar radio astronomy at centimeter wavelengths: The temporal variability of the 10.7-cm flux", J. Geophys. Res., 92, 829-838. [DOI], [ADS] (Cited on page 12.)

Tapping, K.F. and Charrois, D.P., 1994, "Limits to the Accuracy of the 10.7 cm Flux", Solar Phys., 150, 305-315. [DOI], [ADS] (Cited on page 12.)

Taricco, C., Bhandari, N., Cane, D., Colombetti, P. and Verma, N., 2006, "Galactic cosmic ray flux decline and periodicities in the interplanetary space during the last 3 centuries revealed by 44Ti in meteorites", J. Geophys. Res., 111(A10), A08102. [DOI], [ADS] (Cited on page 46.)

Temmer, M., Veronig, A. and Hanslmeier, A., 2002, "Hemispheric Sunspot Numbers R„ and Rs: Catalogue andN-S asymmetry analysis", Astron. Astrophys., 390, 707-715. [DOI], [ADS], [arXiv:astro-ph/0208436] (Cited on page 11.)

Tinsley, B.A. and Zhou, L., 2006, "Initial results of a global circuit model with variable stratospheric and tropospheric aerosols", J. Geophys. Res., 111(D10), D16205. [DOI], [ADS] (Cited on page 59.)

Tobias, S., Hughes, D. and Weiss, N., 2006, "Unpredictable Sun leaves researchers in the dark", Nature, 442, 26. [DOI], [ADS] (Cited on pages 18 and 19.)

Tobias, S.M., Weiss, N.O. and Kirk, V., 1995, "Chaotically modulated stellar dynamos", Mon. Not. R. Astron. Soc., 273, 1150-1166. [ADS] (Cited on pages 54 and 58.)

Toptygin, I.N., 1985, Cosmic Rays in Interplanetary Magnetic Fields, Geophysics and Astrophysics Monographs, Kluwer Academic Publishers, Dordrecht (Cited on page 21.)

Traversi, R., Usoskin, I.G., Solanki, S.K., Becagli, S., Frezzotti, M., Severi, M., Stenni, B. and Udisti, R., 2012, "Nitrate in Polar Ice: A New Tracer of Solar Variability", Solar Phys., 280, 237-254. [DOI], [ADS] (Cited on page 38.)

Tylka, A. and Dietrich, W., 2009, "A new and comprehensive analysis of proton spectra in ground-level enhanced (GLE) solar particle events", in 31th International Cosmic Ray Conference, 31st ICRC, held on July 7-15, 2009, LodZ, Poland, ICRC0273, Universal Academy Press, LodZ (Cited on page 61.)

Tylka, A.J., Dietrich, W.F. and Boberg, P.R., 1997, "Probability distribution of high-energy solar-heavy-ion fluxes from IMP-8: 1973-1996", IEEE Trans. Nucl. Sci., 44, 2140-2149. [DOI], [ADS] (Cited on page 61.)

Usoskin, I.G., 2011, "Cosmic rays and climate forcing", Mem. Soc. Astron. Ital., 82, 937-942. [ADS] (Cited on page 26.)

Usoskin, I.G. and Kovaltsov, G.A., 2004, "Long-Term Solar Activity: Direct and Indirect Study", Solar Phys., 224, 37-47. [DOI], [ADS] (Cited on pages 39 and 51.)

Usoskin, I.G. and Kovaltsov, G.A., 2006, "Cosmic ray induced ionization in the atmosphere: Full modeling and practical applications", J. Geophys. Res., 111, D21206. [DOI], [ADS] (Cited on page 59.)

Usoskin, I.G. and Kovaltsov, G.A., 2008a, "Production of cosmogenic 7Be isotope in the atmosphere: Full 3D modelling", J. Geophys. Res., 113, D12107. [DOI] (Cited on pages 25, 35, and 36.)

Usoskin, I.G. and Kovaltsov, G.A., 2008b, "Cosmic rays and climate of the Earth: Possible connection", C. R. Geosci., 340, 441-450. [DOI], [ADS] (Cited on page 59.)

Usoskin, I.G. and Kovaltsov, G.A., 2012, "Occurrence of extreme solar particle events: Assessment from historical proxy data", Astrophys. J., 757, 92. [DOI], [ADS], [arXiv:1207.5932 [astro-ph.SR]] (Cited on pages 62, 63, 64, and 66.)

Usoskin, I.G. and Kromer, B., 2005, "Reconstruction of the 14C production rate from measured relative abundance", Radiocarbon, 47, 31-37 (Cited on pages 28 and 32.)

Usoskin, I.G. and Mursula, K., 2003, "Long-Term Solar Cycle Evolution: Review of Recent Developments", Solar Phys., 218, 319-343. [DOI], [ADS] (Cited on pages 19 and 52.)

Usoskin, I.G., Kananen, H., Mursula, K., Tanskanen, P. and Kovaltsov, G.A., 1998, "Correlative study of solar activity and cosmic ray intensity", J. Geophys. Res., 103(A5), 9567-9574. [DOI], [ADS] (Cited on page 22.)

Usoskin, I.G., Mursula, K. and Kovaltsov, G.A., 2000, "Cyclic behaviour of sunspot activity during the Maunder minimum", Astron. Astrophys., 354, L33-L36. [ADS] (Cited on pages 52 and 53.)

Usoskin, I.G., Mursula, K. and Kovaltsov, G.A., 2001a, "Heliospheric modulation of cosmic rays and solar activity during the Maunder minimum", J. Geophys. Res., 106, 16,039-16,046. [DOI], [ADS] (Cited on pages 52 and 53.)

Usoskin, I.G., Mursula, K. and Kovaltsov, G.A., 2001b, "Was one sunspot cycle lost in late XVIII century?", Astron. Astrophys., 370, L31-L34. [DOI], [ADS] (Cited on page 17.)

Usoskin, I.G., Mursula, K. and Kovaltsov, G.A., 2002a, "Lost sunspot cycle in the beginning of Dalton minimum: New evidence and consequences", Geophys. Res. Lett., 29, 36-1. [DOI], [ADS] (Cited on page 17.)

Usoskin, I.G., Mursula, K., Solanki, S.K., Schüssler, M. and Kovaltsov, G.A., 2002b, "A physical reconstruction of cosmic ray intensity since 1610", J. Geophys. Res., 107, 1374. [DOI], [ADS] (Cited on pages 23, 40, 41, and 42.)

Usoskin, I.G., Mursula, K. and Kovaltsov, G.A., 2003a, "Reconstruction of monthly and yearly group sunspot numbers from sparse daily observations", Solar Phys., 218, 295-305. [DOI] (Cited on page 11.)

Usoskin, I.G., Mursula, K. and Kovaltsov, G.A., 2003b, "The lost sunspot cycle: Reanalysis of sunspot statistics", Astron. Astrophys., 403, 743-748. [DOI], [ADS] (Cited on page 17.)

Usoskin, I.G., Solanki, S.K., Schussler, M., Mursula, K. and Alanko, K., 2003c, "Millennium-Scale Sunspot Number Reconstruction: Evidence for an Unusually Active Sun since the 1940s", Phys. Rev. Lett., 91, 211101. [DOI], [ADS], [arXiv:astro-ph/0310823] (Cited on pages 19, 36, 38, 41, 43, 44, 45, 47, and 56.)

Usoskin, I.G., Mursula, K., Solanki, S.K., Schussler, M. and Alanko, K., 2004, "Reconstruction of solar activity for the last millennium using 10Be data", Astron. Astrophys., 413, 745-751. [DOI], [ADS], [arXiv:astro-ph/0309556] (Cited on pages 17, 35, 42, 43, and 56.)

Usoskin, I.G., Alanko-Huotari, K., Kovaltsov, G.A. and Mursula, K., 2005, "Heliospheric modulation of cosmic rays: Monthly reconstruction for 1951-2004", J. Geophys. Res., 110(A9), A12108. [DOI], [ADS] (Cited on pages 23, 30, 40, and 41.)

Usoskin, I.G., Solanki, S.K. and Korte, M., 2006a, "Solar activity reconstructed over the last 7000 years: The influence of geomagnetic field changes", Geophys. Res. Lett., 33, 8103. [DOI], [ADS] (Cited on pages 34, 44, and 56.)

Usoskin, I.G., Solanki, S.K., Kovaltsov, G.A., Beer, J. and Kromer, B., 2006b, "Solar proton events in cosmogenic isotope data", Geophys. Res. Lett., 33, L08107. [DOI], [ADS] (Cited on pages 31, 35, and 62.)

Usoskin, I.G., Solanki, S.K., Taricco, C., Bhandari, N. and Kovaltsov, G.A., 2006c, "Long-term solar activity reconstructions: direct test by cosmogenic 44Ti in meteorites", Astron. Astrophys., 457, L25-L28. [DOI], [ADS] (Cited on page 46.)

Usoskin, I.G., Solanki, S.K. and Kovaltsov, G.A., 2007, "Grand minima and maxima of solar activity: new observational constraints", Astron. Astrophys., 471, 301-309. [DOI], [ADS], [arXiv:0706.0385] (Cited on pages 38, 42, 43, 44, 50, 51, 52, 54, 55, 57, and 58.)

Usoskin, I.G., Korte, M. and Kovaltsov, G.A., 2008, "Role of centennial geomagnetic changes in local atmospheric ionization", Geophys. Res. Lett., 35, L05811. [DOI], [ADS] (Cited on pages 24 and 59.)

Usoskin, I.G., Field, C.V., Schmidt, G.A., Leppanen, A.-P., Aldahan, A., Kovaltsov, G.A., Possnert, G. and Ungar, R.K., 2009a, "Short-term production and synoptic influences on atmospheric 7Be concentrations", J. Geophys. Res., 114(D13), D06108. [DOI], [ADS] (Cited on pages 35 and 36.)

Usoskin, I.G., Horiuchi, K., Solanki, S.K., Kovaltsov, G.A. and Bard, E., 2009b, "On the common solar signal in different cosmogenic isotope data sets", J. Geophys. Res., 114(A13), A03112. [DOI], [ADS] (Cited on page 47.)

Usoskin, I.G., Mursula, K., Arlt, R. and Kovaltsov, G.A., 2009c, "A Solar Cycle Lost in 1793-1800: Early Sunspot Observations Resolve the Old Mystery", Astrophys. J. Lett., 700, L154-L157. [DOI], [ADS], [arXiv:0907.0063 [astro-ph.SR]] (Cited on pages 12, 17, and 53.)

Usoskin, I.G., Sokoloff, D. and Moss, D., 2009d, "Grand Minima of Solar Activity and the Mean-Field Dynamo", Solar Phys., 254, 345-355. [DOI], [ADS] (Cited on pages 52 and 55.)

Usoskin, I.G., Mironova, I.A., Korte, M. and Kovaltsov, G.A., 2010, "Regional millennial trend in the cosmic ray induced ionization of the troposphere", J. Atmos. Sol.-Terr. Phys., 72, 19-25. [DOI], [ADS] (Cited on page 24.)

Usoskin, I.G., Bazilevskaya, G.A. and Kovaltsov, G.A., 2011, "Solar modulation parameter for cosmic rays since 1936 reconstructed from ground-based neutron monitors and ionization chambers", J. Geophys. Res., 116, A02104. [DOI], [ADS] (Cited on pages 23 and 41.)

Usoskin, I.G., Kromer, B., Ludlow, F., Beer, J., Friedrich, M., Kovaltsov, G.A., Solanki, S.K. and Wacker, L., 2013, "The AD775 cosmic event revisited: the Sun is to blame", Astron. Astrophys., 552, L3. [DOI], [ADS], [arXiv:1302.6897 [astro-ph.SR]] (Cited on pages 62, 63, and 64.)

Vaquero, J.M., 2007, "Historical sunspot observations: A review", Adv. Space Res., 40, 929-941. [DOI], [ADS], [arXiv:astro-ph/0702068] (Cited on page 14.)

Vaquero, J.M. and Vazquez, M., 2009, The Sun Recorded Through History: Scientific Data Extracted from Historical Documents, Astrophysics and Space Science Library, 361, Springer, Berlin; New York. [Google Books] (Cited on pages 11 and 14.)

Vaquero, J.M., Gallego, M.C. and Garcia, J.A., 2002, "A 250-year cycle in naked-eye observations of sunspots", Geophys. Res. Lett., 29, 58-1. [DOI], [ADS] (Cited on page 15.)

Vaquero, J.M., Gallego, M.C. and Sanchez-Bajo, F., 2004, "Reconstruction of a Monthly Homogeneous Sunspot Area Series Since 1832", Solar Phys., 221, 179-189. [DOI], [ADS] (Cited on page 13.)

Vaquero, J.M., Trigo, R.M. and Gallego, M.C., 2005, "A 'lost' sunspot observation in 1785", Astron. Nachr., 326, 112-114. [DOI], [ADS] (Cited on page 11.)

Vaquero, J.M., Trigo, R.M., Gallego, M.C. and Moreno-Corral, M.A., 2007, "Two Early Sunspots Observers: Teodoro de Almeida and Jose Antonio Alzate", Solar Phys., 240, 165-175. [DOI], [ADS] (Cited on page 11.)

Vaquero, J.M., Gallego, M.C., Usoskin, I.G. and Kovaltsov, G.A., 2011, "Revisited Sunspot Data: A New Scenario for the Onset of the Maunder Minimum", Astrophys. J. Lett., 731, L24. [DOI], [ADS], [arXiv:1103.1520 [astro-ph.SR]] (Cited on pages 12 and 53.)

Vaquero, J.M., Trigo, R.M. and Gallego, M.C., 2012, "A Simple Method to Check the Reliability of Annual Sunspot Number in the Historical Period 1610-1847", Solar Phys., 277, 389-395. [DOI], [ADS], [arXiv:1111.2633 [astro-ph.SR]] (Cited on page 11.)

Vasiliev, S.S. and Dergachev, V.A., 2002, "The ^2400-year cycle in atmospheric radiocarbon concentration: bispectrum of 14C data over the last 8000 years", Ann. Geophys., 20, 115-120. [ADS]. URL (accessed

14 October 2008): (Cited on pages 17 and 51.)

Vazquez, M. and Vaquero, J.M., 2010, "Aurorae Observed at the Canary Islands", Solar Phys., 267, 431-444. [DOI], [ADS] (Cited on page 15.)

Versteegh, G.J.M., 2005, "Solar Forcing of Climate. 2: Evidence from the Past", Space Sci. Rev., 120, 243-286. [DOI], [ADS] (Cited on page 59.)

Vieira, L.E.A., Solanki, S.K., Krivova, N.A. and Usoskin, I.G., 2011, "Evolution of the solar irradiance during the Holocene", Astron. Astrophys., 531, A6. [DOI], [ADS], [arXiv:1103.4958 [astro-ph.SR]] (Cited on page 43.)

Viereck, R.A. and Puga, L.C., 1999, "The NOAA Mg II core-to-wing solar index: Construction of a 20-year time series of chromospheric variability from multiple satellites", J. Geophys. Res., 104, 9995-10 006. [DOI], [ADS] (Cited on page 13.)

Vitinsky, Y.I., 1965, Solar Activity Forecasting, Israel Program for Scientific Translations, Jerusalem (Cited on pages 16, 17, and 18.)

Vitinsky, Y.I., Kopecky, M. and Kuklin, G.V., 1986, Statistics of Sunspot Activity (in Russian), Nauka, Moscow (Cited on pages 10, 11, 16, 17, 51, 52, and 53.)

Vogt, S., Herzog, G.F. and Reedy, R.C., 1990, "Cosmogenic nuclides in extraterrestrial materials", Rev. Geophys., 28, 253-275. [DOI], [ADS] (Cited on page 45.)

Vonmoos, M.V., Beer, J. and Muscheler, R., 2006, "Large variations in Holocene solar activity: Constraints from 10Be in the Greenland Ice Core Project ice core", J. Geophys. Res., 111 (A10), A10105. [DOI], [ADS] (Cited on pages 34, 36, 38, 41, 43, and 47.)

Voss, H., Kurths, J. and Schwarz, U., 1996, "Reconstruction of grand minima of solar activity from A14C data: Linear and nonlinear signal analysis", J. Geophys. Res., 101, 15,637-15,644. [DOI], [ADS] (Cited on pages 34 and 53.)

Waldmeier, M., 1961, The Sunspot-Activity in the Years 1610-1960, Schulthess u. Co. / Swiss Federal Observatory, Zurich. [ADS] (Cited on page 10.)

Wang, Y.-M. and Sheeley Jr, N.R., 2013, "The Solar Wind and Interplanetary Field during Very Low Amplitude Sunspot Cycles", Astrophys. J., 764, 90. [DOI] (Cited on page 52.)

Wang, Y.-M., Lean, J.L. and Sheeley Jr, N.R., 2005, "Modeling the Sun's Magnetic Field and Irradiance since 1713", Astrophys. J., 625, 522-538. [DOI], [ADS] (Cited on pages 13, 46, and 59.)

Webber, W.R. and Higbie, P.R., 2003, "Production of cosmogenic Be nuclei in the Earth's atmosphere by cosmic rays: Its dependence on solar modulation and the interstellar cosmic ray spectrum", J. Geophys. Res., 108, 1355. [DOI], [ADS] (Cited on pages 23, 25, 29, 35, and 44.)

Webber, W.R. and Higbie, P.R., 2009, "Galactic propagation of cosmic ray nuclei in a model with an increasing diffusion coefficient at low rigidities: A comparison of the new interstellar spectra with Voyager data in the outer heliosphere", J. Geophys. Res., 114(A13), A02103. [DOI], [ADS] (Cited on page 23.)

Webber, W.R., Higbie, P.R. and McCracken, K.G., 2007, "Production of the cosmogenic isotopes 3H, 7Be,

10Be, and 36Cl in the Earth's atmosphere by solar and galactic cosmic rays", J. Geophys. Res., 112, A10106. [DOI], [ADS] (Cited on pages 25, 35, and 62.)

Weiss, N.O. and Tobias, S.M., 2000, "Physical Causes of Solar Activity", Space Sci. Rev., 94, 99-112. [DOI], [ADS] (Cited on pages 17, 52, and 58.)

Weiss, N.O., Cattaneo, F. and Jones, C.A., 1984, "Periodic and aperiodic dynamo waves", Geophys. Astrophys. Fluid Dyn, 30, 305-341. [DOI], [ADS] (Cited on page 18.)

Wheatland, M.S., 2003, "The Coronal Mass Ejection Waiting-Time Distribution", Solar Phys., 214, 361373. [DOI], [ADS], [arXiv:astro-ph/0303019] (Cited on page 58.)

Willis, D.M. and Stephenson, F.R., 2001, "Solar and auroral evidence for an intense recurrent geomagnetic storm during December in AD 1128", Ann. Geophys., 19, 289-302. [ADS]. URL (accessed 14 October 2008): (Cited on page 14.)

Willis, D.M., Easterbrook, M.G. and Stephenson, F.R., 1980, "Seasonal variation of oriental sunspot sightings", Nature, 287, 617-619. [DOI], [ADS] (Cited on page 14.)

Willis, D.M., Davda, V.N. and Stephenson, F.R., 1996, "Comparison between Oriental and Occidental Sunspot Observations", Quart. J. R. Astron. Soc., 37, 189-229. [ADS] (Cited on page 15.)

Wilson, P.R., 1994, Solar and Stellar Activity Cycles, Cambridge Astrophysics Series, 24, Cambridge University Press, Cambridge; New York (Cited on page 17.)

Wilson, R.M., 1988a, "Bimodality and the Hale cycle", Solar Phys., 117, 269-278. [DOI], [ADS] (Cited on page 16.)

Wilson, R.M., 1988b, "On the long-term secular increase in sunspot number", Solar Phys., 115, 397-408. [DOI], [ADS] (Cited on page 56.)

Wilson, R.M., 1998, "A Comparison of Wolf's Reconstructed Record of Annual Sunspot Number with Schwabe's Observed Record of 'Clusters of Spots' for the Interval of 1826-1868", Solar Phys., 182, 217-230. [DOI], [ADS] (Cited on page 11.)

Wittmann, A.D. and Xu, Z.T., 1987, "A catalogue of sunspot observations from 165 BC to AD 1684", Astron. Astrophys. Suppl., 70, 83-94. [ADS] (Cited on page 14.)

Wolff, E.W., Bigler, M., Curran, M.A.J., Dibb, J.E., Frey, M.M., Legrand, M. and McConnell, J.R., 2012, "The Carrington event not observed in most ice core nitrate records", Geophys. Res. Lett., 39, L08503. [DOI], [ADS] (Cited on pages 64 and 66.)

Yang, S., Odah, H. and Shaw, J., 2000, "Variations in the geomagnetic dipole moment over the last 12000 years", Geophys. J. Int., 140, 158-162. [DOI], [ADS] (Cited on pages 25, 43, and 50.)

Yau, K.K.C. and Stephenson, F.R., 1988, "A Revised Catalog of Far Eastern Observations of Sunspots (165 BC to AD 1918)", Quart. J. R. Astron. Soc., 29, 175-197. [ADS] (Cited on page 14.)

Yiou, F., Raisbeck, G.M., Baumgartner, S., Beer, J., Hammer, C., Johnsen, S., Jouzel, J., Kubik, P.W., Lestringuez, J., Stievenard, M., Suter, M. and Yiou, P., 1997, "Beryllium 10 in the greenland ice core project ice core at Summit, Greenland", J. Geophys. Res., 102, 26,783-26,794. [DOI], [ADS] (Cited on pages 34 and 63.)

Yule, G.U., 1927, "On a Method of Investigating Periodicities in Disturbed Series, with Special Reference to Wolfer's Sunspot Numbers", Philos. Trans. R. Soc. London, Ser. A, 226, 267-298. Online version (accessed 25 February 2008): (Cited on page 18.)

Zeller, E.J. and Dreschhoff, G.A.M., 1995, "Anomalous Nitrate Concentrations in Polar Ice Cores - do they Result from Solar Particle Injections into the Polar Atmosphere?", Geophys. Res. Lett., 22, 2521-2524. [DOI], [ADS] (Cited on page 66.)

Zhentao, X., 1990, "Solar Observations in Ancient China and Solar Variability", Philos. Trans. R. Soc. London, Ser. A, 330, 513-516. [ADS] (Cited on page 17.)

Zolotova, N.V. and Ponyavin, D.I., 2007, "Was the unusual solar cycle at the end of the XVIII century a result of phase asynchronization?", Astron. Astrophys., 470, L17-L20. [DOI], [ADS] (Cited on page 17.)

Copyright of Living Reviews in Solar Physics is the property of Max-Planck-Institut fur Sonnensystemforschung and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use.