The local effects of the construction site play an important role in the seismic motion characteristics and thus the movements of the earth's surface, because it can change the seismic movements to the intensity of large or small before reaching the Earth's surface. Ground response analysis with the goal of predicting the movements of the earth's surface and developing site-specific response spectra to assess the risks arising from the earthquake liquefaction and calculating the forces that can cause instability of earth retaining structures. In the present research, a range of accelerograms recorded in the seismic bedrock level in two fields (near and far) have been chosen for the purpose of deriving site-specific response spectra. The study for the area located at the city of Anzali(Kapoorchal). In the following, shear wave velocity profiles are derived using the relationships between shear wave velocity and number of SPT in the area. Then the dynamic response in surface is done with linear, equivalent linear and nonlinear methods by applying the accelerograms on the bedrock to soil column profile. Then, it used as in the surface for deriving site response spectra. The results show that the linear equivalent model offers satisfying results, but this method is not recommended with increasing the acceleration of the input mtion because dynamic characteristics of the soil is strongly nonlinear and dependant to the strain and stresses in high strain levels.