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.

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

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

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.

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

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

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

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.