We present numerical simulations on the propagation of UHE protons with energies of $(10^{19.5}-10^{22})$ eV in extragalactic magnetic fields over 1 Gpc. We use the ORS galaxy sample, which allow us to accurately quantify the contribution of nearby sources to the energy spectrum and the arrival distribution, as a source model. We calculate three observable quantities, cosmic ray spectrum, harmonic amplitude, and two point correlation function from our data of numerical simulations. With these quantities, we compare the results of our numerical calculations with the observation. We show that the three observable quantities including the GZK cutoff of the energy spectrum can be reproduced in the case that the number fraction $\sim 10^{-1.7}$ of the ORS galaxies more luminous than -20.5 mag is selected as UHECR sources. In terms of the source number density, this constraint corresponds to $10^{-6}$ Mpc$^{-3}$. However, since mean number of sources within the GZK sphere is only $\sim 0.5$ in this case, the AGASA 8 events above $10^{20.0}$ eV, which do not constitute the clustered events with each other, can not be reproduced. On the other hand, if the cosmic ray flux measured by the HiRes, which is consistent with the GZK cutoff, is correct and observational features about the arrival distribution of UHECRs are same as the AGASA, our source model can explain both the arrival distribution and the flux at the same time. Thus, we conclude that large fraction of the AGASA 8 events above $10^{20}$ eV might originate in the topdown scenarios, or that the cosmic ray flux measured by the HiRes experiment might be better. We also discuss the origin of UHECRs below $10^{20.0}$ eV through comparisons between the number density of astrophysical source candidates and our result ($\sim 10^{-6}$ Mpc$^{-3}$).

We have performed the detailed numerical simulations on the propagation of the UHE protons in the energy range $E=(10^{19.5} - 10^{22.0}$) eV in the relatively strong extra-galactic magnetic field with strength $B= (10, 100)$ nG within about 40 Mpc. In this case, the deflection angles of UHECRs become so large that the no counterparts problem is simply solved. As for the source distribution, we assumed that it is proportional to the number distribution of galaxies within the GZK sphere. We have found many clusters, which mean the small-scale anisotropy, in our simulations. It has been also shown that the observed energy spectrum is well reproduced in our models without any fine-tuned parameter. We have used the correlation value in order to investigate statistically the similarity between the distribution of arrival directions of UHECRs and that of galaxies. We have found that each correlation value for each parameter set begins to converge when the number of the detected events becomes ${\large O}(10^3)$. Since the expected number counts by the experiment of the next generation such as TA, HiRes, Auger, and EUSO are thought to be the order of $10^3$, we will be able to determine the source distribution and values of the parameters in this study in the very near future. Compared with the AGASA data, the significant anisotropy on the arrival directions of UHECRs are found in the analysis of first and second harmonics. This may originate from the incompleteness of the ORS database. This problem may be also solved if the source distribution is slightly changed. For example, this problem may be solved if we assume that UHECRs come from some of the galaxies such as AGNs and radio galaxies.