Next Article in Journal
Improvements and Extension of the Linear Carbon Sink Model
Previous Article in Journal
Climate Classification for Major Cities in China Using Cluster Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Precursory Signs of Large Forbush Decreases in Relation to Cosmic Rays Equatorial Anisotropy Variation

by
Maria-Christina Papailiou
1,
Maria Abunina
2,
Helen Mavromichalaki
1,*,
Nataly Shlyk
2,
Semyon Belov
2,
Artem Abunin
2,
Maria Gerontidou
1,
Anatoly Belov
2,
Victor Yanke
2 and
Amalia Triantou
1
1
Nuclear and Particle Physics Section, Physics Department, National and Kapodistrian University of Athens, 15784 Athens, Greece
2
Pushkov Institute of Terrestrial Magnetism, Ionosphere and Radio Wave Propagation, Russian Academy of Sciences (IZMIRAN), Moscow 108840, Russia
*
Author to whom correspondence should be addressed.
Submission received: 27 May 2024 / Revised: 18 June 2024 / Accepted: 20 June 2024 / Published: 21 June 2024
(This article belongs to the Special Issue Cosmic Rays, Ozone Depletion and Climate Change)

Abstract

:
Forbush decreases are usually characterized by increased values of cosmic ray anisotropy. The precursory signs, i.e., pre-increases and especially pre-decreases of the cosmic ray intensity, are highly anisotropic phenomena that ordinarily forewarn of such events. Two Cosmic Ray Groups from the National and Kapodistrian University of Athens (NKUA) and the Pushkov Institute of Terrestrial Magnetism, Ionosphere, and Radiowave Propagation of the Russian Academy of Sciences (IZMIRAN) have been investigating the existence of precursory signs preceding Forbush decreases in relation to different solar phenomena, interplanetary parameters, and geomagnetic conditions. In this study, large Forbush decreases (magnitude > 5%) accompanied by geomagnetic storms (i.e., geomagnetic index Dst < −100 nT and 5 ≤ Kp-index ≤ 9) and characterized by an equatorial anisotropy 1 h before the onset of the event (Axyb, %) less than 0.8% were examined regarding precursors. In total, 50 events with the aforementioned features were selected and analyzed from the IZMIRAN’s Forbush Effects and Interplanetary Disturbances database concerning the time period from 1969 until 2023. The Ring of Stations method, which depicts the cosmic ray variations for various asymptotic longitudes in relation to time, was applied on each event. The results revealed that clear signs of pre-decreases were not present for the majority of the events. Since particularly strong events were considered, most of them still showed some precursory signs, albeit mainly weak. Despite this, the value of Axyb = 0.8% proves to be a good threshold for the manual selection of FDs with well-expressed precursors.

1. Introduction

The Earth is bombarded by cosmic rays, i.e., relativistic charged particles that originate from extraterrestrial sources [1,2,3,4]. Specifically, galactic cosmic rays are composed of protons, α-particles (7–10%), and heavier nuclei (1%) [5,6,7,8], come from unknown sources outside the solar system [9], and the solar activity determines their flux in the heliosphere [10]. Furthermore, solar cosmic rays are related to solar flares [11,12] and coronal mass ejections [13] and mainly consist of protons (>89%), α-particles (10%), and heavier nuclei (<1%). The typical energy values of these particles are a few hundred MeV/nucleon, but sometimes reach GeV/nucleon [4,14]. A strong dependence on solar activity is also observed, since at the maximum of solar activity and during the descending phase of the solar cycle, an increased flux of low-energy particles is recorded [10].
One of the most important phenomena of galactic cosmic rays is the Forbush effect or Forbush decrease (FD). Decreases in the recording rate of cosmic rays, which last about a week, were first observed by Forbush in 1937 [15,16,17] using ionization chambers [18]. In particular, Ref. [16] analyzed the measurements from ionization chambers distributed around the Earth for 17 consecutive months to conclude that these intensity variations observed at stations around the world were related to each other. Furthermore, these intensity variations with time were irrespective of atmospheric effects. Initially, it was supposed that these variations were the direct or indirect result of geomagnetic variations, such as disturbances of the geomagnetic field during geomagnetic storms. That is why it was assumed that their origin was terrestrial [19].
Later on, in the early 1950s, it was concluded that these temporal variations were not owing to geomagnetic field disturbances alone [18]. Specifically, Refs. [20,21,22,23,24] recorded cosmic ray intensity (CRI) decreases at the geomagnetic pole. Ref. [23] continued his accurate research using neutron monitors to study the temporal intensity variations in a part of the lower-energy primary particles spectrum. Finally, he showed that the source of these variations was not related to geomagnetic field variations [19] but to the interplanetary medium [25,26] and, specifically, to solar activity phenomena [27,28]. Moreover, he singled out the sharp, non-periodic CRI decreases in magnitude greater than 6%, which would later be defined as FDs. These initial findings are considered to have laid the foundation for the later development of the field of cosmic ray modulation.
Over the years, many definitions have been proposed for describing FDs. The definition as ‘a cosmic ray variation during a geomagnetic storm’ mentioned in [29] is considered to be out of date, since FDs are neither always accompanied by geomagnetic storms nor are always recorded near the Earth. Moreover, the Glossary of Solar-Terrestrial Terms NOAA (https://www.ngdc.noaa.gov/stp/glossary/glossary.html, accessed on 12 May 2024) defines FDs as ‘the sudden variation of about 10% in the galactic cosmic rays intensity, as recorded by the neutron monitors’. However, this also is not a precise definition, since an FD is a phenomenon that evolves gradually and usually has a magnitude of less than 10%, in some rare cases. In addition, this definition assumes that these phenomena can only be observed with neutron monitors, which is not accurate since other detectors are capable of recording FDs, i.e., ionization chambers, ground-based or underground muon detectors, as well as other detectors used in the study of cosmic rays in space. Essentially, the FD could be described as a cosmic ray storm, during which the galactic cosmic ray flux becomes anisotropic. Therefore, the FD is a cosmic ray storm as well as a heliospheric event [30]. According to their origin, FDs could be defined as ‘the effect of coronal mass ejections (CMEs and ICMEs) and/or fast-moving solar wind streams from coronal holes on cosmic rays’ [31,32,33,34]. After having mentioned all the above, the most contemporary definitions describe FDs as a short and abrupt decrease in the galactic CRI that is accompanied by a relatively slow recovery (lasting almost one week) [35,36,37,38].
One of the main subjects when studying FDs is the existence of precursory signs preceding the main phase of the CRI decrease. Refs. [39,40] were the first to introduce and study the decreases and/or increases in the CRI before the evolution of an FD. Nowadays, the investigation of precursory signs has developed into an essential chapter of space weather research [41,42,43,44], etc. and an important contributor to the forewarning of space weather phenomena [45,46,47,48].
Responsible for the occurrence of precursors are the kinetic interactions between cosmic ray particles and the approaching interplanetary shock [42,49]. More specifically, the pre-decrease is recorded when cosmic ray particles exit the FD zone along the magnetic field lines. This happens because the Earth and the cosmic ray-depleted region behind the shock front are magnetically connected, i.e., the “loss cone” effect [42,50]. On the other hand, cosmic ray particles reflect from the approaching interplanetary disturbance, and furthermore, their acceleration at the front of the approaching shock results in pre-increase recording [51,52]. The two precursory signs, i.e., pre-decreases and/or pre-increases in the CRI, differ regarding the mechanism that creates them, but coexist in the area in front of the shock wave [50]. Both types of precursors can last from some hours up to one day before the main phase of the decrease.
The study of precursory signs has concerned many scientific groups around the world. The Athens Cosmic Ray Group of the National and Kapodistrian University of Athens (NKUA) and the Cosmic Ray Group of the Pushkov Institute of Terrestrial Magnetism, Ionosphere, and Radio Wave Propagation (IZMIRAN) of the Russian Academy of Sciences are conducting an ongoing investigation through which different solar parameters, interplanetary conditions, and geomagnetic disturbances are examined in relation to their effectiveness concerning the existence of precursors. The Ring of Stations method is used for studying the precursory signs of FDs. A more detailed description of the method used is included in the next section (see Section 2). These findings are collected and expressed in the form of different criteria.
Moreover, two more attempts in monitoring precursory signs preceding FDs are described herein. First, the global network of neutron monitors called ‘Spaceship Earth’ of the Bartol Research Institute is mentioned. Ref. [53] proposed to set up this neutron monitor network, since neutron monitor stations located at high latitudes provide an excellent opportunity to study cosmic ray events. As a result, a multinational cooperation (USA, Russia, Australia, Canada) took place (https://neutronm.bartol.udel.edu/catch/sse1.html, accessed on 12 May 2024). This network includes 11 stations (Inuvik—Canada, Fort Smith—Canada, Peawanuck—Canada, Nain—Canada, Mawson—Antarctica, Apatity—Russia, Norilsk—Russia, Tixie Bay—Russia, Cape Schmidt—Russia, Thule—Greenland, McMurdo—Antarctica) which are located at high latitudes on four continents and are distinguished by the most favorable directional sensitivity.
The Bartol Research Institute runs an efficient precursory signs monitoring system (http://neutronm.bartol.udel.edu/, accessed on 12 May 2024). It is based on the cosmic rays distribution regarding pitch angle, and even though it uses a limited number of neutron monitor stations, they are located uniformly in regard to longitude. Amongst other interplanetary parameters, the hourly average of the CRI measured by one Spaceship Earth station in relation to the cosmic ray density (for neutron monitor data) or a single directional channel in relation to the cosmic rays density (for muon detector data) is depicted with a circle. The red color describes intensity deficit, while the blue color describes intensity excess. The deficit or excess’s magnitude corresponds to the size of the circle. The vertical and horizontal axes display the pitch angle of the station or the directional channel and the time (last 3 days), respectively (https://neutronm.bartol.udel.edu/spaceweather/welcome.html, accessed on 12 May 2024).
However, this system can potentially have two drawbacks, since a possible operational failure of one of the stations can be responsible for an almost 50° gap in the scan of the celestial sphere, and the uninterrupted interplanetary magnetic field data flow is essential, since these data are used for pitch angle calculations [54].
The search for precursory signs with muon detectors has also been developed and rather successfully [43,52,55,56]. Two investigations related to precursory signs of cosmic rays were conducted using ground-based detectors. As mentioned in [57], 14 major geomagnetic storms were studied using the neutron monitor network. These storms were listed by [58] and were characterized by Kp 8– or greater. Another 22 large geomagnetic storms were studied using data from ground-based muon detectors [55], i.e., a large storm is defined as one characterized by Kp 7– or greater [59]. It should be noted here that neutron monitors as well as muon detectors record the secondary particles’ intensity produced when primary cosmic rays, mainly protons, interact with atoms in the atmosphere of the Earth. The typical energy of the primary cosmic rays that will produce the secondary particles (which will, in turn, be modulated by FDs) is ∼10 GeV and ∼30 GeV for neutron monitors and muon detectors, respectively [42].
It was concluded that the majority of large or major geomagnetic storms are accompanied by precursors. In 68% of events (15 out of 22 large storms) registered by muon detectors and in 79% of events (11 out of 14 major storms) registered by neutron monitors, the precursory signs were clear. As reported by [55], if the muon detectors investigation included only the largest storms (with a Kp value of 8 or greater), then the corresponding percent increases to 89%. It was noted that the precursor signs preceded the sudden storm commencement (SSC) usually by a few to 12 h. This time interval is sufficient to be used in space weather forewarning. Moreover, there is a difference between the duration of the precursory signs provided by neutron monitors and muon detectors. Thus, precursory signs are recorded 8 and 4 h prior to the main event for muon detectors and neutron monitors, respectively [42].
It is known that cosmic rays present stability on a galactic scale and show signs of isotropy in their spatial distribution [60,61]. The mean statistical value of the equatorial component of the first harmonic of cosmic ray anisotropy, Axy, is less than 0.6% [62], and more precisely, as mentioned in [63], this value is almost 0.51–0.53% during quiet conditions, i.e., ambient solar wind. However, disturbances created in the Sun and moving through interplanetary space may affect cosmic rays. An increase in the solar wind speed, the interplanetary magnetic field, or the cosmic ray density gradient may result in CRI variations and enhanced anisotropies, i.e., an increase in the value of Axy and significant variations in the vector’s direction. Moreover, the anisotropy 1 h before the FD onset, Axyb, increases to almost 0.7%, but may jump to ~1% for large (>4%) FDs [64].
As it is understood, the precursory signs preceding FDs are a highly anisotropic phenomenon. FDs with an Axyb greater than 1.2%, regardless of their magnitude, were examined for precursors in [62] and provided clear and irrefutable results. Yet, because this value exceeded the expected Axyb before an FD by far, it was imperative to study precursory signs regarding smaller values of anisotropy. Indeed, in the studies that followed, the threshold value for Axyb was set at 0.8%. Precursors were evident in large FDs (magnitude > 5%), accompanied by geomagnetic storms (geomagnetic index Dst < −100 nT and 5 ≤ Kp-index ≤ 9), and characterized by Axyb greater than 0.8% [65].
In this work, large FDs (i.e., magnitude > 5%) with an equatorial component of the first harmonic of cosmic ray anisotropy 1 h before the development of the event (Axyb, %) less than 0.8% that are accompanied by geomagnetic storms (i.e., geomagnetic index Dst < −100 nT and 5 ≤ Kp-index ≤ 9) are being examined for precursory signs. These events cover the period from 1969 until 2023.

2. Data and Method

The events under study were selected from the IZMIRAN’s Forbush Effects and Interplanetary Disturbances (FEID) database [66], which can be found online (https://tools.izmiran.ru/feid, accessed on 12 May 2024), and met the conditions that were mentioned above. Cosmic ray variation parameters in this database are computed using the Global Survey Method (GSM) [67] for the 10 GV rigidity particles with the neutron monitor network data. For example, the zero spherical harmonic or isotropic part of the cosmic ray flux and the first spherical harmonic (i.e., density and anisotropy, respectively), both calculated in percentage, are provided by this method. The aforementioned anisotropy is the hourly equatorial component Axy of the cosmic ray anisotropy, i.e., the diurnal anisotropy Axy = A x 2 + A y 2 , where Ax and Ay are components of the cosmic ray anisotropy vector.
For the analysis and description of each event, a significant number of data, inserted in the FEID, were required. For example, data on solar flares, i.e., type, class, location, time (ftp://ftp.swpc.noaa.gov/pub/indices/events/, accessed on 12 May 2024), interplanetary conditions (solar wind speed and interplanetary magnetic field) and geomagnetic activity (geomagnetic indices Dst, Ap, and Kp) provided by the OMNI database (http://omniweb.gsfc.nasa.gov, accessed on 12 May 2024), SSCs (http://isgi.unistra.fr/data_download.php, accessed on 12 May 2024), etc. were assessed.
The FEID is an important tool in the study of cosmic ray phenomena since it includes all CRI variations taking place from 1957 (the beginning of neutron monitors operation) onwards. It provides the opportunity for retrospective studies of almost seven solar cycles (beginning from solar cycle 19 until the current solar cycle 25).
For the monitoring of the precursory signs of FDs, the Ring of Stations method (RSM) was applied (https://tools.izmiran.ru/ros, accessed on 12 May 2024). This method uses hourly CRI data derived from the neutron monitor stations located in different parts of the world that meet certain criteria. In particular, stations must have a geomagnetic cut-off rigidity Ri less than 4 GV, a coupling coefficient for the North–South component of the cosmic ray anisotropy C 10 i less than 0.55, and, finally, be located in altitude hi less than 1200 m. In total, the method derives data from 36 neutron monitor stations, which are presented in detail in [50]. Nevertheless, due to operational issues that the stations can sometimes have, this number may vary occasionally.
Furthermore, the hourly values of the CRI variations at each station are calculated in regard to a quiet period. They are, then, plotted for various asymptotic longitudes in relation to time. It should be mentioned that each neutron monitor station records particles from certain asymptotic directions. The sufficient number of neutron monitor stations that are used by this method ensures as much as possible the most complete scanning of the celestial sphere.

3. Results

In this study, large FDs (with magnitude > 5%), accompanied by geomagnetic storms (geomagnetic index Dst < −100 nT and 5 ≤ Kp-index ≤ 9) and characterized by Axyb less than 0.8%, were investigated regarding precursory signs. In total, 50 events covering the period from 1969 to 2023 were examined. These events are presented in Table 1, along with their magnitude (%) and the equatorial component of the first harmonic of cosmic ray anisotropy 1 h before the shock arrival (%). As it is noticed, some events break the rule that dictates increased anisotropy during a large disturbance. For instance, the FD on 13 March 1989 had a magnitude of 20.9% but the Axyb was only 0.11%. Another example of a large event (magnitude of 20.9%) with a small Axyb (0.38%) was the FD on 14 February 1978. This is most likely associated with the characteristics of the solar sources, the geometry of particle acceleration, and the particles’ asymptotic directions of arrival. It should be noted that around half of the events in question happened in 24 h or less after another FD. Though this situation is to be expected in periods of high solar activity, it makes precursors and anisotropy study in general much harder, because a phenomenon which may be regarded as a precursor of an event that may as well turn out to be related to the development of the previous one.
Figure 1a shows the CRI variations for various asymptotic longitudes with time for the FDs on 14 February 1978. In these plots, which were derived from the RSM, magenta-colored circles depict CRI decreases and blue colored circles depict CRI increases (both estimated with respect to an undisturbed base period). Each neutron monitor station’s recording for 24 h is demonstrated diagonally (i.e., line of circles). The circle’s size depends on the size of the variation; the vertical line denotes the FDs’ onset. The asymptotic longitude (in o) and the time are marked in the vertical and horizontal axes, respectively. Additionally, for the event on 14 February 1978, data on cosmic ray variations calculated by the GSM are presented (Figure 1b): green curve (A0, left scale)—cosmic ray density variations; pink bars (Axy, right scale)—changes in hourly values of the equatorial component of cosmic ray vector anisotropy. The Axyb parameter used to select events in this study is the bar before the event start hour (vertical line).
The event in Figure 1a happened on 14 February 1978, and as was already mentioned, while being of great magnitude (20.9%), it shows a small Axyb (0.38%), and indeed does not show any clear precursory signs, which is an unexpected behavior for an event of such size. The second event on 6 August 1982 (Af = 6%, Dstmin = −155 nT, Kpmax = 8–, and Axyb = 0.42%) is a good example of the Axyb criterion. The Axy is quite low before this event, and no clear precursors are visible in the picture (Figure 2). Though the weak increase between longitudes 180° and 320° may be related to the shock wave, it is most likely just related to background solar diurnal anisotropy variation.
Note that sometimes anisotropy may be visible on the RSM pictures, while the value of Axy computed with GSM is very low. This can be explained by the nature of both methods: while RSM uses count variations directly, the GSM uses more stations and fits a spherical model to their data, which inevitably leads to some smoothing effect.
However, due to the nature of the sample, which includes large FDs, many of them are expected to at least have a pre-increase. Since pre-increase is caused by particle acceleration on the interplanetary shockwave, it thus depends on the shockwave speed, and for the discussed sample of events the average solar wind speed is very high (740 km/s). Two FDs that still show precursory signs while having relatively low Axyb are presented in Figure 3.
The precursory signs for the FD on 9 January 1983 (Af = 9.4%, Dstmin = −213 nT, Kpmax = 8+, and Axyb = 0.32%) indicate a pre-increase on all longitudes, lasting about 8 h (Figure 3a). For the FD on 24 November 2001 (Af = 9.8%, Dstmin = −221 nT, Kpmax = 8.33, and Axyb = 0.67%), a pre-increase is also noticed in the longitudinal zone 90°–270°, lasting for almost 6 h (Figure 3b). As can be concluded, by the size of the circles, these precursory signs are relatively weak and last a few hours until the development of the main phase of the FD.
Another situation where the Axyb criterion does not work as expected is observed in the relatively recent event of 23 April 2023, shown in Figure 4. A very significant pre-increase is observable as blue circles on all longitudes. This effect may be explained by the cosmic ray particles being trapped in a compressing space between two shockwaves, the second one moving faster than the first. Such a large and almost isotropic increase led to very low values of Axyb.

4. Conclusions and Discussion

Cosmic rays are an important part of space weather. In particular, the study of precursory signs of Forbush decreases can be useful for forecasting large interplanetary disturbances, causing significant geomagnetic storms. These in turn cause changes in the upper atmosphere and other effects, which can be harmful to satellites and other technology [68].
To be able to forewarn of FD events, the NKUA and the IZMIRAN’s Cosmic Ray Groups, in this study, examined whether or not precursory signs were present before the development of large FDs (magnitude > 5%), accompanied by geomagnetic storms (Dst-index < −100 nT and 5 ≤ Kp-index ≤ 9), and characterized by an equatorial component of the first harmonic of cosmic ray anisotropy 1 h before the shock arrival less than 0.8%. In total, 50 events, which occurred from 1969 until 2023, were analyzed.
Events with an equatorial component of the first harmonic of cosmic ray anisotropy 1 h before the shock arrival greater than 0.8% have also been examined and analyzed in a different study [65]. However, this does not affect the novelty of this particular work. The two studies are complementary and not repetitive or even contradictory. The threshold (0.8%) of anisotropy is fully investigated (for greater or lesser values) and, as is shown, the chosen threshold is correctly chosen since the results differ according to this.
It is expected for large events to have great values of anisotropy 1 h before the shock arrival. However, herein, large events with small values of anisotropy 1 h before the shock arrival are studied. This on its own is a very unique situation, indicating specific environmental conditions. Moreover, the fact that the selected events presented signs of a pre-increase and not clear signs of pre-decrease is an important result which helps us understand the physical mechanisms related to these events and evaluate the methods used to analyze them.
The results of this work are summarized as follows:
(1) Sixteen out of the 50 events under study show some signs of a pre-increase, due to the fast-moving interplanetary shock waves. These are the events on 24 April 1979, 19 December 1980, 5 September 1982, 9 January 1983, 13 March 1989, 8 June 1989, 12 March 1990, 28 July 1990, 4 June 1991, 20 February 1992, 9 May 1992, 24 November 2001, 17 April 2002, 29 May 2003, 14 December 2006, and 14 July 2012;
(2) For the majority of the events under study, clear signs of pre-decrease were not observed. Only three events out of 50 showed signs of pre-decrease (7 November 2004, 11 September 2005, and 26 September 2011);
(3) Only two events (1 March 1982 and 22 January 2004) have signs of pre-decrease and pre-increase;
(4) The Ring of Stations method is more sensitive to precursor-like anomalies, which explains some inconsistencies with the Global Survey Method 10 GV anisotropy approximations;
(5) Significant FDs often appear consecutively, which makes the study of the precursors harder;
(6) Setting the threshold value of Axyb at 0.8% helps to select events with clearly expressed precursors, but ideally more criteria should be used along with it.

Author Contributions

Software, A.A., A.B. and S.B.; data curation, A.B.; data analysis, A.T.; writing—original draft preparation, M.-C.P.; writing—review and editing, M.A., N.S., S.B. and M.G.; supervision, H.M., A.B. and V.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This work received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are available at https://tools.izmiran.ru/feid (accessed on 12 May 2024).

Acknowledgments

This research is supported by ESA SSA SWE Space Radiation Expert Service Centre activities (ESA contract number 4000113187/15/D/MRP) and the ESA Space Safety Program’s network of space weather service development and pre-operational activities (ESA contract number 4000134036/21/D/MRP). The high-resolution Neutron Monitor Database (NMDB) is funded by the European Union. A.Ne.Mo.S is supported by the Special Research Account of Athens University (70/4/5803).

Conflicts of Interest

The authors declare that there are no conflicts of interest.

References

  1. Grieder, P.K.F. Cosmic Rays at Earth; Elsevier: Amsterdam, The Netherlands, 2001. [Google Scholar]
  2. Vainio, R.; Desorgher, L.; Heynderickx, D.; Storini, M.; Flückiger, E.; Horne, R.B.; Kovaltsov, G.A.; Kudela, K.; Laurenza, M.; McKenna-Lawlor, S.; et al. Dynamics of the Earth’s particle radiation environment. Space Sci. Rev. 2009, 147, 187–231. [Google Scholar] [CrossRef]
  3. Drury, L.O. Origin of cosmic rays. Astropart. Phys. 2012, 39–40, 52–60. [Google Scholar] [CrossRef]
  4. Mironova, I.A.; Aplin, K.L.; Arnold, F.; Bazilevskaya, G.A.; Harrison, R.G.; Krivolutsky, A.A.; Nicoll, K.A.; Rozanov, E.V.; Turunen, E.; Usoskin, I.G. Energetic Particle Influence on the Earth’s Atmosphere. Space Sci. Rev. 2015, 194, 1–96. [Google Scholar] [CrossRef]
  5. Baker, D.N. What is Space Weather. Adv. Space Res. 1998, 22, 7–16. [Google Scholar] [CrossRef]
  6. Masarik, J.; Beer, J. Simulation of particle fluxes and cosmogenic nuclide production in the Earth’s atmosphere. J. Geophys. Res. 1999, 104, 12099–12111. [Google Scholar] [CrossRef]
  7. Daglis, I.A. Effects of Space Weather on Technology Infrastructure; NATO ASI Series; Wiley: Hoboken, NJ, USA, 2004. [Google Scholar]
  8. Singh, A.K.; Singh, D.; Singh, R.P. Impact of galactic cosmic rays on Earth’s atmosphere and human health. Atmos. Environ. 2011, 45, 3806–3818. [Google Scholar] [CrossRef]
  9. Li, W.-H.; DAMPE Collaboration. A study of Forbush Decreases effects with DAMPE experiment. In Proceedings of the 38th International Cosmic Ray Conference (ICRC2023), Nagoya, Japan, 26 July–3 August 2023; p. 1311. [Google Scholar] [CrossRef]
  10. Daglis, I.A. Space Storms and Space Weather Hazards; NATO ASI Series; Wiley: Hoboken, NJ, USA, 2001. [Google Scholar]
  11. Miroshnichenko, L. Solar Cosmic Rays; Kluwer Academic Press: Norwell, MA, USA, 2001. [Google Scholar]
  12. Miroshnichenko, L.I. Solar cosmic rays in the system of solar–terrestrial relations. J. Atmos. Sol. Terr. Phys. 2008, 70, 450–466. [Google Scholar] [CrossRef]
  13. Anastasiadis, A.; Lario, D.; Papaioannou, A.; Kouloumvakos, A.; Vourlidas, A. Solar energetic particles in the inner heliosphere: Status and open questions. Phil. Trans. R. Soc. 2019, A377, 20180100. [Google Scholar] [CrossRef] [PubMed]
  14. Bazilevskaya, G.A. Once again about origin of the solar cosmic rays. J. Phys. Conf. Ser. 2017, 798, 012034. [Google Scholar] [CrossRef]
  15. Forbush, S.E. On the Effects in the Cosmic-Ray Intensity Observed During the Recent Magnetic Storm. Phys. Rev. 1937, 51, 1108–1109. [Google Scholar] [CrossRef]
  16. Forbush, S.E. On the World-Wide Changes in Cosmic-Ray Intensity. Phys. Rev. 1938, 54, 975. [Google Scholar] [CrossRef]
  17. Hess, V.F.; Demmelmair, A. World-wide Effect in Cosmic Ray Intensity, as Observed during a Recent Magnetic Storm. Nature 1937, 140, 316. [Google Scholar] [CrossRef]
  18. Kane, R.P. Severe geomagnetic storms and Forbush decreases: Interplanetary relationships re-examined. Ann. Geophys. 2010, 28, 479–489. [Google Scholar] [CrossRef]
  19. Lockwood, J.A. Forbush decreases in the cosmic radiation. Space Sci. Rev. 1971, 12, 658–715. [Google Scholar] [CrossRef]
  20. Singer, S.F. Storm decrease of the cosmic radiation at the pole. Phys. Rev. 1954, 95, 647. [Google Scholar]
  21. Singer, S.F. Observations of Cosmic Ray Decreases at the Pole. Nuovo C. 1958, 8, 326–333. [Google Scholar] [CrossRef]
  22. Simpson, J.A. Neutrons Produced in the Atmosphere by the Cosmic Radiations. Phys. Rev. 1951, 83, 1175–1188. [Google Scholar] [CrossRef]
  23. Simpson, J.A. Cosmic-Radiation Intensity-Time Variations and Their Origin. III The Origin of 27-Day Variations. Phys. Rev. 1954, 94, 426–440. [Google Scholar] [CrossRef]
  24. Simpson, J.A.; Fonger, W.; Treiman, S. Cosmic Radiation Intensity-Time Variations and Their Origin. I. Neutron Intensity Variation Method and Meteorological Factors. Phys. Rev. 1953, 90, 934–950. [Google Scholar] [CrossRef]
  25. Cane, H.V. Coronal mass ejections and Forbush decreases. Space Sci. Rev. 2000, 93, 55–77. [Google Scholar] [CrossRef]
  26. Qin, G.; Wu, S.-S. Magnetic Cloud and Sheath in the Ground-level Enhancement Event of 2000 July 14. II. Effects on the Forbush Decrease. Astrophys. J. 2021, 908, 236. [Google Scholar] [CrossRef]
  27. Cane, H.V.; Richardson, I.G.; Wibberenz, G. The Response of Energetic Particles to the Presence of Ejecta Material. Int. Cosm. Ray Conf. 1995, 4, 377. [Google Scholar]
  28. Baral, R.; Adhikari, B.; Calabia, A.; Shah, M.; Mishra, R.K.; Silwal, A.; Bohara, S.; Manandhar, R.; del Peral, L.; Frías, M.D.R. Spectral Features of Forbush Decreases during Geomagnetic Storms. J. Atmos. Sol. Terr. Phys. 2023, 242, 10598. [Google Scholar] [CrossRef]
  29. Dorman, L.I. Cosmic Ray Variations and Space Research; Nauka: Moscow, Russia, 1963; p. 1027. [Google Scholar]
  30. Belov, A.V. Forbush Effects and Their Connection with Solar, Interplanetary and Geomagnetic Phenomena; Gopalswamy, N., Webb, D.F., Eds.; Proceedings of the International Astronomical Union, IAU Symposium S257: Universal Heliophysical Processes; Cambridge University Press: Cambridge, UK, 2008; pp. 439–450. [Google Scholar] [CrossRef]
  31. Richardson, I.G. Energetic Particles and Corotating Interaction Regions in the Solar Wind. Space Sci. Rev. 2004, 111, 267. [Google Scholar] [CrossRef]
  32. Richardson, I.G. Solar wind stream interaction regions throughout the heliosphere. Living Rev. Sol. Phys. 2018, 15, 1. [Google Scholar] [CrossRef] [PubMed]
  33. Dumbovic, M.; Vrsnak, B.; Calogovic, J.; Karlica, M. Cosmic ray modulation by solar wind disturbances. Astron. Astrophys. 2011, 531, A91. [Google Scholar] [CrossRef]
  34. Dumbovic, M.; Vrsnak, B.; Calogovic, J.; Zupan, R. Cosmic ray modulation by different types of solar wind disturbances. Astron. Astrophys. 2012, 538, A28. [Google Scholar] [CrossRef]
  35. Wawrzynczak, A.; Alania, M. Modeling and data analysis of a Forbush decrease. Adv. Space Res. 2010, 45, 622–631. [Google Scholar] [CrossRef]
  36. Arunbabu, K.P.; Antia, H.M.; Dugad, S.R.; Gupta, S.K.; Hayashi, Y.; Kawakami, S.; Mohanty, P.K.; Oshima, A.; Subramanian, P. How are Forbush decreases related to interplanetary magnetic field enhancements? Astron. Astrophys. 2015, 580, A41. [Google Scholar] [CrossRef]
  37. Zhao, L.L.; Zhang, H. Transient Galactic Cosmic-Ray Modulation During Solar Cycle 24: A Comparative Study of Two Prominent Forbush Decrease Events. Astrophys. J. 2016, 827, 13. [Google Scholar] [CrossRef]
  38. Dumbović, M.; Kramarić, L.; Benko, I.; Heber, B.; Vršnak, B. A new method of measuring Forbush decreases. arXiv 2024, arXiv:2401.09000. [Google Scholar] [CrossRef]
  39. Blokh, Y.L.; Dorman, L.I.; Kaminer, N.S. Proceedings of the Moscow Cosmic Ray Conference, July 6–11, 1959; International Conference of Cosmic Radiation: Moscow, Russia, 1960. [Google Scholar]
  40. Fenton, A.G.; McCracken, K.G.; Rose, D.C.; Wilson, B.G. The onset times of Forbush-type cosmic ray intensity decreases. Can. J. Phys. 1959, 37, 970–982. [Google Scholar] [CrossRef]
  41. Nagashima, K.; Sakakibara, S.; Fujimoto, K.; Fujji, Z.; Ueno, H. Local-time- dependent precursory decrease of cosmic rays, in front of Forbush-decrease-associated IMF shock wave, observed on September 9, 1992. In Proceedings of the 23rd ICRC, Calgary, AB, Canada, 19–30 July 1993; Volume 3, pp. 711–714. [Google Scholar]
  42. Leerungnavarat, K.; Ruffolo, D.; Bieber, J.W. Loss cone precursors to Forbush decreases and advance warning of space weather effects. Astrophys. J. 2003, 593, 587–596. [Google Scholar] [CrossRef]
  43. Munakata, K.; Kuwabara, T.; Yasue, S.; Kato, C.; Akahane, S.; Koyama, M.; Ohashi, Y.; Okada, A.; Aoki, T.; Mitsui, K.; et al. A “Loss-Cone” Precursor of an Approaching Shock Observed by a Cosmic-Ray Muon Hodoscope on October 28, 2003. Geophys. Res. Lett. 2005, 32, L03S04. [Google Scholar] [CrossRef]
  44. Fushishita, A.; Kuwabara, T.; Kato, C.; Yasue, S.; Bieber, J.W.; Evenson, P.; Da Silva, M.R.; Dal Lago, A.; Schuch, N.J.; Tokumaru, M.; et al. Precursors of the Forbush Decrease on 2006 December 14 Observed with the Global Muon Detector Network (GMDN). Astrophys. J. 2010, 715, 1239–1247. [Google Scholar] [CrossRef]
  45. Dorman, L.I. Space weather and dangerous phenomena on the Earth: Principles of great geomagnetic storms forcasting by online cosmic ray data. Ann. Geophys. 2005, 23, 2997–3002. [Google Scholar] [CrossRef]
  46. Badruddin. Transient perturbations and their effects in the heliosphere, the geo-magnetosphere, and the Earth’s atmosphere: Space weather perspective. J. Astrophys. Astron. 2006, 27, 209–2017. [Google Scholar] [CrossRef]
  47. Pulkkinen, T. Space Weather: Terrestrial Perspective. Living Rev. Sol. Phys. 2007, 4, 1. Available online: http://www.livingreviews.org/lrsp-2007-1 (accessed on 12 May 2024). [CrossRef]
  48. Badruddin, B.; Aslam, O.P.M.; Derouich, M.; Asiri, H.; Kudela, K. Forbush decreases and geomagnetic storms during a highly disturbed solar and interplanetary period, 4–10 September 2017. Space Weather 2019, 17, 487–496. [Google Scholar] [CrossRef]
  49. Ruffolo, D.; Bieber, J.W.; Evenson, P.; Pyle, R. Precursors to Forbush decreases and space weather prediction. In Proceedings of the 26th International Cosmic-Ray Conference, Salt Lake City, UT, USA, 17–25 August 1999; Kieda, D., Salamon, M., Dingus, B., Eds.; International Union of Pure and Applied Physics: Salt Lake City, UT, USA, 1999; Volume 6, pp. 440–443. [Google Scholar]
  50. Abunina, M.A.; Belov, A.V.; Eroshenko, E.A.; Abunin, A.A.; Yanke, V.G.; Melkumyan, A.A.; Shlyk, N.S.; Pryamushkina, I. Ring of Stations Method in Cosmic Rays Variations Research. Sol. Phys. 2020, 295, 69. [Google Scholar] [CrossRef]
  51. Dorman, L.I.; Iucci, N.; Villoresi, G. The nature of cosmic ray Forbush decrease and precursory effects. In Proceedings of the 24th International Cosmic-Ray Conference, Rome, Italy, 28 August–8 September 1995; Iucci, N., Lamanna, E., Eds.; International Union of Pure and Applied Physics: Rome, Italy, 1995; Volume 4, pp. 892–895. [Google Scholar]
  52. Kudela, K.; Storini, M. Possible tools for space weather issues from cosmic ray continuous records. Adv. Space Res. 2006, 37, 1443–1449. [Google Scholar] [CrossRef]
  53. Bieber, J.W.; Evenson, P. Spaceship Earth—An optimized network of neutron monitors. In Proceedings of the 24th International Cosmic Ray Conference, Rome, Italy, 28 August–8 September 1995; Volume 4, pp. 1316–1319. [Google Scholar]
  54. Asipenka, A.; Belov, A.V.; Eroshenko, E.; Mavromichalaki, H.; Papailiou, M.; Papaioannou, A.; Oleneva, V.; Yanke, V.G. Asymptotic longitudinal distribution of cosmic ray variations in real time as the method of interplanetary space diagnostic. In Proceedings of the 31st International Cosmic Ray Conference, Lodz, Poland, 7–15 July 2009; Available online: http://cosray.phys.uoa.gr/index.php/publications-menu/publications (accessed on 12 May 2024).
  55. Munakata, K.; Bieber, J.W.; Yasue, S.; Kato, C.; Koyama, M.; Akahane, S.; Fujimoto, K.; Fujii, Z.; Humble, J.E.; Duldig, M.L. Precursors of geomagnetic storms observed by muon detector network. J. Geophys. Res. 2000, 105, 27457–27468. [Google Scholar] [CrossRef]
  56. Munakata, K.; Yasue, S.; Kato, C.; Kota, J.; Tokumaru, M.; Kojima, M.; Darwish, A.A.; Kuwabara, T.; Bieber, J.W. On the cross-field diffusion of galactic cosmic rays into an ICME, WSPC/SPI-B368. Adv. Geosci. 2006, 2, 115–124. [Google Scholar]
  57. Belov, A.V.; Bieber, J.W.; Eroshenko, E.A.; Evenson, P.; Pyle, R.; Yanke, V.G. Pitch-angle features in cosmic rays in advance of severe magnetic storms: Neutron monitor observations. In Proceedings of the 27th International Cosmic Ray Conference (ICRC 2001), Hamburg, Germany, 7–15 August 2001; Droege, W., Kunow, H., Scholer, M., Eds.; International Union of Pure and Applied Physics: Hamburg, Germany, 2001; Volume 9, pp. 3507–3510. [Google Scholar]
  58. Gosling, J.T.; Bame, S.J.; McComas, D.J.; Phillips, J.L. Coronal mass ejections and large geomagnetic storms. Geophys. Res. Lett. 1990, 17, 901. [Google Scholar] [CrossRef]
  59. Richardson, I.G.; Cane, H.V. Solar wind drivers of geomagnetic storms during more than four solar cycles. J. Space Weather Space Clim. 2012, 2, A01. [Google Scholar] [CrossRef]
  60. Calisto, M.; Usoskin, I.; Rozanov, E.; Peter, T. Influence of Galactic Cosmic Rays on atmospheric composition and dynamics. Atmos. Chem. Phys. 2011, 11, 4547–4556. [Google Scholar] [CrossRef]
  61. Armano, M.; Audley, H.; Baird, J.; Benella, S.; Binetruy, P.; Born, M.; Bortoluzzi, D.; Castelli, E.; Cavalleri, A.; Cesarini, A.; et al. Forbush Decreases and <2 Day GCR Flux Non-recurrent Variations Studied with LISA Pathfinder. Astrophys. J. 2019, 874, 167. [Google Scholar] [CrossRef]
  62. Papailiou, M.; Mavromichalaki, H.; Belov, A.; Eroshenko, E.; Yanke, V. Precursor Effects in Different Cases of Forbush Decreases. Sol. Phys. 2012, 276, 337–350. [Google Scholar] [CrossRef]
  63. Belov, A.V.; Abunina, M.A.; Abunin, A.A.; Eroshenko, E.A.; Oleneva, V.A.; Yanke, V.G. Cosmic ray vector anisotropy and local characteristics of the interplanetary medium. Geomagn. Aeron. 2017, 57, 389. [Google Scholar] [CrossRef]
  64. Papailiou, M.; Abunina, M.; Belov, A.; Eroshenko, E.; Yanke, V. Large Forbush Decreases and their Solar Sources: Features and Characteristics. Sol. Phys. 2020, 295, 164. [Google Scholar] [CrossRef]
  65. Papailiou, M.; Abunina, M.; Mavromichalaki, H.; Abunin, A.; Belov SGerontidou, M.; Belov, A.; Shlyk, N.; Yanke, V. Precursory signs of Large Forbush decreases: The criterion of anisotropy. Sol. Phys. 2024; submitted. [Google Scholar]
  66. Abunin, A.; Abunina, M.; Belov, A.; Eroshenko, E.; Oleneva, V.; Yanke, V. Forbush-decreases in 19th solar cycle. J. Phys. Conf. Ser. 2013, 409, 012165. [Google Scholar] [CrossRef]
  67. Belov, A.; Eroshenko, E.; Yanke, V.; Oleneva, V.; Abunin, A.; Abunina, M.; Papaioannou, A.; Mavromichalaki, H. The Global Survey Method Applied to Ground-level Cosmic Ray Measurements. Sol. Phys. 2018, 293, 68. [Google Scholar] [CrossRef]
  68. Hanslmeier, A. The SUN and Space Weather; Astrophysics and Space Science Library: Toledo, Spain, 2007; Volume 347. [Google Scholar]
Figure 1. (a) The asymptotic longitude–time diagram for the event on 14 February 1978 and (b) data on cosmic ray variations calculated by the GSM for the same event.
Figure 1. (a) The asymptotic longitude–time diagram for the event on 14 February 1978 and (b) data on cosmic ray variations calculated by the GSM for the same event.
Atmosphere 15 00742 g001
Figure 2. The asymptotic longitude–time diagram for the event on 6 August 1982.
Figure 2. The asymptotic longitude–time diagram for the event on 6 August 1982.
Atmosphere 15 00742 g002
Figure 3. The asymptotic longitude–time diagram for the event on (a) 9 January 1983 and (b) 24 November 2001.
Figure 3. The asymptotic longitude–time diagram for the event on (a) 9 January 1983 and (b) 24 November 2001.
Atmosphere 15 00742 g003aAtmosphere 15 00742 g003b
Figure 4. The asymptotic longitude–time diagram for the event on 23 April 2023.
Figure 4. The asymptotic longitude–time diagram for the event on 23 April 2023.
Atmosphere 15 00742 g004
Table 1. The FDs under examination and the corresponding magnitude and equatorial component of the first harmonic of cosmic ray anisotropy 1 h before the shock arrival.
Table 1. The FDs under examination and the corresponding magnitude and equatorial component of the first harmonic of cosmic ray anisotropy 1 h before the shock arrival.
a/aForbush
Decrease
Magnitude,
%
Axyb, %
123 March 19699.20.79
214 May 19698.70.65
329 September 19695.70.59
424 July 19705.60.43
514 December 19705.80.67
617 June 19727.90.44
714 February 197820.90.38
824 April 19795.30.67
919 December 198080.33
1025 July 19818.80.51
1110 October 19816.10.61
1213 October 198180.54
131 March 19827.20.64
146 August 19826.10.42
155 September 19825.30.63
1621 September 19828.50.57
1724 November 19827.40.45
189 January 19839.40.32
1921 February 19888.70.41
2013 March 198920.90.11
218 June 198960.30
2218 September 19898.20.69
2317 November 19895.90.08
2412 March 199070.63
2520 March 19905.60.18
2628 July 19905.30.52
2724 March 199122.10.60
284 June 19918.80.71
2918 August 19918.60.38
308 February 19925.80.32
3120 February 19925.10.23
329 May 199210.50.32
3322 August 19926.30.50
3418 February 19996.10.27
356 November 20008.70.37
3624 November 20019.80.67
3717 April 200270.63
3823 May 20027.50.34
3918 August 20025.10.76
4029 May 20036.90.47
4122 January 20049.40.35
427 November 20048.60.51
439 November 20048.60.41
4411 September 200513.20.41
4514 December 20069.60.74
4626 September 20115.10.22
4714 July 20129.60.53
4817 March 20155.60.43
497 September 20177.70.42
5023 April 20237.20.16
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Papailiou, M.-C.; Abunina, M.; Mavromichalaki, H.; Shlyk, N.; Belov, S.; Abunin, A.; Gerontidou, M.; Belov, A.; Yanke, V.; Triantou, A. Precursory Signs of Large Forbush Decreases in Relation to Cosmic Rays Equatorial Anisotropy Variation. Atmosphere 2024, 15, 742. https://0-doi-org.brum.beds.ac.uk/10.3390/atmos15070742

AMA Style

Papailiou M-C, Abunina M, Mavromichalaki H, Shlyk N, Belov S, Abunin A, Gerontidou M, Belov A, Yanke V, Triantou A. Precursory Signs of Large Forbush Decreases in Relation to Cosmic Rays Equatorial Anisotropy Variation. Atmosphere. 2024; 15(7):742. https://0-doi-org.brum.beds.ac.uk/10.3390/atmos15070742

Chicago/Turabian Style

Papailiou, Maria-Christina, Maria Abunina, Helen Mavromichalaki, Nataly Shlyk, Semyon Belov, Artem Abunin, Maria Gerontidou, Anatoly Belov, Victor Yanke, and Amalia Triantou. 2024. "Precursory Signs of Large Forbush Decreases in Relation to Cosmic Rays Equatorial Anisotropy Variation" Atmosphere 15, no. 7: 742. https://0-doi-org.brum.beds.ac.uk/10.3390/atmos15070742

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop