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.

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.

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.

The paper is a script of a lecture given at the ISAPP-Baikal summer school in 2018. The lecture gives an overview of the Tunka Advanced Instrument for cosmic rays and Gamma Astronomy (TAIGA) facility including historical introduction, description of existing and future setups, and outreach and open data activities.

Lateral distributions of electrons in air showers initiated by photons of ultra high energies ($10^{18}-10^{22}$ eV) obtained on the basis of numerical solution of adjoint cascade equations are presented. An extended analysis has been made considering separately the Landau-Pomeranchuk-Migdal (LPM) effect and the interaction of photons and electrons with the geomagnetic field (GMF) with respect to the scaling formalism for lateral distributions. It is shown that one-parametric scaling description of the lateral distribution of electrons remains valid up to the highest energies considering the LPM and GMF effects, that allow effective primary particle type discrimination using the surface detectors data of largest ground-based air shower arrays.

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.

Scaling properties of nuclei- and photon-initiated air showers are examined in wide primary energy range ($10^{14}\div 10^{22}$~eV) taking into account Landau-Pomeranchuk-Migdal and geomagnetic field effects. It is shown that the invariance in lateral distribution and also the universal dependence between shower age and root mean square radius of electron component exist up to the highest energies. The implications of shower universality for reanalysis, cross-checking and validation of the results of different experiments together with decreasing of the influence of hadronic model uncertainties are discussed in detail.

It is shown, that primary proton spectrum, reconstructed from sea-level and underground data on muon spectrum with the use of QGSJET 01, QGSJET II, NEXUS 3.97 and SIBYLL 2.1 interaction models, demonstrates not only model-dependent intensity, but also model-dependent form. For correct reproduction of muon spectrum shape primary proton flux should have non-constant power index for all considered models, except SIBYLL 2.1, with break at energies around 10-15 TeV and value of exponent before break close to that obtained in ATIC-2 experiment. To validate presence of this break understanding of inclusive spectra behavior in fragmentation region in p-air collisions should be improved, but we show, that it is impossible to do on the basis of the existing experimental data on primary nuclei, atmospheric muon and hadron fluxes.

Influence of shower fluctuations on the shape of lateral distribution of electrons in EAS of fixed size measured by scintillation counters is analyzed in framework of scaling formalism. Correction factors for the mean square radius of electrons are calculated for the experimental conditions of KASCADE array. Possible improvement of the primary mass discrimination by analysis of lateral distribution of EAS electrons is discussed in detail.

Self-consistency of interaction models QGSJET 01, SIBYLL 2.1, NEXUS 3.97 and QGSJET II is checked in terms of their ability to reproduce simultaneously experimental data on fluxes of muons and hadrons. From this point of view SIBYLL 2.1 gives the most acceptable, though not quite satisfactory, results. Analysis of the situation for muons supports our previous conclusions, that high-energy muon deficit is due both to underestimation of primary light nuclei fluxes in direct emulsion chamber experiments and to softness of $p+A\to\pi^\pm,K^\pm+X$ inclusive spectra in fragmentation region, especially prominent in case of QGSJET 01 model.

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.

The problem of the origin of cosmic rays in the shocks produced by supernova explosions at energies below the so called 'knee' (at ~3*10$^6$ GeV) in the energy spectrum is addressed, with special attention to the propagation of the particles through the inhomogenious interstellar medium and the need to explain recent anisotropy results, [1]. It is shown that the fractal character of the matter density and magnetic field distribution leads to the likelihood of a substantial increase of spatial fluctuations in the cosmic ray energy spectra. While the spatial distribution of cosmic rays in the vicinity of their sources (eg. inside the Galactic disk) does not depend much on the character of propagation and is largely determined by the distribution of their sources, the distribution at large distances from the Galactic disk depends strongly on the character of the propagation. In particular, the fractal character of the ISM leads to what is known as 'anomalous diffusion' and such diffusion helps us to understand the formation of Cosmic Ray Halo. Anomalous diffusion allows an explanation of the recent important result from the Chacaltaya extensive air shower experiment [1], viz. a Galactic Plane Enhancement of cosmic ray intensity in the Outer Galaxy, which is otherwise absent for the case of the so-called 'normal' diffusion. All these effects are for just one reason: anomalous diffusion emphasizes the role of local phenomena in the formation of cosmic ray characteristics in our Galaxy and elsewhere.

We have made a new calculation of the cosmic ray electron spectrum using an anomalous diffusion model to describe the propagation of electrons in the Galaxy. The parameters defining the anomalous diffusion have been recently determined from the study of nuclei propagation. The predicted electron spectrum is in a good agreement with the measurements. The source spectral index, found from experimental data, in this approach turns out to be equal to 2.95.

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.

We consider the propagation of galactic cosmic rays under assumption that the interstellar medium is a fractal one. An anomalous diffusion equation in terms of fractional derivatives is used to describe of cosmic ray propagation. The anomaly in used model results from large free paths ("Levy flights") of particles between galactic inhomogeneities and long time a particle stays in a trap. An asymptotical solution of the anomalous diffusion equation for point instantaneous and impulse sources with inverse power spectrum relating to supernova bursts is found. It covers both subdiffusive and superdiffusive regimes and is expressed in terms of the stable distributions. The energy dependence of spectral exponent of observed particles in different regimes is discussed.

Theoretical predictions for lateral distribution function of electrons in extensive air showers based on scaling formalism are presented. Our results are tested by comparison with AGASA experimental data taking into account the contribution of low-energy muons and simultaneously the effect of scintillation detectors response, according to recent simulation results. The possibility for ultrahigh energy cosmic ray primary composition deduction from the shape of LDF is discussed in detail.

Detailed theoretical study of the mean square radius of extensive air shower electrons has been made in connection with further development of scaling formalism for electron lateral distribution function. A very simple approximation formula, which allows joint description of all our results obtained in wide primary energy range and for different observation depths is presented. The sensitivity of the mean square radius to variations of basic parameters of hadronic interaction model is discussed.