A simplified methodology is presented for the computation of elastic response spectra at liquefiable sites, based on the widely used equivalent linear site response analysis method. The target response spectrum is obtained by frequency dependent interpolation between the response spectra for “non-liquefied” and for “liquefied” ground, the former computed for the native soil properties, without excess pore pressure buildup, and the latter for the properties of the liquefied soil. Apart from the structural frequency, the interpolation coefficient is also related to the factor of safety against liquefaction. In this way, the proposed methodology takes indirectly into account the initial segment of the seismic excitation time-history, until the onset of liquefaction, which is of cornerstone importance for the overall ground response at sites with low and moderate liquefaction potential. Results of parametric nonlinear numerical analyses and three field case studies are used to calibrate the methodology and evaluate its accuracy. In addition, the potential of the proposed methodology for defining code-oriented design spectra for liquefiable sites is examined.