- We apply an algorithm based on the modal summation method to theoretically estimate the site effect at selected locations underlain by different geological formations within the city of Thessaloniki (Greece). Complete strong motion synthetics are constructed for all components of motion at each site, for a maximum frequency of 10Hz. The anelastic, local 1-D velocity models are based on cross-hole data. Four point sources with different azimuths and distances from the city are used to compute the input signals. The theoretical amplification is estimated through spectral ratios of accelerograms obtained by the local 1-D over those obtained by the regional 1-D velocity model. The results from the numerical modeling are compared with those derived from experimental techniques, such as of Standard Spectral Ratio and Horizontal-to-Vertical Spectral Ratio, which had been applied to acceleration data recorded at the same sites. The comparison demonstrates that the theoretical amplifications based on known and simple subsurface geology can be used as a first-order estimate, while for cases of more complex geometries the use of at least 2-D modeling in site effects estimation is mandatory.