AbstractFollowing the 2016 Chauk earthquake in Bagan (Myanmar), numerous old pagodas and temples suffered severe damage. This research presents a study on the seismic hazard analysis of the Bagan city in Myanmar, based on a probabilistic framework focussing on analysing 43 temples with their associated local soil information. To this end, two seismic source models are developed based on the tectonic setting of the region and information available. Instrumental and historical records are compiled from both literature and international earthquake catalogues while conducting catalogue completeness. This study uses state-of-the-art ground motion models to perform probabilistic seismic hazard analysis and develop seismic hazard maps for different return periods in the region. Results are also expressed for selected temples in the region in terms of site-specific uniform hazard spectra. The findings indicate significant seismic activity, with peak ground acceleration in the region ranging from 0.25 to 0.36 g for a return period of 2475 years, 0.22–0.32 g for a return period of 975 years, and 0.18–0.24 g for a return period of 475 years. The updated hazard levels indicate that the literature slightly underestimates hazard in the region under study.