A method of simulation of particle cascades induced by ultra-high (>10$^{15}$ eV) and extremely high (>10$^{18}$ eV) energy neutrinos in water or other dense medium has been elaborated. The lateral spread of high-energy particles in cascades due to Coulomb scattering is negligible. So it is possible to use an approximation of the one-dimensional development of the UHE (or EHE) cascade if energies of particles in it higher than 10-100 GeV. An original program of the one-dimensional development of UHE and EHE cascades in dense media taking into account fluctuations and the Landau-Pomeranchuk-Migdal effect has been elaborated. The GEANT4 package has been used when particle energies are below 1000 GeV, for example, or less to calculate correctly the longitudinal distribution of an energy deposition by charge particles in a cascade. In advance the library of longitudinal characteristics of cascades was calculated. When a particle with the energy E which is below the threshold energy (E<1000 GeV or E<10$^6$ GeV) appears in the basic Monte Carlo program, the longitudinal distribution of an energy deposition induced by it was taken from the library and was added into the corresponding places of the basic distribution. The results of simulation of cascades in a dense medium can be used for calculations of acoustical or radio signals.

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

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

An attractive technique to explore for super-high-energy cosmic neutrino fluxes, via deep underwater acoustic detection, is discussed. Acoustic signals emitted by the neutrino induced cascades at large distances (10-50 km) from cascades are considered. It is argued that an existing hydroacoustic array of 2400 hydrophones, which is available in the Great Ocean near Kamchatka Peninsula, could be used as a base for an exploratory acoustic neutrino telescope SADCO (Sea Acoustic Detector of Cosmic Objects). The detection volume for registration of cascades with energies in the range of $10^{20-21} eV$ is estimated to be hundreds of cubic kilometers. Some models of extremely high energy elementary particle production in the Universe (for example the topological defect model) may be examined by such a detector. Tests of this technique are hoped for within a year.