Essential features of long-term changes of areas and diameters of sunspot groups in solar activity cycles 12–24

We analyze the Greenwich catalog data on areas of sunspot groups of last thirteen solar cycles. Various parameters of sunspots are considered, namely: average monthly smoothed areas, maximum area for each year and equivalent diameters of groups of sunspots. The first parameter shows an exceptional power of the 19th cycle of solar activity, which appears here more contrastively than in the numbers of spots (that is, in Wolf’s numbers). It was found that in the maximum areas of sunspot groups for a year there is a unique phenomenon: a short and high jump in the 18th cycle (in 1946–1947) that has no analogues in other cycles. We also studied the integral distributions for equivalent diameters and found the following: (a) the average value of the index of power-law approximation is 5.4 for the last 13 cycles and (b) there is reliable evidence of Hale's double cycle (about 44?years). Since this indicator reflects the dispersion of sunspot group diameters, the results obtained show that the convective zone of the Sun generates embryos of active regions in different statistical regimes which change with a cycle of about 44?years.     

Specific features of the high-speed plasma stream cycles

The high-speed plasma streams in the solar wind are investigated during the solar cycles nos. 20–22 (1964–1996), separately on the two types of streams according to their solar origin: the HSPS produced by coronal holes (co-rotating) and the flare-generated, in keeping with the classification made in different catalogues. The analysis is performed taking into account the following high-speed stream parameters: the durations (in days), the maximum velocities, the velocity gradients and, the importance of the streams. The time variation of these parameters and the high-speed plasma streams occurrence rate show an 11-year periodicity with some differences between the solar cycles considered. A detailed analysis of the high-speed stream 11-year cycles is made by comparison with the “standard” cycles of the sunspot relative number (Wolf number). The different behaviour of the high-speed stream parameters between even and odd solar cycles could be due to the 22-year solar magnetic cycle. The increased activity of the high-speed plasma streams on the descendant phases of the cycles, regardless of their solar sources, proves the existence of some special local conditions of the solar plasma and the magnetic field on a large scale that allow the ejection of the high energy plasma streams. This fact has led us to the analysis the stream parameters during the different phases of the solar cycles (minimum, ascendant, maximum and, descendant) as well as during the polar magnetic field reversal intervals. The differences between the phases considered are pointed out. The solar cycles 20 and 22 reveal very similar dynamics of the flare-generated and also co-rotating stream parameters during the maximum, descendant and reversal intervals. This fact could be due to their position in a Hale Cycle (the first component of the 22-year solar magnetic cycle). The 21st solar cycle dominance of all co-rotating stream parameters against the 20th and 22nd solar cycle ones, during almost all phases, could be due to the same structure of a Hale Cycle – solar cycle 21 is the second component in a 22-year SC. During the reversal intervals, all high-speed stream parameters have comparable values with the ones of the maximum phases of the cycles even if this interval contains a small part of the descendant branch (solar cycles 20 and 22).     

Solar terrestrial effects of two distinct types

The occurrence frequencies or fluxes of most of the solar phenomena show a 11-year cycle like that of sunspots. However, the average characteristics of these phenomena may not show a 11-year cycle. Among the terrestrial parameters, some related directly to the occurrence frequencies of solar phenomena (for example, ionospheric number densities related to solar EUV fluxes which show 11-year cycle like sunspots) show 11-year cycles, including the double-peak structures near sunspot maxima. Other terrestrial parameters related to average characteristics may not show 11-year sunspot cycles. For example, long-term geomagnetic activity (Ap or Dst indices) is related to the average interplanetary solar wind speed V and the total magnetic field B. The average values of V depend not on the occurrence frequency of ICMEs and/or CIRs as such, but on the relative proportion of slow and high-speed events in them. Hence, V values (and Ap values) in any year could be low, normal or high irrespective of the phase of the 11-year cycle, except that during sunspot minimum, V (and Ap) values are also low. However, 2–3 years after the solar minimum (well before sunspot maximum), V values increase, oscillate near a high level for several years, and may even increase further during the declining phase of sunspot activity, due to increased influence of high-speed CIRs (corotating interplanetary regions). Thus, Ap would have no fixed relationship with sunspot activity. If some terrestrial parameter shows a 11-year cycle, chances are that the solar connection is through the occurrence frequencies (and not average characteristics) of some solar parameter.     

利用EMD方法提取太阳活动周期成分

EMD(经验模态分解)方法在处理非线性及非平稳时间序列时表现出了很大的优势和应用潜力.利用EMD方法研究太阳活动周期,对110年(1894-2003)和55年(1949-2003)的太阳黑子数月均值进行分解,分别得到一系列模式和一个趋势项,其中都可能包含有1.3至1.4年周期分量,25至30个月QBO(准双年振荡)分量,11年太阳周分量和22年Hale周分量.其中11年周期分量幅度最大,变化特征与太阳黑子数原始数据具有很高的相似性.不同于传统方法,EMD方法给出了太阳活动在不同时间尺度上各自分离的变化特征.      

1996-2002年太阳耀斑的统计分析

分析了1996-2002年南北半球的太阳黑子相对数和南北半球太阳X射线耀斑级别(简称Imp)≥M1．0的太阳X射线耀斑的特征和不对称性．分析结果表明，南北半球的太阳耀斑活动的程度交替上升，在2001年7月以前北半球的太阳耀斑活动强于南半球，2001年7月开始耀斑活动逐渐以南半球为主．本文还逐月分析了1996—2001年南北半球的耀斑指数．2000年7月为第23周太阳指数最大的一个月，与第23周太阳黑子相对数最大月均值吻合．     

中等磁暴太阳周分布的统计研究

依据实际观测的中等磁暴数据,统计分析了中等磁暴的太阳周分布.分析结果表明,在一个太阳活动周内,每年中等磁暴随时间的变化出现多个峰值,其中,最大峰值均出现在太阳活动周的下降段,即中等磁暴的峰值比太阳黑子数平滑年均值的峰值要滞后,滞后的时间为2～3年.超过70%的中等磁暴出现在太阳活动周的下降段,这表明绝大多数中等磁暴出现在太阳活动周的下降段.通过对中等磁暴平滑月均值与太阳黑子数平滑月均值相位差的计算分析发现,中等磁暴峰值出现的时间比太阳黑子数峰值出现的时间要滞后,不同太阳活动周中等磁暴峰值出现的时间与太阳黑子数峰值时间滞后的程度不同.      

太阳10.7 cm射电辐射流量预报方法初探

引进“相似周”方法,通过两种方式对第23周太阳10.7cm辐射流量（F10.7）月均值进行预报和预报,由“相似周”方法得到的第23周太阳黑子数月均平滑预测值来预测F10.7月均值和直接采用“相似周”方法对F10.7月均平滑值进行预测,通过对两种预报试验结果分析,得到以下结论。（1）两种预测结果与实际都比较吻合,都得到双峰结构。（2）直接通过相似周预测的F10.7的月均值结果较间接通过太阳黑子数的预测结果所推断的F10.7预测结果更接近实际观测结果。（3）使用“相似周”,预报方法,可以给出F10.7比较精细的剖面结构,这是其他普通预报方法很难做到的。     

Comparison of heliospheric conditions near the earth during four recent solar maxima

Statistical properties of the daily averaged values of the solar activity (sunspot numbers, total solar irradiance and 10.7 cm radio emission indices), the solar wind plasma and the interplanetary magnetic field parameters near the Earth’s orbit are investigated for a period from 1964 to 2002 covering the maxima of four solar cycles from 20th to 23rd. Running half-year averages show significant solar cycle variations in the solar activity indices but only marginal and insignificant changes in comparison with background fluctuations for heliospheric bulk plasma and magnetic field parameters. The current 23rd cycle maximum is weaker than 21st and 22nd maxima, but slightly stronger than 20th cycle in most of solar and heliospheric manifestations.     

Study of the influence of solar variability on a regional (Indian) climate: 1901–2007

We use Indian temperature data of more than 100 years to study the influence of solar activity on climate. We study the Sun–climate relationship by averaging solar and climate data at various time scales; decadal, solar activity and solar magnetic cycles. We also consider the minimum and maximum values of sunspot number (SSN) during each solar cycle. This parameter SSN is correlated better with Indian temperature when these data are averaged over solar magnetic polarity epochs (SSN maximum to maximum). Our results indicate that the solar variability may still be contributing to ongoing climate change and suggest for more investigations.     

利用相似周方法预测第24活动周开始时间

利用相似周方法对第24活动周的开始时间与第23活动周下降相后期的太阳黑子数进行了预报.根据第23周已经出现的特征参量和下降相的形态特征,选取9,10,11,15,17和20等六个太阳活动周作为第23周下降相的相似周,对第24周开始时间进行预报.预报结果显示,第24活动周的开始时间为2007年5±1月,黑子数平滑月均极小值为7.1±2.6,第23太阳活动周长度为11.1年.与其他研究者的预报结果相比较,本文给出的结果与文献[11]和[12]及MSFC的结果比较一致.通过对相似周方法在下降相预测太阳活动周结束时间的研究讨论,及对第23周上升阶段的太阳黑子数和F10.7平滑月均值预报结果的评估,可以看出,相似周预报方法在太阳活动周长期预报中是很有应用价值的.      

Forecasting the 23 and 24 solar cycles on the basis of MGM spectrum

A new method of nonlinear spectral analysis (called the method of global minimum: MGM), based on the best presentation (in sense of minimal squares) of a given time data set as a sum of sinusoids whose frequencies, amplitudes and phases are to be determined, has been used to find periodicities in annual Wolf sunspot numbers (W) during the period 1700–1995. The possible future behaviour of the 11-year solar cycle (based on an extrapolation of the calculated model) is also presented. The main characteristics of the 23^rd solar cycle are as follows: the W maximum occurs about 2004, with a peak of nearly 220. An unusually large value of W will occur during the 23^rd cycle, which should be characterised by the longest maximum, specifically, W will be greater than 100 during the 11-year period from 1997 to 2007. The first sharp rise will occur during the period 1996–1998, the second sharp during 2002–2004. The main features of the 24 year cycle are as follows: the next minimum in W, associated with the 24^th solar cycle, should occur in the year 2008 and the maximum in 2014. W is expected to peak at about 180. The minimum value for the 25^th year cycle is expected to occur in the year 2019. It is shown that the accuracy of these predictions depends, first of all, on the extrapolation of the hyperlong harmonic of the calculated polyharmonic model fit of observed annual sunspot numbers during the period 1700–1995. The error bars in the definition of the maximum and minimum epochs can be as large as two years.     

Sunspot cycle 23 descent to an unusual minimum and forecasts for cycle 24 activity

The decay phase of the sunspot cycle 23 exhibited two unusual features. First, it lasted too long. Second, the interplanetary magnetic field intensity at earth orbit reached the lowest value since in situ measurements in space began in October 1963. These physical anomalies significantly altered the early forecasts for the sunspot activity parameters for cycle 24, made by several colleagues. We note that there was a significant change in the solar behavior during cycle 22. We discuss the observed trends and their effect on our empirical solar activity forecast technique, leading to our prediction for cycle 24 parameters; cycle 24 will be only half as active as cycle 23, reaching its peak in May 2013. We speculate on the possible implications of this outcome on future earth climate change and the ensuing socio-economic consequences.     

The asymmetry of the solar cycle: A result of non-linearity

The paper presents a numerical analysis of Wolf sunspot numbers, with emphasis being laid on the asymmetry of the cyclic variation. To that purpose we have used the standard tables of monthly numbers and, in addition to the Fourier transform, we have done an overall analysis of the trend around each maximum. Many of these maxima present an asymmetry, and sometimes the presence of two maxima is evident. The non-linear RLC van der Pol model suggested by Polygiannakis and Moussas [Polygiannakis, J.M., Moussas, X. A nonlinear RLC solar cycle model. Solar Physics 163, 193–203, 1996] can explain many features of the observed asymmetries. Our analysis shows that a consistent deconvolution in two Gaussian curves is possible for each maximum. We may presume that the observed sunspot time series includes a hidden complex structure. This could give some hints of a behavior typical for coupled non-linear oscillators. It is a matter of further interpretations. whether such “oscillators” are just a simple approximation of a much complex phenomenon, or are a sign of another more physically based model like the dynamo model (or other models).     

Solar radiation variability in the last 1400 years recorded in the carbon isotope ratio of a mediterranean sea core

We present measurements and data analysis of the carbon stable isotopes (δ¹³C) in the planktonic Globigerinoides ruber extracted from the GT90/3 shallow water Ionian sea core, dated with high precision. It is commonly accepted that δ¹³C variations in symbiontic foraminifera mainly record the effects of productivity and of photosynthetic activity, varying with the ambient light level. Therefore from this time series we can deduce information on the sea surface illumination at the time of the planktonic foraminifera growth. The profile (359 points) covers the period 590–1979 AD, with a resolution of 3.87 years and it is an extension of the time series (215 points) previously published in this journal. The spectral analysis of the longer time series confirms the presence of the 11 y signal, with amplitude 0.08‰ (peak-to-trough), found in the shorter time series in phase with the sunspot solar cycle; furthermore it shows the presence of two centennial cycles of 100 and 200 years, with amplitude 0.08‰ and 0.02‰ respectively. These components are identified at high significance level by Monte Carlo singular spectrum analysis (MC-SSA). A comparison between the δ¹³C profile and the historical aurorae series (600–1500 AD) shows that the long-term δ¹³C variations are at least partially generated by the solar activity modulation and in phase with the solar output, as represented by the solar wind interaction with the magnetosphere.     

Sunspot numbers during 1736–1739 revisited

Some periods in the sunspot number reconstruction composed by Hoyt and Schatten [Hoyt, D.V., Schatten, K.H. Group Sunspot Numbers: a new solar activity reconstruction. Sol. Phys. 179, 189–219, 1998. Reprinted with figures in Sol. Phys. 181, 491–512, 1998], are based on very few records. For example, there are only a few solar observations during the years 1736–1739. In this paper we intend to improve the reliability of the sunspot numbers reconstruction developed by Hoyt and Schatten for this 4-years period based on information about solar activity published in three journals of that epoch: “Philosophical Transactions”, “Histoire de l’Académie Royale des Sciences”, and “Nova Acta Eruditorum”. We were able to identify 42 papers with solar observations, including 30 with relevant information on sunspots. Based upon this new outlook, a reconstruction of the monthly solar activity for these years is proposed.     

Latitudinal distribution of solar coronal active regions

We studied the cyclic evolution of the latitudinal distribution of solar coronal active regions based on daily images from SOHO EIT for the period 1995–2017. Fully automated software was used, which included the following steps: initial preparation of images in the data series, normalization of histograms and correction of limb brightening, segmentation of images using threshold intensity values obtained from their histograms, scanning of segmented images in heliographic coordinates and obtaining profiles of latitudinal distribution of coronal active regions for each image of the data series. From the output data, we obtained a temporary change in the latitudinal distribution profiles and the migration of activity centers on the solar disk. From the period of minimum activity to the next minimum in both hemispheres, activity centers begin to migrate from high latitudes towards the equator. At the same time, the general center of activity repeatedly changes the direction of migration. The latitudinal distribution of the so-called presence factor of coronal active regions closely resembles the magnetic butterfly diagram, which proves their direct causal relationships. Variations in the presence factor of coronal active regions are correlated with cyclic variations in the sunspot daily numbers.     

Sunspot activity and cosmic ray modulation at 1 a.u. for 1900–2013

The descent of sunspot cycle 23 to an unprecedented minimum of long duration in 2006–2009 led to a prolonged galactic cosmic ray (GCR) recovery to the highest level observed in the instrumental era for a variety of energetic charged particle species on Earth, over a wide range of rigidities. The remarkable GCR increase measured by several ground-based, balloon-borne, and detectors on a satellite is described and discussed. It is accompanied by a decrease in solar wind velocity and interplanetary magnetic field at 1 a.u., reaching the lowest values since measurements of the solar wind began in October 1963; the solar polar field strength (μT) measured at the Wilcox Solar Observatory (WSO) is also significantly reduced compared to prior cycles since the start of the program in 1976, the polar field in the northern hemisphere reversed in June 2012 and again in February 2014, that in the southern hemisphere reversed in July 2013. If updates of WSO data confirm the second reversal in northern solar hemisphere, it would pose a serious challenge to the Dynamo Theory. The long-term change in solar behavior may have begun in 1992, perhaps earlier. The physical underpinnings of these solar changes need to be understood and their effect on GCR modulation processes clarified. The study discusses the recent phenomena in the context of GCR modulation since 1900. These happenings affected our empirical predictions for the key parameters for the next two sunspot cycles (they may be progressively less active than sunspot cycle 24) but it enhanced support for our prediction that solar activity is descending into a Dalton-like grand minimum in the middle of the twentyfirst century, reducing the frequency of the coronal mass ejections; they determine the space weather affecting the quality of life on Earth, radiation dose for hardware and human activities in space as well as the frequency of large Forbush decreases at 1 a.u.     

Sunspot cycle prediction using Warped Gaussian process regression

Solar cycle prediction is a key activity in space weather research. Several techniques have been employed in recent decades in order to try to forecast the next sunspot-cycle maxima and time. In this work, the Gaussian process, a machine-learning technique, is used to make a prediction for the solar cycle 25 based on the annual sunspot number 2.0 data from 1700 to 2018. A variation known as Warped Gaussian process is employed in order to deal with the non-negativity constraint and asymmetrical data distribution. Tests using holdout data yielded a root mean square error of 10.0 within 5 years and 25.0–35.0 within 10 years. Simulations using the predictive distribution were performed to account for the uncertainty in the prediction. Cycle 25 is expected to last from 2019 to 2029, with a peak sunspot number about 117 (110 by the median) occurring most likely in 2024. Thus our method predicts that solar Cycle 25 will be weaker than previous ones, implying a continuing trend of declining solar activity as observed in the past two cycles.     

Impact of rogue active regions on hemispheric asymmetry

The solar dipole moment at activity minimum is a good predictor of the strength of the subsequent solar cycle. Through a systematic analysis using a state-of-the-art $2\times 2$D solar dynamo model, we found that bipolar magnetic regions (BMR) with atypical characteristics can modify the strength of the next cycle via their impact on the buildup of the dipole moment as a sunspot cycle unfolds. In addition to summarizing these results, we present further effects of such “rogue” BMRs. These have the ability to generate hemispheric asymmetry in the subsequent sunspot cycle, since they modify the polar cap flux asymmetry of the ongoing cycle. We found strong correlation between the polar cap flux asymmetry of cycle i and the total pseudo sunspot number asymmetry of cycle $i+1$. Good correlation also appears in the case of the time lag of the hemispheres of cycle $i+1$.