Probabilistic seismic hazard analysis (PSHA) is needed to understand the regions where strong ground motions occur and its occurrence probability. In PSHA for the specific area, spatial correlation of seismic ground motion should be considered. Meanwhile, the return period of some large earthquake around Japan are estimated as a few hundreds of years, and the probability of occurrence rise significantly over several decades with the passage of time. In areas affected by such earthquakes, it is necessary to consider the temporal changes of occurrence probability during the period to be considered in PSHA. In this study, we proposed the spatial probabilistic seismic hazard analysis considering the spatial correlation of earthquake ground motion with the temporal change of the occurrence probabilities of large earthquakes. We modeled spatial correlation of seismic ground motion from observation data by statistical processing and calculate the occurrence probability of a large earthquake that changes with time taking into consideration the uncertainty of data. We calculated the relationship between the maximum seismic intensity and probability of occurrence in a specific area by trial calculation for Kanagawa prefecture and discussed its temporal change. Furthermore, application of such analysis results to regional disaster prevention planning was discussed.