Cosmic Ray Intensity Variations in Relation with Solar Flare Index and Sunspot Numbers

From the monthly data of cosmic ray intensity (CRI), sunspot numbers (SSN) and solar flare index (SFI), an attempt has been made to study the relationship between CRI and solar activity (SA) parameters SSN and SFI. The correlation between SA parameters and CRI for different neutron monitoring stations having low, middle and high cut-off rigidity has been investigated. The anti-correlation between SA and CRI is found to exist with some time lag. Based on the method of minimizing correlation coefficient and time-delayed component method, the observed time-lag between SA parameters (SSN and SFI) and CRI has been found to be large for odd solar cycles in comparison to even solar cycles. The results of time-lag analysis between CRI and SSN and between CRI-SFI have also been compared. The findings of correlative study between CRI and SSN are in agreement with earlier results, while the CRI-SFI relationship provides new insights to understand the solar modulation of cosmic rays.     

Long-term modulation of cosmic ray intensity in relation to sunspot numbers and tilt angle

A detailed correlative analysis between sunspot numbers (SSN) and tilt angle (TA) with cosmic ray intensity (CRI) in the neutron monitor energy range has been performed for the solar cycles 21, 22 and 23. It is found that solar activity parameters (SSN and TA) are highly (positive) correlated with each other and have inverse correlation with cosmic ray intensity (CRI). The ‘running cross correlation coefficient’ between cosmic ray intensity and tilt angle has also been calculated and it is found that the correlation is positive during the maxima of odd cycles 21 and 23. Moreover, the time lag analysis between CRI and SSN, and between CRI and TA has also been performed and is supported by hysteresis curves, which are wide for odd cycles and narrow for even cycles.     

Long-Term Modulation of Cosmic-Ray Intensity and Correlation Analysis Using Solar and Heliospheric Parameters

Based on the monthly sunspot numbers (SSNs), the solar-flare index (SFI), grouped solar flares (GSFs), the tilt angle of heliospheric current sheet (HCS), and cosmic-ray intensity (CRI) for Solar Cycles 21?–?24, a detailed correlation study has been performed using the cycle-wise average correlation (with and without time lag) method as well as by the “running cross-correlation” method. It is found that the slope of regression lines between SSN and SFI, as well as between SSN and GSF, is continuously decreasing from Solar Cycle 21 to 24. The length of regression lines has significantly decreased during Cycles 23 and 24 in comparison to Cycles 21 and 22. The cross-correlation coefficient (without time lag) between SSN–CRI, SFI–CRI, and GSF–CRI has been found to be almost the same during Cycles 21 and 22, while during Cycles 23 and 24 it is significantly higher between SSN–CRI and HCS–CRI than for SFI–CRI and GSF–CRI. Considering time lags of 1 to 20 months, the maximum correlation coefficient (negative) amongst all of the sets of solar parameters is observed with almost the same time lags during Cycles 21?–?23, whereas exceptional behaviour of the time lag has been observed during Cycle 24, as the correlation coefficient attains its maximum value with two time lags (four and ten months) in the case of the SSN–CRI relationship. A remarkably large time lag (22 months) between HCS and CRI has been observed during the odd-numbered Cycle 21, whereas during another odd cycle, Cycle 23, the lag is small (nine months) in comparison to that for other solar/flare parameters (13?–?15 months). On the other hand, the time lag between SSN–CRI and HCS–CRI has been found to be almost the same during even-numbered Solar Cycles 22 and 24. A similar analysis has been performed between SFI and CRI, and it is found that the correlation coefficient is maximum at zero time lag during the present solar cycle. The GSFs have shown better maximum correlation with CRI as compared to SFI during Cycles 21 to 23, indicating that GSF could also be used as a significant solar parameter to study the cosmic-ray modulation. Furthermore, the running cross-correlation coefficient between SSN–CRI and HCS–CRI, as well as between solar-flare activity parameters (SFI and GSF) and CRI is observed to be strong during the ascending and descending phases of solar cycles. The level of cosmic-ray modulation during the period of investigation shows the appropriateness of different parameters in different cycles, and even during the different phases of a particular solar cycle. We have also studied the galactic cosmic-ray modulation in relation to combined solar and heliospheric parameters using the empirical model suggested by Paouris et al. (Solar Phys.280, 255, 2012). The proposed model for the calculation of the modulated cosmic-ray intensity obtained from the combination of solar and heliospheric parameter gives a very satisfactory value of standard deviation as well as \(R^{2}\) (the coefficient of determination) for Solar Cycles 21?–?24.     

Aspects of the long-term cosmic-ray modulation

To investigate the long-term modulation of galactic cosmic rays at the ground-based detector energies, the monthly values of the neutron monitor (Climax, Mt. Washington, Deep River, and Huancayo) and ionization chamber (Cheltenham/Fredericksburg, Huancayo, and Yakutsk) intensities have been correlated with the sunspot numbers (used as a proxy index for transient solar activity) for each phase of sunspot cycles 18 to 22. Systematic differences are found for results concerning odd and even sunspot cycles. During odd cycles (19 and 21) the onset time of cosmic-ray modulation is delayed when compared with the onset time of the sunspot cycle, while they are more similar during even (18, 20, and 22) cycles. Checking the green corona data, on a half-year basis, we found typical heliolatitudinal differences during ascending phases of consecutive sunspot cycles. This finding suggests a significant role of the latitudinal coronal behaviour in the heliospherical dynamics during a Hale cycle. Such effectiveness concerns not only the transient interplanetary perturbations but also the recurrent ones. In fact, when lag between cosmic-ray data and sunspot numbers is considered, the anticorrelation between both parameters is very high (correlation coefficient |r| > 0.9) for all the phases considered, except for the declining ones of cycles 20 and 21, when high-speed solar wind streams coming from coronal holes affect the cosmic-ray propagation, and theRz parameter is no longer the right proxy index for solar-induced effects in the interplanetary medium.     

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

Prediction for the amplitude of solar cycle 24 from the pattern of activity near the cycle minimum

Correlations are investigated between the pattern of solar activity described by the smoothed monthly relative sunspot numbers (Wolf numbers) near the minimum of a solar cycle and the cycle amplitude. The closest correlation is found between the amplitude of a solar cycle and the sum of the decrease in activity over two years prior to the cycle minimum and the increase in activity over two years after the minimum; the correlation coefficient between these parameters is 0.92. This parameter is used as a precursor to predict the amplitude of solar cycle 24, which is expected to reach its maximum amplitude (85 ± 12) in February 2014. Based on the correlations between the mean parameters of solar cycles, cycle 24 is expected to last for approximately 11.3 years and the minimum of the next cycle 25 is predicted for May 2020.     

Analysis of Sunspot Activity Cycles

A nonlinear analysis of the daily sunspot number for each of cycles 10 to 23 is used to indicate whether the convective turbulence is stochastic or chaotic. There is a short review of recent papers considering sunspot statistics and solar activity cycles. The differences in the three possible regimes – deterministic laminar flow, chaotic flow, and stochastic flow – are discussed. The length of data sets necessary to analyze the regimes is investigated. Chaos is described and a chronology of recent results that utilize chaos and fractals to analyze sunspot numbers follows. The parameters necessary to describe chaos – time lag, phase space, embedding dimension, local dimension, correlation dimension, and the Lyapunov exponents – are determined for the attractor for each cycle. Assuming the laminar regime is unlikely if chaos is not indicated in a cycle by the calculations, the regime must be stochastic. The sunspot numbers in each of cycles 10 to 19 indicate stochastic behavior. There is a transition from stochastic to chaotic behavior of the sunspot numbers in cycles 20, 21, 22, and 23. These changes in cycles 20 – 23 may indicate a change in the scale of turbulence in the convection zone that could result in a change in the convective heat transfer and a change in the size of the convection region for these four cycles.     

Hysteresis,time lag,and relation between solar activity and cosmic rays during solar cycle 24

We study galactic cosmic ray (GCR) modulation during solar cycle 24. For this study we utilize neutron monitor data together with sunspot number (SSN) and 10.7 cm solar radio flux (SRF) data. We plot hysteresis curve between the GCR intensity and SSN, and GCR intensity and SRF. We performed time-lag correlation analysis to determine the time lag between GCR intensity and solar activity parameters. The time lag is determined not only for the whole solar cycle, but also during the two polarity states of the heliosphere (A<0 and A>0) in solar cycle 24. We notice differences in time lags during two polarity epochs of the solar cycle. We discuss these differences in the light of existing modulation models. We compare the results of this very weak solar activity cycle with the corresponding results reported for the previous comparatively more active solar cycles.     

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.     

Area of polar coronal holes and sunspot activity: Years 1939–1993

The correlation between the size of polar coronal holes and sunspot numbers has been investigated for the last five solar cycles. The area of polar coronal holes over the period from 1939 to 1993 was derived from ground-based observations of the green coronal line at 530.3 nm (Fe xiv). Correlation analysis revealed that there is no general shift in the maxima of the curves of these two solar indices. The analysis showed the same shift in months in cycles 21 and 22 when the best correlation between the indices is reached; the time shift found in cycle 20 is slightly different from that in cycle 18; in cycle 19, there is found a shift with a value between the values in cycles 18, 20 and 21, 22. The time between succesive peaks of smoothed polar hole size and smoothed sunspot number is different in each cycle.     

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.     

Longterm Prediction of Solar Activity Using the Combined Method

The Combined Method is a non-parametric regression technique for long-term prediction of smoothed monthly sunspot numbers. Starting from a solar minimum, a prediction of the succeeding maximum is obtained by using a dynamo-based relation between the geomagnetic aa index and succeeding solar maxima. Then a series of predictions is calculated by computing the weighted average of past cycles of similar level. This technique leads to a good prediction performance, particularly in the ascending phase of the solar cycle where purely statistical methods tend to be inaccurate. For cycle 23 the combined method predicts a maximum of 160 (in terms of smoothed sunspot number) early in the year 2000.     

Properties of the Current 23rd Solar-Activity Cycle

The main properties of the current cycle match almost completely those of average-magnitude solar cycles, and some of the features of the current cycle may indicate a change in the generation mode of magnetic fields in the solar convection zone. In this case, the Sun enters a period of intermediate and weak cycles of solar activity (SA) in terms of the Wolf numbers, which may last for 50 to 100 years. This change may result in further pollution of the Earth's environment (near-Earth space) due to the unfavorable regime of removing cosmic garbage from low-Earth orbit, the substantial increase of the radiation background in near space (the weakening of interplanetary magnetic fields will result in an increased concentration of galactic cosmic rays in the heliosphere), and other, possibly unfavorable, consequences. The main development stages of the 23rd solar-activity cycle are the following: the minimum of the 22nd solar cycle, May 1996 (W*=8.0); the beginning of the growth phase, September 1997; the maximum of the smoothed relative sunspot number, April, 2000; the global polarity reversal of the general solar magnetic field, July to December 2000; the secondary maximum of the relative sunspot number, November 2001; the maximum of the 10.7-cm radio flux, February 2002; the phase of the cycle maximum, October 1999 to June 2002; the beginning of the decrease phase, July 2002; the most powerful flare events of the current cycle, October to November 2003; and the likely point of minimum of the current SA cycle, November to December 2006.     

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.     

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

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

Galactic Cosmic Ray Modulation and the Last Solar Minimum

In this work the galactic cosmic ray modulation in relation to solar activity indices and heliospheric parameters during the years 1996??C?2010 covering solar cycle 23 and the solar minimum between cycles 23 and 24 is studied. A new perspective of this contribution is that cosmic ray data with a rigidity of 10 GV at the top of the atmosphere obtained from many ground-based neutron monitors were used. The proposed empirical relation gave much better results than those in previous works concerning the hysteresis effect. The proposed models obtained from a combination of solar activity indices and heliospheric parameters give a standard deviation <?10?% for all the cases. The correlation coefficient between the cosmic ray variations of 10?GV and the sunspot number reached a value of r=?0.89 with a time lag of 13.6±0.4 months. The best reproduction of the cosmic ray intensity is obtained by taking into account solar and interplanetary indices such as sunspot number, interplanetary magnetic field, CME index, and heliospheric current sheet tilt. The standard deviation between the observed and calculated values is about 7.15?% for all of solar cycle 23; it also works very well during the different phases of the cycle. Moreover, the use of the cosmic ray intensity of 10?GV during the long minimum period between cycles 23 and 24 is of special interest and is discussed in terms of cosmic ray intensity modulation.     

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.     

Long-Term Sunspot Number Prediction based on EMD Analysis and AR Model

The Empirical Mode Decomposition (EMD) and Auto-Regressive model (AR) are applied to a long-term prediction of sunspot numbers. With the sample data of sunspot numbers from 1848 to 1992, the method is evaluated by examining the measured data of the solar cycle 23 with the prediction: different time scale components are obtained by the EMD method and multi-step predicted values are combined to reconstruct the sunspot number time series. The result is remarkably good in comparison to the predictions made by the solar dynamo and precursor approaches for cycle 23. Sunspot numbers of the coming solar cycle 24 are obtained with the data from 1848 to 2007, the maximum amplitude of the next solar cycle is predicted to be about 112 in 2011-2012.     

INTER-CYCLE VARIATIONS OF SOLAR IRRADIANCE: SUNSPOT AREAS AS A POINTER

Most of the present models and reconstructions of solar irradiance use the concept of Photometric Sunspot Index (PSI) to account for the influence of sunspots on solar brightness. Since PSI is based on measured sunspot areas a firm database of such areas is essential. We show, however, that a significant disagreement exists between the data provided by the Royal Greenwich Observatory (from 1874 to 1976) and newer measurements provided by the observatories of Rome, Yunnan, Catania, and the US Air Force. The overlap of the time intervals over which sunspot areas were measured at Greenwich and Rome allows us to quantify the difference between the Greenwich and other data sets. We find that the various data sets differ, at least in a statistical sense, mainly by a correction factor of between 1.15 and 1.25.The revised time series of sunspot areas correlates well with the Zürich sunspot relative numbers over the last 120 years, with the relationship between sunspot areas and sunspot numbers changing only slightly from one cycle to the next. In particular, no indication exists for any extraordinary magnetic behavior of the Sun during the last 2 decades, as might falsely be concluded if the various sunspot area data sets are uncritically combined. There are, however, some indications that cycles 15 and 16 deviate from the rest. We expect that our results should have a significant influence on the reconstruction of the historical solar irradiance.     

Do polar faculae on the sun predict a sunspot cycle?

The paper reports the results of the analysis of the data on polar faculae for three solar cycles (1960–1986) at the Kislovodsk Station of the Pulkovo Observatory and on polar bright points in Ca ii K line for two solar cycles (1940–1957) at the Kodaikanal Station of the Indian Institute of Astrophysics. We have noticed that the monthly numbers of polar faculae and polar bright points in Ca ii K line and monthly sunspot areas in each hemisphere of the following solar cycle have a correlation with each other. A new cycle of polar faculae and polar bright points in the Ca ii K line begins after the polar magnetic field reversal. We find that the smaller the period between the ending of the polar field reversal and the beginning of a new sunspot cycle is, the more intense is the cycle itself. The intensity of the forthcoming solar cycle (cycle 22) and the periods of strong fluctuations in activity expected in this cycle are also discussed.