A new project is developed with the implementation of a relatively new method of studying the primary cosmic ray -- the registration of extensive air showers' optical Vavilov-Cherenkov radiation (Cherenkov light) reflected from the snow surface. The aim of the project is the study of the cosmic ray mass composition in the energy range of 1-100 PeV by detecting the reflected extensive air showers' Cherenkov light. Silicon photomultipliers are planned to be used as the main photosensitive element of the detector and an unmanned aerial vehicle will is planned to lift the measuring equipment over the snow-covered ground.

The SPHERE-2 balloon-borne detector designed for extensive air shower (EAS) observations using EAS optical Vavilov-Cherenkov radiation (``Cherenkov light''), reflected from the snow-covered surface of Lake Baikal is described. We briefly discuss the concept behind the reflected Cherenkov light method, characterize the conditions at the experimental site and overview the construction of the tethered balloon used to lift the SPHERE-2 telescope above the surface. This paper is mainly dedicated to a detailed technical description of the detector, including its optical system, sensitive elements, electronics, and data acquisition system (DAQ). The results of some laboratory and field tests of the optical system are presented.

A compact device lifted over the ground surface might be used to observe optical radiation of extensive air showers (EAS). Here we consider spatial and temporal characteristics of Vavilov-Cherenkov radiation ("Cherenkov light") reflected from the snow surface of Lake Baikal, as registered by the SPHERE-2 detector. We perform detailed full direct Monte Carlo simulations of EAS development and present a dedicated highly modular code intended for detector response simulations. Detector response properties are illustrated by example of several model EAS events. The instrumental acceptance of the SPHERE-2 detector was calculated for a range of observation conditions. We introduce the concept of "composite model quantities", calculated for detector responses averaged over photoelectron count fluctuations, but retaining EAS development fluctuations. The distortions of EAS Cherenkov light lateral distribution function (LDF) introduced by the SPHERE-2 telescope are understood by comparing composite model LDF with the corresponding function as would be recorded by an ideal detector situated at the ground surface. We show that the uncertainty of snow optical properties does not change our conclusions, and, moreover, that the expected performance of the SPHERE experiment in the task of cosmic ray mass composition study in the energy region $\sim$10 PeV is comparable with other contemporary experiments. Finally, we compare the reflected Cherenkov light method with other experimental techniques and briefly discuss its prospects.

Preliminary results on the development of a separation method for Cerenkov (CL) and fluorescence (FL) light from EAS are shown. The results are based on the measurement of attenuation coefficients of CL and FL for different filters. A total of six optical filters were investigated: filters from optical glass UFS-1, UFS-5, FS6 (analogue BG3) and interference filters SL 360-50, SL 280-380, FF01-375/110. The measurements were performed using silicon photomultipliers (SiPM). To improve existing fluorescent light detectors, a segment of 7 SiPM was developed, which would be able to separate both components of the light flux from EAS at the level of primary data processing.

Local muon density spectra (LMDS) at various zenith angles have been reconstructed from the data of two detectors of the Experimental complex NEVOD. The inclined muon bundles at the ground level were detected with the coordinate detector DECOR, and for the near-vertical direction with the calibration telescope system (CTS) of the Cherenkov water detector. In comparison with the earlier DECOR results, the experimental statistics has been increased by 2-3 times for different ranges of zenith angle and muon bundle multiplicity and is now based on about 40,000 h of the setup operation. The live time of measurements with CTS is about 12,000 h. It is found with both setups that the slope of LMDS is increasing above the primary energy of about 1017 eV. The details of the experiment and data analysis are presented.

The spatial distribution of Cherenkov radiation from cascade showers generated by muons in water has been measured with Cherenkov water calorimeter (CWC) NEVOD. This result allowed to improve the techniques of treating cascade showers with unknown axes by means of CWC response analysis. The techniques of selecting the events with high energy cascade showers and reconstructing their parameters are discussed. Preliminary results of measurements of the spectrum of cascade showers in the energy range 100 GeV - 20 TeV generated by cosmic ray muons at large zenith angles and their comparison with expectation are presented.

Results of investigations of the near-horizontal muons in the range of zenith angles of 85-95 degrees are presented. In this range, so-called "albedo" muons (atmospheric muons scattered in the ground into the upper hemisphere) are detected. Albedo muons are one of the main sources of the background in neutrino experiments. Experimental data of two series of measurements conducted at the experimental complex NEVOD-DECOR with the duration of about 30 thousand hours "live" time are analyzed. The results of measurements of the muon flux intensity are compared with simulation results using Monte-Carlo on the basis of two multiple Coulomb scattering models: model of point-like nuclei and model taking into account finite size of nuclei.

Change of the energy characteristics of muon bundles with an increase of the primary cosmic ray particles energy can be a key to solving the problem of muon excess in the extensive air showers (EAS) observed in a number of experiments. In this work the data on the energy deposit of multi-muon events in a wide range of zenith angles (and as a consequence in a wide range of primary particles energies) obtained with NEVOD-DECOR setup over a long time period are presented. The experimental data are compared with the results of simulations of EAS muon component performed using CORSIKA code.

We present an event-by-event study of cosmic ray (CR) composition with the reflected Cherenkov light method. The fraction of CR light component above 5 PeV was reconstructed using the 2013 run data of the SPHERE experiment which observed optical Vavilov-Cherenkov radiation of extensive air showers, reflected from snow surface of Lake Baikal. Additionally, we discuss a possibility to improve the elemental groups separability by means of multidimensional criteria.

We give an overview of the SPHERE experiment based on detection of reflected Vavilov-Cherenkov radiation (Cherenkov light) from extensive air showers in the energy region E>10^{15} eV. A brief history of the reflected Cherenkov light technique is given; the observations carried out with the SPHERE-2 detector are summarized; the methods of the experimental datasample analysis are described. The first results on the primary cosmic ray all-nuclei energy spectrum and mass composition are presented. Finally, the prospects of the SPHERE experiment and the reflected Cherenkov light technique are given.

Here is presented the current state of the SPHERE-2 balloon-borne experiment. The detector is elevated up to 1 km above the snow surface and registers the reflected Vavilov-Cherenkov radiation from extensive air showers. This method has good sensitivity to the mass-composition of the primary cosmic rays due to its high resolution near the shower axis. The detector consists of a 1500 mm spherical mirror with a 109 PMT cluster in its focus. The electronics record a signal pulse profile in each PMT. In the last 2 years the detector was upgraded: time resolution of pulse registration was enhanced up to 12.5 ns, channel sensitivity was increased by a factor of 3, a new LED-based relative PMT calibration method was introduced, and new hardware and etc. was installed.

First preliminary results of the balloon-borne experiment SPHERE-2 on the all-nuclei primary cosmic rays (PCR) spectrum and primary composition are presented. The primary spectrum in the energy range $10^{16}$--$5\cdot10^{17}$ eV was reconstructed using characteristics of Vavilov-Cherenkov radiation of extensive air showers (EAS), reflected from a snow surface. Several sources of systematic uncertainties of the spectrum were analysed. A method for separation of the primary nuclei' groups based on the lateral distribution function' (LDF) steepness parameter is presented. Preliminary estimate of the mean light nuclei' fraction $f_{30-150}$ at energies $3\cdot10^{16}$--$1.5\cdot10^{17}$ eV was performed and yielded $f_{30-150}$= (21$\pm$11)%.

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.

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.

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 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.