A methodology to model seismic microzonation maps is required in the hazard mitigation decision plans of the earthquake prone areas. The stage of disaster preparedness for new residential places is of great importance for detailed seismic microzonation models. The effects of local geological and geotechnical site conditions were considered in order to establish site characterization as the initial stage of the models in this study. Dynamic soil properties based on the empirical correlations between shear wave velocity (V (s)) and standard penetration test blow counts were taken into account in order to define representative soil profiles extending down to the engineering bedrock. One-dimensional site response analyses were performed to analyze earthquake characteristics on the ground surface. The layers for soil classification, geology, depth to groundwater level, amplification, distance to fault, slope and aspect, and liquefaction-induced ground deformation potential of the study area were prepared in seismic microzonation models. The study area, Erbaa, is placed along the seismically active North Anatolian Fault Zone. Final seismic microzonation map of the study area was evaluated applying different GIS-based Multi-Criteria Decision Analysis (MCDA) techniques. Two of the MCDA techniques, simple additive weighting and analytical hierarchical process (AHP), are considered during the evaluation step of the final seismic microzonation map. The comparison is made in order to distinguish two different maps based on these MCDA techniques. Eventually, AHP-based seismic microzonation map is more preferable for the seismic design purposes in this study.