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.

Turbulence is ubiquitous in the interstellar medium (ISM) of the Milky Way and other spiral galaxies. The energy source for this turbulence has been much debated with many possible origins proposed. The universality of turbulence, its reported large-scale driving, and that it occurs also in starless molecular clouds, challenges models invoking any stellar source. A more general process is needed to explain the observations. In this work we study the role of galactic spiral arms. This is accomplished by means of three-dimensional hydrodynamical simulations which follow the dynamical evolution of interstellar diffuse clouds (100cm-3) interacting with the gravitational potential field of the spiral pattern. We find that the tidal effects of the arm's potential on the cloud result in internal vorticity, fragmentation and hydrodynamical instabilities. The triggered turbulence result in large-scale driving, on sizes of the ISM inhomogeneities, i.e. as large as 100pc, and efficiencies in converting potential energy into turbulence in the range 10 to 25 percent per arm crossing. This efficiency is much higher than those found in previous models. The statistics of the turbulence in our simulations are strikingly similar to the observed power spectrum and Larson scaling relations of molecular clouds and the general ISM. The dependency found from different models indicate that the ISM turbulence is mainly related to local spiral arm properties, such as its mass density and width. This correlation seems in agreement with recent high angular resolution observations of spiral galaxies, e.g. M51 and M33.

We propose a new, more realistic, description of the perturbed gravitational potential of spiral galaxies, with spiral arms having Gaussian-shaped groove profiles. We investigate the stable stellar orbits in galactic disks, using the new perturbed potential. The influence of the bulge mass on the stellar orbits in the inner regions of a disk is also investigated. The new description offers the advantage of easy control of the parameters of the Gaussian profile of its potential. We find a range of values for the perturbation amplitude from 400 to 800 km^2 s^{-2} kpc^{-1} which implies a maximum ratio of the tangential force to the axisymmetric force between 3% and 6%, approximately. Good self-consistency of arm shapes is obtained between the Inner Lindblad resonance (ILR) and the 4:1 resonance. Near the 4:1 resonance the response density starts to deviate from the imposed logarithmic spiral form. This creates bifurcations that appear as short arms. Therefore the deviation from a perfect logarithmic spiral in galaxies can be understood as a natural effect of the 4:1 resonance. Beyond the 4:1 resonance we find closed orbits which have similarities with the arms observed in our Galaxy. In regions near the center, in the presence of a massive bulge, elongated stellar orbits appear naturally, without imposing any bar-shaped potential, but only extending the spiral perturbation a little inward of the ILR. This suggests that a bar is formed with a half-size around 3 kpc by a mechanism similar to that of the spiral arms. The potential energy perturbation that we adopted represents an important step in the direction of self-consistency, compared to previous sine function descriptions of the potential. Our model produces a realistic description of the spiral structure, able to explain several details that were not yet understood.