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.

The NUCLEON satellite experiment is designed to directly investigate the energy spectra of cosmic-ray nuclei and the chemical composition (Z=1-30) in the energy range of 2-500 TeV. The experimental results are presented, including the energy spectra of different abundant nuclei measured using the new Kinematic Lightweight Energy Meter (KLEM) technique. The primary energy is reconstructed by registration of spatial density of the secondary particles. The particles are generated by the first hadronic inelastic interaction in a carbon target. Then additional particles are produced in a thin tungsten converter, by electromagnetic and hadronic interactions.

The NUCLEON experiment was designed to study the chemical composition and energy spectra of galactic cosmic ray nuclei from protons to zinc at energies of $\sim10^{11}$--$10^{15}$\,eV per particle. The research was carried out with the NUCLEON scientific equipment installed on the Russian satellite "Resurs-P" No.\,2 as an additional payload. This article presents the results for the measured nuclei spectra related to the first approximately 250 days of the scientific data collection during 2015 and 2016. The all-particle spectrum and the spectra of p, He, C, O, Ne, Mg, Si and Fe are presented. Some interesting ratios of the spectra are also presented and discussed. The experiment is now in its beginning stage and the data still have a preliminary character, but they already give numerous indications of the existence of various non-canonical phenomena in the physics of cosmic rays, which are expressed in the violation of a simple universal power law of the energy spectra. These features of the the data are briefly discussed.

In this paper we propose a 'knee-like' approximation of the lateral distribution of the Cherenkov light from extensive air showers in the energy range 30-3000 TeV and study a possibility of its practical application in high energy ground-based gamma-ray astronomy experiments (in particular, in TAIGA-HiSCORE). The approximation has a very good accuracy for individual showers and can be easily simplified for practical application in the HiSCORE wide angle timing array in the condition of a limited number of triggered stations.

This work is a methodical study of another option of the hybrid method originally aimed at gamma/hadron separation in the TAIGA experiment. In the present paper this technique was performed to distinguish between different mass groups of cosmic rays in the energy range 200 TeV - 500 TeV. The study was based on simulation data of TAIGA prototype and included analysis of geometrical form of images produced by different nuclei in the IACT simulation as well as shower core parameters reconstructed using timing array simulation. We show that the hybrid method can be sufficiently effective to precisely distinguish between mass groups of cosmic rays.

This work is a methodical study on hybrid reconstruction techniques for hybrid imaging/timing Cherenkov observations. This type of hybrid array is to be realized at the gamma-observatory TAIGA intended for very high energy gamma-ray astronomy (>30 TeV). It aims at combining the cost-effective timing-array technique with imaging telescopes. Hybrid operation of both of these techniques can lead to a relatively cheap way of development of a large area array. The joint approach of gamma event selection was investigated on both types of simulated data: the image parameters from the telescopes, and the shower parameters reconstructed from the timing array. The optimal set of imaging parameters and shower parameters to be combined is revealed. The cosmic ray background suppression factor depending on distance and energy is calculated. The optimal selection technique leads to cosmic ray background suppression of about 2 orders of magnitude on distances up to 450 m for energies greater than 50 TeV.

A 'knee-like' approximation of Cherenkov light Lateral Distribution Functions, which we developed earlier, now is used for the actual tasks of background rejection methods for high energy (tens and hundreds of TeV) gamma-ray astronomy. In this work we implement this technique to the HiSCORE wide angle timing array consisting of Cherenkov light detectors with spacing of 100 m covering 0.2 km$^2$ presently and up to 5 km$^2$ in future. However, it can be applied to other similar arrays. We also show that the application of a multivariable approach (where 3 parameters of the knee-like approximation are used) allows us to reach a high level of background rejection, but it strongly depends on the number of hit detectors.

A simple 'knee-like' approximation of the Lateral Distribution Function (LDF) of Cherenkov light emitted by EAS (extensive air showers) in the atmosphere is proposed for solving various tasks of data analysis in HiSCORE and other wide angle ground-based experiments designed to detect gamma rays and cosmic rays with the energy above tens of TeV. Simulation-based parametric analysis of individual LDF curves revealed that on the radial distance 20-500 m the 5-parameter 'knee-like' approximation fits individual LDFs as well as the mean LDF with a very good accuracy. In this paper we demonstrate the efficiency and flexibility of the 'knee-like' LDF approximation for various primary particles and shower parameters and the advantages of its application to suppressing proton background and selecting primary gamma rays.

The role of nearby galactic sources, the supernova remnants, in formation of observed energy spectrum and anisotropy of high-energy cosmic rays is studied. The list of these sources is made up based on radio, X-ray and gamma-ray catalogues. The distant sources are treated statistically as ensemble of sources with random positions and ages. The source spectra are defined based on the modern theory of cosmic ray acceleration in supernova remnants while the propagation of cosmic rays in the interstellar medium is described in the frameworks of galactic diffusion model. Calculations of anisotropy are made to reproduce the experimental procedure of getting the two-dimensional anisotropy. The detailed comparison of calculations with cosmic ray data is made. We explained simultaneously the new cosmic ray data on a fine structure of all particle spectra around the knee and the amplitude and direction of the dipole component of anisotropy in the wide energy range 1 TeV - 1 EeV. Made assumptions do not look exotic, and they confirm the modern understanding of cosmic ray origin.

In this paper, a probable interpretation of a remarkable fine structure of all particle spectra between the knee and the ankle, as well as a high content of heavy nuclei around 10^17 eV measured recently in Tunka-133 and KASCADE Grande experiments, is presented in the model where production of cosmic rays in Type Ia SNRs provides observed cosmic rays (CR) flux around the knee. Subtracting the contribution of these sources from all particle spectrum we obtained the residual flux of CR in the transition region between Galactic and Extragalactic CR in the range 10^17-5 10^18 eV. The obtained spectrum also has a pronounced knee at energy (2-3)10^17 eV and slopes g1 ~ 1.8+-0.3 below it and g2 ~3.4 above it. We analyzed the possible contribution from known Galactic sources and found that formally the best candidate to contribute significantly to the transition region is Cas A. This source posses a number of unusual properties and considered as Type IIb SNR. However we showed that the hypothesis of extragalactic origin of CR starting from 10^17 eV seems more realistic.

A new EAS Cherenkov light array, Tunka-133, with ~1 km^2 geometrical area has been installed at the Tunka Valley (50 km from Lake Baikal) in 2009. The array permits a detailed study of cosmic ray energy spectrum and mass composition in the energy range 10^16 - 10^18 eV with a uniform method. We describe the array construction, DAQ and methods of the array calibration.The method of energy reconstruction and absolute calibration of measurements are discussed. The analysis of spatial and time structure of EAS Cherenkov light allows to estimate the depth of the EAS maximum X_max. The results on the all particles energy spectrum and the mean depth of the EAS maximum X_max vs. primary energy derived from the data of two winter seasons (2009 -- 2011), are presented. Preliminary results of joint operation of the Cherenkov array with antennas for detection of EAS radio signals are shown. Plans for future upgrades -- deployment of remote clusters, radioantennas and a scintillator detector network and a prototype of the HiSCORE gamma-telescope -- are discussed.

The simple leaky-box model of propagation of cosmic rays in the Galaxy is quite suitable for handling data on cosmic ray nuclei energy spectra and composition at E ? 1 GeV [1,2]. In the leakybox model a full information about cosmic ray propagation in Galaxy is compressed to the single parameter - escape length, Xe, characterizing mean grams of a matter passed by cosmic rays from sources to the Earth. In this paper we analyze the world data on proton and iron cosmic ray spectra collected in the past (HEAO, CRN et al.) and in series of recent electronic experiments (ATIC, CREAM, MS, BESS, Tracer et al.) and obtain the rigidity dependence of escape length, Xe(R) = R^(-0.47\pm-0.03), from the measured rigidity dependence of the protons/iron ratio. It quite agrees with the one stimated in standard manner from the secondary/primary nuclei ratio. But at R > 300 GV the behavior of Xe(R)distinctly changes, that can (variant of explanation)point out to the change of proton/iron ratio in cosmic ray sources.

We confirm the result obtained many years ago at Tien-Shan mountain station with the large 36-m2 lead calorimeter that in extensive air showers (EAS) with energies of few PeV the attenuation of core energy deposit in lead becomes slower than it could be predicted by modern codes. It is shown that this effect is connected with the appearance of the excess of abnormal cores with a large ionisation released in lower layers of the calorimeter and these abnormal EAS cores are most probably produced by high-energy muon groups. A few hypotheses of the excess of muon reach EAS cores are considered. To study the absorption of EAS hadrons and muons in a lead ionization calorimeter, the EAS development in the atmosphere was simulated in the framework of the CORSIKA+QGSJET code whereupon the passage of hadrons and muons through a lead calorimeter has been modeled with using the FLUKA transport code.

The present paper deals with calculations based on a standard model of cosmic ray acceleration in Supernova remnants (SNR) and on the latest astronomical data concerning the variety in Supernovae types, energies as well as the sites of their explosions. The results obtained allow explaination of the main features of all particle cosmic ray spectrum, including the 'knee'. This conclusion is based essentially on the premise that energy of accelerated CR in the given SNR is proportional to the kinetic energy of this explosion, that stresses the input of high energy explosions and suppresses the contribution of low energy explosions to the total CR flux.