Group Sunspot Numbers: A New Solar Activity Reconstruction

In this paper, we construct a time series known as the Group Sunspot Number. The Group Sunspot Number is designed to be more internally self-consistent (i.e., less dependent upon seeing the tiniest spots) and less noisy than the Wolf Sunspot Number. It uses the number of sunspot groups observed, rather than groups and individual sunspots. Daily, monthly, and yearly means are derived from 1610 to the present. The Group Sunspot Numbers use 65941 observations from 117 observers active before 1874 that were not used by Wolf in constructing his time series. Hence, we have calculated daily values of solar activity on 111358 days for 1610–1995, compared to 66168 days for the Wolf Sunspot Numbers. The Group Sunspot Numbers also have estimates of their random and systematic errors tabulated. The generation and preliminary analysis of the Group Sunspot Numbers allow us to make several conclusions: (1) Solar activity before 1882 is lower than generally assumed and consequently solar activity in the last few decades is higher than it has been for several centuries. (2) There was a solar activity peak in 1801 and not 1805 so there is no long anomalous cycle of 17 years as reported in the Wolf Sunspot Numbers. The longest cycle now lasts no more than 15 years. (3) The Wolf Sunspot Numbers have many inhomogeneities in them arising from observer noise and this noise affects the daily, monthly, and yearly means. The Group Sunspot Numbers also have observer noise, but it is considerably less than the noise in the Wolf Sunspot Numbers. The Group Sunspot Number is designed to be similar to the Wolf Sunspot Number, but, even if both indices had perfect inputs, some differences are expected, primarily in the daily values.     

Group Sunspot Numbers: Sunspot Cycle Characteristics

We examine the `Group' sunspot numbers constructed by Hoyt and Schatten to determine their utility in characterizing the solar activity cycle. We compare smoothed monthly Group sunspot numbers to Zürich (International) sunspot numbers, 10.7-cm radio flux, and total sunspot area. We find that the Zürich numbers follow the 10.7-cm radio flux and total sunspot area measurements only slightly better than the Group numbers. We examine several significant characteristics of the sunspot cycle using both Group numbers and Zürich numbers. We find that the `Waldmeier Effect' – the anti-correlation between cycle amplitude and the elapsed time between minimum and maximum of a cycle – is much more apparent in the Zürich numbers. The `Amplitude–Period Effect' – the anti-correlation between cycle amplitude and the length of the previous cycle from minimum to minimum – is also much more apparent in the Zürich numbers. The `Amplitude–Minimum Effect' – the correlation between cycle amplitude and the activity level at the previous (onset) minimum is equally apparent in both the Zürich numbers and the Group numbers. The `Even–Odd Effect' – in which odd-numbered cycles are larger than their even-numbered precursors – is somewhat stronger in the Group numbers but with a tighter relationship in the Zürich numbers. The `Secular Trend' – the increase in cycle amplitudes since the Maunder Minimum – is much stronger in Group numbers. After removing this trend we find little evidence for multi-cycle periodicities like the 80-year Gleissberg cycle or the two- and three-cycle periodicities. We also find little evidence for a correlation between the amplitude of a cycle and its period or for a bimodal distribution of cycle periods. We conclude that the Group numbers are most useful for extending the sunspot cycle data further back in time and thereby adding more cycles and improving the statistics. However, the Zürich numbers are slightly more useful for characterizing the on-going levels of solar activity.     

Reconstruction of Monthly and Yearly Group Sunspot Numbers From Sparse Daily Observations

Some periods before 1820 are poorly covered by sunspot observations. In addition to apparent, long observational gaps, there are also periods when there are only few sparse daily sunspot observations during a long time. It is important to estimate the reliability of the monthly and yearly mean sunspot values obtained from such sparse daily data. Here we suggest a new method to estimate the reliability of individual monthly means. The method is based on comparing the actual sparse data (sample population) to the well-measured sunspot data in 1850–1996 (reference population), and assumes that the statistical properties of sunspot activity remain similar throughout the entire period. For each sample population we first found those months in the reference population that contain the same data set, and constructed the statistical distribution of the corresponding monthly means. The mean and standard error of this distribution represent the mean and uncertainty of a monthly mean sunspot number reconstructed from sparse daily observations. The simple arithmetic mean of daily values can be adequately applied for months which contain more than 4–5 evenly distributed daily observations. However, the reliability of monthly means for less covered months has to be estimated more carefully. Using the estimated, new monthly values, we have also calculated the weighted annual sunspot numbers.     

The Schwabe and Gleissberg Periods in the Wolf Sunspot Numbers and the Group Sunspot Numbers

Three wavelet functions: the Morlet wavelet, the Paul wavelet, and the DOG wavelet have been respectively performed on both the monthly Wolf sunspot numbers (R_z) from January 1749 to May 2004 and the monthly group sunspot numbers (R_g) from June 1795 to December 1995 to study the evolution of the Gleissberg and Schwabe periods of solar activity. The main results obtained are (1) the two most obvious periods in both the R_z and R_g are the Schwabe and Gleissberg periods. The Schwabe period oscillated during the second half of the eighteenth century and was steady from the 1850s onward. No obvious drifting trend of the Schwabe period exists. (2) The Gleissberg period obviously drifts to longer periods the whole consideration time, and the drifting speed of the Gleissberg period is larger for R_z than for R_g. (3) Although the Schwabe-period values for R_z and R_g are about 10.7 years, the value for R_z seems slightly larger than that for R_g. The Schwabe period of R_z is highly significant after the 1820s, and the Schwabe period of R_g is highly significant over almost the whole consideration time except for about 20 years around the 1800s. The evolution of the Schwabe period for both R_z and R_g in time is similar to each other. (4) The Gleissberg period in R_z and R_g is highly significant during the whole consideration time, but this result is unreliable at the two ends of each of the time series of the data. The evolution of the Gleissberg period in R_z is similar to that in R_g.     

Low Dimensional Chaos from the Group Sunspot Numbers

We examine the nonlinear dynamical properties of the monthly smoothed group sunspot number Rg and find that the solar activity underlying the time series of Rg is globally governed by a low-dimensional chaotic attractor. This finding is consistent with the nonlinear study results of the monthly Wolf sunspot numbers. We estimate the maximal Lyaponuv exponent (MLE) for the Rg series to be positive and to equal approximately 0.0187 ± 0.0023 (month-1). Thus, the Lyaponuv time or predictability time of the chaotic motion is obtained to be about 4.46 ± 0.5 years, which is slightly different with the predictability time obtained from Rz. However, they both indicate that solar activity forecast should be done only for a short to medium term due to the intrinsic complexity of the time behavior concerned.     

Reconstruction of the Heliospheric Current Sheet Tilts Using Sunspot Numbers

We have investigated the correlation between the relative sunspot number and tilt of the heliospheric current sheet (HCS) in solar cycles 21–23. Strong and highly significant positive correlation (r > 0.8, P < 0.001) was found for corresponding data in the time interval from May 1976 through December 2004. Cross-correlation analysis does not reveal any time shift between the data sets. Reconstructed values of the HCS tilt, for the time interval before 1976, are found using sunspot numbers. To take different amplitude of solar cycles into account they were then normalized to zero in the minima of the solar activity and to average in solar cycles 21–23 maximal calculated HCS tilt in the maxima. These normalized reconstructed HCS data are compared with the angular positions of the brightest coronal streamers observed during total solar eclipses in 1870–2002, and their agreement is better for the minima of the solar activity than for the maxima.     

Synchronization of Sunspot Numbers and Sunspot Areas

Sunspot activity is usually described by either sunspot numbers or sunspot areas. The smoothed monthly mean sunspot numbers (SNs) and the smoothed monthly mean areas (SAs) in the time interval from November 1874 to September 2007 are used to analyze their phase synchronization. Both the linear method (fast Fourier transform) and some nonlinear approaches (continuous wavelet transform, cross-wavelet transform, wavelet coherence, cross-recurrence plot, and line of synchronization) are utilized to show the phase relation between the two series. There is a high level of phase synchronization between SNs and SAs, but the phase synchronization is detected only in their low-frequency components, corresponding to time scales of about 7 to 12 years. Their high-frequency components show a noisy behavior with strong phase mixing. Coherent phase variables should exist only for a frequency band with periodicities around the dominating 11-year cycle for SNs and SAs. There are some small phase differences between them. SNs lag SAs during most of the considered time interval, and they are in general more asynchronous around the minimum and maximum times of a cycle than at the ascending and descending phases.     

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

The method of active-day fraction (ADF) was proposed recently to calibrate different solar observers to 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-to-moderate activity. The ADF method becomes inapplicable 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.     

大黑子群对太阳总辐照度的影响

本文用云南天文台在第２２周太阳活动峰年期间拍摄到的大太阳黑子群照相资料，太阳黑子目视描述资料，以及Ｎｉｍｂｕｓ—７卫星上辐射计测量的太阳总辐照度，分别计算了太阳总辐射照度与大黑子群的本影视面积，大黑子群全群视面积和日面上全部黑子的总视面积的相关系数。结果表明，太阳总辐射照度与这三种视面积均存在强的负相关。其中与大黑子群本影视面积的相关最强，其次是与全群视面积的相关，最后是与日面上全部黑子的总视面积的相关。并对以上结果和其它有关结果进行了分析和讨论。     

Long-Term Variations of Solar Differential Rotation and Sunspot Activity: Revisited

Long-term variations of solar differential rotation and sunspot activity are investigated through re-analyzing the data on parameters of the differential-rotation law obtained by Makarov, Tlatov, and Callebaut (Solar Phys. 170, 373, 1997), Javaraiah, Bertello, and Ulrich (Astrophys. J. 626, 579, 2005a; Solar Phys. 232, 25, 2005b), and Javaraiah et al. (Solar Phys. 257, 61, 2009). Our results indicate that the solar-surface-rotation rate at the Equator (indicated by the A-parameter of the standard solar-rotation law) shows a secular decrease since Cycle 12 onwards, given by about 1?–?1.5×10^?3 (deg?day^?1?year^?1). The B-parameter of the standard differential-rotation law seems to also show a secular decrease since Cycle 12 onwards, but of weak statistical significance. The rotation rate averaged over latitudes 0^°?–?40^° does not show a secular trend of statistical significance. Moreover, the average sunspot area shows a secular increase of statistical significance since Cycle 12 onwards, while a negative correlation is found between the level of sunspot activity (indicated by the average sunspot area) and the solar equatorial rotation on long-term scales.     

Long-Term Variations in Solar Differential Rotation and Sunspot Activity

The solar equatorial rotation rate, determined from sunspot group data during the period 1879–2004, decreased over the last century, whereas the level of activity has increased considerably. The latitude gradient term of the solar rotation shows a significant modulation of about 79 year, which is consistent with what is expected for the existence of the Gleissberg cycle. Our analysis indicates that the level of activity will remain almost the same as the present cycle during the next few solar cycles (i.e., during the current double Hale cycle), while the length of the next double Hale cycle in sunspot activity is predicted to be longer than the current one. We find evidence for the existence of a weak linear relationship between the equatorial rotation rate and the length of sunspot cycle. Finally, we find that the length of the current cycle will be as short as that of cycle 22, indicating that the present Hale cycle may be a combination of two shorter cycles. Presently working for the Mt. Wilson Solar Archive Digitization Project at UCLA.     

Short-Term Correlation of Solar Activity and Sunspot: Evidence of Lifetime Increase

We perform a nonlinear study of the short-term correlation properties of the solar activity (daily range) in order to reveal their long-life variations. We estimate the lifetime of the high-frequency component of a Markov-type signal when the high-frequency component is modulated by a slowly varying multiplicative factor. This treatment is applied to different series of solar activity: Wolf Sunspot numbers (WSN), Sunspot Group numbers (SGN), and Royal Greenwich Observatory (RGO) sunspot group series. We obtain that all the lifetime estimates exhibit similar temporal variations that agree with the variations of the sunspot lifetimes directly measured from the RGO data and those of the sunspot areas. An increase of lifetimes by a factor 1.4 is observed from 1915 to 1940. At the same time, a stable ratio is observed between the sunspot group’s maximal area and the lifetime, confirming the Gnevyshev–Waldmeier-type relationship. The analysis identifies also time intervals where the homogeneity of the different time series may be questioned.     

Cosmic Ray Intensity Variations in Relation with Solar Flare Index and Sunspot Numbers

From the monthly data of cosmic ray intensity (CRI), sunspot numbers (SSN) and solar flare index (SFI), an attempt has been made to study the relationship between CRI and solar activity (SA) parameters SSN and SFI. The correlation between SA parameters and CRI for different neutron monitoring stations having low, middle and high cut-off rigidity has been investigated. The anti-correlation between SA and CRI is found to exist with some time lag. Based on the method of minimizing correlation coefficient and time-delayed component method, the observed time-lag between SA parameters (SSN and SFI) and CRI has been found to be large for odd solar cycles in comparison to even solar cycles. The results of time-lag analysis between CRI and SSN and between CRI-SFI have also been compared. The findings of correlative study between CRI and SSN are in agreement with earlier results, while the CRI-SFI relationship provides new insights to understand the solar modulation of cosmic rays.     

A Relationship Between the Solar Rotation and Activity Analysed by Tracing Sunspot Groups

The sunspot position published in the data bases of the Greenwich Photoheliographic Results (GPR), the US Air Force Solar Optical Observing Network and National Oceanic and Atmospheric Administration (USAF/NOAA), and of the Debrecen Photoheliographic Data (DPD) in the period 1874 to 2016 were used to calculate yearly values of the solar differential-rotation parameters \(A\) and \(B\). These differential-rotation parameters were compared with the solar-activity level. We found that the Sun rotates more differentially at the minimum than at the maximum of activity during the epoch 1977?–?2016. An inverse correlation between equatorial rotation and solar activity was found using the recently revised sunspot number. The secular decrease of the equatorial rotation rate that accompanies the increase in activity stopped in the last part of the twentieth century. It was noted that when a significant peak in equatorial rotation velocity is observed during activity minimum, the next maximum is weaker than the previous one.     

太阳黑子数的子波分析

本文讨论了子波变换用于信号突变检测的原理，用它分析了１７００－１９９３年间的太阳黑子数的年均值．精确地检测到了太阳活动的突变点，用相邻两个突变点的时间长度求得了不同尺度下太阳黑子变化的周期．结果表明：利用子波变换检测太阳黑子周期与传统方法相比具有独到之处．