Structure of the crust and African slab beneath the central Anatolian plateau from receiver functions: New insights on isostatic compensation and slab dynamics

The central Anatolian plateau in Turkey is a region with a long history of subduction, continental collision, accretion of continental fragments, and slab tearing and/or breakoff and tectonic escape. Central Anatolia is currently characterized as a nascent plateau with widespread Neogene volcanism and predominantly transtensional deformation. To elucidate the present-day crustal and upper mantle structure of this region, teleseismic receiver functions were calculated from 500 seismic events recorded on 92 temporary and permanent broadband seismic stations. Overall, we see a good correlation between crustal thickness and elevation throughout central Anatolia, indicating that the crust may be well compensated throughout the region. We observe the thickest crust beneath the Taurus Mountains (>40 km); it thins rapidly to the south in the Adana Basin and Arabian plate and to the northwest across the Inner Tauride suture beneath the Tuz Gölü Basin and Kırşehir block. Within the Central Anatolian Volcanic Province, we observe several low seismic velocity layers ranging from 15 to 25 km depth that spatially correlate with the Neogene volcanism in the region, and may represent crustal magma reservoirs. Beneath the central Taurus Mountains, we observe a positive amplitude, subhorizontal receiver function arrival below the Anatolian continental Moho at ∼50–80 km that we interpret as the gently dipping Moho of the subducting African lithosphere abruptly ending near the northernmost extent of the central Taurus Mountains. We suggest that the uplift of the central Taurus Mountains (∼2 km since 8 Ma), which are capped by flat-lying carbonates of late Miocene marine units, can be explained by an isostatic uplift during the late Miocene–Pliocene followed by slab breakoff and subsequent rebound coeval with the onset of faster uplift rates during the late Pliocene–early Pleistocene. The Moho signature of the subducting African lithosphere terminates near the southernmost extent of the Central Anatolian Volcanic Province, where geochemical signatures in the Quaternary volcanics indicate that asthenospheric material is rising to shallow mantle depths

Pasif sismik veriler kullanılarak Türkiye boyunca yerkabuğu ve mantonun analizi.

Anatolia is one of the most tectonically active regions in the world which is amalgamation of different terranes and structures. This region has a complextectonic evolution which includes various stages of subduction and collision. This complexity demands more detailed observations of the crust and upper mantle to generate a comprehensive geodynamic model. However, the knowledge related to crust and mantle structure beneath Anatolia is still very limited so the fates of the structures in deeper interiors are enigmatic. During the last decade, the number of broadband seismic stations operating in Turkey has been expanded significantly which provides a unique opportunity to analyze the crust and mantle in detail. Recently, additional temporary broadband stations are deployed with our collaboration across central Turkey in the scope of an international project funded by National Science Foundation (NSF). In this study, we analyzed the broadband seismic data recorded during earthquakes by national networks and recently deployed temporary stations including the CAT and also the NAF experiments. Our aim is to image the inner structure of crust and upper mantle, detect thicknesses of the crust and analyze the multi-layered seismic anisotropy along crust and upper mantle. At the end, our results are interpreted and modeled along with previous geological and geophysical data to identify the factors controlling the tectonic settings and test plausible geodynamic models.Ph.D. - Doctoral Progra

M-Split: a graphical user interface to analyze multilayered anisotropy from shear wave splitting

Shear wave splitting analysis are commonly used to infer deep anisotropic structure. For simple cases, obtained delay times and fast-axis orientations are averaged from reliable results to define anisotropy beneath recording seismic stations. However, splitting parameters show systematic variations with back azimuth in the presence of complex anisotropy and cannot be represented by average time delay and fast axis orientation. Previous researchers had identified anisotropic complexities at different tectonic settings and applied various approaches to model them. Most commonly, such complexities are modeled by using multiple anisotropic layers with priori constraints from geologic data. In this study, a graphical user interface called M-Split is developed to easily process and model multilayered anisotropy with capabilities to properly address the inherited non-uniqueness. M-Split program runs user defined grid searches through the model parameter space for two-layer anisotropy using formulation of Silver and Savage (1994) and creates sensitivity contour plots to locate local maximas and analyze all possible models with parameter tradeoffs. In order to minimize model ambiguity and identify the robust model parameters, various misfit calculation procedures are also developed and embedded to M-Split which can be used depending on the quality of the observations and their back-azimuthal coverage. Case studies carried out to evaluate the reliability of the program using real noisy data and for this purpose stations from two different networks are utilized. First seismic network is the Kandilli Observatory and Earthquake research institute (KOERI) which includes long term running permanent stations and second network comprises seismic stations deployed temporary as part of the “Continental Dynamics-Central Anatolian Tectonics (CD-CAT)” project funded by NSF. It is also worth to note that M-Split is designed as open source program which can be modified by users for additional capabilities or for other applications

Subsurface signature of North Anatolian Fault Zone and its relation with old sutures: New insight from receiver function analysis

The North Anatolian Fault Zone (NAFZ) is an active continental transform plate boundary that accommodates the westward extrusion of the Anatolian plate. The central segment of NAFZ displays northward convex surface trace which coincides partly with the Paleo-Tethyan suture formed during the early Cenozoic. The depth extent and detailed structure of the actively deforming crust along the NAF is still under much debate and processes responsible from rapid uplift are enigmatic. In this study, over five thousand high quality P receiver functions are computed using teleseismic earthquakes recorded by permanent stations of national agencies and temporary North Anatolian Fault Passive Seismic experiment (2005-2008). In order to map the crustal thickness and Vp/Vs variations accurately, the study area is divided into grids with 20 km spacing and along each grid line Moho phase and its multiples are picked through constructed common conversion point (CCP) profiles. According to our results, nature of discontinuities and crustal thickness display sharp changes across the main strand of NAFZ supporting a lithospheric scale faulting that offsets Moho discontinuity. In the southern block, crust is relatively thin in the west (∼35 km) and becomes thicker gradually towards east (∼40 km). In contrast, the northern block displays a strong lateral change in crustal thickness reaching up to 10 km across a narrow roughly N-S oriented zone which is interpreted as the subsurface signature of the ambiguous boundary between Istanbul Block and Pontides located further west at the surface

Water content estimated from point scale to plot scale

Soil moisture controls the portioning of rainfall into infiltration and runoff. Here we investigate measurements of soil moisture using a range of techniques spanning different spatial scales. In order to understand soil water content in a test basin, 512 km2 in area, in the south of Turkey, a Cosmic Ray CRS200B soil moisture probe was installed at elevation of 1459 m and an ML3 ThetaProbe (CS 616) soil moisture sensor was established at 5cm depth used to get continuous soil moisture. Neutron count measurements were corrected for the changes in atmospheric pressure, atmospheric water vapour and intensity of incoming neutron flux. The calibration of the volumetric soil moisture was performed, from the laboratory analysis, the bulk density varies between 1.719 (g/cm3) -1.390 (g/cm3), and the dominant soil texture is silty clay loam and silt loamThe water content reflectometer was calibrated for soil-specific conditions and soil moisture estimates were also corrected with respect to soil temperature. In order to characterize the subsurface, soil electrical resistivity tomography was used. Wenner and Schlumberger array geometries were used with electrode spacing varied from 1m- 5 m along 40 m and 200 m profiles. From the inversions of ERT data it is apparent that within 50 m distance from the CRS200B, the soil is moderately resistive to a depth of 2m and more conductive at greater depths. At greater distances from the CRS200B, the ERT results indicate more resistive soils. In addition to the ERT surveys, ground penetrating radar surveys using a common mid-point configuration was used with 200MHz antennas. The volumetric soil moisture obtained from GPR appears to overestimate those based on TDR observations. The values obtained from CS616 (at a point scale) and CRS200B (at a mesoscale) are compared with the values obtained at a plot scale. For the field study dates (20-22.06.2017) the volumetric moisture content obtained from CS616 were 25.14%, 25.22% and 25.96% respectively. The values obtained from CRS200B were 23.23%, 22.81% and 23.26% for the same dates. Whereas the values obtained from GPR were between 32%-44%. Soil moisture observed by CRS200B is promising to monitor the water content in the soil at the mesoscale and ERT surveys help to understand the spatial variability of the soil water content within the footprint of CRS200B

The effects of subduction termination on the continental lithosphere: Linking volcanism, deformation, surface uplift, and slab tearing in central Anatolia

Subduction beneath central Anatolia represents the transition between continuous subduction along the Aegean trench in the west and slab break-off and/or subduction termination at the Arabian-Eurasian collision zone in the east. Using recently collected seismic data from the Continental Dynamics-Central Anatolian Tectonics project alongside a newly developed approach to the creation of a 3D shear-velocity model from the joint inversion of receiver functions and surface-wave dispersion data, we can gain important insights into the character of the downgoing, segmenting African lithosphere and its relationship with the overriding Central Anatolian plate