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 matrices formed of silicon photomultipliers (SiPMs) are very promising photosensors for the scintillation detectors. The use of SiPM matrices with appropriate optical collector gives, in principle, a possibility to do a snapshot of glowing track of charged particle traversing a scintillator. The prototype of such scintillation detector is under development now in INR RAS. The preliminary results of characterization study of the matrix ArrayC-60035-64P-PCB (SensL company) for the prototype of such detector 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.

According to the Baksan underground scintillation telescope data the barometric coefficients for six intervals of zenith angles of particle registration were obtained for ten years (2003-2012). The obtained barometric coefficients exceed an order of magnitude the theoretical values. Nevertheless, these results are in good agreement with the data obtained for the barometric effect of high-energy muons in a number papers of other authors. And there is a noticeable increase of the absolute values of the measured barometric coefficients with increase of the threshold energy of muons.

The multichannel data acquisition system is intended to be used in the EMMA experiment studying cosmic rays. The array will be in the Pihasalmi mine (central Finland) at a depth of about 85 m. The scintillator counters (SC-1) of the array are cast plastic scintillators with a wavelength of each SC-1 detector is 12.2{\times}12.2{\times}3.0 cm^3. 16 SC-1 detectors are placed in the metal case of 50.0{\times}50.0{\times}13.0 cm^3 dimension. Each case, called SC-16 detector, contains electronics of preliminary processing of signals and operating mode stabilization. The whole of the array will contain 96 detectors SC-16. It will make 1536 channels placed in three planes (48+24+24 detectors). The array will allow us to measure the time of flight of particles between SC-16 detectors and the coordinates of SC-1 fired detectors. This paper presents the function diagram of data acquisition system that includes electronic of detectors, the hodoscope pulse channels, the trigger block and VME blocks.

Current status and results of the experiment on recording neutrino bursts are presented. The observation time (since 1980) is 24.7 years. The upper bound of collapse frequency in our Galaxy is 0.093 $y^{-1}$ (90% CL).

Experimental data of arrays "Andyrchy" and "Carpet-2" of Baksan Neutrino Observatory (Institute for Nuclear Research), obtained in the regime of a detection of the single cosmic-ray component, are used for a search of the bursts of cosmic gamma rays from evaporating primordial black holes. Different theoretical models of the evaporation process are used for the analysis. Distributions of the counting rate fluctuations on both arrays agree with the expectations from the cosmic ray background. The new constraints on the concentration of evaporating primordial black holes in the local region of Galaxy are obtained. The comparison of the results of different experiments is given.

The technical procedure of a search of bursts of very high energy gamma rays from evaporation of primordial black holes on air-shower array "Andyrchy" of Baksan Neutrino Observatory of Institute for Nuclear Research is described. The theoretical model used in the present work assumes that the chromosphere around the evaporating black hole does not form. For minimization of the cosmic ray background the method of multidimensional analysis of modelled as well as experimentally detected events is applied. The new upper limit on the concentration of evaporating primordial black holes in the local region of Galaxy is obtained. The comparison of the results of different experiments is given.

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.

Temporal and energy characteristics of the very-high-energy gamma-ray bursts from evaporating primordial black holes have been calculated by assuming that the photospheric and chromospheric effects are negligible. The technique of searching for such bursts on shower arrays is described. We show that the burst time profile and the array dead time should be taken into account to interpret experimental data. Based on data from the Andyrchy array of the Baksan Neutrino Observatory (Institute for Nuclear Research, Russian Academy of Sciences), we have obtained an upper limit on the number density of evaporating primordial black holes in a local region of space with a scale size of ~10^{-3} pc. Comparison with the results of previous experiments is made.