Recent balloon-borne and satellite experiments have established new features in the behavior of the spectra of cosmic rays. An analysis of all the data showed that hardening of most abundant primary cosmic ray nuclei spectra with increasing rigidity is observed in $>200$ GV region. At the same time, rigidity dependencies of secondary cosmic rays are distinctly different. The AMS-02 data show that above 200 GV the secondary cosmic rays Li, Be, B harden more than the primary He, C, O. In this paper we discuss a new scenario that self-consistently describes these new features of the cosmic ray spectra. We demonstrate that the measured by the AMS-02 changes in the slope of the cosmic ray spectra caused by the transition from the contribution of multiple distant Galactic sources, including the nuclear interactions of the particles accelerated by these sources with an interstellar medium during their wandering in the Galaxy, to the contribution of mainly local ones. We also found that the spectral universality is observed at rigidity $R> 10^5$ GV.

We consider the problem of the cosmic ray spectrum formation assuming that cosmic rays are produced by Galactic sources. The anomalous diffusion equation proposed in our recent papers is used to describe cosmic ray propagation in the interstellar medium. We show that in the framework of this approach and with generation spectrum exponent $\gamma=2.85$ it is possible to reproduce locally observed basic features of cosmic rays in the energy region of $10^{10}\div 10^{20}$ eV: difference between spectral exponents of protons and other nuclei, mass composition variation, `knee' problem, flattening of the primary spectrum at $E \geq 10^{18}\div 10^{19}$ eV. The crucial model predictions for the mass composition behaviour in the ultra-high energy region are discussed.

We present the results of new calculations of the energy spectra of cosmic ray electrons, positrons and also positron fraction under assumption that both electrons and positrons are generated by the same Galactic sources, which accelerate particles with same power-law injection spectral index $p$. The value of the injection index $p\approx 2.85$ retrieved in our works from the analysis of the observed CR spectra has been used. The propagation of particles through the highly non-homogeneous interstellar medium has been described by the anomalous diffusion model. We show that proposed approach allows the self-consistent description of the recent AMS-02 data. In contrast to the conclusion made by AMS-02 collaboration about the different origin of the high-energy positrons and electrons, our results demonstrate that the differing behavior of spectral indices of electrons and positrons with energy can be described under the assumption that both positrons and electrons have common Galactic sources. Our prediction (Lagutin et al. Bull. Russian Academy of Sciences. Physics, 2013, Vol. 77, No. 11, p. 1312) that the $e^{+}$ to total $e^{-} + e^{+}$ ratio reaches a constant value for E $\gg 200$ GeV is confirmed.

The possibility for a self-consistent description of all the basic features of the observed cosmic ray spectra and primary composition variations in the energy range of $10^{15}\div 10^{20}$ eV within the Galactic origin scenario is examined. We assume the existence of Galactic sources that accelerate particles up to $\sim 3\cdot 10^{18}Z$ eV and take into account a highly inhomogeneous (fractal-like) distribution of matter and magnetic fields in the Galaxy that leads to extremely large free paths of particles ("L\'{e}vy flights"), along with an overwhelming contribution to the cosmic ray fluxes observed above $\sim 10^{18}$ eV from particles reaching the Solar System without scattering. Our scenario was refined on the basis of recent experimental results on primary mass composition. Model predictions, which could be verified with the improved high-precision measurements in the nearest future are discussed.

Analysis of the experimental data on cosmic ray spectra in the framework of the proposed model with two types of sources leads to conclusion, that sources with particle generation spectral exponent $p\sim 2.85$ give the major contribution to the all-particle spectrum in the energy range $10^5-10^7$ GeV. `Fine structure' of spectrum around the `knee' may arise due to presence of nearby supernova type source, accelerating particles up to the energies $\sim3\cdot 10^4 Z$ GeV, if the energy output of such source is $\sim2\cdot10^{48}$ erg/source.

Primary proton spectrum in the energy range $5-10^3$ TeV is reconstructed from the sea level muon spectrum with the use of QGSJET01 and SYBILL2.1 interaction models. Heavier nuclei are taken in accordance with the direct measurements data, 100% uncertainty in helium flux is accounted for. The obtained proton intensity strongly contradicts to the available data of balloon experiments, exceeding them at the least by 100% for QGSJET01. This discrepancy is due to the combined effect of primary nucleon flux underestimation in the direct measurements and incorrect description of extensive air shower development. In the latter case it is required earlier shower development and harder spectra of secondary pions and kaons in comparison with QGSJET01. This conclusion is in agreement with the obtained by the KASCADE group on the basis of events rate study.

In this paper we discuss the problem, why the use of the direct data on primary nuclei spectra together with the modern hadronic interaction models leads to significant deficit of computed vertical muon flux at sea level for energies $>100$ GeV. We suggest, that to find out the source of this inconsistency it is necessary to perform an analysis of sensitivity of emulsion chamber data to variations of hadron-nucleus interaction characteristics. Such analysis will give more ground for discussion of adequacy of the up-to-date interaction models and of mutual compatibility of primary nuclei spectra, obtained in direct and EAS experiments.

We discuss the problem of cosmic ray mass composition variation in wide energy region 10^{14}-10^{20} eV. The mass composition predicted in the framework of recently developed anomalous diffusion model is tested using the results of CORSIKA calculations and experimental data on the depth of maximum of extensive air showers. We show that the model predictions for the mass composition are consistent with observations of different experiments.

We consider the propagation of galactic cosmic rays in the fractal interstellar medium. Steady state solution of the fractional diffusion equation, describing cosmic ray propagation, is found. We show that the exponent of the steady state spectrum turns out to be equal to the exponent of the cosmic ray spectrum above the "knee", that is approximately equal to 3.1.

We calculate the energy spectra and mass composition of cosmic rays in energy region 1 - 10^8 GeV/particle under the assumption that cosmic rays propagation in the Galaxy is described by anomalous diffusion equation. Our results and comparisons with experimental data are presented.