Earth Planets Space, 57, 1019–1038, 2005 Three-dimensional structure of Vp, Vs and Vp/Vs in the upper crust of the Marmara region, NW Turkey

We applied a 3-D seismic tomography inversion algorithm to arrival-time data obtained, during 18 years from 1985 to 2002, from local seismic networks and aftershock studies in the Marmara region, in order to better understand the upper crustal structure of the complex tectonic region. We integrated all the available data set into a common data set and relocated the events, using a 1-D velocity model. We then selected 3,949 earthquakes and obtained 92355 arrival times, in total, consisting of 59,313 P-wave and 33,042 S-wave arrival times. In this paper we present detailed crustal structures for V p, V s and V p/V s ratios from the surface down to 15 km depth with good resolution in terms of the hit count analysis of seismic rays, the checkerboard and restoring resolution tests for the studied region. The results obtained from the inversion suggest that the western part of the North Anatolian Fault Zone shows strong lateral heterogeneity. We concluded that no clear pattern exists between the distribution of microearthquakes and aftershocks and the velocity perturbations presented in this paper. This is probably due to complex tectonic and geological structures. Large coseismic slip associated with the two recent strong earthquakes was found to correspond to higher velocity anomalies, as was often found recently. On the other hand, low velocity values correspond to the sedimentary units or the alluvium regions, as supported by the low resistivity and gravity values. Thus the results presented in this paper are quite consistent with the other data such as gravity, resistivity and magnetic anomalies, indicating that our model is reliable and efficient and i