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.     

Study of cosmic ray intensity in relation to the interplanetary magnetic field and geomagnetic storms for solar cycle 24

In the present study, we investigate the association of cosmic ray intensity (CRI) with various solar wind parameters (i.e. solar wind speed V, plasma proton temperature, plasma proton density), interplanetary magnetic field (IMF B), geomagnetic storms (GSs), averaged planetary A-index (Ap index) and sun spot number (SSN) for the period 2009–2016 (solar cycle 24) by using their daily mean average. To find the association of CRI with various solar wind parameters, GSs, IMF B, Ap index and SSN, we incorporate the analysis technique by superposed-epoch method. We have observed that CRI decreases with the increase in IMF B. Moreover the time-lag analysis has been performed by the method of correlation coefficient and observed a time lag of 0 to 2 day between the decrease in CRI and increase in IMF B. In addition, we show that the CRI is found to decrease in a similar pattern to disturbance storm time (Dst index) for most of the period of solar cycle 24. The high and positive correlation is found between CRI and Dst index. The CRI and Ap index are better anti-correlated to each other than CRI and IMF. CRI and SSN are positively correlated with each other. Solar wind parameters such as solar wind speed V is a CR-effective parameter while plasma proton temperature and plasma proton density are not CR-effective parameters. The indicated parameters such as Dst index, Ap index, IMF B and solar wind parameters such as solar wind speed V, plasma proton temperature, plasma proton density shows a kind of irregular variations for solar cycle 23 and 24 while CRI and SSN shows distinct behaviour for the two cycle.     

Multi-parametric effect of solar activity on cosmic rays

The long-term modulation of cosmic ray intensity (CRI) by different solar activity (SA) parameters and an inverse correlation between individual SA parameter and CRI is well known. Earlier, it has been suggested that the concept of multi-parametric modulation of CRI may play an important role in the study of long-term modulation of CRI. In the present study, we have tried to investigate the combined effect of a set of two SA parameters in the long-term modulation of CRI. For this purpose, we have used a new statistical technique called “Running multiple correlation method”, based on the “Running cross correlation method”. The running multiple correlation functions among different sets of two SA parameters (e.g., sunspot numbers and solar flux, sunspot numbers and coronal index, sunspot numbers and grouped solar flares, etc.) and CRI have been correlated separately. It is found that the strength of multiple correlation (among two SA parameters and CRI) and cross correlation (between individual SA parameter and CRI) is almost similar throughout the period of investigation (1955–2005). It is also found that the multiple correlations among various SA parameters and CRI is stronger during ascending and descending phases of the solar cycles and it becomes weaker during maxima and minima of the solar cycles, which is in accordance with the linear relationship between SA parameters and CRI. The values of multiple correlation functions among different sets of SA parameters and CRI fall well within the 95% confidence interval. In the view of odd-even hypothesis of solar cycles, the strange behaviour of present cycle 23 (odd cycle), as this is characterized by many peculiarities with double peaks and many quiet periods (Gnevyshev gaps) interrupted the solar activity (for example April 2001, October–November 2003 and January 2005), leads us to speculate that the solar cycle 24 (even cycle) might be of exceptional nature.     

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.     

Study of Geomagnetic Storms and Solar and Interplanetary Parameters for Solar Cycles 22 and 24

The aim of this paper is to investigate the association of geomagnetic storms with the component of the interplanetary magnetic field (IMF) perpendicular to the ecliptic (\(Bz\)), the solar wind speed (\(V\)), the product of solar wind speed and \(Bz\) (VBz), the Kp index, and the sunspot number (SSN) for two consecutive even solar cycles, Solar Cycles 22 (1986?–?1995) and 24 (2009?–?2017). A comparative study has been done using the superposed epoch method (Chree analysis). The results of the present analysis show that \(Bz\) is a geoeffective parameter. The correlation coefficient between Dst and \(Bz\) is found to be 0.8 for both Solar Cycles 22 and 24, which indicates that these two parameters are highly correlated. Statistical relationships between Dst and Kp are established and it is shown that for the two consecutive even solar cycles, Solar Cycles 22 and 24, the patterns are strikingly similar. The correlation coefficient between Dst and Kp is found to be the same for the two solar cycles (?0.8), which clearly indicates that these parameters are well anti-correlated. For the same studied period we found that the SSN does not show any relationship with Dst and Kp, while there exists an inverse relation between Dst and the solar wind speed, with some time lag. We have also found that VBz is a more relevant parameter for the production of geomagnetic storms, as compared to \(V\) and \(Bz\) separately. In addition, we have found that in Solar Cycles 22 and 24 this combined parameter is more relevant during the descending phase as compared to the ascending phase.     

Solar modulation of galactic cosmic rays during 19–23 solar cycles

We make a detailed analysis of cross-correlation and time-lag between monthly data of galactic cosmic rays (GCRs) intensity and different solar activity indices (e.g., sunspot number, sunspot area, green coronal Fe line and 10.7 cm solar radio flux) during 19–23 solar cycles. GCRs time-series data from Kiel neutron monitor station and solar data from the last 50 years period, covering five solar cycles (19–23), and alternating solar polarity states (i.e., five A < 0 and four A > 0) have been investigated. We find a clear asymmetry in the cross-correlation between GCRs and solar activity indicators for both odd and even-numbered solar cycles. The time-lags between GCRs and solar parameters are found different in different solar cycles as well as in the opposite polarity states (A < 0 and A > 0) within the same solar cycle. Possible explanations of the observed results are discussed in light of modulation models, including drift effects.     

Similarities and Distinctions in Cosmic-Ray Modulation During Different Phases of Solar and Magnetic Activity Cycles

We study the solar-activity and solar-polarity dependence of galactic cosmic-ray intensity (CRI) on the solar and heliospheric parameters playing a significant role in solar modulation. We utilize the data for cosmic-ray intensity as measured by neutron monitors, solar activity as measured by sunspot number (SSN), interplanetary plasma/field parameters, solar-wind velocity [V] and magnetic field [B], as well as the tilt of the heliospheric current sheet [Λ], and we analyze these data for Solar Cycles 20?–?24 (1965?–?2011). We divide individual solar cycles into four phases, i.e. low, high, increasing, and decreasing solar activity. We perform regression analysis to calculate and compare the CRI-response to changes in different solar/interplanetary parameters during

1. different phases of solar activity and
2. similar activity phases but different polarity states.

We find that the CRI-response is different during negative (A<0) as compared to positive (A>0) polarity states not only with SSN and Λ but also with B and V. The relative CRI-response to changes in various parameters, in negative (A<0) as compared to positive (A>0) state, is solar-activity dependent; it is ≈?2 to 3 times higher in low solar activity, ≈?1.5 to 2 times higher in moderate (increasing/decreasing) activity, and it is nearly equal in high solar-activity conditions. Although our results can be ascribed to the preferential entry of charged particles via the equatorial/polar regions of the heliosphere as predicted by drift models, these results also suggest that we should look for any polarity-dependent response of solar-wind and transport parameters in modulating CRI in the heliosphere.     

Influence of the 11-Year Solar Cycle on Variations of Cosmic Ray Intensity

The monthly cosmic ray intensity (CRI) time series from Climax, Huancayo, Moscow, Kiel, and Calgary are used to investigate the presence of the 11-year periodic component with special attention paid to the solar influence on these variations. The results show obvious 11-year temporal characteristics in CRI variations. We also find a close anticorrelation between the 11-year solar cycle and CRI variations and time delays of the CRI relative to solar activity.     

On the relationship of cosmic ray intensity with solar,interplanetary, and geophysical parameters

The flux rate of cosmic rays incident on the Earth’s upper atmosphere is modulated by the solar wind and the Earth’s magnetic field. The amount of solar wind is not constant due to changes in solar activity in each solar cycle, and hence the level of cosmic ray modulation varies with solar activity. In this context, we have investigated the variability and the relationship of cosmic ray intensity with solar, interplanetary, and geophysical parameters from January 1982 through December 2008. Simultaneous observations have been made to quantify the exact relationship between the cosmic ray intensity and those parameters during the solar maxima and minima, respectively. It is found that the stronger the interplanetary magnetic field, solar wind plasma velocity, and solar wind plasma temperature, the weaker the cosmic ray intensity. Hence, the lowest cosmic ray intensity has good correlations with simultaneous solar parameters, while the highest cosmic ray intensity does not. Our results show that higher solar activity is responsible for a higher geomagnetic effect and vice versa.     

Mid-term periodicities and heliospheric modulation of coronal index and solar flare index during solar cycles 22–23

We applied fast Fourier transform techniques and Morlet wavelet transform on the time series data of coronal index, solar flare index, and galactic cosmic ray, for the period 1986–2008, in order to investigate the long- and mid-term periodicities including the Rieger (\({\sim }130\) to \({\sim }190\) days), quasi-period (\({\sim }200\) to \({\sim }374\) days), and quasi-biennial periodicities (\({\sim }1.20\) to \({\sim }3.27\) years) during the combined solar cycles 22–23. We emphasize the fact that a lesser number of periodicities are found in the range of low frequencies, while the higher frequencies show a greater number of periodicities. The rotation rates at the base of convection zone have periods for coronal index of \({\sim }1.43\) years and for solar flare index of \({\sim }1.41\) year, and galactic cosmic ray, \({\sim }1.35\) year, during combined solar cycles 22–23. In relation to these two solar parameters (coronal index and solar flare index), for the solar cycles 22–23, we found that galactic cosmic ray modulation at mid cut-off rigidity (\(\hbox {Rc} = 2.43\hbox {GV}\)) is anti-correlated with time-lag of few months.     

Modeling the Relationship Between Neutron Counting Rates and Sunspot Numbers Using the Hysteresis Effect

Several studies show that temporal variations in the Galactic cosmic ray (GCR) intensity display a distinct 11-year periodicity due to solar modulation of the galactic cosmic rays in the heliosphere. The 11-year periodicity of GCRs is inversely proportional to, but out of phase with, the 11-year solar cycle, implying that there is a time lag between actual solar cycle and the GCR intensity, which is known as the hysteresis effect. In this study, we use the hysteresis effect to model the relationship between neutron counting rates (NCRs), an indicator of the GCR intensity, and sunspot numbers (SSNs) over the period that covers the last four solar cycles (20, 21, 22, and 23). Both linear and ellipse models were applied to SSNs during odd and even cycles in order to calculate temporal variations of NCRs. We find that ellipse modeling provides higher correlation coefficients for odd cycles compared to linear models, e.g. 0.97, 0.97, 0.92, and 0.97 compared to 0.69, 0.72, 0.53, and 0.68 for data from McMurdo, Swarthmore, South Pole, and Thule neutron monitors, respectively, during solar cycle 21 with overall improvement of 31 % for odd cycles. When combined to a continuous model, the better correlation observed for the odd cycles increases the overall correlation between observed and modeled NCRs. The new empirical model therefore provides a better representation of the relationship between NCRs and SSNs. A major goal of the ongoing research is to use the new non-linear empirical model to reconstruct SSNs on annual time scales prior to 1610, where we do not have observational records of SSNs, based on changes in NCRs reconstructed from ¹⁰Be in ice cores.     

Solar Modulation of Cosmic Rays during the Declining and Minimum Phases of Solar Cycle 23: Comparison with Past Three Solar Cycles

We study solar modulation of galactic cosmic rays (GCRs) during the deep solar minimum, including the declining phase, of solar cycle 23 and compare the results of this unusual period with the results obtained during similar phases of the previous solar cycles 20, 21, and 22. These periods consist of two epochs each of negative and positive polarities of the heliospheric magnetic field from the north polar region of the Sun. In addition to cosmic-ray data, we utilize simultaneous solar and interplanetary plasma/field data including the tilt angle of the heliospheric current sheet. We study the relation between simultaneous variations in cosmic ray intensity and solar/interplanetary parameters during the declining and the minimum phases of cycle 23. We compare these relations with those obtained for the same phases in the three previous solar cycles. We observe certain peculiar features in cosmic ray modulation during the minimum of solar cycle 23 including the record high GCR intensity. We find, during this unusual minimum, that the correlation of GCR intensity is poor with sunspot number (correlation coefficient R=?0.41), better with interplanetary magnetic field (R=?0.66), still better with solar wind velocity (R=?0.80) and much better with the tilt angle of the heliospheric current sheet (R=?0.92). In our view, it is not the diffusion or the drift alone, but the solar wind convection that is the most likely additional effect responsible for the record high GCR intensity observed during the deep minimum of solar cycle 23.     

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.     

Studies of the correlation between the global radiation,ultraviolet radiation,cosmic radiation,sunspot number and meteorological parameters

The correlation between the ratio of the global irradiation to the extraterrestrial solar radiation (H/H ₀), and the ratio of the ultraviolet solar irradiation to the extraterrestrial solar radiation (H _u /H ₀) on a horizontal surface at Bahrain (=26°), and some terrestrial and solar parameters (the monthly average relative humidity, temperature, relative sunshine duration, cosmic radiation intensity, and sunspot number) have been studied. Moreover, the role of the solar effects and the terrestrial effects on the global and the solar ultraviolet radiation has been studied. A detailed investigation has been carried between the level of the cosmic radiation received at Bahrain and the sunspot number. It was concluded that as the solar activity increases, cosmic radiation and sunspot number play a predominant effect on the correlation of (H/H ₀) and (H _u /H ₀). Furthermore, the correlation between cosmic radiation and sunspot number also increases.     

One Possible Reason for Double-Peaked Maxima in Solar Cycles: Is a Second Maximum of Solar Cycle 24 Expected?

We investigate solar activity by focusing on double maxima in solar cycles and try to estimate the shape of the current solar cycle (Cycle 24) during its maximum. We analyzed data for Solar Cycle 24 by using Learmonth Solar Observatory sunspot-group data collected since 2008. All sunspot groups (SGs) recorded during this time interval were separated into two groups: The first group includes small SGs [A, B, C, and H classes according to the Zurich classification], the second group consists of large SGs [D, E, and F]. We then calculated how many small and large sunspot groups occurred, their sunspot numbers [SSN], and the Zurich numbers [Rz] from their daily mean numbers as observed on the solar disk during a given month. We found that the temporal variations for these three different separations behave similarly. We also analyzed the general shape of solar cycles from Cycle 1 to 23 by using monthly International Sunspot Number [ISSN] data and found that the durations of maxima were about 2.9 years. Finally, we used the ascending time and SSN relationship and found that the maximum of Solar Cycle 24 is expected to occur later than 2011. Thus, we conclude that i) one possible reason for a double maximum in solar cycles is the different behavior of large and small sunspot groups, and ii) a double maximum is expected for Solar Cycle 24.     

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.     

Cosmic-ray Periodicity at 170 Days

Power-spectral analysis of cosmic-ray indices (CRI) data for the years 1989–1991 shows a 170-day periodicity of cosmic rays. The periodicity is related to a strong magnetic field. Power-spectral analysis of the long-term periodicity (11 years) of the CRI data for the years 1953–1997 shows that the period 1989–1991 is a unique one in the sense of the discussed pronounced periodicity. The 170-day periodicity of cosmic rays was interpreted in the base of six solar rotations (1 SR = 28.3-day periodicity of 10.7 cm solar radio flux) and may be connected to the instability of the solar core.     

Evidence for the effect of sunspot activity on the El Niño/Southern Oscillation

The El Niño No. 3 area index (5°S∼ 5°N, 150°W∼ 90°W) and yearly sunspot number (SSN) from 1408 to 1978 are used to investigate the influence of solar activity on the El Niño/Southern Oscillation (ENSO), through periodicity analysis, cross wavelet transform (XWT), cross correlation and ensemble empirical mode decomposition (EEMD) analyses. The solar activity period, the Hale period, and the Gleissberg period are determined in the El Niño index time series, but of weak statistical significance. Cross correlation analysis of the index with SSN, and that of its low-frequency components decomposed by EEMD clearly indicate that solar activity may take effect on the ENSO, and such an impact should undergo an accumulation procedure (phase delay). XWT also indicates the existence of the impact. It is found that the index is negatively correlated with SSN when SSN is large during a certain long-term interval, and positively when SSN is small. Strong El Niño is inferred to be taken place in decade(s) to come.     

Li Be B Energetics and Cosmic Ray Origin

Three different models have been proposed for LiBeB production bycosmic rays: the CRI model in which the cosmic rays areaccelerated out of an ISM of solar composition scaled withmetallicity; the CRS model in which cosmic rays with compositionsimilar to that of the current epoch cosmic rays are acceleratedout of fresh supernova ejecta; and the LECR model in which adistinct low energy component coexists with the postulated cosmicrays of the CRI model. These models are usually distinguished bytheir predictions concerning the evolution of the Be and Babundances. Here we emphasize the energetics which favor the CRSmodel. This model is also favored by observations showing that thebulk (80 to 90%) of all supernovae occur in hot, low densitysuperbubbles, where supernova shocks can accelerate the cosmicrays from supernova ejecta enriched matter. This revised version was published online in July 2006 with corrections to the Cover Date.