A Transfer Function Model for the Sunspot Cycle

The mean annual sunspot record for the time interval 1700–2002 can be considered as a sequence of independent, partly overlapping events, triggered quasi-periodically at intervals of the order of 11 years. The individual cycles are approximated by the step response of a band-pass dynamical system and the resulting model consists of the superposition of the response to the independent pulses. The simulated sunspot data explain 98.4% of the cycle peak height variance and the residual standard deviation is 8.2 mean annual sunspots. An empirical linear relationship is found between the amplitude of the transfer function model for each cycle and the pulse interval of the preceding cycle that can be used as a tool of short-term forecasting of solar activity. A peak height of 112 for the solar cycle 23 occurring in 2000 is predicted, whereas the next cycle would start at about 2007 and will have a maximum around 110 in 2011. Cycle 24 is expected to have an annual mean peak value in the range 95 to 125. The model reproduces the high level of amplitude modulation in the interval 1950–2000 with a decrease afterwards, but the peak values for the cycles 18, 19, 21, and 22 are fairly underestimated. The semi-empirical model also recreates recurring sunspot minima and is linked to the phenomenon of the reversal of the solar magnetic field.     

Periodicities in solar activity

The techniques of power spectral analysis are used to determine significant periodicities in the annual mean relative sunspot numbers. The main conclusion is that a period of 10.45 yr is very basic and can be associated with an excitation of new solar cycles. When combined with a period of 11.8 yr, associated here with the free-running length of a solar cycle, the mean cycle length of 11.06 yr and a phase variation of 190 yr are explained. Similarly the amplitude variations with periods 88 and 59 yr (previously described as the 80-yr cycle) are due to an amplitude modulation of the solar cycle by a period of 11.9±0.3 yr. The results dispute several associations of planetary position and solar activity.Radiophysics Publication RPP 1647, January, 1973.     

Periodicities in Irradiance and in other Solar Activity Indices During Cycle 23

Magnetic fields give rise to distinctive features in different solar atmospheric regimes. To study this, time variations of the flare index, sunspot number and sunspot area, each index arising from different physical conditions, were compared with the solar composite irradiance throughout cycle 23. Rieger-type periodicities in these time series were calculated using Fourier and wavelet transforms (WTs). The peaks of the wavelet power of these periodicities appeared between the years 1999 and 2002. We found that the solar irradiance oscillations are less significant than those in the other indices during this cycle. The irradiance shows non-periodic fluctuations during this time interval. The peaks of the flare index, sunspot number and sunspot total area were seen around 2000.4, 1999.9 and 2001.0, respectively. These periodicities appeared intermittently and were not simultaneous in different solar activity indices during the three years of the maximum phase of solar cycle 23.     

The thermoluminescence profile of a recent sea sediment core and the solar variability

The analysis of the thermoluminescence (TL) profile of the GT14 recent sea sedimentary core shows the existence of four main periodicities of 137.7, 59,12.06, and 10.8 years. Here we discuss the affinity of these waves to the known cycles of solar variability. The beats of the two high frequency components produce a modulated wavetrain with a carrier wave of 11.4 years and an amplitude modulation with period 206 years. The minima of this squared amplitude modulation fall in 1810 and 1913 A.D. and closely correspond to the periods of lowest solar activity as indicated by the sunspot series. The sum of the two low frequency waves can in turn be rewritten as a component with period 82.6 years which is amplitude modulated by a second component with period of 206 years. The 82.6-yr wave has the period commonly attributed to the Gleissberg cycle of solar activity. The maxima of the 82.6-yr wave occur in agreement with the dates of maximum solar radius as suggested by Gilliland (1981).     

A nonlinear solar cycle model with potential for forecasting on a decadal time scale

This paper describes a novel non-linear oscillator model of the sunspot cycle which accurately reproduces several of the observed qualitative and quantitative characteristics of the real cycle including the long term amplitude modulation pattern. The model accounts for 96% of cycle peak height variance over the period 1859 to 1980. The aim of this work is to assess the potential of such models for forecasting solar activity on decadal and possibly longer time scales. Longer term forecasts may have practical economic significance because of the growing evidence for relationships between solar cycle variations and terrestrial weather and climatic variations (Bandeen and Moran, 1975; Currie, 1980; Williams, 1981). The model predicts that cycle 22 will have an annual mean peak amplitude in the range 25 to 45, the lowest peak activity for 260 yr.     

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.     

On Mid-Term Periodicities in Sunspot Groups and Flare Index

In this work we study the mid-term periodicities (MTPs), between 1 and 2 years, of the sunspot groups and the flare index (FI), by separating the data into hemispheres and spectral bands (SBs) according to the most significant periodicities presented by these phenomena. We found that the MTP of sunspot groups has a diminished power during the Modern Minimum and an increased power during the Modern Maximum, with the exception of cycle 20. For flares, the MTP has a diminished power during the low activity cycle 20, and an increased power during cycles 21 and 22. Therefore, for both sunspot groups and FI, cycle 20 shows a very diminished power followed by the active and higher-power cycles 21 and 22; cycle 23 shows a weaker power than cycles 21 and 22. It is uncertain whether MTP can be a precursor of a long-term minimum of solar activity or not, as has been previously suggested. Also, there is no one-to-one correlation between the cycle intensity and the importance of MTP. Concerning the quasi-biennial periodicities and the theory of two kinds of dynamos, we notice the tendency that higher-power cycles mean weaker coupling in the model. Concerning the hemispheric north-south asymmetry, for sunspot groups the southern hemisphere dominates in most of the SBs, while for FI the northern hemisphere dominates for all the SBs. Additionally, the time lag found between the two hemispheres indicates that the degrees of coupling in the photosphere for sunspot groups and in the corona for flares are between moderate and strong. Finally, the modulation shown by the MTP time series suggests that these periodicities are the product of chaotic quasi-periodic processes and not of stochastic processes.     

Reconstruction of Wolf Sunspot Numbers on the Basis of Spectral Characteristics and Estimates of Associated Radio Flux and Solar Wind Parameters for the Last Millennium

A reconstruction of sunspot numbers for the last 1000 years was obtained using a sum of sine waves derived from spectral analysis of the time series of sunspot number R _z for the period 1700–1999. The time series was decomposed in frequency levels using the wavelet transform, and an iterative regression model (ARIST) was used to identify the amplitude and phase of the main periodicities. The 1000-year reconstructed sunspot number reproduces well the great maximums and minimums in solar activity, identified in cosmonuclides variation records, and, specifically, the epochs of the Oort, Wolf, Spörer, Maunder, and Dalton Minimums as well the Medieval and Modern Maximums. The average sunspot number activity in each anomalous period was used in linear equations to obtain estimates of the solar radio flux F _10.7, solar wind velocity, and the southward component of the interplanetary magnetic field.     

Coronal Activity During the 22-Year Solar Magnetic Cycle

The existence of a 22-year heliomagnetic cycle was inferred long ago not only from direct measurements of the solar magnetic field but also from a cyclic variability of a number of the solar activity phenomena. In particular, it was stated (a rule derived after Gnevyshev and Ohl (1948) findings and referenced as the G–O rule in the following) that if sunspot number Rz cycles are organized in pairs of even–odd numbered cycles, then the height of the peak in the curve of the yearly-averaged sunspot numbers Rz-y is always lower for a given even cycle in comparison with the corresponding height of the following odd cycle. Exceptions to this rule are only cycles 4 and 8 which, at the same time, are the nearest even cycles to the limits of the so-called Dalton minimum of solar activity (i.e., the 1795–1823 time interval). In the present paper, we are looking for traces of the mentioned G–O rule in green corona brightness (measured in terms of the Fexiv 530.3 nm emission line intensity), using data covering almost five solar cycles (1943–1994). It was found that the G–O rule seems to work within the green-line corona brightness, namely, when coronal intensity measured in an extended solar middle-latitude zone is considered separately from the rest of the solar surface. On the other hand, the same G–O rule is valid at the photospheric level, as the heliographic latitudinal dependence of sunspot numbers (1947–1984) shows.     

The phase variations of the solar cycle

It has previously been shown that the statistics of the phase fluctuation of the sunspot cycle are compatible with the assumption that the solar magnetic field is generated deep in the Sun by a frequency stable oscillator and that the observed substantial phase fluctuation in the sunspot cycle is due to variation in the time required for the magnetic field to move to the solar surface (Dicke, 1978, 1979). It was shown that the observed phase shifts are strongly correlated with the amplitude of the solar cycle. It is shown here that of two empirical models for the transport of magnetic flux to the surface, the best fit to the data is obtained with a model for which the magnetic flux is carried to the surface by convection with the convection velocity proportional to a function of the solar cycle amplitude. The best fit of this model to the data is obtained for a 12-yr transit time. The period obtained for the solar cycle is T = 22.219 ± 0.032 yr. It is shown that the great solar anomaly of 1760–1800 is most likely real and not due to poor data.     

Seasonal Variation of 5577 å Line Emission At Calcutta and its Variation with Different Solar Parameters

The paper presents the variation of 5577 Å line intensity with relative sunspot number, and 10.7 cm solar flux. The study has obtained the following important results.[(i)] The 5577 Å line intensity at Calcutta is plotted against relative sunspot number, and the variable component of 10.7 cm solar flux during 1984–1985, which is the secondary peak of the descending phase of the 21st solar cycle. The intensity curves show periodic variation with different solar parameters.[(ii)] The 5577 Å line intensity at Mt. Abu also shows periodic variation with solar parameters during the period 1965–1968 when there was a peak phase of the 20th solar cycle.[(iii)] A possible explanation for such variation is also presented.     

The Effect of Relative Sunspot Numbers,Solar Flare Numbers and Variable Component of 10.7 cm Solar Flux on The Seasonal Variation of 6300 Å Line Intensity at Calcutta

The paper presents the seasonal variation of 6300 Å line intensity at Calcutta with relative sunspot number, solar flare number and variable component of 10.7 cm solar flux. A study has been made and important results have been obtained which are as follows. (i) Intensity of 6300 Å line shows periodic variation with relative sunspot number, solar flare number and variable component of 10.7 cm solar flux during the period 1984–1986 which is the secondary peak of the descending phase of 21st solar cycle. (ii) 6300 Å line intensity at Cachoeira Paulista station, taken by Sahai et al. (1988), also shows periodic variation with solar parameters during the period 1978–1980 which is the peak phase of the solar cycle. (iii) A possible explanation of such a type of variation is also presented.     

Short-Term Periodicities in Solar Indices

The purpose of the present communication is to identify the short-term (few tens of months) periodicities of several solar indices (sunspot number, Caii area and K index, Lyman , 2800 MHz radio emission, coronal green-line index, solar magnetic field). The procedure used was: from the 3-month running means (3m) the 37-month running means (37m) were subtracted, and the factor (3m – 37m) was examined for several parameters. For solar indices, considerable fluctuations were seen during the ± 4 years around sunspot maxima of cycles 18–23, and virtually no fluctuations were seen in the ± 2 years around sunspot minima. The spacings between successive peaks were irregular but common for various solar indices. Assuming that there are stationary periodicities, a spectral analysis was carried out which indicated periodicities of months: 5.1–5.7, 6.2–7.0, 7.6–7.9, 8.9–9.6, 10.4–12.0, 12.8–13.4, 14.5–17.5, 22–25, 28 (QBO), 31–36 (QBO), 41–47 (QTO). The periodicities of 1.3 year (15.6 months) and 1.7 years (20.4 months) often mentioned in the literature were seen neither often nor prominently. Other periodicities occurred more often and more prominently. For the open magnetic flux estimated by Wang, Lean, and Sheeley (2000) and Wang and Sheeley (2002), it was noticed that the variations were radically different at different solar latitudes. The open flux for < 45 solar latitudes had variations very similar (parallel) to the sunspot cycle, while open flux for > 45 solar latitudes had variations anti-parallel to the sunspot cycle. The open fluxes, interplanetary magnetic field and cosmic rays, all showed periodicities similar to those of solar indices. Many peaks (but not all) matched, indicating that the open flux for < 45 solar latitudes was at least partially an adequate carrier of the solar characteristics to the interplanetary space and thence for galactic cosmic ray modulation.     

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.     

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.     

Periodicities of solar irradiance and solar activity indices,II

Using a dynamic power spectral analysis technique, the time-varying nature of solar periodicities is investigated for background X-ray flux, 10.7 cm flux, several indices to UV chromospheric flux, total solar irradiance, projected sunspot areas, and a sunspot blocking function. Many prior studies by a host of authors have differed over a wide range on solar periodicities. This investigation was designed to help resolve the differences by examining how periodicities change over time, and how the power spectra of solar data depend on the layer of the solar atmosphere. Using contour diagrams that show the percent of total power over time for periods ranging from 8 to 400 days, the transitory nature of solar periodicities is demonstrated, including periods at 12–14, 26–28, 51–52, and approximately 154 days. Results indicate that indices related to strong magnetic fields show the greatest variation in the number of periodicities, seldom persist for more than three solar rotations, and are highly variable in their frequency and amplitude. Periodicities found in the chromospheric indices are fewer, persist for up to 8–12 solar rotations, and are more stable in their frequency and amplitude. An additional result, found in all indices to varying degrees and related to the combined effects of solar rotation and active region evolution, is the fashion in which periodicities vary from about 20 to 36 days. I conclude that the solar data examined here are both quasi-periodic and quasistationary, with chromospheric indices showing the longest intervals of stationarity, and data representing strong magnetic fields showing the least stationarity. These results may have important implications to the results of linear statistical analysis techniques that assume stationarity, and in the interpretation of time series studies of solar variability.     

The 10-year time delay of the solar cycle luminosity modulation of the Sun behind the solar cycle magnetic oscillation

We found an evidence that the solar cycle luminosity modulation of the Sun deduced from the total irradiance modulation which was measured by the Earth Radiation Budget (ERB) experiment on board of Nimbus 7 from November 16, 1978 to December 13, 1993 was not in phase with the solar cycle magnetic oscillation when we used the sunspot relative number as its index. The modulation was delayed in time behind the solar cycle magnetic oscillation by an amount of about 10.3 years on the order of length of one solar cycle. In order to quantitatively evaluate the correlation between the two quantities, we devised a method to extract characteristics which were proper to a particular solar cycle by defining a new index of the correlation called multiplied correlation index (MCI). We found that the characteristics of the ERB data time profile between solar cycles 21 and 22 were more similar to those of the solar cycle magnetic oscillation between solar cycles 20 and 21 than those between solar cycles 21 and 22 and thus the time profile of the luminosity modulation from the maximum phase of solar cycle 21 to the declining phase of the solar cycle 22 corresponded to the solar cycle magnetic oscillation from the maximum phase of solar cycle 20 to the declining phase of solar cycle 21. We interpret this phenomenon as an evidence that main features of the modulation is not caused by dark sunspots and bright faculae and plages on the surface of the Sun that should instantaneously affect the luminosity modulation but is caused by time-delayed modulation of global convection by the Lorentz force of the magnetic field of the solar cycle. The delay time of about 10.3 years is the time needed for the force to modify the flows of the convection and to modulate heat flow. Thus the delay time is a function of the strength of the magnetic field oscillation of the solar cycle which is represented by amplitude of the solar cycle. Accordingly, the delay time for other time intervals of the solar cycle magnetic oscillation with different amplitudes can be different from 10.3 years for the interval of the present analysis.     

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.     

Predicting Solar Cycle 24 Using a Geomagnetic Precursor Pair

We describe using Ap and F_10.7 as a geomagnetic-precursor pair to predict the amplitude of Solar Cycle 24. The precursor is created by using F_10.7 to remove the direct solar-activity component of Ap. Four peaks are seen in the precursor function during the decline of Solar Cycle 23. A recurrence index that is generated by a local correlation of Ap is then used to determine which peak is the correct precursor. The earliest peak is the most prominent but coincides with high levels of non-recurrent solar activity associated with the intense solar activity of October and November 2003. The second and third peaks coincide with some recurrent activity on the Sun and show that a weak cycle precursor closely following a period of strong solar activity may be difficult to resolve. A fourth peak, which appears in early 2008 and has recurrent activity similar to precursors of earlier solar cycles, appears to be the “true” precursor peak for Solar Cycle 24 and predicts the smallest amplitude for Solar Cycle 24. To determine the timing of peak activity it is noted that the average time between the precursor peak and the following maximum is ≈?6.4 years. Hence, Solar Cycle 24 would peak during 2014. Several effects contribute to the smaller prediction when compared with other geomagnetic-precursor predictions. During Solar Cycle 23 the correlation between sunspot number and F_10.7 shows that F_10.7 is higher than the equivalent sunspot number over most of the cycle, implying that the sunspot number underestimates the solar-activity component described by F_10.7. During 2003 the correlation between aa and Ap shows that aa is 10 % higher than the value predicted from Ap, leading to an overestimate of the aa precursor for that year. However, the most important difference is the lack of recurrent activity in the first three peaks and the presence of significant recurrent activity in the fourth. While the prediction is for an amplitude of Solar Cycle 24 of 65±20 in smoothed sunspot number, a below-average amplitude for Solar Cycle 24, with maximum at 2014.5±0.5, we conclude that Solar Cycle 24 will be no stronger than average and could be much weaker than average.