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

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

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)

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)

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)

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

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)

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)

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

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

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

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)

Evidence for distinct modes of solar activity

Authors: I. G. Usoskin, G. Hulot, Y. Gallet, R. Roth, A. Licht, F. Joos, G. A. Kovaltsov, E. Thebault, A. Khokhlov

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

Submitted 19 February, 2014; originally announced February 2014.

Comments: Astron. Astrophys. Letters (in press)

Journal ref: Astron. Astrophys., 562, L10, 2014

Inconsistency of the Wolf sunspot number series around 1848

Authors: R. Leussu, I. G. Usoskin, R. Arlt, K. Mursula

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

Submitted 1 November, 2013; v1 submitted 31 October, 2013; originally announced October 2013.

Comments: Published as: Leussu, R., I.G. Usoskin, R. Arlt and K. Mursula, Inconsistency of the Wolf sunspot number series around 1848, Astron. Astrophys., 559, A28, 2013

Journal ref: Astron. Astrophys., 559, A28, 2013

Sunspot positions and sizes for 1825-1867 from the observations by Samuel Heinrich Schwabe

Authors: R. Arlt, R. Leussu, N. Giese, K. Mursula, I. G. Usoskin

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

Submitted 31 May, 2013; originally announced May 2013.

Comments: 8 pages, 4 figures, Mon. Not. R. Astron. Soc., accepted

Occurrence Probability of Large Solar Energetic Particle Events: Assessment from Data on Cosmogenic Radionuclides in Lunar Rocks

Authors: Gennady A. Kovaltsov, Ilya G. Usoskin

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

Submitted 22 May, 2013; originally announced May 2013.

Comments: Solar Physics, in press

Journal ref: Solar Phys., 289, 211-220 (2014)

The AD775 cosmic event revisited: the Sun is to blame

Authors: I. G. Usoskin, B. Kromer, F. Ludlow, J. Beer, M. Friedrich, G. A. Kovaltsov, S. K. Solanki, L. Wacker

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

Submitted 22 May, 2013; v1 submitted 27 February, 2013; originally announced February 2013.

Comments: Published by Astron. Astrophys. Lett

Journal ref: Astron. Astrophys. Lett., 552, L3, 2013

The First Ground Level Enhancement Event of Solar Cycle 24: Direct Observation of Shock Formation and Particle Release Heights

Authors: N. Gopalswamy, H. Xie, S. Akiyama, S. Yashiro, I. G. Usoskin, J. M. Davila

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

Submitted 6 February, 2013; originally announced February 2013.

Comments: 16 pages, 4 figures, 1 table

Waldmeier relations and the solar cycle dynamics by the mean-field dynamos

Authors: V. V. Pipin, D. D. Sokoloff, I. G. Usoskin

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

Submitted 11 November, 2012; originally announced November 2012.

Comments: 2 pages, IAU 294 Proceedings

Occurrence of extreme solar particle events: Assessment from historical proxy data

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

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

Submitted 25 July, 2012; originally announced July 2012.

Comments: accepted to Astrophys. J

Journal ref: Astrophys. J., 757, 92, 2012

A new model of cosmogenic production of radiocarbon 14C in the atmosphere

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

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

Submitted 29 June, 2012; originally announced June 2012.

Comments: Published in EPSL, 337, 114, 2012

Journal ref: Earth Planet. Sci. Lett., 337, 114-120, 2012

Estimating the frequency of extremely energetic solar events, based on solar, stellar, lunar, and terrestrial records

Authors: C. J. Schrijver, J. Beer, U. Baltensperger, E. W. Cliver, M. Guedel, H. S. Hudson, K. G. McCracken, R. A. Osten, Th. Peter, D. R. Soderblom, I. G. Usoskin, E. W. Wolff

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

Submitted 21 June, 2012; originally announced June 2012.

Comments: 14 pages, 3 figures (in press as of 2012/06/18); Journal of Geophysical Research (Space Physics), 2012

Properties of Ground Level Enhancement Events and the Associated Solar Eruptions during Solar Cycle 23

Authors: N. Gopalswamy, H. Xie, S. Yashiro, S. Akiyama, P. Mäkelä, I. G. Usoskin

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

Submitted 3 May, 2012; originally announced May 2012.

Comments: 17 figures, 56 pages, based on two LWS/CDAWs on GLE events, part of Topical Issue on GLE events in Space Science Reviews

Variations of the solar cycle profile in a solar dynamo with fluctuating dynamo governing parameters

Authors: V. V. Pipin, D. D. Sokoloff, I. G. Usoskin

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

Submitted 23 May, 2012; v1 submitted 28 December, 2011; originally announced December 2011.

Comments: 12 pages, 9 figures, to appear in A&A

Revisited sunspot data: A new scenario for the onset of the Maunder minimum

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

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

Submitted 8 March, 2011; originally announced March 2011.

Comments: 8 pages, 2 figures, accepted for publication in ApJ Letters

A solar cycle lost in 1793--1800: Early sunspot observations resolve the old mystery

Authors: Ilya G. Usoskin, Kalevi Mursula, Rainer Arlt, Gennady A. Kovaltsov

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

Submitted 19 July, 2009; v1 submitted 1 July, 2009; originally announced July 2009.

Comments: Published by Astrophys. J. Lett

Journal ref: Astrophys.J.700:L154-L157,2009

A History of Solar Activity over Millennia

Authors: Ilya G. Usoskin

Abstract: Presented here is a review of present knowledge of the long-term behavior of solar activity on a multi-millennial timescale, as reconstructed using the indirect proxy method. The concept of solar activity is discussed along with an overview of the special indices used to quantify different aspects of variable solar activity, with special emphasis upon sunspot number. Over long timescales, quantita… ▽ More Presented here is a review of present knowledge of the long-term behavior of solar activity on a multi-millennial timescale, as reconstructed using the indirect proxy method. The concept of solar activity is discussed along with an overview of the special indices used to quantify different aspects of variable solar activity, with special emphasis upon sunspot number. Over long timescales, quantitative information about past solar activity can only be obtained using a method based upon indirect proxies, such as the cosmogenic isotopes \super{14}C and \super{10}Be in natural stratified archives (e.g., tree rings or ice cores). We give an historical overview of the development of the proxy-based method for past solar-activity reconstruction over millennia, as well as a description of the modern state. Special attention is paid to the verification and cross-calibration of reconstructions. It is argued that this method of cosmogenic isotopes makes a solid basis for studies of solar variability in the past on a long timescale (centuries to millennia) during the Holocene. A separate section is devoted to reconstructions of strong solar energetic-particle (SEP) events in the past, that suggest that the present-day average SEP flux is broadly consistent with estimates on longer timescales, and that the occurrence of extra-strong events is unlikely. Finally, the main features of the long-term evolution of solar magnetic activity, including the statistics of grand minima and maxima occurrence, are summarized and their possible implications, especially for solar/stellar dynamo theory, are discussed. △ Less

Submitted 12 April, 2013; v1 submitted 22 October, 2008; originally announced October 2008.

Comments: A review, 91 pages, 28 figures. available online at http://solarphysics.livingreviews.org/Articles/lrsp-2013-1/

Journal ref: Usoskin, I.G., A History of Solar Activity over Millennia, Living. Rev. Solar Phys., 10, 1, 2013