Mathematical modelling of the sunspot cycle

The sunspot record for the time interval 1749–1977 can be represented conveniently by an harmonic model comprising a relatively large number of lines. Solar activity can otherwise be considered as a sequence of partly overlapping events, triggered periodically at intervals of the order of 11 years. Each individual cycle is approximated by a function of the Maxwell distribution type; the resulting impulse model consists of the superposition of the independent pulses. Application of these two models for the prediction of annual values of the Wolf sunspot numbers leads to controversial results. Mathematical modelling of the sunspot time series does not give an unambiguous result.     

Relative sunspot numbers in the first half of eighteenth century

We revised relative sunspot numbers in the time interval 1700–1748 for which Wolf derived their annual means. The frequency of daily observations, counting simultaneously the number of sunspots and the number of sunspot groups necessary for determinating Wolf's relative sunspot numbers, is in this time interval very low and covers, on average, 4.8% of the number of all days only. There also exist incomplete observations not convenient to determine relative sunspot numbers. To enlarge the number of daily relative sunspot numbers we used the nonlinear, two-step interpolation method derived earlier by Letfus (1996, 1999). After interpolation, the mean value increased to 13.8%. Waldmeier (1968) found that the scaling factor k can be derived directly from the observed number of spots f and from the number of sunspot groups g. From the observations made at Zürich (Wolf and his assistants, Wolfer), at Peckeloh, and at Moncalieri during the years 1861–1928, we derived a new, more correct empirical relation. The resulting annual relative sunspot numbers are given in Table II. However, only for 26 years (53.0%) from the total number of 49 years was it possible to derive annual relative sunspot numbers. The observations were missing for the other years. This corresponds with results of Wolf, which gives the annual relative sunspot numbers for all 49 years. For the years when the data were missing, he marked these values as interpolated or very uncertain ones. Most of the observations originate from two data series (Kirch, Plantade), for which Wolf derived a higher scaling factor (k=2.0) than followed from the newly derived relation (k=1.40). The investigated time interval covers four solar cycles. After our results, the height of the first cycle (No. –4), given by Wolf, should be lowered by about two-thirds, the following two cycles (Nos. –3 and –2) lowered by one-third, as given by Wolf, and only the height of the fourth one (No. –1) should be unchanged. The activity levels of the cycles, as represented by group sunspot numbers, are lower by about one-fourth and, in the case of the first one (No. –4) even by two-thirds of the levels derived by us. The group sunspot numbers, derived from a much greater number of observations, have also greater credibility than other estimates. The shapes of the cycles, as given by Wolf, can be considered only as their more or less idealized form.     

Relative Sunspot Numbers in the First Half of the Seventeenth Century

We derived daily relative sunspot numbers and their monthly and annual means in the first half of the seventeenth century. The series of observations collected by Wolf were recorded in the years 1611–1613 and 1642–1644. We used a nonlinear two-step interpolation method derived earlier (Letfus, 1996, 1999) to enlarge the number of daily data. Before interpolation the relative monthly frequency of observations in 24 months of the first time interval 1611–1613 was 49.4% and in 22 months of the second interval 1642–1644 was 49.9%. After interpolation the relative frequency increased in the first time interval to 91.3%, in the second time interval to 82.6%. Most data series in the years 1611–1613 overlap one another and also overlap with a series, for which Wolf estimated a scaling factor converting relative sunspot numbers on the Zürich scale. We derived the scaling factors of all individual series of observations also from the ratios of observed numbers of sunspots to the numbers of sunspot groups (Letfus, 2000). The differences between almost all scaling factors derived in one and the other way are not substantial. All data series were homogenized by application of scaling factors and parallel data in the overlapping parts of data series were averaged. Resulting daily relative sunspot numbers and their monthly and annual means in the years l61l–1613 are given in Table I and those in the years 1642–1644 in Table II. The annual means of these data are compared with analogous data obtained otherwise.     

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.     

Variations in the solar rotation rate derived from Ca+ K plage areas

Daily calcium plage areas for the period 1951–1981 (which include the solar cycle 19 and 20) have been used to derive the rotation period of the Sun at latitude belts 10–15 ° N, 15–20 ° N, 10–15 ° S, and 15–20 ° S and also for the entire visible solar disk. The mean rotation periods derived from 10–20 ° S and N, total active area and sunspot numbers were 27.5, 27.9, and 27.8 days (synodic), respectively. A power spectral analysis of the derived rotation rate as a function of time indicates that the rotation rate in each latitude belt varies over time scales ranging from the solar activity cycle, down to about 2 years. Variations in adjacent latitude belts are in phase, whereas those in different hemispheres are not correlated. The rotation rates derived from sunspot numbers also behave similarly though the dependence over the solar cycle are not very apparent. The total plage areas, integrated over the entire visible hemisphere of the Sun shows a dominant periodicity of 7 years in rotation rate, while the other time scales are also discernible.     

Rotational Modulation of Microwave Solar Flux

Time series data of 10.7 cm solar flux for one solar cycle (1985–1995 years) was processed through autocorrelation. Rotation modulation with varying persistence and period was quite evident. The persistence of modulation seems to have no relation with sunspot numbers. The persistence of modulation is more noticeable during 1985–1986, 1989–1990, and 1990–1991. In other years the modulation is seen, but its persistence is less. The sidereal rotation period varies from 24.07 days to 26.44 days with no systematic relation with sunspot numbers. The results indicate that the solar corona rotates slightly faster than photospheric features. The solar flux was split into two parts, i.e., background emission which remains unaffected by solar rotation and the localized emission which produces the observed rotational modulation. Both these parts show a direct relation with the sunspot numbers. The magnitude of localized emission almost diminishes during the period of low sunspot number, whereas background emission remains at a 33% level even when almost no sunspots may be present. The localized regions appear to shift on the solar surface in heliolongitudes.     

A Fractal Structure of the Time Series of Global Indices of Solar Activity

The structure of time series of daily global indices of solar activity is investigated: the sunspot numbers for the time interval between the years 1854 and 1996, the Greenwich total sunspot area for 1874–1983, the radio-flux at 10.7 cm (F10.7) for 1964–1996, and the Stanford mean solar magnetic field for 1975–1996. The fractal dimensions are determined by two fractal and spectral methods. The identified three time-scale ranges, 2 days–2 months, 2 months–2 years, 2 years–8 and more years, with the fractal dimensions 1.4–1.6, 2, 1.2–1.6, respectively, show perhaps some fractal structure of time series of global indices. The first time-scale range may correspond to ordinary brownian noise and the second to flicker noise. The solar rotation influence of the value of the fractal dimensions at the time range close to the rotational period is studied.     

Relationship between sunspot numbers during years of sunspot maximum and sunspot minimum

A correlation analysis shows that the sunspot numbers at the peaks of the last eight solar cycles are well-correlated with the sunspot numbers in heliolatitudes 20°–40° (specially in the southern hemisphere) occurring in the solar minimum years immediately preceding the solar maximum years.On leave from Physical Research Laboratory, Ahmedabad, India.     

Prediction of the maximum annual mean sunspot number in the coming solar maximum epoch

Using an earlier correlation analysis between the annual sunspot numbers at sunspot maximum epochs and the minimum annual aa index in the immediately preceding years, the minimum annual aa index (21.6) during 1985–86 implies a maximum annual sunspot number of about 190±40 in the coming solar maximum epoch, in about 1988–89.     

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.     

Prediction of the sunspot maximum of solar cycle 23 by extrapolation of spectral components

A simple method MEM-MRA, where spectral peaks are located by MEM (Maximum Entropy Method) and about a dozen most prominent ones are used in MRA (Multiple Regression Analysis) to estimate their amplitudes and phases, was applied to the sunspot number (Rz) series of 1748–1996. Spectral characteristics were different in the successive 3 intervals of 83 years each. Hence, for predictions, only data for the recent 83 years were considered relevant. From the spectra for 1914–1996, the most significant peaks at 5.3, 8.3, 10.5, 12.2, 47 years were used for reconstruction. The match between observed and reconstructed values was good (correlation +0.90). When extrapolated, the reconstructed values indicate a sunspot number maximum for the present solar cycle 23 as 140±9, to occur in year 2000 and for the next solar cycle 24 as 105±9, to occur in year 2010–2011.     

Daily relative sunspot numbers 1749–1848: reconstruction of missing observations

A great part of missing daily relative sunspot numbers in the time interval 1749–1848 was reconstructed by nonlinear two-step method of interpolation. In the first step gaps of missing observations not longer than five days were directly interpolated. In the second step data were sorted to so-called Bartels scheme, i.e., to rows of the length of 27 days subsequently ranged in a matrix. In this step the missing value at any position was interpolated from the data at the same position of preceding and following rows. The interpolation was limited to sequences of no more than four missing data. The procedure enables to interpolate long gaps and simultaneously to respect the 27-day variation of solar activity. Monthly and annual means of relative sunspot numbers are presented. The differences between monthly and annual means of the primary observations and of the data completed by interpolation fluctuate around zero. The amplitude of fluctuations depends inversely on the frequency of observations. Most conspicuous are the deviations in the time interval 1776–1795 where the frequency of observations is very low or almost zero. The average dispersion of monthly differences is ±11.5 R and that of annual differences is ±7.8 R. The two-step method of interpolation was tested on the series of daily data in the time interval 1918–1948. The sequence of missing daily data in the years 1818–1848 represents a masking function. This function was applied to the continuous data series in the time interval 1918–1948 and then the modified series was reconstructed. The differences between the monthly and annual means of primary and reconstructed data are small with fluctuations around zero and with dispersion for monthly differences ±2.7 R a for annual differences ±0.6 R. Corresponding dispersions of the data differences for monthly means in the time interval 1818–1848 are ±4.3 R and ±1.1 R for annual means. The small dispersion values and small differences among them give evidence about the applicability and the effectiveness of the nonlinear two-step method of interpolation and also about high credibility of relative sunspot numbers after reconstruction.     

New Evidence for Long-Term Persistence in the Sun's Activity

Possible persistence of sunspot activity was studied using rescaled range and detrended fluctuation analyses. In addition to actual Wolf numbers (1700–2000 A.D.), two solar proxies were used in this research, viz., an annual sunspot proxy obtained for 1090–1700 A.D. and sunspot numbers reconstructed from the decadal radiocarbon series (8005 B.C. – 1895 A.D). The reconstruction was made using a five-box carbon exchange model. Analyses showed that in all cases the scaling exponent is significantly higher than 0.5 in the range of scales from 25 yr up to 3000 yr. This indicates the existence of a long-term memory in solar activity, in agreement with results obtained for other solar indices.     

Digitization of Sunspot Drawings by Staudacher in 1749 – 1796

Original drawings by J.C. Staudacher made in the period of 1749 – 1796 were digitized. The drawings provide information about the size of the sunspots and are therefore useful for analyses sensitive to sunspot area rather than Wolf numbers. The total sunspot area as a function of time is shown for the observing period. The sunspot areas measured do not support the proposition of a weak, “lost” cycle between cycles 4 and 5. We also evaluate the usefulness of the drawings for the determination of sunspot positions for future studies.     

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.     

Solar neutrino in relation to solar activity

Here we have carried out a power-spectrum analysis of solar nuclear gamma-ray (NGR) flares observed by SMM and HINOTORI satellites. The solar NGR flares show a periodicity of 152 days, confirming the existence of a 152–158 days periodicity in the occurrence of solar activity phenomena and also indicating that the NGR flares are a separate class of solar flares. The power-spectrum analysis of the daily sunspot areas on the Sun for the period 1980–1982 shows a peak around 159 days while sunspot number data do not show any periodicity (Verma and Joshi, 1987). Therefore, only sunspot area data should be treated as an indicator of solar activity and not the daily sunspot number data.     

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.     

Modeling of Sunspot Numbers by a Modified Binary Mixture of Laplace Distribution Functions

This paper presents a new approach for describing the shape of 11-year sunspot cycles by considering the monthly averaged values. This paper also brings out a prediction model based on the analysis of 22 sunspot cycles from the year 1749 onward. It is found that the shape of the sunspot cycles with monthly averaged values can be described by a functional form of modified binary mixture of Laplace density functions, modified suitably by introducing two additional parameters in the standard functional form. The six parameters, namely two locations, two scales, and two area parameters, characterize this model. The nature of the estimated parameters for the sunspot cycles from 1749 onward has been analyzed and finally we arrived at a sufficient set of the parameters for the proposed model. It is seen that this model picks up the sunspot peaks more closely than any other model without losing the match at other places at the same time. The goodness of fit for the proposed model is also computed with the Hathaway – Wilson – Reichmann measure, which shows, on average, that the fitted model passes within 0.47 standard deviations of the actual averaged monthly sunspot numbers.     

Sunspot Plumes and Flow Channels

It is well known that sunspots are dark. This statement is not correct in the sunspot atmosphere between the chromosphere and the corona, where sunspots often are brighter than their surroundings. The brightest feature in the sunspot transition region is called a sunspot plume. Not all sunspots contain a plume. We find that 20 out of 21 sunspots show a plume when one magnetic polarity dominates the sunspot region out to a distance of 50 from the sunspot. Most sunspots show downflows that exceed 25 km s^–1 in the sunspot plumes at temperatures close to 250000 K. This downflow is not maintained by inflow from the corona, but by gas at transition region temperatures, streaming in flow channels from locations well outside the sunspot. We suggest that this inflow is a necessary requirement for the sunspot plume to occur and present a working hypothesis for the origin of sunspot plumes. This paper is the first thorough spectral analysis of sunspot plumes. It is based on simultaneous observations of ten or six EUV emission lines in 42 sunspot regions with the Coronal Diagnostic Spectrometer – CDS on the Solar and Heliospheric Observatory – SOHO. The line profiles are studied in detail with another SOHO instrument, the Solar Ultraviolet Measurements of Emitted Radiation – SUMER.     

Rotation of Plasma in Sunspots

In this paper we present the results of a sunspot rotation study using Abastumani Astrophysical Observatory photoheliogram data for 324 sunspots. The rotation amplitudes vary in theinebreak 2–64° range (with maximum at 12–14°), and the periods around 0–20 days (with maximum atinebreak 4–6 days). It could be concluded that sunspot rotations are rather inhomogeneous and asymmetric, but several types of sunspots are distinguished by their rotational parameters.During solar activity maximum, sunspot average rotation periods and amplitudes slightly increase. This can be affected by the increase of sunspot magnetic flux tube depth. So we can suppose that sunspot formation during solar activity is connected to a rise of magnetic tubes from deeper layers of the solar photosphere, strengthening the processes within the tube and causing variations in rotation.There is a linear relation between tilt-angle oscillation periods and amplitudes, showing higher amplitudes for large periods. The variations of those periods and especially amplitudes have a periodical shape for all types of sunspots and correlate well with the solar activity maxima with a phase delay of about 1–2 years.