Does Machine Learning reconstruct missing sunspots and forecast a new solar minimum?

Abstract The retrodiction and prediction of solar activity are two closely-related problems in dynamo theory. We applied Machine Learning (ML) algorithms and analyses to the World Data Center’s newly constructed annual sunspot time series (1700-2019; Version 2.0). This provides a unique model that gives insights into the various patterns of the Sun’s magnetic dynamo that drives solar activity maxima and minima. We found that the variability in the ∼ 11 -year Sunspot Cycle is closely connected with 120-year oscillatory magnetic activity variations. We also identified a previously under-eported 5.5 year periodicity in the sunspot record. This 5.5-year pattern is co-modulated by the 120-year oscillation and appears to influence the shape and energy/power content of individual 11-year cycles. Our ML algorithm was trained to recognize such underlying patterns and provides a convincing hindcast of the full sunspot record from 1700 to 2019. It also suggests the possibility of missing sunspots during Sunspot Cycles -1, 0, and 1 (ca. 1730s-1760s). In addition, our ML model forecasts a new phase of extended solar minima that began prior to Sunspot Cycle 24 (ca. 2008-2019) and will persist until Sunspot Cycle 27 (ca. 2050 or so). Our ML Bayesian model forecasts a peak annual sunspot number (SSN) of 95 with a probable range of 80 to 115 for Cycle 25 between 2023 and 2025.

[1]  H. Yoshimura The solar-cycle period-amplitude relation as evidence of hysteresis of the solar-cycle nonlinear magnetic oscillation and the long-term (55 year) cyclic modulation , 1979 .

[2]  Z. Kolláth,et al.  Multiple and changing cycles of active stars - I. Methods of analysis and application to the solar cycles , 2009, 0904.1725.

[3]  J. V. Lincoln,et al.  Prediction of sunspot numbers , 1949 .

[4]  Stephen A. Billings,et al.  Modeling and prediction of global magnetic disturbance in near‐Earth space: A case study for Kp index using NARX models , 2016 .

[5]  L. Lefèvre,et al.  Uncertainties in the Sunspot Numbers: Estimation and Implications , 2016, 1608.05261.

[6]  W. Gleissberg Probability Laws of Sunspot Variations. , 1942 .

[7]  K. M. Hiremath Prediction of solar cycle 24 and beyond , 2008 .

[8]  D. Sokoloff,et al.  The Maunder minimum: a mixed-parity dynamo mode? , 1994 .

[9]  Aslak Grinsted,et al.  Nonlinear Processes in Geophysics Application of the Cross Wavelet Transform and Wavelet Coherence to Geophysical Time Series , 2022 .

[10]  J. M. Mitchell,et al.  On the Power Spectrum of “Red Noise” , 1963 .

[11]  I. Usoskin,et al.  On the Claimed 5.5-year Periodicity in Solar Activity , 1997 .

[12]  S. Solanki,et al.  Analyzing Solar Cycles , 2011, Science.

[13]  W. Soon,et al.  Covariations of chromospheric and photometric variability of the young Sun analogue HD 30495: evidence for and interpretation of mid-term periodicities , 2018, Monthly Notices of the Royal Astronomical Society.

[14]  Joan Feynman,et al.  Period and phase of the 88-year solar cycle and the Maunder minimum: Evidence for a chaotic sun , 1990 .

[15]  F. Lopes,et al.  Characteristic Time Scales of Decadal to Centennial Changes in Global Surface Temperatures Over the Past 150 Years , 2020, Earth and Space Science.

[16]  C. Torrence,et al.  A Practical Guide to Wavelet Analysis. , 1998 .

[17]  F. Lopes,et al.  A Solar Signature in Many Climate Indices , 2019, Journal of Geophysical Research: Atmospheres.

[18]  P. Preka-Papadema,et al.  On signal–noise decomposition of time‐series using the continuous wavelet transform: application to sunspot index , 2003 .

[19]  W. Soon,et al.  A review of Holocene solar-linked climatic variation on centennial to millennial timescales: Physical processes, interpretative frameworks and a new multiple cross-wavelet transform algorithm , 2014 .

[20]  R. Cameron,et al.  Solar activity: periodicities beyond 11 years are consistent with random forcing , 2019, Astronomy & Astrophysics.

[21]  B. Mendoza,et al.  On Mid-Term Periodicities in Sunspot Groups and Flare Index , 2011 .

[22]  D. Legates,et al.  A refined index of model performance: a rejoinder , 2013 .

[23]  N. Weiss,et al.  For how long will the current grand maximum of solar activity persist? , 2008 .

[24]  M. Kremliovsky Can we understand time scales of solar activity? , 1994 .

[25]  J. M. Vaquero Historical sunspot observations: A review , 2007 .

[26]  S. Silverman,et al.  Power spectral analysis of auroral occurrence frequency , 1983 .

[27]  Judit Pap,et al.  Periodicities of solar irradiance and solar activity indices, I , 1990 .

[28]  A. Grossmann,et al.  DECOMPOSITION OF HARDY FUNCTIONS INTO SQUARE INTEGRABLE WAVELETS OF CONSTANT SHAPE , 1984 .

[29]  W. Pesnell Effects of Version 2 of the International Sunspot Number on Naïve Predictions of Solar Cycle 25 , 2018, Space Weather.

[30]  D. Nandy,et al.  Prediction of the strength and timing of sunspot cycle 25 reveal decadal-scale space environmental conditions , 2018, Nature Communications.

[31]  Colin Jacobs,et al.  Surveying the reach and maturity of machine learning and artificial intelligence in astronomy , 2019, WIREs Data Mining Knowl. Discov..

[32]  V. Letfus Relative sunspot numbers in the first half of eighteenth century , 2000 .

[33]  F. Steinhilber,et al.  Prediction of solar activity for the next 500 years , 2013 .

[34]  W. Gleissberg A long-periodic fluctuation of the sun-spot numbers , 1939 .

[35]  W. Feng,et al.  WHY IS THE SOLAR CONSTANT NOT A CONSTANT? , 2012, 1201.0819.

[36]  S. D. Bouwer,et al.  Periodicities of solar irradiance and solar activity indices, II , 1990 .

[37]  Paul E. Damon,et al.  Persistence of the Gleissberg 88‐year solar cycle over the last ∼12,000 years: Evidence from cosmogenic isotopes , 2003 .

[38]  A. Veronig,et al.  Sunspot Number Second Differences as a Precursor of the Following 11-year Sunspot Cycle , 2017, 1712.05782.

[39]  Xenophon Moussas,et al.  A nonlinear RLC solar cycle model , 1996 .

[40]  S. Ananthakrishnan,et al.  Another Mini Solar Maximum in the Offing: A Prediction for the Amplitude of Solar Cycle 25 , 2019, Journal of Geophysical Research: Space Physics.

[41]  Edward W. Cliver,et al.  Revisiting the Sunspot Number , 2014, 1407.3231.

[42]  D. Legates,et al.  Evaluating the use of “goodness‐of‐fit” Measures in hydrologic and hydroclimatic model validation , 1999 .

[43]  M. Kremliovsky Limits of predictability of solar activity , 1995 .

[44]  J. Xie,et al.  Phase Relations between the Sunspot Numbers and Total Solar Irradiance , 2017 .

[45]  Enrico Camporeale,et al.  The Challenge of Machine Learning in Space Weather: Nowcasting and Forecasting , 2019, Space Weather.

[46]  P. Paquet,et al.  The common oscillations of solar activity, the geomagnetic field, and the Earth's rotation , 1996 .

[47]  A. Andersen,et al.  The analysis and forecasting of the Wolf sunspot numbers , 1980 .

[48]  N. D. Lewis,et al.  Deep Learning Made Easy with R: A Gentle Introduction For Data Science , 2016 .

[49]  N. Watkins,et al.  Overlapping Magnetic Activity Cycles and the Sunspot Number: Forecasting Sunspot Cycle 25 Amplitude , 2020, Solar Physics.

[50]  T. Weier,et al.  Phase coherence and phase jumps in the Schwabe cycle , 2020, Astronomische Nachrichten.

[51]  Heinrich Schwabe,et al.  Sonnen — Beobachtungen im Jahre 1843 , 1844 .

[52]  P. Charbonneau Dynamo models of the solar cycle , 2020, Living Reviews in Solar Physics.

[53]  R. Shah,et al.  Least Squares Support Vector Machines , 2022 .

[54]  T. Bayes LII. An essay towards solving a problem in the doctrine of chances. By the late Rev. Mr. Bayes, F. R. S. communicated by Mr. Price, in a letter to John Canton, A. M. F. R. S , 1763, Philosophical Transactions of the Royal Society of London.

[55]  S. Solanki,et al.  Solar activity reconstructed over the last 7000 years: The influence of geomagnetic field changes , 2006 .

[56]  Blanca Mendoza,et al.  Reconstruction and prediction of the total solar irradiance: From the Medieval Warm Period to the 21st century , 2015 .

[57]  Douglas V. Hoyt,et al.  Group Sunspot Numbers: A New Solar Activity Reconstruction , 1998 .

[58]  T. Weier,et al.  Schwabe, Gleissberg, Suess-de Vries: Towards a consistent model of planetary synchronization of solar cycles , 2019, 1910.10383.

[59]  D. Callebaut,et al.  Large-Scale Magnetic Field and Sunspot Cycles , 2001 .

[60]  R. P. Kane Some Implications Using the Group Sunspot Number Reconstruction , 2002 .

[61]  D. Hathaway The Solar Cycle , 2010, Living reviews in solar physics.

[62]  S. Sello,et al.  Long-Term Solar Variability and the Solar Cycle in the 21st Century , 2004 .

[63]  J Kugler [Power spectrum]. , 1978, Deutsche medizinische Wochenschrift.