Experimental Investigation of the Delay Time in Galactic Cosmic Ray Flux in Different Epochs of Solar Magnetic Cycles: 1959 – 2014

In the present article, we analyze long-term changes in the intensity of galactic cosmic rays (GCRs) in different polarity epochs of the solar magnetic cycles from 1959 to 2014. Our purpose is to carry out a study of the delay time (DT) between the changes of the GCR intensity and various parameters characterizing the conditions in the heliosphere. We prove the existence of varying DTs between the changes of GCR intensity and the parameters characterizing solar activity, such as sunspot number and tilt angle. Based on our investigation, we obtained different DTs in epochs with different global solar magnetic field polarities. We conclude that the observed DTs are very important parameters for the study of GCR transport in the heliosphere.


Introduction
The 11-year variation in the intensity of galactic cosmic rays (GCRs) is generally associated with a similar variation of the solar activity (SA); there is an anticorrelation between noteworthy. The increase of the TA to the maximum value in each of these three cycles is remarkably similar, while the decay of the TA to low values is more gradual in odd cycles. Inceoglu et al. (2014) used the hysteresis effect to model the relationship between neutron counting rates (NCRs), an indicator of the GCR intensity, and SSN over the period that covers the four Solar Cycles 20, 21, 22, and 23. On the basis of a cross-correlation coefficient between monthly NCRs and SSN, time lags of about 5 -6 months were established and calculated for NCRs in the period from 1957 to 2012. Singh, Saxena, and Tiwari (2018) analyzed the Solar Cycles: 22 and 23 (1986 -2008). In their article, the SA is characterized by the coronal index (CI) and the solar flare index (SFI). These authors used the fast Fourier transform and wavelet techniques. Sierra-Porta (2018) studied Solar Cycles 21, 22, and23 (1976 -2008). GCR modulation is controlled by SA, which depends on the conditions of the heliosphere. In the latter article, SA is represented by three parameters: SSN, the flare index (FI), and the Ap index. The author used a cross-correlation analysis to find the maximum negative correlation needed to find DTs. Chowdhury, Kudela, and Dwivedi (2013) and Chowdhury and Kudela (2018) stated that there is no DT between the GCR intensity and the IMF magnitude, nor between the GCR intensity and the Ap index. The above considerations show that DTs are a complex subject that depends on many parameters. DTs are influenced, not only by particle drifts, but also by the structure of the magnetic turbulence at various periods of SA characterized by many parameters, which determine the character of diffusion and of coronal mass ejection (CME) activity.

Experimental Results
Unlike the case of previous works dealing with the DT problem, the purpose of this article is to carry out a study of the DTs between changes of GCR intensity and various parameters characterizing the conditions in the heliosphere such as SSN, TA, and IMF intensity, B. From 1959 -2014 we took monthly data of the GCR flux from the neutron monitor (NM) at Kiel and SSN, TA, and IMF magnitude. This time period is divided into five 11-year periods, each spanning a period between the SA maximum and minimum with different polarities (positive or negative) of the solar magnetic cycle (SMC).
These five 11-year periods are further divided into ten sub-periods for analysis in five different cases. In Figure 1, we present time profiles of the inverted intensity of the GCR flux and from the NM at Kiel in the 1959 -2014, divided into five sub-periods corresponding to different polarities of the SMC, and the SSN.
The correlation coefficient, r, between the GCR intensity, I (GCR), and each of the SA parameters, SSN, TA, and B, was calculated to examine the DTs. This was done for each of the five 11-year time periods, each one having a different GSMF polarity. In Figures 2 -6, time profiles of monthly data of the SSN and inverted intensity of GCR from the Kiel NM in five periods are presented.

Methodology
The correlation coefficient, r, between monthly data of the GCR flux and each of the SA parameters, SSN, TA, and B, was calculated for the considered periods, with various monthly shifts: t = 0, 1, 2, . . . , 18 months, to study the DT.
For instance, we calculated the correlation coefficient, for two independent series of input data of SSN and GCR. Next, I (GCR) data was shifted in 0, 1, 2, . . . , 18 monthly steps relative to the SSN input data and again we calculated r after every step. In each of the five periods, the DT was determined at each step when the correlation coefficient, r, reached the maximum value, r SSN, I i (GCR) → max DT, where i = 0, 1, 2, 3, . . . , 18 is the number of shifted months.
In Table 1, and Figures 7a, and 7b, we present sample correlation coefficients r(SSN, I i (GCR)) for the periods 1959 -1969 (A < 0) and 1991 -1999 (A > 0). We used the same method for the parameters I (GCR), SSN, TA, and B (see Table 2). Table 2 presents the correlation coefficients and DTs in five different periods.
In Table 2 and Figures 2 -6, one can see that the DTs between SSN and I (GCR) and between SSN and B are larger in A < 0 periods than in A > 0 periods. At the same time, we do not observe significant DTs between B and I (GCR). We also see significant correlations between the remaining parameters with various DTs; however, we do not observe a clear dependence of the DTs on the SMC polarity (A > 0 and A < 0). For detailed examination of the DTs, each of the five periods with a given polarity of SMC has been divided into two sub-periods: ascending and descending (marked with an upward and downward arrow, respectively) regarding the GCR intensity. Table 3 shows r and DT between each of the six pairs calculated for these ten sub-periods.
In Table 3, we observe a 22-year periodicity of the DT in the correlation between SNN and I (GCR), but we do not observe such periodicity in the correlation between other parameters, despite the high correlation between them. On the contrary, from DT (SSN, I (GCR)) in Table 3, one can see that in A < 0, SMC is larger in the GCR ascending period than in the descending one. For SMC where A > 0, there is no significant DT observed in the ascending  period, but a small DT is observed in the descending period. In the next section, we discuss these results.

Discussion and Interpretation of the Results
Analyzing the correlations between changes of the GCR intensity and various parameters (see Table 3) representing the conditions in the heliosphere, we find that there is a certain regularity between SSN changes and I (GCR); particularly, we find the 22-year variation of the DT, while we do not recognize anything similar in the correlations between other parameters. Each cycle of SA contains two periods with different polarity of the GSMF. In Cycles 20, 22, and 24, we observe a small DT (one-five months) in periods with A < 0 and with descending GCR intensity, while no significant DT is observed in periods with A > 0 and ascending GCR intensity, which consequently results in small DTs in even SA cycles. In turn, in the odd SA Cycles 19, 21, and 23, we observe a large DT in periods with A < 0 and with ascending I (GCR), while in periods with A > 0 and descending I (GCR), we observe a small DTs (one-five months), which consequently results in large DTs in odd SA cycles. The DT between I (GCR) and each of the other parameters is very important regarding different periods of SA. The mechanism for the transport of GCRs in the heliosphere is governed by four major processes: diffusion, convection, adiabatic energy changes, and drift caused by the gradient and curvature of IMF. Spatial diffusion is caused by pitch angle scattering of GCRs by IMF turbulence. All these processes are vital, but their contribution in the modulation is different and depends on the SA. Many studies of the DT problem showed that the DTs are more pronounced in odd SA cycles (Burlaga, McDonald, and Ness, 1993;Parker, 1963) than in even SA cycles (Dorman, 1959;Jokipii and Thomas, 1981) and interpreted this phenomenon in terms of drift effects (Cliver and Ling, 2001;Chowdhury, Kudela, and Dwivedi, 2013;Chowdhury and Kudela, 2018;Mavromichalaki, Belehaki, and Rafios, 1998).
Based on the drift theory of Jokipii (1971) of the modulation of GCR in the A > 0 epoch, a drift stream of GCRs caused by the gradient and curvature of the IMF preferentially enters the heliosphere from the polar region and is ejected outward along the equatorial current sheet.
The reverse situation occurs in the A < 0 epoch, where the drift stream of GCRs enters the heliosphere along the HCS and exits the polar region. In the A < 0 epoch, GCR particles measured at the Earth are more affected by propagating diffusive barriers associated with SA and the waviness of HCS, resulting in large DTs.
The opposite situation is observed in the A > 0 epoch. GCR particles measured at the Earth are less affected by propagating diffusive barriers which causes small or insignificant DTs. Cliver and Ling (2001) showed 22-year patterns in the relationship of SSN with TA and I (GCR). Based on the analysis of the correlations, the authors noted that in odd cycles (19 and 21), DTs between SSN and I (GCR) and between TA and I (GCR) are ≈ 13 and ≈ 15 months, respectively. During even cycles, the correlations between SSN and I (GCR) and TA and I (GCR) are essentially similar. Chowdhury, Kudela, and Dwivedi (2013) noted that the DT between I (GCR) and some solar indices, such as SSN or area and SFI, in Cycle 23 is remarkably large, reaching ≈ 18 months. Previous investigations (e.g. Nagashima and Morishita, 1980a;Mavromichalaki, Belehaki, and Rafios, 1998;Mavromichalaki, Paouris, and Karalidi, 2007) of Cycles 17 to 22 showed an average DT of ≈ 2.4 months in even cycles and of ≈ 12.4 months for odd ones. According to the interpretation of 22-year patterns in the correlations between SSN and I (GCR) and TA and I (GCR) found by Cliver and Ling (2001) incoming GCRs during the A > 0 epoch will be less affected by the drift effects due to increasing of TA at the beginning of an odd cycle or due to the diffusion associated with enhanced coronal mass ejection (CME) activity. At the beginning of even SA cycles for A < 0, when GCRs approach the Sun along the HCS, they will be easily affected by the changes in the TA and low-latitude CMEs. Also, the interpretation by Van Allen (2000) is consistent with the drift effect, although the author rightly emphasizes the importance of the irregular fields at high latitudes (e.g. Balogh et al., 1999) and the small latitudinal gradients (Simpson et al., 1995) measured by Ulysses.   Table 3 Correlation coefficients, r, and DTs for 10 sub-periods. In the cycles column, upward and downward arrows represent ascending and descending GCR intensity, respectively (see text and Table 2).
No Period  Chowdhury, Kudela, and Dwivedi (2013) and Chowdhury and Kudela (2018) discussed in detail the phenomena of GCR modulation and the DT problem between I (GCR) and various parameters characterizing conditions in the heliosphere. Chowdhury, Kudela, and Dwivedi (2013) mentioned that this is a very complex phenomenon occurring throughout the heliosphere that depends on several factors. No single solar parameter can account for the GCR intensity variations. The IMF magnitude, B, plays a vital role in GCR modulation, because the Larmor radius (an important parameter determining particle transport in space) is inversely proportional to the strength of the IMF. An increase in the IMF strength decreases the Larmor radius and the diffusion coefficient, leading to an increase in GCR modulation. It is remarkable, however, that the highest correlation between B and I (GCR) occurs at zero DT in Cycle 23 (see DT (B, I (GCR)) in Table 2). Chowdhury, Kudela, and Dwivedi (2013) proposed that there is no significant DT between B and I (GCR), because the local disturbances such as CMEs, traveling shocks, etc., which are injected into the inner heliosphere, dominate over the effects of merged interaction regions (MIRs) and globally merged interaction regions (GMIRs) operating at large distances in the outer heliosphere in this interval.
It is possible that some indices, such as SSN and SFI, represent global effects, whereas others, such as B and Ap, represent local effects, limited to near the Earth on the ecliptic plane. This is consistent with Usoskin's et al. (1998) suggestion that the GCR modulation appears to be clearly correlated with only the global indices because of the complicated transport of GCRs in the heliosphere. Cane et al. (1999) reported that the physical mechanisms responsible for the 22-year modulation of GCRs are time-dependent heliospheric drifts and outward propagating diffusive barriers. The relative importance of these two, fundamentally different, mechanisms changes with the SMC. Diffusive barriers are supposedly formed by CMEs, shocks, and high-speed flows at 10 -15 AU, the above-mentioned MIRs. Drifts are expected to be most influential in the years when SA is lowest and the large-scale heliospheric magnetic field changes slowly. In such a case, diffusive barriers or GMIRs tend to be absent Cane et al. (1999).
In the present article, we examine the whole period from 1959 to 2014 which includes two and a half 22-year solar magnetic cycles, and we divided it into five periods according to their polarity (see Table 2). To examine the DTs in a more detailed way, each of these five periods were further divided into two sub-periods: ascending and descending periods of I (GCR) (see Table 3). We observed the 22-year periodicity of the DT between SSN and I (GCR) and SSN and B but not between the remaining parameters (see Table 2). In particular, we observed the 22-year periodicity of the DT only between SSN and I (GCR), but we did not observe a significant DT between B and I (GCR), despite the high correlation between them (see Table 3). We think that the main reason for this phenomenon is essentially the temporal rearrangement of the structure of the IMF turbulence from the minimum to the maximum epoch of SA, which causes the structure of the IMF turbulence to show polarity dependence .

Conclusions
1. Solar activity characterized by the global sunspot number (SSN) and its change determines the conditions prevailing in the heliosphere, i.e. variable IMF and its turbulence, HCS, CME, MIR, shocks. These changing conditions affect the GCR propagation in the heliosphere. 2. The structure of the IMF turbulence significantly changes from the minimum to the maximum epoch of SA and shows polarity dependence. 3. A 22-year variation of the DT between SSN and I (GCR) was observed. DTs for periods A < 0 are greater than DTs for periods A > 0. 4. We also found that DTs in the A < 0 epoch are larger in the ascending I (GCR) period than in the descending I (GCR) period. In the A > 0 epoch, we did not find DTs in the ascending I (GCR) period, but small DTs in the descending I (GCR) period. 5. The main reasons for this phenomenon are essential temporal rearrangements of the structure of the IMF turbulence. The drift of GCRs caused by the gradient and curvature of the IMF is an additional factor, which strengthens this phenomenon. 6. To establish the properties of DTs from various cycles of SA, it is very important to model and experimentally analyze the propagation of cosmic ray particles in the heliosphere.