It has been shown that muon flux intensities calculated in terms of the EPOS LHC and EPOS 1.99 models at the energy of 10^4 GeV exceed the data of the classical experiments L3+Cosmic, MACRO and LVD on the spectra of atmospheric muons by a factor of 1.9 and below these data at the same energy by a factor of 1.8 in case of the QGSJET II-03 model. It has been concluded that these tested models overestimate (underestimate in case of QGSJET II-03 model) the production of secondary particles with the highest energies in interactions of hadrons by a factor of ~1.5. The LHCf and TOTEM accelerator experiments show also this type of disagreements with these model predictions at highest energies of secondary particles.

A simple method of the vertical muon energy spectrum simulations have been suggested. These calculations have been carried out in terms of various models of hadronic interactions. The most energetic $ \pi^\pm $-mesons and K$^\pm $-mesons produced in hadron interactions contribute mainly in to this energy spectrum of muons due to the very steep energy spectrum of the primary particles. So, some constraints on the hadronic interaction models may be set from a comparison of calculated results with the cosmic data on the vertical muon energy spectrum. This comparison showed that the most energetic secondary particles production is too high in case of the QGSJET II-04 model and rather low in case of the QGSJET II-03 model. These conclusion have been supported by the LHC data.

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.

Here the results of calculations of pulses in Cherenkov light detectors for the Yakutsk array are presented. As long as the Vavilov-Cherenkov light is used to calibrate signals in scintillation detectors at the Yakutsk array it is vital for these measurements to be precise. The validity of measurements of the signal Q(400) used as the estimation parameter at the Yakutsk array has been confirmed. Our calculations show that the width of time pulses increases from nearly 50 ns at a distance of 100 m from the shower axis up to 700 ns at 1000 m.

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 primary cosmic rays particles with energies above 10**20 eV have been observed at many extensive air shower arrays since the beginning of observations over 40 years ago. The validity of measurements of signal s(600) used as energy estimation parameter at the Yakutsk array has been confirmed. Our calculations show that the width of the time pulses increases from nearly 100 ns at a distance of 100 m from the shower axis up to 4 - 5 $\mu$s at 1500 m. The calculated estimate of energy of extensive air shower is ~ 1.7 times smaller than the experimental estimate for the same value s(600). The pointing directions of extensive air showers observed at the Pierre Auger Observatory were fitted within +-3.1**o with positions of the nearby active galactic nuclei from the Veron-Cetty and P. Veron catalog. The cosmic ray luminosity of the active galactic nuclei which happened to be a source of the particular cosmic ray event constitutes a fraction ~10**-4 of the optical one if only cosmic ray particles with energies above 6*10**19 eV are produced. If produced cosmic ray particles have a spectrum up to ~ 100 GeV then the cosmic ray luminosity of the active galactic nuclei should be much higher than the optical one.

The pointing directions of extensive air showers observed at the Pierre Auger Observatory were fitted within 3.1 degree with positions of the nearby active galactic nuclei from the Veron-Cetty and P. Veron catalog. The cosmic ray luminosity of the active galactic nuclei which happened to be a source of the particular cosmic ray event constitutes a fraction ~0.0001 of the optical one if only cosmic ray particles with energies above 60 EeV are produced. If produced cosmic ray particles have a spectrum dE/E^3 up to ~100 GeV then the cosmic ray luminosity would be much higher than the optical one of the active galactic nuclei.

In terms of the quark-gluon string model the analysis of the classic procedure to estimate the energy of giant air showers with help of the parameter s(600) (a density of energy deposition in the scintillator at a distance of 600 m from the shower core) have been carried out. Simulations of the signal s(600) with help of the CORSIKA code in terms of the hybrid scheme show energy estimates which are approximately a factor of 1.6 times lower than adopted at the Yakutsk array. The energy estimates calculated with the help of the Cherenkov radiation coincide with the experimental data. Simulations of deposited energy distributions in the atmosphere with help of the GEANT4 code and the CORSIKA code show that more than 20% of this energy may be deposited at distances above 100 m from the shower axis.

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.