Arrival directions of ultrahigh energy extensive air showers (EAS) by world arrays data are analyzed. It is shown that arrival directions of ultrahigh energy EAS by Yakutsk and P. Auger data are some correlated with pulsars of Local Arm Orion of Galaxy. It is found large amplitude of first harmony of Fourier row by data Yakutsk array. The problem of cosmic rays origin is discussed.

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.

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.