High Resolution Velocity Structure in Eastern Turkey

We investigate the crust and upper mantle structure of eastern Turkey where the Anatolian, Arabian and Eurasian Plates meet, forming a complex tectonic regime. The Bitlis suture is a continental collision zone between the Anatolian plateau and the Arabian plate. Broadband data available through the Eastern Turkey Seismic Experiment (ETSE) provide a unique opportunity for studying the high resolution velocity structure of the region. Zor et al. (2003) found an average 46 km thick crust in the Anatolian plateau using a six-layered grid search inversion of the ETSE receiver functions. Receiver functions are sensitive to the velocity contrast of interfaces and the relative travel time of converted and reverberated waves between those interfaces. The interpretation of receiver functions alone, however, may result in an apparent depth-velocity trade-off [Ammon et al., 1990]. In order to improve upon this velocity model, we have combined the receiver functions with surface wave data using the joint inversion method of Julia et al. (2000). In this technique, the two sets of observations are combined into a single algebraic equation and each data set is weighted by an estimate of the uncertainty in the observations. The receiver functions are calculated using an iterative time-domain deconvolution technique. We also consider azimuthal changes in the receiver functions and have stacked them into different groups accordingly. We are improving our surface wave model by making Love and Rayleigh dispersion measurements at the ETSE stations and incorporating them into a regional group velocity model for periods between 10 and 100 seconds. Preliminary results indicate a strong trend in the long period group velocities toward the northeast, indicating slow upper mantle velocities in the area consistent with Pn, Sn and receiver function results. Starting models used for the joint inversions include both a 1-D model from a 12-ton dam shot recorded by ETSE [Gurbuz et al., 2004] and the models from the original receiver function inversions. We observe that the joint inversion results are independent of the starting model and converge to the same final model, with some differences compared to the original profiles. While we don't observe significant changes in the first order discontinuities of the model, such as Moho depth, we are better able to resolve features in the crust.