We present a description of the new complex installation for the study of extensive air showers which was created at the Tien Shan mountain cosmic ray station, as well as the results of the test measurements made there in 2014-2016. At present, the system for registration of electromagnetic shower component consists of $\sim$100 detector points built on the basis of plastic scintillator plates with the sensitive area of 0.25m$^2$ and 1m$^2$, spread equidistantly over $\sim$10$^4$m$^2$ space. The dynamic range of scintillation amplitude measurements is currently about $(3-7)\cdot 10^4$, and there is a prospect of it being extended up to $\sim$10$^6$. The direction of shower arrival is defined by signal delays from a number of the scintillators placed cross-wise at the periphery of the detector system. For the investigation of nuclear active shower components there was created a multi-tier 55m$^2$ ionization-neutron calorimeter with a sum absorber thickness of $\sim$1000g/cm$^2$, typical spatial resolution of the order of 10cm, and dynamic range of ionization measurement channel about $\sim$10$^5$. Also, the use of saturation-free neutron detectors is anticipated for registration of the high- and low-energy hadron components in the region of shower core. A complex of underground detectors is designed for the study of muonic and penetrative nuclear-active components of the shower. The full stack of data acquisition, detector calibration, and shower parameters restoration procedures are now completed, and the newly obtained shower size spectrum and lateral distribution of shower particles occur in agreement with conventional data. Future studies in the field of $10^{14}-10^{17}$eV cosmic ray physics to be held at the new shower installation are discussed.

In the most works which deal with EAS the CR energy spectrum is deduced by means of the model defined dependence $E_0=a\cdot{N_e}^{\alpha}$. An electron total number $N_e$ is evaluated by integration of the NKG-function f(r). This algorithm breaks down for young EAS with age parameter $s\sim0$. This work shows, that part of the young EAS becomes large in the range $N_e\geq10^7$, it cause to divergency of the $N_e$ integral for them and distorts the shape of EAS (CR) spectrum. A final analysis of the experimental data permits to conclude that EAS spectrum has local maximum at $N_e\sim10^9$, which results in a decrease of the EAS spectrum slope for $N_e\geq10^7$ (inverse break or "knee"). A local maximum can arise because of the additional CR component in the range $E_0\geq10$ PeV.

We present a study of the high energy spectra of hadrons in cores of extensive air showers. These data were obtained for the first time in the hybrid {\it HADRON} experiment (Tien-Shan) by means of a large X-ray emulsion chamber combined with the shower array. In the local energy interval 3--100 PeV an increase in the energy of hadrons was found, which means the appearance of a penetrating component. This component in our experiment was observed in the atmosphere that indicates the presence of a penetrating strongly interacting component in primary cosmic rays. Along with that, it is worth emphasising that the region where this component is observed coincides with the region of the so-called knee in the spectrum of cosmic rays. On this basis, a new hypothesis of knee formation can be put forward.

The experimental gamma-families in which events called "halos" were detected using X-ray emulsion chambers (XRECs) of the PAMIR experiment are described. The halo nature is explained within the standard model of nuclear interactions. Based on events with halo, the fraction of p+He in primary cosmic rays (PCR) was estimated. It is shown that the PCR mass composition remains mixed in the range E0 = 1-100 PeV.

First preliminary results of the balloon-borne experiment SPHERE-2 on the all-nuclei primary cosmic rays (PCR) spectrum and primary composition are presented. The primary spectrum in the energy range $10^{16}$--$5\cdot10^{17}$ eV was reconstructed using characteristics of Vavilov-Cherenkov radiation of extensive air showers (EAS), reflected from a snow surface. Several sources of systematic uncertainties of the spectrum were analysed. A method for separation of the primary nuclei' groups based on the lateral distribution function' (LDF) steepness parameter is presented. Preliminary estimate of the mean light nuclei' fraction $f_{30-150}$ at energies $3\cdot10^{16}$--$1.5\cdot10^{17}$ eV was performed and yielded $f_{30-150}$= (21$\pm$11)%.

Here is presented the current state of the SPHERE-2 balloon-borne experiment. The detector is elevated up to 1 km above the snow surface and registers the reflected Vavilov-Cherenkov radiation from extensive air showers. This method has good sensitivity to the mass-composition of the primary cosmic rays due to its high resolution near the shower axis. The detector consists of a 1500 mm spherical mirror with a 109 PMT cluster in its focus. The electronics record a signal pulse profile in each PMT. In the last 2 years the detector was upgraded: time resolution of pulse registration was enhanced up to 12.5 ns, channel sensitivity was increased by a factor of 3, a new LED-based relative PMT calibration method was introduced, and new hardware and etc. was installed.