Three-dimensional thermal structure of the Chinese continental crust and upper mantle

We invert S-wave velocities for the 3D upper-mantle temperatures, in which the position with a temperature crossing the 1300℃ adiabat is corresponding to the top of the seismic low velocity zone. The temperatures down to the depth of 80 km are then calculated by solving steady-state thermal conduction equation with the constraints of the inverted upper-mantle temperatures and the surface temperatures, and then surface heat flows are calculated from the crustal temperatures. The misfit between the calculated and observed surface heat flow is smaller than 20% for most regions. The result shows that, at a depth of 25 km, the crustal temperature of eastern China (500―600℃) is higher than that of western China (<500℃). At a depth of 100 km, temperatures beneath eastern and southeastern China are higher than the adiabatic temperature of 1300℃, while that beneath west China is lower. The Tarim craton and the Sichuan basin show generally low temperature. At a depth of 150 km, temperatures beneath south China, eastern Yangtze craton, North China craton and around the Qiangtang terrane are higher than the adiabatic temperature of 1300℃, but is the lowest beneath the Sichuan basin and the regions near the Indian-Eurasian collision zone. At a depth of 200 km, very low temperature occurs beneath the Qinghai-Tibet Plateau and the south to the Tarim craton.     

S-wave crustal and upper mantle’s velocity structure in the eastern Tibetan Plateau — Deep environment of lower crustal flow

A teleseismic profile consisting of 26 stations was deployed along 30°N latitude in the eastern Tibetan Plateau. By use of the inversion of P-wave receiver function, the S-wave velocity structures at depth from surface to 80 km beneath the profile have been determined. The inversion results reveal that there is significant lateral variation of the crustal structure between the tectonic blocks on the profile. From Linzhi north of the eastern Himalayan Syntaxis, the crust is gradually thickened in NE direction; the crustal thickness reaches to the maximum value (∼72 km) at the Bangong-Nujiang suture, and then decreased to 65 km in the Qiangtang block, to 57–64 km in the Bayan Har block, and to 40–45 km in the Sichuan Basin. The eastern segment of the teleseismic profile (to the east of Batang) coincides geographically with the Zhubalong-Zizhong deep seismic sounding profile carried out in 2000, and the S-wave velocity structure determined from receiver functions is consistent with the P-wave velocity structure obtained by deep seismic sounding in respect of the depths of Moho and major crustal interfaces. In the Qiangtang and the Bayan Har blocks, the lower velocity layer is widespread in the lower crust (at depth of 30–60 km) along the profile, while there is a normal velocity distribution in lower crust in the Sichuan Basin. On an average, the crustal velocity ratio (Poisson ratio) in tectonic blocks on the profile is 1.73 (σ = 0.247) in the Lhasa block, 1.78 (σ = 0.269) in the Banggong-Nujiang suture, 1.80 (σ = 0.275) in the Qiangtang block, 1.86 (σ = 0.294) in the Bayan Har blocks, and 1.77 (σ = 0.265) in the Yangtze block, respectively. The Qiangtang and the Bayan Har blocks are characterized by lower S-wave velocity anomaly in lower crust, complicated Moho transition, and higher crustal Poisson ratio, indicating that there is a hot and weak medium in lower crust. These are considered as the deep environment of lower crustal flow in the eastern Tibetan Plateau. Flowage of the ductile material in lower crust may be attributable to the variation of the gravitational potential energy in upper crust from higher on the plateau to lower off plateau. Supported by the National Natural Science Foundation of China (Grants No. 40334041 and 40774037) and the International Cooperation Program of the Ministry of Science and Technology of China (Grant No. 2003DF000011)     

Crustal structure in Dabieshan UHP metamorphic belt and its tectonic implication

IntroductionGeologistsfirstlyfoundcoesite-bearingecologitesattheendof1980'sandthenthemicrodiamond(Xu,elal,1992)inDabieshanarea.Theultra-highpressure(UHP)metamorphismandthegeodynamicprocessesofDabieorogenhaveattractedmanygeoscientists(Wang,etal,1989,O...     

2-D crustal Poisson’s ratio from seismic travel time inversion in Changbaishan Tianchi volcanic region

Based on the inversion method of 2D velocity structure and interface, the crustal velocity structures of P-wave and S-wave along the profile L ₁ are determined simultaneously with deep seismic sounding data in Changbaishan Tianchi volcanic region, and then its Poisson’s ratio is obtained. Calculated results show that this technique overcomes some defects of traditional forward calculation method, and it is also very effective to determine Poisson’s ratio distribution of deep seismic sounding profile, especially useful for study on volcanic magma and crustal fault zone. Study result indicates that there is an abnormally high Poisson’s ratio body that is about 30 km wide and 12 km high in the low velocity region under Tianchi crater. Its value of Poisson’s ratio is 8% higher than that of surrounding medium and it should be the magma chamber formed from melted rock with high temperature. There is a high Poisson’s ratio zone ranging from magma chamber to the top of crust, which may be the uprise passage of hot substance. The lower part with high Poisson’s ratio, which stretches downward to Moho, is possibly the extrusion way of hot substance from the uppermost mantle. The conclusions above are consistent with the study results of both tomographic determination of 3D crustal structure and magnetotelluric survey in this region. Foundation item: Key Project from China Earthquake Administration and the Project (95-11-02-01) from Ministry of Science and Technology (2001DIA10003). Contribution No. RCEG200401, Geophysical Exploration Center, China Earthquake Administration.     

Structure of the crust and uppermost mantle from broadband seismic array in the western Dabie Mountains,east central China

A broadband seismic array of 7 stations was set up in the western Dabie Mountains (31°20′-31°50′N, 114°30′-115°E). Teleseismic events from May 2001 to November 2001 were collected and analyzed by radial receiver function to determine the S-wave velocity structure of the crust and uppermost mantle. The crustal thickness is 32-38 km beneath the array. The crust-mantle boundary appears as a gently north-dipping velocity discontinuity, but turns to be a velocity gradient beneath a station near the Qiliping shea...     

Deep tectonic setting of the 2008 Wenchuan <Emphasis Type="Italic">M</Emphasis><Subscript>s</Subscript>8.0 earthquake in southwestern China

Teleseismic P-wave receiver functions at 20 broadband seismic stations in the Longmenshan fault zone (LMFZ) and its vicinity were extracted, and the crustal thickness and the P- and S-wave velocity ratio were calculated by use of the H-k stacking algorithm. With the results as constraints, the S-wave velocity structures beneath each station were determined by the inversion of receiver functions. The crustal structure of the Rear-range zone is similar to that of the Songpan-Garze Block, whereas the velocity structure of the Fore-range zone resembles that of Sichuan Basin, implying that the Central Principal Fault of LMFZ is the boundary between the eastern Tibetan Plateau and the Yangtze Block. Lower velocity zone exists in lower crust of the Songpan-Garze Block and the central-southern segment of the Rear-range zone, which facilitates the detachment of the material in upper and middle crust. Joint analysis of the receiver functions and the Bouguer gravity anomalies supports the thesis on the detachment-thrust mode of the LMFZ. A double-detachment pattern is suggested to the tectonic setting in the Songpan-Garze Block. The upper detachment occurs at the depth of 10-15 km, and represents a high-temperature ductile shear zone. There is a lower detachment at the depth of about 30 km, below which the lower crust flow exists in the eastern Tibetan Plateau. Interpretation of the Bouguer gravity anomalies indicates that the Sichuan Basin is of higher density in upper and middle crust in comparison with that of the Songpan-Garze Block. The LMFZ with higher density is the result from the thrusting of the Songpan-Garze Block over the Sichuan Basin. In the lower crust, higher P velocity and higher density in the Sichuan Basin are related to more rigid material, while lower S velocity and lower density in the Songpan-Garze Block are related to the softened and weakened material. The higher density block beneath the Sichuan Basin obstructs the eastward flow of lower crustal material from the Tibetan Plateau, which is driven by the compression of northward movement of Indian Plate. The eastward movement of upper and middle crustal material is also obstructed by the rigid Yangtze Block, resulting in the stress concentrated and accumulated along the LMFZ. When the stress releases sharply, the Wenchuan M _s8.0 earthquake occurs. Supported by the National Natural Science Foundation of China (Grant Nos. 40334041, 40774037) and Joint Foundation of Earthquake Science (Grant No. 1040062)     

Crustal structure beneath the Songpan—Garze orogenic belt

The Benzilan-Tangke deepseismic sounding profile in the western Sichuan region passes through the Song-pan-Garze orogenic belt with trend of NNE.Based on the travel times and the related amplitudes of phases in the record sections,the 2-D P-wave crustal structure was ascertained in this paper.The velocity structure has quite strong lateral variation along the profile.The crust is divided into 5layers,where the first,second and third layer belong to the upper crust,the forth and fifth layer belong to the lower crust.The low velocity anomaly zone gener-ally exists in the central part of the upper crust on the profile,and it integrates into the overlying low velocity basement in the area to the north of Ma‘erkang.The crustal structure in the section can be divided into 4parts:in the south of Garze-litang fault,between Garze-Litang fault and Xianshuihe fault,between Xianshuihe fault and Longriba fault and in the north of Longriba fault,which are basically coincided with the regional tectonics division.The crustal thickness decreases from southwest to northeast along the profile,that is ,from62km in the region of the Jinshajiang River to 52km in the region of the Yellow River.The Moho discontinuity does not obviously change across the Xianshuihe fault basesd on the PmP phase analysis.The crustal average velocity along the profile is lower,about 6.30 km/s.The Benzilan-Tangke profile reveals that the crust in the study area is orogenic.The Xianshuihe fault belt is located in the central part of the profile,and the velocity is positive anomaly on the upper crust,and negative anomaly on the lower crust and upper mantle.It is considered as a deep tectonhic setting in favor of strong earthquake‘s accumulation and occurrence.     

Deep electrical structure of the Sulu orogen and neighboring areas

Because of the discovery of ultrahigh pressure metamorphic (UHPM) belt beneath the Sulu (Jiangsu Province-Shandong Province) orogen, this area has become a focused subject of current geoscience, as it has a close relationship with the evolution of the orogen and the neighboring North China craton. Probing the deep structure beneath this area would be of great significance for the geological interpretation of this issue. In this study, we make an analysis of magnetotelluric (MT) data along a profile across the Sulu orogen to provide evidence of deep structure below this region. The profile begins in west from the North China block, extending in S129°E, across the Tan-Lu fault, Sulu UHPM zone, and Sulu high pressure metamorphic (HPM) zone, and terminates in the Yangtze block in east. We use the nonlinear conjugate gradient method and TE-TM combined mode to perform inversion and interpretation of the MT data, and obtain an electrical structure image above depth of 150 km along the profile. It shows that the structure can be divided into seven sections in lateral direction, between which the electric boundaries coincide well with the major faults, such as the Tan-Lu, Haizhou-Siyang, and Jiashan-Xiangshui faults. In vertical direction the electrical structure can be subdivided into six layers of different resistivities. It is noted that there exist high-conductivity areas in crust below the North China block and Yangtze block, while such a feature is not present beneath the Sulu orogen, which is very different from the Dabie orogen. It is also observed that a fairly continuous zone of relatively low-resistivity exists at depths of 50–90 km of the electrical structure image, which is presumably a weak zone in the uppermost mantle. Just below this low-resistivity zone are the relatively high- resistivity layer of the North China block, relatively low-resistivity layer of the Sulu orogen, and relatively high-resistivity layer of the Yangtze block, all in the shallow upper mantle, respectively. From the whole 2D electrical structure image, there is no abnormally low-resistivity layer in the shallow upper mantle beneath the Sulu orogen and neighboring areas, indicating that no hot asthenoshperic material associated with lithospheric thinning exists at present. Supported by National Natural Science Foundation of China (Grant No. 40534023) and Director Foundation of Institute of Geology, China Earthquake Administration (Grant No. DF-IGCEA-0608-2-16)     

云南思茅—中甸地震剖面的地壳结构

云南思茅—中甸宽角反射/折射地震剖面切割松潘—甘孜、扬子和华南三个构造单元的部分区域. 我们利用初至波和壳内反射波走时层析成像获得地壳纵波速度结构. 在获得新的地壳速度结构模型基础上，利用地震散射成像思想和低叠加次数的叠前深度偏移方法重建了研究区的地壳、上地幔反射结构. 综合分析研究区地壳P波速度模型和壳内地震反射剖面发现：沿测线从北至南地壳厚度从约50 km减薄至35 km左右，地壳厚度的减薄量主要体现在下地壳，剖面北段下地壳厚度约为30 km，剖面南段下地壳厚度仅为15 km左右；上地幔顶部局部位置P波速度值偏低，一般为76～78 km/s，反映出云南地区是典型的构造活动区的特点.剖面沿线地壳内地震反射发育，其中莫霍强反射出现在景云桥下方；在景云桥弧形断裂带8～10 km深处出现宽约50 km的强反射带.     

Geometric form of Haiyuan fault zone in the crustal interior and dynamics implications

The deep seismic reflection data on profile HY2 are reprocessed by the method of simultaneous inversion of velocity distribution and interface position. By the travel-time inversion with the data of the diving wave Pg and fault plane reflection wave, we determine the geometric form and velocity of Haiyuan fault zone interior and surrounding rock down to 10 km depth. The measured data show that the amplitudes have strong attenuation in the range of stake number 37–39 km, suggesting the fault zone has considerable width in the crustal interior. The results of this paper indicate that to the north of the fault zone the crystalline basement interface upheaves gradually from southwest to northeast and becomes shallow gradually towards northeast, and that to the south of the fault zone, within the basin between Xihua and Nanhua mountains, the folded basement becomes shallow gradually towards southwest. The obliquity of the fault zone is about 70° above the 3 km depth, about 60° in the range of the 3–10 km depths. From the results of this paper and other various citations, we believe that Haiyuan fault zone is in steep state from the Earth’s surface to the depth of 10 km. Foundation item: Joint Seismological Science Foundation of China (201001) and State Key Basic Research Development and Programming Project (95-13-02-02). Contribution No. RCEG200308, Exploration Geophysical Center, China Earthquake Administration.     

Crustal electric structure of Haiyuan arcuate tectonic region in the northeastern margin of Qinghai-Xizang Plateau, China

Through the analysis and 2-D inversion for the 5 profiles in Haiyuan arcuate tectonic region (105°～107°E,36°～37.5°N) in the northeastern margin of Qinghai-Xizang Plateau, we have obtained the electric structure within a range of 160 km in width (east-west) and 60 km in depth in the studied area. The results show that the crustal electric structure can be divided into 6 sections, corresponding respectively to Xiji basin (Ⅰ), Xihuashan-Nanhuashan uplift (Ⅱ), Xingrenbu-Haiyuan basin (Ⅲ), Zhongwei-Qingshuihe basin (Ⅳ), Zhongning-Hongsibu basin (Ⅴ) and west-margin zone of Ordos (Ⅵ) from the southwest to the northeast. The crustal electric structure is characterized by a broom-shaped pattern, which scatters to the northwest and shrinks to the southeast. The structures in the top part of Haiyuan arcuate tectonic region are complete and large, however, they diminish from the arc top to the northwest and southeast ends. In the depth from 0 km to 10 km, the resistivity is high in the sections Ⅱ and Ⅵ, but relatively low in the other four sections, showing a similar pattern of basin depression. The electrical basement in the section Ⅲ is the deepest, displaying a "dustpan" shape that is deep in the southwest and shallow in the northeast. A series of discontinuous zones with high conductivity exist in the middle-lower crust in Haiyuan arcuate tectonic region, which is possibly related to the moderate and strong earthquakes in the region. The resistivity distribution in the focal area of the 1920 Haiyuan earthquake is significantly heterogeneous with an obviously high conductivity zone near the hypocenter regime.     

Crustal structure using receiver function in the east part of A＇nyemaqen suture belt

Twenty broadband seismographs were deployed along Hongyuan, Sichuan to Wuwei, Gansu. 81 teleseismic events were recorded in one year. We computed receiver functions from teleseismic waveform data and obtained S wave velocity structure beneath each station along the profile by using receiver function inversion method. The results revealed that the crustal structure is very complex and crustal average S wave velocity is to be on the low side. Low velocity structure generally exists in the depth range of 10～40 km in the crust between Aba arc fault and northern edge fault of Qinling earth＇s axis and it is a tectonic feature of complex geological process such as ancient A＇nyemaqen Tethys ocean from closing and side colliding to subducted plate exhumed or thrust rock slice lifted. The Moho is about 50 km depth along the profile and is slightly deeper in the south than in the north.     

Correction to the crustal thickness of the northwestern China using the data of seis-mic tomography

Introduction The gravity anomaly is an indicator of the density distribution of the underground material. Therefore the gravity anomalies have been important data used for studying the deep crustal struc-ture for a long time. Many people have made detailed researches on the regional crustal structure inverted by Bouguer anomalies. In particular some empirical formulae and practical algorithms about the crustal thickness were brought forward, and a series of results were obtained (MENG, 1996)…     

Crust and uppermost mantle structure of the Ailaoshan-Red River fault from receiver function analysis

S-wave velocity structure beneath the Ailaoshan-Red River fault was obtained from receiver functions by using teleseismic body wave records of broadband digital seismic stations. The average crustal thickness, Vp/Vs ratio and Poisson’s ratio were also estimated. The results indicate that the interface of crust and mantle beneath the Ailaoshan-Red River fault is not a sharp velocity discontinuity but a characteristic transition zone. The velocity increases relatively fast at the depth of Moho and then increases slowly in the uppermost mantle. The average crustal thickness across the fault is 36―37 km on the southwest side and 40―42 km on the northeast side, indicating that the fault cuts the crust. The relatively high Poisson’s ratio (0.26―0.28) of the crust implies a high content of mafic materials in the lower crust. Moreover, the lower crust with low velocity could be an ideal position for decoupling between the crust and upper mantle.     

Crust and upper mantle structures beneath Northeast China from receiver function studies

P-wave and S-wave receiver function analyses have been performed along a profile consisted of 27 broadband seismic stations to image the crustal and upper mantle discontinuities across Northeast China. The results show that the average Moho depth varies from about 37 km beneath the Daxing’anling orogenic belt in the west to about 33 km beneath the Songliao Basin, and to about 35 km beneath the Changbai mountain region in the east. Our results reveal that the Moho is generally flat beneath the Daxing’anling region and a remarkable Moho offset (about 4 km) exists beneath the basin-mountain boundary, the Daxing’anling-Taihang Gravity Line. Beneath the Tanlu faults zone, which seperates the Songliao Basin and Changbai region, the Moho is uplift and the crustal thickness changes rapidly. We interpret this feature as that the Tanlu faults might deeply penetrate into the upper mantle, and facilitate the mantle upwelling along the faults during the Cenozoic era. The average depth of the lithosphere-asthenosphere boundary (LAB) is ~80 km along the profile which is thinner than an average thickness of a continental lithosphere. The LAB shows an arc-like shape in the basin, with the shallowest part approximately beneath the center of the basin. The uplift LAB beneath the basin might be related to the extensive lithospheric stretching in the Mesozoic. In the mantle transition zone, a structurally complicated 660 km discontinuity with a maximum 35 km depression beneath the Changbai region is observed. The 35 km depression is roughly coincident with the location of the stagnant western pacific slab on top of the 660 km discontinuity revealed by the recent P wave tomography.     

玉树震区及邻近地区壳幔速度结构特征

对获得的玉树震区的人工地震测深资料,利用正演拟合地震测深方法进行处理与计算,获得了玉树震区的地壳速度结构和构造的基本特征。结果表明,不同区域的地壳速度结构在纵向和横向上均具有明显的非均匀性。 地壳呈层状结构,研究区结晶基底界面起伏较大,玉树附近下方较厚达8km左右,在400.0km桩号向北至温泉方向基底结构逐渐减薄为2.5km。基底界面的凹陷、隆起与不同的地质构造单元具有较好的对应关系;地壳厚度由囊谦、玉树附近的72km向南、北逐渐减薄为62km。在200.0～400.0km桩号之间二维速度等值线及壳内界面起伏变化较大,在玉树附近下方莫霍面有一弧形凹陷。     

大别造山带上部地壳结构的有限差分层析成像

将有限差分层析反演方法用于大别造山带人工地震测深剖面上部地壳初至波的走时层析成像，得到这一剖面上部地壳横向不均匀结构图像表明，在南大别地区的超高压变质带发现榴辉岩相岩石，其下方3km深度以内的基底仍具有大陆地壳构造中正常的结晶基底速度（6．00km／s左右）；北大别地区在整个上地壳保持正常的速度，而在超高压带下方3km深度以下为6．20-630km／s的相对高速异常区，这一现象可能与超高压变质岩的含量增大有关．同时还表明，南大别和北大别之间至少在上地壳已有明显差异，它们之间的水吼一五河断裂可能是大别造山带内部的主要构造分界线．     

Study of crustal thickness and poisson’s ratio of the south of Erenhot area by teleseismic receiver function

The Xing’an Mongolian Orogenic Belt (XMOB) and the northern margin of North China Craton (NCC) have undergone multistage tectonic superimposition and the tectonic evolution is extremely complicated. We collect the teleseismic data of 44 temporary broadband seismic stations deployed in the XMOB and the northern margin of NCC to calculate the P wave receiver functions. The crustal thickness and average crustal ratio as well as the Poisson’s ratios beneath 33 stations are estimated using the H-κ stacking method. The results show: (1) the crustal thickness of the study area ranges from 38.7 to 42.7 km, with an average thickness of 41.2 km. There is a great difference in crustal thickness on both sides of Solonker suture zone. The characteristics of crustal thickness support the geodynamic model that the Paleo-Asian Ocean subducted and closed at the Solonker suture zone. (2) The Poisson’s ratios in the study area are low, ranging from 0.215 to 0.277, with an average value of 0.243, suggesting that the rock composition of the area is dominated by felsic-acid rocks. (3) There exists a negative correlation between the Poisson’s ratio and the crustal thickness. Combined with the lower values of Poisson’s ratio, we speculate that the delamination is the major mechanism in crustal extension and thinning in the study area.     

Crustal structure beneath the Dabie orogenic belt from ambient noise tomography

The Qinling–Dabie–Sulu orogenic belt in east-central China is the largest high and ultrahigh pressure (HP and UHP) metamorphic zone in the world. The Dabie Mountains are the central segment of this orogenic belt between the North China and Yangtze cratons. This work studies the nature of the crustal structure beneath the Dabie orogenic belt to better understand the orogeny. To do that, we apply ambient noise tomography to the Dabie orogenic belt using ambient noise data from 40 stations of the China National Seismic Network (CNSN) between January 2008 and December 2009. We retrieve high signal noise ratio (SNR) Rayleigh waves by cross-correlating ambient noise data between most of the station pairs and then extract phase velocity dispersion measurements from those cross-correlations using a spectral method. Taking those dispersion measurements, we obtain high-resolution phase velocity maps at 8–35 second periods. By inverting Rayleigh wave phase velocity maps, we construct a high-resolution 3D shear velocity model of the crust in the Dabie orogenic belt.The resulting 3D model reveals interesting crustal features related to the orogeny. High shear wave velocities are imaged beneath the HP/UHP metaphoric zones at depths shallower than 9 km, suggesting that HP/UHP metaphoric rocks are primarily concentrated in the upper crust. Underlying the high velocity HP/UHP metamorphic zones, low shear velocities are observed in the middle crust, probably representing ductile shear zones and/or brittle fracture zones developed during the exhumation of the HP/UHP metamorphic rocks. Strong high velocities are present beneath the Northern Dabie complex unit in the middle crust, possibly related to cooling and crystallization of intrusive igneous rocks in the middle crust resulting from the post-collisional lithosphere delamination and subsequent magmatism. A north-dipping Moho is revealed in the eastern Dabie with the deepest Moho appearing beneath the Northern Dabie complex unit, consistent with the model of Triassic northward subduction of the Yangtze Craton beneath the North China Craton.     

Insight into the Crustal Structure of the Eastern Marmara Region, NW Turkey

In order to investigate crustal structure beneath the eastern Marmara region, a seismic refraction survey was conducted across the North Anatolian Fault (NAF) zone in north west Turkey. Two reversed profiles across two strands of the NAF zone were recorded in the Armutlu Highland where a tectonically active region was formed by different continents. We used land explosions in boreholes and quarry blasts as seismic sources. A reliable crustal velocity and depth model is obtained from the inversion of first arrival travel times. The velocity-depth model will improve the positioning of the earthquake activities in this active portion of the NAF. A high velocity anomaly (5.6–5.8 km s⁻¹) in the central highland of Armutlu block and the low velocity (4.90 km s⁻¹) pattern north of Iznik Lake are the two dominant features. The crustal thickness is about 26 ± 2 km in the north and increases to about 32 ± 2 km beneath the central Armutlu block in the south. P-wave velocities are about 3.95 km s⁻¹ to 4.70 km s⁻¹ for the depth range between about 1 km and 5 km in the upper crust. The eastern Marmara region has different units of upper crust with velocities varying with depth to almost 8 km. The high upper crust velocities are associated with Armutlu metamorphic rocks, while the low velocity anomalies are due to unconsolidated sedimentary sequences. The western side of Armutlu block has complex tectonics and is well known for geothermal sources. If these sources are continuous throughout the portions of the crust, it may be associated with a granitic intrusion and deformation along the NAF zone. That is, the geothermal sources associated with the low velocity may be due to the occurrence of widespread shear heating, even shear melting. The presence of shear melting may indicate the presence of crustal fluid imposed by two blocks of the NAF system.