

By scanning the amplitude spectra of the Ps and PpPs phases, we can obtain precise arrival times of these phases for interfaces above 250 km (errors < 0.13 s) and derive the average velocities and thicknesses of the layers between interfaces in the crust and upper mantle. In this study, we propose a receiver function velocity analysis technique (RFVAT) that can constrain the velocity structure above discontinuities and remove multiple reflections arising from target discontinuities. Converted phases of crustal and upper mantle discontinuities are often interfered by multiple reflections from shallower structures, adding further complexity to the interpretation of P wave RF images. Results of P wave receiver function (RF) analysis (e.g., H‐κ stacking, common conversion point (CCP) stacking, and migration) depend on the velocity models employed. As a result, the findings of this study provide useful new constraints to the crustal thickness and composition in the study region. The prominent crustal low-velocity zones especially around the Nemrut and Ararat volcanoes are identified below ∼20 km depth and likely imply the presence of partial melt zones or active magma chambers in the region. A high average crustal Vp/Vs ratio of 1.82–1.87 is found near volcanic centers, characterized by predominantly a mafic crustal composition. The southwest of study region is characterized overall by a low average crustal Vp/Vs ratio (less than 1.71), implying a felsic crustal composition.

Inversion results indicate that Moho depth varies from at most ∼40 km near the Bitlis-Zagros thrust zone to 45 km in the northeastward of Quaternary volcanic centers. For slant stacking method, radial receiver functions are calculated using deconvolution based on frequency domain. Further, slant stacked method, employed for the first time in the study region, is applied to identify the arrival times of P to S conversion and reverberations (Ps and PpPms) that give Vp/Vs ratio beneath each station. Teleseismic receiver functions obtained through deconvolution in the time domain are used to estimate crustal discontinuities and Moho interface. This study aims to explore the crustal seismic properties across the East Anatolian volcanic belt.
