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.

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.