-
Large amplitude bidirectional anisotropy of cosmic-ray intensity observed with world-wide networks of ground-based neutron monitors and muon detectors in November, 2021
Authors:
K. Munakata,
M. Kozai,
C. Kato,
Y. Hayashi,
R. Kataoka,
A. Kadokura,
M. Tokumaru,
R. R. S. Mendonça,
E. Echer,
A. Dal Lago,
M. Rockenbach,
N. J. Schuch,
J. V. Bageston,
C. R. Braga,
H. K. Al Jassar,
M. M. Sharma,
M. L. Duldig,
J. E. Humble,
I. Sabbah,
P. Evenson,
P. -S. Mangeard,
T. Kuwabara,
D. Ruffolo,
A. Sáiz,
W. Mitthumsiri
, et al. (2 additional authors not shown)
Abstract:
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 fir…
▽ More
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.
△ Less
Submitted 13 September, 2022;
originally announced September 2022.
-
A Peculiar ICME Event in August 2018 Observed with the Global Muon Detector Network
Authors:
W. Kihara,
K. Munakata,
C. Kato,
R. Kataoka,
A. Kadokura,
S. Miyake,
M. Kozai,
T. Kuwabara,
M. Tokumaru,
R. R. S. Mendonça,
E. Echer,
A. Dal Lago,
M. Rockenbach,
N. J. Schuch,
J. V. Bageston,
C. R. Braga,
H. K. Al Jassar,
M. M. Sharma,
M. L. Duldig,
J. E. Humble,
P. Evenson,
I. Sabbah,
J. Kóta
Abstract:
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 acco…
▽ More
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.
△ Less
Submitted 19 January, 2021;
originally announced January 2021.
-
Cosmic ray short burst observed with the Global Muon Detector Network (GMDN) on June 22, 2015
Authors:
K. Munakata,
M. Kozai,
P. Evenson,
T. Kuwabara,
C. Kato,
M. Tokumaru,
M. Rockenbach,
A. Dal Lago,
R. R. S. Mendonca,
C. R. Braga,
N. J. Schuch,
H. K. Al Jassar,
M. M. Sharma,
M. L. Duldig,
J. E. Humble,
I. Sabbah,
J. Kota
Abstract:
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 appe…
▽ More
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.
△ Less
Submitted 26 June, 2018;
originally announced June 2018.
-
Average spatial distribution of cosmic rays behind the interplanetary shock -Global Muon Detector Network observations-
Authors:
M. Kozai,
K. Munakata,
C. Kato,
T. Kuwabara,
M. Rockenbach,
A. Dal Lago,
N. J. Schuch,
C. R. Braga,
R. R. S. Mendonça,
H. K. Al Jassar,
M. M. Sharma,
M. L. Duldig,
J. E. Humble,
P. Evenson,
I. Sabbah,
M. Tokumaru
Abstract:
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…
▽ More
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.
△ Less
Submitted 21 May, 2016;
originally announced May 2016.
-
First Very Low Frequency detection of short repeated bursts from magnetar SGR J1550-5418
Authors:
Y. T. Tanaka,
J. P. Raulin,
F. C. P. Bertoni,
P. R. Fagundes,
J. Chau,
N. J. Schuch,
M. Hayakawa,
Y. Hobara,
T. Terasawa,
T. Takahashi
Abstract:
We report on the first detection of ionospheric disturbances caused by short repeated gamma-ray bursts from the magnetar SGR J1550-5418. Very low frequency (VLF) radio wave data obtained in South America clearly show sudden amplitude and phase changes at the corresponding times of eight SGR bursts. Maximum amplitude and phase changes of the VLF signals appear to be correlated with the gamma-ray fl…
▽ More
We report on the first detection of ionospheric disturbances caused by short repeated gamma-ray bursts from the magnetar SGR J1550-5418. Very low frequency (VLF) radio wave data obtained in South America clearly show sudden amplitude and phase changes at the corresponding times of eight SGR bursts. Maximum amplitude and phase changes of the VLF signals appear to be correlated with the gamma-ray fluence. On the other hand, VLF recovery timescales do not show any significant correlation with the fluence, possibly suggesting that the bursts' spectra are not similar to each other. In summary, the Earth's ionosphere can be used as a very large gamma-ray detector and the VLF observations provide us with a new method to monitor high energy astrophysical phenomena without interruption such as Earth Occultation.
△ Less
Submitted 21 August, 2010;
originally announced August 2010.
-
Drift effects and the average features of cosmic ray density gradient in CIRs during successive two solar minimum periods
Authors:
A. Fushishita,
Y. Okazaki,
T. Narumi,
C. Kato,
S. Yasue,
T. Kuwabara,
J. W. Bieber,
P. Evenson,
M. R. Da Silva,
A. Dal Lago,
N. J. Schuch,
M. Tokumaru,
M. L. Duldig,
J. E. Humble,
I. Sabbah,
J. Kota,
And K. Munakata
Abstract:
We deduce on hourly basis the spatial gradient of the cosmic ray density in three dimensions from the directional anisotropy of high-energy (~50 GeV) galactic cosmic ray (GCR) intensity observed with a global network of muon detectors on the Earth's surface. By analyzing the average features of the gradient in the corotational interaction regions (CIRs) recorded in successive two solar activity…
▽ More
We deduce on hourly basis the spatial gradient of the cosmic ray density in three dimensions from the directional anisotropy of high-energy (~50 GeV) galactic cosmic ray (GCR) intensity observed with a global network of muon detectors on the Earth's surface. By analyzing the average features of the gradient in the corotational interaction regions (CIRs) recorded in successive two solar activity minimum periods, we find that the observed latitudinal gradient (Gz) changes its sign from negative to positive on the Earth's heliospheric current sheet (HCS) crossing from the northern to the southern hemisphere in A<0 epoch, while it changes from positive to negative in A>0 epoch. This is in accordance with the drift prediction. We also find a negative enhancement in Gx after the HCS crossing in both A<0 and A>0 epochs, but not in Gy. This asymmetrical feature of Gx and Gy indicates significant contributions from the parallel and perpendicular diffusions to the the gradient in CIRs in addition to the contribution from the drift effect.
△ Less
Submitted 24 December, 2009;
originally announced December 2009.
-
Precursors of the Forbush Decrease on December 14, 2006 observed with the Global Muon Detector Network (GMDN)
Authors:
A. Fushishita,
T. Kuwabara,
C. Kato,
S. Yasue,
J. W. Bieber,
P. Evenson,
M. R. Da Silva,
A. Dal Lago,
N. J. Schuch,
M. Tokumaru,
M. L. Duldig,
J. E. Humble,
I. Sabbah,
H. K. Al Jassar,
M. M. Sharma,
K. Munakata
Abstract:
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 prec…
▽ More
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.
△ Less
Submitted 1 February, 2010; v1 submitted 5 September, 2009;
originally announced September 2009.
-
Drift effects and the cosmic ray density gradient in a solar rotation period: First observation with the Global Muon Detector Network (GMDN)
Authors:
Y. Okazaki,
A. Fushishita,
T. Narumi,
C. Kato,
S. Yasue,
T. Kuwabara,
J. W. Bieber,
P. Evenson,
M. R. Da Silva,
A. Dal Lago,
N. J. Schuch,
Z. Fujii,
M. L. Duldig,
J. E. Humble,
I. Sabbah,
J. Kóta,
K. Munakata
Abstract:
We present for the first time hourly variations of the spatial density gradient of 50 GeV cosmic rays within a sample solar rotation period in 2006. By inversely solving the transport equation, including diffusion, we deduce the gradient from the anisotropy that is derived from the observation made by the Global Muon Detector Network (GMDN). The anisotropy obtained by applying a new analysis met…
▽ More
We present for the first time hourly variations of the spatial density gradient of 50 GeV cosmic rays within a sample solar rotation period in 2006. By inversely solving the transport equation, including diffusion, we deduce the gradient from the anisotropy that is derived from the observation made by the Global Muon Detector Network (GMDN). The anisotropy obtained by applying a new analysis method to the GMDN data is precise and free from atmospheric temperature effects on the muon count rate recorded by ground based detectors. We find the derived north-south gradient perpendicular to the ecliptic plane is oriented toward the Helioshperic Current Sheet (HCS) (i.e. southward in the toward sector of the Interplanetary Magnetic Field (IMF) and northward in the away sector). The orientation of the gradient component parallel to the ecliptic plane remains similar in both sectors with an enhancement of its magnitude seen after the Earth crosses the HCS. These temporal features are interpreted in terms of a local maximum of the cosmic ray density at the HCS. This is consistent with the prediction of the drift model for the $A<0$ epoch. By comparing the observed gradient with the numerical prediction of a simple drift model, we conclude that particle drifts in the large-scale magnetic field play an important role in organizing the density gradient, at least in the present $A<0$ epoch. We also found that corotating interaction regions did not have such a notable effect. Observations with the GMDN provide us with a new tool for investigating cosmic ray transport in the IMF.
△ Less
Submitted 16 February, 2008;
originally announced February 2008.
