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.

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.

Appreciable hard X-ray (HXR) and gamma-ray emissions in the 0.04 to 150 MeV energy range associated with the 2003 October 29 solar flare (X10/3B) were observed at 20:38 to 20:58 UT by the SONG instrument onboard the CORONAS-F mission. To restore flare gamma-ray spectra we fitted the SONG energy loss spectra with a three-component model of the incident spectrum: (1) a power law in energy, assumed to be due to electron bremsstrahlung; (2) a broad continuum produced by prompt nuclear de-excitation gamma-lines; and (3) a broad gamma-line generated from pion decay. We also restored spectra from the RHESSI data, compared them with the SONG spectra and found a reasonable agreement between these spectra in the 0.1 to 10 MeV energy range. The pion-decay emission was observed from 20:44:20 UT and had its maximum at 20:48 to 20:51 UT. The power-law spectral index of accelerated protons estimated from the ratio between intensities of different components of gamma rays changed with time. The hardest spectrum with a power-law index S = - 3.5 to - 3.6 was observed at 20:48 to 20:51 UT. Time histories of the pion-decay emission and proton spectrum were compared with changes of the locations of flare energy release as shown by RHESSI hard X-ray images and remote H alpha brightening. An apparent temporal correlation between processes of particle acceleration and restructuring of flare magnetic field was found. In particular, the protons were accelerated to sub-relativistic energies after radical change of the character of footpoint motion from a converging motion to a separation motion.

Basic ideas of muon tracker technique for the solution of primary cosmic ray (PCR) composition problem in the energy range $10^{17}-10^{18}$ eV are presented. The approach uses MC simulation data made with CORSIKA6.990 for "Pamir-XXI" site conditions. Similar technology can certainly be developed for other observation levels and interaction models. One can probably extend it to much higher primary energies.

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.

In the present paper we consider neutrino events due to quasi-elastic scattering (QEL) as the most reliable events among various candidate events to be analyzed, and have carried out the first step of an L/E analysis which aims to confirm the survival probability with a Numerical Computer Experiment. The most important factor in the survival probability is Lnu and Enu, but this cannot be measured for such neutral particles. Instead, Lmu and Emu is utilized in the L/E analysis, where Lnu, Lmu, Enu and Emu denote the flight path lengths of the incident neutrinos, those of the emitted leptons, the energies of the incident neutrinos and those of the emitted leptons, respectively. According to our Computer Numerical Experiment, the relation of Lnu/Enu is nearly equal to Lmu/Emu doesn't hold. In subsequent papers, we show the results on an L/E analysis with the Computer Numerical Experiment based on our results obtained in the present paper.

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

We retrieved the total content of the atmospheric water vapor (or Integrated Water Vapor, IWV) from extensive sets of photometric data obtained since 1995 at Lindenberg Meteorological Observatory with star and sun photometers. Different methods of determination of the empirical parameters that are necessary for the retrieval are discussed. The instruments were independently calibrated using laboratory measurements made at Pulkovo Observatory with the VKM-100 multi-pass vacuum cell. The empirical parameters were also calculated by the simulation of the atmospheric absorption by water vapor, using the MODRAN-4 program package for different model atmospheres. The results are compared to those presented in the literature, obtained with different instruments and methods of the retrieval. The reliability of the empirical parameters, used for the power approximation that links the water vapor content with the observed absorption, is analyzed. Currently, the total (from measurements, calibration, and calculations) errors yield the standard uncertainty of about 10% in the total column water vapor. We discuss the possibilities for improving the accuracy of calibration to ~1% as indispensable condition in order to make it possible to use data obtained by optical photometry as an independent reference for other methods (GPS, MW-radiometers, lidar, etc).

We describe the laboratory complex for the calibration of photometers that are used in weather service to measure the water vapor content in the Earth atmosphere. The complex was built up in Pulkovo Observatory and developed within the framework of collaboration between Pulkovo Observatory and Lindenberg Meteorological Observatory (Meteorologisches Observatorium Lindenberg - Richard-A{\ss}mann-Observatorium, Lindenberg, Germany). It is used to obtain calibration dependences for individual devices, and also to develop and compare various methods of construction of calibration dependences. These techniques are based on direct calibration of the photometers, on the use of spectral laboratory transmission functions for water vapor, on calculation methods using spectroscopical databases for individual lines. We hope that when the parameters of the equipment are taken into account in detail and new results for the absorptive power of water vapor are used, the accuracy of determination of the water vapor content in the atmosphere of 1-2% may be attained.

In spectral studies of water vapor under laboratory conditions (determination of molecular constants, measurement for spectral transmission functions), the amount of water vapor in the time of the measurements is one of the most essential parameters, which should be determined accurately. We discuss the application for this purpose of polymer sensors of humidity manufactured by Praktik-NC (Moscow) and used in the Pulkovo VKM-100 multipass vacuum cell. These sensors were examined in the laboratory of Lindenberg Meteorological observatory (Germany) by comparison between their readings and those of standard measuring devices for various values of relative humidity, pressure, and temperature. We also carried out measurements of relative humidity in boxes with saline solution, in which the relative humidity that corresponds to a given solution is guaranteed with the accuracy of several tenths of percent. The analysis of the results of the laboratory examination of the sensors and extended sets of measurements made with the Pulkovo cell made it possible to conclude that in measurements in the interval of relative humidity 40-80%, the ~5% accuracy of the measurements for the water vapor content is reached. Further paths are indicated for the increase of the accuracy of measurements and extending the interval of the relative humidity, in which accurate measurements may be carried out.

We obtained laboratory spectra of absorption by water vapor at the wavelengths 6500-10500 {\AA} with the multipass cell. The water vapor content along the line of view varied from 0.1 to 3.0 cm of precipitated water, the pressure from 0.1 to 1.0 atm. The spectra were taken with the width of the exit slit of the spectrophotometer 25, 50, 100, and 150 {\AA}. To match these spectra, we selected empirical functions, which approximate the observed absorption within the indicated interval of water vapor content and pressure with the accuracy about 1%. For the water vapor band at the wavelengths regions 7200, 8200, and 9300 {\AA}, with the step 25 {\AA}, we determined the parameters necessary for the calculation of empirical transmission functions. The presented data make it possible to select the parameters for taking into account the radiation attenuation in the spectral region of telluric water vapor under the conditions of real astronomical observations for a specific place and spectrophotometer. The suggested set of empirical parameters may provide correction of observed stellar spectra for the extinction in the atmosphere with the accuracy 0.m01-0.m02.

On the basis of experience acquired at creation of the Pulkovo Spectrophotometric Catalog the method of investigation of a terrestrial atmospheric components (aerosols and water vapor) in night time are designed. For these purposes the small-sized photometers were created. Carried out in 1995-1999{\Gamma}.{\Gamma}. series of night and daily monitoring of the atmospheric condition in Pulkovo, in MGO by A.I.Voejkov., in Germany (complex experiments LITFASS 98 and LACE 98) confirmed suitability of devices, techniques of observations and their reduction designed in Pulkovo Observatory for the solution of geophysical and ecological problems. A final aim of this work - creation of small-sized automatic complexes (telescope + photometer), which would be rightful component of meteorological observatories. Such complexes will work without the help of the observer and would provide the daily monitoring of a terrestrial atmosphere.

During summer of 1998 two large-scale complex campaigns, LITFASS98 (May 25th to June 22nd) and LACE98 (July 13th to August 14th), took place at the Meteorological Observatory Lindenberg (MOL). The aim of both experiments focus on the intensive daily observations of atmospheric conditions and the determination of their fundamental meteorological parameters in the vertical column over Lindenberg (Lindenberg's Column). About 20 German research institutions and addition one from the Netherlands, Austria and Russia participated at the experiments. A wide variety of ground-based instruments was operated in Lindenberg and Falkenberg, including LIDARs, microwave radiometer and radiosondes complemented by tethered balloons and aircraft measurements. For the first time the star- and sunphotometer of MOL were used together with other geophysical tools. The observations with both photometers were carried out practically every day and night except during absolutely overcast conditions. The observed data were processed immediately by a series of programs developed at Pulkovo Observatory (Russia), and the results (daily variations of aerosol optical depths and water vapour column content) were presented at daily briefings. The comparison of these results with radiosonde and microwave radiometer data demonstrated the usefulness of photometer data for the calibration of other ground-based observations and satellite measurements.

This paper describes the accuracy and the errors of water vapour content measurements in the atmosphere using optical methods, especially starphotometer. After the general explanations of the used expressions for the star-magnitude observations of the water vapour absorption in section 3 the absorption model for the water vapour band will be discussed. Sections 4 and 5 give an overview on the technique to determine the model parameters both from spectroscopic laboratory and radiosonde observation data. Finally, the sections 6 and 7 are dealing with the details of the errors; that means errors of observable magnitude, of instrumental extraterrestrial magnitude, of atmospheric extinction determination and of water vapour content determination by radiosonde humidity measurements. The main conclusion is: Because of the high precision of the results the optical methods for water vapour observation are suited to validate and calibrate alternative methods (GPS, LIDAR, MICROWAVE) which are making constant progress world-wide in these days.

In the previous paper (Part~1), we have verified that the SK assumption on the direction does not hold in the analysis of neutrino events occurred inside the SK detector. We have made four possible L/E analyses, L_nu/E_nu, L_nu/E_mu, L_mu/E_nu and L_mu/E_mu. Among four kinds of L/E analyses, we have shown that only L_nu/E_nu analysis can give the signature of maximum oscillations clearly, while the L_mu/E_mu analysis which are really done by Super-Kamiokande Collaboration cannot give the maximum oscillation at all. It is thus concluded that Super-Kamiokande type experiment cannot find the maximum oscillation from L/E analysis. Therefore, we would suggest Super-Kamiokande Collaboration to re-analyze the zenith angle distribution of the neutrino events which occur inside the detector carefully.

Quasi Elastic Scattering (QEL) is the dominant mechanism for producing both Fully Contained Events and Partially Contained Events in the SuperKamiokande (SK) detector for atmospheric neutrinos in the range 0.1GeV to 10GeV. In the analysis of SK events, it is assumed that the zenith angle of the incident neutrino is the same as that of the detected charged lepton. In the present paper, we derive the distribution function for the scattering angle of charged leptons and derive the average scattering angles as well as their standard deviation due to QEL. Thus, it is shown that the SK assumption for the scattering angle of the charged lepton in QEL is not valid. Further, for a real analysis of the experimental data, we should examine the correlations between the zenith angle of the charged leptons and those of their incident neutrinos, taking into account the effect of the azimuthal angle of the charged particle over its zenith angle. We conclude that it is not possible to reliably assign the zenith angle of the incident neutrino to that of the charged lepton, and so the zenith angle distribution of charged leptons does not reflect that of the incident neutrinos under the SK assumption. This result has clear implication for attempts to detect neutrino oscillations from an analysis of Fully Contained Events and Partially Contained Events in SuperKamiokande.

Quasi Elastic Scattering (QEL) is the dominant source for producing both Fully Contained Events and Partially Contained Events in the Superkamiokande(SK) detector for the atmospheric neutrinos, in the range 0.1 GeV to 10 GeV. In the analysis of SK events, it is assumed that the direction of the incident neutrino is the same as that of the detected charged lepton. In the present letter, we derive the distribution function for the scattering angle of the charged leptons, their averaged scattering angle and their standard deviation due to QEL. Then, it is shown that the SK assumption for the scattering angle of the charged leptons in the QEL is not valid. Further, we examine the influence of the azimuthal angle of the charged leptons over their zenith angle. As the result, we conclude that the zenith angle distribution of the neutrino under the SK assumption does not reflect the real zenith angle distribution of the atmospheric neutrino which produces Fully Contained Events and Partially Contained Events. This result has clear implication for attempts to detect neutrino oscillations from the analyses of Fully Contained Events and Partially Contained Events in Superkamiokande.