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.

We report on a measurement of thermal neutrons, generated by the hadronic component of extensive air showers (EAS), by means of a small array of EN-detectors developed for the PRISMA project (PRImary Spectrum Measurement Array), novel devices based on a compound alloy of ZnS(Ag) and $^{6}$LiF. This array has been operated within the ARGO-YBJ experiment at the high altitude Cosmic Ray Observatory in Yangbajing (Tibet, 4300 m a.s.l.). Due to the tight correlation between the air shower hadrons and thermal neutrons, this technique can be envisaged as a simple way to estimate the number of high energy hadrons in EAS. Coincident events generated by primary cosmic rays of energies greater than 100 TeV have been selected and analyzed. The EN-detectors have been used to record simultaneously thermal neutrons and the air shower electromagnetic component. The density distributions of both components and the total number of thermal neutrons have been measured. The correlation of these data with the measurements carried out by ARGO-YBJ confirms the excellent performance of the EN-detector.

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.

The measurement of large scale anisotropies in cosmic ray arrival directions at energies above 10^13 eV is performed through the detection of Extensive Air Showers produced by cosmic ray interactions in the atmosphere. The observed anisotropies are small, so accurate measurements require small statistical uncertainties, i.e. large datasets. These can be obtained by employing ground detector arrays with large extensions (from 10^4 to 10^9 m^2) and long operation time (up to 20 years). The control of such arrays is challenging and spurious variations in the counting rate due to instrumental effects (e.g. data taking interruptions or changes in the acceptance) and atmospheric effects (e.g. air temperature and pressure effects on EAS development) are usually present. These modulations must be corrected very precisely before performing standard anisotropy analyses, i.e. harmonic analysis of the counting rate versus local sidereal time. In this paper we discuss an alternative method to measure large scale anisotropies, the "East-West method", originally proposed by Nagashima in 1989. It is a differential method, as it is based on the analysis of the difference of the counting rates in the East and West directions. Besides explaining the principle, we present here its mathematical derivation, showing that the method is largely independent of experimental effects, that is, it does not require corrections for acceptance and/or for atmospheric effects. We explain the use of the method to derive the amplitude and phase of the anisotropy and we demonstrate its power under different conditions of detector operation.

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.

The amplitude and phase of the cosmic ray anisotropy are well established experimentally between 10^{11} eV and 10^{14} eV. The study of their evolution into the energy region 10^{14}-10^{16} eV can provide a significant tool for the understanding of the steepening ("knee") of the primary spectrum. In this letter we extend the EAS-TOP measurement performed at E_0 around 10^{14} eV, to higher energies by using the full data set (8 years of data taking). Results derived at about 10^{14} and 4x10^{14} eV are compared and discussed. Hints of increasing amplitude and change of phase above 10^{14} eV are reported. The significance of the observation for the understanding of cosmic ray propagation is discussed.