In this study, we aimed to explore the S-wave velocity structure of shallow soils using microtremors in order to estimate site responses in Tekirdag and surrounding areas (NW Turkey). We collected microtremor array data at 44 sites in Tekirdag, Marmara Ereglisi, Corlu, and Muratli. The phase velocities of Rayleigh waves were estimated from the microtremor data using a Spatial Autocorrelation method. Then, we applied a hybrid genetic simulated annealing algorithm to obtain a 1D S-wave velocity structure at each site. Comparison between the horizontal-to-vertical ratio of microtremors and computed ellipticities of the fundamental mode Rayleigh waves showed good agreement with validation models. The depth of the engineering bedrock changed from 20 to 50 m in the Tekirdag city center and along the coastline with a velocity range of 700-930 m/s, and it ranged between 10 and 65 m in Marmara Ereglisi. The average S-wave velocity of the engineering bedrock was 780 m/s in the region. We obtained average S-wave velocities in the upper 30 m to compare site amplifications. Empirical relationships between the AVs30, the site amplifications, and also average topographic slopes were established for use in future site effects microzonation studies in the region.