Neutrino air showers can be formed in any part of the atmosphere passing a long was in the matter due to its physical properties. In general, air showers produced by neutrinos are highly inclined and formed near the ground level, i.e. young showers. Therefore, one should expect a large number of peaks in the signal of such air showers [1, 2, 3]. The goal of our work is to search for air shower candidates produced by neutrino. For this purposes, we analyzed large amount of data from scintillation detectors with different area and energy thresholds [1, 2]. Preliminary analysis of Yakutsk array data indicated the absence of air shower produced by neutrino, but it does not mean that such air showers does not exist. It is going to need a further analysis of highly inclined showers. In order to do that improved methodology for recording and processing of air showers required.

We discuss the lateral distribution of charged particles in extensive air showers with energy above $10^{17}$ eV measured by surface scintillation detectors of Yakutsk EAS array. The analysis covers the data obtained during the period from 1977 to 2013. Experimental values are compared to theoretical predictions obtained with the use of CORSIKA code within frameworks of different hadron interaction models. The best agreement between theory and experiment is observed for QGSJet01 and QGSJet-II-04 models. A change in the cosmic ray mass composition towards proton is observed in the energy range $(1-20) \times 10^{17}$ eV.

Responses of surface and underground scintillation detectors of Yakutsk array are calculated for showers initiated by primary particles with energy E0>=1.0E17 eV within the frameworks of QGSJet01 QGSIIJet-II-04, SIBYLL-2.1 and EPOS-LHC hadron interaction models. A new estimation of E0 is obtained with the use of various methods. The resulting energy is lower compared to the obtained with earlier method by factor ~1.33.

We present the results of searches for anisotropy in the right ascension ($RA$) distribution of arrival directions of cosmic rays (CRs) detected with the Yakutsk array during the 1974--2008 observational period in the energy range above $10^{18}$ eV. Two methods of analysis are applied to two sub-samples of the data. Particularly, estimations of the first and second harmonic amplitudes are given, as well as the first harmonic phase in adjacent energy intervals. Analysis of variance demonstrates a significant contraction of the minimal width of the $RA$ distribution in the energy bin $(10^{19},1.78\times10^{19})$ eV with respect to the isotropic distribution, which may be attributed to a possible source of CRs within the interval $RA\in(15^0,45^0)$.

The lateral distribution of muons with ~1.0x\sec(\theta) GeV in extensive air showers within ~10^{17}-10^{19} eV energy region obtained during different observational periods from November 1987 to June 2013 has been analyzed. Experimental data have been compared to predictions of various hadron interaction models. The best agreement is observed with QGSJETII-04. Until 1996, the mass composition of cosmic rays with energy below 2x10^{18} eV was significantly lighter than in later periods.

Here we consider EAS events with energy above 1E17 eV with recorded pulses delayed by t>=5 mcs in scintillation detectors with different thresholds: 10, 5 and 1.8 MeV. In order to identify pulses from electrons, muons and neutrons, experimental data were compared to computational results performed within the framework of QGSJET01d model. Preliminary, one may speculate of registration of low-energy electrons arisen from moderation of neutrons in a detector or a medium surrounding a detector or in the snow cover and frozen crust (albedo particles). The fact that such pulses were registered mostly in low-threshold detectors confirms this hypothesis.

The lateral distribution of muons with threshold energy e ~ 1*sec(theta) GeV have been studied in showers with energy E>=10^17 eV. The data considered in the analysis have been collected from November 2011 to June 2013. Experimental values are compared to predictions obtained with the use of CORSIKA code within the framework of various hadron interaction models. The best agreement between experiment and theory is observed for QGSJETII-04 model. At E~10^17 eV it complies to a mixed cosmic ray composition with the mean atomic number <lnA>~3.0+/-0.5. At E>= 4*10^17 eV the composition varies around the value <lnA>~0.5.

One of several working groups established for this workshop was charged with examining results and methods associated with the UHECR energy spectrum. We summarize the results of our discussions, which include a better understanding of the analysis choices made by groups and their motivation. We find that the energy spectra determined by the larger experiments are consistent in normalization and shape after energy scaling factors are applied. Those scaling factors are within systematic uncertainties in the energy scale, and we discuss future work aimed at reducing these systematics.

The aim of this report of the Working Group on Hadronic Interactions and Air Shower Simulation is to give an overview of the status of the field, emphasizing open questions and a comparison of relevant results of the different experiments. It is shown that an approximate overall understanding of extensive air showers and the corresponding hadronic interactions has been reached. The simulations provide a qualitative description of the bulk of the air shower observables. Discrepancies are however found when the correlation between measurements of the longitudinal shower profile are compared to that of the lateral particle distributions at ground. The report concludes with a list of important problems that should be addressed to make progress in understanding hadronic interactions and, hence, improve the reliability of air shower simulations.

We present estimation of proton-air inelastic interaction cross-section obtained for different energy values are. Results are compared with different hadron interaction models.

A spectrum of cosmic rays within energy range 10^15 - 3x10^17 eV was derived from the data of the small Cherenkov setup, which is a part of the Yakutsk complex EAS array. In this, work a new series of observation is covered. These observations lasted from 2000 till 2010 and resulted in increased number of registered events within interval 10^16 - 10^18 eV, which in turn made it possible to reproduce cosmic ray spectrum in this energy domain with better precision. A sign of a thin structure is observed in the shape of the spectrum. It could be related to the escape of heavy nuclei from our Galaxy. Cosmic ray mass composition was obtained for the energy region 10^16 - 10^18 eV. A joint analysis of spectrum and mass composition of cosmic rays was performed. Obtained results are considered in the context of theoretical computations that were performed with the use of hypothesis of galactic and meta-galactic origin of cosmic rays.

Characteristics of the muon component in EAS are analyzed together with their fluctuations. The aim of this analysis -- a comparison of experimental data with computational results obtained within frameworks of various hadron interaction models for protons and iron nuclei and an estimation of cosmic ray mass composition in the ultra-high energy region.

We present new data on Cherenkov light observations obtained during the period 1994-2009, after a modernization of the Yakutsk EAS array. A complex analysis of $\xmax$ and its fluctuations $\sigma(\xmax)$ was performed over a wide energy range. With the new data, according to QGSJet II model, an estimation was made of the cosmic ray mass composition for $\E \sim 10^{17} - 3 \times 10^{19}$ eV. The result points towards a mixed composition with a large portion of heavy nuclei at $\E \sim 10^{17}$ eV and the dominance of light nuclei at $\E \sim 10^{19}$ eV. The analysis of $\sigma(\xmax)$ energy dependence for the same energies qualitatively confirms this result. The shape of the $\xmax$ distribution at fixed energy $10^{18}$ eV is analysed to make more precise conclusions on cosmic ray mass composition.

The characteristics relating to the lateral and longitudinal development of EAS in the energy region of 10^15-10^19eV have been analyzed in the framework of the QGSJET model and of mass composition of primary cosmic rays. It is found that at E(0) >= 5*10^15eV the mean mass composition of primary cosmic rays begins to vary as indicated by a rise of <ln A> with increasing energy. The maximum value of <ln A> is observed at E(0) ~ (5-50)*10^16eV. It is confirmed by data of many compact EAS arrays and does not contradict an anomalous diffusion model of cosmic ray propagation in our Galaxy. In the superhigh energy region (>=10^18eV) the value <ln A> begins to decrease, i.e. the mass composition becomes lighter and consists of protons and nuclei of He and C. It does not contradict our earlier estimations for the mass composition and points to a growing role of the metagalactic component of cosmic rays in the superhigh energy region.

The differential energy spectrum of cosmic rays from Cherenkov radiation measurements in EAS in the energy range of 10^15-10^20eV has been compared with an anomalous diffusion model for the particles in interstellar space having fractal properties (Lagutin et al, 2001). The close association between experimental data and calculated "all particle" spectra in form at E(0) (10^15-10^18)eV is found. In this case, the average mass composition of cosmic rays calculated by five components does not contradict the average mass composition from experimental data which was obtained by several of EAS characteristics in that energy region.

The energy dependence of such characteristics as a ratio of the total number of charged particles to the total flux of EAS Cherenkov radiation, a ratio of E(thr)>=1GeV muon flux density at the distance of 600m from a shower core to charged particle flux density, a ratio of the energy transferred to the electromagnetic component of EAS to the primary particle energy is presented. Their comparison with two-component mass composition of cosmic rays (p-Fe) in the framework of calculations by a QGSJET model is given.

A shape of lateral distribution for charged particles in events with energy above 10^19eV is considered. Two methods were used for individual LDF parametrization. In the first approach, the index of power was determined for generalized Greisen-Linsley approximation. In second, mean square radius of the shower was determined for approximation proposed by Lagutin et al. Comparison of resulted parameters is presented for individual events arrived from different celestial regions -- Galactic planes and the region with increased flux of particles with E(0)>=10^19eV (according to Yakutsk array): 1.7h-3.7h right ascension; 45-60 degrees declination.

Experimental data obtained at the Yakutsk array after the modernization in 1993 are analyzed. The characteristics of EAS longitudinal and radial development found from the charged particle flux and EAS Cherenkov light registered at the Yakutsk complex array are presented. The energy spectrum of EAS obtained from Cherenkov light and an estimate of the PCR mass composition are presented.

The energy spectrum of primary cosmic rays with ultra high energies based on the Yakutsk EAS Array data is presented. For the largest events values of S600 and axis coordinates have been obtained using revised lateral distribution function. The affect of the arrival time distribution at several axis distance on estimated density for Yakutsk and AGASA is considered.

There are some discrepancies in the results on energy spectrum from Yakutsk, AGASA, and HiRes experiments. In this work differential energy spectrum of primary cosmic rays based on the Yakutsk EAS Array data is presented. For the largest events values of $S_{600}$ and axes coordinates have been obtained using revised lateral distribution function. Simulation of converter's response at large distances showed no considerable underestimation of particle density. Complex shape of spectrum in region of $\E > 10^{17}$ eV is confirmed. After adjustment of parameters and additional exposition at the Yakutsk array there are three events with energy $\E > 10^{20}$ eV.

The energy spectrum of UHECRs is discussed on the basis of the Yakutsk array database analysis. In the region E=0.1 to 30 EeV the showers are detected under trigger-500, while at energies above 30 EeV the whole acceptance area for trigger-1000 is used in order to utilize all the data available in the region of GZK cutoff.