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

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.

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.

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.

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.