The Large High Altitude Air Shower Observatory (LHAASO) project is a new generation multi-component instrument, to be built at 4410 meters of altitude in the Sichuan province of China, with the aim to study with unprecedented sensitivity the spec trum, the composition and the anisotropy of cosmic rays in the energy range between 10$^{12}$ and 10$^{18}$ eV, as well as to act simultaneously as a wide aperture (one stereoradiant), continuously-operated gamma ray telescope in the energy range between 10$^{11}$ and $10^{15}$ eV. The experiment will be able of continuously surveying the TeV sky for steady and transient sources from 100 GeV to 1 PeV, t hus opening for the first time the 100-1000 TeV range to the direct observations of the high energy cosmic ray sources. In addition, the different observables (electronic, muonic and Cherenkov/fluorescence components) that will be measured in LHAASO will allow to investigate origin, acceleration and propagation of the radiation through a measurement of energy spec trum, elemental composition and anisotropy with unprecedented resolution. The remarkable sensitivity of LHAASO in cosmic rays physics and gamma astronomy would play a key-role in the comprehensive general program to explore the High Energy Universe. LHAASO will allow important studies of fundamental physics (such as indirect dark matter search, Lorentz invariance violation, quantum gravity) and solar and heliospheric physics. In this document we introduce the concept of LHAASO and the main science goals, providing an overview of the project.

Measuring radon-due neutron flux at various altitude (100, 1000, 1700, 4300 m above sea level) we found an evidence of significant increase of radon concentration with altitude. It was also conirmed by direct radon measurements at high altitude. This allowed us to assume cosmic rays could take part in process of soil activation: they transform long-lived nuclei of uranium and thorium to nuclei with shorter life-time through specific nuclear reactions. If the resulting nuclei belong to the U-238 radioactive chain they can lead to producion of Ra-226 and then to Rn-222, thus significantly increasing its production at high altitudes where cosmic ray flux is high.

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.

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.

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.

It is shown that a "second (iron) knee" in cosmic ray spectrum expected at energy about 100 PeV could not probably be found there. The reason is very simple: the position of the "iron knee" depends on an answer to the questions: "What do we see at 3-5 PeV? Is the knee seen at this energy associated with proton or iron primaries?"

Preliminary results of one year anisotropy measurement in the energy range 10^{13} -10^{14} eV as a function of energy are presented. The results are compared for two methods of data analysis: the standard one with meteo correction approach in use and another one so-called "East minus West" method. Amplitudes and phases of anisotropy for three median energies E = 25 TeV, E = 75 TeV and E = 120 TeV are reported. Brief consideration of amplitude-phase dependence of anisotropy on energy is expounded.

A novel type of Extensive Air Shower (EAS) array is proposed and described. It is shown that only new approaches to the so called "knee problem" could solve this complicated and old problem.

It is shown that Extensive Air Shower (EAS) longitudinal development has a critical point where an equilibrium between the main hadronic component and the secondary electromagnetic one exhibits a brake. This results in a change of slope in quasi-power law function $N_{e}(Eo)$. The latter leads to a knee in the EAS size spectrum at primary energy of about 100 TeV/nucleon. Many ``strange'' experimental results can be successfully explained in the frames of current approach.

Examples of the power law spectra observations in the Nature are given. It is shown that the observed meteoroid mass spectrum is very similar to the observed cosmic ray energy spectrum. It is concluded that the "knees" observed in both spectra have a similar origin: the Earth's atmosphere. Both spectra are studied by indirect methods through secondary components and specific behavior of these components in the atmosphere could produce the "knee" even in a case when primary spectrum follows a pure power law.

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.

Two approaches to the old problem of the so-called "knee" in cosmic ray spectrum are compared with the existing experimental data. It is shown that "standard astrophysical model" of the knee is not supported by some experimental data, while the alternative model agrees better with the data.

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.

All astrophysical explanations of the knee in cosmic ray spectrum accept the hypothesis of its existence {\it a priori} without any doubt. But, there exist experimental evidences against this hypothesis. Experimental data on the knee in various secondary cosmic ray components obtained by KASCADE array as well as by many other arrays reveal the existence of the knee ($\Delta\gamma\approx0.4\div0.5$) in electromagnetic component (e and $\gamma$ including Cherenkov light) and the absence of it or a very small knee ($\Delta\gamma\approx0.1\div0.2$) in muonic and hadronic components. But, one could expect the inverse relationship in case of the astrophysical knee origin. A brief review of experimental data concerning this problem is given.