Simulation of plasmas in the electromagnetic fields requires to solve numerically a kinetic equation, describing the time evolution of the particle distribution function. Here, we propose a finite volume scheme based on the integral relation for the Poisson bracket to solve the most fundamental kinetic equation, namely, the Liouville equation. The proposed scheme conserves the number of particles, maintains the total-variation-diminishing (TVD) property, and provides high-quality numerical results. Some other types of kinetic equations may be also formulated in terms of the Poisson brackets and solved with the proposed method. Among them is the focused transport equation describing the acceleration and propagation of the Solar Energetic Particles (SEPs), which is of practical importance, since the high energy SEPs produce radiation hazards. The newly proposed scheme is demonstrated to be accurate and efficient, which makes it applicable to global simulation systems analysing the space weather. We also discuss a role of focused transport and the accuracy of the diffusive approximation, in application to the SEPs

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 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.

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.

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 $\kappa_\perp$ is generally much smaller than parallel diffusion coefficient $\kappa_\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.

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.

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.

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.

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.

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.

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.