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.

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.

It has been shown that muon flux intensities calculated in terms of the EPOS LHC and EPOS 1.99 models at the energy of 10^4 GeV exceed the data of the classical experiments L3+Cosmic, MACRO and LVD on the spectra of atmospheric muons by a factor of 1.9 and below these data at the same energy by a factor of 1.8 in case of the QGSJET II-03 model. It has been concluded that these tested models overestimate (underestimate in case of QGSJET II-03 model) the production of secondary particles with the highest energies in interactions of hadrons by a factor of ~1.5. The LHCf and TOTEM accelerator experiments show also this type of disagreements with these model predictions at highest energies of secondary particles.

A simple method of the vertical muon energy spectrum simulations have been suggested. These calculations have been carried out in terms of various models of hadronic interactions. The most energetic $ \pi^\pm $-mesons and K$^\pm $-mesons produced in hadron interactions contribute mainly in to this energy spectrum of muons due to the very steep energy spectrum of the primary particles. So, some constraints on the hadronic interaction models may be set from a comparison of calculated results with the cosmic data on the vertical muon energy spectrum. This comparison showed that the most energetic secondary particles production is too high in case of the QGSJET II-04 model and rather low in case of the QGSJET II-03 model. These conclusion have been supported by the LHC data.

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.

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

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

We present the interpretation of the muon and scintillation signals of ultra-high-energy air showers observed by AGASA and Yakutsk extensive air shower array experiments. We consider case-by-case ten highest energy events with known muon content and conclude that at the 95% confidence level (C.L.) none of them was induced by a primary photon. Taking into account statistical fluctuations and differences in the energy estimation of proton and photon primaries, we derive an upper limit of 36% at 95% C.L. on the fraction of primary photons in the cosmic-ray flux above 10^20 eV. This result disfavors the Z-burst and superheavy dark-matter solutions to the GZK-cutoff problem.