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.

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.

We analyze the precursor of a Forbush Decrease (FD) observed with the Global Muon Detector Network on December 14, 2006. An intense geomagnetic storm is also recorded during this FD with the peak Kp index of 8+. By using the "two-dimensional map" of the cosmic ray intensity produced after removing the contribution from the diurnal anisotropy, we succeed in extracting clear signatures of the precursor. A striking feature of this event is that a weak loss-cone signature is first recorded more than a day prior to the Storm Sudden Commencement (SSC) onset. This suggests that the loss-cone precursor appeared only 7 hours after the Coronal Mass Ejection (CME) eruption from the Sun, when the Interplanetary (IP) shock driven by the Interplanetary Coronal Mass Ejection (ICME) located at 0.4 AU from the Sun. We find the precursor being successively observed with multiple detectors in the network according to the Earth's spin and confirmed that the precursor continuously exists in space. The long lead time (15.6 hours) of this precursor which is almost twice the typical value indicates that the IMF was more quiet in this event than a typical power spectrum assumed for the IMF turbulence. The amplitude (-6.45 %) of the loss-cone anisotropy at the SSC onset is more than twice the FD size, indicating that the maximum intensity depression behind the IP shock is much larger than the FD size recorded at the Earth in this event. We also find the excess intensity from the sunward IMF direction clearly observed during ~10 hours preceding the SSC onset. It is shown that this excess intensity is consistent with the measurement of the particles accelerated by the head-on collisions with the approaching shock. This is the first detailed observation of the precursor due to the shock reflected particles with muon detectors.