## 1 Introduction  
The state of Odisha is famous for its cultural heritage, historical monuments, and hot springs, and is considered as the center for economic activity in central-eastern India. However, Odisha is seismically active due to ancient faults developed in the Mahanadi rift basins. Earthquakes of intra-cratonic events have occurred in western and northeastern Odisha [1]. Low-to moderate-events have been reported by the United States Geological Survey (USGS) and National Center for Seismology (NCS). The main reasons for earthquake risk in this region could be the intra-cratonic seismicity and Bay of Bengal events. Odisha has experienced several low to moderate events with magnitudes in excess of 5.8 at a distance of 310 km southeast of Bhubaneswar. The Mahanadi rift basin elongated in a NW-SE orientation and divides the land into two parts: the North Odisha Boundary Fault (NOBF) and Mahanadi shear zone (MSZ) [2; 3]. Several sediment-filled basins such as Talcher, IB River, and Jajpur-Chandbali basins were developed in the northwest while Athgarh and Salepur basins are in the southern part of Mahanadi rift valley. However, major faults have indicated reactivation during the Holocene epoch, including the Brahmani Fault. A graben structure has formed in the Mahanadi valley, where many deep-seated faults are situated. [PERSON] et al. [4] described that proximity to faults does not affect the hazard, as earthquake damage depends on several factors, i.e., subsurface geology. The intensity at a particular location depends on amplification values of near- and sub-surface lithology and associated tectonic structures [4]. Several studies on soil amplification have been conducted showing local amplification with unconsolidated sedimentary deposits [5; 6]. The surface seismic waves, such as R and L waves, were usually amplified due to alluvium deposits [5]. [PERSON] et al. [6] studied earthquake probability (EPA) by using several indicators, such as sub-surface lithology, i.e., the amplification factor. They achieved an accuracy of 84.5% based on the artificial neural network (ANN) technique. Earth's crust was convoluted due to complicated tectonic activity, where lateral and vertical heterogeneity in soil and rocks could be observed in terms of thickness, composition, rock density, and strength [7]. Therefore, nature's heterogeneity could be the cause for energy trapping and shear wave interference, generating a response to shallow subsurfaces in terms of ground motion [7]. [PERSON] [8] interpreted that grabens, depressions, and basins carry immature loose sediments with causal negative curvature factors and a 10\({}^{\circ}\) slope gradient. Therefore, wave amplification in several basins including the Mahanadi basin can increase damages even by a moderate magnitude event [9; 10]. Thus, India's existing seismic zonation map did not count site-specific peak ground acceleration (PGA) parameters nor response spectra-based amplification factors on rock types [9]. To deal with this problem, [PERSON] et al. [11] assessed the ground motion seismic hazard using a probabilistic seismic hazard assessment (PSHA) approach in India to report PGA values using the global seismic hazard assessment (GSHAP). Eventually, [PERSON] et al. [12] conducted the deterministic seismic hazard assessment (DSHA) for India, using seismogenic zones, historical inventory, infrastructure model, active tectonics, and focal mechanisms. Consequently, several other FSHA-based studies were conducted in Bangalore and Mumbai [13; 14; 15]. Lack of relevant local attenuation relations became an issue that can cause ground motion uncertainty. However, in the FSHA method based on attenuation relation [16; 17; 18], it can be used for ground motion study. [PERSON] and [PERSON] [13] investigated and observed high seismic hazards for a magnitude of 8.0 Mw in relation to the existing code IS 1893-2002. They observed that PGA values were overestimated and did not comply with the earthquake inventory.