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.

Direct Dark Matter searches are nowadays one of the most fervid research topics with many experimental efforts devoted to the search for nuclear recoils induced by the scattering of Weakly Interactive Massive Particles (WIMPs). Detectors able to reconstruct the direction of the nucleus recoiling against the scattering WIMP are opening a new frontier to possibly extend Dark Matter searches beyond the neutrino background. Exploiting directionality would also prove the galactic origin of Dark Matter with an unambiguous signal-to-background separation. Indeed, the angular distribution of recoiled nuclei is centered around the direction of the Cygnus constellation, while the background distribution is expected to be isotropic. Current directional experiments are based on gas TPC whose sensitivity is limited by the small achievable detector mass. In this paper we present the discovery potential of a directional experiment based on the use of a solid target made of newly developed nuclear emulsions and of optical read-out systems reaching unprecedented nanometric resolution.

Nowadays there is compelling evidence for the existence of dark matter in the Universe. A general consensus has been expressed on the need for a directional sensitive detector to confirm, with a complementary approach, the candidates found in conventional searches and to finally extend their sensitivity beyond the limit of neutrino-induced background. We propose here the use of a detector based on nuclear emulsions to measure the direction of WIMP-induced nuclear recoils. The production of nuclear emulsion films with nanometric grains is established. Several measurement campaigns have demonstrated the capability of detecting sub-micrometric tracks left by low energy ions in such emulsion films. Innovative analysis technologies with fully automated optical microscopes have made it possible to achieve the track reconstruction for path lengths down to one hundred nanometers and there are good prospects to further exceed this limit. The detector concept we propose foresees the use of a bulk of nuclear emulsion films surrounded by a shield from environmental radioactivity, to be placed on an equatorial telescope in order to cancel out the effect of the Earth rotation, thus keeping the detector at a fixed orientation toward the expected direction of galactic WIMPs. We report the schedule and cost estimate for a one-kilogram mass pilot experiment, aiming at delivering the first results on the time scale of six years.

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

Here the results of calculations of pulses in Cherenkov light detectors for the Yakutsk array are presented. As long as the Vavilov-Cherenkov light is used to calibrate signals in scintillation detectors at the Yakutsk array it is vital for these measurements to be precise. The validity of measurements of the signal Q(400) used as the estimation parameter at the Yakutsk array has been confirmed. Our calculations show that the width of time pulses increases from nearly 50 ns at a distance of 100 m from the shower axis up to 700 ns at 1000 m.

The fluxes of electrons, positrons, gammas, Cherenkov photons and muons in individual extensive air showers induced by the primary protons and helium, oxygen and iron nuclei at the level of observation have been estimated with help of the code CORSICA 6.616. The comparison show that the values of the function Xi**2 per one degree of freedom changes from 1.1 for iron nuclei to 0.9 for primary protons. As this difference is small all readings of detectors of the Vavilov-Cherenkov radiation have been used. At last, readings of underground detectors of muons with energies above 1 GeV have been exploited to make definite conclusion about chemical composition.

The energies of the most energetic extensive air showers observed at the Yakutsk array have been estimated with help of the all detectors readings instead of using of the standard procedure with a parameter s(600). The energy of the most energetic extensive air shower observed at the Yakutsk array happened to be 200, 200, 180 and 165 EeV with the values of the Xi**2 function per one degree of freedom 0.9, 1., 0.9 and 1.1 for the primary protons and helium, oxygen and iron nuclei accordingly.

The primary cosmic rays particles with energies above 10**20 eV have been observed at many extensive air shower arrays since the beginning of observations over 40 years ago. The validity of measurements of signal s(600) used as energy estimation parameter at the Yakutsk array has been confirmed. Our calculations show that the width of the time pulses increases from nearly 100 ns at a distance of 100 m from the shower axis up to 4 - 5 $\mu$s at 1500 m. The calculated estimate of energy of extensive air shower is ~ 1.7 times smaller than the experimental estimate for the same value s(600). The pointing directions of extensive air showers observed at the Pierre Auger Observatory were fitted within +-3.1**o with positions of the nearby active galactic nuclei from the Veron-Cetty and P. Veron catalog. The cosmic ray luminosity of the active galactic nuclei which happened to be a source of the particular cosmic ray event constitutes a fraction ~10**-4 of the optical one if only cosmic ray particles with energies above 6*10**19 eV are produced. If produced cosmic ray particles have a spectrum up to ~ 100 GeV then the cosmic ray luminosity of the active galactic nuclei should be much higher than the optical one.

The pointing directions of extensive air showers observed at the Pierre Auger Observatory were fitted within 3.1 degree with positions of the nearby active galactic nuclei from the Veron-Cetty and P. Veron catalog. The cosmic ray luminosity of the active galactic nuclei which happened to be a source of the particular cosmic ray event constitutes a fraction ~0.0001 of the optical one if only cosmic ray particles with energies above 60 EeV are produced. If produced cosmic ray particles have a spectrum dE/E^3 up to ~100 GeV then the cosmic ray luminosity would be much higher than the optical one of the active galactic nuclei.

In terms of the quark-gluon string model the analysis of the classic procedure to estimate the energy of giant air showers with help of the parameter s(600) (a density of energy deposition in the scintillator at a distance of 600 m from the shower core) have been carried out. Simulations of the signal s(600) with help of the CORSIKA code in terms of the hybrid scheme show energy estimates which are approximately a factor of 1.6 times lower than adopted at the Yakutsk array. The energy estimates calculated with the help of the Cherenkov radiation coincide with the experimental data. Simulations of deposited energy distributions in the atmosphere with help of the GEANT4 code and the CORSIKA code show that more than 20% of this energy may be deposited at distances above 100 m from the shower axis.