Seismic intensity maps are employed globally in the aftermaths of earthquakes for rapid response purposes. These maps involve correlations between intensity and peak ground-motion values. In this study, we focus on eastern sections of the North Anatolian fault zone (NAFZ). The eastern segments of the NAFZ are less investigated and have sparse seismic networks compared with western ones. In particular, we study Erzincan, which is a small city in eastern Turkey, located in the conjunction of three active faults: North Anatolian, North East Anatolian, and Ovacik faults. Erzincan city center is in a pull-apart basin underlain by soft sediments, which significantly amplify the ground motions. Combination of the tectonic and geological settings of the region have led to destructive earthquakes such as the 27 December 1939 (Ms = 8.0) and the 13 March 1992 (Mw = 6.6) events resulting in extensive losses. In this study, we initially perform ground-motion simulations for a set of scenario events as well as the 1992 Erzincan earthquake. We then use local relationships between Modified Mercalli Intensity (MMI) and Peak Ground Acceleration (PGA) as well as Peak Ground Velocity (PGV) to obtain the corresponding MMI values. We present our results in the form of synthetic intensity maps for the 1992 event and the scenario earthquakes