Crustal structure of the Peninsular India

Summary Earthquake parameters for the forty aftershocks of the main Koyna earthquake of 10 December, 1967, have been determined. Depths of the foci of the earthquakes have been found to vary between 2 to 17 km. The velocities for the phasesP _g ,P ^*,P _n have been observed to be 5.78±0.00, 6.58±0.04, 8.19±0.02 km/sec, and forS _g ,S ^*,S _n to be 3.42±0.00, 3.92±0.01 and 4.62±0.01 km/sec respectively. A two-layered crustal model has been interpreted for the Peninsular shield with the average thickness of the granitic layer as 20 km and that for the basaltic layer as 18.7 km. A plot of the epicenters suggests a NNE to SSW orientation of the fault.     

不同地区人工爆炸与天然地震记录特征及识别研究

本文分析了河北怀来多次爆炸、河北三河采石场多次爆炸和低震级天然地震事件的记录特征和时频差异。结果显示:河北怀来爆炸的P波能量强、衰减快、S波发育弱;河北三河采石场爆炸的P波、S波主频均低于怀来爆炸,S波与面波混淆,不同震中距的台站记录低频发育明显;而天然地震的有效频带更宽,频率成分更为复杂。将Pg/Sg谱比判据应用于小震级地震与爆炸的识别中,探索交叉频带谱比对不同地区爆炸的识别。结果表明:高频（＞5 Hz）Pg/Sg谱比判据可将研究数据中的爆炸与小震级地震完全区分;与Sg低频（0—2 Hz）有关的交叉频带谱比可对两个不同地区的爆炸进行识别,交叉频带的谱比判据较传统的单一频带谱比判据能够更好地反映出不同类型事件的特征差异。      

A crust-mantle model for the NORSAR area

Summary Elastic waves from explosions were recorded at NORSAR and at a number of field stations, and the data were used for determining a crust-mantle model under the array. The number of explosions was eleven distributed on seven shot points. The total number of recording points was fifty-one, and the interpretation was based on 350 individual records.The velocities obtained for the crustal phases were 6.2, 6.6 and 8.2 km/sec for theP _g ,P _g andP _n waves respectively. A deep crustal phase with a velocity of about 7.4 km/sec was observed. The mean depths to the discontinuities within the crust were determined to be 17 and 26 km. The depth to Moho varied greatly across the array from 31.5 km in the central part to 38 km under the C-ring. The maximum dip observed for the Moho was 12^o.Contribution No. 57 to Norwegian Geotraverse Project.     

北京、天津、唐山和张家口地区的地壳结构

本文主要叙述利用近年来北京及其附近地区八次工业爆破的观测资料,同时也利用了一部份天然地震的资料来研究北京-天津-唐山-张家口地区地壳结构的结果,采用了工业爆破中的58个观测点和天然地震的54个观测点的数据,得到了本区的三层地壳模型和地壳厚度。 地壳的平均厚度分布各地区不很一致;北京及其东南地区为35±1.5公里,北部地区为37±1.0公里,西部地区为39±1.5公里,地壳介质中的纵波和横波的平均速度分别为6.23±0.03公里/秒和3.55±0.05公里/秒,上地幔顶部介质的平均速度为7.98±0.13公里/秒和4.60±0.03公里/秒。     

Crustal Architecture of the Inverted Central Lapland Rift Along the HUKKA 2007 Profile

We have studied the lateral velocity variations along a partly buried inverted paleo–rift in Central Lapland, Northern Europe with a 2D wide-angle reflection and refraction experiment, HUKKA 2007. The experiment was designed to use seven chemical explosions from commercial and military sites as sources of seismic energy. The shots were recorded by 102 stations with an average spacing of 3.45 km. Two-dimensional crustal models of variations in P-wave velocity and Vp/Vs-ratio were calculated using the ray tracing forward modeling technique. The HUKKA 2007 experiment comprises a 455 km long profile that runs NNW–SSE parallel to the Kittilä Shear Zone, a major deformation zone hosting gold deposits in the area. The profile crosses Paleoproterozoic and reactivated Archean terranes of Central Lapland. The velocity model shows a significant difference in crustal velocity structure between the northern (distances 0–120 km) and southern parts of the profile. The difference in P-wave velocities and Vp/Vs ratio can be followed through the whole crust down to the Moho boundary indicating major tectonic boundaries. Upper crustal velocities seem to vary with the terranes/compositional differences mapped at the surface. The lower layer of the upper crust displays velocities of 6.0–6.1 km/s. Both Paleoproterozoic and Archean terranes are associated with high velocity bodies (6.30–6.35 km/s) at 100 and 200–350 km distances. The Central Lapland greenstone belt and Central Lapland Granitoid complex are associated with a 4 km-thick zone of unusually low velocities (<6.0 km/s) at distances between 120 and 220 km. We interpret the HUKKA 2007 profile to image an old, partly buried, inverted continental rift zone that has been closed and modified by younger tectonic events. It has structural features typical of rifts: inward dipping rift shoulders, undulating thickness of the middle crust, high velocity lower crust and a rather uniform crustal thickness of 48 km.     

Crustal structure of Tushka Region, Abu-Simbel, Egypt, inferred from spectral ratios of P waves of local earthquakes

The crustal structure of North Abu-Simbel area was studied using spectral ratios of short-period P waves. Three-component short period seismograms from the Masmas seismic station of the Egyptian National Seismic Network Stations were used. The Thomson-Haskell matrix formulation was applied for linearly elastic, homogeneous crustal layers. The obtained model suggests that the crust under the study region consists of a thin (0.8 km) superficial top layer with a P-wave velocity of 3.8±0.7 km/s and three distinct layers with a mean P-wave velocity of 6.6 km/s, overlaying the upper mantle with a P-wave velocity of 8.3 km/s (fixed). The results were obtained for 14 different earthquakes. The P-wave velocities of the three layers are: 5.8±0.6 km/s, 6.5±0.4 km/s and 7.2±0.3 km/s. The total depth to the Moho interface is 32±2 km. The crustal velocity model estimated using observations is relatively simple, being characterized by smooth velocity variations through the middle and lower crust and normal crustal thickness. The resultant crustal model is consistent with the model obtained from previous deep seismic soundings along the northern part of Aswan lake zone.     

A study of the crustal structure beneath the Himalayas from body waves

Summary The crustal structure beneath the Himalayas has been investigated using body wave data from near earthquakes having epicentres over the Himalayas and recorded by the observatories situated over, or very near, the foothills of the mountains. A three-layered crustal model, without the top sedimentary layer, with velocities for theP wave group in Granite I, Granite II and the Basaltic layer as 5.48, 6.00 and 6.45 and for theS wave group as 3.33, 3.56 and 3.90 km/sec respectively, has been interpreted. The upper mantle velocity for theP wave has been observed to be 8.07 km/sec and for theS wave as 4.57 km/sec. Average thickness for the Granite I layer has been computed as 22.7 km, for the Granite II layer as 16.3 km and for the Basaltic layer as 18.7 km. Crustal and sub-crustal velocities indicate a lower trend under the mountain. A thicker crust has been obtained beneath the Himalayas.     

Ray parameters of diffracted long period P and S waves and the velocities at the base of the mantle

The derivation of P and S velocities at the core-mantle boundary (CMB) from long-period diffracted waves by the use of the simple ray-theoretical formulav _CMB=r _c /p (v _CMB=velocity at the CMB;r _c =core radius;p=ray parameter) yields apparent velocity values which differ from the true velocities. Using a dominant period of about 20 sec for calculating theoretical seismograms, we found a linear relation between the apparent velocity and the average velocity in a transition zone at the base of the mantle with fixed velocity on top.The ray parameters determined from long-period earthquake data are found to be 4.540±0.035 and 8.427±0.072 sec/deg for P_diff and S_diff, respectively. These values yield apparent velocities of 13.378±0.103 for P and 7.207±0.062 km/sec for S waves. By means of the theoretical relation between apparent and average velocity and under the assumption of linear variation of velocity with depth, one can invert the apparent velocities into true CMB velocities of 13.736±0.170 and 7.320±0.124 km/sec. These results imply positive velocity gradients at the base of the mantle and hence no significant departures from adiabaticity and homogeneity.Contribution No. 211 of the Geophysical Institute, University of Karlsruhe.     

西藏高原当雄-亚东地带地壳与上地幔结构和速度分布的爆炸地震研究

本文概述了在西藏高原长达450公里左右的南北向测线上取得九次湖中水下爆炸地震记录的处理结果。通过数字处理、拟合和反演等计算,得出了该地区地壳与上地幔的成层结构和速度分布。 结果表明,该区整个沉积岩层厚约3-5公里,雅鲁藏布江以北到当雄地带,地壳巨厚达70-73公里;江南地区为68-45公里,并逐渐向南翘起。在成层地壳介质中发现下地壳中存在低速层,厚约10公里,速度为5.64公里/秒。分析认为,高原地形与巨厚地壳的形成是印度洋板块与欧亚板块碰撞以及长期挤压和内部物质运移的结果。     

Calibration of the Tibetan Plateau Using Regional Seismic Waveforms

We use the recordings from 51 earthquakes produced by a PASSCAL deployment in Tibet to develop a two-layer crustal model for the region. Starting with their ISC locations, we iteratively fit the P-arrival times to relocate the earthquakes and estimate mantle and crustal seismic parameters. An average crustal P velocity of 6.2–6.3 km/s is obtained for a crustal thickness of 65 km while the P velocity of the uppermost mantle is 8.1 km/s. The upper layer of the model is further fine-tuned by obtaining the best synthetic SH waveform match to an observed waveform for a well-located event. Green's functions from this model are then used to estimate the source parameters for those events using a grid search procedure. Average event relocation relative to the ISC locations, excluding two poorly located earthquakes, is 16 km. All but one earthquake are determined by the waveform inversion to be at depths between 5 and 15 km. This is 15 km shallower, on average, than depths reported by the ISC. The shallow seismicity cut-off depth and low crustal velocities suggest high temperatures in the lower crust. Thrust faulting source mechanisms dominate at the margins of the plateau. Within the plateau, at locations with surface elevations less than 5 km, source mechanisms are a mixture of strike-slip and thrust. Most events occurring in the high plateau where elevations are above 5 km show normal faulting. This indicates that a large portion of the plateau is under EW extension.     

On some problems of seismic crustal phase

ＯｎｓｏｍｅｐｒｏｂｌｅｍｓｏｆｓｅｉｓｍｉｃｃｒｕｓｔａｌｐｈａｓｅＨｕａｎ－ＣｈｅｎｇＧＥ（葛焕称）（ＳｅｉｓｍｏｌｏｇｉｃａｌＢｕｒｅａｕｏｆＪｉａｎｇｓｕＰｒｏｖｉｎｃｅ,Ｎａｎｊｉｎｇ２１００１４,Ｃｈｉｎａ）Ａｂｓｔｒａｃｔ：Ｉｎｔｈｉｓ...     

Crustal seismic structure and depth distribution of earthquakes in the Archean Kuusamo region,Fennoscandian Shield

Two-dimensional crustal velocity models are derived from passive seismic observations for the Archean Karelian bedrock of north-eastern Finland. In addition, an updated Moho depth map is constructed by integrating the results of this study with previous data sets. The structural models image a typical three-layer Archean crust, with thickness varying between 40 and 52 km. P wave velocities within the 12–20 km thick upper crust range from 6.1 to 6.4 km/s. The relatively high velocities are related to layered mafic intrusive and volcanic rocks. The middle crust is a fairly homogeneous layer associated with velocities of 6.5–6.8 km/s. The boundary between middle and lower crust is located at depths between 28 and 38 km. The thickness of the lower crust increases from 5–15 km in the Archean part to 15–22 km in the Archean–Proterozoic transition zone. In the lower crust and uppermost mantle, P wave velocities vary between 6.9–7.3 km/s and 7.9–8.2 km/s. The average V_p/V_s ratio increases from 1.71 in the upper crust to 1.76 in the lower crust.The crust attains its maximum thickness in the south-east, where the Archean crust is both over- and underthrust by the Proterozoic crust. A crustal depression bulging out from that zone to the N–NE towards Kuusamo is linked to a collision between major Archean blocks. Further north, crustal thickening under the Salla and Kittilä greenstone belts is tentatively associated with a NW–SE-oriented collision zone or major shear zone. Elevated Moho beneath the Pudasjärvi block is primarily explained with rift-related extension and crustal thinning at ∼2.4–2.1 Ga.The new crustal velocity models and synthetic waveform modelling are used to outline the thickness of the seismogenic layer beneath the temporary Kuusamo seismic network. Lack of seismic activity within the mafic high-velocity body in the uppermost 8 km of crust and relative abundance of mid-crustal, i.e., 14–30 km deep earthquakes are characteristic features of the Kuusamo seismicity. The upper limit of seismicity is attributed to the excess of strong mafic material in the uppermost crust. Comparison with the rheological profiles of the lithosphere, calculated at nearby locations, indicates that the base of the seismogenic layer correlates best with the onset of brittle to ductile transition at about 30 km depth.We found no evidence on microearthquake activity in the lower crust beneath the Archean Karelian craton. However, a data set of relatively well-constrained events extracted from the regional earthquake catalogue implies a deeper cut-off depth for earthquakes in the Norrbotten tectonic province of northern Sweden.     

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.     

Shallow seismicity, stress distribution and crustal deformation pattern in the Andaman-West Sunda arc and Andaman Sea, northeastern Indian Ocean

Shallow seismicity and available source mechanisms in the Andaman–westSunda arc and Andaman sea region suggest distinct variation in stressdistribution pattern both along and across the arc in the overriding plate.Seismotectonic regionalisation indicates that the region could be dividedinto eight broad seismogenic sources of relatively homogeneousdeformation. Crustal deformation rates have been determined for each oneof these sources based on the summation of moment tensors. The analysisshowed that the entire fore arc region is dominated by compressive stresseswith compression in a mean direction of N23^°, and the rates ofseismic deformation velocities in this belt decrease northward from 5.2± 0.65 mm/yr near Nias island off Sumatra and 1.12 ±0.13 mm/yr near Great Nicobar islands to as much as 0.4 ±0.04 mm/yr north of 8^°N along Andaman–Nicobar islandsregion. The deformation velocities indicate, extension of 0.83 ±0.05 mm/yr along N343^° and compression of 0.19 ±0.01 mm/yr along N73^° in the Andaman back arc spreadingregion, extension of 0.18 ± 0.01 mm/yr along N125^° andcompression of 0.16 ± 0.01 mm/yr along N35^° in NicobarDeep and west Andaman fault zone, compression of 0.84 ±0.12 mm/yr N341^° and extension of 0.77 ± 0.11 mm/yralong N72^° within the transverse tectonic zone in the Andamantrench, N-S compression of 3.19 ± 0.29 mm/yr and an E-Wextension of 1.24 ± 0.11 mm/yr in the Semangko fault zone ofnorth Sumatra. The vertical deformation suggests crustal thinning in theAndaman sea and crustal thickening in the fore arc and Semangko faultzones. The apparent stresses calculated for all major events range between0.1–10 bars and the values increase with increasing seismic moment.However, the apparent stress estimates neither indicate any significantvariation with faulting type nor display any variation across the arc, incontrast to the general observation that the fore arc thrust events showhigher stress levels in the shallow subduction zones. It is inferred that theoblique plate convergence, partial subduction of 90^°E Ridge innorth below the Andaman trench and the active back arc spreading are themain contributing factors for the observed stress field within the overridingplate in this region.     

Structure of the crust in the Black Sea and adjoining regions from surface wave data

Group velocities of Rayleigh and Love waves along the paths across the Black Sea and partly Asia Minor and the Balkan Peninsula are used to estimate lateral variations of the crustal structure in the region. As a first step, lateral variations of group velocities for periods in the range 10–20 s are determined using a 2D tomography method. Since the paths are oriented predominantly in NE–SW or N–S direction, the resolution is estimated as a function of azimuth. The local dispersion curves are actually averaged over the extended areas stretched in the predominant direction of the paths. The size of the averaging area in the direction of the best resolution is approximately 200 km. As a second step, the local averaged dispersion curves are inverted to vertical sections of S-wave velocities. Since the dispersion curves in the 10–20 s period range are mostly affected by the upper crustal structure, the velocities are estimated to a depth of approximately 25 km. Velocity sections along 43° N latitude are determined separately from Rayleigh and Love wave data. It is shown that the crust under the sea contains a low-velocity sedimentary layer of 2–3 km thickness, localized in the eastern and western deeps, as found earlier from DSS data. Beneath the sedimentary layer, two layers are present with velocity values lying between those of granite and consolidated sediments. Velocities in these layers are slightly lower in the deeps, and the boundaries of the layers are lowered. S-wave velocities obtained from Love wave data are found to be larger than those from Rayleigh wave data, the difference being most pronounced in the basaltic layer. If this difference is attributed to anisotropy, the anisotropy coefficient = (SH - SV)/Smean is reasonable (2–3%) in the upper layers, and exceeds 9% in the basaltic layer.     

2-D Crustal Velocity Structure Along Hirapur-Mandla Profile in Central India: An Update

The 2-D crustal velocity model along the Hirapur-Mandla DSS profile across the Narmada-Son lineament in central India (Murty et al., 1998) has been updated based on the analysis of some short and discontinuous seismic wide-angle reflection phases. Three layers, with seismic velocities of 6.5–6.7, 6.35–6.40 and 6.8 km s^–1, and upper boundaries located approximately at 8, 17 and 22 km depth respectively, have been identified between the basement (velocity 5.9 km s^–1) and the uppermost mantle (velocity 7.8 km s^–1). The layer with 6.5–6.7 km s^–1 velocity is thin and is less than 2-km deep between the Narmada north (at Katangi) and south (at Jabalpur) faults. The upper crust shows a horst feature between these faults, which indicates that the Narmada zone acts as a ridge between two pockets of mafic intrusion in the upper crust. The Moho boundary, at 40–44 km depth and the intra-crustal layers exhibit an upwarp suggesting that the Narmada faults have deep origins, involving deep-seated tectonics. A smaller intrusive thickness between the Narmada faults, as compared to those beyond these faults, suggests that the intrusive activities on the two sides are independent. This further suggests that the two Narmada faults may have been active at different geological times. The seismic model is constrained by 2-D gravity modeling. The gravity highs on either side of the Narmada zone are due to the effect of the high velocity/high density mafic intrusion at upper crustal level.     

唐山震区上地壳平均波速比分析

收集河北地震台网记录的唐山震区2012-2016年ML≥1.0地震,选取震源深度5-10km、波形信噪比较高的波形数据,拾取直达P波、S波震相到时,利用和达法计算平均波速比,并分析研究区平均波速比水平分布特征与时变特征,结果表明:唐山震区上地壳平均波速比约为1.71,且水平分布呈区域化特征,反映了研究区复杂的地壳结构。研究区内3次ML≥4.0地震前,波速比曲线均呈较显著的下降-低值-回升形态,表明该区平均波速比变化可能具有一定地球物理异常特征。     

Explosion seismic reflections from the Earth’s core

Recent seismological studies have presented evidence for the existence of a layer with ultra-low seismic velocities at the core-mantle boundary at ca. 2900 km depth. We report high-amplitude, high-frequency, and laterally coherent seismic arrivals from three nuclear explosions in Siberia. With recording station intervals of 15 km, the seismic phases are readily correlated and show the presence of a thin, ultra-low velocity zone in a region where it was not previously reported. The duration and complexity of the arrivals are inconsistent with a simple core-mantle boundary and require a hitherto unidentified, kilometre-scale, fine structure in the ultra-low velocity zone. The observations may be explained by a ca. 7 km thick, two-layer, ultra-low velocity zone with exceptional low velocities, which indicate the presence of high percentages of melt (>15%), in particular in the lower part of the zone. Waveform variation implies lateral change in the thickness and physical properties of the ultra-low velocity zone with a wavelength of less than 100 km.     

P-wave velocity anomalies in the earth's mantle from the Uppsala array observations

Summary Distribution of compressional-wave velocities in the mantle is determined fromdT/d measurements using the Uppsala seismograph array station (UPSAS). Short-period vertical-component seismograms from 181 events in the epicentral distance range 16°–100° have been used. The velocity distribution shows anomalous variations at depths of 750, 1500, 1800, 2300 and 2550 km. Evidence of lateral heterogeneity beneath the northern part of the Asian continent, in the depth range 1700–2300 km, is discussed. Computed travel times, based on this velocity-depth relation, are tested by an examination of travel-time residuals from the Long Shot and Milrow explosions on Amchitka, Aleutian Islands.     

Crustal thermal regime inferred from magnetic anomaly data and its relationship to seismogenic layer thickness: The Japanese islands case study

The seismogenic layer thickness correlates with surface heat flow beneath the Japanese islands. However, this correlation is shown at restricted area, where seismic activity is high. In order to overcome this spatial limitation, we used another approach to estimate the regional thermal structure in the crust beneath the Japanese islands with more uniform coverage. The bottom depths of the magnetized crust determined from the spectral analysis of residual magnetic anomalies is generally interpreted as the level of the Curie point isotherm. We applied this method to estimate the crustal thermal structure in square windows of 2.125° × 2.125°. The obtained depths ranging from 11 to 30 km with average value of 18 km, correlate with the seismogenic layer thickness. It suggests that the Curie point depth is a useful indicator of the crustal thermal structure in these regions.