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.

The paper is devoted to the results of the EAS neutron component investigations by means of the PRISMA-32 array. The array consists of 32 en-detectors and enables to record delayed thermal neutrons accompanying showers. For registration of thermal neutrons, the scintillator based on $^{6}Li$ isotope as a target is used in the detectors. Some results of the processing of data accumulated over a long period of time are presented: the lateral distribution function of neutrons in EAS and preliminary results on EAS neutron multiplicity spectrum and distribution of showers in e/n ratio.

Air radon concentration measurement is useful for research on geophysical effects, but it is strongly sensitive to site geology and many geophysical and microclimatic processes such as wind, ventilation, air humidity and so on that induce very big fluctuations on the concentration of radon in air. On the contrary, monitoring the radon concentration in soil by measuring the thermal neutron flux reduces environmental effects. In this paper we report some experimental results on the natural thermal neutron flux as well as the concentration of air radon and its variations at 4300 m a.s.l. These results were obtained with unshielded thermal neutron scintillation detectors (en-detectors) and radon monitors located inside the ARGO-YBJ experimental hall. The correlation of these variations with the lunar month and 1-year period is undoubtedly confirmed. A method for earthquakes prediction provided by a global net of the en-detectors is currently under study.

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.

The results of Monte-Carlo simulations of Extensive Air Shower are presented to show the difference of hadronic component content at various altitudes with the aim to choose an optimal altitude for the PRISMA-like experiment. CORSIKA program for EAS simulations with QGSJET and GHEISHA models was used to calculate the number of hadrons reaching the observational level inside a ring of 50 m radius around the EAS axis. Then the number of neutrons produced by the hadronic component was calculated using an empirical relationship between the two components. We have tested the results with the ProtoPRISMA array at sea level, and recorded neutrons are close to the simulation results.

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.