Upper mantle P wave velocity structure of the northern part of China and Mongolia

ｎｔｒｏｄｕｃｔｉｏｎＴｈｅｄｅｔｅｒｍｉｎａｔｉｏｎｏｆｆｉｎｅｒａｄｉａｌｖｅｌｏｃｉｔｙｓｔｒｕｃｔｕｒｅｏｆｕｐｐｅｒｍａｎｔｌｅｐｌａｙｓａｎｉｍｐｏｒｔａｎｔｒｏｌｅｉｎｉｎｖｅｓｔｉｇａｔｉｏｎｏｆｍａｎｔｌｅｃｏｍｐｏｓｉｔｉｏｎａ...     

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.     

Crustal structure beneath the eastern part of the Coast Ranges (Diablo Range) of central California from explosion seismic and near earthquake data

Seismic refraction and near earthquake data of the U.S. Geological Survey for central California have been compiled into record sections along profiles and interpreted in terms of crustal structure. The profiles are located northeast of the San Andreas fault of central California and run parallel to the general structures. For the explosion seismic line through the centre of the Diablo Range, an uppermost layer (Franciscan formation) with P velocities of 3.6–5.0 km s^?1 decreases in thickness towards the northwest. The lower boundaries of layers with constant velocities of 5.75 and 6.8 km s^?1 are found at almost constant depths of 12 and 21 km, respectively. Between 21 and 26 km depth a well-defined low-velocity zone appears whose velocity is estimated as ～ 5.3 km s^?1 with the aid of a hedgehog inversion and the calculation of amplitudes. This zone is underlain by a layer 3–5 km thick with a velocity of 7.6 km s^?1. The upper-mantle velocity beneath the Moho at 29–30 km depth is 8.2 km s^?1. The near earthquake profiles, located ～ 20 km southwest and parallel to the explosion seismic line, follow more or less the Hayward and Calaveras fault systems. The velocity-depth distribution derived for the earthquake data is very similar to that found beneath the Diablo Range. However, the low-velocity zone at 21–26 km depth does not seem to exist everywhere along the line. The Moho is not disturbed beneath the Calaveras, Hayward and Silver Creek faults; it rises slightly from the Diablo Range towards the southwest.     

Delineation of Subtrappean Mesozoic Sediments in Deccan Syneclise, India, Using Travel-Time Inversion of Seismic Refraction and Wide-Angle Reflection Data

2-D shallow velocity structure is derived by travel-time inversion of the first arrival seismic refraction and wide-angle reflection data along the E–W trending Narayanpur–Nandurbar and N–S Kothar–Sakri profiles, located in the Narmada–Tapti region of the Deccan syneclise. Deccan volcanic (Trap) rocks are exposed along the two profiles. Inversion of seismic data reveals two layered velocity structures above the basement along the two profiles. The first layer with a P-wave velocity of 5.15–5.25 km s^?1 and thickness varying from 0.7–1.5 km represents the Deccan Trap formation along the Narayanpur–Nandurbar profile. The Trap layer velocity ranges from 4.5 to 5.20 km s^?1 and the thickness varies from 0.95 to 1.5 km along the Kothar–Sakri profile. The second layer represents the low velocity Mesozoic sediments with a P-wave velocity of 3.5 km s^?1 and thickness ranging from about 0.70 to 1.6 km and 0.55 to 1.1 km along the E–W and N–S profiles, respectively. Presence of a low-velocity zone (LVZ) below the volcanic rocks in the study area is inferred from the travel-time ‘skip’ and amplitude decay of the first arrival refraction data together with the prominent wide-angle reflection phase immediately after the first arrivals from the Deccan Trap formation. The basement with a P-wave velocity of 5.8–6.05 km s^?1 lies at a depth ranging from 1.5 to 2.45 km along the profiles. The velocity models of the profiles are similar to each other at the intersection point. The results indicate the existence of a Mesozoic basin in the Narmada–Tapti region of the Deccan syneclise.     

Deep structure of the East European platform according to seismic data

The paper presents a review and analysis of new seismic data related to the structure of the mantle beneath the East European platform. Analysis of observations of long-range profiles revealed pronounced differences in the structure of the lower lithosphere beneath the Russian plate and the North Caspian coastal depression. The highest P-velocities found at depths around 100 km are in the range 8.4–8.5 km s^?1. Deep structure of the Baltic shield is different from the structures of both these regions. No evidence of azimuthal anisotropy in the upper mantle was found. A distribution of P-velocity in the upper mantle and in the transition zone consistent with accurate travel-time data was determined. The model involves several zones of small and large positive velocity gradients in the upper mantle, rapid increases of velocity near 400 and 640 km depths and an almost constant positive velocity gradient between the 400 and 640 km discontinuities. The depth of the 640 km discontinuity was determined from observations of waves converted from P to SV in the mantle.     

Seismic investigations of crustal structure in east-central Minnesota

Relative arrival times from 120 mine blasts recorded on the Central Minnesota Seismic Array (CMSA) out to distances of 350 km have been statistically analysed and interpreted to yield a crustal velocity-depth function for the Minnesota area. Velocities increase continuously from 5.9 km s^?1 at the surface to 7.4 km s^?1 at 40 km depth, with a steepening of the velocity gradient at about 20 km. A long-range refraction profile extending from the Mesabi Iron Range in Northern Minnesota to the seismic array, has also been recorded. A striking feature of the profile is the absence of direct P waves in the distance range 60–120 km. Our preferred explanation for this shadow zone is a fault or fracture zone near x = 60 km. Reflections within the shadow zone have been analysed to give a Moho depth of 42 km under northeastern Minnesota. Teleseismic residuals have been computed for 85 earthquakes recorded on the CMSA. No systematic dependence of the residuals upon azimuth or distance was found. The average residuals agree quantitatively with time delays predicted from earlier upper crustal refraction studies (Mooney et al., 1970).     

Crustal P-wave velocity structure in the northeastern margin of the Qinghai-Tibetan Plateau and insights into crustal deformation

The transitional area between the northeastern margin of the Qinghai-Tibetan Plateau, Ordos Block and Alxa Block,also being the northern segment of the North-South Seismic Belt, is characterized by considerably high seismicity level and high risk of strong earthquakes. In view of the special tectonic environment and deep tectonic setting in this area, this study used two seismic wide-angle reflection/refraction cross profiles for double constraining, so as to more reliably obtain the fine-scale velocity structure characteristics in both the shallow and deep crust of individual blocks and their boundaries in the study area,and further discuss the seismogenic environment in seismic zones with strong historical earthquakes. In this paper, the P-wave data from the two profiles are processed and interpreted, and two-dimensional crustal velocity structure models along the two profiles are constructed by travel time forward modeling. The results show that there are great differences in velocity structure,shape of intra-crustal interfaces and crustal thickness among different blocks sampled by the two seismic profiles. The crustal thickness along the Lanzhou-Huianbu-Yulin seismic sounding profile(L1) increases from ~43 km in the western margin of Ordos Block to ~56 km in the Qilian Block to the west. In the Ordos Block, the velocity contours vary gently, and the average velocity of the crust is about 6.30 km s-1; On the other hand, the velocity structures in the crust of the Qilian Block and the arclike tectonic zone vary dramatically, and the average crustal velocities in these areas are about 0.10 km s-1 lower than that of the Ordos Block. In addition, discontinuous low-velocity bodies(LVZ1 and LVZ2) are identified in the crust of the Qilian Block and the arc-like tectonic zone, the velocity of which is 0.10–0.20 km s-1 lower than that of the surroundings. The average crustal thickness of the Ordos Block is consistently estimated to be around 43 km along both Profile L2(Tongchuan-Huianbu-Alashan left banner seismic sounding profile) and Profile L1. In contrast to the gently varying intra-crustal interfaces and velocity contours in the Ordos Block along Profile L1, which is a typical structure characteristic of stable cratons, the crustal structure in the Ordos Block along Profile L2 exhibits rather complex variations. This indicates the presence of significant structural differences in the crust within the Ordos Block. The crustal structure of the Helan Mountain Qilian Block and the Yinchuan Basin is featured by "uplift and depression" undulations, showing the characteristics of localized compressional deformation.Moreover, there are low-velocity zones with alternative high and low velocities in the middle and lower crust beneath the Helan Mountain, where the velocity is about 0.15–0.25 km s-1 lower than that of the surrounding areas. The crustal thickness of the Alxa Block is about 49 km, and the velocity contours in the upper and middle-lower crust of the block vary significantly. The complex crustal velocity structure images along the two seismic sounding profiles L1 and L2 reveal considerable structural differences among different tectonic blocks, their coupling relationships and velocity structural features in the seismic zones where strong historical earthquakes occurred. The imaging result of this study provides fine-scale crustal structure information for further understanding the seismogenic environment and mechanism in the study area.     

Three-dimensional structure of mesoscale eddies in the western tropical Pacific as revealed by a high-resolution ocean simulation

The three-dimensional structure of mesoscale eddies in the western tropical Pacific(6°S–20°N, 120°E–150°E)is investigated using a high-resolution ocean model simulation. Eddy detection and eddy tracking algorithms are applied to simulated horizontal velocity vectors, and the anticyclonic and cyclonic eddies identified are composited to obtain their three-dimensional structures. The mean lifetime of all long-lived eddies is about 52 days, and their mean diameter is 147 km. Two typical characteristics of mesoscale eddies are revealed and possible dynamic explanations are analyzed. One typical characteristic is that surface eddies are generally separated from subthermocline eddies along the bifurcation latitude(~13°N) of the North Equatorial Current in the western tropical Pacific, which may be associated with different eddy energy sources and vertical eddy energy fluxes in subtropical and tropical gyres. Surface eddies have maximum swirl velocities of 8–9 cm s~(-1) and can extend to about 1500 m depth. Subthermocline eddies occur below 200 m, with their cores at about 400–600 m depth, and their maximum swirl velocities can reach 10 cm s~(-1). The other typical characteristic is that the meridional velocity component of the eddy is much larger than the zonal component. This characteristic might be due to more zonal eddy pairs(two eddies at the same latitude),which is also supported by the zonal wavelength(about 200 km) in the high-frequency meridional velocity component of the horizontal velocity.     

Velocity constraints for the inner core inferred from long-period PKP amplitudes

Long-period PKP amplitudes from 16 earthquakes in the distance range 110– 170° are compared with theoretical amplitudes which are derived from synthetic seismograms calculated for 56 systematic modifications of Earth model 1066B in the inner core. A suitable normalization procedure allows for the common representation of all observed amplitudes as a function of epicentral distance. Using the theoretical amplitude distributions it can be shown that the parameters of a regression line through the logarithmic and normalized amplitudes between 110 and 134° are related to the velocity and density jump at the inner-core boundary (ICB). The analysis shows that the dominant influence on the PKP amplitudes is the P-velocity jump which can be restricted to 0.64 ± 0.05 km s^?1. There exists a trade-off between the S-velocity jump and the density jump. Restricting the latter to the reasonable range 0–1.2 g cm^?3 the S-velocity jump at the ICB can be inferred to be 2.5–3.0 km s^?1. A rather strong S-velocity gradient below the ICB follows from the condition that the S-wave travel-time through the inner core agrees with that implied by free oscillation observations. This leads to central S-velocities between 3.81 and 4.15 km s^?1, assuming a parabolic velocity law.     

Studies of mantle structure of U.S.S.R. territory on long-range seismic profiles

Long-range seismic sounding carried out during the last few years on the territory of the U.S.S.R. has shown a basic inhomogeneity of the uppermost mantle, as well as evidence of regularities in the distribution of its seismic parameters. The following data were used: times and apparent velocities of P- and S-waves for investigation of mantle velocities, converted waves for seismic discontinuity model studies and wave attenuation for Q-factor estimation. Strong regularities were distinguished in the distribution of average seismic velocities for the uppermost mantle, in their dependence on the age and type of geostructure and on their position relative to the central part of the continent. Old platforms and the inner part of the continent are marked by velocities under the Mohorovi?i? discontinuity of more than 8.2–8.3 km s^?1, young platforms and outer parts of the continent by 8.0–8.2 km s^?1, and orogenic and rift zones by 7.8–8.0 km s^?1. The difference becomes more pronounced at a depth of about 100–200 km: for the old platform mantle velocities of 8.5–8.6 km s^?1 are typical; beneath the orogenic and rift areas, inversion zones with velocities less than 7.8 km s^?1 are observed.The converted waves show fine inhomogeneities of the crust and uppermost mantle, the presence of many discontinuities with positive and negative changes of velocity, and anisotropy of seismic waves in some of the layers. Wave attenuation allowed the determination of the Q-factor in the mantle. It varied from one region to another but a close relation between Q and P-wave velocity is the main cause of its variation.     

Velocity Structure of the West-Bengal Sedimentary Basin, India along the Palashi-Kandi Profile Using a Travel-time Inversion of Wide-angle Seismic data and Gravity Modeling – An Update

The 2-D shallow velocity structure along the north-south Palashi-Kandi profile in the West Bengal sedimentary basin has been updated by travel-time inversion of seismic refraction, wide-angle reflection and gravity data. A six-layer shallow model up to a depth of about 7 km has been derived. The first layer, which has an average velocity of 2.0 kms^?1, represents the alluvium deposit, which rests over the shale formation with average velocity of 3.0 kms^?1. The thin (200 m) Sylhet limestone, observed at a nearby Palashi well, remains hidden in the present data set. Hence a 200-m thin layer with a velocity of 3.7 kms^?1, corresponding to the Sylhet limestone, has been assumed to be present throughout the profile. The fourth layer with a velocity of 4.5–4.7 kms^?1 at a depth of 1.7–2.4 km represents the Rajmahal traps. The ‘skip’ phenomenon and rapid amplitude decay of first arrivals indicate a low-velocity zone (LVZ) in the study area. Using the ‘skip’ phenomena and wide-angle reflection data, identified on seismograms, the LVZ with a velocity of 4.0 kms^?1, indicating the Gondwana sediments, has been delineated below the Rajmahal traps. The next layer with a velocity 5.4–5.6 kms^?1 overlying the crystalline basement (5.8–6.25 kms^?1) may be associated with the Singhbhum group of meta volcanic rock that has been exposed in the western part of the basin. The basement lies at a variable depth of 4.9 to 6.8 km. The overall uncertainties of various velocity and boundary nodes are ± 0.15 kms^?1 and ± 0.5 km, respectively. The elevated basement feature in the north might have acted as a structural barrier for the deposition of Sylhet limestone during the Eocene epoch. The seismically derived shallow structure correctly explains the observed Bouguer gravity anomaly along the profile. The addition of reflections in the present analysis provides a stronger control on the depths and velocities of basement and overlying sedimentary formations, compared to the earlier model derived mainly by the first arrival seismic data.     

Upper mantle heterogeneity in the northern part of Eurasia

A two dimensional velocity model of the upper mantle has been compiled from a long-range seismic profile crossing the West Siberian young plate and the old Siberian platform. It revealed considerable horizontal and vertical heterogeneity of the mantle. A sharp seismic boundary at a depth of 400 km outlines the high-velocity gradient transition zone, its base lying at a depth of 650 km. Several layers with different velocities, velocity gradients and wave attenuation are distinguished in the upper mantle. They likewise differ in their inner structure. For instance, the uppermost 50–70 km of the mantle are divided into blocks with velocities from 7.9–8.1 to 8.4–8.6 km s^?1.Comparison of the travel-time curves for the Siberian long-range profile with those compiled from seismological data for Europe distinguished large-scale upper mantle inhomogeneities of the Eurasian continent and allowed for the correlation of tectonic features and geophysical fields. The velocity heterogeneity of the uppermost 50–100 km of the mantle correlates with the platform age and heat flow, i.e., the young plates of Western Europe and Western Siberia have slightly lower velocities and higher heat flows than the ancient East European and Siberian platforms. At greater depths (150–250 km) the upper mantle velocities increase from the ocean to the inner parts of the continent. The structure of the transition zone differs significantly beneath Western Europe and the other parts of Eurasia. The sharp boundary at a depth of 400 km, traced throughout the whole continent as the boundary reflecting intensive waves, transforms beneath Western Europe into a gradient zone. This transition zone feature correlates with positions of the North Atlantic-west Europe geoid and heat-flow anomalies.     

Crustal Structure Along the Lawrencepur-Astor Profile in the Northwest Himalayas

During the Pamir Himalayan project in the year 1975 seismic refraction and wide-angle reflection data were recorded along a 270 km long Lawrencepur-Astor (Sango Sar) profile in the northwest Himalayas. The profile starts in the Indus plains and crosses the Main Central Thrust (MCT), the Hazara Syntaxis, the Main Mantle Thrust (MMT) and ends to the east of Nanga Parbat. The seismic data, as published by Guerra et al. (1983), are reinterpreted using the travel-time ray inversion method of Zelt and Smith (1992) and the results of inversion are constrained in terms of parameter resolution and uncertainty estimation. The present model shows that the High Himalayan Crystallines (HHC, velocity 5.4 km s⁻¹) overlie the Indian basement (velocity 5.8–6.0 km s⁻¹). The crust consists of four layers of velocity 5.8–6.0, 6.2, 6.4 and 6.8 km s⁻¹ followed by the upper mantle velocity of 8.2 km s⁻¹ at a depth of about 60 km.     

Accurate Establishment of Error Models for the Satellite Gravity Gradiometry Recovery and Requirements Analysis for the Future GOCE Follow-On Mission

Firstly, the new single and combined error models applied to estimate the cumulative geoid height error are efficiently produced by the dominating error sources consisting of the gravity gradient of the satellite-equipped gradiometer and the orbital position of the space-borne GPS/GLONASS receiver using the power spectral principle. At degree 250, the cumulative geoid height error is 1.769 × 10^?1 m based on the new combined error model, which preferably accords with a recovery accuracy of 1.760 ×10^?1 m from the GOCE-only Earth gravity field model GO_CONS_GCF_2_TIM_R2 released in Germany. Therefore, the new combined error model of the cumulative geoid height is correct and reliable in this study. Secondly, the requirements analysis for the future GOCE Follow-On satellite system is carried out in respect of the preferred design of the matching measurement accuracy of key payloads comprising the gravity gradient and orbital position and the optimal selection of the orbital altitude of the satellite. We recommend the gravity gradient with an accuracy of 10^?13?10^?15 /s², the orbital position with a precision of 1-0.1 cm and the orbital altitude of 200-250 km in the future GOCE Follow-On mission.     

The distribution of thermal diffusivity in the Earth's crust

Data in the literature and additional measurements on the thermal diffusivities of granites, granulites and ultrabasic rocks at temperatures up to 1000 K and pressures to 2 GPa, have been used to propose a new model for thermal diffusivity distribution in the crust and upper mantle.The laboratory measurements were made using a pulse method or the Angstroem method with cylindrical heat flow. After making particular assumptions about the pressure and temperature distribution within the top 60 km the pressure and temperature dependencies of diffusivity were transformed into a depth dependence.The model is characterised by a continuous decrease of diffusivity to a depth of ～30 km where there is a small but rapid increase to a nearly constant value of 7.3 × 10^?3 cm² s^?1.     

Crustal velocity structure in the southern edge of the Central Alborz (Iran)

Inversion of local earthquake travel times and joint inversion of receiver functions and Rayleigh wave group velocity measurements were used to derive a simple model for the velocity crustal structure beneath the southern edge of the Central Alborz (Iran), including the seismically active area around the megacity of Tehran. The P and S travel times from 115 well-located earthquakes recorded by a dense local seismic network, operated from June to November 2006, were inverted to determine a 1D velocity model of the upper crust. The limited range of earthquake depths (between 2 km and 26 km) prevents us determining any velocity interfaces deeper than 25 km. The velocity of the lower crust and the depth of the Moho were found by joint inversion of receiver functions and Rayleigh wave group velocity data. The resulting P-wave velocity model comprises an upper crust with 3 km and 4 km thick sedimentary layers with P wave velocities (V_p) of ～5.4 and ～5.8 km s^?1, respectively, above 9 km and 8 km thick layers of upper crystalline crust (V_p ～6.1 and ～6.25 km s^?1 respectively). The lower crystalline crust is ～34 km thick (V_p ～ 6.40 km s^?1). The total crustal thickness beneath this part of the Central Alborz is 58 ± 2 km.     

Geophysical constraints on deep weathering and water storage potential in the Southern Sierra Critical Zone Observatory

The conversion of bedrock to regolith marks the inception of critical zone processes, but the factors that regulate it remain poorly understood. Although the thickness and degree of weathering of regolith are widely thought to be important regulators of the development of regolith and its water‐storage potential, the functional relationships between regolith properties and the processes that generate it remain poorly documented. This is due in part to the fact that regolith is difficult to characterize by direct observations over the broad scales needed for process‐based understanding of the critical zone. Here we use seismic refraction and resistivity imaging techniques to estimate variations in regolith thickness and porosity across a forested slope and swampy meadow in the Southern Sierra Critical Zone Observatory (SSCZO). Inferred seismic velocities and electrical resistivities image a weathering zone ranging in thickness from 10 to 35 m (average = 23 m) along one intensively studied transect. The inferred weathering zone consists of roughly equal thicknesses of saprolite (P‐velocity < 2 km s^?1) and moderately weathered bedrock (P‐velocity = 2–4 km s^?1). A minimum‐porosity model assuming dry pore space shows porosities as high as 50% near the surface, decreasing to near zero at the base of weathered rock. Physical properties of saprolite samples from hand augering and push cores are consistent with our rock physics model when variations in pore saturation are taken into account. Our results indicate that saprolite is a crucial reservoir of water, potentially storing an average of 3 m³ m^?2 of water along a forested slope in the headwaters of the SSCZO. When coupled with published erosion rates from cosmogenic nuclides, our geophysical estimates of weathering zone thickness imply regolith residence times on the order of 10⁵ years. Thus, soils at the surface today may integrate weathering over glacial–interglacial fluctuations in climate. Copyright © 2013 John Wiley & Sons, Ltd.     

Receiver function study in northern Sumatra and the Malaysian peninsula

In this receiver function study, we investigate the structure of the crust beneath six seismic broadband stations close to the Sunda Arc formed by subduction of the Indo-Australian under the Sunda plate. We apply three different methods to analyse receiver functions at single stations. A recently developed algorithm determines absolute shear-wave velocities from observed frequency-dependent apparent incidence angles of P waves. Using waveform inversion of receiver functions and a modified Zhu and Kanamori algorithm, properties of discontinuities such as depth, velocity contrast, and sharpness are determined. The combination of the methods leads to robust results. The approach is validated by synthetic tests. Stations located on Malaysia show high-shear-wave velocities (V _S) near the surface in the range of 3.4–3.6 km s^− 1 attributed to crystalline rocks and 3.6–4.0 km s^− 1 in the lower crust. Upper and lower crust are clearly separated, the Moho is found at normal depths of 30–34 km where it forms a sharp discontinuity at station KUM or a gradient at stations IPM and KOM. For stations close to the subduction zone (BSI, GSI and PSI) complexity within the crust is high. Near the surface low V _S of 2.6–2.9 km s^− 1 indicate sediment layers. High V _S of 4.2 km s^− 1 are found at depth greater than 6 and 2 km at BSI and PSI, respectively. There, the Moho is located at 37 and 40 km depth. At station GSI, situated closest to the trench, the subducting slab is imaged as a north-east dipping structure separated from the sediment layer by a 10 km wide gradient in V _S between 10 and 20 km depth. Within the subducting slab V _S ≈ 4.7 km s^− 1. At station BSI, the subducting slab is found at depth between 90 and 110 km dipping 20° ± 8° in approximately N 60° E. A velocity increase in similar depth is indicated at station PSI, however no evidence for a dipping layer is found.     

One-Dimensional Thermal Modeling of the Eastern Pontides Orogenic Belt (NE Turkey)

The eastern Pontides orogenic belt is one of the most complex geodynamic settings in the Alpine–Himalayan belt due to the lack of systematical geological and geophysical data. In this study, 1D crustal structure and P-wave velocity distribution obtained from gravity modeling and seismological data in the area has been used for the development of the thermal model of the eastern Pontides orogenic belt. The computed temperature-depth profiles suggest a temperature of 590?±?60°C at a Moho depth of 35?km indicates the presence of a brittle-ductile transition zone. This temperature value might be related to water in the subducted crust of the Tethys oceanic lithosphere. The Curie temperature depth value of 29?km, which may correspond to the crustal magma chambers, is found 5–7?km below the Moho depth. Surface heat flow density values vary from 66.5 and 104.7?mW?m^?2. High mantle heat flow density value of 48?mW?m^?2 is obtained for the area should be related to melting of the lithospheric mantle caused by upwelling of asthenosphere.