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.

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 novel type of EAS array (PRISMA-32) has been constructed on the base of NEVOD-DECOR experiment (MEPhI,Moscow) and is now taking data. It consists of 32 specially designed scintillator en-detectors able to measure two main EAS components: hadrons (n) and electrons (e). First results on thermal neutron lateral as well as temporal distributions are presented. Obtained exponential neutron lateral distributions are consistent with that expected for normal hadron production with exponential transverse momentum distribution. As there are no other experimental data on thermal neutron distributions and so, to compare results with other measurements, we additionally obtained electron lateral distribution function (using the same detectors) and compared it with NKG - function. Recorded neutron temporal distributions are very close to that obtained with data of our previous prototypes.