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 and upper mantle structure of stable continental regions in North America and northern Europe

From an analysis of many seismic profiles across the stable continental regions of North America and northern Europe, the crustal and upper mantle velocity structure is determined. Analysis procedures include ray theory calculations and synthetic seismograms computed using reflectivity techniques. TheP wave velocity structure beneath the Canadian Shield is virtually identical to that beneath the Baltic Shield to a depth of at least 800 km. Two major layers with a total thickness of about 42 km characterize the crust of these shield regions. Features of the upper mantle of these region include velocity discontinuities at depths of about 74 km, 330 km, 430 km and 700 km. A 13 km thickP wave low velocity channel beginning at a depth of about 94 km is also present.A number of problems associated with record section interpretation are identified and a generalized approach to seismic profile analysis using many record sections is described. TheS wave velocity structure beneath the Canadian Shield is derived from constrained surface wave data. The thickness of the lithosphere beneath the Canadian and Baltic Shields is determined to be 95–100 km. The continental plate thickness may be the same as the lithospheric thickness, although available data do not exclude the possibility of the continental plate being thicker than the lithosphere.     

中国境内天山地壳上地幔结构的地震层析成像

根据横跨中国境内天山的库车—奎屯宽频带流动地震台阵和区域地震台网记录的近震和远震P波走时数据,利用地震层析成像方法重建了沿该地震台阵剖面下方400 km深度范围内地壳上地幔的P波速度结构.结果表明：沿新疆库车—奎屯剖面,天山地壳具有明显的横向分块结构,且南、北天山地壳显示了较为强烈的横向变形特征,表明塔里木地块对天山地壳具有强烈的侧向挤压作用;在塔里木和准噶尔地块上地幔顶部有厚度约60～90 km的高速异常体,塔里木—南天山下方的高速异常体产生了较为明显的弯曲变形,而准噶尔—北天山下方的高速异常体向南一直俯冲到中天山南侧边界下方300 km的深度,两者形成了不对称对冲构造;在塔里木和准噶尔地块下方150～400 km深度存在上地幔低速体,其中塔里木地块一侧的上地幔低速物质上涌到南天山地块的下方;在塔里木—南天山200～300 km深度范围的上地幔存在高速异常体,它可能是地幔热物质向上迁移过程融断的塔里木岩石圈的拆离体. 上述结果表明,塔里木地块的俯冲可能涉及整个岩石圈深度,但其前缘仅限于南天山的北缘;青藏高原隆升的远程效应可能不但驱动塔里木岩石圈向北俯冲,同时还造成天山造山带南侧上地幔物质的涌入;天山造山带上地幔广泛存在的低速异常有助于其上地幔的变形,而上地幔物质的强烈非均匀性应有助于推动天山造山带上地幔小尺度地幔对流的形成;根据研究区地壳上地幔速度结构特征推断,新近纪以来天山快速隆升的主要力源来自青藏高原快速隆升的远程效应,相对软弱的上地幔为加速天山造山带的变形和隆升创造了必要条件.     

Seismic evidence for a metastable olivine wedge in the subducting Pacific slab under Japan Sea

We apply a forward-modeling approach to high-quality arrival time data from 23 deep earthquakes greater than 400 km depth to investigate the detailed structure of the subducting Pacific slab beneath the Japan Sea. Our results show that a finger-like anomaly exists within the subducting Pacific slab below 400 km depth, which has a P-wave velocity 5% lower than the surrounding slab velocity (or 3% lower than that of the normal mantle), suggesting the existence of a metastable olivine wedge (MOW) in the slab. The MOW top and bottom depths are 400 and 560 km, respectively. The MOW is estimated to be about 50 km wide at 400 km depth and close to the slab upper boundary. At 560 km depth the MOW is located at about 25 km below the slab upper boundary. Most of the deep earthquakes are located in the MOW. Our results favor transformational faulting as the mechanism for deep earthquakes.     

Siberian traps: Hypotheses and seismology data

Siberian traps are the result of huge basalt eruptions which took place about 250 Ma ago over a vast territory of Siberia. The genesis of Siberian traps is attributed to a mantle plume with a center in the region of Iceland or beneath the central Urals in terms of their present coordinates. The eruption mechanism is associated with delamination—replacement of the mantle lithosphere by the deep magma material. The receiver function analysis of the records from the Norilsk seismic station (NRIL) allows comparing these hypotheses with the factual data on the depth structure of the region of Siberian traps. The S-wave velocity section place the seismic lithosphere/asthenosphere boundary (LAB) at a depth of 155–190 km, commensurate with the data for the other cratons. The mantle lithosphere has a high S-wave velocity characteristic of cratons (4.6–4.8 km/s instead of the typical value 4.5 km/s). The seismic boundary, which is located at a depth around 410 km beneath the continents is depressed by ~10 km in the region of the NRIL station. The phase diagram of olivine/wadsleyite transformation accounts for this depression by a 50–100°С increase in temperature. At the depths of 350–400 km, the S-wave velocity drops due to partial melting. A new reduction in the S-wave velocities is observed at a depth of 460 km. The similar anomalies (deepening of the 410-km seismic boundary and low shear wave velocity at depths of 350–400 and 460–500 km, respectively) were previously revealed in the other regions of the Meso-Cenozoic volcanism. In the case of a differently directed drift of the Siberian lithosphere and underlying mantle at depths down to 500 km, these anomalies are barely accountable. In particular, if the mantle at a depth ranging from 200 to 500 km is fixed, the anomalies should be observed at the original locations where they emerged 250 Ma ago, i.e. thousands of km from the Siberian traps. Our seismic data suggest that despite the low viscosity of the asthenosphere, the mantle drift at depths ranging from 200 to 500 km is correlated with the drift of the Siberian lithospheric plate. Furthermore, the position of the mantle plume beneath the Urals is easier to reconcile with the seismic data than its position beneath Iceland because of the Siberian traps being less remote from the Urals.     

Upper mantle structure in the ocean-continent transition zone: Kamchatka

We consider results from modeling the crustal and upper mantle velocity structure in Kamchatka by seismic tomography and compare these with gravity data and present-day tectonics. We found a well-pronounced (in the physical fields) vertical and lateral variation for the upper mantle and found that it is controlled by fault tectonics. Not only are individual lithosphere blocks moving along faults, but also parts of the Benioff zone. The East Kamchatka volcanic belt (EKVB) is confined to the asthenospheric layer (the asthenosphere lens) at a depth of 70–80 km; this lens is 10–20 km thick and seismic velocity in it is lower by 2–4%. The top of the asthenosphere lens has the shape of a dome uplift beneath the Klyuchevskoi group of volcanoes and its thickness is appreciably greater; overall, the upper mantle in this region is appreciably stratified. A low-velocity heterogeneity (asthenolith) at least 100 km thick has been identified beneath the Central Kamchatka depression; we have determined its extent in the upper mantle and how it is related to the EKVB heterogeneities. Gravity data suggest the development of a rift structure under the Sredinnyi Range volcanic belt. The Benioff zone was found to exhibit velocity inhomogeneity; the anomalous zones that have been identified within it are related to asthenosphere inhomogeneities in the continental and oceanic blocks of the mantle.     

Phase transformations in the earth's mantle and subducting slabs: Implications for their compositions, seismic velocity and density structures and dynamics

Abstract Phase transformations in model mantle compositions and those in subducting slabs have been reviewed to a depth of 800 km on the basis of recent high-pressure experimental data. Seismic velocity and density profiles in these compositions have also been calculated using these and other mineral physics data. The nature of the seismic velocity and density profiles calculated for a pyrolite composition was found to generally agree with those determined by seismic observations (e.g. PREM). The locations of the seismic discontinuities at 400 and 670 km correspond almost exactly to the depths where the transformations of the olivine component to denser phases take place. Moreover, the steep gradients in the seismic velocity/density profiles observed between these depths are qualitatively consistent with those expected from the successive transformations in the complementary pyroxene-garnet component in the pyrolite composition. Further, the calculated seismic velocity and density values agree well with those observed in the upper mantle and mantle transition region within the uncertainties attached to these calculations and observations. Pyrolite or peridotite compositions are thus most likely to represent the composition of the mantle above 670 km depth, although some degrees of chemical heterogeneity may exist in the transition region. The observed sharp discontinuous increases of seismic velocities and density at this depth may be attributed either to the phase transformation to a perovskite-bearing assemblage in pyrolite or to chemical composition changes. Density profiles in subducted slabs have been calculated along adequate geotherms assuming that the slabs are composed of the former oceanic crust underlain by a thicker harzburgitic layer. It is shown that the former oceanic crust is substantially less dense than the surrounding pyrolite mantle at depths below 670 km, while it is denser than pyrolite in the upper mantle and the transition region. The subducted former oceanic crust may be trapped in this region, forming a geochemically enriched layer at the upper mantle-lower mantle boundary. Thick and cool slabs may penetrate into the lower mantle, but the chemically derived buoyancy may result in strong deformation and formation of megalith structures around the 670 km seismic discontinuity. These structures are consistent with those detected by recent seismic tomography studies for subduction zones.     

Seismic structure of Iceland from Rayleigh wave inversions and geodynamic implications

We have constrained the shear-wave structure of crust and upper mantle beneath Iceland by analyzing fundamental mode Rayleigh waves recorded at the ICEMELT and HOTSPOT seismic stations in Iceland. The crust varies in thickness from 20 to 28 km in western and northern Iceland and from 26 to 34 km in eastern Iceland. The thickest crust of 34–40 km lies in central Iceland, roughly 100 km west to the current location of the Iceland hotspot. The crust at the hotspot is ∼32 km thick and is underlain by low shear-wave velocities of 4.0–4.1 km/s in the uppermost mantle, indicating that the Moho at the hotspot is probably a weak discontinuity. This low velocity anomaly beneath the hotspot could be associated with partial melting and hot temperature. The lithosphere in Iceland is confined above 60 km and a low velocity zone (LVZ) is imaged at depths of 60 to 120 km. Shear wave velocity in the LVZ is up to 10% lower than a global reference model, indicating the influence of the Mid-Atlantic Ridge and the hotspot in Iceland. The lowest velocities in the LVZ are found beneath the rift zones, suggesting that plume material is channeled along the Mid-Atlantic Ridge. At depths of 100 to 200 km, low velocity anomalies appear at the Tjornes fracture zone to the north of Iceland and beneath the western volcanic zone in southwestern Iceland. Interestingly, a relatively fast anomaly is imaged beneath the hotspot with its center at ∼135 km depth, which could be due to radial anisotropy associated with the strong upwelling within the plume stem or an Mg-enriched mantle residual caused by the extensive extraction of melts.     

Heterogeneities in the mantle transition zone from observations of P-to-SV converted waves

The method of detection of P-to-SV converted waves from distant earthquakes (Vinnik, 1977) was applied to sets of long-period records from a few seismograph stations in Europe and the west of North America. The results obtained suggest that the converted phases related to the major boundaries in the mantle can be reliably detected and the depths of conversion evaluated with an accuracy of a few kilometres. The depth of the olivine-spinel transition is close to 400 km and no difference between the estimates for the north of Europe and the west of North America is found. The depth of the boundary separating the upper and lower mantle is close to 640 km, which is 30 km less than in the recent Earth-reference models. Fine S velocity stratification of this transition changes laterally from a high-gradient layer 50 km thick, terminated at the bottom by a sharp discontinuity, to a gradient layer 100 km or more thick without the discontinuity. A striking anomaly of the mantle transition zone is found in the Rio Grande rift area where a well pronounced boundary is found at 510 km depth.     

An upper mantleP-wave velocity model for eastern and southern Africa

AP-wave velocity model for the upper mantle beneath eastern and southern Africa is proposed. The top 250 km of the model is characterized by relatively low velocities similar to those deduced for the upper mantle beneath the western United States of America. At greater depths, the velocities gradually change to normal mantle values.     

中国南北带地壳和上地幔的三维速度图象

本文采用作者提出的地震层析成象法得到了中国南北带地壳和上地幔的三维速度图象。通过误差和分辨分析,以及同爆炸地震测深剖面的比较,证明了成象结果的可靠性。 成象结果表明:1.南北带的地壳和上地幔存在显著的横向不均匀性,深达450km这种不均匀性还依然存在;2.地壳上部的速度图象与地表的已知地质特征明显相关:四川盆地显著低速,康滇地轴显著高速;3.中地壳在很大范围内存在低速层,其最低速度值达5.60km/s;4.在25°N-38°N和100.0°E-103.2°E的长条带内,上地幔顶部出现低速异常,异常速度值约为7.49km/s。 成象结果还清晰地勾划了各块体间的焊接边界。120km深度的速度图象表明,扬子准地台自秦岭以南以龙门-大巴和盐源-丽江台缘褶带为其西部边界;西南以哀劳山褶皱带为界;东南则以右江褶皱带内的南盘江为界。 统计表明,地震活动与南北带的速度结构相关:从20km以上的速度图象发现,大地震大都发生在高速和低速间的过渡条带上。     

华北地区上地幔温度及岩石圈厚度分布研究

上地幔温度是控制地幔流变性和动力学过程的关键参数之一.本文利用高分辨率S波地震层析成像速度结果,基于岩石温度与地震波速度的关系,研究了华北地区上地幔50~300 km深度范围内的温度分布和"热"岩石圈厚度.为了验证结果的可靠性,本文用计算的上地幔60 km深度处的温度作为底面约束,得到了相应的地表热流.计算地表热流与观测地表热流之间符合程度较好,相对误差大部分都在地表热流观测误差范围之内.通过对上地幔的温度分布进行分析,我们研究发现:(1)在上地幔浅部,温度与地表构造之间有很好的对应关系.在小于170 km的深度上,温度呈现东高西低的分布态势.温度较高的区域集中在东部的河淮盆地、渤海湾盆地、华北平原和中部陆块的交界处、西部鄂尔多斯高原北缘的银川―河套地堑以及阴山地区,同时,这些地区的岩石圈厚度也相应较薄,大约为80~100 km;(2)西部的鄂尔多斯高原是整个华北地区岩石圈地幔温度最低的地区,比东部地区低200~400 ℃,岩石圈厚度相应最厚,平均岩石圈厚度达到140~150 km,最厚处超过160 km.(3)在170 km以下的软流圈地幔部分,温度分布发生反转,西部温度高于东部,表明东、西部陆块在地质历史时期经历了不同的深部地幔动力学过程.     

重震反演中国东北地壳上地幔三维密度结构

本文利用重力和地震P波到时数据反演得到了中国东北地区地壳上地幔三维密度结构.与单一的重力或地震反演相比,重震反演一方面有效地克服了重力反演结果垂向分辨率低的问题,另一方面也提高了地震反演结果的可靠性.结果显示:中国东北地区的地壳及上地幔剩余密度异常分布与构造单元具有明显的相关性,造山带对应低密度异常,盆地对应高密度异常;区域内火山下方有明显的低密度体存在,可能是由于太平洋板块俯冲进入上地幔并部分滞留,在滞留板块深部脱水和软流圈热物质共同作用下产生了上涌岩浆,喷发后形成了火山.     

华北克拉通及东邻西太平洋活动大陆边缘地区的P波速度结构：对岩石圈减薄动力学过程的探讨

本文通过地震层析成像研究获得了华北克拉通及其东邻地区（30°N-50°N,95°E -145°E）1°×1°的P波速度扰动图像.结果显示,在西太平洋俯冲带地区,上地幔中西倾的板片状高速异常体与其上方的低速异常区构成俯冲带与上覆地幔楔的典型速度结构式样.俯冲板片高速体在约300~400 km深度范围内被低速物质充填,暗示俯冲板片可能发生了断离.在华北克拉通地区的上地幔中发现三个东倾排列的高速异常带.在此基础上,本文构建了华北克拉通及其东邻西太平洋活动大陆边缘地区的上地幔速度结构模式图,并据此探讨克拉通岩石圈减薄与西太平洋活动大陆边缘的深部动力学联系.本文认为,太平洋板片的俯冲（断离）,触发热地幔物质上涌并在上覆地幔楔中形成对流,使克拉通岩石圈受到改造（底侵与弱化）.随着俯冲板片后撤,地幔楔中的对流场以及对岩石圈改造的影响范围均随之东移,最终导致华北克拉通岩石圈自下而上、从西向东分三个阶段依次拆沉减薄.这一模式能很好地解释现今克拉通岩石圈自西向东呈台阶状减薄的深部现象.     

基于Love波相速度反演南北地震带地壳上地幔结构

收集了南北地震带区域地震台网中292个地震台站2008年1月至2011年3月期间的地震波形数据,由频时分析方法提取了Love波相速度频散曲线,经过反演得到了研究区内的Love波相速度分布.根据Love波纯路径频散,采用线性反演方法对0.25°×0.25°的网格点进行了一维S波速度结构反演,利用线性插值获取了南北地震带地区的三维S波速度结构.结果显示了松潘—甘孜地体和川滇菱形块体地区的下地壳具有明显的S波低速层分布,该异常分布特征支持解释青藏高原隆升及其地壳物质运移的下地壳流模型.在100至120km深度上,川滇菱形块体西北部呈现较强的S波高速异常,这可能是印度岩石圈板块沿喜马拉雅东构造结下插至该区域所致,该区域下地壳的低速软弱物质与上地幔的高速强硬物质形成了鲜明对比,暗示了地壳和上地幔可能具有不同的构造运动和变形方式,这为该区域的壳幔动力学解耦提供了条件.     

华南地区上地幔P波三维速度结构和动力学意义:来自有限频层析成像的证据

利用155个宽频带流动地震仪记录的连续地震波形数据,通过有限频层析成像技术,反演获得了华南地区上地幔的高分辨率P波三维速度模型.结果显示,大致以江南造山带为界,研究区域南部的华夏块体的大部分区域上地幔存在一个清晰的低速异常构造,而研究区北部的扬子克拉通的大部分区域上地幔则存在高速异常结构,并且这些速度异常体都向下延伸到地幔转换带.一个重要的结果是在(27°N,118°E)处观测到通过410-km界面的上涌流,并且在上升的过程中逐渐向西和向北扩展,显示为华夏块体深部200～400 km深度的大范围低速异常,可能为华夏块体广泛分布的新生代岩浆活动提供深部来源.更重要的是华夏块体通过410-km界面上涌流在上涌的过程中向北延伸,越过江南造山带"侵入"到扬子克拉通的南部地区,造成了扬子克拉通较厚的岩石圈对应的高速异常体向南倾斜的假象.最后,位于117°E(郯庐断裂的南端)以东的扬子克拉通岩石圈已经被"活化",即被来自南部热的软流圈物质替而代之.同时,推断在华夏块体下方地幔转换带内低速异常体可能是与海南地幔柱有关.海南地幔柱和(27°N,118°E)410-km界面上涌流的关系还有待于今后更大范...     

青藏高原东北缘远震P波走时层析成像研究

利用青藏高原东北缘地区固定地震台网2010年4月至2015年3月期间记录的远震事件,采用多道波形互相关方法（Multi-Channel Cross-Correlation）拾取了10697个有效P波相对走时残差数据,进而采用FMTT （Fast Marching Teleseismic Tomography）方法获得了青藏高原东北缘上地幔400 km深度范围内的P波速度结构,结果显示：秦祁地块下面存在深达70 km的高速异常,阻断了青藏高原块体中下地壳低速层向东北方向的延伸;40~140 km深度范围内,四川西南部存在一个低速区,该低速区穿过龙门山断裂带进入到四川盆地内部;祁连山造山带东部低速异常区从地壳一直延伸到上地幔400 km处,表明这里可能存在一个上地幔到地壳间的热流通道;松潘-甘孜地块分布大面积的低速异常区,而鄂尔多斯地块西南缘相对速度较高,这与青藏高原为软块体、介质密度低和鄂尔多斯块体为硬块体、介质密度高相吻合.     

Seismic image of the Mount Spurr magmatic system

 The three-dimensional P-wave velocity structure of Mount Spurr is determined to depths of 10 km by tomographic inversion of 3,754 first-arriving P-wave times from local earthquakes recorded by a permanent network of 11 seismographs. Results show a prominent low-velocity zone extending from the surface to 3–4 km below sea level beneath the southeastern flank of Crater Peak, spatially coincident with a geothermal system. P-wave velocities in this low-velocity zone are approximately 20% slower than those in the shallow crystalline basement rocks. Beneath Crater Peak an approximately 3-km-wide zone of relative low velocities correlates with a near-vertical band of seismicity, suggestive of a magmatic conduit. No large low-velocity zone indicative of a magma chamber occurs within the upper 10 km of the crust. These observations are consistent with petrologic and geochemical studies suggesting that Crater Peak magmas originate in the lower crust or upper mantle and have a short residence time in the shallow crust. Earthquakes relocated using the three-dimensional velocity structure correlate well with surface geology and other geophysical observations; thus, they provide additional constraints on the kinematics of the Mount Spurr magmatic system. Received: 4 December 1997 / Accepted: 27 February 1998     

The application of synthetic seismograms to the interpretation of the upper mantle P-wave velocity structure in northern Australia

During the period 1974 to 1977, a long range seismic refraction project was conducted in Central Australia, along a profile extending south from Darwin. Earthquakes from the Banda Sea region were used as seismic sources for this experiment. An analysis by Hales and co-workers of the resulting data based on travel times, and using geometric raytracing techniques, has resulted in the construction of an upper mantle velocity model. Using synthetic seismograms to model amplitudes, it is shown that additional constraints can be placed on the derived velocity profile. The low velocity zone beneath the “200 km” discontinuity is found to have a more abrupt onset than was previously suggested. A smaller discontinuity at 325 km depth is now implied. The analysis suggests that the “400 km” discontinuity is a first order velocity increase, whereas all other observed upper mantle discontinuities are more satisfactorily modelled as second order type structures.     

根据接收函数反演得到的首都圈地壳上地幔三维S波速度结构

利用2002~2003年中国地震局地质研究所台阵实验室以唐山大震区为中心布设的40个流动宽频带地震台站和首都圈数字台网的33个宽频带台站的远震数据,采用接收函数非线性反演方法得到其中72个宽频带台站下方60 km深度范围内的S波速度结构.根据得到的各台站下方地壳上地幔的S波速度结构,并综合刘启元等（1997）用接收函数非线性反演方法得到的延怀盆地15个宽频带流动台站下方的地壳上地幔S波速度结构模型,给出了39°N～41°N,114°E~119.5°E区域内沿不同走向、不同深度S波速度分布.由于综合了利用首都圈数字地震台网的宽频带台站以及流动地震台阵的观测数据,本文给出了较前人同类研究空间分辨率更好的结果.结果表明: (1)研究区的速度结构,特别是怀来以东的速度结构十分复杂.在10~20 km深度范围内,研究区地壳具有高速和低速异常块体的交错结构.研究区中上地壳速度结构主要被与张渤地震带大体重合的NW向高速条带和穿越唐山大震区的NE向高速条带所控制,而其中下地壳的速度结构主要为延怀—三河—唐山地区上地幔隆起所控制.(2)研究区内存在若干壳内S波低速体,它们主要分布在唐山,三河及延怀盆地等地区.在这些地区,壳内低速体伴随着壳幔界面的隆起和上地幔顶部速度结构的横向变化.(3)地表断层分布与地壳速度结构分区有较好的相关性,表明断层对不同块体有明显的控制作用.其中,宝坻断裂,香河断裂和唐山断裂均为超壳断裂.(4)首都圈内大地震的分布与壳内低速体及上地幔顶部的速度结构有密切关系.对于唐山大地震的成因,仅考虑板块作用引起的水平应力场是不够的,有必要充分重视由于上地幔变形引起的地壳垂直变形和上地幔物质侵入造成的热效应.