Search

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 report on the 22-23 June 2015 geomagnetic storm. There has been a shortage of intense geomagnetic storms during the current solar cycle 24 in relation to the previous cycle. This situation changed after mid-June 2015 when one of the biggest solar active regions (AR 2371) of current solar cycle 24, close to the central meridian produced several coronal mass ejections (CMEs) associated with M-class flares. The CMEs impact on the Earth's magnetosphere resulted in a moderately-severe G4-class geomagnetic storm on 22-23 June 2015 and a G2 (moderate) geomagnetic storms on 24 June. The G4 solstice storm was the second biggest (so far) geomagnetic storms of cycle 24. We highlight the ground level observations made by New-Tupi, Muonca and the CARPET El Leoncito cosmic ray detectors that are located within the South Atlantic Anomaly (SAA) region. These observations are studied in correlation with data obtained by space-borne detectors and other ground-based experiments. The CME designations are from the Computer Aided CME Tracking (CACTus) automated catalog. As expected, Forbush Decreases (FD) associated with the passing CMEs were recorded by these detectors. We noticed a peculiar feature linked to a severe geomagnetic storm event. The 21 June 2015 CME 0091 was likely associated with the 22 June summer solstice FD event. The angular width of CME 0091 was very narrow and measured 56 degrees seen from Earth. In most cases, only CME halos and partial halos, lead to severe geomagnetic storms. We performed a cross-check analysis of the FD events detected during the rise phase of the current solar cycle 24, the geomagnetic parameters, and the CACTus CME catalog. Our study suggests that narrow angular-width CMEs that erupt in the western region of the ecliptic plane can lead to moderate and severe geomagnetic storms., Comment: 20 pages, 14 figures, accepted in Solar Physics

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