In current practice, the seismic hazard curves are generated by combining all possible choices for seismotectonic model (SM), earthquake magnitude recurrence models, and ground motion prediction equation (GMPE) models using logic trees. However, studies have shown that significant uncertainties exist in the evaluation of the seismic hazard, especially in low-to-moderate seismicity regions, such as Central \& Eastern United States and Europe, due to the scarcity of the recorded earthquake data and lack of precise knowledge concerning the earthquake generation and propagation mechanisms. Specifically, these uncertainties lead to conservatisms in the seismic hazard which poses a considerable risk for nuclear facilities. Designing or upgrading facilities with such conservative estimates of seismic hazard results in large capital costs. In standard probabilistic seismic hazard assessment (PSHA) practices, the estimate of branch weights for the GMPE models is based on expert judgment and there is no holistic approach for estimating the weights under uncertainty. In this research, we focus on the GMPE part of the PSHA and we propose a statistical methodology that is based on Bayesian inference, to evaluate the GMPE model weights from a list of GMPE models in the logic tree. In this study, Bayesian model averaging (BMA) approach with Bayesian linear models (BLM) is employed to combine GMPE models and the performance of the BMA model is tested against European Strong-Motion (ESM) database. It is shown that the proposed methodology can calculate the weights at different time periods without any approximations using analytical formulas and it is computationally inexpensive.