A&A 392, 301-307 (2002)
DOI: 10.1051/0004-6361:20020616
K. J. Li^{1,2} - L. S. Zhan ^{1,2,3,4} - J. X. Wang^{2} - X. H. Liu^{5} - H. S. Yun^{6} - S. Y. Xiong^{7} - H. F. Liang^{1,2} - H. Z. Zhao^{1,2}
1 - Yunnan Observatory, YN650011, China
2 - National Astronomical Observatories, CAS, China
3 - Academy of Graduates, CAS, Beijing, China
4 - Department of Physics, Jingdezhen Comprehensive College, Jiangxi, China
5 - Department of Physics, Jiaozuo Educational College, Henan, China
6 - Astronomy Program, SEES, Seoul National University, Seoul, Korea
7 - Library, Yunnan Observatory, Yunnan, China
Received 28 May 2001 / Accepted 18 April 2002
Abstract
In this paper, we investigate the prospect of using previously occurring
sunspot cycle signatures to determine future behavior in an ongoing cycle,
with specific application to cycle 23, the current sunspot cycle. We find
that the gross level of solar activity (i.e., the sum of the total number
of sunspots over the course of a sunspot cycle) associated with cycle 23,
based on a comparison of its first several years of activity against similar
periods of preceding cycles, is such that cycle 23 best compares to cycle 2.
Compared to cycles 2 and 22, respectively, cycle 23 appears 1.08 times larger
and 0.75 times as large. Because cycle 2 was of shorter period, we infer that
cycle 23 also might be of shorter length (period less than 11 years), ending
sometime in late 2006 or early 2007.
Key words: Sun: activity - Sun: sunspots - Sun: general
The well-known 11-year period associated with sunspot cycles was first discovered by Schwabe (1843), then reported again by Wolf (1852) from his examination of the Zurich Observatory sunspot records. Waldmeier (1935, 1939) noted the asymmetry between the rise to maximum and fall to minimum as well as the importance of the amplitude in determining the shape and length of the cycle (Hathaway et al. 1994). Wilson (1988) and Wilson et al. (1998) provides a useful description of other relationships found by examining the historical records of previous sunspot cycle. These findings are useful for predicting of solar activity. However, unfortunately, sunspot cycle has no strictly period. Cycles are observed to vary both in size and length, and periods of inactivity like the Maunder Minimum are known to have occurred (Hathaway et al. 1994). Since the year 1755, the beginning of solar cycle 1, the sunspot cycle varies markedly in amplitude by at least a factor of 4. Also the cycle length (counted from one sunspot minimum to the next) varies appreciably, between approximately 8 and 15 yr (Schrijver & Zwaan 2000). Therefore, prediction of the sunspot number cycle is of great difficulty but deserves to be done due to the extreme importance for many applications, including telecommunications and satellite orbitography (Lantos 2000; Li et al. 2001). Various methods for estimating the relative timing and strength of cycle 23 are examined and summarized by Hathaway et al. (1999, 2001), Li et al. (2001), and Kane (2001). In the present work, investigated is the possibility of using the known sunspot numbers to predict the forthcoming solar activity of a progressing cycle.
Figure 1: Sunspot numbers of cycles 1 to 23. | |
Open with DEXTER |
Next, we want to know whether the set of the two parametric values can be
determined at the early stage of a Part I cycle.
We delete the last one value of the Part I cycle
and repeat the above procedure to get a best-fit result of Part I to Part II
and obtain a set of parametric values.
Then again, we delete the last two values of Part I and
repeat the above procedure to get a best-fit result and obtain a set of
parametric values, and so on until only twelve values of sunspot
numbers left in Part I.
Afterwards, we know the variation of
the expected values of parameters C and T with the points of the used
data of Part I when moving a Part I cycle toward the other Part II cycles.
Figures 2 and 3 show the variation of the determined parameters C and T with the used data points, respectively.
There are 21 lines in each panel of the figures, corresponding to
moving a Part I cycle toward the other 21 complete cycle, called as
Part II cycles here.
Figure 2: Variation of the determined parameter C with the used data points of a Part I cycle when moving the Part I cycle into its corresponding 21 Part II cycles. | |
Open with DEXTER |
Figure 3: Variation of the determined parameters T with the used data points of a Part I cycle when moving the Part I cycle into its corresponding 21 Part II cycles. | |
Open with DEXTER |
Figure 4: Determinations of the parameter C at one-month interval into each cycle of cycles 17 to 22 when moving them respectively into their similar cycles. | |
Open with DEXTER |
Here, in fact, we choose a cycle as the similar cycle of a moved cycle according to the two criteria: the two cycles have a minimum value of the square difference (min ) with each other when the given moved cycle progresses 3 to 4 years into cycle, and meantime, the determined parameters C and T have very small change, variation of parameter C is less than 10, and parameter Tvaries within 4 points (months).
Figures 4 and 5 respectively shows the determinations of parameters Cand T at one-month time interval
into the cycle of Part I when moving the cycle into its similar cycle.
Concerning the time shift T, we have to point out that:
at first we put the series of the sunspot numbers of the moved
cycle at the end of a certain cycle of Part II, at this time the time shift Tis counted as zero, then the given cycle is moved into the certain cycle
of Part II, and the time shift begins to count.
Figures 4 and 5 indicate that the parameters C and T can
actually be quite accurately determined within the first 2 to 3 years
following the start of the given cycle (parameter T can
accurately determined about 0.5 year earlier than parameter C).
This suggests that fitting the behavior of
the given solar cycle with the known solar cycles after
it progresses 3 to 4 years into cycle could find its similar cycle,
and then give a good prediction of the behavior of the given solar cycle over the
remaining 6 to 11 years of the cycle through its similar cycle.
Figure 5: Determinations of the parameter T at one-month interval into each cycle of cycles 17 to 22 when moving them respectively into their similar cycles. | |
Open with DEXTER |
Figure 6: Comparison of the moved cycles (the dashed lines) respectively with their similar cycles (the solid lines) which are multiplied by their determined parametric values C. | |
Open with DEXTER |
Here, we give a new index, the "weight center'' of sunspot numbers of a cycle, which is defined as the time position: the total sunspot number before the time position is equal to that after. It is found that the time shift T can accurately give the distant (D) between two weight centers of the two involved cycles when moving one solar cycle toward the other. For cycles 17 to 22, the error between T and its corresponding Dis not larger than 2 months. For each of cycles 17 to 22, we have calculated the ratio of total sunspot number of its similar cycle to that of its own cycle, resultingly, the ratios for cycles 17 to 22 are 1.10, 1.26, 1.15, 1.87, 1.10, and 1.03, respectively. Values of the parameter C for cycles 17 to 22 are determined to be 1.15, 1.30, 1.17, 1.77, 1.09, and 1.05, respectively, so we know that the parameter C can reflect the ratio of the two involved cycles with error less than 6, it means that we can accurately estimate the gross level (the total sunspot number) of solar activity of a cycle after 3 years following the beginning of the cycle. After understanding the implication of the two parameters, we can easily estimate an approximate range of their possible values: and 0 < C < 8.
Among the 22 complete cycles,
we found that cycle 2 is the similar cycle of cycle 23.
Figures 9 and 10 respectively show the determined parameters Cand T at one-month time interval
into cycle 23 when moving the cycle into its similar cycle.
The gross level of solar activity of cycle 23 is estimated to be 1.08
times as much as cycle 2 and 0.75 times as much as its neighboring cycle,
cycle 22. Comparison of cycle 23 and its similar cycle is shown in Fig. 11,
which indicates that
the similar cycle after multiplied by a constant
can accurately describe the known part of solar activity of
cycle 23, so here we predict that solar activity of cycle 23 will develop
into future along the trace of solar activity of the similar cycle, namely,
the rest part of the similar cycle after March 2001 may be
regarded as the developing tendency of solar cycle 23. MSFC Solar Physics
Branch members Wilson, Hathaway, and Reichmann have studied the sunspot record
for characteristic behavior that might help in predicting future
sunspot activity and given a prediction of sunspot numbers of the
current cycle. Their predicted numbers, which are available through
web site http://www.science.nasa.gov/ssl/pad/solar/predict.htm, are also
plotted in Fig. 11. The figure shows that solar activity predicted by us
fits well to the solar activity predicted by these authors.
Because cycle 2 was of shorter period (Wilson et al. 1996),
we infer here that cycle 23 also might be of shorter length (period less than
11 years), ending sometime in later 2006 or early 2007, well matching our
result (Li et al. 2002).
Figure 7: Variation of the determined parameters C with the used data points of cycle 23 when moving cycle 23 into the 22 complete cycles. | |
Open with DEXTER |
Figure 8: Variation of the determined parameters T with the used data points of cycle 23 when moving cycle 23 into the 22 complete cycles. | |
Open with DEXTER |
Figure 9: Determinations of the parameter C at one-month interval into cycle 23 when moving it into its similar cycle, cycle 2. | |
Open with DEXTER |
Figure 10: Determinations of the parameter T at one-month interval into cycle 23 when moving it into its similar cycle, cycle 2. | |
Open with DEXTER |
Figure 11: Comparison of sunspot numbers of solar cycle 23 (the dashed line), those of its similar cycle (the solid line), and those predicted by MSFC Solar Physics Branch members (the dashed and dotted line). | |
Open with DEXTER |
Up to now, sunspot numbers have been available of 22 complete cycles. They are enough to allow us to utilise these 22 complete cycles, namely 22 "profiles'' of solar activity to predict the forthcoming solar activity of a progressing cycle. On the other hand, as we know, the early record of sunspot sunspots (cycles 1 to 11) is of less accuracy and reliability than those that follow, it means that shapes of the early cycles are somewhat untrue. An analysis on comparisons against these early cycles certainly seems to invite criticism. We believe that, the more complete cycles of large reliability are, the larger is the possibility of using sunspot numbers of the known solar cycles to predict the forthcoming solar activity of a progressing cycle, and the more accurate the predicted solar activity of the progressing cycle should be.
We have not known exactly the physics behind our proposed method. But two empirical principles are justified by the current investigation. First, for a given progressing solar cycle, there are always one previous cycle, whose behavior is mostly similar to the given cycle after its time profile been shifted and/or linearly multiplied by a constant. Secondly, as stated by Hathaway et al. (1994), the first 3 to 4 years behavior of a given solar cycle seems very indicative to the behavior of remainders of the cycle^{*}. These two empirical principles may be used for the predictions of cyclic behavior of other solar and terrestrial activity indices^{}.
Acknowledgements
We would like to thank the referee for very useful and helpful comments. One of authors (Li K. J.) is grateful to the staff of SEES, Seoul National University for their hospitality. This work is supported by the National Science Foundation of China, the State Ministry of Science and Technology of China (NKBRSF G20000784), the Chinese Academy of Sciences, the Science Foundation of Yunnan Province, and the BK21 project of the Korean Government.