A new set of empirical equations for prediction of displacement response spectral ordinates from 20 Hz to T = 20 s is illustrated. The coefficients of the equations were obtained by regressing a dataset based on 1,155 tri-axial digital and 9 analog accelerometer records from 60 earthquakes worldwide. Long period disturbances in the accelerograms were evaluated and removed using a very recent method, aimed at preserving the long-period spectral content of the records. Analysis of variance has disclosed only little evidence for regional dependence of ground motions, while a carefully conducted evaluation of site effects resulted in clearly differentiated spectral amplification bands associated to the main ground types B, C, and D of Eurocode 8. Spectral ordinates for vibration periods >5 s were found to scale with magnitude quite consistently with theoretical scaling from Brune’s model. On the other hand, comparison of results with those yielded by recent prediction models in Europe and the United States (NGA), indicated that the latter may not be uniformly reliable at long periods. The proposed empirical equations are easily implemented in computer programs for seismic hazard assessment, being characterized by a simple functional form and a restricted number of predictor variables.