Energy Spectrum of the Recurrent Variation of Galactic Cosmic Rays During the Solar Minimum of Cycles 23/24

The Sun during the recent epoch of solar activity operated in a different way than during the last 60 years, being less active. We study temporal changes of the energy spectrum of the first three harmonics of the 27-day variation of the galactic cosmic rays (GCR) intensity during the unusual, recent solar minimum, between Solar Cycles 23 and 24 (SC 23/24) and compare with four previous minima. We show that the energy spectrum of the amplitudes of the recurrent variation of the GCR intensity is hard in the maximum epochs and is soft in the minimum epochs during Solar Cycles 20 – 24, but with peculiarities during the Solar Minimum 23/24. In particular, while the energy/rigidity spectrum of the amplitudes of the first harmonic of the recurrent variation of the GCR intensity behaves practically the same as for previous epochs, the energy/rigidity spectrum of the amplitudes of the second and the third harmonics demonstrates a pronounced softening. We attribute this phenomenon to the decrease of the extension of the heliosphere caused by the decrease of the solar-wind dynamic pressure during the unusual Solar Minimum 23/24.


Introduction
Galactic cosmic rays are highly energetic nuclei of extra-solar origin. When they are propagating through the heliosphere, they are partially modulated by the solar wind and heliospheric magnetic field (HMF). Moreover, the heliosphere changes over time, and those changes influence the GCR stream on different time scales (e.g. Kudela, 2013;Yeeram et al., 2014). The recurrent variation of the galactic cosmic-ray intensity was described for the  (Forbush, 1938) and is still of interest. Among recent publications devoted to this issue one can mention, e.g., Richardson (2004); Dunzlaff et al. (2008); Gil, Modzelewska, and Alania (2012); Guo and Florinski (2014). However, changes of the amplitudes of the recurrent variation of the GCR intensity versus the rigidity [R = T 2 + 2T T r , where T is kinetic energy and T r = 0.938 GeV is a proton's rest-mass energy] were studied very rarely until the work of Gil and Alania (2008). Alania (2010, 2013b) demonstrated that during the period 1965 -2002 the rigidity spectrum of the recurrent variation of the GCR intensity is hard in the maximum epochs, and it is soft in the minimum epochs of solar activity (SA). We relate this phenomenon to the changes of the effective size [L] of the modulation area of the recurrent variation of the GCR intensity in different epochs of SA. The effective size [L] of the modulation region of the GCR recurrence is smaller in the minimum epochs than in the maximum epochs. We base this explanation on the relation of the GCR particle's rigidity [R] with the product of Bρ (R ∝ Bρ), where B is the strength of the regular part of the HMF in the vicinity of heliolongitudinal asymmetry of the solar wind (being the source of the recurrent variation of the GCR intensity), and ρ is the cosmic-ray particle's Larmor radius. When other conditions are equal, the larger are ρ and B, the higher-energy are the cosmic-ray particles that are involved in the modulation, and the energy/rigidity spectrum of the recurrent variation of the GCR intensity, consequently, becomes harder.
The maximum epoch of the current SC 24 was about 50 % weaker than the maxima of the four previous SCs (e.g. Bazilevskaya et al., 2014;Balogh et al., 2014), similarly as the unusual Solar Minimum 23/24 (e.g. Dikpati, 2013) was weaker comparing with the previous minima. The sunspot number (SSN) is treated as a key indicator of solar activity. During SC 23 the lowest value of smoothed SSN was 1.7 in December 2008, while the lowest SSN in SC 22 was 12.6 (e.g. Dikpati, 2013). McComas et al. (2008) indicated a huge drop (≈ 36 %) in the average strength of the HMF. The lowest monthly value of HMF strength during the last minimum was 3.5 nT in April 2009. Sheeley (2010) concluded that the polar fields were weaker than they have been in the last 100 years. According to Hoeksema (2010) that behavior could signal a weaker maximum of Solar Cycle 24, which indeed took place (e.g. Balogh et al., 2014). These low values of the polar field throughout the last solar cycle caused a much slower decreasing of the tilt angle of the heliospheric current sheet (up to 2008), in contrast to Cycles 21 and 22. The average solar-wind speed (SWs) during this 23/24 solar minimum was 388 km s −1 , being 8 % slower than in the previous minimum, and 15 % slower than earlier ones (e.g. Jian, Russell, and Luhmann, 2011). This might be caused by the fact that the polar coronal holes were smaller than in previous cycles. Abramenko et al. (2010) reported that the area occupied by the near-equatorial coronal holes (±40 • ) was larger during the Minimum 23/24 in comparison with the previous solar minimum. Kirk et al. (2009) stated that the north and south polar-hole areas were reduced by ≈ 15 % in 2007 compared with 1996. Mason et al. (2009) found that during the recent Minimum 23/24 in comparison with the previous one (22/23) there were many more co-rotating interaction regions observed, connected with the existence of more high-speed streams of solar wind. Moreover, a low solar-wind dynamic pressure, ≈ 40 % in comparison with the previous minimum, was observed (Mewaldt et al., 2010) and this was caused by the drop in the average proton density, which was ≈ 5.5 cm −3 , being less dense in comparison with the previous three minima by ≈ 70 % (Jian, Russell, and Luhmann, 2011). All of the phenomena and processes mentioned caused diminishing of the solar-modulation influence on GCR intensity and its variations. McComas et al. (2008) recognized a shrinkage of the size of the heliosphere (caused by the drop of the solar-wind dynamic pressure), which should result in an increase in the GCR intensity at 1 AU (Mewaldt et al., 2010). Indeed, the GCR intensity values measured by neutron monitors were the highest ever recorded (e.g. Moraal and Stoker, 2010;Oh et al., 2013). In 2007, neutron monitors (NMs) recorded values last measured during the last A < 0 solar minimum, 1985 -1987, and the values continued to grow, exceeding the previous levels by a few percent, throughout the end of 2009. Stratospheric balloon measurements showed an even higher increase: ≈ 20 % (Stozhkov et al., 2009). The Payload for Antimatter Matter Exploration and Light-nuclei Astrophysics (PAMELA) experiment observed that from mid-2006 to the end of 2009 the cosmic-ray spectra became softer (up to the intensity of 2.7 particles(m 2 s sr MeV) −1 (e.g. Vos and Potgieter, 2015) since gradually more protons with lower energy were able to reach the Earth (Potgieter et al., 2014). Mewaldt et al. (2010) found ≈ 44 % growth of the parallel mean free path of GCR particles [λ ] from the minimum in 1997 to 2009 and its rapid reduction in 2010. Cliver, Richardson, and Ling (2013) stated that particle drifts were not the main processes in the recent minimum, but the GCR modulation was determined by diffusion and convection, contradicting the results obtained by Alania, Modzelewska, and Wawrzynczak (2014), i.e. that a drift related to the delay time between tilt angles of heliospheric current sheet and cosmic-ray intensity was an important reason for an increase of cosmic-ray intensity in 2009. Potgieter et al. (2014) indicated that diffusion and particle-drift contributions were both ≈ 50 % for the GCR protons observed at the Earth. The average mass and kinetic energy of coronal mass ejections were around five times lower in comparison with the previous minimum (Vourlidas et al., 2010), diminishing one of the barriers for entrance of GCR particles into the inner heliosphere (Mewaldt, 2013). Zhao et al. (2014) also stated that a low magnetic turbulence, increasing the parallel diffusion and reducing the perpendicular diffusion in the polar direction, might be another phenomenon giving rise to the observed higher GCR intensity.
The amplitudes of a new type of quasi-periodicity of the GCR intensity, found using the amplitudes of the 27-day variation, called the three to four Carrington-rotation period (3 -4 CRP) recurrence during the last minimum epoch had similar values as in the previous minimum, much lower than during the solar minimum one Hale cycle earlier (Gil and Alania, 2013a). It was ascribed to the anomaly increasing the diffusion of cosmic rays (Moraal and Stoker, 2010) during the Minimum 23/24. Dunzlaff et al. (2008) reported, based on Ulysses observational data, that in SC 22 the decrease in the component of the recurrent GCR arising from the fast solar wind showed a visible maximum from 25 • to 40 • and beyond, but in SC 23 there was no such substantial change. The authors suggested that the absence of large, stable, coronal-hole structures in SC 23 might be a reason for this dissimilarity (e.g. Kirk et al., 2009). Guo and Florinski (2014) stated that the GCR recurrent-modulation effect of co-rotating interaction regions played a significant role. Gil, Modzelewska, and Alania (2012) reported that the average amplitude of the recurrent variation of the GCR intensity calculated for the mid-latitude NMs was as high as during the previous minimum epoch 1996 (A > 0) and much higher than during the minimum one Hale cycle earlier; but at the same time, the average amplitude of the recurrent periodicity of the GCR anisotropy was at the same level as in the previous negative-polarity cycle. Modzelewska and Alania (2013) showed that the recurrent variation of the GCR intensity was stable during 2007 -2008 and evolved to longer periods (up to 33 -36 days) during 2009. Aslam and Badruddin (2015) showed, using wavelet analysis, that the recurrent variation during the last minimum remained steady compared with the four previous minima.
Because of all of those facts, it is important to check that the strength of the solar activity influenced the spectrum of the recurrence of the GCR intensity. Hence, the main aim of this article is to study the behavior of the rigidity spectrum of the recurrent variation of the galactic cosmic rays during the recent epochs of SA, particularly in the period covering the last solar-minimum and solar-maximum epochs.

Data and Mathematical Background
We present the computed rigidity spectrum of the amplitudes of the first three harmonics of the recurrent variation of the GCR intensity during the maximum of SC 23 up to the ongoing SC 24. For this purpose we used daily data from eight neutron monitors from a worldwide network, with different cut-off rigidities, in the period 01 January 2000 to 31 December 2014. Table 1 presents details, such as location and effective cut-off rigidity, of the neutron monitor (NM) stations used in this article.
Using harmonic analysis (e.g. Kincaid and Cheney, 2002) we computed the amplitudes of the first three harmonics of the 27-day variation of the GCR intensity by means of daily data from all NMs for each consecutive solar rotation during the period considered. For this purpose we used the formula where a n = 2 , 2, 3} is the amplitude of the j th harmonic of the 27-day variation of the GCR intensity.
The variability of the primary spectrum of the GCR flux in the heliosphere can be defined by (e.g. Dorman, 2004) where D i (R) means the primary spectrum of GCR during the ith rotation and D 0 (R) the primary spectrum for the considered period. A typical procedure of reconstructing the primary cosmic-ray characteristics from NMs measurements is by the usage of the neutron-monitor yield function Y (R, h) [in units of m 2 sr] dependent on the primary particle rigidity/energy and observational level h (e.g. Dorman, 1974;Clem and Dorman, 2000;Mishev, Usoskin, and Kovaltsov, 2013): where N means count rates of specific NM, D(R, t) primary spectrum (in units of particles(m 2 sr s GeV nucleon −1 ) −1 ) of GCR at the orbit of the Earth depending on rigid- ity and time. Thus, the amplitudes of the recurrence of the GCR intensity for the kth neutron monitor with characteristic cut-off rigidity R ck can be calculated as follows (e.g. Dorman, 2004): where [δD(R)/D(R)] i means the rigidity spectrum for the ith rotation. With the aid of the broadly used approximation, we define the rigidity spectrum of the recurrent variation as (e.g. Dorman, 1974;Yasue et al., 1982): where R max is a lower limit of the rigidity values above which the recurrence does not exist, in our article R max = 200 GV (e.g. Alania and Wawrzynczak, 2012), a is a constant, for R 0 being equal to δD(R)/D(R), and γ is the spectral index for the amplitudes of the harmonics of the 27-day variation of the GCR intensity. Using Equations (1) and (2) one can write Owing to the yield function, the 27-day variation amplitudes for each NM should be the same (to the accuracy available); hence in order to find spectral index γ for a particular solar rotation one has to minimize the expression (using the least-square method): where a av is an average amplitude of the 27-day variation of the GCR intensity for all NMs.
To compare the behavior of the spectral index γ of the 27-day variation of the GCR intensity during the ongoing SC 24 and especially the recent, unusual Solar Minimum 23/24 with its behavior during four previous SCs 20 -23 we need to have continuous values of γ for a long period. For this purpose, neutron monitors with substantially different cut-off rigidity, with well-established, continuous measurements over a long period, are essential. This requirement is reasonably well fulfilled by the Kiel and Rome NMs.
Using those NMs with the aid of the above-described method we compute values of the spectral index γ of the first three harmonics of the recurrent variation of the GCR intensity for the period of 01 January 1964 to 31 December 2014.

Behavior of the Rigidity Spectrum of the 27-Day Variation of the GCR Intensity During Solar Cycles 23/24
To study precisely the behavior of the energy/rigidity spectrum of the recurrent variation of the galactic cosmic-ray intensity during the current Solar Cycle 24 and the unusual Solar Minimum 23/24, we use eight various neutron monitors (Table 1) with considerably different cut-off rigidity during the period of 01 January 2000 to 31 December 2014. Figure 1 presents the results of the calculations. To emphasize the properties of the long-term changes of the rigidity spectra, we smoothed the values of the spectral index [γ ] with a running average over 39 solar rotations. The upper panel of Figure 1 shows that during the recent, unusual minimum epoch 23/24 the value of the spectral index [γ ] of the amplitudes of the first harmonic of the recurrent variation of the GCR intensity was the highest, being equal to 1.67 ± 0.08. This could be a signal that more particles with lower energy, being modulated by the recurrent solar behavior, were able to reach the Earth. Moreover, the lowest value of γ during the recent maximum was 0.96 ± 0.05, being 35 % higher in comparison with the value of γ during the maximum epoch of SC 23. The middle panel of Figure 1 shows that the rigidity spectrum of the amplitudes of the second harmonic of the GCR recurrence was extremely soft, with the highest value, γ ≈ 1.79 ± 0.09, being even larger than the γ of the first harmonic. The lower panel of Figure 1 illustrates that the lowest value of γ for the third harmonics during the recent maximum was 1.02 ± 0.05, being almost 60 % higher in comparison with the γ during the maximum of SC 23. This lower panel of Figure 1 shows that the rigidity spectrum of the amplitudes of the third harmonic of this variability of the GCR intensity during the recent minimum was soft, with the highest value being γ ≈ 1.43 ± 0.07. This could be caused by the fact that during the recent minimum epoch, at the end of SC 23, the second and the third harmonics of this variability were indeed enhanced, in addition to the enhancement in the solar-wind and solar-activity parameters (e.g. Modzelewska and Alania, 2013). To see changes of the spectral index [γ ] of the amplitudes of the GCR recurrent variation on the background of solar activity and cosmic-ray intensity, the upper panel of Figure 2 displays the 27-day averages of the Oulu NM count rates during the same period of 01 January 2000 to 31 December 2014. One can see a strong compatibility between the NM  measurements and the spectral index [γ ]. The correlation coefficient between the spectral index [γ ] and cosmic rays measured by the Oulu station is high: 0.89 ± 0.02.
One excellent proxy of solar activity is the solar radio flux (SRF or F 10.7 index). Its wavelength equals 10.7 cm and this has been measured continuously since 1947. SRF is derived from the solar atmosphere, high in the chromosphere and low in the solar corona. The middle panel of Figure 2 displays the 27-day averages of SRF during the period of 01 January 2000 to 31 December 2014. There also exists a very clear negative correlation between solar activity illustrated by SRF and the spectral index [γ ]. The anti-correlation between the spectral index γ and SRF is high: −0.89 ± 0.03.
Moreover, we have calculated the rigidity spectrum of the 11-year variation of the GCR intensity during the period of 01 January 2000 to 31 December 2014, using the same eight neutron monitors (Table 1). The lower panel of Figure 2 presents the spectral index [γ ] of the 11-year variation of the GCR intensity, illustrating the fact that there is a strong negative relationship between the changes of the spectral index [γ ] of the 11-year variation and the spectral index [γ ] of the recurrent variation of the GCR intensity, with the correlation coefficient being equal to −0.82 ± 0.03. This result is in good agreement with our findings for previous solar-activity cycles (Gil and Alania, 2011).
In order to make a quantitative, not only qualitative, comparison of the rigidity spectrum [δD(R)/D(R)] of the amplitudes of the first three harmonics of the recurrent variation of the GCR intensity during the ongoing SC 24 with previous SCs, we have calculated the value of the spectral index [γ ] of those three harmonics for the period of 01 January 1964 to 31 December 2014 by the method described above, using the Kiel and Rome NMs. Figure 3 shows that the rigidity spectrum of recurrent variation harmonics is, indeed, hard in the maximum and soft in the minimum epochs of solar activity, as was found by Alania (2010, 2013b). Figure 3 confirms that during the SC 24 the rule, concerning the hardening of the rigidity spectrum of the recurrent variation of the GCR intensity during the solarmaximum epochs, is maintained. Furthermore, Figure 3 endorses the idea that during the ongoing SC 24 values of the spectral index [γ ] for the first harmonic of this variability of the GCR intensity (upper panel of Figure 3) were at a comparable level as during previous SCs 20 -23, but values of γ for the amplitudes of the second (middle panel of Figure 3) and the third (lower panel of Figure 3) harmonics were much higher, from a few up to several dozen percent, than during the previous four solar-activity cycles. We have explained (Gil and Alania, 2011) the hardening of the rigidity spectrum of the recurrent variation of the GCR intensity in the maximum epochs of SA by the larger average zone of the co-rotating interaction regions (causing the recurrent variation of the GCR intensity) in comparison with the minimum epochs of SA. Hence, a shrinkage of the size of the heliosphere caused by the drop of the solar-wind dynamic pressure during the last Solar Minimum 23/24 (McComas et al., 2008), together with the enhancement of the second and the third harmonics of the recurrent variability effect, could be a significant reason for the exceptional softening of the rigidity spectrum of the second and third harmonics of this recurrence of the GCR intensity during the last Solar Minimum 23/24. Also, the peculiar temporal changes of the coronal-hole structure during SC 23/24 could affect the rigidity spectrum of the recurrent variation.

Conclusions
i) The energy/rigidity spectrum of the amplitudes of the first three harmonics of the recurrent variation of the galactic cosmic-ray intensity was hard during the maximum epoch of the lasting SC 24 and soft during the recent, unusually prolonged, minimum epoch between SCs 23 and 24. We connected this phenomenon with the changes of the effective size of the modulation region of the recurrent variation of the GCR intensity in different epochs of solar activity. ii) The energy/rigidity spectrum of the amplitudes of the second and the third harmonics of the 27-day recurrence of the GCR intensity was extremely soft during the last Solar Minimum 23/24 in comparison with the previous four solar minima. We ascribe this phenomenon to the reduction of the heliosphere's extension affected by the decrease of the solar-wind dynamic pressure combined with the enhancement of the second and the third harmonics of the 27-day variation during the Solar Minimum 23/24. In addition, the unusual temporal changes of the coronal-hole structure during SC 23/24 could affect the rigidity spectrum of the recurrent variation.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.