The Thessaly basin, the largest basin in central Greece, has a well-known history of large earthquakes, with mainshocks having typical magnitudes between 6.0 and 7.0. During the 20thcentury, eight major seismic sequences with mainshock magnitudes equal or larger than 6.0 occurred in the area (1905, 1911, 1930, 1941, 1954, 1955, 1957, 1980). The purpose of the present study is to simulate the seismic motions of the M=7.0 1954 Sofades earthquake, which had observed macroseismic intensities up to 9.5, resulting in heavy damages of the town and villages of the broader southern Thessaly area. This event marked the beginning of a series of earthquakes along the southern Thessaly rupture zone (Central Greece), that also involved the destructive events of 1955 and 1957 which mostly affected the city of Volos and the whole eastern part of the southern Thessaly basin (Papastamatiou and Mouyaris, 1986). The 1980 Almyros earthquake is another notable, recent seismic event with similar damaging consequences.