Showing 1–18 of 18 results for author: Kota, J

Role of magnetic arcades in explaining the puzzle of the gamma-ray emission from the solar disk

Authors: Eleonora Puzzoni, Federico Fraschetti, József Kóta, Joe Giacalone

Abstract: The interpretation of gamma-ray emission originating from the solar disk ($0.5^\circ$ in angular size) as due to the interaction of Galactic Cosmic Rays (GCRs) with the solar atmosphere has remained a central challenge in solar physics. After the seminal work by Seckel, Stanev, and Gaisser (SSG91) based on GCRs magnetic mirroring, discrepancies between models and observations persist, indicating t… ▽ More The interpretation of gamma-ray emission originating from the solar disk ($0.5^\circ$ in angular size) as due to the interaction of Galactic Cosmic Rays (GCRs) with the solar atmosphere has remained a central challenge in solar physics. After the seminal work by Seckel, Stanev, and Gaisser (SSG91) based on GCRs magnetic mirroring, discrepancies between models and observations persist, indicating the need for a novel approach. The present work focuses on exploring the impact of a closed magnetic field geometry in the low photosphere on the observed gamma-ray flux. We track numerically with the PLUTO code the trajectories of test-particle protons within a static $\sim 20$ Mm scale height magnetic arcade adjacent to jets. By making use of numerical vertical density profiles we inject particles at distinct chromospheric/photospheric altitudes, mimicking the migration of GCRs from neighboring flux tubes into closed arcades. Remarkably, our model reproduces a flat gamma-ray spectrum below $\sim 33$ GeV, a nearly-isotropic emission at $\sim 10$ GeV, both consistent with Fermi-LAT observations, and a near-limb emission at $\sim 1$ TeV. Our model can also reproduce the flux-drop detected by HAWC ($\sim 1$ TeV). Finally, we argue that the spectral dip observed at $\sim$ 40 GeV may result from the flux suppression at low energy due to the cross-field diffusion, which would produce a cutoff. These findings underscore the pivotal role of closed magnetic field structures in shaping the solar disk gamma-ray emission. △ Less

Submitted 23 July, 2024; originally announced July 2024.

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.

Comments: Accepted for publication in the Astrophysical Journal

On the Energy Dependence of Galactic Cosmic Ray Anisotropies in the Very Local Interstellar Medium

Authors: Romina Nikoukar, Matthew E. Hill, Lawrence Brown, Stamatios M. Krimigis, Robert B. Decker, Konstantinos Dialynas, Jozsef Kota, Edmond C. Roelof, Scott Lasley, Douglas C. Hamilton, Vladimir Florinski, Joe Giacalone, John Richardson, Merav Opher

Abstract: We report on the energy dependence of galactic cosmic rays (GCRs) in the very local interstellar medium (VLISM) as measured by the Low Energy Charged Particle (LECP) instrument on the Voyager 1 (V1) spacecraft. The LECP instrument includes a dual-ended solid state detector particle telescope mechanically scanning through 360 deg across eight equally-spaced angular sectors. As reported previously,… ▽ More We report on the energy dependence of galactic cosmic rays (GCRs) in the very local interstellar medium (VLISM) as measured by the Low Energy Charged Particle (LECP) instrument on the Voyager 1 (V1) spacecraft. The LECP instrument includes a dual-ended solid state detector particle telescope mechanically scanning through 360 deg across eight equally-spaced angular sectors. As reported previously, LECP measurements showed a dramatic increase in GCR intensities for all sectors of the >=211 MeV count rate (CH31) at the V1 heliopause (HP) crossing in 2012, however, since then the count rate data have demonstrated systematic episodes of intensity decrease for particles around 90° pitch angle. To shed light on the energy dependence of these GCR anisotropies over a wide range of energies, we use V1 LECP count rate and pulse height analyzer (PHA) data from >=211 MeV channel together with lower energy LECP channels. Our analysis shows that while GCR anisotropies are present over a wide range of energies, there is a decreasing trend in the amplitude of second-order anisotropy with increasing energy during anisotropy episodes. A stronger pitch-angle scattering at the higher velocities is argued as a potential cause for this energy dependence. A possible cause for this velocity dependence arising from weak rigidity dependence of the scattering mean free path and resulting velocity-dominated scattering rate is discussed. This interpretation is consistent with a recently reported lack of corresponding GCR electron anisotropies. △ Less

Submitted 19 January, 2022; originally announced January 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.

Comments: 19 pages, 3 figures, accepted for publication in the Space Weather

No Stagnation Region Before the Heliopause at Voyager 1? Inferences From New Voyager 2 Results

Authors: A. C. Cummings, E. C. Stone, J. D. Richardson, B. C. Heikkila, N. Lal, J. Kóta

Abstract: We present anisotropy results for anomalous cosmic-ray (ACR) protons in the energy range $\sim$0.5-35 MeV from Cosmic Ray Subsytem (CRS) data collected during calibration roll maneuvers for the magnetometer instrument when Voyager 2 (V2) was in the inner heliosheath. We use a new technique to derive for the first time the radial component of the anisotropy vector from CRS data. We find that the CR… ▽ More We present anisotropy results for anomalous cosmic-ray (ACR) protons in the energy range $\sim$0.5-35 MeV from Cosmic Ray Subsytem (CRS) data collected during calibration roll maneuvers for the magnetometer instrument when Voyager 2 (V2) was in the inner heliosheath. We use a new technique to derive for the first time the radial component of the anisotropy vector from CRS data. We find that the CRS-derived radial solar wind speeds, when converted from the radial components of the anisotropy vectors via the Compton-Getting (C-G) effect, generally agree with those similarly-derived speeds from the Low-Energy Charged Particle experiment using 28-43 keV data. However, they often differ significantly from the radial solar wind speeds measured directly by the Plasma Science (PLS) instrument. There are both periods when the C-G-derived radial solar wind speeds are significantly higher than those measured by PLS and times when they are significantly lower. The differences are not expected nor explained, but it appears that after a few years in the heliosheath the V2 radial solar wind speeds derived from the C-G method underestimate the true speeds as the spacecraft approaches the heliopause. We discuss the implications of this observation for the stagnation region reported along the Voyager 1 trajectory as it approached the heliopause inferred using the C-G method. △ Less

Submitted 25 November, 2020; originally announced November 2020.

High Resolution Finite Volume Method for Kinetic Equations with Poisson Brackets

Authors: Igor V. Sokolov, Haomin Sun, Gabor Toth, Zhenguang Huang, Valeriy Tenishev, Lulu Zhao, Jozsef Kota, Ofer Cohen, Tamas Gombosi

Abstract: Simulation of plasmas in electromagnetic fields requires numerical solution of a kinetic equation that describes the time evolution of the particle distribution function. In this paper we propose a finite volume scheme based on integral relation for Poisson brackets to solve the Liouville equation, the most fundamental kinetic equation. The proposed scheme conserves the number of particles, mainta… ▽ More Simulation of plasmas in electromagnetic fields requires numerical solution of a kinetic equation that describes the time evolution of the particle distribution function. In this paper we propose a finite volume scheme based on integral relation for Poisson brackets to solve the Liouville equation, the most fundamental kinetic equation. The proposed scheme conserves the number of particles, maintains the total-variation-diminishing (TVD) property, and provides high-quality numerical results. Other types of kinetic equations may be also formulated in terms of Poisson brackets and solved with the proposed method including the transport equations describing the acceleration and propagation of Solar Energetic Particles (SEPs), which is of practical importance, since the high energy SEPs produce radiation hazards. The proposed scheme is demonstrated to be accurate and efficient, which makes it applicable to global simulation systems analyzing space weather. △ Less

Submitted 26 November, 2022; v1 submitted 25 October, 2019; originally announced October 2019.

Comments: 50 pages, 9 figures

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.

Comments: accepted for the publication in the Astrophysical Journal

Journal ref: Astrophys. J. 862:170 2018

On the causes and mechanisms of the long-term variations in the GCR characteristics

Authors: M. Krainev, J. Kóta, M. S. Potgieter

Abstract: We argue that the degree of understanding the causes and mechanisms of the long-term variations (11-year and 22-year) in the galactic cosmic rays (GCR) characteristics is still insufficient and to improve it we need new approaches and methods. For the time being there is a long-lasting controversy on how these long-term variations, observed for more than 50 years in the inner heliosphere, are form… ▽ More We argue that the degree of understanding the causes and mechanisms of the long-term variations (11-year and 22-year) in the galactic cosmic rays (GCR) characteristics is still insufficient and to improve it we need new approaches and methods. For the time being there is a long-lasting controversy on how these long-term variations, observed for more than 50 years in the inner heliosphere, are formed. It is widely believed that the 11-year variation is due entirely to the toroidal branch of solar activity (the area and number of sunspots, the strength of the heliospheric magnetic field etc) because of the diffusion, convection and adiabatic energy loss, while the much smaller 22-year variation is caused by the particle drifts connected with the poloidal branch of solar activity (the high-latitude solar magnetic fields). At the same time, both past and more recent numerical simulations indicate that the contribution of particle drifts could be significant for both 22- and 11-year variations in the GCR intensities. However, even those who agree on the significant influence of drifts appear to have different perceptions on the mechanisms of this influence. In this paper, we present an analysis of the possible causes of the first point of view (small role of drifts in the 11-year GCR variation) and the reasons why one can expect the significant contri- bution of the processes connected with the poloidal branch of solar activity in both types of the long-term variations of the GCR characteristics. Then we briefly discuss some numerical methods suggested in the past and recently and the approaches and perspectives for the sought-for methods are considered. △ Less

Submitted 2 September, 2015; originally announced September 2015.

Comments: 34th International Cosmic Ray Conference, The Hague, The Netherlands, 30 July - 6 August, 2015, paper ID 198, 8 pages. in Proceedings of Science, 2015

Long term variation of the solar diurnal anisotropy of galactic cosmic rays observed with the Nagoya multi-directional muon detector

Authors: K. Munakata, M. Kozai, C. Kato, J. Kota

Abstract: We analyze the three dimensional anisotropy of the galactic cosmic ray (GCR) intensities observed independently with a muon detector at Nagoya in Japan and neutron monitors over four solar activity cycles. We clearly see the phase of the free-space diurnal anisotropy shifting toward earlier hours around solar activity minima in A>0 epochs, due to the reduced anisotropy component parallel to the me… ▽ More We analyze the three dimensional anisotropy of the galactic cosmic ray (GCR) intensities observed independently with a muon detector at Nagoya in Japan and neutron monitors over four solar activity cycles. We clearly see the phase of the free-space diurnal anisotropy shifting toward earlier hours around solar activity minima in A>0 epochs, due to the reduced anisotropy component parallel to the mean magnetic field. The average parallel component is consistent with a rigidity independent spectrum, while the perpendicular component increases with GCR rigidity. We suggest that this harder spectrum of the perpendicular component is due to contribution from the drift streaming. We find that the bidirectional latitudinal density gradient is positive in A>0 epoch, while it is negative in A<0 epoch, in accord with the drift model prediction. The radial density gradient, on the other hand, varies with ~11-year cycle with maxima (minima) in solar maximum (minimum) periods, but there is no significant difference seen between average radial gradients in A>0 and A<0 epochs. The average parallel mean free path is larger in A<0 than in A>0. We also find, however, that parallel mean free path (radial gradient) appears to persistently increase (decreasing) in the last three cycles of weakening solar activity. We suggest that simple differences between these parameters in A>0 and A<0 epochs are seriously biased by these long-term trends. △ Less

Submitted 23 June, 2014; v1 submitted 7 April, 2014; originally announced April 2014.

Comments: accepted for the publication in the Astrophysical Journal

Journal ref: Astrophys. J. 791:22 2014

The Cosmic Ray Intensity Near the Archean Earth

Authors: O. Cohen, J. J. Drake, J. Kota

Abstract: We employ three-dimensional state of the art magnetohydrodynamic models of the early solar wind and heliosphere and a two-dimensional model for cosmic ray transport to investigate the cosmic ray spectrum and flux near the Archean Earth. We assess how sensitive the cosmic ray spectrum is to changes in the sunspot placement and magnetic field strength, the large scale dipole magnetic field strength,… ▽ More We employ three-dimensional state of the art magnetohydrodynamic models of the early solar wind and heliosphere and a two-dimensional model for cosmic ray transport to investigate the cosmic ray spectrum and flux near the Archean Earth. We assess how sensitive the cosmic ray spectrum is to changes in the sunspot placement and magnetic field strength, the large scale dipole magnetic field strength, the wind ram pressure, and the Sun's rotation period. Overall, our results confirm earlier work that suggested the Archean Earth would have experienced a greatly reduced cosmic ray flux than is the case today. The cosmic ray reduction for the early Sun is mainly due to the shorter solar rotation period and tighter winding of the Parker spiral, and to the different surface distribution of the more active solar magnetic field. These effects lead to a global reduction of the cosmic ray flux at 1AU by up to two orders of magnitude or more. Variations in the sunspot magnetic field have more effect on the flux than variations in the dipole field component. The wind ram pressure affects the cosmic ray flux through its influence on the size of the heliosphere via the pressure balance with the ambient interstellar medium. Variations in the interstellar medium pressure experienced by the solar system in orbit through Galaxy could lead to order of magnitude changes in the cosmic ray flux at Earth on timescales of a few million years. △ Less

Submitted 26 September, 2012; originally announced September 2012.

Comments: 7 pages, 4 figures,accepted to ApJ

Particle acceleration by collisionless shocks containing large-scale magnetic-field variations

Authors: F. Guo, J. R. Jokipii, J. Kota

Abstract: Diffusive shock acceleration at collisionless shocks is thought to be the source of many of the energetic particles observed in space. Large-scale spatial variations of the magnetic field has been shown to be important in understanding observations. The effects are complex, so here we consider a simple, illustrative model. Here, we solve numerically the Parker transport equation for a shock in the… ▽ More Diffusive shock acceleration at collisionless shocks is thought to be the source of many of the energetic particles observed in space. Large-scale spatial variations of the magnetic field has been shown to be important in understanding observations. The effects are complex, so here we consider a simple, illustrative model. Here, we solve numerically the Parker transport equation for a shock in the presence of large-scale sinusoidal magnetic-field variations. We demonstrate that the familiar planar-shock results can be significantly altered as a consequence of large-scale, meandering magnetic lines of force. Because perpendicular diffusion coefficient $κ_\perp$ is generally much smaller than parallel diffusion coefficient $κ_\parallel$, the energetic charged particles are trapped and preferentially accelerated along the shock front in the regions where the connection points of magnetic field lines intersecting the shock surface converge, and thus create the "hot spots" of the accelerated particles. For the regions where the connection points separate from each other, the acceleration to high energies will be suppressed. Further, the particles diffuse away from the "hot spot" regions and modify the spectra of downstream particle distribution. These features are qualitatively similar to the recent Voyager's observation in the Heliosheath. These results are potentially important for particle acceleration at shocks propagating in turbulent magnetized plasmas as well as those which contain large-scale nonplanar structures. Examples include anomalous cosmic rays accelerated by the solar wind termination shock, energetic particles observed in propagating heliospheric shocks, and galactic cosmic rays accelerated by supernova blast waves, etc. △ Less

Submitted 2 October, 2010; v1 submitted 28 September, 2010; originally announced September 2010.

Comments: accepted to ApJ

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.

Comments: 11 pages, 1 table, 5 figures, Proc. of AOGS 2008, accepted for the publication on the Advances in Geoscience, 2010

Journal ref: Advances in Geosciences 21 199 2008

Solar cycle dependence of the diurnal anisotropy of 0.6 TeV cosmic ray intensity observed with the Matsushiro underground muon detector

Authors: K. Munakata, Y. Mizoguchi, C. Kato, S. Yasue, S. Mori, M. Takita, J. Kota

Abstract: We analyze the temporal variation of the diurnal anisotropy of sub-TeV cosmic ray intensity observed with the Matsushiro (Japan) underground muon detector over two full solar activity cycles in 1985-2008. The average sidereal amplitude over the entire period is 0.034+-0.003 %, which is roughly one third of the amplitude reported from AS and deep-underground muon experiments monitoring multi-TeV… ▽ More We analyze the temporal variation of the diurnal anisotropy of sub-TeV cosmic ray intensity observed with the Matsushiro (Japan) underground muon detector over two full solar activity cycles in 1985-2008. The average sidereal amplitude over the entire period is 0.034+-0.003 %, which is roughly one third of the amplitude reported from AS and deep-underground muon experiments monitoring multi-TeV GCR intensity suggesting a significant attenuation of the anisotropy due to the solar modulation. We find, on the other hand, only weak correlations either with the solar activity- or magnetic-cycles. We examine the temporal variation of the "single-band valley depth" (SBVD) quoted by the Milagro experiment and, by contrast with recent Milagro reports, we find no steady increase in the Matsushiro observations in a 7-year period between 2000 and 2007. We suggest, therefore, that the steady increase of the SBVD reported by the Milagro experiment is not caused by the decreasing solar modulation in the declining phase of the 23rd solar activity cycle. △ Less

Submitted 24 January, 2010; v1 submitted 5 November, 2009; originally announced November 2009.

Comments: submitted to the Astrophysical Journal

Journal ref: Astrophys.J.712:1100-1106,2010

The sidereal anisotropy of multi-TeV cosmic rays in an expanding Local Interstellar Cloud

Authors: Y. Mizoguchi, K. Munakata, M. Takita, J. Kota

Abstract: The sidereal anisotropy of galactic cosmic ray (GCR) intensity observed with the Tibet Air Shower (AS) experiment still awaits theoretical interpretation. The observed global feature of the anisotropy is well reproduced by a superposition of the bi-directional and uni-directional flows (BDF and UDF, respectively) of GCRs. If the orientation of the deduced BDF represents the orientation of the lo… ▽ More The sidereal anisotropy of galactic cosmic ray (GCR) intensity observed with the Tibet Air Shower (AS) experiment still awaits theoretical interpretation. The observed global feature of the anisotropy is well reproduced by a superposition of the bi-directional and uni-directional flows (BDF and UDF, respectively) of GCRs. If the orientation of the deduced BDF represents the orientation of the local interstellar magnetic field (LISMF), as indicated by best-fitting a model to the data, the UDF deviating from the BDF orientation implies a significant contribution from the streaming perpendicular to the LISMF. This perpendicular streaming is probably due to the drift anisotropy, because the contribution from the perpendicular diffusion is expected to be much smaller than the drift effect. The large amplitude deduced for the UDF indicates a large spatial gradient of the GCR density. We suggest that such a density gradient can be expected at the heliosphere sitting close to the boundary of the Local Interstellar Cloud (LIC), if the LIC is expanding. The spatial distribution of GCR density in the LIC reaches a stationary state because of the balance between the inward cross-field diffusion and the adiabatic cooling due to the expansion. We derive the steady-state distribution of GCR density in the LIC based on radial transport of GCRs in a spherical LIC expanding at a constant rate. By comparing the expected gradient with the observation by Tibet experiment, we estimate the perpendicular diffusion coefficient of multi-TeV GCRs in the local interstellar space. △ Less

Submitted 5 September, 2009; originally announced September 2009.

Comments: 4 pages, 2 figures, 31st International Cosmic Ray Conference (Lodz, Poland), 2009

Large-scale sidereal anisotropy of multi-TeV galactic cosmic rays and the heliosphere

Authors: M. Amenomori, X. J. Bi, D. Chen, S. W. Cui, Danzengluobu, L. K. Ding, X. H. Ding, C. Fan, C. F. Feng, Zhaoyang Feng, Z. Y. Feng, X. Y. Gao, Q. X. Geng, Q. B. Gou, H. W. Guo, H. H. He, M. He, K. Hibino, N. Hotta, Haibing Hu, H. B. Hu, J. Huang, Q. Huang, H. Y. Jia, L. Jiang , et al. (61 additional authors not shown)

Abstract: We develop a model anisotropy best-fitting to the two-dimensional sky-map of multi-TeV galactic cosmic ray (GCR) intensity observed with the Tibet III air shower (AS) array. By incorporating a pair of intensity excesses in the hydrogen deflection plane (HDP) suggested by Gurnett et al., together with the uni-directional and bi-directional flows for reproducing the observed global feature, this m… ▽ More We develop a model anisotropy best-fitting to the two-dimensional sky-map of multi-TeV galactic cosmic ray (GCR) intensity observed with the Tibet III air shower (AS) array. By incorporating a pair of intensity excesses in the hydrogen deflection plane (HDP) suggested by Gurnett et al., together with the uni-directional and bi-directional flows for reproducing the observed global feature, this model successfully reproduces the observed sky-map including the "skewed" feature of the excess intensity from the heliotail direction, whose physical origin has long remained unknown. These additional excesses are modeled by a pair of the northern and southern Gaussian distributions, each placed ~50 degree away from the heliotail direction. The amplitude of the southern excess is as large as ~0.2 %, more than twice the amplitude of the northern excess. This implies that the Tibet AS experiment discovered for the first time a clear evidence of the significant modulation of GCR intensity in the heliotail and the asymmetric heliosphere. △ Less

Submitted 5 September, 2009; originally announced September 2009.

Comments: 4 pages, 1 figure, 31st International Cosmic Ray Conference (Lodz, Poland), 2009

A Re-interpretation of the STEREO/STE Observations and it's Consequences

Authors: K. C. Hsieh, P. C. Frisch, J. Giacalone, J. R. Jokipii, J. Kota, D. E. Larson, R. P Lin, J. G. Luhmann, L. Wang

Abstract: We present an alternate interpretation of recent STEREO/STE observations that were originally attributed to energetic neutral atoms (ENA) from the heliosheath. The signal attributed to the diffuse ENA source instead shows the characteristics of a point source. We point out that the peak intensity seen by STEREO/STE is centered at the ecliptic longitude of the bright X-ray source Sco X-1. The obs… ▽ More We present an alternate interpretation of recent STEREO/STE observations that were originally attributed to energetic neutral atoms (ENA) from the heliosheath. The signal attributed to the diffuse ENA source instead shows the characteristics of a point source. We point out that the peak intensity seen by STEREO/STE is centered at the ecliptic longitude of the bright X-ray source Sco X-1. The observed energy spectrum and intensity are also consistent with the X-rays from Sco X-1. The problem of energy dissipation at the solar wind termination shock remains unsolved while current understanding of the interaction between the solar wind and interstellar wind awaits future observations. △ Less

Submitted 20 February, 2009; originally announced February 2009.

Comments: Accepted by ApJL

Sky-maps of the sidereal anisotropy of galactic cosmic ray intensity and its energy dependence

Authors: K. Munakata, N. Matsumoto, S. Yasue, C. Kato, S. Mori, M. Takita, M. L. Duldig, J. E. Humble, J. Kota

Abstract: We analyze the sidereal daily variations observed between 1985 and 2006 at Matsushiro, Japan (MAT) and between 1993 and 2005 at Liapootah, Tasmania (LPT). These stations comprise the two hemisphere network (THN) of underground muon detectors in Japan and Australia. Yearly mean harmonic vectors at MAT and LPT are more or less stable without any significant change in phase and amplitude in correla… ▽ More We analyze the sidereal daily variations observed between 1985 and 2006 at Matsushiro, Japan (MAT) and between 1993 and 2005 at Liapootah, Tasmania (LPT). These stations comprise the two hemisphere network (THN) of underground muon detectors in Japan and Australia. Yearly mean harmonic vectors at MAT and LPT are more or less stable without any significant change in phase and amplitude in correlation with either the solar activity- or magnetic-cycles. In this paper, therefore, we analyze the average anisotropy over the entire observation periods, i.e. 1985-2006 for MAT and 1993-2005 for LPT. We apply to the THN data a best-fitting analysis based on a model anisotropy in space identical to that adopted by Amenomori et al. (2007) for Tibet III data. The median energies of primary cosmic rays recorded are ~0.5 TeV for THN and ~5 TeV for the Tibet III experiment. It is shown that the intensity distribution of the best-fit anisotropy is quite similar to that derived from Tibet III data, regardless of the order of magnitude difference in energies of primary particles. This, together with the THN observations, confirms that the analysis by Amenomori et al. (2007) based on the Tibet III experiment in the northern hemisphere is not seriously biased. The best-fit amplitudes of the anisotropy, on the other hand, are only one third or less of those reported by the Tibet III experiment, indicating attenuation due to solar modulation. The rigidity dependence of the anisotropy amplitude in the sub-TeV region is consistent with the spectrum reported by Hall et al. (1999), smoothly extending to the Tibet III result in the multi-TeV region. The amplitude at higher energies appears almost constant or gradually decreasing with increasing rigidity. △ Less

Submitted 4 November, 2008; v1 submitted 4 November, 2008; originally announced November 2008.

Comments: 13 pages, 1 table, 4 figures, submitted to the Advances in Geoschiences

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.

Comments: 35 pages, 10 figures, submitted to the Astrophysical Journal

Journal ref: Astrophys. J. 681:693-707 2008