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.

All information about primary cosmic rays above the knee has been obtained from results of EAS investigations. At that, two alternative approaches exist: cosmophysical and nuclear physical. In the frame of the first one, all changes in measured EAS characteristics are explained by the changes in energy spectrum and mass composition of primary cosmic rays. In this paper, the second approach is considered, in frame of which corresponding changes in EAS parameters are explained by changes of interaction model above the knee. Some experimental possibilities of proof of the correctness of the nuclear physical approach are considered.

Local muon density spectra (LMDS) at various zenith angles have been reconstructed from the data of two detectors of the Experimental complex NEVOD. The inclined muon bundles at the ground level were detected with the coordinate detector DECOR, and for the near-vertical direction with the calibration telescope system (CTS) of the Cherenkov water detector. In comparison with the earlier DECOR results, the experimental statistics has been increased by 2-3 times for different ranges of zenith angle and muon bundle multiplicity and is now based on about 40,000 h of the setup operation. The live time of measurements with CTS is about 12,000 h. It is found with both setups that the slope of LMDS is increasing above the primary energy of about 1017 eV. The details of the experiment and data analysis are presented.

The spatial distribution of Cherenkov radiation from cascade showers generated by muons in water has been measured with Cherenkov water calorimeter (CWC) NEVOD. This result allowed to improve the techniques of treating cascade showers with unknown axes by means of CWC response analysis. The techniques of selecting the events with high energy cascade showers and reconstructing their parameters are discussed. Preliminary results of measurements of the spectrum of cascade showers in the energy range 100 GeV - 20 TeV generated by cosmic ray muons at large zenith angles and their comparison with expectation are presented.

Results of investigations of the near-horizontal muons in the range of zenith angles of 85-95 degrees are presented. In this range, so-called "albedo" muons (atmospheric muons scattered in the ground into the upper hemisphere) are detected. Albedo muons are one of the main sources of the background in neutrino experiments. Experimental data of two series of measurements conducted at the experimental complex NEVOD-DECOR with the duration of about 30 thousand hours "live" time are analyzed. The results of measurements of the muon flux intensity are compared with simulation results using Monte-Carlo on the basis of two multiple Coulomb scattering models: model of point-like nuclei and model taking into account finite size of nuclei.

Experimental complex (EC) NEVOD includes a number of unique experimental facilities for studies of main components of cosmic rays on the Earth's surface. The complex is used for the basic research of CR flux characteristics and their interactions in the energy range 10^15 - 10^19 eV, and for applied investigations directed to the development of methods of the muon diagnostics of the atmosphere and the Earth's magnetosphere and near-terrestrial space. To extend the experimental capabilities and raising the status of the installation to the Mega Science level, nowadays new large-scale detectors: array for the EAS registration - NEVOD-EAS, detector of atmospheric neutrons - URAN, and large-area coordinate-tracking detector - TREK, are being deployed around EC NEVOD. The description of new detectors and a common trigger system to ensure the joint operation together with other detectors of EC NEVOD are presented.

A new cluster type shower array NEVOD-EAS is designed for estimating the size, axis position and arrival direction of extensive air showers registered by the Cherenkov water detector NEVOD and coordinate detector DECOR (Moscow, Russia). In 2015-2016, the central part of the array was deployed and started its operation. It includes 4 independent clusters of scintillation detector stations located around the NEVOD-DECOR experimental complex on the area of about $10^{4} m^{2}$. This article presents the results of studying amplitude and timing characteristics of the array clusters which are critical for EAS parameters reconstruction, as well as the examples of registered events.

Change of the energy characteristics of muon bundles with an increase of the primary cosmic ray particles energy can be a key to solving the problem of muon excess in the extensive air showers (EAS) observed in a number of experiments. In this work the data on the energy deposit of multi-muon events in a wide range of zenith angles (and as a consequence in a wide range of primary particles energies) obtained with NEVOD-DECOR setup over a long time period are presented. The experimental data are compared with the results of simulations of EAS muon component performed using CORSIKA code.

Differential and integral energy spectra of cosmic ray muons in the energy range from several TeV to ~ 1 PeV obtained by means of the analysis of multiple interactions of muons (pair meter technique) in the Baksan underground scintillation telescope (BUST) are presented. The results are compared with preceding BUST data on muon energy spectrum based on electromagnetic cascade shower measurements and depth-intensity curve analysis, with calculations for different muon spectrum models, and also with data of other experiments.

A new approach to investigations of ultra-high energy cosmic rays based on the ground-level measurements of the spectra of local density of EAS muons at various zenith angles is considered. Basic features of the local muon density phenomenology are illustrated using a simple semi-analytical model. It is shown that muon density spectra are sensitive to the spectrum slope, primary composition, and to the features of hadronic interaction. New experimental data on muon bundles at zenith angles from 30 degrees to horizon obtained with the coordinate detector DECOR are compared with CORSIKA-based simulations. It is found that measurements of muon density spectra in inclined EAS give possibility to study characteristics of primary cosmic ray flux in a very wide energy range from 10^15 to 10^19 eV.

Muon rate variations during Forbush effects registered by means of muon detectors DECOR, TEMP and URAGAN operated in the experimental complex NEVOD (MEPhI, Moscow) have been studied, and comparative analysis with neutron monitor data has been performed. The ratio of values of Forbush decreases in muon and in neutron fluxes is about one third, and preliminary rigidity dependence of Forbush decrease amplitude using muon data has been also obtained (for 2.4 GV cut-off rigidity). The detection of muon flux in the hodoscopic mode allows to study the dynamics of muon flux anisotropy related with magnetic field perturbations. Results of analysis of data from the new unique muon detector URAGAN indicate the change of muon flux asymmetry direction during the Forbush decrease. This phenomenon is related with the motion of solar plasma cloud and hodoscope acceptance cone relative to each other.