Search

This paper presents the global analysis of two extended decreases of the galactic cosmic ray intensity observed by world-wide networks of ground-based detectors in 2012. This analysis is capable of separately deriving the cosmic ray density (or omnidirectional intensity) and anisotropy each as a function of time and rigidity. A simple diffusion model along the spiral field line between Earth and a cosmic-ray barrier indicates the long duration of these events resulting from about 190$^\circ$ eastern extension of a barrier such as an IP-shock followed by the sheath region and/or the corotating interaction region (CIR). It is suggested that the coronal mass ejection merging and compressing the preexisting CIR at its flank can produce such the extended barrier. The derived rigidity spectra of the density and anisotropy both vary in time during each event period. In particular we find that the temporal feature of the ``phantom Forbush decrease'' reported in an analyzed period is dependent on rigidity, looking quite different at different rigidities. From these rigidity spectra of the density and anisotropy, we derive the rigidity spectrum of the average parallel mean-free-path of pitch angle scattering along the spiral field line and infer the power spectrum of the magnetic fluctuation and its temporal variation. Possible physical cause of the strong rigidity dependence of the ``phantom Forbush decrease'' is also discussed. These results demonstrate the high-energy cosmic rays observed at Earth responding to remote space weather., Comment: Accepted for publication in ApJ

We analyze the cosmic-ray variations during a significant Forbush decrease observed with world-wide networks of ground-based neutron monitors and muon detectors during November 3-5, 2021. Utilizing the difference between primary cosmic-ray rigidities monitored by neutron monitors and muon detectors, we deduce the rigidity spectra of the cosmic-ray density (or omnidirectional intensity) and the first- and second-order anisotropies separately, for each hour of data. A clear two-step decrease is seen in the cosmic-ray density with the first $\sim2\%$ decrease after the interplanetary shock arrival followed by the second $\sim5\%$ decrease inside the magnetic flux rope (MFR) at 15 GV. Most strikingly, a large bidirectional streaming along the magnetic field is observed in the MFR with a peak amplitude of $\sim5\%$ at 15 GV which is comparable to the total density decrease inside the MFR. The bidirectional streaming could be explained by adiabatic deceleration and/or focusing in the expanding MFR, which have stronger effects for pitch angles near 90$^\circ$, or by selective entry of GCRs along a leg of the MFR. The peak anisotropy and density depression in the flux rope both decrease with increasing rigidity. The spectra vary dynamically indicating that the temporal variations of density and anisotropy appear different in neutron monitor and muon detector data., Comment: Accepted for publication in the Astrophysical Journal

We demonstrate that global observations of high-energy cosmic rays contribute to understanding unique characteristics of a large-scale magnetic flux rope causing a magnetic storm in August 2018. Following a weak interplanetary shock on 25 August 2018, a magnetic flux rope caused an unexpectedly large geomagnetic storm. It is likely that this event became geoeffective because the flux rope was accompanied by a corotating interaction region and compressed by high-speed solar wind following the flux rope. In fact, a Forbush decrease was observed in cosmic-ray data inside the flux rope as expected, and a significant cosmic-ray density increase exceeding the unmodulated level before the shock was also observed near the trailing edge of the flux rope. The cosmic-ray density increase can be interpreted in terms of the adiabatic heating of cosmic rays near the trailing edge of the flux rope, as the corotating interaction region prevents free expansion of the flux rope and results in the compression near the trailing edge. A northeast-directed spatial gradient in the cosmic-ray density was also derived during the cosmic-ray density increase, suggesting that the center of the heating near the trailing edge is located northeast of Earth. This is one of the best examples demonstrating that the observation of high-energy cosmic rays provides us with information that can only be derived from the cosmic ray measurements to observationally constrain the three-dimensional macroscopic picture of the interaction between coronal mass ejections and the ambient solar wind, which is essential for prediction of large magnetic storms., Comment: 19 pages, 3 figures, accepted for publication in the Space Weather

Cosmic rays are charged particles whose flux observed at Earth shows temporal variations related to space weather phenomena and may be an important tool to study them. The cosmic ray intensity recorded with ground-based detectors also shows temporal variations arising from atmospheric variations. In the case of muon detectors, the main atmospheric effects are related to pressure and temperature changes. In this work, we analyze both effects using data recorded by the Global Muon Detector Network (GMDN), consisting of four multidirectional muon detectors at different locations, in the period between 2007 and 2016. For each GMDN directional channel, we obtain coefficients that describe the pressure and temperature effects. We then analyze how these coefficients can be related to the geomagnetic cutoff rigidity and zenith angle associated with cosmic-ray particles observed by each channel. In the pressure effect analysis, we found that the observed barometric coefficients show a very clear logarithmic correlation with the cutoff rigidity divided by the zenith angle cosine. On the other hand, the temperature coefficients show a good logarithmic correlation with the product of the cutoff and zenith angle cosine after adding a term proportional to the sine of geographical latitude of the observation site. This additional term implies that the temperature effect measured in the northern hemisphere detectors is stronger than that observed in the southern hemisphere. The physical origin of this term and of the good correlations found in this analysis should be studied in detail in future works.

We analyze the short cosmic ray intensity increase ("cosmic ray burst": CRB) on June 22, 2015 utilizing a global network of muon detectors and derive the global anisotropy of cosmic ray intensity and the density (i.e. the omnidirectional intensity) with 10-minute time resolution. We find that the CRB was caused by a local density maximum and an enhanced anisotropy of cosmic rays both of which appeared in association with Earth's crossing of the heliospheric current sheet (HCS). This enhanced anisotropy was normal to the HCS and consistent with a diamagnetic drift arising from the spatial gradient of cosmic ray density, which indicates that cosmic rays were drifting along the HCS from the north of Earth. We also find a significant anisotropy along the HCS, lasting a few hours after the HCS crossing, indicating that cosmic rays penetrated into the inner heliosphere along the HCS. Based on the latest geomagnetic field model, we quantitatively evaluate the reduction of the geomagnetic cut-off rigidity and the variation of the asymptotic viewing direction of cosmic rays due to a major geomagnetic storm which occurred during the CRB and conclude that the CRB is not caused by the geomagnetic storm, but by a rapid change in the cosmic ray anisotropy and density outside the magnetosphere., Comment: accepted for the publication in the Astrophysical Journal

We analyze the galactic cosmic ray (GCR) density and its spatial gradient in Forbush Decreases (FDs) observed with the Global Muon Detector Network (GMDN) and neutron monitors (NMs). By superposing the GCR density and density gradient observed in FDs following 45 interplanetary shocks (IP-shocks), each associated with an identified eruption on the sun, we infer the average spatial distribution of GCRs behind IP-shocks. We find two distinct modulations of GCR density in FDs, one in the magnetic sheath and the other in the coronal mass ejection (CME) behind the sheath. The density modulation in the sheath is dominant in the western flank of the shock, while the modulation in the CME ejecta stands out in the eastern flank. This east-west asymmetry is more prominent in GMDN data responding to $\sim$ 60 GV GCRs than in NM data responding to $\sim$ 10 GV GCRs, because of softer rigidity spectrum of the modulation in the CME ejecta than in the sheath. The GSE-y component of the density gradient, $G_y$ shows a negative (positive) enhancement in FDs caused by the eastern (western) eruptions, while $G_z$ shows a negative (positive) enhancement in FDs by the northern (southern) eruptions. This implies the GCR density minimum being located behind the central flank of IP-shock and propagating radially outward from location of the solar eruption. We also confirmed the average $G_z$ changing its sign above and below the heliospheric current sheet, in accord with the prediction of the drift model for the large-scale GCR transport in the heliosphere., Comment: 54 pages, 11 figures, 1 table, accepted for publication in the Astrophysical Journal

Energy coupling between the solar wind and the Earth's magnetosphere can affect the electron population in the outer radiation belt. However, the precise role of different internal and external mechanisms that leads to changes of the relativistic electron population is not entirely known. This paper describes how ultralow frequency (ULF) wave activity during the passage of Alfvenic solar wind streams contributes to the global recovery of the relativistic electron population in the outer radiation belt. To investigate the contribution of the ULF waves, we searched the Van Allen Probes data for a period in which we can clearly distinguish the enhancement of electron fluxes from the background. We found that the global recovery that started on 22 September 2014, which coincides with the corotating interaction region preceding a high‐speed stream and the occurrence of persistent substorm activity, provides an excellent scenario to explore the contribution of ULF waves. To support our analyses, we employed ground‐ and space‐based observational data and global magnetohydrodynamic simulations and calculated the ULF wave radial diffusion coefficients employing an empirical model. Observations show a gradual increase of electron fluxes in the outer radiation belt and a concomitant enhancement of ULF activity that spreads from higher to lower L‐shells. Magnetohydrodynamic simulation results agree with observed ULF wave activity in the magnetotail, which leads to both fast and Alfven modes in the magnetospheric nightside sector. The observations agree with the empirical model and are confirmed by phase space density calculations for this global recovery period.

Analysis of particle pitch angle distributions (PADs) has been used as a means to comprehend a multitude of different physical mechanisms that lead to flux variations in the Van Allen belts and also to particle precipitation into the upper atmosphere. In this work we developed a neural network-based data clustering methodology that automatically identifies distinct PAD types in an unsupervised way using particle flux data. One can promptly identify and locate three well-known PAD types in both time and radial distance, namely, 90deg peaked, butterfly, and flattop distributions. In order to illustrate the applicability of our methodology, we used relativistic electron flux data from the whole month of November 2014, acquired from the Relativistic Electron-Proton Telescope instrument on board the Van Allen Probes, but it is emphasized that our approach can also be used with multiplatform spacecraft data. Our PAD classification results are in reasonably good agreement with those obtained by standard statistical fitting algorithms. The proposed methodology has a potential use for Van Allen belt's monitoring.

Magnetopause shadowing and wave-particle interactions are recognized as the two primary mechanisms for losses of electrons from the outer radiation belt. We investigate these mechanisms, sing satellite observations both in interplanetary space and within the magnetosphere and particle drift modeling. Two interplanetary shocks sheaths impinged upon the magnetopause causing a relativistic electron flux dropout. The magnetic cloud (C) and interplanetary structure sunward of the MC had primarily northward magnetic field, perhaps leading to a concomitant lack of substorm activity and a 10 day long quiescent period. The arrival of two shocks caused an unusual electron flux dropout. Test-particle simulations have shown 2 to 5 MeV energy, equatorially mirroring electrons with initial values of L 5.5can be lost to the magnetosheath via magnetopause shadowing alone. For electron losses at lower L-shells, coherent chorus wave-driven pitch angle scattering and ULF wave-driven radial transport have been shownto be viable mechanisms.