Search

A neutron monitor (NM) is, since the 1950s, a standard ground-based detector whose count rate reflects cosmic-ray variability. The worldwide network of NMs forms a rough spectrometer for cosmic rays. Recently, a brand-new data acquisition (DAQ) system has been installed on the DOMC and DOMB NMs, located at the Concordia research station on the Central Antarctic plateau. The new DAQ system digitizes, at a 2-MHz sampling rate, and records all individual pulses corresponding to secondary particles in the detector. An analysis of the pulse characteristics (viz. shape, magnitude, duration, waiting time) has been performed, and several clearly distinguishable branches were identified: (A) corresponding to signal from individual secondary neutrons; (B) representing the detector's noise; (C) double pulses corresponding to the {shortly separated nucleons of the same} atmospheric cascades; (D) very-high multiple pulses which are likely caused by atmospheric muons; and (E) double pulses potentially caused by contamination of the neighbouring detector. An analysis of the waiting time distributions has revealed two clearly distinguishable peaks: peak (I) at about 1 msec being related to the intra-cascade diffusion and thermalisation of secondary atmospheric neutrons; and peak (II) at 30--1000 msec corresponding to individual atmospheric cascades. This opens a new possibility to study spectra of cosmic-ray particles in a single location as well as details of the cosmic-ray induced atmospheric cascades, using the same dataset., Comment: Accepted to JGR Space Physics (doi: 10.1029/2020JA028959)

Fluences of solar energetic particles (SEPs) are not easy to evaluate, especially for high-energy events (i.e. ground-level enhancements, GLEs). Earlier estimates of event-integrated SEP fluences for GLEs were based on partly outdated assumptions and data, and they required revisions. Here, we present the results of a full revision of the spectral fluences for most major SEP events (GLEs) for the period from 1956 -- 2017 using updated low-energy flux estimates along with greatly revisited high-energy flux data and applying the newly invented reconstruction method including an improved neutron-monitor yield function. Low- and high-energy parts of the SEP fluence were estimated using a revised space-borne/ionospheric data and ground-based neutron monitors, respectively. The measured data were fitted by the modified Band function spectral shape. The best-fit parameters and their uncertainties were assessed using a direct Monte Carlo method. As a result, a full reconstruction of the event-integrated spectral fluences was performed in the energy range above 30 MeV, parametrised, and tabulated for easy use along with estimates of the 68% confidence intervals. This forms a solid basis for more precise studies of the physics of solar eruptive events and the transport of energetic particles in the interplanetary medium, as well as the related applications., Comment: 19 pages, 3 figures, to be published in Astronomy and Astrophysics

Empirical rules of solar cycle evolution form important observational constraints for the solar dynamo theory. This includes the Waldmeier rule relating the magnitude of a solar cycle to the length of its ascending phase, and the Gnevyshev--Ohl rule clustering cycles to pairs of an even-numbered cycle followed by a stronger odd-numbered cycle. These rules were established as based on the "classical" Wolf sunspot number series, which has been essentially revisited recently, with several revised sets released by the research community. Here we test the robustness of these empirical rules for different sunspot (group) series for the period 1749--1996, using four classical and revised international sunspot numbers and group sunspot-number series. We also provide an update of the sunspot group series based on the active-day fraction (ADF) method, using the new database of solar observations. We show that the Waldmeier rule is robust and independent of the exact sunspot (group) series: its classical and $n+1$ (relating the length of $n$-th cycle to the magnitude of ($n+1$)-th cycle) formulations are significant or highly significant for all series, while its simplified formulation (relating the magnitude of a cycle to its full length) is insignificant for all series. The Gnevyshev--Ohl rule was found robust for all analyzed series for Cycles 8--21, but unstable across the Dalton minimum and before it., Comment: Solar Physics, accepted

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., Comment: J. Geophys. Res. in press

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., Comment: Solar Physics, accepted

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., Comment: Published in Astron. Astrophys

A new reconstruction of the heliospheric modulation potential for galactic cosmic rays is presented for the neutron monitor era, since 1951. The new reconstruction is based on an updated methodology in comparison to previous reconstructions: (1) the use of the new-generation neutron monitor yield function, (2) the use of the new model of the local interstellar spectrum, employing in particular direct data from the distant missions, and (3) the calibration of the neutron monitor responses to direct measurements of the cosmic ray spectrum performed by the PAMELA space-borne spectrometer over 47 time intervals during 2006{2010. The reconstruction is based on data from six standard NM64-type neutron monitors (Apatity, Inuvik, Kergulen, Moscow, Newark and Oulu) since 1965, and two IGY-type ground-based detectors (Climax and Mt.Washington) for 1951-1964. The new reconstruction, along with the estimated uncertainties is tabulated in the paper. The presented series forms a benchmark record of the cosmic ray variability (in the energy range between 1-30 GeV) for the last 60 years, and can be used in long-term studies in the fields of solar, heliospheric and solar-terrestrial physics., Comment: Published in JGR

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., Comment: 19 pages, 12 figures, Astronomy and Astrophysics, accepted

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., Comment: in press (J. Geophys. Res. (Atmosphere), doi: 10.1002/2016JD025034)

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., Comment: Accepted to Astron. Astrophys

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., Comment: Icarus (accepted)

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., Comment: Solar Physics, in press

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., Comment: Published by Astrophys. J. Lett