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.

We present a summary of recent tests and measurements of hadronic interaction properties with air showers. This report has a special focus on muon density measurements. Several experiments reported deviations between simulated and recorded muon densities in extensive air showers, while others reported no discrepancies. We combine data from eight leading air shower experiments to cover shower energies from PeV to tens of EeV. Data are combined using the z-scale, a unified reference scale based on simulated air showers. Energy-scales of experiments are cross-calibrated. Above 10 PeV, we find a muon deficit in simulated air showers for each of the six considered hadronic interaction models. The deficit is increasing with shower energy. For the models EPOS-LHC and QGSJet-II.04, the slope is found significant at 8 sigma.

The Sydney University Giant Air-shower Recorder (SUGAR) measured the energy spectrum of ultra-high-energy cosmic rays reconstructed from muon-detector readings, while the Pierre Auger Observatory, looking at the same Southern sky, uses the calorimetric fluorescence method for the same purpose. Comparison of their two spectra allows us to reconstruct the empirical dependence of the number of muons in the shower on the primary energy for energies between $10^{17}$ and $10^{18.5}$ eV. We compare this dependence with the predictions of hadronic interaction models \mbox{QGSJET-II-04} and \mbox{EPOS-LHC}. The empirically determined number of muons with energies above 0.75 GeV exceeds the simulated one by the factors $\sim$1.67 and $\sim$1.28 for $10^{17}$ eV proton and iron primaries, respectively. The muon excess grows moderately with the primary energy, increasing by an additional factor of $\sim 1.2$ for $10^{18.5}$ eV primaries.

The origin and nature of extreme energy cosmic rays (EECRs), which have energies above the 50 EeV, the Greisen-Zatsepin-Kuzmin (GZK) energy limit, is one of the most interesting and complicated problems in modern cosmic-ray physics. Existing ground-based detectors have helped to obtain remarkable results in studying cosmic rays before and after the GZK limit, but have also produced some contradictions in our understanding of cosmic ray mass composition. Moreover, each of these detectors covers only a part of the celestial sphere, which poses problems for studying the arrival directions of EECRs and identifying their sources. As a new generation of EECR space detectors, TUS (Tracking Ultraviolet Set-up), KLYPVE and JEM-EUSO, are intended to study the most energetic cosmic-ray particles, providing larger, uniform exposures of the entire celestial sphere. The TUS detector, launched on board the Lomonosov satellite on April 28, 2016, from Vostochny Cosmodrome in Russia, is the first of these. It employs a single-mirror optical system and a photomultiplier tube matrix as a photo-detector and will test the fluorescent method of measuring EECRs from space. Utilizing the Earth's atmosphere as a huge calorimeter, it is expected to detect EECRs with energies above 100 EeV. It will also be able to register slower atmospheric transient events: atmospheric fluorescence in electrical discharges of various types including precipitating electrons escaping the magnetosphere and from the radiation of meteors passing through the atmosphere. We describe the design of the TUS detector and present results of different ground-based tests and simulations.

TUS (Tracking Ultraviolet Set-up), the first orbital detector of extreme energy cosmic rays (EECRs), those with energies above 50 EeV, was launched into orbit on April 28, 2016, as a part of the Lomonosov satellite scientific payload. The main aim of the mission is to test a technique of registering fluorescent and Cherenkov radiation of extensive air showers generated by EECRs in the atmosphere with a space telescope. We present preliminary results of its operation in a mode dedicated to registering extensive air showers in the period from August 16, 2016, to November 4, 2016. No EECRs have been conclusively identified in the data yet, but the diversity of ultraviolet emission in the atmosphere was found to be unexpectedly rich. We discuss typical examples of data obtained with TUS and their possible origin. The data is important for obtaining more accurate estimates of the nocturnal ultraviolet glow of the atmosphere, necessary for successful development of more advanced orbital EECR detectors including those of the KLYPVE (K-EUSO) and JEM-EUSO missions.

TUS (Track Ultraviolet Setup) is the first space experiment aimed to check the possibility of registering extreme energy cosmic rays (EECRs) at E>50 EeV by measuring the fluorescence signal of extensive air showers in the atmosphere. The detector operates as a part of the scientific payload of the Lomonosov satellite for more than a year. We describe an algorithm of searching for EECR events in the TUS data and briefly discuss a number of candidates selected by formal criteria.

Some discrepancies have been reported between observed and simulated muon content of extensive air showers: the number of observed muons exceeded the expectations in HiRes-MIA, Yakutsk and Pierre Auger Observatory data. Here, we analyze the data of the Moscow State University Extensive Air Shower (EAS-MSU) array on E_mu>~10 GeV muons in showers caused by ~(10^17-10^18) eV primary particles and demonstrate that they agree with simulations (QGSJET-II-04 hadronic interaction model) once the primary composition inferred from the surface-detector data is assumed.

On April 28, 2016, the Lomonosov satellite, equipped with a number of scientific instruments, was launched into orbit. Here we present briefly some of the results obtained with the first orbital telescope of extreme energy cosmic rays TUS and by a group of detectors aimed at multi-messenger observations of gamma-ray bursts.

Results of the search for $\sim (10^{16} - 10^{17.5})$ eV primary cosmic-ray photons with the data of the Moscow State University (MSU) Extensive Air Shower (EAS) array are reported. The full-scale reanalysis of the data with modern simulations of the installation does not confirm previous indications of the excess of gamma-ray candidate events. Upper limits on the corresponding gamma-ray flux are presented. The limits are the most stringent published ones at energies $\sim 10^{17}$ eV.

The Moscow State University Extensive Air Shower (EAS-MSU) array studied high-energy cosmic rays with primary energies ~(1-500) PeV in the Northern hemisphere. The EAS-MSU data are being revisited following recently found indications to an excess of muonless showers, which may be interpreted as the first observation of cosmic gamma rays at ~100 PeV. In this paper, we present a complete Monte-Carlo model of the surface detector which results in a good agreement between data and simulations. The model allows us to study the performance of the detector and will be used to obtain physical results in further studies.

Results of a blind search for localized regions of an excessive flux of cosmic rays in the energy range from 50 TeV to 20 PeV with the data of the FIAN KLARA-Chronotron experiment, the EAS MSU array and the Prototype of the EAS-1000 array are presented. A number of regions with a significant excess of the registered flux over an expected isotropic background are found. Some of the regions are present in at least two of the data sets considered.

Archival EAS-MSU data are searched for anomalous muonless events which may be caused by primary gamma rays with energies between 10^15 eV and 10^18 eV. We consider a refined sample of high-quality data and confirm the previously reported detection of a non-zero gamma-ray flux at ~5x10^16 eV with a similar flux value but at somewhat lower statistical significance, corresponding to the depletion of the sample. We present upper limits on the flux below and above these energies, including the first constraints in the range (10^17-10^18) eV never studied by any other experiment.

We reanalyze archival EAS-MSU data in order to search for events with anomalously low content of muons with energies E_mu>10 GeV in extensive air showers with number of particles N_e>2*10^7. We confirm the first evidence for nonzero flux of primary cosmic gamma rays at energies E~10^17 eV. The estimated fraction of primary gamma rays in the flux of cosmic particles with energies E>5.4*10^16 eV is (0.43 +0.12-0.11)%, which corresponds to the intensity of (1.2 +0.4-0.3)*10^{-16} cm^{-2} s^{-1} sr^{-1}. The study of arrival directions does not favour any particular mechanism of the origin of the photon-like events.

The origin of ultra-high-energy cosmic rays (UHECRs) is one of the most intriguing problems of modern cosmic ray physics. We briefly review the main astrophysical models of their origin and the forthcoming orbital experiments TUS and JEM-EUSO, and discuss how the new data can help one solve the long-standing puzzle.

We present spectra and mass composition of cosmic rays incoming to the Earth in the energy range (0.5-2)10^20 eV. As their sources we consider Seyfert galaxies located at distances < 40 Mpc, following an acceleration model for such moderate-power objects. Mass composition of the particles at sources is assumed to be mixed. Generation spectra are described by a function E^(-g0), where g0 is an arbitrary parameter. It is shown that the assumptions adopted make it possible to describe experimental data provided by HiRes and Pierre Auger Observatory, using different values of g0.

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.

We present computing results of the Cherenkov light propagation in air and water from extensive air showers developing over the ocean. Limits on zenith angles of the showers, at which the registration of flashes of reflected Cherenkov photons by the satellite-based detector TUS is possible, are analyzed with consideration for waves on the ocean surface.

We have estimated the reflected component of Cherenkov radiation, which arises in developing of an extensive air shower with primary energy of 10^20 eV over the ocean surface. It has been shown that, under conditions of the TUS experiment, a flash of the reflected Cherenkov photons at the end of the fluorescence track can be identified in showers with zenith angles up to 20 degrees.

The new EAS Cherenkov array Tunka-133 with about 1 km**2 geometric acceptance area is installed in the Tunka Valley (50 km from Lake Baikal). The array willpermit a detailed study of cosmic ray energy spectrum and mass composition in the energy range of 10**15 - 10**18 eV with uniform method. The array consistsof 19 clusters, each composed of 7 optical detectors with 20 cm PMTs. Since November 2008, the first 12 clusters are in operation, commissioning of the whole array is planned for September 2009 (At the time of submission of this paperto electronic arXiv(February 2010) the completed Tunka-133 array is already taking data). We describe the array construction and DAQ, preliminary results and plans for the future development: deployment of radio-antennas and muon detectors network.

Angular distributions of extensive air showers with different number of charged particles in the range 2.5x10^5--4x10^7 are derived using the experimental data obtained with the EAS MSU array. Possible approximations of the obtained distributions with different empiric functions available in literature, are analysed. It is shown that the exponential function provides the best approximation of the angular distributions in the sense of the chi-squared criterion.

The final analysis of the Extensive Air Shower (EAS) maximum X_max depth distribution derived from the data of Tunka-25 atmospheric Cherenkov light array in the energy range 3.10^15 - 3.10^16 eV is presented. The perspectives of X_max studies with the new Cherenkov light array Tunka-133 of 1 km^2 area, extending the measurements up to 10^18 eV, are discussed.

The new EAS Cherenkov array TUNKA-133, with about 1 km**2 sensitive area, is being installed in the Tunka Valley. The investigated energy range is 10**15-10**18 eV. It will consist of 133 optical detectors based on EMI9350 PMTs. Optical detectors are grouped into 19 clusters with 7 detectors each. The detectors are connected to the cluster box with RG-58 cables. Every PMT signal is digitized in the cluster box with 200 MHz FADC. The cluster boxes are connected to the data acquisition center with a 1 Gb/s optical link. A detailed description of the data acquisition system (DAQ) is presented.

The new EAS Cherenkov array Tunka-133, with about 1 km**2 sensitive area, is being installed in the Tunka Valley since the end of 2005. This array will permit a detailed study of the cosmic ray energy spectrum and the mass composition in the energy range of 10**15-10**18 eV with a unique method. The array will consist of 19 clusters, each composed of 7 optical detectors. The first cluster started operation in October 2006. We describe the data acquisition system and present preliminary results from data taken with the first cluster.

The propagation of high-energy cosmic rays in the Galaxy is investigated. Solutions of a diffusion model are combined with numerically calculated trajectories of particles. The resulting escape path length and interaction path length are presented and energy spectra obtained at Earth are discussed. It is shown that the energy spectra for heavy elements should be flatter as compared to light ones due to nuclear interactions during the propagation process. The obtained propagation properties of ultra-heavy elements indicate that these elements could play an important role for the explanation of the second knee in the cosmic-ray energy spectrum around 400 PeV.

An efficient scheme for one-dimensional extensive air shower simulation and its implementation in the program CONEX are presented. Explicit Monte Carlo simulation of the high-energy part of hadronic and electromagnetic cascades in the atmosphere is combined with a numeric solution of cascade equations for smaller energy sub-showers to obtain accurate shower predictions. The developed scheme allows us to calculate not only observables related to the number of particles (shower size) but also ionization energy deposit profiles which are needed for the interpretation of data of experiments employing the fluorescence light technique. We discuss in detail the basic algorithms developed and illustrate the power of the method. It is shown that Monte Carlo, numerical, and hybrid air shower calculations give consistent results which agree very well with those obtained within the CORSIKA program.

The project of an EAS Cherenkov array in the Tunka valley/Siberia with an area of about 1 km2 is presented. The new array will have a ten times bigger area than the existing Tunka-25 array and will permit a detailed study of the cosmic ray energy spectrum and the mass composition in the energy range from 10^15 to 10^18 eV.

Two structures in the all-particle energy spectrum of cosmic rays, the knee at 4 PeV and the second knee around 400 PeV are proposed to be explained by a phenomenological model, the poly gonato model, connecting direct and indirect measurements. Within this approach the knee is caused by a successive cut-off of the flux for individual elements starting with protons at 4.5 PeV. The second knee is interpreted as the end of the stable nuclei of the periodic table. To check some key features of this model calculations of the cosmic ray energy spectrum and the propagation path length at energies from 10^14 to 10^19 eV have been performed within the framework of a combined approach based on the diffusion model of cosmic rays and a direct simulation of charged-particle trajectories in the Galaxy.

A hybrid simulation code is developed that is suited for fast one-dimensional simulations of shower profiles, including fluctuations. It combines the Monte Carlo simulation of high energy interactions with a fast numerical solution of cascade equations for the resulting distributions of secondary particles. Results obtained with this new code, called CONEX, are presented and compared to CORSIKA predictions.

The project of an EAS Cherenkov array in the Tunka valley/Siberia with an area of about 1 km^2 is presented. The new array will have a ten times bigger area than the existing Tunka-25 array and will permit a detailed study of the cosmic ray energy spectrum and the mass composition in the energy range from 10^15 to 10^18 eV.

We present an introductory lecture on general features of cosmic rays, for non-experts, and some recent developments concerning cascade equations for air shower developments.

We have analysed arrival times of extensive air showers (EAS) registered with the EAS-1000 prototype array during the period from August, 1997 till February, 1999. Our analysis has revealed that though the vast majority of samples of consecutive time intervals between EAS arrival times obey the exponential distribution, there are sequences of showers that have another distribution and thus violate the homogeneity hypothesis. The search for correlation between such events and clusters of showers and events with big delays between arrival times was also carried out.

The purpose of this paper is twofold: first we want to introduce a new type of hadronic interaction model (NEXUS), which has a much more solid theoretical basis as, for example, presently used models like QGSJET and VENUS, and ensures therefore a much more reliable extrapolation towards high energies. Secondly, we want to promote an extensive air shower (EAS) calculation scheme, based on cascade equations rather than explicit Monte Carlo simulations, which is very accurate in calculations of main EAS characteristics and extremely fast concerning computing time. We employ the NEXUS model to provide the necessary data on particle production in hadron-air collisions and present the average EAS characteristics for energies 10^14 - 10^17 eV. The experimental data of the casa-blanka group are analyzed in the framework of the new model.