Geodynamic model of the evolution of the Pacific Ocean

The present Pacific Ocean differs significantly in its structure and evolution from the expanding Atlantic Ocean. The Pacific is asymmetric. Its mid-ocean ridge is located not along its median line but is closer to South America and adjoins North America. The Pacific is surrounded by a ring of subduction zones but has marginal seas only at its Eurasian margins. After the breakup of Pangea, the Atlantic began to open and the Pacific began to close. This paper examines the evolution of the Pacific Ocean and, in particular, the formation mechanisms of its present structures. Numerical modeling of the long-term drift of a large continent is performed, with the initial position of the continent corresponding to the state after the breakup of the supercontinent. At first the continent, driven by the nearest descending mantle flow, begins to approach a subduction zone. Since the mantle flows beneath a large continent have different directions, its velocity is a few times lower than that of the mantle flows near the subduction zone. As a result, a zone of extension arises at the active continental margin and a fragment is broken off from the continent; this fragment rapidly moves away and stops above the descending mantle flow as in a trap. A marginal sea forms at the active continental margin. The continent continues its slow movement toward the subduction zone. The oceanic lithosphere, which earlier sank vertically, begins to descend obliquely. This evolutionary stage corresponds to the present position of Eurasia. The modeling shows how the interaction of the continent with the mantle causes the subduction zone to roll back toward the ocean. Subsequently, the continent nevertheless catches up with the subduction zone, and they move together for a while. The marginal sea then closes and high compressive stresses arise at the active continental margin. This state corresponds to the present position of South America. During the subsequent drift, the continent together with the subduction zone reaches the mid-ocean ridge and partially overrides it. This state corresponds to North America, which was the first to break off from Pangea and passed through the stages of both Eurasia and South America. The large and slowly moving Eurasia, which formed only at the time of Pangea, is still in the first evolutionary stage of the Pacific Ocean closure.     

Horizontal stresses in the mantle and in the moving continent for the model of two-dimensional convection with varying viscosity

Numerical experiments on studying the spatial fields and evolution of viscous overlithostatic horizontal stresses and pressure in the mantle and in the moving continent are carried out. The continent moves consistently with time-dependent forces, which act from the viscous mantle. By introducing the varying viscosity, we gain the possibility for taking into account the oceanic lithosphere and the difference between the viscosity of the upper and the lower mantle in the context of a purely viscous model. The typical overlithostatic horizontal stresses in the main part of the mantle are ±(7–9) MPa (70–90 bar); in the highly viscous regions and, particularly, in the subduction zones they are at least three times larger. The descending mantle flows in the depth interval from approximately 50 km to about 300 km are more sharply pronounced in the pressure field than in the field of horizontal stresses. At the considered stages of motion and in different parts, the continent is characterized by the following typical values of stresses: the overlithostatic pressure ranges from ?5 to +15 MPa; the horizontal overlithostatic tensile stress amounts up to ?4MPa (?40 bar); and the compressive stress in case of the overriding of the subduction zone attains +35 MPa (350 bar).     

Subduction of continental margins and the uplift of high-pressure metamorphic rocks

The mechanism by which high-pressure metamorphosed continental material is emplaced at high structural levels is a major unsolved problem of collisional orogenesis. We suggest that the emplacement results from partial subduction of the continental margin which, because of its high flexural rigidity, produces a rapid change in the trajectory of the descending slab. We assume a two-fold increase in effective elastic thickness of the lithosphere as the continental margin approaches the subduction zone, and calculate the flexural profile of a thin plate for progressive downward migration of the zone of increased rigidity. We assess the effect of changes in the flexural profile on the overlying accretionary prism and mantle wedge as the continent approaches by estimating the extra stresses that are imposed on the wedge due to the bending moment exerted by the continental part of the plate. The wedges overlying the subduction zones, and the subducting slab itself, experience substantial extra compressional stress at depths of around 100 km, and extensional stress at shallower depths, as the continental margin passes through the zone of maximum curvature. The magnitudes of such extra stresses are probably adequate to effect significant deformation of the wedge and/or the descending plate, and are experienced in a time interval of less than 5 m.y. for typical subduction rates. The spatial variation of yield stresses in the region of the wedge and descending slab indicates that much of this deformation may be taken up in the crustal part of the descending slab, which is the weakest region in the deeper parts of the subduction zone. This may result in rapid upward migration of the crust of the partially subducted continental margin, against the flow of subduction. High-pressure metamorphosed terranes emplaced by the mechanism envisaged in this paper would be bounded by thrust faults below and normal faults above. Movement on the faults would have been coeval, and would have resulted in rapid unroofing of the high-pressure terranes, synchronous with arrival of the continental margin at the subduction zone and, therefore, relatively early in the history of a collisional orogen.     

大陆岩石圈、地幔底部异常体与地幔对流相互作用的数值模拟

在地球表层存在着占地表面积约30%的具有低固有密度、高黏度的大陆岩石圈.由于其特殊的物理化学性质,大陆岩石圈通常不直接参与下方的地幔对流,但其与地幔对流格局有着重要的相互影响.大量研究显示,在中太平洋和非洲的下地幔底部,存在着两块占核幔边界（CMB）面积约20%的高密度热化学异常体（由于其剪切波速度较低,常称作低剪切波速度省（LSVPs））.LSVPs的演化既受地幔对流的影响,同时也影响地幔物质运动的格局和动力学过程.本文系统研究了存在大陆岩石圈,下地幔LSVPs的地幔对流模型.模拟结果显示：（1）当大陆体积较小时,其边缘常伴随着俯冲,大陆区域地幔常处于下涌状态,其上地幔温度较低,大陆岩石圈在水平方向处于压应力状态.随着大陆体积的增大,大陆边缘的俯冲逐渐减弱,大陆区域地幔由下涌转为上涌,其上地幔温度较高,大陆岩石圈水平方向处于拉应力状态.（2） 岩石圈与软流圈边界（LAB）在大陆下方较深,温度较低;在海洋区域较浅,温度较高.随着大陆体积的增大,陆洋之间LAB深度、温度的差异逐渐减小.（3）大陆区域地幔底部LSVPs物质的丰度与大陆的体积呈正相关.当大陆体积较小时,大陆下方的LSVPs丰度比海洋区域少.随着大陆体积的增大,大陆下方LSVPs的丰度逐渐增大.（4）海洋地区地表热流高,且随时间波动大,大陆地区地表热流低,随时间波动较小;LSVPs区域的核幔边界热流低.     

三维板块几何形态对大陆深俯冲动力学的制约

大陆深俯冲及超高压变质作用是大陆动力学的重要研究内容,前人进行了系统的地质、地球物理观测以及数值模拟研究.然而,自然界中大陆板块的俯冲、碰撞及造山过程大部分具有明显的沿走向的差异性,这种典型的三维特征可能很大程度上依赖于会聚大陆板块的初始几何学和运动学特征.本文采用三维高分辨率的动力学数值模拟方法,建立了方形大陆板块和楔形大陆板块两种不同的俯冲-碰撞模型,并且俯冲大陆板块侧面与大洋俯冲带相邻.数值模拟结果揭示大洋板块可以持续地俯冲到地幔之中,而大陆板块俯冲到一定深度处,其前端的俯冲板块将发生断离,并进而造成残余的大陆板块俯冲角度的减小.方形大陆俯冲板块的断离深度约为150km,而楔形大陆俯冲板块的断离深度较大,约250~300km,这很大程度上取决于俯冲带中大洋板块的牵引力和大陆板块的负浮力之间的竞争关系.同时,无论方形还是楔形大陆板块俯冲模型中,板块断离后,侧向的大洋俯冲板块仍可以拖曳约60~70km宽的大陆边缘岩石圈持续向下俯冲,揭示了新西兰东部的洋-陆空间转换俯冲带的动力学机制.并且,数值模型与喜马拉雅造山带和秦岭—大别—苏鲁造山带进行了对比,进而对其高压-超高压岩石空间展布沿走向的差异性特征和机制提供了一定的启示.     

Genetic algorithm-finite element method inversion of the factors determining the recent tectonic stress field of part of East Asia area

ＧｅｎｅｔｉｃａｌｇｏｒｉｔｈｍｆｉｎｉｔｅｅｌｅｍｅｎｔｍｅｔｈｏｄｉｎｖｅｒｓｉｏｎｏｆｔｈｅｆａｃｔｏｒｓｄｅｔｅｒｍｉｎｉｎｇｔｈｅｒｅｃｅｎｔｔｅｃｔｏｎｉｃｓｔｒｅｓｆｉｅｌｄｏｆｐａｒｔｏｆＥａｓｔＡｓｉａａｒｅａＭＥＩＪＩＡＮ...     

中国东部海域及邻区岩石层地幔的P波速度结构与构造分析

利用中国大陆东部及台湾地区、日本和琉球群岛的地震观测数据,通过体波地震层析技术反演了中国东部海域及其邻近地区的P波速度结构.以此为依据分析了不同地区的岩石层性质和深部动力学条件,探讨了中朝与扬子块体、扬子与华夏块体在海区的深部边界及其构造属性,揭示出菲律宾海板块与欧亚大陆的碰撞以及板片俯冲下沉、弧后扩张作用对中国东部海域岩石层结构的影响.结果表明,中国东部海域的岩石层地幔存在明显的横向非均匀性,它们与区域构造的形成演化有一定的联系.中国大陆东部的五莲－青岛断裂与朝鲜半岛西缘断裂、济州岛南缘断裂共同构成中朝和扬子块体的边界,江绍断裂向东延伸至朝鲜半岛南端成为分隔扬子和华夏块体的边界;东海陆架与冲绳海槽的岩石层结构差异明显,东海陆架具有中国东部地区的岩石层特征,属于欧亚大陆向海域的延伸;冲绳海槽的岩石层强烈减薄,为大陆向大洋过渡的区域;沿着日本－琉球－台湾俯冲带,菲律宾海板块俯冲下沉引起的地幔扰动对中国东部海域产生了较大的影响,欧亚大陆与菲律宾海板块之间的相互碰撞导致台湾地区岩石层明显增厚.     

西太平洋俯冲带的研究及其动力学意义

讨论了西太平洋俯冲带的分布及特征、西太平洋Ｗａｄａｔｉ－Ｂｅｎｉｏｆｆ带的形态及俯冲带上的应力状态及太平洋板块、菲律宾海板块与欧亚板块之间的相互作用；总结了地震层析成像结果；计算了俯冲板块在地幔中引起的Ｐ波速度异常，提出了俯冲板块与６６０ｋｍ间断面相互作用的４种可能；研究了俯冲板块物理性质的变化、俯冲板块产生的负浮力及其影响因素；提出需要开展俯冲带对东亚大陆构造运动和演化的影响、俯冲带相互关系及演化的研究．     

基于重力数据特征的东北地区深部动力机制解析

中国东北地区处于古亚洲构造域、蒙古—鄂霍茨克构造域和环太平洋构造域叠加作用最为显著的地区,是地学研究的热点区域.为了探析欧亚大陆下西太平洋板片的俯冲形态以及揭示该区域深部地球动力学机制,利用卫星重力数据通过预处理共轭梯度快速密度反演算法获得了包含东北地区、华北部分地区及日本海海域在内的研究区域上地幔三维密度结构,结合天然地震三维层析成像结果共同揭示太平洋板片的俯冲形态和深部动力机制.俯冲的太平洋板片在日本海沟处呈高密度异常,向西横向持续扩张,深度方向上有逐渐增加趋势.不连续的高密度体俯冲至地幔转换带(410～660km)后继续水平西向俯冲,部分滞留板片可能进入下地幔;在大兴安岭断裂带下面转换带中同样发现水平分布的高密度体,推断是大兴安岭断裂带下方地幔岩石圈拆沉的结果,横向不均匀分布的俯冲板片边缘已抵至大兴安岭造山带附近,这对于研究东北地区深部动力学机制具有重要的意义.     

地幔对流对全球岩石圈应力产生与分布的作用

利用动力学模拟方法研究地幔对流对于大尺度岩石圈内部应力场形成的作用. 地幔物质内部的密度横向非均匀及表面板块运动引起地幔流动,并在岩石圈底部产生一个应力场. 该应力场作为面力将造成岩石圈本身变形,从而产生岩石圈内部的应力分布. 模拟计算结果表明,大部分俯冲带及大陆碰撞带区域应力均呈现挤压特征,如环太平洋俯冲带及印度－欧亚碰撞带等;而东太平洋洋脊、大西洋洋脊及东非裂谷处应力状态均表现为拉张;并且绝大多数热点位置处于应力拉张区域,这与目前对全球构造应力状态的理解是一致的. 计算的岩石圈内部最大水平主压应力的方向与观测表现出相当的一致,其结果总体上吻合得较好,然而在局部区域（例如西北太平洋的俯冲带、青藏高原等地区）存在着较大的差异. 研究表明,地幔对流是造成岩石圈内部大尺度应力状态及分布的一个重要因素.     

The structure of thermal convection and geotherms beneath a continent and ocean

The formation of the thermal cross section of the lithosphere and mantle upon the interaction between the mantle convection and the immobile continent surrounded by the oceanic lithosphere is studied by numerical modeling. The convective temperature and velocity fields and then the averaged geotherms for subcontinental and suboceanic regions up to the boundary with the core are calculated from the solution of convection equations with a jump in viscosity in the continental zone. Using the experimental data on the solidus temperature in the rocks of the upper mantle, the average thickness of the continental and oceanic lithosphere is estimated at 190 and 30 km, respectively. The effect of a hot spot formed in the subcontinental upper mantle at a depth of 250–500 km, which has not been previously noted, is revealed. Although the temperature in this zone is typically assumed to be close to adiabatic, the calculations show that it is actually higher than adiabatic by up to 200°C. The physical mechanism responsible for this effect is associated with the accumulation of convective heat beneath the thermally insulating layer of the continental lithosphere. The revealed anomalies can be important in studying the phase and mineral transformations at the base of the lithosphere and in the regional geodynamical reconstructions.     

The mantle convection model with non-Newtonian rheology and phase transitions: The flow structure and stress fields

The mantle convection model with phase transitions, non-Newtonian viscosity, and internal heat sources is calculated for two-dimensional (2D) Cartesian geometry. The temperature dependence of viscosity is described by the Arrhenius law with a viscosity step of 50 at the boundary between the upper and lower mantle. The viscosity in the model ranges within 4.5 orders of magnitude. The use of the non-Newtonian rheology enabled us to model the processes of softening in the zone of bending and subduction of the oceanic plates. The yield stress in the model is assumed to be 50 MPa. Based on the obtained model, the structure of the mantle flows and the spatial fields of the stresses σ_xz and σ_xx in the Earth’s mantle are studied. The model demonstrates a stepwise migration of the subduction zones and reveals the sharp changes in the stress fields depending on the stage of the slab detachment. In contrast to the previous model (Bobrov and Baranov, 2014), the self-consistent appearance of the rigid moving lithospheric plates on the surface is observed. Here, the intense flows in the upper mantle cause the drift and bending of the top segments of the slabs and the displacement of the plumes. It is established that when the upwelling plume intersects the boundary between the lower and upper mantle, it assumes a characteristic two-level structure: in the upper mantle, the ascending jet of the mantle material gets thinner, whereas its velocity increases. This effect is caused by the jump in the viscosity at the boundary and is enhanced by the effect of the endothermic phase boundary which impedes the penetration of the plume material from the lower mantle to the upper mantle. The values and distribution of the shear stresses σ_xz and superlithostatic horizontal stresses σ_xx are calculated. In the model area of the subducting slabs the stresses are 60–80 MPa, which is by about an order of magnitude higher than in the other mantle regions. The character of the stress fields in the transition region of the phase boundaries and viscosity step by the plumes and slabs is analyzed. It is established that the viscosity step and endothermic phase boundary at a depth of 660 km induce heterogeneities in the stress fields at the upper/lower mantle boundary. With the assumed model parameters, the exothermic phase transition at 410 km barely affects the stress fields. The slab regions manifest themselves in the stress fields much stronger than the plume regions. This numerically demonstrates that it is the slabs, not the plumes that are the main drivers of the convection. The plumes partly drive the convection and are partly passively involved into the convection stirred by the sinking slabs.     

Viscous mantle tangential stresses in a numerical model of a supercontinental cycle

The values of maximum tangential stresses, as well as their spatial orientations, are calculated at each point of the mantle. These calculations are performed for successive stages of the supercontinental cycle in terms of a numerical 2-D model of assembling and diverging continents driven by mantle flows; in turn, the flows themselves are formed under the thermal and mechanical action of continents. Zones of maximum tangential stresses, within which their values are equal to or exceed 30 MPa (300 bar), are shown to be located in places where ascending or descending mantle flows approach the horizontal mantle boundaries. These zones have the shape of a bowl with the broad side toward the boundary, and their maximum tangential stresses are oriented at an angle of about 45° to the horizontal. The dimensions of these zones amount approximately to 300 km laterally and 300 km in depth. Their subhorizontal and subvertical tangential stresses are close to zero. It is also shown that the stress values vary significantly during the evolution of the suboceanic mantle under a new opening ocean (after the breakup of a supercontinent). Thus, the maximum tangential stresses in this vast mantle region at the stage of fast ocean opening are, on average, more than two times higher than at the subsequent slow stage. Note that, at all stages, the axes of the maximum tangential stresses are oriented nearly horizontally or vertically almost throughout the aforementioned area.     

The upper mantle flow beneath the North China Platform

In this paper we establish an upper mantle convection model which is constrained by regional isostatic gravity anomalies. Comparing the computed convection patterns with the tectonic features of the North China Platform we find that there are two positive anomaly centers connected with upward flows. These anomalies belong to the tectonic units of the Shan-Xi geoanticline and the Lu-Xi geoanticline. The centers of downward flows are connected with the tectonic units of the Liao-Ji geosyncline. It is reasonable to suggest that the upward mantle flows push the lithosphere upward and generate the observed positive isostatic gravity anomaly. The downward mantle flows pull the lithosphere down and generate the negative anomaly. However, the use of simple analysis makes it difficult to explain the complex lithospheric dynamics of this region. In order to understand lithospheric structures and tectonic features we must investigate the mechanical properties of the lithosphere and the relationship between the lithosphere and the mantle. These problems are discussed in the last section of this paper.     

Three-dimensional numerical simulations of crustal deformation and subcontinental mantle convection

3-D simulations of mantle convection allowing for continental crust are explored to study the effects of crustal thickening on lithosphere stability and of continents on large-scale mantle flow. Simulations begin with a crustal layer within the upper thermal boundary layer of a mantle convection roll in a 1 × 1 × 1 Cartesian domain. Convective stresses cause crust to thicken above a sheet-like mantle downwelling. For mild convective vigor an initial crustal thickness variation is required to induce 3-D lithospheric instability below the zone of crustal convergence. The amplitude of the required variation decreases with increasing convective vigor. Morphologically, instability is manifest in formation of drip-like thermals that exist within the large-scale roll associated with initial crustal thickening. A strong surface signature of the drips is their ability to cause deviations from local Airy compensation of topography. After the initial thickening phase, the crustal accumulation that forms serves as a model analog to a continent. Its presence leads to mantle flow patterns distinctly different from the steady-state roll that results in its absence. Large lateral thermal gradients are generated at its edge allowing this region to be the initiation site for continued small-scale thermal instabilities. Eventually these instabilities induce a restructuring of large-scale mantle flow, with the roll pattern being replaced by a square cell. Although preliminary and idealized, the simulations do show the fluid dynamical plausibility behind the idea that significant mantle variations can be generated along the strike of a largely 2-D mountain chain by the formation of the chain itself. The ability of a model continent to cause a change in fundamental convective planform also suggests that the effects of continental crust on mantle convection may be low-order despite the seemingly trivial volume of crust relative to mantle.     

汤加—克马德克俯冲带现今非均匀应力场特征及其动力学意义

汤加—克马德克俯冲带是太平洋板块向澳大利亚板块俯冲碰撞的动力作用区,是全球俯冲带动力学研究的热点区域.本研究基于EHB地震目录,对汤加—克马德克俯冲带(18.5°S—28.5°S)区域进行平面拟合,得到该范围内俯冲带走向约为196°,倾角约为48°;利用该俯冲带研究区域内Global CMT目录,对不同位置、不同深度进行区域应力张量反演,得到汤加—克马德克俯冲带研究区内精细的应力图像.结果显示:(1)俯冲带浅部(60~300km)应力结构非均匀特征明显,主应力轴倾伏角变化多样,并且最大主压应力轴方位在24°S左右发生明显偏转,我们推测这可能与洋底构造路易斯维尔海链俯冲有关;(2)中部(300~500km)最大主压、主张应力轴由北向南逐渐发生偏转,这可能与由北向南流动的地幔流对俯冲板片产生推挤作用有关,并且这种推挤作用向南逐渐减弱;(3)深部(500~700km)最大主压应力轴沿俯冲方向分布;(4)本文的结果还发现了主俯冲带深部西侧"偏移"板片与主俯冲带应力结构不同,表明"偏移"板片与主俯冲带是分离的.     

西北太平洋俯冲带日本本州至中国东北段应力场反演

本研究基于Global CMT提供的1196个1976年11月—2017年1月MW4.6地震矩心矩张量解,对西北太平洋俯冲带日本本州至中国东北段的应力场进行反演计算,得到了从浅表到深部俯冲带应力状态的完整分布.结果显示:俯冲带浅表陆壳一侧应力场呈现水平挤压、垂向拉伸状态,洋壳一侧的应力状态则相反,即近水平拉张、近垂向压缩.沿着俯冲板片向下,应力主轴逐渐向俯冲板片轮廓靠拢,其中位于双地震层(120km深度附近)之上的部分,主张应力轴沿俯冲板片轮廓展布而又比其更为陡倾;双地震层内的应力模式同典型I型双层地震带内的应力模式一致,即上层沿俯冲板片轮廓压缩、下层沿俯冲板片轮廓拉伸;双地震层之下,应力模式逐步转变为主压应力轴平行于俯冲板片轮廓.通观所研究的整个俯冲系统,水平面内主压和主张应力轴基本保持了与西北太平洋板片俯冲方向上的一致性,同经典俯冲板片的应力导管模型所预言的俯冲带应力模式相符;而主张应力轴在俯冲板片表面之下的中源地震深度范围内转向海沟走向,或许同研究区域横跨日本海沟与千岛海沟结合带,改变的浅部海沟形态致使完整俯冲板片下部产生横向变形有关.     

唐山地震热力源的研究

本文从应力场角度分析了1976年唐山地震发生的可能力源。震前唐山地区除受有太平洋板块和印度板块的作用外,主要受有附近的热力源的作用。以渤中为中心的下辽河—渤中—黄骅上地幔隆起的高温区和沧东断裂水平温度梯度带的存在,使唐山地区位于热应力集中区。本文根据已有的观测结果,提出了将温度变化区简化为圆形、狭长矩形和椭圆形等三个力学模式,分别求得它们的解析解,并计算了它们的热应力场。结果表明:在增温区以内,正应力都是压应力,剪应力很小。在增温区以外,法向应力为压应力,周向应力为张应力,水平剪应力较大,但各个应力分量随着距离增加而衰减,只在增温区的边缘地带达到最大值,其数值有几百巴的数量级。唐山恰好位于该增温区的边缘地带,表明热力源是促使唐山地震发生的重要原因。 除唐山地震外,华北地区近年发生的几个大地震如邢台、河间、渤海、海城地震,以及1983年发生的菏泽地震,均位于热应力的高值地带,这些地震的孕育和发生显然也与热力源的作用有关。      

Continental subduction and exhumation of high-pressure rocks: insights from thermo-mechanical laboratory modelling

Thermo-mechanical physical modelling of continental subduction is performed to investigate the exhumation of deeply subducted continental crust. The model consists of two lithospheric plates made of new temperature sensitive analogue materials. The lithosphere is underlain by liquid asthenosphere. The continental lithosphere contains three layers: the weak sedimentary layer, the crust made of a stronger material, and of a still stronger lithospheric mantle. The whole model is subjected to a constant vertical thermal gradient, causing the strength reduction with depth in each lithospheric layer. Subduction is driven by both push force and pull force. During subduction, the subducting lithosphere is heating and the strength of its layers reduces. The weakening continental crust reaches maximal depth of about 120 km and cannot subduct deeper because its frontal part starts to flow up. The subducted crust undergoes complex deformation, including indicated upward ductile flow of the most deeply subducted portions and localised failure of the subducted upper crust at about 50-km depth. This failure results in the formation of the first crustal slice which rises up between the plates under the buoyancy force. This process is accompanied by the delamination of the crustal and mantle layers of the subducting lithosphere. The delamination front propagates upwards into the interplate zone resulting in the formation of two other crustal slices that also rise up between the plates. Average equivalent exhumation rate of the crustal material during delamination is about 1 cm/year. The crust-asthenosphere boundary near the interplate zone is uplifted. The subducted mantle layer then breaks off, removing the pull force and thereby stopping the delamination and increasing horizontal compression of the lithosphere. The latter produces shortening of the formed orogen and the growth of relief. The modelling reveals an interesting burial/exhumation evolution of the sedimentary cover. During initial stages of continental subduction the sediments of the continental margin are dragged to the overriding plate base and are partially accreted at the deep part of the interplate zone (at 60-70 km-depth). These sediments remain there until the beginning of delamination during which the pressure between the subducted crust and the overriding plate increases. This results in squeezing the underplated sediments out. Part of them is extruded upwards along the interplate zone to about 30-km depth at an equivalent rate of 5-10 cm/year.     

中国大陆西部及周边地区地震活动特征的研究

探讨了中国周缘板块的联合作用对中国大陆地震的控制和影响，进一步研究了中国大陆西部及周边地区的地震活动特征。结果发现，该区域的地震活动除了有高潮和低潮的轮回特征外，还有相互消长的关系，并且地震相互消长有一个特定的比例。这一特点，对于中国大陆强震活动主体区的预测，尤其是对地震高潮期的结束时间提供了一个判据。