Astrometric observations performed by the Gaia Follow-Up Network for Solar System Objects (Gaia-FUN-SSO) play a key role in ensuring that moving objects first detected by ESA's Gaia mission remain recoverable after their discovery. An observation campaign on the potentially hazardous asteroid (99 942) Apophis was conducted during the asteroid's latest period of visibility, from 12/21/2012 to 5/2/2013, to test the coordination and evaluate the overall performance of the Gaia-FUN-SSO . The 2732 high quality astrometric observations acquired during the Gaia-FUN-SSO campaign were reduced with the Platform for Reduction of Astronomical Images Automatically (PRAIA), using the USNO CCD Astrograph Catalogue 4 (UCAC4) as a reference. The astrometric reduction process and the precision of the newly obtained measurements are discussed. We compare the residuals of astrometric observations that we obtained using this reduction process to data sets that were individually reduced by observers and accepted by the Minor Planet Center. We obtained 2103 previously unpublished astrometric positions and provide these to the scientific community. Using these data we show that our reduction of this astrometric campaign with a reliable stellar catalog substantially improves the quality of the astrometric results. We present evidence that the new data will help to reduce the orbit uncertainty of Apophis during its close approach in 2029. We show that uncertainties due to geolocations of observing stations, as well as rounding of astrometric data can introduce an unnecessary degradation in the quality of the resulting astrometric positions. Finally, we discuss the impact of our campaign reduction on the recovery process of newly discovered asteroids.

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.

We report a measurement of the neutrino-electron elastic scattering rate from 8B solar neutrinos based on a 123 kton-day exposure of KamLAND. The background-subtracted electron recoil rate, above a 5.5 MeV analysis threshold is 1.49+/-0.14(stat)+/-0.17(syst) events per kton-day. Interpreted as due to a pure electron flavor flux with a 8B neutrino spectrum, this corresponds to a spectrum integrated flux of 2.77+/-0.26(stat)+/-0.32(syst) x 10^6 cm^-2s^-1. The analysis threshold is driven by 208Tl present in the liquid scintillator, and the main source of systematic uncertainty is due to background from cosmogenic 11Be. The measured rate is consistent with existing measurements and with Standard Solar Model predictions which include matter enhanced neutrino oscillation.

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.

We present lightcurve observations and multiband photometry for 107P/Wilson-Harrington using five small- and medium-sized telescopes. The lightcurve has shown a periodicity of 0.2979 day (7.15 hour) and 0.0993 day (2.38 hour), which has a commensurability of 3:1. The physical properties of the lightcurve indicate two models: (1) 107P/Wilson-Harrington is a tumbling object with a sidereal rotation period of 0.2979 day and a precession period of 0.0993 day. The shape has a long axis mode (LAM) of L1:L2:L3 = 1.0:1.0:1.6. The direction of the total rotational angular momentum is around {\lambda} = 310{\deg}, {\beta} = -10{\deg}, or {\lambda} = 132{\deg}, {\beta} = -17{\deg}. The nutation angle is approximately constant at 65{\deg}. (2) 107P/Wilson-Harrington is not a tumbler. The sidereal rotation period is 0.2979 day. The shape is nearly spherical but slightly hexagonal with a short axis mode (SAM) of L1:L2:L3 = 1.5:1.5:1.0. The pole orientation is around {\lambda} = 330{\deg}, {\beta} = -27{\deg}. In addition, the model includes the possibility of binary hosting. For both models, the sense of rotation is retrograde. Furthermore, multiband photometry indicates that the taxonomy class of 107P/Wilson-Harrington is C-type. No clear rotational color variations are confirmed on the surface.

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.

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.

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.