SUMMARYSoil e ects on peak ground acceleration, velocity and elastic response spectra (5% damping) are expressed by simple approximate relations in terms of ÿve key parameters: (a) the fundamental vibration period of the non-linear soil, T S , (b) the period of a bedrock site of equal thickness, T b , (c) the predominant excitation period, Te, (d) the peak seismic acceleration at outcropping bedrock, a b max , and (e) the number of signiÿcant excitation cycles, n. Furthermore, another relation is proposed for the estimation of T S in terms of the soil thickness H , the average shear wave velocity of the soil V S; o and a b max . The aforementioned parameters were ÿrst identiÿed through a simpliÿed analytical simulation of the site excitation. The multivariable approximate relations were then formulated via a statistical analysis of relevant data from more than 700 one-dimensional equivalent-linear seismic ground response analyses, for actual seismic excitations and natural soil conditions. Use of these relations to back-calculate the numerical results in the database gives an estimate of their error margin, which is found to be relatively small and unbiased. The proposed relations are also independently veriÿed through a detailed comparison with strong motion recordings from seven well-documented case studies: (a) two sites in the San Fernando valley during the Northridge earthquake, and (b) ÿve di erent seismic events recorded at the SMART-1 accelerometer array in Taiwan. It is deduced that the accuracy of the relations is comparable to that of the equivalent-linear method. Hence, they can be readily used as a quick alternative for routine applications, as well as for spreadsheet computations (e.g. GIS-aided seismic microzonation studies) where numerical methods are cumbersome to implement.