The paper presents the results of the longitudinal development of extensive air showers (X$_{max}$) of ultra-high energies and mass composition of cosmic rays. The measurements of X$_{max}$ are based on data from observations of the Cherenkov radiation at the Yakutsk array for the period 1974-2014. The cascade curves of individual showers and the depth of maximum X$_{max}$ were reconstructed over the energy range 10$^{16}$-5.7$\cdot$10$^{19}$ eV. It is shown that the displacement rate of the parameter dX$_{max}$ / dE in the atmosphere is nonlinear and depends on the energy. Such a feature indicates a change in mass composition, which is confirmed by fluctuations of X$_{max}$ in this energy region. The composition of cosmic rays was determined by interpolation using the QGSJetII-04 model.

The Yakutsk Array is designed to study cosmic rays at energy 10$^{15}$-10$^{20}$ eV. It consists several independent arrays that register charged particles, muons with energy E$\geq$1 GeV, Cherenkov light and radio emission. The paper presents a technical description of the Yakutsk Radio Array and some preliminary results obtained from measurements of radio emission at 30-35 MHz frequency induced by air shower particles with energy $\varepsilon$ $\geq$ 1$\cdot$10$^{17}$ eV. The data obtained at the Yakutsk array in 1986-1989 (first set of measurements) and 2009-2014 (new set of measurements). Based on the obtained results we determined: Lateral distribution function (LDF) of air showers radio emission with energy $\geq$ 10$^{17}$ eV. Radio emission amplitude empirical connection with air shower energy. Determination of depth of maximum by the ratio of amplitude at different distances from the shower axis. For the first time, at the Yakutsk array, radio emission from the air shower with energy $>$ 10$^{19}$ eV was registered including the shower with the highest energy ever registered at the Yakutsk array with energy $\sim$2 $\cdot$ 10$^{20}$ eV.

Long-term series of measurement of spectral transparency of the atmosphere (\lambda = 430 nm) and atmospheric optical thickness (AOT) measured by multimode photometer CE 318 in the region of Yakutsk array are analyzed. Correlation of AOT with intensity of air showers with small energies 10^{15} - 10^{16} eV is found. The variability of aerosol composition of the atmosphere during the registration period of the Cherenkov light should be taken into account since it may affect the quality of determining characteristics of air showers.

This paper presents data on the optical thickness AOT in the context of long-term observations, including in times of large-scale forest fires in Yakutia. Sudden changes in such features as the AOT point directly to the impact of fires on aerosol structure of the atmosphere during the summer.

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.

In this paper, we present results obtained from the measurements of radio emission at frequency of 32 MHz with energy more than 10$^{19}$ eV. Generalized formula that describe lateral distribution and depends on main characteristic of the air showers: energy E$_0$ and depth of maximum X$_{max}$ was derived. The formula has a good agreement with data at average and large distances from shower axis. Employing the ratio of radio emission amplitude at distances 175 m and 725 m we determined the depth of maximum X$_{max}$ for air shower with energy 3.7$\cdot$10$^{19}$ eV, which in our case is equal to X$_{max}$ = 769$\pm$34g$\cdot$cm$^{-2}$.

The paper presents short technical description of Yakutsk Radio Array and some preliminary results obtained from measurements of radio emission at 32 MHz frequency induced by air shower particles with energy {\epsilon}$\geq$1$\cdot$10$^{19}$ eV. The data obtained at Yakutsk array in 1987-1989 (first set of measurements) and 2009-2014 (new set of measurements). For the first time, at Yakutsk array radio emission from air shower with energy > 10$^{19}$ eV was registered including the shower with highest energy ever registered at Yakutsk array with energy $\sim$2$\cdot$10$^{20}$ eV.

The paper presents long-term observations of the atmosphere in Yakutsk region. Analysis of the data for 40 year period indicates a gradual strengthening of cyclonic activity in the region and hence the increase of the average winter temperature, increase variations of the rest atmosphere, which greatly softens the continental climate of Central Yakutia.

Calibration of a PMT matrix is crucial for the treatment of the data obtained with Cherenkov tracking detector. Furthermore, due to high variability of the aerosol abundance in the atmosphere depending on season, weather etc. A constant monitoring of the atmospheric transparency is required during the measurements. For this purpose, besides traditional methods, a station for laser atmospheric probing is used.

The paper present correlation of radio signal with air shower parameters: shower energy E0 and depth of maximum Xmax. It is shown that from radio emission measurements of air showers one can obtain individual showers parameters and mass composition of cosmic rays. We also derived generalized formula for calculating energy of the air showers.

The Yakutsk array group is developing a wide field-of-view Cherenkov telescope to be operated in coincidence with the surface detectors of the extensive air shower array. Currently, the engineering prototype of the reflecting telescope with the front-end electronics is designed, assembled, and tested to demonstrate the feasibility of the conceived instrument. The status and specifications of the prototype telescope are presented, as well as the modernization program of the already existing Cherenkov light detectors subset of the array measuring ultra-high energy cosmic rays.

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.

From the data on Cherenkov light from EASs collected over the period of more than 15 years, the spectrum of cosmic rays was obtained in the energy range 1.0E15-1.0E18 eV. This spectrum has two features: at ~3.0E15 eV (first knee) and at ~1.0E17 eV (second knee). The first knee is characterized by the indexes \gamma_1=2.7+/-0.03 and \gamma_2 = 3.12+/-0.03 and the second knee - by \gamma_1 = 2.92+/-0.03 and \gamma_2=3.24+/-0.04. In the second case the difference amounts to \Delta\gamma_23=0.32{\pm}0.03{\pm}0.05 which is less than the case of the first knee \Delta\gamma_12=0.42+/-0.03+/-0.05 (here dual errors represent statistics and systematics correspondingly). The lesser difference \Delta\gamma_23 can be explained with the influx of cosmic rays from Meta-galaxy and, hence, with some increase of the cosmic rays intensity in the energy range 5.0E16-1.0E18 eV, which compensates the escape of heavier nuclei from the Galaxy. The presence of the second knee could be confirmed by abrupt change in mass composition from from lnA~3 at ~1.0E17 eV to lnA~1.5 at 1.0E18eV found from the analysis of longitudinal development of EAS.

The paper describes the techniques and method of registration of air shower radio emission at the Yakutsk array of extensive air showers at a frequency of 32 MHz. At this stage, emission registration involves two set of antennas, the distance between them is 500m. One set involves 8 antennas, second - 4 antennas. The antennas are perpendicularly crossed dipoles with radiation pattern North South,West East and raised 1.5 m above the ground. Each set of antennas connected to an industrial PC. The registration requires one of two triggers. First trigger are generated by scintillation detectors of Yakutsk array. Scintillation detectors cover area of 12 km^2 and registers air showers with energy more than 10^17 eV. The second trigger is generated by Small Cherenkov Array that covers area of 1 km^2 and registers air showers with energy 10^15 - 5*10^17 eV. Small Cherenkov Array is part of Yakutsk array and involve Cherenkov detectors located at a distance of 50, 100, 250 m. For further selection we are using an additional criterion the radio pulse must be localized in the area corresponding to the delay time on first and second triggers. In addition, descriptions of the algorithm and the flowcharts of the program for the air shower selection and further analysis are given. This method registers EAS radio emission with energy 10^16 - 10^19 eV. With the absolute calibration, the amplitudes of all antennas converted to a single value. Air shower radio emission dependences from zenith angle and shower energy are plotted.

We present of measurements of radio emission from EAS at the frequency 32 MHz performed in 2008-2012. Showers were selected by zenith and azimuth directions, sampled by energy. A mean lateral distribution function was fitted to the data each interval. The paper presents unied approximation formula for all three energy intervals of lateral distribution, using connection between equation (1) and (2) coefficients with energy and geomagnetic angle. Here we analyze the connection between the shape of lateral distribution and the depth of the maximum shower development Xmax too.

This paper presents the set of measurements of ultra-high energy air shower radio emission at frequency 32 MHz in period of 2008-2012. The showers are selected by geomagnetic and azimuth angles and then by the energy in three intervals: 3*10^16 3*10^17 eV, 3*10^17 6*10^17 eV and 6*10^17 5*10^18 eV. In each energy interval average lateral distribution function using mathematically averaged data from antennas with different directions are plotted. In the paper, using experimental data the dependence of radio signal averaged amplitude from geomagnetic angle, the shower axis distance and the energy are determined. Depth of maximum of cosmic ray showers Xmax for the given energy range is evaluated. The evaluation is made according QGSJET model calculations and average lateral distribution function shape.

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 a summary of the measurements of mass sensitive parameters at the highest cosmic ray energies done by several experiments. The Xmax distribution as a function of energy has been measured with fluorescence telescopes by the HiRes, TA and Auger experiments and with Cherenkov light detectors by Yakutsk. The <Xmax> or the average mass (<lnA>) has been also inferred using ground detectors, such as muon and water Cherenkov detectors. We discuss the different data analyses elaborated by each collaboration in order to extract the relevant information. Special attention is given to the different approaches used in the analysis of the data measured by fluorescence detectors in order to take into account detector biases. We present a careful analysis of the stability and performance of each analysis. The results of the different experiments will be compared and the discrepancies or agreements will be quantified.

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.

The aim of the Yakutsk array enhancement project is to create an instrument to study the highest-energy galactic cosmic rays (CRs) -- their sources, energy spectrum, and mass composition. Additionally, there will be unique capabilities for investigations in the transition region between galactic and extragalactic components of CRs. Using the well-developed imaging atmospheric Cherenkov telescope technique adapted to the energy region $E>10^{15}$ eV, we plan to measure the longitudinal structure parameters of the shower, e.g., angular and temporal distributions of the Cherenkov signal related to $X_{max}$ and the mass composition of CRs. The main advantages of the Yakutsk array, such as its multi-component measurements of extensive air showers, and model-independent CR energy estimation based on Cherenkov light measurements, will be inherited by the instrument to be created.

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.

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.

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 energy spectrum of cosmic rays in the range 10^15 eV to 6*10^19 eV has been studied using the air Cherenkov light detectors of the Yakutsk array. The total flux of photons produced by relativistic electrons (including positrons as well, hereafter) of extensive air showers in the atmosphere is used as the energy estimator of the primary particle initiating a shower. The resultant differential flux of cosmic rays exhibits, in accordance with previous measurements, a knee and ankle features at energies 3*10^15 and ~10^19 eV, respectively. A comparison of observational data with simulations is made in the knee and ankle regions in order to choose the models of galactic and extragalactic components of cosmic rays which describe better the energy spectrum measured.

Signal from scintillation detectors in showers with energies 10**17-10**20eV has been analyzed. It was determined, that for core distances >=200m the time scale of a particle arrival time at the detector constitutes 40-5000ns and longer. Events with fine structure of the signal have been observed. The pulse consists of several peaks (for vertical showers) and a single pulse (for inclined showers). There are some anomalous events with two pulses wherein such a pulse shape is presented in two, three and more detectors. A period between these pulses makes 60-350ns.

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.

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.

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.

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.

At the Yakutsk EAS array E_em is determined by using measurements of EAS Cherenkov light flux and charged particle flux. It is known from calculations that these characteristics depend on a sort of primary particle and, therefore, the estimation of E_em depends on a primary particle mass. In the work the dependence of the E_em/E_0 ratio on the energy is given and experimental data are compared with calculations by the QGSJET model. The calculations have been carried out for the primary proton and iron nucleus. The average calculated meaning of the value of E_em/E_0 ratio (between the proton and iron nucleus) within experimental errors is in agreement with experimental data that doesnt contradict to the mixed mass composition of primary cosmic radiation.

The ratio of the muon flux density to charged particle flux density at distances of 300 and 600 m from the shower axis ($\rhom(300)/\rhos(300)$ and $\rhom(600)/\rhos(600)$) is measured. In addition, the energy dependence of $\rhom(1000)$ is analysed for showers with energies above $10^{18}$ eV. A comparison between the experimental data and calculations performed with the QGSJET model is given for the cases of primary proton, iron nucleus and gamma- ray. We conclude that the showers with $\E\ge3\times10^{18}$ eV can be formed by light nuclei with a pronounced fraction of protons and helium nuclei. It is not excluded however that a small part of showers with energies above $10^{19}$ eV could be initiated by primary gamma-rays.

The experimental distributions of $\Xmax$ obtained with the Yakutsk EAS array at fixed energies of $5 \times10^{17}$, $1\times10^{18}$ and $5\times10^{18}$ eV are analysed. A recent version of the QGSJET model is used as a tool of our analysis. In the framework of this model, the most adequate mass composition of primary particles satisfying the experimental data on $\Xmax$ is selected.

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.

In the work the data of the Yakutsk complex EAS array and their comparison with calculation in the case of primary nuclei of different chemical elements are presented. The calculation by QGSJET model have been used interpreting experimental data.

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.

The Yakutsk array observes the Cherenkov light emitted by UHECR in atmosphere. Recently, an autonomous subarray is added consisting of photomultipliers to measure the showers in the knee region. Our aim is to analyze the combined data set in order to derive the cosmic ray spectrum in the energy range as wide as possible using the same technique. The advantage of the air Cherenkov light measurement is the model independent estimation of the EAS primary energy using the total light flux emitted in the atmosphere. A set of the light lateral distributions observed in the extended energy range is presented also.