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

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 consider the problem of night-time atmosphere monitoring at the locations of the Yakutsk EAS array and the TAIGA observatory for the annual October-April period of observations at these facilities. It is shown that the use of data from instruments aboard Terra, Aqua, Suomi-NPP and NOAA-20 satellites, received by the ground stations of the Altai State University, Russia, and processed to Level 2 (retrieval of geophysical parameters of the atmosphere), gives the possibility of the near-real-time monitoring observations of the cloud cover structure and the temperature profiles at night with a spatial resolution and a frequency sufficient for meteorological correction of the detector readings.

New analysis of lateral distributions (LD) of electrons, muons and charged particles in extensive air showers generated by $10^{15}-10^{19}$ eV cosmic rays has been made with respect to the scaling formalism for LD and air shower universality paradigm. It is found that lateral distributions of electrons and muons, when considered separately, are well described by one-parametric scaling functions, while for the charged particles mixture distributions the scaling description is not accurate in radial distance range where electrons and muons give a comparable contribution to the local particle densities. The use of scaling formalism enables enhancing the electron/muon separation capabilities and, when combined with theoretically motivated air shower universality concept, provides robust mass composition estimates. The proposed approach could be implemented for the present and future (multi-) hybrid air shower observations, especially in realizing the potential of Auger and Telescope Array observatories upgrades, as well as for re-analysis and cross-calibration of the data collected from different air shower arrays in a broad primary energy range.

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.

The energy spectrum of high energy cosmic rays emitted by sources change during propagation in the Galaxy. Using the spectrum observed at the Earth and two different transport models, based on anomalous and normal diffusion equations, we have retrieved an average value of the cosmic ray injection exponent at Galactic sources taking into account the inhomogeneity of the interstellar medium. We have shown that the average value of the injection index $p$, obtained in the framework of both transport models, equals to $p\sim (2.8-3.0)$.

On the basis of the scaling approach and CORSIKA simulations data the radial scale factor of lateral distribution of electrons in extensive air showers is confirmed as potentially effective primary mass estimator, and its sensitivity to hadronic interaction model is investigated. It is shown that improved composition results both on average and event-by-event basis can be achieved taking into account the universality property of air shower development expressed by the relation between radial scale factor and longitudinal age parameter. The enhancements of such a theoretically motivated tool for the unbiased cosmic ray composition deduction in a wide primary energy range from (multi-)hybrid air shower measurements of nearest future are discussed.