Our direct knowledge of solar eruptive events is limited to several decades and does not include extreme events, which can only be studied by the indirect proxy method over millennia, or by a large number of sun-like stars. There is a gap, spanning 1--2 orders of magnitude, in the strength of events between directly observed and reconstructed ones. Here, we study the proxy-method sensitivity to identify extreme solar particle events (SPEs). First, the strongest directly observed SPE (23-Feb-1956), used as a reference for proxy-based reconstructions, was revisited using the newly developed method. Next, the sensitivity of the cosmogenic-isotope method to detect a reference SPE was assessed against the precision and number of individual isotopic records, showing that it is too weak by a factor $\approx$30 to be reliably identified in a single record. Uncertainties of 10Be and 14C data are shown to be dominated by local/regional patterns and measurement errors, respectively. By combining several proxy records, an SPE 4--5 times stronger than the reference one can be potentially detected, increasing the present-day sensitivity by an order of magnitude. This will allow filling the observational gap in SPE strength distribution, thus enriching statistics of extreme events from 3--4 presently known ones to several tens. This will provide a solid basis for research in the field of extreme events, both for fundamental science, viz. solar and stellar physics, and practical applications, such as the risk assessments of severe space-based hazards for the modern technological society.

Numerous analyses suggest the existence of various quasi-periodicities in solar activity. The power spectrum of solar activity recorded in sunspot data is dominated by the $\sim$11-year quasi-periodicity, known as the Schwabe cycle. In the mid-term range (1 month -- 11 years) a pronounced variability known as a quasi-biennial oscillation (QBO) is widely discussed. In the shorter time scale a pronounced peak, corresponding to the synodic solar rotation period ($\sim$ 27 days) is observed. Here we revisited the mid-term solar variability in terms of statistical dynamic of fully turbulent systems, where solid arguments are required to accept an isolated dominant frequency in a continuous (smooth) spectrum. For that, we first undertook an unbiased analysis of the standard solar data, sunspot numbers and the F10.7 solar radioflux index, by applying a wavelet tool, which allows one to perform a frequency-time analysis of the signal. Considering the spectral dynamics of solar activity cycle by cycle, we showed that no single periodicity can be separated, in a statistically significant manner, in the specified range of periods. We examine whether a model of solar dynamo can reproduce the mid-term oscillation pattern observed in solar data. We found that a realistically observed spectrum can be explained if small spatial (but not temporal) scales are effectively smoothed. This result is important because solar activity is a it global feature, although monitored via small-scale tracers like sunspots.

The newly published spectra of protons and helium over time directly measured in space by the AMS-02 experiment for the period 2011--2017 provide a unique opportunity to calibrate ground-based neutron monitors (NMs). Here, calibration of several stable sealevel NMs (Inuvik, Apatity, Oulu, Newark, Moscow, Hermanus, Athens) was performed using these spectra. Four modern NM yield functions were verified: Mi13 (Mishev et al., 2013), Ma16 (Mangeard et al., 2016), CM12 (Caballero-Lopez & Moraal, 2012) and CD00 (Clem & Dorman, 2000), on the basis of the cosmic-ray spectra measured by AMS-02. The Mi13 yield function was found to realistically represent the NM response to galactic cosmic rays. CM12 yield function leads to a small skew in the solar cycle dependence of the scaling factor. In contrast, Ma16 and CD00 yield functions tend to overestimate the NM sensitivity to low-rigidity (<10 GV) cosmic rays. This effect may be important for an analysis of ground level enhancements, leading to a potential underestimate of fluxes of solar energetic particles as based on NM data. The Mi13 yield function is recommended for quantitative analyses of NM data, especially for ground-level enhancements. The validity the force-field approximation was studied, and it was found that it fits well the directly measured proton spectra, within a few % for periods of low to moderate activity and up to ~10% for active periods. The results of this work strengthen and validate the method of the cosmic-ray variability analysis based on the NM data and yield-function formalism, and improves its accuracy.

Over the past decades and even centuries, the astronomical community has accumulated a signif-icant heritage of recorded observations of a great many astronomical objects. Those records con-tain irreplaceable information about long-term evolutionary and non-evolutionary changes in our Universe, and their preservation and digitization is vital. Unfortunately, most of those data risk becoming degraded and thence totally lost. We hereby call upon the astronomical community and US funding agencies to recognize the gravity of the situation, and to commit to an interna-tional preservation and digitization efforts through comprehensive long-term planning supported by adequate resources, prioritizing where the expected scientific gains, vulnerability of the origi-nals and availability of relevant infrastructure so dictates. The importance and urgency of this issue has been recognized recently by General Assembly XXX of the International Astronomical Union (IAU) in its Resolution B3: "on preservation, digitization and scientific exploration of his-torical astronomical data". We outline the rationale of this promotion, provide examples of new science through successful recovery efforts, and review the potential losses to science if nothing it done.

Changes in solar irradiance and in its spectral distribution are among the main natural drivers of the climate on Earth. However, irradiance measurements are only available for less than four decades, while assessment of solar influence on Earth requires much longer records. The aim of this work is to provide the most up-to-date physics-based reconstruction of the solar total and spectral irradiance (TSI/SSI) over the last nine millennia. The concentrations of the cosmogenic isotopes 14C and 10Be in natural archives have been converted to decadally averaged sunspot numbers through a chain of physics-based models. TSI and SSI are reconstructed with an updated SATIRE model. Reconstructions are carried out for each isotope record separately, as well as for their composite. We present the first ever SSI reconstruction over the last 9000 years from the individual 14C and 10Be records as well as from their newest composite. The reconstruction employs physics-based models to describe the involved processes at each step of the procedure. Irradiance reconstructions based on two different cosmogenic isotope records, those of 14C and 10Be, agree well with each other in their long-term trends despite their different geochemical paths in the atmosphere of Earth. Over the last 9000 years, the reconstructed secular variability in TSI is of the order of 0.11%, or 1.5 W/m2. After the Maunder minimum, the reconstruction from the cosmogenic isotopes is consistent with that from the direct sunspot number observation. Furthermore, over the nineteenth century, the agreement of irradiance reconstructions using isotope records with the reconstruction from the sunspot number by Chatzistergos et al. (2017) is better than that with the reconstruction from the WDC-SILSO series (Clette et al. 2014), with a lower chi-square-value.

Using neutron monitor and space-borne data we performed an analysis of the second ground level enhancement of solar cycle 24, namely the event of 10 September 2017 (GLE 72) and derive the spectral and angular characteristics of GLE particles. We employ new neutron monitor yield function and a recently proposed model based on optimization procedure. The method consists of simulation of particle propagation in a model magnetosphere in order to derive the cut-off rigidity and neutron monitor asymptotic directions. Subsequently the rigidity spectrum and anisotropy of GLE particles are obtained in their dynamical evolution during the event on the basis of inverse problem solution. The derived angular distribution and spectra are briefly discussed.

Aims. Lunar soil and rocks are not protected by a magnetic field or an atmosphere and are continuously irradiated by energetic particles that can produce cosmogenic radioisotopes directly inside rocks at different depths depending on the particle's energy. This allows the mean fluxes of solar and galactic cosmic rays to be assessed on the very long timescales of millions of years. Methods. Here we show that lunar rocks can serve as a very good particle integral spectrometer in the energy range 20-80 MeV. We have developed a new method based on precise modeling, that is applied to measurements of $^{26}$Al (half-life ~0.7 megayears) in lunar samples from the Apollo mission, and present the first direct reconstruction (i.e., without any a priori assumptions) of the mean energy spectrum of solar and galactic energetic particles over a million of years. Results. We show that the reconstructed spectrum of solar energetic particles is totally consistent with that over the last decades, despite the very different levels of solar modulation of galactic cosmic rays ($\phi=496\pm 40$ MV over a million years versus $\phi= 660\pm 20$ MV for the modern epoch). We also estimated the occurrence probability of extreme solar events and argue that no events with the F(>30 MeV) fluence exceeding $5*10^{10}$ and $10^{11}$ cm$^2$ are expected on timescales of a thousand and million years, respectively. Conclusions. We conclude that the mean flux of solar energetic particles hardly depends on the level of solar activity, in contrast to the solar modulation of galactic cosmic rays. This puts new observational constraints on solar physics and becomes important for assessing radiation hazards for the planned space missions.

Solar activity in the past millennia can only be reconstructed from cosmogenic radionuclide records in terrestrial archives. However, because of the diversity of the proxy archives, it is difficult to build a homogeneous reconstruction. Here we provide a new consistent multiproxy reconstruction of the solar activity over the last 9000 years, using available long-span datasets of 10Be and 14C in terrestrial archives. A new method, based on a Bayesian approach, was applied for the first time to solar activity reconstruction. A Monte Carlo search for the most probable value of the modulation potential was performed to match data from different datasets for a given time. We used six 10Be series from Greenland and Antarctica, and the global 14C production series. The 10Be series were resampled to match wiggles related to the grand minima in the 14C reference dataset. The GRIP and the EDML 10Be series diverge from each other during the second half of the Holocene, while the 14C series lies between them. A likely reason for this is the insufficiently precise beryllium transport and deposition model for Greenland. A slow 6-millennia variability with lows at ca. 5500 BC and 1500 AD of solar activity is found. Two components of solar activity can be statistically distinguished: the main 'normal' component and a component corresponding to grand minima. A possible existence of a component representing grand maxima is indicated, but it cannot be separated from the main component in a statistically significant manner. A new consistent reconstruction of solar activity over the last nine millennia is presented with the most probable values of decadal sunspot numbers and their realistic uncertainties. Independent components of solar activity corresponding to the main moderate activity and the grand-minimum state are identified; they may be related to different operation modes of the dynamo.

Context. There is no consensus on the amplitude of the historical solar forcing. The estimated magnitude of the total solar irradiance difference between Maunder minimum and present time ranges from 0.1 to 6 W/m2 making uncertain the simulation of the past and future climate. One reason for this disagreement is the applied evolution of the quiet Sun brightness in the solar irradiance reconstruction models. This work addresses the role of the quiet Sun model choice and updated solar magnetic activity proxies on the solar forcing reconstruction. Aims. We aim to establish a plausible range of the solar irradiance variability on decadal to millennial time scales. Methods. The spectral solar irradiance (SSI) is calculated as a weighted sum of the contributions from sunspot umbra/penumbra, fac- ulae and quiet Sun, which are pre-calculated with the spectral synthesis code NESSY. We introduce activity belts of the contributions from sunspots and faculae and a new structure model for the quietest state of the Sun. We assume that the brightness of the quiet Sun varies in time proportionally to the secular (22-year smoothed) variation of the solar modulation potential. Results. A new reconstruction of the TSI and SSI covering the period 6000 BCE - 2015 CE is presented. The model simulates solar irradiance variability during the satellite era well. The TSI change between the Maunder and recent minima ranges between 3.7 and 4.5 W/m2 depending on the applied solar modulation potential. The implementation of a new quietest Sun model reduces, by approximately a factor of two, the relative solar forcing compared to the largest previous estimation, while the application of updated solar modulation potential increases the forcing difference between Maunder minimum and the present by 25-40 %.

The method of active day fraction (ADF) was proposed recently to calibrate different solar observers to the standard observational conditions. The result of the calibration may depend on the overall level of solar activity during the observational period. This dependency is studied quantitatively using data of the Royal Greenwich Observatory, by formally calibrating synthetic pseudo-observers to the full reference dataset. It is shown that the sunspot group number is precisely estimated by the ADF method for periods of moderate activity, may be slightly underestimated by 0.5\,--\,1.5 groups ($\leq$10\,\%) for strong and very strong activity, and is strongly overestimated by up to 2.5 groups ($\leq$30\,\%) for weak--moderate activity. The ADF method becomes unapplicable for the periods of grand minima of activity. In general, the ADF method tends to overestimate the overall level of activity and to reduce the long-term trends.

In this work, we use two decadal sunspot number series reconstructed from cosmogenic radionuclide data (14C in tree trunks, SN-14C and 10Be in polar ice, SN-10Be) and the Extreme Value Theory to study variability of solar activity during the last 9 millennia. The peaks-over-threshold technique was used to compute, in particular, the shape parameter of the generalized Pareto distribution for different thresholds. Its negative value implies an upper bound of the extreme SN-10Be and SN-14C time series. The return level for 1000 and 10000 years were estimated leading to values lower than the maximum observed values, expected for the 1000-year, but not for the 10000-year return levels, for both series. A comparison of these results with those obtained using the observed sunspot numbers from telescopic observations during the last four centuries suggest that the main characteristics of solar activity have already been recorded in the telescopic period (from 1610 to nowadays) which covers the full range of solar variability from a Grand minimum to a Grand maximum.

The conventional definition of ground-level enhancement (GLE) events requires a detection of solar energetic particles (SEP) by at least two differently located neutron monitors. Some places are exceptionally well suitable for ground-based detection of SEP - high-elevation polar regions with negligible geomagnetic and reduced atmospheric energy/rigidity cutoffs. At present, there are two neutron-monitor stations in such locations on the Antarctic plateau: SOPO/SOPB (at Amundsen-Scott station, 2835 m elevation), and DOMC/DOMB (at Concordia station, 3233 m elevation). Since 2015, when the DOMC/DOMB station started continuous operation, a relatively weak SEP event that was not detected by sea-level neutron-monitor stations was registered by both SOPO/SOPB and DOMC/DOMB, and it was accordingly classified as a GLE. This would lead to a distortion of the homogeneity of the historic GLE list and the corresponding statistics. To address this issue, we propose to modify the GLE definition so that it maintains the homogeneity: A GLE event is registered when there are near-time coincident and statistically significant enhancements of the count rates of at least two differently located neutron monitors, including at least one neutron monitor near sea level and a corresponding enhancement in the proton flux measured by a space-borne instrument(s). Relatively weak SEP events registered only by high-altitude polar neutron monitors, but with no response from cosmic-ray stations at sea level, can be classified as sub-GLEs.

The paper by Popova et al. presents an oversimplified mathematical model of solar activity with a claim of predicting/postdicting it for several millennia ahead/backwards. The work contains several flaws devaluating the results: (1) the method is unreliable from the point of view of signal processing (it is impossible to make harmonic predictions for thousands of years based on only 35 years of data) and lacks quality control, (2) the result of post-diction apparently contradicts the observational data. (3) theoretical speculations make little sense, To summarize, a multi-harmonic mathematical model, hardly related to full solar dynamo theory, is presented, which is not applicable to realistic solar conditions because of the significant chaotic/stochastic intrinsic component and strong non-stationarity of solar activity. The obtained result is apparently inconsistent with the data in the past and thus cannot be trusted for the future predictions.

Recent data from the Kepler mission has revealed the occurrence of superflares in sun-like stars which exceed by far any observed solar flares in release of energy. A natural idea is that the dynamo mechanism in superflaring stars differs in some respect from that in the Sun. We search for a difference in the dynamo-related parameters between superflaring stars and the Sun to suggest a dynamo-mechanism as close as possible to the conventional solar/stellar dynamo but capable of providing much higher magnetic energy. Dynamo based on joint action of differential rotation and mirror asymmetric motions can in principle result in excitation of two types of magnetic fields. First of all, it is well-known in solar physics dynamo waves. The point is that another magnetic configuration with initial growth and further stabilisation is also possible for excitation. For comparable conditions, magnetic field strength of second configuration is much larger rather of the first one just because dynamo do not spend its efforts for periodic magnetic field inversions but use its for magnetic field growth. We analysed available data from the Kepler mission concerning the superflaring stars in order to find tracers of anomalous magnetic activity. Starting from the recent paper, we find that anti-solar differential rotation or anti-solar sign of the mirror-asymmetry of stellar convection can provide the desired strong magnetic field in dynamo models. We confirm this concept by numerical models of stellar dynamos with corresponding governing parameters. We conclude that the proposed mechanism can plausibly explain the superflaring events at least for some cool stars, including binaries, subgiants and, possibly, low-mass stars and young rapid rotators.

Sunspot number series are composed from observations of hundreds of different observers that requires careful normalization of the observers to the standard conditions. Here we present a new normalized series of the number of sunspot groups for the period 1749-1996. The reconstruction is based on the active day fraction (ADF) method, which is slightly updated with respect to the previous works, and a revised database of sunspot group observations. Stability of some key solar observers has been evaluated against the composite series. The Royal Greenwich Observatory dataset appears fairly stable since the 1890s but is about 10% too low before that. A declining trend of 10-15% in the quality of Wolfer's observation is found between the 1880s and 1920s, suggesting that using him as the reference observer may lead to additional uncertainties. Wolf (small telescope) appears fairly stable between the 1860s and 1890s, without any obvious trend. The new reconstruction reflects the centennial variability of solar activity as evaluated using the singular spectrum analysis method. It depicts a highly significant feature of the Modern grand maximum of solar activity in the second half of the 20th century, being a factor 1.33-1.77 higher than during the 18-19th centuries. The new series of the sunspot group numbers with monthly and annual resolution, available also in the electronic format, is provided forming a basis for new studies of the solar variability and solar dynamo for the last 250 years.

Group sunspot number (GSN) series constitute the longest instrumental astronomical database providing information on solar activity. It is a compilation of observations by many individual observers, and their inter-calibration has usually been performed using linear rescaling. There are multiple published series that show different long-term trends for solar activity. We aim at producing a GSN series, with a non-linear non-parametric calibration. The only underlying assumptions are that the differences between the various series are due to different acuity thresholds of the observers, and that the threshold of each observer remains constant throughout the observing period. We use a daisy chain process with backbone (BB) observers and calibrate all overlapping observers to them. The calibration of each individual observer is performed with a probability distribution function (PDF) matrix, constructed considering all daily values for the overlapping period with the BB. The calibration of the BBs is done in a similar manner. Propagation of errors is modelled with Monte Carlo simulations. The final series extends back to 1739 and includes data from 314 observers. It suggests moderate activity during the 18th and 19th century, which is significantly lower than the high level of solar activity predicted by other recent reconstructions applying linear regressions. The new series provides a robust reconstruction, based on modern and non-parametric methods, of sunspot group numbers since 1739, and it confirms the existence of the Modern grand maximum of solar activity in the second half of the 20th century.

An important topic in the field of space weather is the precise assessment of the contribution of galactic cosmic rays and solar energetic particles on air crew exposure, specifically during eruptive events on the Sun. Here we present a model, a full chain analysis based on ground based measurements of cosmic rays with neutron monitors, subsequent derivation of particle spectral and angular characteristics and computation of dose rate. The model uses method for ground level enhancement analysis and newly numerically computed yield functions for conversion of secondary particle fluence to effective dose and/or the ambient dose equivalent. The precise an adequate information about the solar energetic particle spectra (SEPs) is the basis of the model. Since SEPs possess an essential isotropic part, specifically during the event onset, the angular characteristics should be also derived with good precision. This can be achieved using neutron monitor data during a special class of SEP events the ground level enhancements (GLEs). On the basis of the method representing a sequence of consecutive steps: computation of the NM asymptotic cones, NM rigidity cut-off and application of convenient optimization procedure, we derive the rigidity spectra and anisotropy characteristics of GLE particles. For the computation we use newly computed yield function of the standard sea-level 6NM64 neutron monitor for primary proton and alpha CR nuclei as well as 6NM64 yield function at altitudes ranging from the sea level up to 5000 m above the sea level. We derive the SEP spectra and pitch angle distributions in their dynamical development throughout the event. Subsequently on the basis of the derived spectra and angular characteristics and previously computed yield functions we calculate the effective dose and/or ambient dose equivalent during the GLE. Several examples are shown.

It has been recently claimed (Zolotova and Ponyavin, Solar Phys., 291, 2869, 2016, ZP16 henceforth) that a mid-latitude optical phenomenon, which took place over the city of Astrakhan in July 1670, according to Russian chronicles, was a strong aurora borealis. If this was true, it would imply a very strong or even severe geomagnetic storm during the quietest part of the Maunder minimum. However, as we argue in this article, this conclusion is erroneous and caused by a misinterpretation of the chronicle record. As a result of a thorough analysis of the chronicle text, we show that the described phenomenon occurred during the daylight period of the day ("the last morning hour"), in the south direction ("towards noon"), and its description does not match that of an aurora. The date of the event was also incorrectly interpreted. We conclude that this phenomenon was not a mid-latitude aurora but an atmospheric phenomenon, the so-called sundog (or parhelion) which is a particular type of solar halo. Accordingly, the claim about a strong mid-latitude aurora during the deep Maunder minimum is not correct and should be dismissed.

A ground-based neutron monitor is a standard tool to measure cosmic ray variability near Earth, and it is crucially important to know its yield function for primary cosmic rays. Although there are several earlier theoretically calculated yield functions, none of them agrees with experimental data of latitude surveys of sea-level neutron monitors, thus suggesting for an inconsistency. A newly computed yield function of the standard sea-level 6NM64 neutron monitor is presented here separately for primary cosmic ray protons and $\alpha-$particles, the latter representing also heavier species of cosmic rays. The computations have been done using the GEANT-4 Planetocosmics Monte-Carlo tool and a realistic curved atmospheric model. For the first time, an effect of the geometrical correction of the neutron monitor effective area, related to the finite lateral expansion of the cosmic ray induced atmospheric cascade, is considered, that was neglected in the previous studies. This correction slightly enhances the relative impact of higher-energy cosmic rays (energy above 5--10 GeV/nucleon) in neutron monitor count rate. The new computation finally resolves the long-standing problem of disagreement between the theoretically calculated spatial variability of cosmic rays over the globe and experimental latitude surveys. The newly calculated yield function, corrected for this geometrical factor, appears fully consistent with the experimental latitude surveys of neutron monitors performed during three consecutive solar minima in 1976--77, 1986--87 and 1996--97. Thus, we provide a new yield function of the standard sea-level neutron monitor 6NM64 that is validated against experimental data.

Nearly 70 Ground Level Enhancements (GLEs) of cosmic rays have been recorded by the worldwide neutron monitor network since the 1950s depicting a big variety of energy spectra of solar energetic particles (SEP). Here we studied a statistical relation between the event-integrated intensity of GLEs (calculated as count-rate relative excess, averaged over all available polar neutron monitors, and expressed in percent-hours) and the hardness of the solar particle energy spectra. For each event the integral omnidirectional event-integrated fluences of particles with energy above 30 MeV ($F_{30}$) and above 200 MeV ($F_{200}$) were computed using the reconstructed spectra, and the ratio between the two fluences was considered as a simple index of the event's hardness. We also provided a justification of the spectrum estimate in the form of the Band-function, using direct PAMELA data for GLE 71 (17-May-2012). We found that, while there is no clear relation between the intensity and the hardness for weak events, all strong events with the intensity greater 100 \%*hr are characterized by a very hard spectrum. This implies that a hard spectrum can be securely assumed for all extreme GLE events, e.g., those studied using cosmogenic isotope data in the past.

The distribution of the sunspot group size (area) and its dependence on the level of solar activity is studied. It is shown that the fraction of small groups is not constant but decreases with the level of solar activity so that high solar activity is largely defined by big groups. We study the possible influence of solar activity on the ability of a realistic observer to see and report the daily number of sunspot groups. It is shown that the relation between the number of sunspot groups as seen by different observers with different observational acuity thresholds is strongly non-linear and cannot be approximated by the traditionally used linear scaling ($k-$factors). The observational acuity threshold [$A_{\rm th}$] is considered to quantify the quality of each observer, instead of the traditional relative $k-$factor. A nonlinear $c-$factor based on $A_{\rm th}$ is proposed, which can be used to correct each observer to the reference conditions. The method is tested on a pair of principal solar observers, Wolf and Wolfer, and it is shown that the traditional linear correction, with the constant $k-$factor of 1.66 to scale Wolf to Wolfer, leads to an overestimate of solar activity around solar maxima.

New consistent and precise computations of the production of five cosmogenic radio-isotopes, 7Be, 10Be, 14C, 22Na and 36Cl, in the Earth's atmosphere by cosmic rays are presented in the form of tabulated yield functions. For the first time, a detailed set of the the altitude profiles of the production functions is provided which makes it possible to apply the results directly as input for atmospheric transport models. Good agreement with most of the earlier published works for columnar and global isotopic production rates is shown. Altitude profiles of the production are important, in particular for such tasks as studies of strong solar particle events in the past, precise reconstructions of solar activity on long-term scale, tracing air-mass dynamics using cosmogenic radio-isotopes, etc. As an example, computations of the $^{10}$Be deposition flux in the polar region are shown for the last decades and also for a period around 780 AD and confronted with the actual measurements in Greenland and Antarctic ice cores.

We compare four sunspot-number data sequences against geomagnetic and terrestrial auroral observations. The comparisons are made for the original SIDC composite of Wolf-Zurich-International sunspot number [$R_{ISNv1}$], the group sunspot number [$R_{G}$] by Hoyt and Schatten (Solar Phys., 1998), the new "backbone" group sunspot number [$R_{BB}$] by Svalgaard and Schatten (Solar Phys., 2016), and the "corrected" sunspot number [$R_{C}$] by Lockwood at al. (J.G.R., 2014). Each sunspot number is fitted with terrestrial observations, or parameters derived from terrestrial observations to be linearly proportional to sunspot number, over a 30-year calibration interval of 1982-2012. The fits are then used to compute test sequences, which extend further back in time and which are compared to $R_{ISNv1}$, $R_{G}$, $R_{BB}$, and $R_{C}$. To study the long-term trends, comparisons are made using averages over whole solar cycles (minimum-to-minimum). The test variations are generated in four ways: i) using the IDV(1d) and IDV geomagnetic indices (for 1845-2013) fitted over the calibration interval using the various sunspot numbers and the phase of the solar cycle; ii) from the open solar flux (OSF) generated for 1845 - 2013 from four pairings of geomagnetic indices by Lockwood et al. (Ann. Geophys., 2014) and analysed using the OSF continuity model of Solanki at al. (Nature, 2000) which employs a constant fractional OSF loss rate; iii) the same OSF data analysed using the OSF continuity model of Owens and Lockwood (J.G.R., 2012) in which the fractional loss rate varies with the tilt of the heliospheric current sheet and hence with the phase of the solar cycle; iv) the occurrence frequency of low-latitude aurora for 1780-1980 from the survey of Legrand and Simon (Ann. Geophys., 1987). For all cases, $R_{BB}$ exceeds the test terrestrial series by an amount that increases as one goes back in time.

New sunspot data composites, some of which are radically different in the character of their long-term variation, are evaluated over the interval 1845-2014. The method commonly used to calibrate historic sunspot data, relative to modern-day data, is "daisy-chaining", whereby calibration is passed from one data subset to the neighbouring one, usually using regressions of the data subsets for the intervals of their overlap. Recent studies have illustrated serious pitfalls in these regressions and the resulting errors can be compounded by their repeated use as the data sequence is extended back in time. Hence the recent composite data series by Usoskin et al. (2016), $R_{UEA}$, is a very important advance because it avoids regressions, daisy-chaining and other common, but invalid, assumptions: this is achieved by comparing the statistics of "active day" fractions to those for a single reference dataset. We study six sunspot data series including $R_{UEA}$ and the new "backbone" data series $R_{BB}$, recently generated by Svalgaard and Schatten (2016) by employing both regression and daisy-chaining. We show that all six can be used with a continuity model to reproduce the main features of the open solar flux variation for 1845-2014, as reconstructed from geomagnetic activity data. However, some differences can be identified that are consistent with tests using a basket of other proxies for solar magnetic fields. Using data from a variety of sunspot observers, we illustrate problems with the method employed in $R_{BB}$ which cause it to increasingly overestimate sunspot numbers going back in time and we recommend using $R_{UEA}$ because it employs more robust procedures that avoid such problems.

On the basis of neutron monitor data we estimate the energy spectrum, anisotropy axis direction and pitch-angle distribution of solar energetic particles during two major ground level enhancements (GLE 59 on 14 July 2000 and GLE 70 on 13 December 2006). For the analysis we use a newly computed neutron monitor yield function. The method consists of several consecutive steps: definition of the asymptotic viewing cones of neutron monitor stations considered for the data analysis by computations of cosmic ray particles propagation in a model magnetosphere with the MAGNETOCOSMICS code; computation of the neutron monitor model responses and derivation of the solar energetic particle characteristics on the basis of inverse problem solution. The pitch-angle distribution and rigidity spectrum of high-energy protons are obtained as function of time in the course of ground level enhancements. A comparison with previously reported results is performed and reasonable agreement is achieved. A discussion of the obtained results is included as well their possible application is included.

Cosmogenic isotopes provide the only quantitative proxy for analyzing the long-term solar variability over a centennial timescale. While essential progress has been achieved in both measurements and modeling of the cosmogenic proxy, uncertainties still remain in the determination of the geomagnetic dipole moment evolution. Here we improve the reconstruction of solar activity over the past nine millennia using a multi-proxy approach. We used records of the 14C and 10Be cosmogenic isotopes, current numerical models of the isotope production and transport in Earth's atmosphere, and available geomagnetic field reconstructions, including a new reconstruction relying on an updated archeo-/paleointensity database. The obtained series were analyzed using the singular spectrum analysis (SSA) method to study the millennial-scale trends. A new reconstruction of the geomagnetic dipole field moment, GMAG.9k, is built for the last nine millennia. New reconstructions of solar activity covering the last nine millennia, quantified in sunspot numbers, are presented and analyzed. A conservative list of grand minima and maxima is provided. The primary components of the reconstructed solar activity, as determined using the SSA method, are different for the series based on 14C and 10Be. These primary components can only be ascribed to long-term changes in the terrestrial system and not to the Sun. They have been removed from the reconstructed series. In contrast, the secondary SSA components of the reconstructed solar activity are found to be dominated by a common ~2400-yr quasi-periodicity, the so-called Hallstatt cycle, in both the 14C and 10Be based series. This Hallstatt cycle thus appears to be related to solar activity. Finally, we show that the grand minima and maxima occurred intermittently over the studied period, with clustering near highs and lows of the Hallstatt cycle, respectively.

We use sunspot group observations from the Royal Greenwich Observatory (RGO) to investigate the effects of intercalibrating data from observers with different visual acuities. The tests are made by counting the number of groups $R_B$ above a variable cut-off threshold of observed total whole-spot area (uncorrected for foreshortening) to simulate what a lower acuity observer would have seen. The synthesised annual means of $R_B$ are then re-scaled to the observed RGO group number $R_A$ using a variety of regression techniques. It is found that a very high correlation between $R_A$ and $R_B$ ($r_{AB}$ > 0.98) does not prevent large errors in the intercalibration (e.g. sunspot maximum values can be over 30% too large even for such levels of $r_{AB}$). In generating the backbone sunspot number, Svalgaard and Schatten [2015] force regression fits to pass through the scatter plot origin which generates unreliable fits (the residuals do not form a normal distribution) and causes sunspot cycle amplitudes to be exaggerated in the intercalibrated data. It is demonstrated that the use of Quantile-Quantile (Q-Q) plots to test for a normal distribution is a useful indicator of erroneous and misleading regression fits. Ordinary least squares linear fits, not forced to pass through the origin, are sometimes reliable (although the optimum method used is shown to be different when matching peak and average sunspot group numbers). However other fits are only reliable if non-linear regression is used. From these results it is entirely possible that the inflation of solar cycle amplitudes in the backbone group sunspot number as one goes back in time, relative to related solar-terrestrial parameters, is entirely caused by the use of inappropriate and non-robust regression techniques to calibrate the sunspot data.

A new model CRAC:EPII (Cosmic Ray Atmospheric Cascade: Electron Precipitation Induced Ionization) is presented. The CRAC:EPII is based on Monte Carlo simulation of precipitating electrons propagation and interaction with matter in the Earth atmosphere. It explicitly considers energy deposit: ionization, pair production, Compton scattering, generation of Bremsstrahlung high energy photons, photo-ionization and annihilation of positrons, multiple scattering as physical processes accordingly. The propagation of precipitating electrons and their interactions with atmospheric molecules is carried out with the GEANT4 simulation tool PLANETOCOSMICS code using NRLMSISE 00 atmospheric model. The ionization yields is compared with an analytical parametrization for various energies of incident precipitating electron, using a flux of mono-energetic particles. A good agreement between the two models is achieved. Subsequently, on the basis of balloon-born measured spectra of precipitating electrons at 30.10.2002 and 07.01.2004, the ion production rate in the middle and upper atmosphere is estimated using the CRAC:EPII model

Although the sunspot-number series have existed since the mid-19th century, they are still the subject of intense debate, with the largest uncertainty being related to the "calibration" of the visual acuity of individual observers in the past. Daisy-chain regression methods are applied to inter-calibrate the observers which may lead to significant bias and error accumulation. Here we present a novel method to calibrate the visual acuity of the key observers to the reference data set of Royal Greenwich Observatory sunspot groups for the period 1900-1976, using the statistics of the active-day fraction. For each observer we independently evaluate their observational thresholds [S_S] defined such that the observer is assumed to miss all of the groups with an area smaller than S_S and report all the groups larger than S_S. Next, using a Monte-Carlo method we construct, from the reference data set, a correction matrix for each observer. The correction matrices are significantly non-linear and cannot be approximated by a linear regression or proportionality. We emphasize that corrections based on a linear proportionality between annually averaged data lead to serious biases and distortions of the data. The correction matrices are applied to the original sunspot group records for each day, and finally the composite corrected series is produced for the period since 1748. The corrected series displays secular minima around 1800 (Dalton minimum) and 1900 (Gleissberg minimum), as well as the Modern grand maximum of activity in the second half of the 20th century. The uniqueness of the grand maximum is confirmed for the last 250 years. It is shown that the adoption of a linear relationship between the data of Wolf and Wolfer results in grossly inflated group numbers in the 18th and 19th centuries in some reconstructions.

A two-wave dynamo model was recently proposed by Zharkova et al. (2015, Zh15 henceforth), which aims at long-term predictions of solar activity for millennia ahead and backwards. Here we confront the backward model predictions for the last 800 years with known variability of solar activity, using both direct sunspot observations since 1610 and reconstructions based on cosmogenic radionuclide data. We show that the Zh15 model fails to reproduce the well-established features of the solar activity evolution during the last millennium. This means that the predictive part for the future is not reliable either.

Aims: Although the time of the Maunder minimum (1645--1715) is widely known as a period of extremely low solar activity, claims are still debated that solar activity during that period might still have been moderate, even higher than the current solar cycle #24. We have revisited all the existing pieces of evidence and datasets, both direct and indirect, to assess the level of solar activity during the Maunder minimum. Methods: We discuss the East Asian naked-eye sunspot observations, the telescopic solar observations, the fraction of sunspot active days, the latitudinal extent of sunspot positions, auroral sightings at high latitudes, cosmogenic radionuclide data as well as solar eclipse observations for that period. We also consider peculiar features of the Sun (very strong hemispheric asymmetry of sunspot location, unusual differential rotation and the lack of the K-corona) that imply a special mode of solar activity during the Maunder minimum. Results: The level of solar activity during the Maunder minimum is reassessed on the basis of all available data sets. Conclusions: We conclude that solar activity was indeed at an exceptionally low level during the Maunder minimum. Although the exact level is still unclear, it was definitely below that during the Dalton minimum around 1800 and significantly below that of the current solar cycle #24. Claims of a moderate-to-high level of solar activity during the Maunder minimum are rejected at a high confidence level.

Long-term records of sunspot number and concentrations of cosmogenic radionuclides (10Be and 14C) on the Earth reveal the variation of the Sun's magnetic activity over hundreds and thousands of years. We identify several clear periods in sunspot, 10Be, and 14C data as 1000, 500, 350, 200 and 100 years. We found that the periods of the first five spherical harmonics of the slow magnetic Rossby mode in the presence of a steady toroidal magnetic field of 1200-1300 G in the lower tachocline are in perfect agreement with the time scales of observed variations. The steady toroidal magnetic field can be generated in the lower tachocline either due to the steady dynamo magnetic field for low magnetic diffusivity or due to the action of the latitudinal differential rotation on the weak poloidal primordial magnetic field, which penetrates from the radiative interior. The slow magnetic Rossby waves lead to variations of the steady toroidal magnetic field in the lower tachocline, which modulate the dynamo magnetic field and consequently the solar cycle strength. This result constitutes a key point for long-term prediction of the cycle strength. According to our model, the next deep minimum in solar activity is expected during the first half of this century.

The Maunder minimum (MM) of greatly reduced solar activity took place in 1645-1715, but the exact level of sunspot activity is uncertain as based, to a large extent, on historical generic statements of the absence of spots on the Sun. Here we aim, using a conservative approach, to assess the level and length of solar cycle during the Maunder minimum, on the basis of direct historical records by astronomers of that time. A database of the active and inactive days (days with and without recorded sunspots on the solar disc respectively) is constructed for three models of different levels of conservatism (loose ML, optimum MO and strict MS models) regarding generic no-spot records. We have used the active day fraction to estimate the group sunspot number during the MM. A clear cyclic variability is found throughout the MM with peaks at around 1655--1657, 1675, 1684 and 1705, and possibly 1666, with the active day fraction not exceeding 0.2, 0.3 or 0.4 during the core MM, for the three models. Estimated sunspot numbers are found very low in accordance with a grand minimum of solar activity. We have found, for the core MM (1650-1700), that: (1) A large fraction of no-spot records, corresponding to the solar meridian observations, may be unreliable in the conventional database. (2) The active day fraction remained low (below 0.3-0.4) throughout the MM, indicating the low level of sunspot activity. (3) The solar cycle appears clearly during the core MM. (4) The length of the solar cycle during the core MM appears $9\pm 1$ years, but there is an uncertainty in that. (5) The magnitude of the sunspot cycle during MM is assessed to be below 5-10 in sunspot numbers; A hypothesis of the high solar cycles during the MM is not confirmed.

Context. The relation between solar surface rotation and sunspot activity still remains open. Sunspot activity has dramatically reduced in solar cycle 24 and several solar activity indices and flux measurements experienced unprecedentedly low levels during the last solar minimum. Aims. We aim to reveal the momentary variation of solar surface rotation, especially during the recent years of reducing solar activity. Methods. We used a dynamic, differentially rotating reference system to determine the best-fit annual values of the differential rotation parameters of active longitudes of solar X-ray flares and sunspots in 1977-2012. Results. The evolution of rotation of solar active longitudes obtained with X-ray flares and with sunspots is very similar. Both hemispheres speed up since the late 1990s, with the southern hemisphere rotating slightly faster than the north. Earlier, in 1980s, rotation in the northern hemisphere was considerably faster, but experienced a major decrease in the early 1990s. On the other hand, little change was found in the southern rotation during these decades. This led to a positive asymmetry in north-south rotation rate in the early part of the time interval studied. Conclusions. The rotation of both hemispheres has been speeding up at roughly the same rate since late 1990s, with the southern hemisphere rotating slightly faster than the north. This period coincides with the start of dramatic weakening of solar activity, as observed in sunspots and several other solar, interplanetary and geomagnetic parameters.

Firestone (2014) asserted evidence for numerous (23) nearby (d<300 pc) supernovae within the Middle and Late Pleistocene. If true, this would have strong implications for the irradiation of the Earth; at this rate, mass extinction level events due to supernovae would be more frequent than 100 Myr. However, there are numerous errors in the application of past research. The paper overestimates likely nitrate and 14C production from moderately nearby supernovae by about four orders of magnitude. Moreover, the results are based on wrongly selected (obsolete) nitrate and 14C datasets. The use of correct and up-to-date datasets does not confirm the claimed results. The claims in the paper are invalidated.

A mysterious increase of radiocarbon 14C ca. 775 AD in the Earth's atmosphere has been recently found by Miyake et al. (Nature, 486, 240, 2012). A possible source of this event has been discussed widely, the most likely being an extreme solar energetic particle event. A new exotic hypothesis has been presented recently by Liu et al. (Nature Sci. Rep., 4, 3728, 2014) who proposed that the event was caused by a comet bringing additional 14C to Earth. Here we calculated a realistic mass and size of such a comet to show that it would have been huge (~100 km across and 10^{14}-10^{15} ton of mass) and would have produced a disastrous impact on Earth. Such an impact could not remain unnoticed in the geological records and chronicles. The absence of an evidence for such a dramatic event makes this hypothesis invalid.

Aims. The Sun shows strong variability in its magnetic activity, from Grand minima to Grand maxima, but the nature of the variability is not fully understood, mostly because of the insufficient length of the directly observed solar activity records and of uncertainties related to long-term reconstructions. Here we present a new adjustment-free reconstruction of solar activity over three millennia and study its different modes. Methods. We present a new adjustment-free, physical reconstruction of solar activity over the past three millennia, using the latest verified carbon cycle, 14C production, and archeomagnetic field models. This great improvement allowed us to study different modes of solar activity at an unprecedented level of details. Results. The distribution of solar activity is clearly bi-modal, implying the existence of distinct modes of activity. The main regular activity mode corresponds to moderate activity that varies in a relatively narrow band between sunspot numbers about 20 and 67. The existence of a separate Grand minimum mode with reduced solar activity, which cannot be explained by random fluctuations of the regular mode, is confirmed at a high confidence level. The possible existence of a separate Grand maximum mode is also suggested, but the statistics is too low to reach a confident conclusion. Conclusions. The Sun is shown to operate in distinct modes - a main general mode, a Grand minimum mode corresponding to an inactive Sun, and a possible Grand maximum mode corresponding to an unusually active Sun. These results provide important constraints for both dynamo models of Sun-like stars and investigations of possible solar influence on Earth's climate.

The present work is a critical revision of the hypothesis of the planetary tidal influence on solar activity published by Abreu et al. (Astron. Astrophys. 548, A88, 2012; called A12 here). A12 describes a hypothesis that planets can have an impact on the solar tachocline and therefore on solar activity. We checked the procedure and results of A12, namely the algorithm of planetary tidal torque calculation and the wavelet coherence between torque and heliospheric modulation potential. We found that the claimed peaks in long-period range of the torque spectrum are artefacts caused by the calculation algorithm. Also the statistical significance of the results of the wavelet coherence is found to be overestimated by an incorrect choice of the background assumption of red noise. Using a more conservative non-parametric random phase method, we found that the long-period coherence between planetary torque and heliospheric modulation potential becomes insignificant. Thus we conclude that the considered hypothesis of planetary tidal influence on solar activity is not based on a solid ground.

Aims. Sunspot number is a benchmark series in many studies, but may still contain inhomogeneities and inconsistencies. In particular, an essential discrepancy exists between the two main sunspot number series, Wolf (WSN) and group (GSN) sunspot numbers, before 1848. The source of this discrepancy has so far remained unresolved. However, the recently digitized series of solar observations in 1825-1867 by Samuel Heinrich Schwabe, who was the primary observer of the WSN before 1848, makes such an assessment possible. Methods. We construct sunspot series, similar to WSN and GSN, but using only Schwabe's data. These series, called WSN-S and GSN-S, respectively, were compared with the original WSN and GSN series for the period 1835-1867 to look for possible inhomogeneities. Results. We show that: (1) The GSN series is homogeneous and consistent with the Schwabe data throughout the entire studied period; (2) The WSN series decreases by roughly ~20% around 1848 caused by the change of the primary observer from Schwabe to Wolf and an inappropriate individual correction factor used for Schwabe in the WSN; (3) This implies a major inhomogeneity in the WSN, which needs to be corrected by reducing its values by 20% before 1848; (4) The corrected WSN series is in good agreement with the GSN series. This study supports the earlier conclusions that the GSN series is more consistent and homogeneous in the earlier part than the WSN series.

Samuel Heinrich Schwabe made 8486 drawings of the solar disk with sunspots in the period from November 5, 1825 to December 29, 1867. We have measured sunspot sizes and heliographic positions on digitized images of these drawings. A total of about 135,000 measurements of individual sunspots are available in a data base. Positions are accurate to about 5% of the solar radius or to about three degrees in heliographic coordinates in the solar disk center. Sizes were given in 12 classes as estimated visually with circular cursor shapes on the screen. Most of the drawings show a coordinate grid aligned with the celestial coordinate system. A subset of 1168 drawings have no indication of their orientation. We have used a Bayesian estimator to infer the orientations of the drawings as well as the average heliographic spot positions from a chain of drawings of several days, using the rotation profile of the present Sun. The data base also includes all information available from Schwabe on spotless days.

We revisited assessments of the occurrence probability distribution of large events in solar energetic particles (SEP), based on measurements of cosmogenic radionuclides in lunar rocks. We present a combined cumulative occurrence probability distribution of SEP events based on three time scales: directly measured SEP fluences for the last 60 years; estimates based on terrestrial cosmogenic radionuclides 10Be and 14C for the multi-millennial (Holocene) time scale; and cosmogenic radionuclides measured in lunar rocks on the time scale of up to 1 Myr. All the three time scales yield a consistent distribution. The data suggest a strong rollover of the occurrence probability so that SEP events with the fluence of protons with energy >30 MeV greater than 10^{11} (protons /cm2/yr) are not expected at the Myr time scale.

Miyake et al. (henceforth M12) recently reported, based on 14C data, an extreme cosmic event ca. AD775. Using a simple model, M12 claimed that the event was too strong to be caused by a solar flare within the standard theory. This implied a new paradigm of either an impossibly strong solar flare or a very strong cosmic ray event of unknown origin occurred ca. AD775. We show that the strength of the event was significantly overestimated by M12. Several subsequent works have attempted to find a possible exotic source for such an event, but they are all based on incorrect estimates by M12. We revisit this event with analysis of new datasets and consistent theoretical modelling. We verified the experimental result for the AD775 event using independent datasets including 10Be series and newly measured 14C annual data. We surveyed available historical chronicles for astronomical observations for the AD770s to identify potential sightings of aurorae or supernovae. We interpreted the 14C measurements using an appropriate carbon cycle model. We show that: (1) The reality of the AD775 event is confirmed by new measurements of 14C; (2) by using an inappropriate carbon cycle model, M12 strongly overestimated the event's strength; (3) The revised magnitude of the event is consistent with different independent datasets (14C, 10Be, 36Cl) and can be associated with a strong, but not inexplicably strong, SEP event (or a sequence of events), and provides the first evidence for an event of this magnitude (the fluence >30 MeV was about 4.5*10^{10} /cm2) in multiple datasets; (4) This is in agreement with increased auroral activity identified in historical chronicles. This point to the likely solar origin of the event, which is the greatest solar event on a multi-millennial time scale, placing a strong observational constraint on the theory of explosive energy releases on the Sun and cool stars.

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 \super{14}C and \super{10}Be 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.

We report on the 2012 May 17 Ground Level Enhancement (GLE) event, which is the first of its kind in Solar Cycle 24. This is the first GLE event to be fully observed close to the surface by the Solar Terrestrial Relations Observatory (STEREO) mission. We determine the coronal mass ejection (CME) height at the start of the associated metric type II radio burst (i.e., shock formation height) as 1.38 Rs (from the Sun center). The CME height at the time of GLE particle release was directly measured from a STEREO image as 2.32 Rs, which agrees well with the estimation from CME kinematics. These heights are consistent with those obtained for cycle-23 GLEs using back-extrapolation. By contrasting the 2012 May 17 GLE with six other non-GLE eruptions from well-connected regions with similar or larger flare size and CME speed, we find that the latitudinal distance from the ecliptic is rather large for the non-GLE events due to a combination of non-radial CME motion and unfavorable solar B0 angle, making the connectivity to Earth poorer. We also find that the coronal environment may play a role in deciding the shock strength.

The long-term variability of the sunspot cycle, as recorded by the Wolf numbers, are imprinted in different kinds of statistical relations which relate the cycle amplitudes, duration and shapes. This subject always gets a special attention because it is important for the solar activity forecast. We discuss statistical properties of the mean-field dynamo model with the fluctuating \alpha-effect. Also, we estimate dynamical properties of the model for the long and short time-scale and compare it with the dynamics of the sunspot numbers data sets.

The probability of occurrence of extreme solar particle events (SPEs) with the fluence of (>30 MeV) protons F30>10^{10} cm^{-2} is evaluated based on data of cosmogenic isotopes 14C and 10Be in terrestrial archives centennial-millennial time scales. Four potential candidates with F30=(1-1.5)x10^{10} cm^{-2} and no events with F30>2x10^{10} cm^{-2} are identified since 1400 AD in the annually resolved 10Be data. A strong SPE related to the Carrington flare of 1859 AD is not supported by the data. For the last 11400 years, 19 SPE candidates with F30=(1-3)x10^{10} cm^{-2} are found and clearly no event with F30>5x10^{10} cm^{-2} (50-fold the SPE of 23-Feb-1956) occurring. This values serve as an observational upper limit for the strength of SPE on the time scale of tens of millennia. Two events, ca. 780 and 1460 AD, appear in different data series making them strong candidates to extreme SPEs. We built a distribution of the occurrence probability of extreme SPEs, providing a new strict observational constraint. Practical limits can be set as F30~1x, 2-3x, and 5x10^{10} cm^{-2} for the occurrence probability ~10^{-2}, 10^{-3} and 10^{-4} year^{-1}, respectively. Because of uncertainties, our results should be interpreted as a conservative upper limit of the SPE occurrence near Earth. The mean SEP flux is evaluated as ~40 (cm2 sec)^{-1} in agreement with estimates from the lunar rocks. On average, extreme SPEs contribute about 10% to the total SEP fluence.

We present the results of full new calculation of radiocarbon 14C production in the Earth atmosphere, using a numerical Monte-Carlo model. We provide, for the first time, a tabulated 14C yield function for the energy of primary cosmic ray particles ranging from 0.1 to 1000 GeV/nucleon. We have calculated the global production rate of 14C, which is 1.64 and 1.88 atoms/cm2/s for the modern time and for the pre-industrial epoch, respectively. This is close to the values obtained from the carbon cycle reservoir inventory. We argue that earlier models overestimated the global 14C production rate because of outdated spectra of cosmic ray heavier nuclei. The mean contribution of solar energetic particles to the global 14C is calculated as about 0.25% for the modern epoch. Our model provides a new tool to calculate the 14C production in the Earth's atmosphere, which can be applied, e.g., to reconstructions of solar activity in the past.

The most powerful explosions on the Sun [...] drive the most severe space-weather storms. Proxy records of flare energies based on SEPs in principle may offer the longest time base to study infrequent large events. We conclude that one suggested proxy, nitrate concentrations in polar ice cores, does not map reliably to SEP events. Concentrations of select radionuclides measured in natural archives may prove useful in extending the time interval of direct observations up to ten millennia, but as their calibration to solar flare fluences depends on multiple poorly known properties and processes, these proxies cannot presently be used to help determine the flare energy frequency distribution. Being thus limited to the use of direct flare observations, we evaluate the probabilities of large-energy solar explosions by combining solar flare observations with an ensemble of stellar flare observations. We conclude that solar flare energies form a relatively smooth distribution from small events to large flares, while flares on magnetically-active, young Sun-like stars have energies and frequencies markedly in excess of strong solar flares, even after an empirical scaling with the mean activity level of these stars. In order to empirically quantify the frequency of uncommonly large solar flares extensive surveys of stars of near-solar age need to be obtained, such as is feasible with the Kepler satellite. Because the likelihood of flares larger than approximately X30 remains empirically unconstrained, we present indirect arguments, based on records of sunspots and on statistical arguments, that solar flares in the past four centuries have likely not substantially exceeded the level of the largest flares observed in the space era, and that there is at most about a 10% chance of a flare larger than about X30 in the next 30 years.

Solar cycles vary in their amplitude and shape. There are several empirical relations between various parameters linking cycle's shape and amplitude, in particular the Waldmeier relations. As solar cycle is believed to be a result of the solar dynamo action, these relations require explanation in the framework of this theory.Here we aim to present a possible explanation of such kind. We relate the cycle-to-cycle variability of solar activity to fluctuations of solar dynamo drivers and primarily to fluctuations in the parameter responsible for recovery of the poloidal magnetic field from the toroidal one. To be specific, we develop such a model in the framework of the mean-field dynamo based on the differential rotation and $\alpha$-effect. We demonstrate that the mean-field dynamo based on a realistic rotation curve and nonlinearity associated with the magnetic helicity balance reproduces both qualitatively and quantitatively the Waldmeier relations observed in sunspot data since 1750 (SIDC data). The model also reproduces more or less successfully other relations between the parameters under discussion, in particular, the link between odd and even cycles (Gnevyshev-Ohl rule). We conclude that the contemporary solar dynamo theory provides a way to explain the cycle-to-cycle variability of solar activity as recorded in sunspots. We discuss the importance of the model for stellar activity cycles which, as known from the data of HK project, demonstrate the cycle-to-cycle variability similar to solar cycles.

We present an overview of the observed properties of the GLEs and those of the associated flares and CMEs. The solar eruptions are very intense involving X-class flares and extreme CME speeds (average ~2000 km/s). The active regions in which the GLE events originate are generally large: 1290 msh (median 1010 msh) compared to 934 msh (median: 790 msh) for SEP-producing active regions. The initial acceleration of GLE-associated CMEs is much larger (by a factor of 2) than that of ordinary CMEs (2.3 km/s2 vs.1 km/s2). The GLE particle release is delayed with respect to the onset of all electromagnetic signatures of the eruptions: type II bursts, low frequency type III bursts, soft X-ray flares and CMEs. The presence of metric type II radio bursts some 17 min (median: 16 min; range: 3 to 48 min) before the GLE onset indicates shock formation well before the particle release. The release of GLE particles occurs when the CMEs reach an average height of ~3.09 Rs for well-connected events. For poorly connected events, the average CME height at GLE particle release is ~66% larger (mean: 5.18 Rs). The longitudinal dependence is consistent with shock accelerations because the shocks from poorly connected events need to expand more to cross the field lines connecting to an Earth observer. The CME height at metric type II burst onset is in the narrow range 1.29 to 1.8 Rs, with A mean of 1.53 Rs. The CME heights at metric type II burst onset and GLE particle release correspond to the minimum and maximum in the Alfven speed profile. The CME heights at GLE particle release are in good agreement with those obtained from the velocity dispersion analysis (Reames, 2009a,b) including the source longitude dependence. We also discuss the implications of the delay of GLE particle release with respect to complex type III bursts and hard X-ray emission.

Aims. We present a physically consistent reconstruction of the total solar irradiance for the Holocene. Methods. We extend the SATIRE models to estimate the evolution of the total (and partly spectral) solar irradiance over the Holocene. The basic assumption is that the variations of the solar irradiance are due to the evolution of the dark and bright magnetic features on the solar surface. The evolution of the decadally averaged magnetic flux is computed from decadal values of cosmogenic isotope concentrations recorded in natural archives employing a series of physics-based models connecting the processes from the modulation of the cosmic ray flux in the heliosphere to their record in natural archives. We then compute the total solar irradiance (TSI) as a linear combination of the jth and jth + 1 decadal values of the open magnetic flux. Results. Reconstructions of the TSI over the Holocene, each valid for a di_erent paleomagnetic time series, are presented. Our analysis suggests that major sources of uncertainty in the TSI in this model are the heritage of the uncertainty of the TSI since 1610 reconstructed from sunspot data and the uncertainty of the evolution of the Earth's magnetic dipole moment. The analysis of the distribution functions of the reconstructed irradiance for the last 3000 years indicates that the estimates based on the virtual axial dipole moment are significantly lower at earlier times than the reconstructions based on the virtual dipole moment. Conclusions. We present the first physics-based reconstruction of the total solar irradiance over the Holocene, which will be of interest for studies of climate change over the last 11500 years. The reconstruction indicates that the decadally averaged total solar irradiance ranges over approximately 1.5 W/m2 from grand maxima to grand minima.

Maunder Minimum forms an archetype for the Grand minima, and detailed knowledge of its temporal development has important consequences for the solar dynamo theory dealing with long-term solar activity evolution. Here we reconsider the current paradigm of the Grand minimum general scenario by using newly recovered sunspot observations by G. Marcgraf and revising some earlier uncertain data for the period 1636--1642, i.e., one solar cycle before the beginning of the Maunder Minimum. The new and revised data dramatically change the magnitude of the sunspot cycle just before the Maunder Minimum, from 60--70 down to about 20, implying a possibly gradual onset of the Minimum with reduced activity started two cycles before it. This revised scenario of the Maunder Minimum changes, through the paradigm for Grand solar/stellar activity minima, the observational constraint on the solar/stellar dynamo theories focused on long-term studies and occurrence of Grand minima.

Because of the lack of reliable sunspot observation, the quality of sunspot number series is poor in the late 18th century, leading to the abnormally long solar cycle (1784--1799) before the Dalton minimum. Using the newly recovered solar drawings by the 18--19th century observers Staudacher and Hamilton, we construct the solar butterfly diagram, i.e. the latitudinal distribution of sunspots in the 1790's. The sudden, systematic occurrence of sunspots at high solar latitudes in 1793--1796 unambiguously shows that a new cycle started in 1793, which was lost in traditional Wolf's sunspot series. This finally confirms the existence of the lost cycle that has been proposed earlier, thus resolving an old mystery. This letter brings the attention of the scientific community to the need of revising the sunspot series in the 18th century. The presence of a new short, asymmetric cycle implies changes and constraints to sunspot cycle statistics, solar activity predictions, solar dynamo theories as well as for solar-terrestrial relations.

We consider to what extent the long-term dynamics of cyclic solar activity in the form of Grand Minima can be associated with random fluctuations of the parameters governing the solar dynamo. We consider fluctuations of the alpha-coefficient in the conventional Parker migratory dynamo, and also in slightly more sophisticated dynamo models, and demonstrate that they can mimic the gross features of the phenomenon of the occurrence of Grand Minima over a suitable parameter range. The temporal distribution of these Grand Minima appears chaotic, with a more or less exponential waiting time distribution, typical of Poisson processes. In contrast however, the available reconstruction of Grand Minima statistics based on cosmogenic isotope data demonstrates substantial deviations from this exponential law. We were unable to reproduce the non-Poissonic tail of the waiting time distribution either in the framework of a simple alpha-quenched Parker model, or in its straightforward generalization, nor in simple models with feedback on the differential rotation. We suggest that the disagreement may only be apparent and is plausibly related to the limited observational data, and that the observations and results of numerical modeling can be consistent and represent physically similar dynamo regimes.

Using a reconstruction of sunspot numbers stretching over multiple millennia, we analyze the statistics of the occurrence of grand minima and maxima and set new observational constraints on long-term solar and stellar dynamo models. We present an updated reconstruction of sunspot number over multiple millennia, from $^{14}$C data by means of a physics-based model, using an updated model of the evolution of the solar open magnetic flux. A list of grand minima and maxima of solar activity is presented for the Holocene (since 9500 BC) and the statistics of both the length of individual events as well as the waiting time between them are analyzed. The occurrence of grand minima/maxima is driven not by long-term cyclic variability, but by a stochastic/chaotic process. The waiting time distribution of the occurrence of grand minima/maxima deviates from an exponential distribution, implying that these events tend to cluster together with long event-free periods between the clusters. Two different types of grand minima are observed: short (30--90 years) minima of Maunder type and long ($>$110 years) minima of Sp\"orer type, implying that a deterministic behaviour of the dynamo during a grand minimum defines its length. The duration of grand maxima follows an exponential distribution, suggesting that the duration of a grand maximum is determined by a random process. These results set new observational constraints upon the long-term behaviour of the solar dynamo.

A new type of cosmic-ray experiment is under construction in the Pyh\"asalmi mine in the underground laboratory of the University of Oulu, Finland. It aims to study the composition of cosmic rays at and above the knee region. The experiment, called EMMA, will cover approximately 150 square-metres of detector area. The array is capable of measuring the multiplicity and the lateral distribution of underground muons, and the arrival direction of the air shower. The full-size detector is expected to run by the end of 2007.

As recently found, the distribution of sunspots is non-axisymmetric and spot group formation implies the existence of two persistent active longitudes separated by 180 degrees. Here we quantitatively study the non-axisymmetry of sunspot occurrence. In a dynamic reference frame inferred from the differential rotation law, the raw sunspot data show a clear clustering around the persistent active longitudes. The differential rotation describing the dynamic frame is quantified in terms of the equatorial angular velocity and the differential rotation rate, which appear to be significantly different from those for individual sunspots. This implies that the active longitudes are not linked to the depth of sunspot anchoring. In order to quantify the observed effect, we introduce a measure of the non-axisymmetry of the sunspot distribution. The non-axisymmetric component is found to be highly significant, and the ratio of its strength to that of the axisymmetric one is roughly 1:10. This provides additional constraints for solar dynamo models.

The extension of the sunspot number series backward in time is of considerable interest for dynamo theory, solar, stellar, and climate research. We have used records of the Be-10 concentration in polar ice to reconstruct the average sunspot activity level for the period between the year 850 to the pr esent. Our method uses physical models for processes connecting the Be-10 concentration with the sunspot number. The reconstruction shows reliably that the period of high solar activity during the last 60 years is unique throughout the past 1150 years. This nearly triples the time interval for which such a statement could be made previously

In a recent paper (Usoskin et al., 2002a), we have reconstructed the concentration of the cosmogenic $^{10}$Be isotope in ice cores from the measured sunspot numbers by using physical models for $^{10}$Be production in the Earth's atmosphere, cosmic ray transport in the heliosphere, and evolution of the Sun's open magnetic flux. Here we take the opposite route: starting from the $^{10}$Be concentration measured in ice cores from Antarctica and Greenland, we invert the models in order to reconstruct the 11-year averaged sunspot numbers since 850 AD. The inversion method is validated by comparing the reconstructed sunspot numbers with the directly observed sunspot record since 1610. The reconstructed sunspot record exhibits a prominent period of about 600 years, in agreement with earlier observations based on cosmogenic isotopes. Also, there is evidence for the century scale Gleissberg cycle and a number of shorter quasi-periodicities whose periods seem to fluctuate in the millennium time scale. This invalidates the earlier extrapolation of multi-harmonic representation of sunspot activity over extended time intervals.