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.     

A new look at wolf sunspot numbers in the late 1700's

Long-term homogeneous observations of solar activity or many solar cycles are essential for investigating many problems in solar physics and climatology. The one key parameter used in most long-term studies is the Wolf sunspot number, which is susceptible to observer bias, particularly because it is highly sensitive to the observer's ability to see the smallest sunspots. In this paper we show how the Wolf sunspot number can be derived from the number of sunspot groups alone. We utilize this approach to obtain a Group Wolf number. This technique has advantages over the classical method of determining the Wolf number because corrections for observer differences are reduced and long-term self-consistent time series can be developed. The level of activity can be calculated to an accuracy of ± 5% using this method. Applying the technique to Christian Horrebow's observations of solar cycles 1, 2, and 3 (1761–1777), we find that the standard Wolf numbers are nearly homogeneous with sunspot numbers measured from 1875 to 1976 except the peak of solar cycle 2 is too low by 30%. This result suggests that further analyses of early sunspot observations could lead to significant improvements in the uniformity of the measurements of solar activity. Such improvements could have important impacts upon our understanding of long-term variations in solar activity, such as the Gleissberg cycle, or secular variations in the Earth's climate.     

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.     

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.     

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.     

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.     

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.     

Solar activity in the sixteenth and seventeenth centuries (a revision)

Maximum relative sunspot numbers for the 16th and 17th century were computed by means of the dependence of the maximum relative sunspot numbers on the solar cycle rise time and on the cycle asymmetry. In these dependencies four separate modes of relations, two for odd and two for even cycles, were identified. These modes are coupled two and two in even-odd cycle pairs. The rise times and the asymmetries of solar cycles in the 16th and 17th centuries were taken from cycle extreme estimates by Schove (1979), from auroral and telescopic sunspot observations during this period, but with some necessary corrections. Annual relative sunspot numbers and decade averages were estimated from the cycle maxima and the epochs of extremes. In addition, the efficiency of auroral records in latitudes lower than 55 deg was computed for the time interval 1500–1868. For this purpose the dependence of occurrence numbers of aurorae on the cycle and decade means of the relative sunspot numbers was derived.     

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.     

A possible relationship between spectral bands in sunspot number and the space-time organization of our planetary system

For the particular purpose of this paper, Zürich relative sunspot numbers of the time spans 1749–1982, 1749–1865, and 1866–1982, have been analysed anew by two different methos. It is shown that the spectral bands in the power spectra of sunspot numbers between 1 and 234 years obtained from these analyses can be clearly related to the modified configuration frequencies of the giant planets and their harmonics. In particular, the clearly dominant spectral band in sunspot number, the solar cycle of 10.8 years, is given by the configuration period of Jupiter and Saturn (19.859 yr) times the ratio of their distances from the Sun (0.545).     

Sunspot turning-points and aurorae since A.D. 1510

Dates of solar maxima and minima extending back to c. 1610 were estimated by Wolf and Wolfer at Zürich (Waldmeier, 1961) in the nineteenth century, and those back to c. 1710 have been generally accepted. Slight modifications have already been suggested by the author (Schove, 1967) for the seventeenth century, although, in that century, even the existence of the eleven-year cycle has been questioned (Eddy, 1976). In the course of any sunspot cycle we find a pattern of the aurorae in place and time characteristic of sunspot cycles of the particular amplitude-class. These patterns since c. 1710 can be linked to the precise dates of the Zürich turning-points by a set of empirical rules. A sunspot rule is based on the Gnevyshev gap, the gap in large sunspots near the smoothed maximum. These rules are here applied to the period c. 1510–1710 to give improved determination of earlier turning-points, and approximately confirm the dates given for the seventeenth century by Wolfer and for most of the later sixteenth century by Link (1978). Some turning-points for the fifteenth century and revised sunspot numbers for the period 1700–48 are also given.     

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.     

Fitting the sunspot cycles 10–21 by a modified F-distribution density function

It has been found that sunspot cycles 10–21, represented by quarterly mean values of Zürich sunspot number, can be suitably described by the F-distribution density function provided it is modified by introducing five characteristic parameters, in order to achieve an optimal fitting of each cycle. The average cycle calculated from cycles 10–21 has been used as a basis to forecast time and magnitude of the maximum of each cycle, as a function of various numbers of the first quarterly mean values in the beginning N = 8 to 16 quarters. The standard deviations at a 99% significance level calculated from the observed values depend on N, and vary from 1.6 to 1.1 quarters and 65 to 16 units of sunspot number. A rather sufficient forecast is obtained from N = 12 quarters (with inaccuracy of ± 1.5 quarters and ± 24 units); the forecast for cycle 22 yielded, for N = 12, the values t _m = (15.4 ± 1.5) quarters ( 1990.I) and f(t _m ) = (175 ±24 units).     

Duration of Polar Activity Cycles and Their Relation to Sunspot Activity

We have defined the duration of polar magnetic activity as the time interval between two successive polar reversals. The epochs of the polarity reversals of the magnetic field at the poles of the Sun have been determined (1) by the time of the final disappearance of the polar crown filaments and (2) by the time between the two neighbouring reversals of the magnetic dipole configuration (l=1) from the H synoptic charts covering the period 1870–2001. It is shown that the reversals for the magnetic dipole configuration (l=1) occur on an average 3.3±0.5 years after the sunspot minimum according to the H synoptic charts (Table I) and the Stanford magnetograms (Table III). If we set the time of the final disappearance of the polar crown filaments (determined from the latitude migration of filaments) as the criterion for deciding the epoch of the polarity reversal of the polar fields, then the reversal occurs on an average 5.8±0.6 years from sunspot minimum (last column of Table I). We consider this as the most reliable diagnostic for fixing the epoch of reversals, as the final disappearance of the polar crown filaments can be observed without ambiguity. We show that shorter the duration of the polar activity cycle (i.e., the shorter the duration between two neighbouring reversals), the more intense is the next sunspot cycle. We also notice that the duration of polar activity is always more in even solar cycles than in odd cycles whereas the maximum Wolf numbers W _\max is always higher for odd solar cycles than for even cycles. Furthermore, we assume there is a secular change in the duration of the polar cycle. It has decreased by 1.2 times during the last 120 years.     

Sunspot groups at high latitude

Data of sunspot groups at high latitude (35°), from the year 1874 to the present (2000 January), are collected to show their evolutional behaviour and to investigate features of the yearly number of sunspot groups at high latitude. Subsequently, an evolutional pattern of sunspot group number at high latitude is given in this paper. Results obtained show that the number of sunspot groups of a solar cycle at high latitude rises to a maximum value about 1 yr earlier than the time of the maximum of sunspot relative numbers of the solar cycle, and then falls to zero more rapidly. The results also show that, at the moment, solar activity described by the sunspot relative numbers has not yet reached its minimum. In general, sunspot groups at high latitude have not appeared on the solar disc during the last 3 yr of a Wolf solar cycle. The asymmetry of the high latitude sunspot group number of a Wolf solar cycle can reflect the asymmetry of solar activity in the Wolf solar cycle, and it is suggested that one could further use the high latitude sunspot group number during the rising time of a Wolf solar cycle, maximum year included, to judge the asymmetry of solar activity over the whole solar cycle.     

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

Long-Period Cycles of the Sun's Activity Recorded in Direct Solar Data and Proxies

Different records of solar activity (Wolf and group sunspot number, data on cosmogenic isotopes, historic data) were analyzed by means of modern statistical methods, including one especially developed for this purpose. It was confirmed that two long-term variations in solar activity – the cycles of Gleissberg and Suess – can be distinguished at least during the last millennium. The results also show that the century-type cycle of Gleissberg has a wide frequency band with a double structure consisting of 50–80 years and 90–140 year periodicities. The structure of the Suess cycle is less complex showing a variation with a period of 170–260 years. Strong variability in Gleissberg and Suess frequency bands was found in northern hemisphere temperature multiproxy that confirms the existence of a long-term relationship between solar activity and terrestial climate.     

A Comparison of Wolf's Reconstructed Record of Annual Sunspot Number with Schwabe's Observed Record of ‘Clusters of Spots’ for the Interval of 1826–1868

Samuel Heinrich Schwabe, the discoverer of the sunspot cycle, observed the Sun routinely from Dessau, Germany during the interval of 1826–1868, averaging about 290 observing days per year. His yearly counts of ‘clusters of spots’ (or, more correctly, the yearly number of newly appearing sunspot groups) provided a simple means for describing the overt features of the sunspot cycle (i.e., the timing and relative strengths of cycle minimum and maximum). In 1848, Rudolf Wolf, a Swiss astronomer, having become aware of Schwabe's discovery, introduced his now familiar ‘relative sunspot number’ and established an international cadre of observers for monitoring the future behavior of the sunspot cycle and for reconstructing its past behavior (backwards in time to 1818, based on daily sunspot number estimates). While Wolf's reconstruction is complete (without gaps) only from 1849 (hence, the beginning of the modern era), the immediately preceding interval of 1818–1848 is incomplete, being based on an average of 260 observing days per year. In this investigation, Wolf's reconstructed record of annual sunspot number is compared against Schwabe's actual observing record of yearly counts of clusters of spots. The comparison suggests that Wolf may have misplaced (by about 1–2 yr) and underestimated (by about 16 units of sunspot number) the maximum amplitude for cycle 7. If true, then, cycle 7's ascent and descent durations should measure about 5 years each instead of 7 and 3 years, respectively, the extremes of the distributions, and its maximum amplitude should measure about 86 instead of 70. This study also indicates that cycle 9's maximum amplitude is more reliably determined than cycle 8's and that both appear to be of comparable size (about 130 units of sunspot number) rather than being significantly different. Therefore, caution is urged against the indiscriminate use of the pre-modern era sunspot numbers in long-term studies of the sunspot cycle, since such use may lead to specious results.     

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.