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 energy range above the PCR "knee" (5-1000 PeV), the emulsion chamber data exhibit some new features which are hard to explain in the framework of the Standard Model. In this survey the results of a systematic quest of some new and exotic types of events which are observed in the Pamir emulsion chamber experiment are presented, namely, huge "halo" events, coplanar events, penetrating hadrons with abnormally weak absorption in lead, "Centauro"-type events. Possible theoretical approaches for their explanation are discussed.

The fraction of halo events among gamma-hadron families with SUM(Egamma) > 500 TeV is analyzed. The comparison of experimental data and Monte-Carlo simulations by MC0 code of quark-gluon string model under various assumptions on the PCR mass composition suggests a slowly increasing contribution of nuclei heavier than protons and alpha-particles to PCR mass composition. Effective atomic number of PCR particles changes in the energy interval 10^15-10^17 eV so that lnA grows from 2 to 2.5.

Adjusted data on the fraction of protons in the mass composition of primary cosmic rays (PCRs) in the energy range of 10^15-10^17 eV are presented. Adjustments are made according to detailed calculations of the response of the X-ray emulsion chamber in the PAMIR experiment. It is demonstrated that the fraction of protons in a PCR is 16-18% for E0 = 10^15-10^16 eV and does not change within the error for E0 = 10^16-10^17 eV.

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.

Recently essential efforts were made to improve measurement routine with X-ray films exposed in the X-ray emulsion chambers at the Pamirs. Analysis of X-ray emulsion response upon recorded events show that gamma-family energy and intensity in early publications were over estimated. The main physical results of the new analysis are presented.

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.

A search for aligned events has been done throughout the muon group events measured with Baksan Underground Scintillation Telescope (BUST) during a long period of time. Only groups of multiplicity >3 for muon threshold energies equal to 0.25; 0.85; 2.0 and 3.2 TeV were selected for the analysis. A distribution of the events on the alignment parameter /lambda has been obtained and compared with the results of Monte-Carlo simulation made for this experiment. Upper limit for aligned muon groups flux is given.

A search for aligned events has been done throughout the muon groups events measured by Baksan Underground Scintillation Telescope (BUST) during a period of 7.7 years. Only groups of multiplicity $>3$ for muon threshold energy equal to 0.85 TeV were selected for the analysis. A distribution of the events on alignment parameter $\lambda$ has been obtained and compared with the results of Monte-Carlo simulation made for this experiment. The upper limit for aligned muon event flux as low as $5.3\cdot10^{-15}$ $cm^{-2} sec^{-1} sr^{-1}$ is given.