Two EAS arrays during a day have recorded 3 particles with energies above 30 EeV arriving from the same sky region. Two events were registered by the Yakutsk array and one - by the Telescope Array. Two Yakutsk events were estimated to be same energy and the Telescope Array shower's energy was almost 2 times higher. This indicates the same magnetic rigidity of all three particles, if the charge differs by a factor of 2. The relatively short sequence of all three events and their "monochromaticity" in rigidity can be due to the magnetic separation of particles in the acceleration process and propagation of the beam.

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.

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

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.

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

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.

Comparing the signals measured by the surface and underground scintillator detectors of the Yakutsk Extensive Air Shower Array, we place upper limits on the integral flux and the fraction of primary cosmic-ray photons with energies E > 10^18 eV, E > 2*10^18 eV and E > 4*10^18 eV. The large collected statistics of the showers measured by large-area muon detectors provides a sensitivity to photon fractions < 10^-2, thus achieving precision previously unreachable at ultra-high energies.

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.

We analyse a sample of 33 extensive air showers (EAS) with estimated primary energies above 2\cdot 10^{19} eV and high-quality muon data recorded by the Yakutsk EAS array. We compare, event-by-event, the observed muon density to that expected from CORSIKA simulations for primary protons and iron, using SIBYLL and EPOS hadronic interaction models. The study suggests the presence of two distinct hadronic components, ``light'' and ``heavy''. Simulations with EPOS are in a good agreement with the expected composition in which the light component corresponds to protons and the heavy component to iron-like nuclei. With SYBILL, simulated muon densities for iron primaries are a factor of \sim 1.5 less than those observed for the heavy component, for the same electromagnetic signal. Assuming two-component proton-iron composition and the EPOS model, the fraction of protons with energies E>10^{19} eV is 0.52^{+0.19}_{-0.20} at 95% confidence level.

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.

By making use of the data on the total signal and on the muon component of the air showers detected by the Yakutsk array, we analyze, in the frameworks of the recently suggested event-by-event approach, how large the fraction of primary gamma-rays at ultra-high energies can be. We derive upper limits on the photon fraction in the integral flux of primary cosmic rays. At the 95% confidence level (CL), these limits are 22% for primary energies E_0>4\cdot 10^{19}eV and 12% for E_0>2\cdot 10^{19}eV. Despite the presence of muonless events, the data are consistent with complete absence of photons at least at 95% CL. The sensitivity of the results to systematic uncertainties, in particular to those of the energy determination for non-photon primaries, is discussed.

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.

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.

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.

The Earth's magnetic field effect on the azimuthal distribution of extensive air showers (EAS) of cosmic rays has been evaluated using a bulk of the Yakutsk array data. The uniform azimuthal distribution of the EAS event rate is rejected at the significance level 10^(-14). Amplitude of the first harmonics of observed distribution depends on zenith angle as A1=0.2*sin^2(theta) and is almost independent of the primary energy; the phase coincides with the magnetic meridian. Basing upon the value of measured effect, the correction factor has been derived for the particle density depending on a geomagnetic parameter of a shower.