Probabilistic seismic hazard analysis (PSHA) plays a crucial role in providing critical information about quake-prone regions to design and practicing engineers. In addition, PSHA also helps in identifying possible locations for nuclear sites and deployment of future seismic networks. Ground motion model plays an important role in the PSHA as it characterizes the ground motion subjected to quakes of all magnitudes and distances. In the regions where there is lack of adequake seismic data, ground motions are estimated by modeling the region using the numerical techniques. However, these approaches are computationally expensive to generate ground motion for broad range of frequencies. Hence, traditionally, the empirical GMMs are developed using the available ground motion records. However, in regions where there is lack of adequate data such as stable continental regions very few GMMs are available and the available GMMs are also for the quakes of magntiudes > 4. Neverthless, it is observed that the large quakes influence the seismic hazard relatively compared to small-moderate quakes. However, the hazard associated to the small quakes in other planetary bodies such as Moon or Mars is similar to hazard associated with the large quakes on the Earth. Hence, in the current study, a GMM that can predict the possible ground motion at a site subjected to both small-moderate or large quakes has been developed using the machine learning technqiues which can be further used in the seismic hazard studies of planetary bodies for their upcoming seismic missions. Moreover, the GMM developed for small-moderate quakes is further validated and compared with the ground motion data available from the Mars Insight mission.As the comparisions proved that the dveeloped GMM for stable continental regions is efficient, the seismic hazard in the lunar south pole region (LSPR) . The outcome of this study helped in identifying the potential landing site among the 13 candidate regions provided by NASA-Artemis III mission. Additionally, the present study also provides alternate ground motion parameter i.e, equivalent velocity (energy based parameter) that can better characterize the ground motion in terms of damage. Using the energy based GMM, PSHA analyis has been conducted for stable continental region i.e, Pensinsular India as a further application.