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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.