Prediction of the Maximum Amplitude and Timing of Sunspot Cycle 24

Precursor techniques, in particular those using geomagnetic indices, often are used in the prediction of the maximum amplitude for a sunspot cycle. Here, the year 2008 is taken as being the sunspot minimum year for cycle 24. Based on the average aa index value for the year of the sunspot minimum and the preceding four years, we estimate the expected annual maximum amplitude for cycle 24 to be about 92.8±19.6 (1-sigma accuracy), indicating a somewhat weaker cycle 24 as compared to cycles 21 – 23. Presuming a smoothed monthly mean sunspot number minimum in August 2008, a smoothed monthly mean sunspot number maximum is expected about October 2012±4 months (1-sigma accuracy).     

A Preliminary Estimate of the Size of the Coming Solar Cycle 24, based on Ohl’s Precursor Method

For many purposes (e.g., satellite drag, operation of power grids on Earth, and satellite communication systems), predictions of the strength of a solar cycle are needed. Predictions are made by using different methods, depending upon the characteristics of sunspot cycles. However, the method most successful seems to be the precursor method by Ohl and his group, in which the geomagnetic activity in the declining phase of a sunspot cycle is found to be well correlated with the sunspot maximum of the next cycle. In the present communication, the method is illustrated by plotting the 12-month running means aa(min ) of the geomagnetic disturbance index aa near sunspot minimum versus the 12-month running means of the sunspot number Rz near sunspot maximum [aa(min ) versus Rz(max )], using data for sunspot cycles 9 – 18 to predict the Rz(max ) of cycle 19, using data for cycles 9 – 19 to predict Rz(max ) of cycle 20, and so on, and finally using data for cycles 9 – 23 to predict Rz(max ) of cycle 24, which is expected to occur in 2011 – 2012. The correlations were good (∼+0.90) and our preliminary predicted Rz(max ) for cycle 24 is 142±24, though this can be regarded as an upper limit, since there are indications that solar minimum may occur as late as March 2008. (Some workers have reported that the aa values before 1957 would have an error of 3 nT; if true, the revised estimate would be 124±26.) This result of the precursor method is compared with several other predictions of cycle 24, which are in a very wide range (50 – 200), so that whatever may be the final observed value, some method or other will be discredited, as happened in the case of cycle 23.     

A Prediction of Solar Cycle 24 Using a Modified Precursor Method

Based on cycles 17 – 23, linear correlations are obtained between 12-month moving averages of the number of disturbed days when Ap is greater than or equal to 25, called the Disturbance Index (DI), at thirteen selected times (called variate blocks 1, 2,… , each of six-month duration) during the declining portion of the ongoing sunspot cycle and the maximum amplitude of the following sunspot cycle. In particular, variate block 9, which occurs just prior to subsequent cycle minimum, gives the best correlation (0.94) with a minimum standard error of estimation of ± 13, and hindcasting shows agreement between predicted and observed maximum amplitudes to about 10%. As applied to cycle 24, the modified precursor technique yields maximum amplitude of about 124±23 occurring about 45±4 months after its minimum amplitude occurrence, probably in mid to late 2011.     

Predicting maximum sunspot number in solar cycle 24

A few prediction methods have been developed based on the precursor technique which is found to be successful for forecasting the solar activity. Considering the geomagnetic activity aa indices during the descending phase of the preceding solar cycle as the precursor, we predict the maximum amplitude of annual mean sunspot number in cycle 24 to be 111 ± 21. This suggests that the maximum amplitude of the upcoming cycle 24 will be less than cycles 21–22. Further, we have estimated the annual mean geomagnetic activity aa index for the solar maximum year in cycle 24 to be 20.6 ± 4.7 and the average of the annual mean sunspot number during the descending phase of cycle 24 is estimated to be 48 ± 16.8.     

Solar Cycle Predictions based on Solar Activity at Different Solar Latitudes

Many methods of predictions of sunspot maximum number use data before or at the preceding sunspot minimum to correlate with the following sunspot maximum of the same cycle, which occurs a few years later. Kane and Trivedi (Solar Phys. 68, 135, 1980) found that correlations of R _z(max) (the maximum in the 12-month running means of sunspot number R _z) with R _z(min) (the minimum in the 12-month running means of sunspot number R _z) in the solar latitude belt 20° – 40°, particularly in the southern hemisphere, exceeded 0.6 and was still higher (0.86) for the narrower belt > 30° S. Recently, Javaraiah (Mon. Not. Roy. Astron. Soc. 377, L34, 2007) studied the relationship of sunspot areas at different solar latitudes and reported correlations 0.95 – 0.97 between minima and maxima of sunspot areas at low latitudes and sunspot maxima of the next cycle, and predictions could be made with an antecedence of more than 11 years. For the present study, we selected another parameter, namely, SGN, the sunspot group number (irrespective of their areas) and found that SGN(min) during a sunspot minimum year at latitudes > 30° S had a correlation +0.78±0.11 with the sunspot number R _z(max) of the same cycle. Also, the SGN during a sunspot minimum year in the latitude belt (10° – 30° N) had a correlation +0.87±0.07 with the sunspot number R _z(max) of the next cycle. We obtain an appropriate regression equation, from which our prediction for the coming cycle 24 is R _z(max )=129.7±16.3.     

The Relationship between Prediction Accuracy and Correlation Coefficient

The correlation coefficient (r) between the maximum amplitude (R _m) of a sunspot cycle and the preceding minimum aa geomagnetic index (aa _min), in terms of geomagnetic cycle, can be fitted by a sinusoidal function with a four-cycle periodicity superimposed on a declining trend. The prediction index (χ) of the prediction error relative to its estimated uncertainty based on a geomagnetic precursor method can be fitted by a sinusoidal function with a four-and-half-cycle periodicity. A revised prediction relationship is found between the two quantities: χ<1.2 if r varies in a rising trend, and χ>1.2 if r varies in a declining trend. The prediction accuracy of R _m depends on the long-term variation in the correlation. These results indicate that the prediction for the next cycle inferred from this method, R _m(24)=87±23 regarding the 75% level of confidence (1.2-σ), is likely to fail. When using another predictor of sunspot area instead of the geomagnetic index, similar results can be also obtained. Dynamo models will have better predictive powers when having considered the long-term periodicities.     

Sunspot Time Series – Relations Inferred from the Location of the Longest Spotless Segments

Spotless days (i.e., days when no sunspots are observed on the Sun) occur during the interval between the declining phase of the old sunspot cycle and the rising phase of the new sunspot cycle, being greatest in number and of longest continuous length near a new cycle minimum. In this paper, we introduce the concept of the longest spotless segment (LSS) and examine its statistical relation to selected characteristic points in the sunspot time series (STS), such as the occurrences of first spotless day and sunspot maximum. The analysis has revealed statistically significant relations that appear to be of predictive value. For example, for Cycle 24 the last spotless day during its rising phase should be about August 2012 (± 9.1 months), the daily maximum sunspot number should be about 227 (± 50; occurring about January 2014±9.5 months), and the maximum Gaussian smoothed sunspot number should be about 87 (± 25; occurring about July 2014). Using the Gaussian-filtered values, slightly earlier dates of August 2011 and March 2013 are indicated for the last spotless day and sunspot maximum for Cycle 24, respectively.     

The 27-Day and 22-Year Cycles in Solar and Geomagnetic Activity

We study the evolution of the longitudinal asymmetry in solar activity through the wave packet technique applied to the period domain of 25 – 31 days (centered at the 27-day solar rotation period) for the sunspot number and geomagnetic aa index. We observe the occurrence of alternating smaller and larger amplitudes of the 11-year cycle, resulting in a 22-year periodicity in the 27-day signal. The evolution of the 22-year cycle shows a change of regime around the year 1912 when the 22-year period disappears from the sunspot number series and appears in the aa index. Other changes, such as a change in the correlation between solar and geomagnetic activity, took place at the same time. Splitting the 27-day frequency domain of aa index shows an 11-year cycle for higher frequencies and a pure22-year cycle for lower frequencies, which we attribute to higher latitude coronal holes. This evidence is particularly clear after 1940, which is another benchmark in the evolution of the aa index. We discuss briefly the mechanisms that could account for the observed features of the 22-year cycle evolution.     

On the level of skill in predicting maximum sunspot number: A comparative study of single variate and bivariate precursor techniques

Precursor prediction techniques have generally performed well in predicting the maximum amplitude of sunspot cycles, based on cycles 10–21. Single variate methods based on minimum sunspot amplitude have reliably predicted the size of the sunspot cycle 9 out of 12 times, where a reliable prediction is defined as one having an observed maximum amplitude within the prediction interval (determined from the average error). On the other hand, single variate methods based on the size of the geomagnetic minimum have reliably predicted the size of the sunspot cycle 8 of 10 times (geomagnetic data are only available since about cycle 12). Bivariate prediction methods have, thus far, performed flawlessly, giving reliable predictions 10 out of 10 times (bivariate methods are based on sunspot and geomagnetic data). For cycle 22, single variate methods (based on geomagnetic data) suggest a maximum amplitude of about 170 ± 25, while bivariate methods suggest a maximum amplitude of about 140 ± 15; thus, both techniques suggest that cycle 22 will be of smaller maximum amplitude than that observed during cycle 19, and possibly even smaller than that observed for cycle 21. Compared to the mean cycle, cycle 22 is presently behaving as if it is a + 2.6 cycle (maximum amplitude about 225). It appears then that either cycle 22 will be the first cycle not to be reliably predicted by the combined precursor techniques (i.e., cycle 22 is an anomaly, a statistical outlier) or the deviation of cycle 22 relative to the mean cycle will substantially decrease over the next 18 months. Because cycle 22 is a large amplitude cycle, maximum smoothed sunspot number is expected to occur early in 1990 (between December 1989 and May 1990).     

Predicting the Amplitude of a Solar Cycle Using the North – South Asymmetry in the Previous Cycle: II. An Improved Prediction for Solar Cycle 24

Recently, using Greenwich and Solar Optical Observing Network sunspot group data during the period 1874 – 2006, Javaraiah (Mon. Not. Roy. Astron. Soc. 377, L34, 2007: Paper I), has found that: (1) the sum of the areas of the sunspot groups in 0° – 10° latitude interval of the Sun’s northern hemisphere and in the time-interval of −1.35 year to +2.15 year from the time of the preceding minimum of a solar cycle n correlates well (corr. coeff. r=0.947) with the amplitude (maximum of the smoothed monthly sunspot number) of the next cycle n+1. (2) The sum of the areas of the spot groups in 0° – 10° latitude interval of the southern hemisphere and in the time-interval of 1.0 year to 1.75 year just after the time of the maximum of the cycle n correlates very well (r=0.966) with the amplitude of cycle n+1. Using these relations, (1) and (2), the values 112±13 and 74±10, respectively, were predicted in Paper I for the amplitude of the upcoming cycle 24. Here we found that the north – south asymmetries in the aforementioned area sums have a strong ≈44-year periodicity and from this we can infer that the upcoming cycle 24 will be weaker than cycle 23. In case of (1), the north – south asymmetry in the area sum of a cycle n also has a relationship, say (3), with the amplitude of cycle n+1, which is similar to (1) but more statistically significant (r=0.968) like (2). By using (3) it is possible to predict the amplitude of a cycle with a better accuracy by about 13 years in advance, and we get 103±10 for the amplitude of the upcoming cycle 24. However, we found a similar but a more statistically significant (r=0.983) relationship, say (4), by using the sum of the area sum used in (2) and the north – south difference used in (3). By using (4) it is possible to predict the amplitude of a cycle by about 9 years in advance with a high accuracy and we get 87±7 for the amplitude of cycle 24, which is about 28% less than the amplitude of cycle 23. Our results also indicate that cycle 25 will be stronger than cycle 24. The variations in the mean meridional motions of the spot groups during odd and even numbered cycles suggest that the solar meridional flows may transport magnetic flux across the solar equator and potentially responsible for all the above relationships. The author did a major part of this work at the Department of Physics and Astronomy, UCLA, 430 Portola Plaza, Los Angeles, CA 90095-1547, USA.     

Gnevyshev Peaks and Gaps for Coronal Mass Ejections of Different Widths Originating in Different Solar Position Angles

The sunspot number series at the peak of sunspot activity often has two or three peaks (Gnevyshev peaks; Gnevyshev, Solar Phys. 1, 107, 1967; Solar Phys. 51, 175, 1977). The sunspot group number (SGN) data were examined for 1997 – 2003 (part of cycle 23) and compared with data for coronal mass ejection (CME) events. It was noticed that they exhibited mostly two Gnevyshev peaks in each of the four latitude belts 0° – 10°, 10° – 20°, 20 ° – 30°, and > 30°, in both N (northern) and S (southern) solar hemispheres. The SGN were confined to within latitudes ± 50° around the Equator, mostly around ± 35°, and seemed to occur later in lower latitudes, indicating possible latitudinal migration as in the Maunder butterfly diagrams. In CMEs, less energetic CMEs (of widths < 71°) showed prominent Gnevyshev peaks during sunspot maximum years in almost all latitude belts, including near the poles. The CME activity lasted longer than the SGN activity. However, the CME peaks did not match the SGN peaks and were almost simultaneous at different latitudes, indicating no latitudinal migration. In energetic CMEs including halo CMEs, the Gnevyshev peaks were obscure and ill-defined. The solar polar magnetic fields show polarity reversal during sunspot maximum years, first at the North Pole and, a few months later, at the South Pole. However, the CME peaks and gaps did not match with the magnetic field reversal times, preceding them by several months, rendering any cause – effect relationship doubtful.     

Short-Term Periodicities in Sunspot Activity and Flare Index Data during Solar Cycle 23

The short-term periodicities in sunspot numbers, sunspot areas, and flare index data are investigated in detail using the Date Compensated Discrete Fourier Transform (DCDFT) for the full disk of the Sun separately over the rising, the maximum, and the declining portions of solar cycle 23 (1996 – 2006). While sunspot numbers and areas show several significant periodicities in a wide range between 23.1 and 36.4 days, the flare index data do not exhibit any significant periodicity. The earlier conclusion of Pap, Tobiska, and Bouwer (1990, Solar Phys. 129, 165) and Kane (2003, J. Atmos. Solar-Terr. Phys. 65, 1169), that the 27-day periodicity is more pronounced in the declining portion of a solar cycle than in the rising and maximum ones, seems to be true for sunspot numbers and sunspot area data analyzed here during solar cycle 23.     

The Behavior of Sunspot Contrast during Cycle 23

Results are presented from a study of various sunspot contrast parameters in broadband red (672.3 nm) Cartesian full-disk digital images taken at the San Fernando Observatory (SFO) over eight years, 1997 – 2004, of the twenty-third sunspot cycle. A subset of over 2700 red sunspots was analyzed and values of average and maximum sunspot contrast as well as maximum umbral contrast were compared to various sunspot parameters. Average and maximum sunspot contrasts were found to be significantly correlated with sunspot area (r _s=− 0.623 and r _s=− 0.714, respectively). Maximum umbral contrast was found to be significantly correlated with umbral area (r _s=− 0.535). These results are in agreement with the works of numerous other authors. No significant dependence was detected between average contrast, maximum contrast, or maximum umbral contrast during the rising phase of the solar cycle (r _s=0.024, r _s=0.033, and r _s=0.064, respectively). During the decay phase, no significant correlation was found between average contrast or maximum contrast and time (r _s=− 0.057 and r _s=0.009, respectively), with a weak dependence seen between maximum umbral contrast and cycle (r _s=0.102).     

Study of nuclear decays during a solar eclipse: Thule Greenland 2008

A new prediction technique based on logarithmic values is proposed to predict the maximum amplitude (R _m) of a solar cycle from the preceding minimum aa geomagnetic index (aa _min). The correlation between lnR _m and lnaa _min (r=0.92) is slightly stronger than that between R _m and aa _min (r=0.90). From this method, cycle 24 is predicted to have a peak size of R _m(24)=81.7(1±13.2%). If the suggested error in aa (3 nT) before 1957 is corrected, the correlation coefficient between R _m and aa _min (r=0.94) will be slightly higher, and the peak of cycle 24 is predicted much lower, R _m(24)=52.5±13.1. Therefore, the prediction of R _m based on the relationship between R _m and aa _min depends greatly on the accurate measurement of aa.     

PREDICTIONS OF THE FEATURES FOR SUNSPOT CYCLE 23

Smoothed monthly mean Ap indices are decomposed into two components (Ap) _c and (Ap) _n. The former is directly correlated with the current sunspot numbers, while the latter is shown to achieve its maximum correlation with the sunspot numbers after some time lag. This latter property is used to develop a method for predicting the sunspot maximum based on the observed value of (Ap) _n maximum which occurs during the preceding cycle. The value of R _M for cycle 23 predicted by this method is 149.3 ± 19.9. A method to estimate the rise time (from solar minimum to maximum) has been developed (based on analyses of Hathaway, Wilson, and Reichmann, 1994) and yields a value of 4.2 years. Using an estimate that the minimum between cycles 22 and 23 occurred in May 1996, it is predicted that the sunspot maximum for cycle 23 will occur in July 2000.     

Two Regularities in the Coronal Green-Line Brightness – Magnetic Field Coupling and the Heating of the Corona

To study the quantitative relationship between the brightness of the coronal green line 530.5 nm Fe xiv and the strength of the magnetic field in the corona, we have calculated the cross-correlation of the corresponding synoptic maps for the period 1977 – 2001. The maps of distribution of the green-line brightness I were plotted using every-day monitoring data. The maps of the magnetic field strength B and the tangential B _t and radial B _r field components at the distance 1.1 R _⊙ were calculated under potential approximation from the Wilcox Solar Observatory (WSO) photospheric data. It is shown that the correlation I with the field and its components calculated separately for the sunspot formation zone ±30° and the zone 40 – 70° has a cyclic character, the corresponding correlation coefficients in these zones changing in anti-phase. In the sunspot formation zone, all three coefficients are positive and have the greatest values near the cycle minimum decreasing significantly by the maximum. Above 40°, the coefficients are alternating in sign and reach the greatest positive values at the maximum and the greatest negative values, at the minimum of the cycle. It is inferred that the green-line emission in the zone ±30° is mainly controlled by B _t, probably due to the existence of low arch systems. In the high-latitude zone, particularly at the minimum of the cycle, an essential influence is exerted by B _r, which may be a manifestation of the dominant role of large-scale magnetic fields. Near the activity minimum, when the magnetic field organization is relatively simple, the relation between I and B for the two latitudinal zones under consideration can be represented as a power-law function of the type I ∝ B ^q. In the sunspot formation zone, the power index q is positive and varies from 0.75 to 1.00. In the zone 40 – 70°, it is negative and varies from −0.6 to −0.8. It is found that there is a short time interval approximately at the middle of the ascending branch of the cycle, when the relationship between I and B vanishes. The results obtained are considered in relation to various mechanisms of the corona heating.     

An Estimate for the Size of Sunspot Cycle 24

For the sunspot cycles in the modern era (cycle?10 to the present), the ratio of R _Z(max)/R _Z(36th month) equals 1.26±0.22, where R _Z(max) is the maximum amplitude of the sunspot cycle?using smoothed monthly mean sunspot number and R _Z(36th month) is the smoothed monthly mean sunspot number 36 months after cycle?minimum. For the current sunspot cycle?24, the 36th month following the cycle?minimum occurred in November 2011, measuring?61.1. Hence, cycle?24 likely will have a maximum amplitude of about 77.0±13.4 (the one-sigma prediction interval), a value well below the average R _Z(max) for the modern era sunspot cycles (about 119.7±39.5).     

A Prediction of the Peak Sunspot Number for Solar Cycle 23

We use a precursor technique based on the geomagneticaa index during the decline (last 30%) of solar cycle 22 to predict a peak sunspot number of 158 (± 18) for cycle 23, under the assumption that solar minimum occurred in May 1996. This method appears to be as reliable as those that require a year of data surrounding the geomagnetic minimum, which typically follows the smoothed sunspot minimum by about six months.     

An Early Prediction of Maximum Sunspot Number in Solar Cycle 24

A non-linear coupling function between sunspot maxima and aa minima modulations has been found as a result of a wavelet analysis of geomagnetic index aa and Wolf sunspot number yearly means since 1844. It has been demonstrated that the increase of these modulations for the past 158 years has not been steady, instead, it has occurred in less than 30 years starting around 1923. Otherwise sunspot maxima have oscillated about a constant level of 90 and 141, prior to 1923 and after 1949, respectively. The relevance of these findings regarding the forecasting of solar activity is analyzed here. It is found that if sunspot cycle maxima were still oscillating around the 141 constant value, then the Gnevyshev–Ohl rule would be violated for two consecutive even–odd sunspot pairs (22–23 and 24–25) for the first time in 1700 years. Instead, we present evidence that solar activity is in a declining episode that started about 1993. A value for maximum sunspot number in solar cycle 24 (87.5±23.5) is estimated from our results.     

Latitude Dependence of the “Gnevyshev” Peaks and Gaps

The occurrence of double peaks near the maximum of sunspot activity was first emphasized by Gnevyshev (Solar Phys. 1, 107, 1967) for the peak years of solar cycle 19 (1954 – 1964). In the present analysis, it is shown that double peaks in sunspot numbers were clearly visible in solar latitudes 10 – 30° N but almost absent in the southern latitudes, where some single peaks were observed out of phase by several months from any of the peaks in the northern latitudes. The spacing between the double peaks increased from higher to lower northern latitudes, hinting at latitudinal migration. In the next cycle 20 (1965 – 1976), which was of about half the strength of cycle 19, no clear-cut double peaks were seen, and the prominent peak in the early part of 1967 in the northern latitudes was seen a few months later in the southern latitudes. A direct relationship of Gnevyshev peaks with changes in the solar polar magnetic fields seems to be dubious. The commencements do not match.