The simulation of the Cherenkov light lateral distribution function (LDF) in Extensive Air Showers (EAS) was perfomed with the CORSIKA code in the energy range 10^13-10^16 eV for configuration of the Tunka-25 EAS array. On the basis of this simulation we obtained sets of approximating functions for primary protons, iron nuclei and gamma-quanta for zenith angles Theta<=20. The comparison of the calculated Cherenkov light LDF with that measured with the Tunka-25 array has shown an opportunity of primary particle identification and definition of its energy around the "knee".

We present results of an improved analysis of the experimental data of the EAS Cherenkov array Tunka-25. A new function to fit the Cherenkov light lateral distribution LDF at core distances from 0 to 350 m has been developed on the base of CORSIKA simulations and applied to the analysis of Tunka data. Two methods to estimate the EAS maximum position have been used. The one is based on the pulse FWHM, the other on the light LDF. We present the primary energy spectrum in the energy range 10^15 - 10^17 eV. The use of the depth of the EAS maximum to determine the mean mass composition is discussed.

The project of an EAS Cherenkov array in the Tunka valley/Siberia with an area of about 1 km2 is presented. The new array will have a ten times bigger area than the existing Tunka-25 array and will permit a detailed study of the cosmic ray energy spectrum and the mass composition in the energy range from 10^15 to 10^18 eV.

The project of an EAS Cherenkov array in the Tunka valley/Siberia with an area of about 1 km^2 is presented. The new array will have a ten times bigger area than the existing Tunka-25 array and will permit a detailed study of the cosmic ray energy spectrum and the mass composition in the energy range from 10^15 to 10^18 eV.

New results of 300 hours of operation of the Tunka array are presented. An improved parametrization of the Cherenkov light lateral distribution function (LDF), based on CORSIKA Monte Carlo simulations and the experiment QUEST, has been used for the reconstruction of EAS parameters. The corrected energy spectrum in the knee region is obtained. The mean depth of the EAS maximum has been derived both from the analysis of LDF steepness and the FWHM of Cerenkov light pulse. The mean mass composition around the knee is estimated.

Experiment on the study of fluctuation of night glow is conducted on the base of 13 module atmospheric Cherenkov telescope TUNKA-13.Basic task of TUNKA-13 is registration of light from Extended Air Shower (EAS) initiated by primary cosmic particles. Observation of light flashes radiating in nigght atmosphere is one more experiment which is run for TUNKA-13 to study the background condition for EAS registration. At a period from December 1997 on a May 1998 was registered 149 light flashes of radiwting night atmosphere with duration from 16 ms to 190 seconds. General time of night observing has formed 237 hours. Coincidence of optical flashes with the gamma-bursts from the catalogue of BATSE experiment of observatories "Compton" was analysed.

Supersymmetric models allow for stable non-topological solitons, Q-balls, which can be produced in the early Universe and contribute to dark matter. Experimental signature of electrically neutral Q-balls is, in fact, the same as is expected for superheavy magnetic monopoles catalyzing baryon decay. Here we use the upper limits on monopole flux obtained with deep underwater Cherenkov array "Gyrlyanda" which operated in the Baikal lake in 1984-90 with 267 days of live time to obtain the limit on Q-ball flux. The last has been found to be equal to 3.9 x 10^{-16} cm^{-2} sr^{-1} s^{-1} (90% CL). This result is discussed and compared with other restrictions.