The Harmonization of Seismic Hazard Maps in the Western Balkan Countries Project (BSHAP) was funded for 7 years by NATO Science for Peace and Security Program to support the preparation of new seismic hazard maps of the Western Balkan Region. In the framework of BSHAP, the regional free field strong motion network capacity was significantly increased by deployed recorders (Gülerce et al. in Final report of improvements in the Harmonized Seismic Hazard Maps for the Western Balkan Countries (BSHAP-2) Project, NATO SfP-984374. http://wbalkanseismicmaps.org/ , 2015) and the BSHAP strong motion database that includes both pre-BSHAP (mostly analogue) and post-BSHAP (all digital) recordings was compiled. The BSHAP strong motion database includes uniformly processed strong motions along with the related earthquake metadata and station information; therefore, it provides a solid base for the ground motion characterization studies in the region. The established database is used for selection of the ground motion prediction equations (GMPEs) to be employed in the probabilistic seismic hazard assessment by comparing the compiled strong ground motions with the predictions of candidate global and Euro-Mediterranean GMPEs in a systematic manner. Details of the selection strategy including the analysis of model predictions for ground motion scaling, data testing to evaluate the model performance, and ranking is provided here along with the GMPE logic tree used in the proposed seismic hazard maps by Kuka et al. (Bull Earthq Eng, under review).