This study suggests a prediction model of ground motion spectral shape considering characteristics of earthquake records in Korea. Based on the Graizer and Kalkan’s prediction procedure, a spectral shape model is defined as a continuous function of period in order to improve the complex problems of the conventional models. The approximate spectral shape function is then developed with parameters such as moment magnitude, fault distance, and average shear velocity of independent variables. This paper finally determines estimator coefficients of subfunctions which explain the corelation among the independent variables using the nonlinear optimization. As a result of generating the prediction model of ground motion spectral shape, the ground motion spectral shape well estimates the response spectrum of earthquake recordings in Korea.