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.     

Analysis of the 1.3–1.7 yr Oscillation Relationship between Solar and Geomagnetic Activities

The study on the 1.3–1.7 yr period of the solar and geomagnetic activities is very important for understanding the possible physical processes in the solar-terrestrial coupling system. The sunspot is the most prominent magnetic field structure in the solar photosphere, and the Ap index is an important indicator for the global geomagnetic activity level. The 1.3–1.7 yr period for the sunspot number and the geomagnetic Ap index is obtained by the synchro-squeezing wavelet transform, and the phase relationship between them is studied by the cross-correlation analysis. The main results are as follows: (1) The 1.3–1.7 yr period of the geomagnetic Ap index and sunspot number exhibits an intermittent evolutionary characteristics, and changes continuously with the time; (2) the geomagnetic Ap index has a higher periodic component in the odd solar cycles than the neighboring even solar cycles, which is characterized by fluctuations; (3) the phase relationship between the geomagnetic Ap index and the sunspot number is not always invariant, in most cases the geomagnetic Ap index lags behind the sunspot number, except in the 18th and 22th solar cycles.     

Polar coronal holes and solar cycles

The relationship between the geomagnetic activity of the three years preceding a sunspot minimum and the peak of the next sunspot maximum confirms the polar origin of the solar wind during one part of the solar cycle. Pointing out that the polar holes have a very small size or disappear at the time of the polar field reversal, we suggest a low latitude origin of the solar wind at sunspot maximum and we describe the cycle variation of solar wind and geomagnetic activity. In addition we note a close relationship between the maximum level of the geomagnetic activity reached few years before a solar minimum and its level at the next sunspot maximum. Studying separately the effects of both the low latitude holes and the solar activity, we point out the possibility of predicting both the level of geomagnetic activity and the sunspot number at the next sunspot maximum. As a conclusion we specify the different categories of phenomena contributing to a solar cycle.     

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.     

Preferred longitudes of sunspot groups and high-speed solar wind streams: Evidence for a ‘solar memory’

Correlation analysis of the mean longitude distribution of sunspot groups (taken from the Greenwich Photoheliographic Results) and high-speed solar wind streams (inferred from the C9 index for geomagnetic disturbances) with the Bartels rotation period P = 27.0 days shows anti-correlation for individual cycles.In particular, the longitudes of post-maximum stable streams of cycle 18 and 19 are well anticorrelated with the preferred longitudes of sunspot groups during the maximum activity periods of these cycles. This is further analyzed using the daily Zürich sunspot number, R, between 1932 and 1980, which reveals a conspicuous similarity of cycle 18 and 19 as well as cycle 20 and 21.We conclude that there is a solar memory for preferred longitudes of activity extending at least over one, probably two cycles (i.e. one magnetic cycle of 22 years). We conjecture that this memory extends over longer intervals of time as a long-term feature of solar activity.     

What the Sunspot Record Tells Us About Space Climate

The records concerning the number, sizes, and positions of sunspots provide a direct means of characterizing solar activity over nearly 400 years. Sunspot numbers are strongly correlated with modern measures of solar activity including: 10.7-cm radio flux, total irradiance, X-ray flares, sunspot area, the baseline level of geomagnetic activity, and the flux of galactic cosmic rays. The Group Sunspot Number provides information on 27 sunspot cycles, far more than any of the modern measures of solar activity, and enough to provide important details about long-term variations in solar activity or “Space Climate.” The sunspot record shows: 1) sunspot cycles have periods of 131± 14 months with a normal distribution; 2) sunspot cycles are asymmetric with a fast rise and slow decline; 3) the rise time from minimum to maximum decreases with cycle amplitude; 4) large amplitude cycles are preceded by short period cycles; 5) large amplitude cycles are preceded by high minima; 6) although the two hemispheres remain linked in phase, there are significant asymmetries in the activity in each hemisphere; 7) the rate at which the active latitudes drift toward the equator is anti-correlated with the cycle period; 8) the rate at which the active latitudes drift toward the equator is positively correlated with the amplitude of the cycle after the next; 9) there has been a significant secular increase in the amplitudes of the sunspot cycles since the end of the Maunder Minimum (1715); and 10) there is weak evidence for a quasi-periodic variation in the sunspot cycle amplitudes with a period of about 90 years. These characteristics indicate that the next solar cycle should have a maximum smoothed sunspot number of about 145 ± 30 in 2010 while the following cycle should have a maximum of about 70 ± 30 in 2023.     

Prediction of Solar Cycle 24 Using Geomagnetic Precursors: Validation and Update

In the previous study (Dabas et al. in Solar Phys. 250, 171, 2008), to predict the maximum sunspot number of the current solar cycle 24 based on the geomagnetic activity of the preceding sunspot minimum, the Ap index was used which is available from the last six to seven solar cycles. Since a longer series of the aa index is available for more than the last 10 – 12 cycles, the present study utilizes aa to validate the earlier prediction. Based on the same methodology, the disturbance index (DI), which is the 12-month moving average of the number of disturbed days (aa≥50), is computed at thirteen selected times (called variate blocks 1,2,…,13; each of them in six-month duration) during the declining portion of the ongoing sunspot cycle. Then its correlation with the maximum sunspot number of the following cycle is evaluated. As in the case of Ap, variate block 9, which occurs exactly 48 months after the current cycle maximum, gives the best correlation (R=0.96) with a minimum standard error of estimation (SEE) of ± 9. As applied to cycle 24, the aa index as precursor yields the maximum sunspot number of about 120±16 (the 90% prediction interval), which is within the 90% prediction interval of the earlier prediction (124±23 using Ap). Furthermore, the same method is applied to an expanded range of cycles 11 – 23, and once again variate block 9 gives the best correlation (R=0.95) with a minimum SEE of ± 13. The relation yields the modified maximum amplitude for cycle 24 of about 131±20, which is also close to our earlier prediction and is likely to occur at about 43±4 months after its minimum (December 2008), probably in July 2012 (± 4 months).     

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).     

Ten cycles of solar and geomagnetic activity

Series of 110 years of sunspot numbers and indices of geomagnetic activity are used with 17 years of solar wind data in order to study through solar cycles both stream and shock event solar activity. According to their patterns on Bartels diagrams of geomagnetic indices, stable wind streams and transient solar activities are separated from each other. Two classes of stable streams are identified: equatorial streams occurring sporadically, for several months, during the main phase of sunspot cycles and both polar streams established, for several years, at each cycle, before sunspot minimum. Polar streams are the first activity of solar cycles. For study of the relationship between transient geomagnetic phenomena and sunspot activity, we raise the importance of the contribution, at high spot number, of severe storms and, at low spot number, of short lived and unstable streams. Solar wind data are used to check and complete the above results. As a conclusion, we suggest a unified scheme of solar activity evolution with a starting point every eleventh year, a total duration of 17 years and an overlapping of 6 years between the first and the last phase of both successive series of phenomena: first, from polar field reversal to sunspot minimum, a phase of polar wind activity of the beginning cycle is superimposed on the weak contribution of shock events of the ending cycle; secondly, an equatorial phase mostly of shock events is superimposed on a variable contribution of short lived and sporadic stable equatorial stream activities; and thirdly a phase of low latitude shock events is superimposed on the polar stream interval of the following cycle.     

Search for relationship between duration of the extended solar cycles and amplitude of sunspot cycle

Duration of the extended solar cycles is taken into the consideration. The beginning of cycles is counted from the moment of polarity reversal of large-scale magnetic field in high latitudes, occurring in the sunspot cycle n till the minimum of the cycle n + 2. The connection between cycle duration and its amplitude is established. Duration of the “latent” period of evolution of extended cycle between reversals and a minimum of the current sunspot cycle is entered. It is shown, that the latent period of cycles evolution is connected with the next sunspot cycle amplitude and can be used for the prognosis of a level and time of a sunspot maximum. The 24th activity cycle prognosis is made. The found dependences correspond to transport dynamo model of generation of solar cyclicity, it is possible with various speed of meridional circulation. Long-term behavior of extended cycle's lengths and connection with change of a climate of the Earth is considered. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A pictorial comparison of interplanetary magnetic field polarity,solar wind speed,and geomagnetic disturbance index during the sunspot cycle

Observations of interplanetary magnetic field polarity, solar wind speed, and geomagnetic disturbance index (C9) during the years 1962–1975 are compared in a 27-day pictorial format that emphasizes their associated variations during the sunspot cycle. This display accentuates graphically several recently reported features of solar wind streams including the fact that the streams were faster, wider, and longer-lived during 1962–1964 and 1973–1975 in the declining phase of the sunspot cycle than during intervening years (Bame et al., 1976; Gosling et al., 1976). The display reveals strikingly that these high-speed streams were associated with the major, recurrent patterns of geomagnetic activity that are characteristic of the declining phase of the sunspot cycle. Finally, the display shows that during 1962–1975 the association between long-lived solar wind streams and recurrent geomagnetic disturbances was modulated by the annual variation (Burch, 1973) of the response of the geomagnetic field to solar wind conditions. The phase of this annual variation depends on the polarity of the interplanetary magnetic field in the sense that negative sectors of the interplanetary field have their greatest geomagnetic effect in northern hemisphere spring, and positive sectors have their greatest effect in the fall. During 1965–1972 when the solar wind streams were relatively slow (500 km s^-1), the annual variation strongly influenced the visibility of the corresponding geomagnetic disturbance patterns.Visiting Scientist, Kitt Peak National Observatory, Tucson, Arizona.Operated by the Association of Universities for Research in Astronomy, Inc., under contract with the National Science Foundation.     

Fluctuations of Solar Activity during the Declining Phase of the 11-Year Sunspot Cycle

The number of coronal mass ejections (CMEs) erupting from the Sun follows a trend similar to that of sunspot numbers during the rising and maximum phase of the solar cycle. In the declining phase, the CME number has large fluctuations, dissimilar to those of sunspot numbers. In several studies of solar – interplanetary and solar – terrestrial relationships, the sunspot numbers and the 2800-MHz flux (F10) are used as representative of solar activity. In the rising phase, this may be adequate, but in the declining phase, solar parameters such as CMEs may have a different behaviour. Cosmic-ray Forbush decreases may occur even when sunspot activity is low. Therefore, when studying the solar influence on the Earth, one has to consider that although geomagnetic conditions at solar maximum will be disturbed, conditions at solar minimum may not be necessarily quiet.     

Amplitude and period of the dynamo wave and prediction of the solar cycle

The relation of the solar cycle period and its amplitude is a complex problem as there is no direct correlation between these two quantities. Nevertheless, the period of the cycle is of important influence to the Earth's climate, which has been noted by many authors. The present authors make an attempt to analyse the solar indices data taking into account recent developments of the asymptotic theory of the solar dynamo. The use of the WKB method enables us to estimate the amplitude and the period of the cycle versus dynamo wave parameters in the framework of the nonlinear development of the one-dimensional Parker migratory dynamo. These estimates link the period T and the amplitude a with dynamo number D and thickness of the generation layer of the solar convective zone h. As previous authors, we have not revealed any considerable correlation between the above quantities calculated in the usual way. However, we have found some similar dependences with good confidence using running cycle periods. We have noticed statistically significant dependences between the Wolf numbers and the running period of the magnetic cycle, as well as between maximum sunspot number and duration of the phase of growth of each sunspot cycle. The latter one supports asymptotic estimates of the nonlinear dynamo wave suggested earlier. These dependences may be useful for understanding the mechanism of the solar dynamo wave and prediction of the average maximum amplitude of solar cycles. Besides that, we have noted that the maximum amplitude of the cycle and the temporal derivative of the monthly Wolf numbers at the very beginning of the phase of growth of the cycle have high correlation coefficient of order 0.95. The link between Wolf number data and their derivative taken with a time shift enabled us to predict the dynamics of the sunspot activity. For the current cycle 23 this yields Wolf numbers of order 107±7.     

The Relative Phase Asynchronization between Sunspot Numbers and Polar Faculae

The monthly sunspot numbers compiled by Temmer et al. and the monthly polar faculae from observations of the National Astronomical Observatory of Japan, for the interval of March 1954 to March 1996, are used to investigate the phase relationship between polar faculae and sunspot activity for total solar disk and for both hemispheres in solar cycles 19, 20, 21 and 22. We found that (1) the polar faculae begin earlier than sunspot activity, and the phase difference exhibits a consistent behaviour for different hemispheres in each of the solar cycles, implying that this phenomenon should not be regarded as a stochastic fluctuation; (2) the inverse correlation between polar faculae and sunspot numbers is not only a long-term behaviour, but also exists in short time range; (3) the polar faculae show leads of about 50–71 months relative to sunspot numbers, and the phase difference between them varies with solar cycle; (4) the phase difference value in the northern hemisphere differs from that in the southern hemisphere in a solar cycle, which means that phase difference also existed between the two hemispheres. Moreover, the phase difference between the two hemispheres exhibits a periodical behaviour. Our results seem to support the finding of Hiremath (2010).     

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.     

Evolution of Solar and Geomagnetic Activity Indices, and Their Relationship: 1960?�C?2001

We employ annually averaged solar and geomagnetic activity indices for the period 1960??C?2001 to analyze the relationship between different measures of solar activity as well as the relationship between solar activity and various aspects of geomagnetic activity. In particular, to quantify the solar activity we use the sunspot number R _s, group sunspot number R _g, cumulative sunspot area Cum, solar radio flux F10.7, and interplanetary magnetic field strength IMF. For the geomagnetic activity we employ global indices Ap, Dst and Dcx, as well as the regional geomagnetic index RES, specifically estimated for the European region. In the paper we present the relative evolution of these indices and quantify the correlations between them. Variations have been found in: i) time lag between the solar and geomagnetic indices; ii) relative amplitude of the geomagnetic and solar activity peaks; iii) dual-peak distribution in some of solar and geomagnetic indices. The behavior of geomagnetic indices is correlated the best with IMF variations. Interestingly, among geomagnetic indices, RES shows the highest degree of correlation with solar indices.     

Variations of the radiation signatures of noise storm-emitting sunspot groups during a solar cycle

The flux density of a noise storm continuum is known to depend on importance parameters of the associated sunspot group, e.g., its total area A. A study of the continua at 287, 234, 113, and 64 MHz, however, reveals in case of the two cycles Nos. 20 and 21 that the radiation signatures of sunspot groups, with a value of A kept fixed, vary systematically with time indicating regular changes of relevant parameters of the overlying loop systems with the phase of the solar cycle. A trend of intense continua at high frequencies (for definition, cf. Figures 1, 2(b)) to occur preferably during the first activity maximum of a solar cycle has been obtained in either case suggesting a decrease of the emissivity of sunspot groups with time. Vice versa, intense continua at lower frequencies (for definition, cf. Figures 1, 2(b)) were mainly observed during the later phase of both cycles. The latter effect is shown to be attributed rather to a long-term variation of the spectral characteristics of the type-I continua than to an enhanced number of intense type-III continua. From the result obtained it follows that non-potential loops extending to great heights into the corona or developing at least conditions favourable for the generation of an intense type-I continuum even at the frequencies < 100 MHz tend to occur more frequently above sunspot groups during the later phase of a solar cycle than above the comparable groups of its first activity maximum.Furthermore, characteristic periods have been found for both cycles during which the emissivity, especially of the very large sunspot groups, was significantly diminished with reference to the comparable groups of the adjacent time intervals.     

On the use of ‘first spotless day’ as a predictor for sunspot minimum

Defining the first spotless day of a sunspot cycle as the first day without spots relative to sunspot maximum during the decline of the solar cycle, one finds that the timing of that occurrence can be used as a predictor for the occurrence of solar minimum of the following cycle. For cycle 22, the first spotless day occurred in April 1994, based on the International sunspot number index, although other indices (Boulder and American) indicated the first spotless day to have occurred earlier (September 1993). For cycles 9–14, sunspot minimum followed the first spotless day by about 72 months, having a range of 62–82 months; for cycles 15–21, sunspot minimum followed the first spotless day by about 35 months, having a range of 27–40 months. Similarly, the timing of first spotless day relative to sunspot minimum and maximum for the same cycle reveals that it followed minimum (maximum) by about 69 (18) months during cycles 9–14 and by about 90 (44) months during cycles 15–21. Accepting April 1994 as the month of first spotless day occurrence for cycle 22, one finds that it occurred 91 months into the cycle and 57 months following sunspot maximum. Such values indicate that its behavior more closely matches that found for cycles 15–21 rather than for cycles 9–14. Therefore, one infers that sunspot minimum for cycle 23 will occur in about 2–3 years, or about April 1996 to April 1997. Accepting the earlier date of first spotless day occurrence indicates that sunspot minimum for cycle 23 could come several months earlier, perhaps late 1995.The U.S. Government right to retain a non-exclusive, royalty free licence in and to any copyright is acknowledged.     

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.     

太阳黑子月平均面积数活动周及其特征

我们对第12周至第22周的太阳黑子月平均面积数进行统计分析,并与相应的太阳黑子月平均数相比较,结果表明太阳黑子月平均面积数活动周与太阳黑子月平均数活动周有一定的关系。在多数情况下,太阳黑子出现最大值的时间与太阳黑子面积数出现最大值的时间上不一致;太阳黑子平滑月平均数活动周上升期与太阳黑子平滑月平均面积数上升期在大多数情况下不相同;太阳黑子平滑月平均数活动周平均效果的瓦德迈尔效应(Waldmeiereffect)一般要比太阳黑子平滑平均面积数的活动周明显;文中还对太阳黑子平滑月平均面积数活动周的特征进行了分析。