The D″ region

Two very different types of models are currently being proposed for D″, the lowest region of the earth's mantle: (a) those in which the P and S velocities vary smoothly down to the core-mantle boundary, without any extreme change in gradient; (b) those in which the velocity gradients decrease fairly abruptly at a height of 100 km or so above the core-mantle boundary, and maintain a value close to the critical gradient down to the boundary.Type (a) is represented by model UTD124A′ of Dziewonski and Gilbert (1972) and model B1 of Jordan and Anderson (1974). Both models are in good agreement with most travel time and free oscillation data. Their validity rests on the supposition, supported in part by theoretical studies, that data which suggest the presence of a low velocity zone in D″ result from distortion of seismic waves by the core-mantle boundary.On the other hand, slowness and amplitude data from short period P waves indicate a fairly rapid decrease in velocity gradient at a depth corresponding to an epicentral distance of about 92°, and it is very unlikely that these data can be interpreted as interface phenomena. The measured P and S times at distances beyond about 96° also indicate reduced velocities in D″. The suggestion that the measured velocities are in error as a result of interface effects is weakened by the fact that the results are apparently not wavelength-dependent.Type (b) is represented by model B2 of Jordan (1972), Bolt's (1972) model, and a new model designated as ANU2. All models have high density gradients indicative of inhomogeneity in the region. Model B2 fits the oscillation data reasonably well, but has an unjustifiably low S velocity at the core-mantle boundary. In Bolt's model the P and S velocities at the top of D″ are based on the models of Herrin et al. (1968) and Jeffreys (1939), whereas in ANU2 the values are taken from Hales and Herrin (1972) and Hales and Roberts (1970b). The velocities at the core-mantle boundary in Bolt's model and ANU2 are based on observations of “diffracted” P and S. Both of these models were designed to produce flattening of the P curve at about 92°. Both may require some modification in order to be compatible with free oscillation data.     

The seismic velocity distribution in the vicinity of a mine tunnel at Thabazimbi, South Africa

Analysis of the refracted arrivals on a seismic reflection profile recorded along the wall of a tunnel at an iron mine near Thabazimbi, South Africa, shows variations in P-wave velocity in dolomite away from the de-stressed zone that vary between 4.4 and 7.2 km/s, though values greater than 5.8 km/s predominate along most of the profile. The seismic velocities at the tunnel wall, however, vary between 4.2 and 5.2 km/s. Time–depth terms are in the range from 0.1 to 0.9 ms, and yield thicknesses of the zone disturbed by the tunnel excavations of between 2 and 9 m. The very low seismic velocities away from the tunnel wall in two regions are associated with alcoves or ‘cubbies’ involving offsets in the wall of up to 10 m. The large variations in seismic velocity resolved over distances less than 15 m with signals of wavelength around 6–9 m are attributed to variations in the sizes and concentrations of fracture systems and cracks, and in the degree of groundwater saturation of the fracture systems. The results suggest that seismic velocity variations from reflection surveys may also assist modelling studies of the stress regime in deep mines, particularly if both P and S wave velocity variations can be determined. The seismic velocity variations inferred also show that application of refraction static corrections in the processing of ‘in-mine’ seismic reflection profiles is as important as in surface surveys, because of the higher frequencies of the seismic energy recorded in the deep mine environment.     

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.     

意大利埃特纳火山的三维速度结构与地震活动性

本文的研究工作,得到了埃特纳火山下面的一组新的三维速度模型。作者们用１９８０年以来地方性地震在永久和临时地震台网在４个或更多台记录到的１２４９次地震的Ｐ和Ｓ波,被选来作走时反演。选择了几种判别标志与参数化办法以显示其类似的基本特征。表明在火山的东南象限在浅层,有Ｐ波的高速分怖,它同布格重力高异常有密切关系。在该区存在低的Ｖｐ／Ｖｓ比值;沿着中央火山管道,分布有高速Ｐ波和高的Ｖｐ／Ｖｓ比值;建议该处存在有稠密的侵入的岩浆体,延伸到２０千米左右的深度,该区附近有低速的Ｐ波速度。对震源的重新定位也显示出向外倾斜的易碎区,与中央火山管道延伸出一段Ｐ波高速异常区,位于中央火山口的附近。沿着中央火山口附近,有Ｐ波的低速区,它同次生的火山锥的分布有关,与部份熔融的深部岩浆库有联系。     

The amplitude ratioPcP/P and the core-mantle boundary

Summary Using the Haskell matrix formulation, theoretical reflection coefficient curves have been calculated for a multi-layered core-mantle boundary for comparison with observational data. Two cases are considered, first when the shear velocity in the core is equal to zero and second when the core has a finite rigidity. If the velocity contrast is large between the imbedded layer and the mantle, the reflection coefficient curves for the multi-layered medium are irregular in shape as compared to those for two half-spaces, representing the core and the mantle, respectively. The reflection coefficient curves show an oscillatory character if the imbedded layer is thick and has a high velocity contrast.The observational data consist of short-period vertical-component seismograph records ofP andPcP from nuclear explosions in the Aleutian chain, Nevada, Novaya Zemlya, Kazakh and Sahara. Attenuation and geometrical spreading are taken into consideration. Four different models for the quality factorQ are applied to the observational data. The data are found to be much affected by theQ-model used for the corrections.Based on proposedQ-values, a model for the core-mantle boundary is found, characterized by two low-velocity layers at the bottom of the mantle. The thicknesses are 16.10 km (outer layer) and 19.96 km (inner layer), the compressional wave velocities 12.17 km/sec and 10.94 km/sec and the shear wave velocities are 6.29 km/sec and 5.33 km/sec, respectively. A better fit to this model is found when in addition the shear velocity in the outer core is 2.20 km/sec and the density ratio at the core-mantle boundary is 1.07. In other words, the observations favour a layer of finite rigidity in the outer core rather than a fluid one.     

Inhomogeneities near the core-mantle boundary evidenced from scattered waves: A review

Statistical properties of small-scale inhomogeneities (wavelengths between 20 and 70 km) near the core-mantle boundary are inferred from scattered core waves. Observations of scattered core waves at large seismic arrays and worldwide networks indicate that the inhomogeneities have a global nature with similar characteristics. However, there may exist a few regions having markedly stronger or weaker strengths. Scattering by volumetric inhomogeneities of about 1% inP-wave velocity in the lower mantle or by about 300 m of topographic relief of the core-mantle boundary can explain the observations. At present it is not possible to rule out either of these two alternatives, or a combination of both.     

郯庐断裂带鲁苏皖段及邻区地壳速度结构

郯庐断裂带是我国东部规模最大的深断裂带.为了揭示该断裂带的深部结构,本文利用江苏、安徽、山东、上海和浙江地震台网记录的近震到时资料,对8700个地震事件重新精确定位,进而开展了多震相地震走时成像法反演地壳速度结构.通过分析郯庐断裂带鲁苏皖段及邻区三维地壳速度结构图像,发现(1)研究区内不同构造块体具有差异明显的地壳速度...     

Seismic Imaging of a Bimaterial Interface Along the Hayward Fault,CA, with Fault Zone Head Waves and Direct P Arrivals

We observe fault zone head waves (FZHW) that are generated by and propagate along a roughly 80 km section of the Hayward fault in the San Francisco Bay area. Moveout values between the arrival times of FZHW and direct P waves are used to obtain average P-wave velocity contrasts across different sections of the fault. The results are based on waveforms generated by more than 5,800 earthquakes and recorded at up to 12 stations of the Berkeley digital seismic network (BDSN) and the Northern California seismic network (NCSN). Robust identification of FZHW requires the combination of multiple techniques due to the diverse instrumentation of the BDSN and NCSN. For single-component short-period instruments, FZHW are identified by examining sets of waveforms from both sides of the fault, and finding on one (the slow) side emergent reversed-polarity arrivals before the direct P waves. For three-component broadband and strong-motion instruments, the FZHW are identified with polarization analysis that detects early arrivals from the fault direction before the regular body waves which have polarizations along the source-receiver backazimuth. The results indicate average velocity contrasts of 3–8 % along the Hayward fault, with the southwest side having faster P wave velocities in agreement with tomographic images. A systematic moveout between the FZHW and direct P waves for about a 80 km long fault section suggests a single continuous interface in the seismogenic zone over that distance. We observe some complexities near the junction with the Calaveras fault in the SE-most portion and near the city of Oakland. Regions giving rise to variable FZHW arrival times can be correlated to first order with the presence of lithological complexity such as slivers of high-velocity metamorphic serpentinized rocks and relatively distributed seismicity. The seismic velocity contrast and geological complexity have important implications for earthquake and rupture dynamics of the Hayward fault, including a statistically preferred propagation direction of earthquake ruptures to the SE.     

Three-dimensional velocity structure and seismicity of Mt. Etna Volcano, Italy

A new set of three-dimensional velocity models beneath Mt. Etna volcano is derived in the present work. We have used P- and S-wave arrivals from local earthquakes recorded at permanent and temporary seismic networks installed since 1980. A set of 1249 earthquakes recorded at more than four seismic stations was selected for traveltime inversion. The velocity models obtained by using different data selection criteria and parametrization display similar basic features, showing a high P-wave velocity at shallow depth in the SE quadrant, in close connection with a high gravimetric Bouguer anomaly. This area shares a low V_p/V_s ratio. High P-wave velocities and high V_p/V_s ratios are obtained along the central conduits, suggesting the presence of dense, intrusive magmatic bodies extending to a depth of about 20 km. The central intrusive core is surrounded by lower P-wave velocities. The relocated earthquake hypocenters also display the presence of an outward dipping brittle region, away from the central conduits, surrounding a ductile zone spatially related to the high P-wave velocity anomalies located in proximity to the central craters.     

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

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

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.     

THE 3-D VELOCITY STRUCTURE OF CRUST AND UPPERMOST MANTLE AND ITS TECTONIC IMPLICATIONS IN FUJIAN PROVINCE

It is important to detect the fine velocity structures of the crust and uppermost mantle to understand the regional tectonic evolution, earthquake generation processes, and to conduct earthquake risk assessment. The inversion of uppermost mantle velocity and Moho depth are strongly influenced by crustal velocity heterogeneity. In this study, we collected first arrivals of Pg and Pn and secondary arrivals of Pg wave from the seismograms recorded at Fujian provincial seismic network stations. New 3-D P-wave velocities were inverted by multi-phase joint inversion method in Fujian Province. Our results show that the fault zones in Fujian Province have various velocity patterns. The shallow crust is characterized by high velocity that represents mountains, while the mid-lower crust shows low velocities. The anomalous velocities are correlated closely with tectonic faults in Fujian Province. Velocity anomalies mainly show NE-trending distribution, especially in the mid-lower crust and uppermost mantle, which is consistent with the NE-trending of the regional main fault zones. Meanwhile, a part of velocity patterns show NW trending, which is related to the secondary NW-oriented faults. Such velocity distribution also shows a geological structural pattern of "zoning in east-west direction and blocking in north-south direction" in Fujian area. In the crust, a low velocity zone is found along Zhenghe-Dapu fault zone as mentioned by previous study, however our result shows the low velocity exists at depth of 20~30km in mid-lower crust. Compared with previous study, this low velocity zone is larger and deeper both in range and depth. The crustal thickness of 28~35km from our joint inversion is similar to the results from the receiver functions of previous studies. The thinnest crust(28km)is observed at offshore in the north of Quanzhou; while the thickest crust(35km)is located west of Zhangzhou near the Zhenghe-Dapu fault zone. Generally, thinner crustal thickness is found in offshore of Fujian Province, and thicker crustal thickness is in the mainland. However, we also found that crustal thickness becomes thinner along the east side of Yongan-Jinjiang Fault. The values of Pn velocities in the region vary from 7.71 to 8.26km/s. The velocity distribution of the uppermost mantle presents a large inhomogeneity, which is correlated with the distribution of the fault zone. High Pn velocity anomalies are found mainly along the west side of the Zhenghe-Dapu fault zone(F₂), and the east side of the Shaowu-Heyuan fault zone(F₁), which is strip-shaped throughout the central part of Fujian. Low Pn velocity anomalies are observed along the coast and Taiwan Straits, including the Changle-Zhaoan fault zone, the coastal fault zone, and the Fuzhou Basin. We also found a low Pn velocity anomaly zone, which extends to the coast, in the Shaowu-Heyuan fault zone at the junction of the Fujian, Guangdong and Jiangxi Provinces. In the west of Taiwan Straits, both high and low Pn velocity anomalies are observed. Our results show that the historical strong earthquakes(larger than magnitude 6.0) are mainly distributed between positive and negative anomaly zones at different depth profiles of the crust, and similar anomalies distribution also exists at the uppermost mantle, suggesting that the occurrence of strong earthquakes in the region is not only related to the anomalous crustal velocity structure, but also affected by the velocity anomaly structure from the uppermost mantle.     

The P-wave velocity structure of Deception Island,Antarctica, from two-dimensional seismic tomography

Deception Island is a volcanic island with a flooded caldera that has a complex geological setting in Bransfield Strait, Antarctica. We use P-wave arrivals recorded on land and seafloor seismometers from airgun shots within the caldera and around the island to invert for the P-wave velocity structure along two orthogonal profiles. The results show that there is a sharp increase in velocity to the north of the caldera which coincides with a regional normal fault that defines the northwestern boundary of the Bransfield Strait backarc basin. There is a low-velocity region beneath the caldera extending from the seafloor to > 4 km depth with a maximum negative anomaly of 1 km/s. Refracted arrivals are consistent with a 1.2-km-thick layer of low-velocity sediments and pyroclastites infilling the caldera. Synthetic inversions show that this layer accounts for only a small portion of the velocity anomaly, implying that there is a significant region of low velocities at greater depths. Further synthetic inversions and melt fraction calculations are consistent with, but do not require, the presence of an extensive magma chamber beneath the caldera that extends downwards from ≤ 2 km depth.     

Origin of ULVZs near the African LLSVP: Implications from their distribution and characteristics

Ultra-low velocity zones (ULVZs) provide important information on the composition and dynamics of the core-mantle boundary (CMB). However, their global distribution and characteristics are not well constrained, especially near African large low-shear velocity provinces (LLSVPs). Here, we used ScS precursor (SdS) and postcursor (ScscS) phases recorded by various seismic networks in Africa and South America to investigate the ULVZ characteristics underlying the South Atlantic Ocean. We found no evidence of ULVZs near the SE boundary of South America, but an ULVZ was found within the SW boundary of the African LLSVP, with thicknesses ranging from 11–18 km and reductions in S-wave velocities of 18%–34%. Our results, combined with the global distribution of ULVZs, suggest that thermal activity may be essential to ULVZ formation. Moreover, subducted slab and mantle flow may also play a key role, depending on the location of the ULVZs.     

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.     

SEISMIC AND GRAVITY TWO-DIMENSIONAL INHOMOGENEOUS MODEL OF THE EARTH'S CRUST (ON THE EXAMPLE OF CENTRAL KAZAKHSTAN)*

The interpretation of deep seismic sounding (DSS) data has been made on the basis of a two-dimensional inhomogeneous model. The refracted first arrivals as well as reflected and diffracted waves on the seismic records have been utilized. The seismic section was modeled in the iso-veolcity lines v(x, y) = const, taking into account the zones of diffraction associated with deep faults. Gravity observations have been used to construct a block model of the Earth's crust with vertical boundaries. It is suggested to define the base of the crust as the zone with velocities between 7.8 and 8.2 km/s. The reflecting boundaries of different length occurring in this zone can be conformal or unconformal with the iso-velocity lines near the base of the crust. As an example of our approach to the interpretation of DSS data the folded-blocky structure of the crust with horizontal inhomogeneities of velocity and density is shown in the Kzyl-Orda-Dzheskazgan profile in Central Kazakhstan.     

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 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.     

Imaging 3-D crustal P-wave velocity structure of western Yunnan with bulletin data

Western Yunnan is a region with intensive tectonic activity and serious earthquake risk. It is of significant importance to study three dimensional crustal structure of this region to understand the tectonic setting and disaster mechanism. Densification and digitalization of seismic networks in this region provides an opportunity to study the velocity structure with bulletin data. In this study, we collect P-wave data of 10 403 regional earthquakes recorded by 79 seismic stations from January 2008 to December 2010. In addition to first arrivals data (Pg with epicentral distance less than 200 km and Pn), the Pg (or P) data with epicentral distance more than 200 km are also considered as later direct arrivals in the tomographic inversion. We also compare the quantity and the quality of the seismic data before 2010 and after 2010. The test results show that adding the follow-up Pg phase can effectively improve the inversion ability of crustal imaging, and quantity and the data quality are significantly improved since 2010. The tomographic results show that: (1) The Honghe fault zone, which is the major fault systems in this region, may cut through the entire crust, and the velocity contrasts between two sides at lower crust beneath the Honghe fault are estimated at higher than 10%, while the velocity difference below Nujiang fault zone extends only in the upper crust; (2) Most of the earthquakes in the region occurred at the interface of high-velocity media and low-velocity media, i.e., the areas with high velocity gradient, which has been validated in other areas.     

2017年米林6.9级地震震源区速度结构与余震重定位

利用西藏自治区林芝地区的固定地震台站与南迦巴瓦流动测震台站在2017年11月18日至2017年11月24日记录到的430个余震的直达波走时数据反演得到了震源区的三维P波速度、S波速度结构,并利用三维速度结构对余震进行了重定位.成像结果显示,米林地震震源区在0~5km深度内存在低地震波速度异常;在5~15km深度内,存在高地震波速度异常,该高速异常致使震源区西南侧的地震波速度高于东北侧.重定位结果中,余震呈条带状以NW-SE走向展布,震源深度具有西南方向深、东北方向浅的特征.主震位于11km深度处、高地震波速异常体顶部,余震主要分布在高地震波速度与低地震波速度过渡的区域.对成像结果的分析表明,震源区浅部的低速异常具有低泊松比的特性,与富石英的沉积变质杂岩体-东久杂岩单元的岩性特征有关;深部的速度结构特征则可能反映了发震断层上盘地震波速度高,下盘地震波速度低的介质特性.余震重定位结果与成像结果联合表明:此次地震发震断层从11km深度处,东久杂岩体下方的高地震波速度异常顶部开始破裂,继而在5~15km深度内发生后续破裂,后续破裂的发生区域正处于喜马拉雅构造单元与冈底斯构造单元接触的形变区内.此外,根据地震波速度计算的泊松比反映了震源区持续的低泊松比特征,暗示此次地震与流体活动并无直接关系.