Showing 1–50 of 74 results for author: Usoskin, I

Solar Cycle Observations

Authors: Aimee Norton, Rachel Howe, Lisa Upton, Ilya Usoskin

Abstract: We describe the defining observations of the solar cycle that provide constraints for the dynamo processes operating within the Sun. Specifically, we report on the following topics: historical sunspot numbers and revisions; active region (AR) flux ranges and lifetimes; bipolar magnetic region tilt angles; Hale and Joy's law; the impact of rogue ARs on cycle progression and the amplitude of the fol… ▽ More We describe the defining observations of the solar cycle that provide constraints for the dynamo processes operating within the Sun. Specifically, we report on the following topics: historical sunspot numbers and revisions; active region (AR) flux ranges and lifetimes; bipolar magnetic region tilt angles; Hale and Joy's law; the impact of rogue ARs on cycle progression and the amplitude of the following cycle; the spatio-temporal emergence of ARs that creates the butterfly diagram; polar fields; large-scale flows including zonal, meridional, and AR in-flows; short-term cycle variability; and helioseismic results including mode parameter changes. △ Less

Submitted 28 September, 2023; v1 submitted 31 May, 2023; originally announced May 2023.

Comments: 39 pages, 11 figures, ISSI workshop on Solar and Stellar Dynamos (June 2022). Submitted to Space Science Reviews, May 2023. This is the revised version accepted September 2023

Long-term modulation of solar cycles

Authors: Akash Biswas, Bidya Karak, Ilya Usoskin, Eckhard Weisshaar

Abstract: Solar activity has a cyclic nature with the ~11-year Schwabe cycle dominating its variability on the interannual timescale. However, solar cycles are significantly modulated in length, shape and magnitude, from near-spotless grand minima to very active grand maxima. The ~400-year-long direct sunspot-number series is inhomogeneous in quality and too short to study robust parameters of long-term sol… ▽ More Solar activity has a cyclic nature with the ~11-year Schwabe cycle dominating its variability on the interannual timescale. However, solar cycles are significantly modulated in length, shape and magnitude, from near-spotless grand minima to very active grand maxima. The ~400-year-long direct sunspot-number series is inhomogeneous in quality and too short to study robust parameters of long-term solar variability. The cosmogenic-isotope proxy extends the timescale to twelve millennia and provides crucial observational constraints of the long-term solar dynamo modulation. Here, we present a brief up-to-date overview of the long-term variability of solar activity at centennial--millennial timescales. The occurrence of grand minima and maxima is discussed as well as the existing quasi-periodicities such as centennial Gleissberg, 210-year Suess/de Vries and 2400-year Hallstatt cycles. It is shown that the solar cycles contain an important random component and have no clock-like phase locking implying a lack of long-term memory. A brief yet comprehensive review of the theoretical perspectives to explain the observed features in the framework of the dynamo models is presented, including the nonlinearity and stochastic fluctuations in the dynamo. We keep gaining knowledge of the processes driving solar variability with the new data acquainted and new models developed. △ Less

Submitted 28 February, 2023; originally announced February 2023.

Comments: Accepted for Space Science Review

Re-calibration of the Sunspot Number: Status Report

Authors: F. Clette, L. Lefèvre, T. Chatzistergos, H. Hayakawa, V. M. Carrasco, R. Arlt, E. W. Cliver, T. Dudok de Wit, T. Friedli, N. Karachik, G. Kopp, M. Lockwood, S. Mathieu, A. Muñoz-Jaramillo, M. Owens, D. Pesnell, A. Pevtsov, L. Svalgaard, I. G. Usoskin, L. van Driel-Gesztelyi, J. M. Vaquero

Abstract: We report progress on the ongoing recalibration of the Wolf sunspot number (SN) and Group sunspot number (GN) following the release of version 2.0 of SN in 2015. This report constitutes both an update of the efforts reported in the 2016 Topical Issue of Solar Physics and a summary of work by the International Space Science Institute (ISSI) International Team formed in 2017 to develop optimal SN an… ▽ More We report progress on the ongoing recalibration of the Wolf sunspot number (SN) and Group sunspot number (GN) following the release of version 2.0 of SN in 2015. This report constitutes both an update of the efforts reported in the 2016 Topical Issue of Solar Physics and a summary of work by the International Space Science Institute (ISSI) International Team formed in 2017 to develop optimal SN and GN re-construction methods while continuing to expand the historical sunspot number database. Significant progress has been made on the database side while more work is needed to bring the various proposed SN and (primarily) GN reconstruction methods closer to maturity, after which the new reconstructions (or combinations thereof) can be compared with (a) ``benchmark'' expectations for any normalization scheme (e.g., a general increase in observer normalization factors going back in time), and (b) independent proxy data series such as F10.7 and the daily range of variations of Earth's undisturbed magnetic field. New versions of the underlying databases for SN and GN will shortly become available for years through 2022 and we anticipate the release of next versions of these two time series in 2024. △ Less

Submitted 5 January, 2023; originally announced January 2023.

Comments: 21 figures, 4 tables. To be published in Solar Physics

Exploring the Solar Poles: The Last Great Frontier of the Sun

Authors: Dibyendu Nandy, Dipankar Banerjee, Prantika Bhowmik, Allan Sacha Brun, Robert H. Cameron, S. E. Gibson, Shravan Hanasoge, Louise Harra, Donald M. Hassler, Rekha Jain, Jie Jiang, Laurène Jouve, Duncan H. Mackay, Sushant S. Mahajan, Cristina H. Mandrini, Mathew Owens, Shaonwita Pal, Rui F. Pinto, Chitradeep Saha, Xudong Sun, Durgesh Tripathi, Ilya G. Usoskin

Abstract: Despite investments in multiple space and ground-based solar observatories by the global community, the Sun's polar regions remain unchartered territory - the last great frontier for solar observations. Breaching this frontier is fundamental to understanding the solar cycle - the ultimate driver of short-to-long term solar activity that encompasses space weather and space climate. Magnetohydrodyna… ▽ More Despite investments in multiple space and ground-based solar observatories by the global community, the Sun's polar regions remain unchartered territory - the last great frontier for solar observations. Breaching this frontier is fundamental to understanding the solar cycle - the ultimate driver of short-to-long term solar activity that encompasses space weather and space climate. Magnetohydrodynamic dynamo models and empirically observed relationships have established that the polar field is the primary determinant of the future solar cycle amplitude. Models of solar surface evolution of tilted active regions indicate that the mid to high latitude surges of magnetic flux govern dynamics leading to the reversal and build-up of polar fields. Our theoretical understanding and numerical models of this high latitude magnetic field dynamics and plasma flows - that are a critical component of the sunspot cycle - lack precise observational constraints. This limitation compromises our ability to observe the enigmatic kilo Gauss polar flux patches and constrain the polar field distribution at high latitudes. The lack of these observations handicap our understanding of how high latitude magnetic fields power polar jets, plumes, and the fast solar wind that extend to the boundaries of the heliosphere and modulate solar open flux and cosmic ray flux within the solar system. Accurate observation of the Sun's polar regions, therefore, is the single most outstanding challenge that confronts Heliophysics. This paper argues the scientific case for novel out of ecliptic observations of the Sun's polar regions, in conjunction with existing, or future multi-vantage point heliospheric observatories. Such a mission concept can revolutionize the field of Heliophysics like no other mission concept has - with relevance that transcends spatial regimes from the solar interior to the heliosphere. △ Less

Submitted 30 December, 2022; originally announced January 2023.

Comments: This White Paper was submitted in 2022 to the United States National Academies Solar and Space Physics (Heliophysics) Decadal Survey

Superflares on solar-like stars: A new method for identifying the true flare sources in photometric surveys

Authors: Valeriy Vasilyev, Timo Reinhold, Alexander I. Shapiro, Natalie A. Krivova, Ilya Usoskin, Benjamin T. Montet, Sami K. Solanki, Laurent Gizon

Abstract: Over the past years, thousands of stellar flares have been detected by harvesting data from large photometric surveys. These detections, however, do not account for potential sources of contamination such as background stars appearing in the same aperture as the primary target. We present a new method for identifying the true flare sources in large photometric surveys using data from the Kepler mi… ▽ More Over the past years, thousands of stellar flares have been detected by harvesting data from large photometric surveys. These detections, however, do not account for potential sources of contamination such as background stars appearing in the same aperture as the primary target. We present a new method for identifying the true flare sources in large photometric surveys using data from the Kepler mission. Potential flares are identified in two steps: first, we search the light curves for at least two subsequent data points exceeding a 5σ threshold above the running mean. For these two cadences, we subtract the "quiet" stellar flux from the Kepler pixel data to obtain new images where the potential flare is the main light source. In the second step, we use a Bayesian approach to fit the point spread function of the instrument to determine the most likely location of the flux excess on the detector. We applied our method to 5862 main-sequence stars with near-solar effective temperatures. We found 2274 events exceeding the 5-sigma in at least two consecutive points in the light curves. Applying the second step reduced this number to 342 superflares. Of these, 283 flares happened on 178 target stars, 47 events are associated with fainter background stars, and in 10 cases, the flare location cannot be distinguished between the target and a background star. We also present cases where flares have been reported previously but our technique could not attribute them to the target star. We conclude that 1) identifying outliers in the light curves alone is insufficient to attribute them to stellar flares and 2) flares can only be uniquely attributed to a certain star when the instrument pixel-level data together with the point spread function are taken into account. As a consequence, previous flare statistics are likely contaminated by instrumental effects and unresolved astrophysical sources. △ Less

Submitted 28 September, 2022; originally announced September 2022.

Comments: accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 668, A167 (2022)

Long-term forcing of Sun's coronal field, open flux and cosmic ray modulation potential during grand minima, maxima and regular activity phases by the solar dynamo mechanism

Authors: Soumyaranjan Dash, Dibyendu Nandy, Ilya Usoskin

Abstract: Magnetic fields generated in the Sun's interior by the solar dynamo mechanism drive solar activity over a range of time-scales. While space-based observations of the Sun's corona exist only for few decades, direct sunspot observations exist for a few centuries, solar open flux and cosmic ray flux variations can be reconstructed through studies of cosmogenic isotopes over thousands of years. While… ▽ More Magnetic fields generated in the Sun's interior by the solar dynamo mechanism drive solar activity over a range of time-scales. While space-based observations of the Sun's corona exist only for few decades, direct sunspot observations exist for a few centuries, solar open flux and cosmic ray flux variations can be reconstructed through studies of cosmogenic isotopes over thousands of years. While such reconstructions indicate the presence of extreme solar activity fluctuations in the past, causal links between millennia scale dynamo activity, consequent coronal field, solar open flux and cosmic ray modulation remain elusive. By utilizing a stochastically forced solar dynamo model we perform long-term simulations to illuminate how the dynamo generated magnetic fields govern the structure of the solar corona and the state of the heliosphere -- as indicated by variations in the open flux and cosmic ray modulation potential. We establish differences in the nature of the large-scale structuring of the solar corona during grand maximum, minimum, and regular solar activity phases and simulate how the open flux and cosmic ray modulation potential varies over time scales encompassing these different phases of solar activity. We demonstrate that the power spectrum of simulated and reconstructed solar open flux are consistent with each other. Our study provides the theoretical basis for interpreting long-term solar cycle variability based on reconstructions relying on cosmogenic isotopes and connects solar internal variations to the forcing of the state of the heliosphere. △ Less

Submitted 30 June, 2023; v1 submitted 25 August, 2022; originally announced August 2022.

Comments: 15 Pages, 8 Figures, In Press (MNRAS)

Solar Energetic Particle Ground-Level Enhancements and the Solar Cycle

Authors: Mathew Owens, Luke Barnard, Benjamin Pope, Mike Lockwood, Ilya Usoskin, Eleanna Asvestari

Abstract: Severe geomagnetic storms appear to be ordered by the solar cycle in a number of ways. They occur more frequently close to solar maximum and declining phase, are more common in larger solar cycles and show different patterns of occurrence in odd- and even-numbered solar cycles. Our knowledge of the most extreme space weather events, however, comes from the spikes in cosmogenic-isotope ($^{14}$C,… ▽ More Severe geomagnetic storms appear to be ordered by the solar cycle in a number of ways. They occur more frequently close to solar maximum and declining phase, are more common in larger solar cycles and show different patterns of occurrence in odd- and even-numbered solar cycles. Our knowledge of the most extreme space weather events, however, comes from the spikes in cosmogenic-isotope ($^{14}$C, $^{10}$Be and $^{36}$Cl) records that are attributed to significantly larger solar energetic particle (SEP) events than have been observed during the space age. Despite both storms and SEPs being driven by solar eruptive phenomena, the event-by-event correspondence between extreme storms and extreme SEPs is low. Thus it should not be assumed a priori that the solar cycle patterns found for storms also hold for SEPs and the cosmogenic-isotope events. In this study we investigate the solar cycle trends in the timing and magnitude of the 67 SEP ground-level enhancements (GLEs) recorded by neutron monitors since the mid 1950s. Using a number of models of GLE occurrence probability, we show that GLEs are around a factor four more likely around solar maximum than around solar minimum, and that they preferentially occur earlier in even-numbered solar cycles than in odd-numbered cycles. There are insufficient data to conclusively determine whether larger solar cycles produce more GLEs. Implications for putative space-weather events in the cosmogenic-isotope records are discussed. We find that GLEs tend to cluster within a few tens of days, likely due to particularly productive individual active regions, and with approximately 11-year separations, owing to the solar cycle ordering. But these timescales do not explain cosmogenic-isotope spikes which require multiple extreme SEP events over consecutive years. △ Less

Submitted 27 July, 2022; v1 submitted 26 July, 2022; originally announced July 2022.

Comments: 22 pages, 9 figures

Extreme solar events

Authors: Edward W. Cliver, Carolus J. Schrijver, Kazunari Shibata, Ilya G. Usoskin

Abstract: We trace the evolution of research on extreme solar and solar-terrestrial events from the 1859 Carrington event to the rapid development of the last twenty years. Our focus is on the largest observed/inferred/theoretical cases of sunspot groups, flares on the Sun and Sun-like stars, coronal mass ejections, solar proton events, and geomagnetic storms. The reviewed studies are based on modern observ… ▽ More We trace the evolution of research on extreme solar and solar-terrestrial events from the 1859 Carrington event to the rapid development of the last twenty years. Our focus is on the largest observed/inferred/theoretical cases of sunspot groups, flares on the Sun and Sun-like stars, coronal mass ejections, solar proton events, and geomagnetic storms. The reviewed studies are based on modern observations, historical or long-term data including the auroral and cosmogenic radionuclide record, and Kepler observations of Sun-like stars. We compile a table of 100- and 1000-year events based on occurrence frequency distributions for the space weather phenomena listed above. Questions considered include the Sun-like nature of superflare stars and the existence of impactful but unpredictable solar "black swans" and extreme "dragon king" solar phenomena that can involve different physics from that operating in events which are merely large. △ Less

Submitted 18 May, 2022; originally announced May 2022.

Comments: Invited review article for Living Reviews in Solar Physics

Journal ref: Living Rev Sol Phys 19, 2 (2022)

The First Ground Level Enhancement of Solar Cycle 25 on 28 October 2021

Authors: A. Papaioannou, A. Kouloumvakos, A. Mishev, R. Vainio, I. Usoskin, K. Herbst, A. P. Rouillard, A. Anastasiadis, J. Gieseler, R. Wimmer-Schweingruber, P. Kühl

Abstract: Aims. The first relativistic solar proton event of solar cycle 25 (SC25) was detected on 28 October 2021 by neutron monitors (NMs) on the ground and particle detectors onboard spacecraft in the near-Earth space. This is the first ground level enhancement (GLE) of the current cycle. A detailed reconstruction of the NM response together with the identification of the solar eruption that generated th… ▽ More Aims. The first relativistic solar proton event of solar cycle 25 (SC25) was detected on 28 October 2021 by neutron monitors (NMs) on the ground and particle detectors onboard spacecraft in the near-Earth space. This is the first ground level enhancement (GLE) of the current cycle. A detailed reconstruction of the NM response together with the identification of the solar eruption that generated these particles is investigated based on in-situ and remote-sensing measurements. Methods. In-situ proton observations from a few MeV to $\sim$500 MeV were combined with the detection of a solar flare in soft X-rays (SXRs), a coronal mass ejection (CME), radio bursts and extreme ultraviolet (EUV) observations to identify the solar origin of the GLE. Timing analysis was performed and a relation to the solar sources was outlined. Results. GLE73 reached a maximum particle rigidity of $\sim$2.4 GV and is associated with type III, type II, type IV radio bursts and an EUV wave. A diversity of time profiles recorded by NMs was observed. This points to an anisotropic nature of the event. The peak flux at E$>$10 MeV was only $\sim$30 pfu and remained at this level for several days. The release time of $\geq$1 GV particles was found to be $\sim$15:40 UT. GLE73 had a moderately hard rigidity spectrum at very high energies ($γ$ $\sim$5.5). Comparison of GLE73 to previous GLEs with similar solar drivers is performed △ Less

Submitted 16 February, 2022; originally announced February 2022.

Comments: accepted for publication in A&A Let

Journal ref: A&A 660, L5 (2022)

High-altitude polar NM with the new DAQ system as a tool to study details of the cosmic-ray induced nucleonic cascade

Authors: Markus Similä, Ilya Usoskin, Stepan Poluianov, Alexander Mishev, Gennady A. Kovaltsov, Du Toit Strauss

Abstract: 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 digitize… ▽ More 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. △ Less

Submitted 10 April, 2021; originally announced April 2021.

Comments: Accepted to JGR Space Physics (doi: 10.1029/2020JA028959)

Modelling the evolution of the Sun's open and total magnetic flux

Authors: N. A. Krivova, S. K. Solanki, B. Hofer, C. -J. Wu, I. G. Usoskin, R. Cameron

Abstract: Solar activity in all its varied manifestations is driven by the magnetic field. Particularly important for many purposes are two global quantities, the Sun's total and open magnetic flux, which can be computed from sunspot number records using models. Such sunspot-driven models, however, do not take into account the presence of magnetic flux during grand minima, such as the Maunder minimum. Here… ▽ More Solar activity in all its varied manifestations is driven by the magnetic field. Particularly important for many purposes are two global quantities, the Sun's total and open magnetic flux, which can be computed from sunspot number records using models. Such sunspot-driven models, however, do not take into account the presence of magnetic flux during grand minima, such as the Maunder minimum. Here we present a major update of a widely used simple model, which now takes into account the observation that the distribution of all magnetic features on the Sun follows a single power law. The exponent of the power law changes over the solar cycle. This allows for the emergence of small-scale magnetic flux even when no sunspots are present for multiple decades and leads to non-zero total and open magnetic flux also in the deepest grand minima, such as the Maunder minimum, thus overcoming a major shortcoming of the earlier models. The results of the updated model compare well with the available observations and reconstructions of the solar total and open magnetic flux. This opens up the possibility of improved reconstructions of sunspot number from time series of cosmogenic isotope production rate. △ Less

Submitted 29 March, 2021; originally announced March 2021.

Journal ref: A&A 650, A70 (2021)

Solar cyclic activity over the last millennium reconstructed from annual 14C data

Authors: I. G. Usoskin, S. K. Solanki, N. Krivova, B. Hofer, G. A. Kovaltsov, L. Wacker, N. Brehm, B. Kromer

Abstract: The 11-year solar cycle is the dominant pattern of solar activity reflecting the oscillatory dynamo mechanism in the Sun. Solar cycles were directly observed since 1700, while indirect proxies suggest their existence over a much longer period of time but generally without resolving individual cycles and their continuity. Here we reconstruct individual cycles for the last millennium using recent 14… ▽ More The 11-year solar cycle is the dominant pattern of solar activity reflecting the oscillatory dynamo mechanism in the Sun. Solar cycles were directly observed since 1700, while indirect proxies suggest their existence over a much longer period of time but generally without resolving individual cycles and their continuity. Here we reconstruct individual cycles for the last millennium using recent 14C data and state-of-the-art models. Starting with the 14C production rate determined from the so far most precise measurements of radiocarbon content in tree rings, solar activity is reconstructed in three physics-based steps: (1) Correction of the 14C production rate for the changing geomagnetic field; (2) Computation of the open solar magnetic flux; and (3) Conversion into sunspot numbers outside of grand minima. Solar activity is reconstructed for the period 971-1900 (85 individual cycles). This more than doubles the number of solar cycles known from direct solar observations. We found that lengths and strengths of well-defined cycles outside grand minima are consistent with those obtained from the direct sunspot observations after 1750. The validity of the Waldmeier rule is confirmed at a highly significant level. Solar activity is found to be in a deep grand minimum when the activity is mostly below the sunspot formation threshold, during about 250 years. Therefore, although considerable cyclic variability in 14C is seen even during grand minima, individual solar cycles can hardly be reliably resolved therein. Three potential solar particle events, ca. 994, 1052 and 1279 AD, are shown. A new about 1000-year long solar activity reconstruction, in the form of annual (pseudo) sunspot numbers with full assessment of uncertainties, is presented based on new high-precision 14C measurements and state-of-the-art models, more than doubling the number of individually resolved solar cycles. △ Less

Submitted 28 March, 2021; originally announced March 2021.

Comments: Accepted for publication in Astronomy and Astrophysics

New reconstruction of event-integrated spectra (spectral fluences) for major solar energetic particle events

Authors: Sergey A. Koldobskiy, Osku Raukunen, Rami Vainio, Gennady A. Kovaltsov, Ilya G. Usoskin

Abstract: 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… ▽ More 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. △ Less

Submitted 25 January, 2021; originally announced January 2021.

Comments: 19 pages, 3 figures, to be published in Astronomy and Astrophysics

Journal ref: A&A 647, A132 (2021)

Robustness of Solar-Cycle Empirical Rules Across Different Series Including an Updated ADF Sunspot Group Series

Authors: Ilya Usoskin, Gennady Kovaltsov, Wilma Kiviaho

Abstract: 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 "classica… ▽ More 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. △ Less

Submitted 15 December, 2020; originally announced December 2020.

Comments: Solar Physics, accepted

A New Full 3-D Model of Cosmogenic Tritium $^3$H Production in the Atmosphere (CRAC:3H)

Authors: S. V. Poluianov, G. A. Kovaltsov, I. G. Usoskin

Abstract: A new model of cosmogenic tritium ($^3$H) production in the atmosphere is presented. The model belongs to the CRAC (Cosmic-Ray Atmospheric Cascade) family and is named as CRAC:3H. It is based on a full Monte-Carlo simulation of the cosmic-ray induced atmospheric cascade using the Geant4 toolkit. The CRAC:3H model is able, for the first time, to compute tritium production at any location and time,… ▽ More A new model of cosmogenic tritium ($^3$H) production in the atmosphere is presented. The model belongs to the CRAC (Cosmic-Ray Atmospheric Cascade) family and is named as CRAC:3H. It is based on a full Monte-Carlo simulation of the cosmic-ray induced atmospheric cascade using the Geant4 toolkit. The CRAC:3H model is able, for the first time, to compute tritium production at any location and time, for any given energy spectrum of the primary incident cosmic ray particles, explicitly treating, also for the first time, particles heavier than protons. This model provides a useful tool for the use of $^3$H as a tracer of atmospheric and hydrological circulation. A numerical recipe for practical use of the model is appended. △ Less

Submitted 16 September, 2020; originally announced September 2020.

Journal ref: J. Geophys. Res. Atmos., 125, e2020JD033147 (2020)

Revisited reference solar proton event of 23-Feb-1956: Assessment of the cosmogenic-isotope method sensitivity to extreme solar events

Authors: Ilya G. Usoskin, Sergey A. Koldobskiy, Gennady A. Kovaltsov, Eugene V. Rozanov, Timophei V. Sukhodolov, Alexander L. Mishev, Irina A. Mironova

Abstract: 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 id… ▽ More 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. △ Less

Submitted 21 May, 2020; originally announced May 2020.

Journal ref: Journal of Geophysical Research: Space Physics, 125, e2020JA027921 (2020)

Spectral characteristic of mid-term quasi-periodicities in sunspots data

Authors: P. Frick, D. Sokoloff, R. Stepanov, V. Pipin, I. Usoskin

Abstract: 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 prono… ▽ More 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. △ Less

Submitted 15 November, 2019; originally announced November 2019.

Comments: accepted in MNRAS

Validation of the neutron monitor yield function using data from AMS-02 experiment, 2011--2017

Authors: Sergey A. Koldobskiy, Veronica Bindi, Claudio Corti, Gennady A. Kovaltsov, Ilya G. Usoskin

Abstract: 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 verif… ▽ More 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. △ Less

Submitted 3 April, 2019; originally announced April 2019.

Comments: J. Geophys. Res. in press

Journal ref: J. Geophys. Res. (Space Phys.), 124, 2367-2379, 2019

Historical astronomical data: urgent need for preservation, digitization enabling scientific exploration

Authors: Alexei Pevtsov, Elizabeth Griffin, Jonathan Grindlay, Stella Kafka, Jennifer Lynn Bartlett, Ilya Usoskin, Kalevi Mursula, Sarah Gibson, Valentin M. Pillet, Joan Burkepile, David Webb, Frederic Clette, James Hesser, Peter Stetson, Andres Munoz-Jaramillo, Frank Hill, Rick Bogart, Wayne Osborn, Dana Longcope

Abstract: 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 degrad… ▽ More 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. △ Less

Submitted 12 March, 2019; originally announced March 2019.

Comments: 8 pages, White Paper submitted to ASTRO2020 Decadal Survey

Journal ref: White Paper, ASTRO2020 Decadal Survey, 2020

Solar total and spectral irradiance reconstruction over the last 9000 years

Authors: C. -J. Wu, N. A. Krivova, S. K. Solanki, I. G. Usoskin

Abstract: 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 irradian… ▽ More 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. △ Less

Submitted 8 November, 2018; originally announced November 2018.

First analysis of GLE 72 event on 10 September 2017: Spectral and anisotropy characteristics

Authors: A. L. Mishev, I. G. Usoskin, O. Raukunen, M. Paassilta, E. Valtonen, L. G. Kocharov, R. Vainio

Abstract: 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 pro… ▽ More 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. △ Less

Submitted 24 October, 2018; originally announced October 2018.

Comments: 21 pages, 10 figures, 1 table

Journal ref: Solar Physics 293, 136, 2018

Solar energetic particles and galactic cosmic rays over millions of years as inferred from data on cosmogenic $^{26}$Al in lunar samples

Authors: S. Poluianov, G. A. Kovaltsov, I. G. Usoskin

Abstract: 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… ▽ More 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 ($φ=496\pm 40$ MV over a million years versus $φ= 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. △ Less

Submitted 26 July, 2018; originally announced July 2018.

Comments: Astron. Astrophys., in press

Solar activity over nine millennia: A consistent multi-proxy reconstruction

Authors: Chi Ju Wu, I. G. Usoskin, N. Krivova, G. A. Kovaltsov, M. Baroni, E. Bard, S. K. Solanki

Abstract: 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 te… ▽ More 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. △ Less

Submitted 4 April, 2018; originally announced April 2018.

Comments: A&A, in print (the abstract has been slightly modified w/r to the original paper to match the ArXiv length limitations)

Journal ref: A&A 615, A93 (2018)

Revised historical solar irradiance forcing

Authors: T. Egorova, W. Schmutz, E. Rozanov, A. I. Shapiro, I. Usoskin, J. Beer, R. V. Tagirov, T. Peter

Abstract: 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 reconstructio… ▽ More 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 %. △ Less

Submitted 1 April, 2018; originally announced April 2018.

Comments: Accepted for publication in Astronomy & Astrophysics, 10 pages, 10 figures

Journal ref: A&A 615, A85 (2018)

A Test of the Active Day Fraction Method of Sunspot Group Number Calibration: Dependence on the Level of Solar Activity

Authors: Teemu Willamo, Ilya G. Usoskin, Gennady A. Kovaltsov

Abstract: 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 t… ▽ More 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. △ Less

Submitted 28 March, 2018; originally announced March 2018.

Comments: Solar Physics, accepted

Extreme value theory applied to the millennial sunspot number series

Authors: F. J. Acero, M. C. Gallego, J. A. García, I. G. Usoskin, J. M. Vaquero

Abstract: 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 differ… ▽ More 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. △ Less

Submitted 29 January, 2018; originally announced January 2018.

Comments: 15 pages, 5 figures, 5 tables, accepted for publication in ApJ

GLE and Sub-GLE Redefinition in the Light of High-Altitude Polar Neutron Monitors

Authors: S. V. Poluianov, I. G. Usoskin, A. L. Mishev, M. A. Shea, D. F. Smart

Abstract: 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-m… ▽ More 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. △ Less

Submitted 28 November, 2017; v1 submitted 15 November, 2017; originally announced November 2017.

Journal ref: Poluianov, S., Usoskin, I., Mishev, A., Shea, M., Smart, D., GLE and Sub-GLE Redefinition in the Light of High-Altitude Polar Neutron Monitors, Solar Phys., 292 (11), 2017, doi: 10.1007/s11207-017-1202-4

Comment on the paper by Popova et al. "On a role of quadruple component of magnetic field in defining solar activity in grand cycles"

Authors: Ilya G. Usoskin

Abstract: 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 da… ▽ More 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. △ Less

Submitted 14 October, 2017; originally announced October 2017.

Comments: Accepted to J.Atmos. Solar-Terr. Phys

Can superflares occur on the Sun? A view from dynamo theory

Authors: M. M. Katsova, L. L. Kitchatinov, M. A. Livshits, D. L. Moss, D. D. Sokoloff, I. G. Usoskin

Abstract: 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-mech… ▽ More 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. △ Less

Submitted 29 September, 2017; originally announced October 2017.

Comments: 13 pages, 3 figures and 1 table

Updated sunspot group number reconstruction for 1749-1996 using the active day fraction method

Authors: Teemu Willamo, Ilya G. Usoskin, Gennady A. Kovaltsov

Abstract: 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, an… ▽ More 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. △ Less

Submitted 15 May, 2017; originally announced May 2017.

Comments: Published in Astron. Astrophys

Journal ref: Astron, Astrophys., 601, A109, 2017

New reconstruction of the sunspot group number since 1739 using the direct calibration and "backbone" methods

Authors: Theodosios Chatzistergos, Ilya G. Usoskin, Gennady A. Kovaltsov, Natalie A. Krivova, Sami K. Solanki

Abstract: 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 ser… ▽ More 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. △ Less

Submitted 27 February, 2017; v1 submitted 20 February, 2017; originally announced February 2017.

Comments: 19 pages, 12 figures, Astronomy and Astrophysics, accepted

Journal ref: A&A 602, A69 (2017)

Application of a full chain analysis using neutron monitor data for space weather studies

Authors: A. Mishev, I. Usoskin

Abstract: 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 charac… ▽ More 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. △ Less

Submitted 29 December, 2016; v1 submitted 21 December, 2016; originally announced December 2016.

Comments: "XXV ECRS 2016 Proceedings - eConf C16-09-04.3", highlight talk

Analysis of Ground Level Enhancements (GLE): Extreme solar energetic particle events have hard spectra

Authors: E. Asvestari, T. Willamo, A. Gil, I. G. Usoskin, G. A. Kovaltsov, V. V. Mikhailov, A. Mayorov

Abstract: 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 exp… ▽ More 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. △ Less

Submitted 18 September, 2016; originally announced December 2016.

Comments: Advances in Space Research (in press)

Neutron Monitor Yield Function: New Improved computations

Authors: A. L. Mishev, I. G. Usoskin, G. A. Kovaltsov

Abstract: 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 compu… ▽ More 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 $α-$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. △ Less

Submitted 25 October, 2016; originally announced December 2016.

Comments: Includes also Erratum (submitted to JGR) to Table 1

Journal ref: Mishev, A., I.G. Usoskin, G.A. Kovaltsov, Neutron monitor yield function: New improved computations, J. Geophys. Res., 118, 783-2788, 2013

An Optical Atmospheric Phenomenon Observed in 1670 over the City of Astrakhan Was not a Mid-Latitude Aurora

Authors: I. G. Usoskin, G. A. Kovaltsov, L. N. Mishina, D. D. Sokoloff, J. Vaquero

Abstract: 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, a… ▽ More 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. △ Less

Submitted 28 November, 2016; originally announced December 2016.

Comments: accepted to Solar Physics

Dependence of the Sunspot-group Size on the Level of Solar Activity and its Influence on the Calibration of Solar Observers

Authors: I. G. Usoskin, G. A. Kovaltsov, T. Chatzistergos

Abstract: 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… ▽ More 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. △ Less

Submitted 2 September, 2016; originally announced September 2016.

Comments: Solar Physics, accepted

Production of cosmogenic isotopes 7Be, 10Be, 14C, 22Na and 36Cl in the atmosphere: Altitudinal profiles of yield functions

Authors: Stepan Poluianov, Gennady A. Kovaltsov, Alexander L. Mishev, Ilya G. Usoskin

Abstract: 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 atmosphe… ▽ More 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. △ Less

Submitted 19 June, 2016; originally announced June 2016.

Comments: in press (J. Geophys. Res. (Atmosphere), doi: 10.1002/2016JD025034)

Tests of Sunspot Number Sequences: 2. Using Geomagnetic and Auroral Data

Authors: Mike Lockwood, Mathew J. Owens, Luke A Barnard, Chris J. Scott, Ilya G. Usoskin, Heikki Nevanlinna

Abstract: 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 t… ▽ More 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. △ Less

Submitted 6 May, 2016; originally announced May 2016.

Comments: Solar Physics, in Press. 26 pages, 8 figures

An Assessment of Sunspot Number Data Composites over 1845-2014

Authors: Mike Lockwood, Mathew J. Owens, Luke A. Barnard, Ilya G. Usoskin

Abstract: 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 in… ▽ More 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. △ Less

Submitted 15 April, 2016; originally announced April 2016.

Comments: In press, Astrophys. J., 2016

Analysis of the ground level enhancements on 14 July 2000 and on 13 December 2006 using neutron monitor data

Authors: Alexander Mishev, Ilya Usoskin

Abstract: 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 asymp… ▽ More 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. △ Less

Submitted 29 March, 2016; originally announced March 2016.

Comments: 23 pages, 9 figures, 3 tables; Accepted for publication in Solar Physics

Solar activity during the Holocene: the Hallstatt cycle and its consequence for grand minima and maxim

Authors: I. G. Usoskin, Y. Gallet, F. Lopes, G. A. Kovaltsov, G. Hulot

Abstract: 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 mil… ▽ More 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. △ Less

Submitted 8 February, 2016; originally announced February 2016.

Comments: In press in Astronomy & Astrophysics, doi: 10.1051/0004-6361/201527295

Journal ref: Astron. Astrophys., 587, A150, 2016

Atmospheric ionization induced by precipitating electrons: Comparison of CRAC:EPII model with parametrization model

Authors: A. A. Artamonov, A. L. Mishev, I. G. Usoskin

Abstract: 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-io… ▽ More 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 △ Less

Submitted 22 January, 2016; originally announced January 2016.

Comments: 13 pages, 4 figures. Under review in JASTP

A New Calibrated Sunspot Group Series Since 1749: Statistics of Active Day Fractions

Authors: I. G. Usoskin, G. A. Kovaltsov, M. Lockwood, K. Mursula, M. Owens, S. K. Solanki

Abstract: 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 nov… ▽ More 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. △ Less

Submitted 20 December, 2015; originally announced December 2015.

Comments: Accepted for publication at Solar Physics

A two-wave dynamo model by Zharkova et al. (2015) disagrees with data on long-term solar variability

Authors: I. Usoskin, G. Kovaltsov

Abstract: 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.… ▽ More 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. △ Less

Submitted 17 December, 2015; originally announced December 2015.

Comments: 8 pages, 1 figure, an extended version of a comment sent to Sci. Rep

Tests of sunspot number sequences: 3. Effects of regression procedures on the calibration of historic sunspot data

Authors: M. Lockwood, M. J Owens, L. Barnard, I. G. Usoskin

Abstract: 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 sy… ▽ More 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. △ Less

Submitted 24 January, 2016; v1 submitted 27 October, 2015; originally announced October 2015.

Comments: 14 pages, 1 Table, 11 figures

The Maunder minimum (1645--1715) was indeed a Grand minimum: A reassessment of multiple datasets

Authors: Ilya G. Usoskin, Rainer Arlt, Eleanna Asvestari, Ed Hawkins, Maarit Käpylä, Gennady A. Kovaltsov, Natalie Krivova, Michael Lockwood, Kalevi Mursula, Jezebel O'Reilly, Matthew Owens, Chris J. Scott, Dmitry D. Sokoloff, Sami K. Solanki, Willie Soon, José M. Vaquero

Abstract: 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 durin… ▽ More 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. △ Less

Submitted 18 July, 2015; originally announced July 2015.

Comments: Accepted to Astron. Astrophys

Journal ref: A&A 581, A95 (2015)

Long-term variation in the Sun's activity caused by magnetic Rossby waves in the tachocline

Authors: T. V. Zaqarashvili, R. Oliver, A. Hanslmeier, M. Carbonell, J. L. Ballester, T. Gachechiladze, I. G. Usoskin

Abstract: 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… ▽ More 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. △ Less

Submitted 11 May, 2015; originally announced May 2015.

Comments: 4 pages, 4 figures, accepted in ApJL

Level and length of cyclic solar activity during the Maunder minimum as deduced from the active day statistics

Authors: J. M. Vaquero, G. A. Kovaltsov, I. G. Usoskin, V. M. S. Carrasco, M. C. Gallego

Abstract: 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 as… ▽ More 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. △ Less

Submitted 26 March, 2015; originally announced March 2015.

Comments: Accepted to Astron. Astrophys

Journal ref: Astron. Astrophys., 577, A71, 2015

Solar surface rotation: N-S asymmetry and recent speed-up

Authors: L. Zhang, K. Mursula, I. Usoskin

Abstract: 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 reduc… ▽ More 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. △ Less

Submitted 12 January, 2015; originally announced January 2015.

Comments: Astron. Astrophys. Lett. (accepted)

Journal ref: A&A 575, L2 (2015)

Has the Earth been exposed to numerous supernovae within the last 300 kyr?

Authors: A. L. Melott, I. G. Usoskin, G. A Kovaltsov, C. M. Laird

Abstract: 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 likel… ▽ More 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. △ Less

Submitted 25 September, 2014; v1 submitted 19 June, 2014; originally announced June 2014.

Comments: 8 pages, 1 figure. to be published in the Interenational Journal of Astrobiology

Journal ref: International Journal of Astrobiology 14, 375-378 (2015)