Lateral distributions of electrons in air showers initiated by photons of ultra high energies ($10^{18}-10^{22}$ eV) obtained on the basis of numerical solution of adjoint cascade equations are presented. An extended analysis has been made considering separately the Landau-Pomeranchuk-Migdal (LPM) effect and the interaction of photons and electrons with the geomagnetic field (GMF) with respect to the scaling formalism for lateral distributions. It is shown that one-parametric scaling description of the lateral distribution of electrons remains valid up to the highest energies considering the LPM and GMF effects, that allow effective primary particle type discrimination using the surface detectors data of largest ground-based air shower arrays.

Lipari and Stanev developed a method for range fluctuation of high energy muons, stressing the importance of accounting for the fluctuations of the energy loss in radiative processes in 1991 and, now, their method has become the basement for the energy determination of high energy muons through the measurement of the Cherenkov light yields due to those muons in KM3 physics. Once, Takahashi et al. developed a method for the investigation on the depth intensity relation of high energy muons in which all the stochastic processes concerned are taken into account exactly (1983). Now, we make the method by Takahashi et al. revival for the same purpose of the application to the analysis of future KM3 physics. In the present paper, our concern is restricted to the introduction to the fundamental of our method and some subsequent results thereby in which the real simulated behaviors of high energy muons from $10^{12}$eV to $10^{18}$eV, the survival probabilities of high energy and so on are included. The discussion around the practical application of our method to the KM3 physics is entrusted in the subsequent papers. As far as the survival probability of high energy muons is concerned, our method gives nearly the same results to Lipari and Stanev's in some regions and gives the deviated results from theirs in another ones. Thus, we examine the application limit of their method and clarify the reason why, comparing with our method. The most distinct difference in the both methods may become apparent in the treatment on the Cherenkov light yields spectrum by which one may estimate the energies of the muons concerned. We will mention to them in subsequent papers.

In the present paper we consider neutrino events due to quasi-elastic scattering (QEL) as the most reliable events among various candidate events to be analyzed, and have carried out the first step of an L/E analysis which aims to confirm the survival probability with a Numerical Computer Experiment. The most important factor in the survival probability is Lnu and Enu, but this cannot be measured for such neutral particles. Instead, Lmu and Emu is utilized in the L/E analysis, where Lnu, Lmu, Enu and Emu denote the flight path lengths of the incident neutrinos, those of the emitted leptons, the energies of the incident neutrinos and those of the emitted leptons, respectively. According to our Computer Numerical Experiment, the relation of Lnu/Enu is nearly equal to Lmu/Emu doesn't hold. In subsequent papers, we show the results on an L/E analysis with the Computer Numerical Experiment based on our results obtained in the present paper.

The validity of our Monte Carlo simulation procedure (the time sequential procedure) had been verified by the corresponding analytical procedure which is independent of our method methodologically. Also, the results obtained by our procedure are compared with those obtained by the different Monte Carlo simulation procedure (the v_{cut} procedure) which has been exclusively utilized by different authors and the agreement between ours and theirs are found to be good. By utilizing our Monte Carlo procedure, the validity of which is guaranteed in two different procedures, we investigate not only the fluctuation of high energy muons themselves but also the fluctuation of the various quantities related to the energy losses by the muons, which are difficult to obtain by the v_{cut} procedure. Namely, we obtain the fluctuation on energy losses of muons in the present paper(Part1). In a subsequent paper (Part2), based on the idea developed in Part1, we will examine the correlation between the primary energies of muons and their energy deposits and the correlation between the energy deposits of the muons concerned and their Cherenkov light yields in the KM3 detector in detail.

The validity of our Monte Carlo simulation procedure ({\it the integral method}) had been verified by the corresponding analytical procedure of which is quite independent of our method methodologically. Also, the results obtained by our procedure are compared with those obtained by the different Monte Carlo simulation procedure ({\it the differential method}) which have been exclusively utilized by the different authors and the agreement between them are found to be well. By utilizing our Monte Carlo procedures, the validity of which is guaranteed in two different procedures, we investigate not only the fluctuation of high energy muons themselves but also fluctuation of the various quantities related to the energy losses by the muons, which are difficult to obtain by {\it the differential method}. Namely, we obtain fluctuation on energy losses of the muons, fluctuation on Cherenkov lights due to the accompanied cascade showers initiated by the muon and the correlations between them . Finally, we obtain the transition curves for Cherenkov lights in KM3 detector, taking into account of all possible fluctuations in the stochastic processes and point out the difficulty of the reliable estimation of the energy of the muons which are resultants of muon neutrino events in the KM3 detectors.

We calculate high and (ultra-)high energy upward-going muon neutrino propagation through the Earth and the induced muon energy distribution near the one cubic kilometer detector using the Monte Carlo simulation, due to both charged current interaction and neutral one. The initiated neutrino energies on the surface of the Earth are 1PeV, 1EeV and 1ZeV. The mean free paths of (ultra-)high energy neutrino events generated by the deep inelastic scattering may be comparable with the diameter of the Earth or less than it. Therefore, the induced muon production distribution is influenced by the change of the densities interior to the Earth. Furthermore, in such situation, the contribution from the neutral current neutrino interaction to the induced muon production distribution cannot be neglected. We report several examples of the deep inelastic scattered depth of muon neutrino in the Earth and the induced muon energy distribution near the detector.

By using the integral method in the muon propagation through water, we calculate the range fluctuation of high and ultra high energy muons. Many authors divide all radiative processes into two parts, namely, the continuous part and radiative part in their Monte Carlo simulation in order to consider the fluctuation in the both ranges and energies of the muons, while we treat all stochastic processes as exactly as possible, without the introduction of the continuous parts in all stochastic processes. The validity of our Monte Carlo method is checked by the corresponding analytical method which is methodologically independent on the Monte Carlo procedure. Accompanied cascade showers are generated by the direct electron pair production, bremsstrahlung and photo-nuclear interaction. These showers are calculated by the exact Monte Carlo Method in one dimensional way. We report survival probabilities, range distributions and examples of individual muon behavior.

In the near future the energy region above few hundreds of TeV may really be accessible for measurements of the atmospheric muon spectrum by the IceCube array. Therefore one expects that muon flux uncertainties above 50 TeV, related to a poor knowledge of charm production cross sections and insufficiently examined primary spectra and composition, will be diminished. We give predictions for the very high-energy muon spectrum at sea level, obtained with the three hadronic interaction models, taking into account also the muon contribution due to decays of the charmed hadrons.

In the previous paper (Part~1), we have verified that the SK assumption on the direction does not hold in the analysis of neutrino events occurred inside the SK detector. We have made four possible L/E analyses, L_nu/E_nu, L_nu/E_mu, L_mu/E_nu and L_mu/E_mu. Among four kinds of L/E analyses, we have shown that only L_nu/E_nu analysis can give the signature of maximum oscillations clearly, while the L_mu/E_mu analysis which are really done by Super-Kamiokande Collaboration cannot give the maximum oscillation at all. It is thus concluded that Super-Kamiokande type experiment cannot find the maximum oscillation from L/E analysis. Therefore, we would suggest Super-Kamiokande Collaboration to re-analyze the zenith angle distribution of the neutrino events which occur inside the detector carefully.

The energy spectra of hadron cascade showers produced by the cosmic ray muons travelling through water as well as the muon energy spectra underwater at the depth up to 4 km are calculated with two models of muon inelastic scattering on nuclei, the recent hybrid model (two-component, 2C) and the well-known generalized ector-meson-dominance model for the comparison. The 2C model involves photonuclear interactions at low and moderate virtualities as well as the hard scattering including the weak neutral current processes. For the muon scattering off nuclei substantial uclear effects, shadowing, nuclear binding and Fermi motion of nucleons are taken into account. It is shown that deep nderwater muon energy spectrum calculated with the 2C model are noticeably distorted at energies above 100 TeV as compared to that obtained with the GVMD model.

Extragalactic background of high energy neutrinos arising from the interactions of cosmic ray protons with far-infrared extragalactic background radiation is calculated. The main assumption is that the cosmic ray spectrum at energies higher than $10^8 GeV$ has extragalactic origin and, therefore, is proton dominated. All calculations are performed with taking into account the possible cosmological evolution of extragalactic sources of cosmic ray protons as well as infrared-luminous galaxies.

The SuperKamiokande group assert that they have found an oscillatory signature in atmospheric neutrinos through the analysis of \textit{Fully Contained Events} and \textit{Partially Contained Events}. We have performed an $L/E$ (length/energy) analysis of \textit{Upward Through-Going Muon Events} and \textit{Stopping Muon Events} in a numerical computer simulations both with and without neutrino oscillations but were unable to find an oscillatory signature. We give likely explanations for the absence of the oscillatory signature in our simulations and its apparent presence in the SuperKamiokande data.

Quasi Elastic Scattering (QEL) is the dominant mechanism for producing both Fully Contained Events and Partially Contained Events in the SuperKamiokande (SK) detector for atmospheric neutrinos in the range 0.1GeV to 10GeV. In the analysis of SK events, it is assumed that the zenith angle of the incident neutrino is the same as that of the detected charged lepton. In the present paper, we derive the distribution function for the scattering angle of charged leptons and derive the average scattering angles as well as their standard deviation due to QEL. Thus, it is shown that the SK assumption for the scattering angle of the charged lepton in QEL is not valid. Further, for a real analysis of the experimental data, we should examine the correlations between the zenith angle of the charged leptons and those of their incident neutrinos, taking into account the effect of the azimuthal angle of the charged particle over its zenith angle. We conclude that it is not possible to reliably assign the zenith angle of the incident neutrino to that of the charged lepton, and so the zenith angle distribution of charged leptons does not reflect that of the incident neutrinos under the SK assumption. This result has clear implication for attempts to detect neutrino oscillations from an analysis of Fully Contained Events and Partially Contained Events in SuperKamiokande.

Quasi Elastic Scattering (QEL) is the dominant source for producing both Fully Contained Events and Partially Contained Events in the Superkamiokande(SK) detector for the atmospheric neutrinos, in the range 0.1 GeV to 10 GeV. In the analysis of SK events, it is assumed that the direction of the incident neutrino is the same as that of the detected charged lepton. In the present letter, we derive the distribution function for the scattering angle of the charged leptons, their averaged scattering angle and their standard deviation due to QEL. Then, it is shown that the SK assumption for the scattering angle of the charged leptons in the QEL is not valid. Further, we examine the influence of the azimuthal angle of the charged leptons over their zenith angle. As the result, we conclude that the zenith angle distribution of the neutrino under the SK assumption does not reflect the real zenith angle distribution of the atmospheric neutrino which produces Fully Contained Events and Partially Contained Events. This result has clear implication for attempts to detect neutrino oscillations from the analyses of Fully Contained Events and Partially Contained Events in Superkamiokande.

We discuss the problem of cosmic ray mass composition variation in wide energy region 10^{14}-10^{20} eV. The mass composition predicted in the framework of recently developed anomalous diffusion model is tested using the results of CORSIKA calculations and experimental data on the depth of maximum of extensive air showers. We show that the model predictions for the mass composition are consistent with observations of different experiments.

Detailed theoretical study of the mean square radius of extensive air shower electrons has been made in connection with further development of scaling formalism for electron lateral distribution function. A very simple approximation formula, which allows joint description of all our results obtained in wide primary energy range and for different observation depths is presented. The sensitivity of the mean square radius to variations of basic parameters of hadronic interaction model is discussed.

Theoretical predictions for lateral distribution function of electrons in extensive air showers based on scaling formalism are presented. Our results are tested by comparison with AGASA experimental data taking into account the contribution of low-energy muons and simultaneously the effect of scintillation detectors response, according to recent simulation results. The possibility for ultrahigh energy cosmic ray primary composition deduction from the shape of LDF is discussed in detail.

The vertical sea-level muon spectrum at energies above 1 GeV and the underground/underwater muon intensities at depths up to 18 km w.e. are calculated. The results are particularly collated with a great body of the ground-level, underground, and underwater muon data. In the hadron-cascade calculations, the growth with energy of inelastic cross sections and pion, kaon, and nucleon generation in pion-nucleus collisions are taken into account. For evaluating the prompt muon contribution to the muon flux, we apply two phenomenological approaches to the charm production problem: the recombination quark-parton model and the quark-gluon string model. To solve the muon transport equation at large depths of homogeneous medium, a semi-analytical method is used. The simple fitting formulas describing our numerical results are given. Our analysis shows that, at depths up to 6-7 km w. e., essentially all underground data on the muon intensity correlate with each other and with predicted depth-intensity relation for conventional muons to within 10%. However, the high-energy sea-level data as well as the data at large depths are contradictory and cannot be quantitatively decribed by a single nuclear-cascade model.