Geodynamical basis for crustal deformation under the Tibetan Plateau

This paper uses plate tectonics and satellite-derived gravity data to further discussion of crustal deformation under the Tibetan Plateau. The first of our three contributions is a spherical harmonic analysis of the global plate boundary system. A distribution of 470 Dirac delta functions is applied to describe the generating forces according to the rates of crustal creation and destruction on the plate boundaries. Analysis of the extensional and compressional forces in the spreading and subducting zones shows that the present global plate motion causes compressional stresses in the N-S direction under the Tibetan Plateau. The second contribution is the calculation of the crustal stresses in Tibet as inferred from satellite gravity data. By applying solutions to the problem of the spherical shells, the satellite-determined stresses indicate that the up-welling mantle material under Tibet induce N-S and E-W extension. Finally, a superimposed stress system is constructed. This stress system shows that the present crustal deformation in Tibet does not produce N-S shortening but generates E-W extension.The results of this paper have provided geodynamical explanations for geological field observations in Tibet and fault plane solutions of earthquakes in the Tibetan side of the India-Eurasia collision. The stress patterns reveal that the cold downwelling mantle convection flow beneath southern Tibet pulls the Indian plate down but applies a bending moment on the end of the plate to uplift and support the mass of the Himalayas.     

利用海洋重力场变化分析洋底板块运动

洋底板块运动是地球动力学和全球变化研究的重要内容.本文根据质量迁移与地球外部重力场变化的对应关系,利用不同时期测高资料推算的1995—2019全球海洋重力场变化结果,反演分析全球洋底板块运动特征.研究表明,板块汇聚边界、板块内无震海岭、海山群、断裂带等区域重力异常变化显著,而在板块离散边界无明显变化趋势;西南印度洋中脊、大西洋中脊、中印度洋中脊等地区重力异常垂直梯度变化显著,且在西太平洋俯冲带、部分海岭区域也存在明显变化,其空间分布与地形基本吻合.海洋重力场变化整体上准确反映了全球洋底板块构造运动.相较于重力异常变化反演结果,重力垂直梯度的变化能够更为准确地反映洋底板块运动特征,特别是在洋中脊区域,扩张速率越小,垂直重力梯度变化越显著.此外,详细讨论了测高海洋重力场不确定因素对洋底板块运动分析结果的影响,海面坡度改正是主要因素之一.     

中国大陆及邻区板内应力场的数值模拟及动力机制探讨

中国大陆位于欧亚板块的东南部, 受到印度板块、太平洋板块和菲律宾海板块的碰撞挤压与俯冲作用, 其构造应力场形态和动力学机制相当复杂. 本文采用伪三维有限元方法, 以世界应力图2008年版本数据(WSM2008)的应力方向和应力型两类指标作为主要约束, 对中国大陆及邻区的动力驱动机制进行数值模拟, 给出了中国大陆周边地区板块边界力的大小和方向估计. 同时对3个典型情况的数值模型进行了分析. 结果显示, 软流层静压推力对该区域构造应力场影响相对较小, 板块边界力作用则起主导作用; 印度板块在喜马拉雅造山带对欧亚板块的碰撞控制了中国大陆地区应力场的基本形态, 是形成川滇地区走滑型地震为主的重要原因; 琉球海沟——南海海槽俯冲带边界力显示了挤压-张性的分段特性, 贝加尔裂谷表现为拉张作用. 进一步的分析表明, 中国大陆大部分区域内最大水平剪应力分布图像与该地区地震辐射能量密度的分布存在较好的空间正相关性．      

Origin of tectonic stresses in the Chinese continent and adjacent areas

Based on data of principal stress orientation from focal mechanism and of geological features in China, we made pseudo-3D genetic algorithm finite element (GA-FEM) inversion to investigate the main forces acting on the Chinese continent and adjacent areas which form the Chinese tectonic stress field. The results confirm that plate boundary forces play the dominant role in forming the stress field in China, as noticed by many previous researchers. However, we also find that topographic spreading forces, as well as basal drag forces of the lower crust to the upper crust, make significant contribution to stresses in regional scale. Forces acting on the Chinese continent can be outlined as follows: the collision of the India plate to the NNE is the most important action, whereby forces oriented to the NW by the Philippine plate and forces oriented to the SWW by the Pacific plate are also important. Topographic spreading forces are not negligible at high topographic gradient zones, these forces are perpendicular to edges of the Tibetan Plateau and a topographic gradient belt running in the NNE direction across Eastern China. Basal drag forces applied by the ductile flow of the lower crust to the base of upper crust affect the regional stress field in the Tibetan Plateau remarkably, producing the clockwise rotation around the eastern Himalaya syntax.     

冲绳板块应力场数值模拟及其动力学特征div>

冲绳板块位于菲律宾海板块向欧亚板块俯冲形成的西太平洋边缘活动带上,构造应力场图像及其动力学机制表现得相当复杂．采用伪三维有限元方法,以WSM2008 观测应力场数据的应力取向和应力型两方面指标作为主要约束,对冲绳板块构造应力场进行了数值模拟．通过对计算结果的分析,对模型涉及的各种作用力作出了估计．在此基础上,对冲绳板块岩石层的状态,以及该地区的板块动力学特征进行了探讨,并得到了以下一些初步认识:① 软流层静压推力控制着该地区构造应力场的基本形态;② 冲绳海槽的演化过程,例如该地区的岩石层减薄与其下地幔流的上升等,也在很大程度上影响了该地区的板内应力场空间分布特征;③ 琉球俯冲带边界力的作用是分段的,不同区段作用力对板内应力场的影响有所不同．      

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

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

Hotspot distribution, gravity, mantle tomography: evidence for plumes

Thermal convection is the motor of Earth dynamics and therefore is the link between plate motions, hotspots, seismic velocity variations in the mantle, and anomalies of the gravity field. Small scale mantle anomalies, such as plumes, do, however, generally escape detection by tomographic methods. It is attempted to approach the problem of detection in a somewhat statistical manner. Correlations are sought between spherical harmonic expansions of the fields under study: the hotspot distribution, mantle velocity variations, gravity, heat flow. Using spherical harmonic representations of global fields implies integration and averaging over the whole globe. Thus, although relationships may remain masked in the space domain by a multitude of effects, tendencies may become visible in the spectra or in appropriate averages.The main results are the following: There is a significant long wavelength (n=2,3) negative correlation between the hotspot density and the P-wave velocity variation in the lower mantle. Positive hotspot density of degree 2 to 9 generally correlates with low seismic velocity in all depths of the upper mantle and with positive gravity. This fits well with plume-type convection. These results are also confirmed regionally for a number of individual mid-ocean ridges and hotspots. The hotspot density and the free air anomalies are distinctly positive above regions of low velocity extending to great depth. The effect is not distinct at ridges with shallow velocity anomalies. In a general way, we suggest that the antipodal upwellings (Pacific, Africa) are divided by downwelling currents around the shrinking Pacific. Plate boundaries can easily move away from their past connections with the deeper mantle. Small scale plume currents seem to be depicted in the hotspot expansion. © 1999 Elsevier Science Ltd. All rights reserved.     

亚太网三期GPS观测数据处理结果及地壳形变特征分析

本文对亚洲太平洋地区地球动力学计划(APRGP,简称亚太网)在1999-2001年采集的三期GPS观测数据,利用GAMIT/GLOBK软件进行了分析处理,得出各期的观测站在ITRF97框架下的三维直角坐标,并根据三期观测数据中的重复观测站进行了速度场计算和地壳形变特征分析。结果表明亚洲大陆板块具有明显的顺时针旋转运动特征;印度板块和澳大利亚板块从西南方向对亚洲大陆板块进行挤压;而西太平洋地区各观测站则向西北方向运动。这些特征与国内GPS区域网观测数据处理得出的结论是一致的,而本文则揭示了更大范围的亚洲太平洋地区地壳形变信息。     

Early Tertiary rupture of the Pacific plate: 1700 km of dextral offset along the Emperor trough-Line Islands lineament

Between 67 and ～40 Ma ago a northwest-southeast-trending fracture system over 8000 km long split the Pacific plate and accumulated at least 1700 km of dextral offset between the east and west portions. This system, here named the Emperor fracture zone (EFZ) system, consisted of several segments, one along the present trace of the Emperor trough and another along the Line Islands, joined by short spreading ridges. The EFZ terminated at its northern end against the Kula-Pacific ridge, and at its southern end in a ridge-transform system, called the Emperor spreading system, which extended to the west, north of Australia.The finite angular velocity vector describing the relative motion between the East and West Pacific plates is ～0.6°/Ma about a pole at 36°N, 70°W. This vector, added to the known Early Tertiary motion of the Pacific plate with respect to the global hotspot reference frame, accounts in large part for the NNW trend of the Emperor seamount chain relative to the WNW Hawaiian trend, without violation of the integrity of the Antarctic plate. The Meiji-Emperor and Emperor-Hawaiian bends date, respectively, the initiation (～67 Ma ago) and cessation (～40 Ma ago) of seafloor spreading on the Emperor spreading system.The postulated Early Tertiary relative motion along the EFZ between the East and West Pacific plates explains (1) the present misalignment of the two sets of magnetic bights of the Pacific, (2) the abrupt truncation of eastern Pacific bathymetric lineaments against the Emperor trough and Line Islands, (3) the contrast in paleolatitude between the eastern and western Pacific as indicated by paleomagnetic and sedimentologic studies, and (4) the anomalous gravity signature of the central Hawaiian ridge that indicates that the ridge loaded thin hot lithosphere.     

地幔深部运动,全球构造格局及中国地震应力场背景c

本文探讨了形成岩石层内部应力场两种可能的力源:地幔对流产生的作用于岩石层底部的切向拖曳力和沿板块边界分布力系对岩石层内应力场的影响.发现这两种力同时作用控制了中国大陆应力场的基本格局,它们产生的应力场主压应力方向和用地震震源机制解、钻孔应力测量及地质构造推断的中国应力场分布基本吻合.      

从克拉通破坏到板块动力模型的研究

华北克拉通破坏动力机制研究导致了全球动力系统及“板块动力模型”研究.板块运动最有可能的动力是地幔物质流动,但由于地幔物质流动的成因至今尚未查明,所以板块构造学说研究仍处于“运动模型”阶段,而没有进入“动力模型”阶段.如果地幔密度异常是驱动地幔物质流动的成因,那么就有可能基于重力学方法以“板块动力模型”的形式建立地幔密度异常驱动模式;软流圈中可能存在着动力特性不同的区块,地幔密度正异常代表物质盈余、区块内的物质要向区块外移动,地幔密度负异常代表物质亏损、区块外的物质要向区块内移动.本项目采用重力和地震资料相结合研究地球的整体分层,根据重力大地水准面联合地震波速度结构反演求解地幔密度异常,再根据地球正常密度假说和板块运动重力学机制的观点并与现有“板块运动模型”相结合,分析地幔密度异常动力区块,初步建立“基于重力学机制的板块动力模型”;为最终建立多学科机制的“全球板块动力模型”,迈出重要一步.     

带弹簧支撑的新型泥石流拦挡坝抗冲击性能研究

为了解决普通重力式泥石流拦挡坝在大规模泥石流发生时容易被大块石损毁的问题,提出带弹簧支撑的新型泥石流拦挡坝,并采用钢球代替泥石流中的大块石,运用有限元软件ANSYS/LS-DYNA对这种新型坝体进行泥石流大块石的冲击模拟,从坝体的变形、支反力、Von Mises应力、加速度和应变率5个方面与普通重力式拦挡坝进行对比分析。结果表明：与普通泥石流拦挡坝相比,带弹簧支撑的新型泥石流拦挡坝的支反力可减小3/4左右,钢球的冲击加速度可降到1/3.87,坝体的响应加速度可降到1/2.61,说明带弹簧支撑的新型泥石流拦挡坝具有非常好的抗冲击性能。     

One-dimensional modelling of mantle flow

A one-dimensional model of flow between a fixed boundary at the bottom and a moving one on top with no net flow through vertical sections is tested for geophysically interesting mantle viscosity-depth functions. Such a model, although simplistic, may help in answering the question to what depth the return flow extends, at least in the case of moving plates measuring many thousand kilometers across, such as the Pacific plate.It the viscosity in the asthenosphere is less than three orders of magnitude smaller than that of the mantle below, the return flow extends to great depth and the asthenosphere is a zone of concentrated shear. If the viscosity contrast is greater, the return flow is concentrated in the asthenosphere. For a wide range of model parameters typical flow velocities below the asthenosphere are about one-tenth of the plate velocity. The pressure gradient required by the mantle flow may be manifest in gravity trends across moving plates, but no excessive gravity anomalies are required by the model if the absolute viscosity values conform to those inferred from post-glacial rebound data. A thinner and lower-viscosity layer is favored over a thicker and more viscous layer if both fit glacial rebound evidence. The present model may not be applicable if down to the core the viscosity is as low as about 10²¹ N s m^–2 with a free-slip bottom boundary.     

Two-phase flow properties of a horizontal fracture: The effect of aperture distribution

A systematic numerical method has been presented to investigate the constitutive relationships between two-phase flow properties of horizontal fractures and aperture distributions. Based on fractal geometry, single rough-walled fractures are generated numerically by modified successive random addition (SRA) method and then aperture distributions with truncated Gaussian distribution are formed by shear displacement between lower and upper surfaces. (The truncated Gaussian distribution is used to describe aperture evolution under different normal stresses.) According to the assumption of two-dimensional porous media and local parallel plate model, invasion percolation approach is employed to model the two-phase flow displacement (imbibition) in generated horizontal fractures, in which capillary forces are dominant over viscous and gravity forces. For truncated Gaussian distributions, constitutive relationships from numerical simulation are compared to closed-form relationships and a good agreement is obtained. The simulation results indicate strong phase interference with the sum of two phase relative permeability values being less than one in the intermediate saturations. It is found that fracture properties related to residual saturations depend on spatial correlation of aperture distributions. Based on the simulation results, we proposed an empirical relationship between the fracture residual-saturation-rated parameters and the corresponding aperture distributions.     

东亚地震活动的时空分布及其与区域应力场的关系

来自板块之间相对运动的构造力可以传递到大陆地壳,从而形成了大陆内部的地震应力场.在某一地区,包括小地震在内的地震活动性的分布可以反映该地区的应力场的变化.根据这一观点,本文根据大量的震源机制解的结果以及最近500年的地震活动资料,详细地研究了东亚地区内几个地区的区域应力场的特征.其结果表明,来自太平洋板块相对欧亚板块的俯冲所形成的构造力,控制了从华北地区到南北地震带北段的应力场.本文根据小震的地震活动变化的特征,讨论了日本一部分地区由地震活动性的变化所反映的区域应力场的变化.中国西部以及印度-澳大利亚和欧亚大陆板块边界地区,最近大约100年地震活动性的同步变化表明,来自印度-澳大利亚板块和欧亚板块碰撞所产生的构造力,传递到了中国西部.印度-澳大利亚板块和欧亚大陆板块边界,以及中国西部的地震活动,现在依然处在地震活动高潮期.      

西太平洋板块向我国东北地区深部俯冲的数值模拟

本文采用依赖温度的黏度结构以及考虑海洋板块和大陆板块厚度差异等特征,以太平洋板块向欧亚板块会聚速率作为板块速度的主要约束,通过变化海沟后撤速度模型,数值模拟西太平洋板块向中国东北的俯冲过程.结果表明,要产生类似于中国东北之下低角度的板片俯冲,海沟后撤是重要条件;而上下地幔黏度的较大差异是决定俯冲板片不穿透660 km相变面的决定因素;西太平洋板块向欧亚板块的俯冲应早于70 Ma B.P.,海沟后撤速度可能小于一些地质学家估计的45 mm/a, 而且可能是分阶段变化的;速度场表明运动学模型的反过程：大陆岩石圈之下物质的不断水平向东的流动和推挤可能成为海沟后撤的力源之一,地幔物质的这种东向流动可能与印度板块挤压碰撞欧亚板块有关,沿欧亚板块东缘的扩张构造可能是太平洋-欧亚板块运动和印度-欧亚板块运动的综合效应.     

中国海—西太平洋莫霍面深度分布特征及其地质意义

中国海—西太平洋位于欧亚板块、印澳板块和太平洋板块的交汇处,构造运动剧烈,地质情况复杂,是认识板块运动、洋陆相互作用、物质交换和能量传递不可多得的窗口,而莫霍面深度对于研究壳幔结构以及深部动力过程有着重要的意义.本文使用最新的覆盖全球的重力和地形数据,收集了深地震测深、多道地震测深等剖面183条,数字化得到2982个控制点,使用带控制点的三维界面反演方法来约束反演过程,得到中国海—西太平洋莫霍面深度,由莫霍面形态分析可知大洋板块的俯冲和印澳板块与欧亚板块的碰撞对西太平洋边缘海的形成演化有着重要作用.结合地热、岩石圈厚度、地震活动等地质地球物理资料,分析得知研究区内各个海域莫霍深度和地壳性质的变化是处于不同构造演化阶段的表现.并在马里亚纳沟弧盆拟合一条重力2.5维剖面,结果表明热物质上涌导致了马里亚纳海槽处地幔密度减小,马里亚纳海槽以及帕里西维拉海盆到西马里亚纳海岭的下地壳高密度异常是由残留的岩浆岩引起的.     

Mantle convection pattern and subcrustal stress field under Asia

Gravitational field models derived from satellite tracking and surface gravity data have been used to derive the forces in the earth's mantle under Asia. Based on studies of tectonic forces from these models, a subcrustal stress field under China has been obtained. The stresses are due to mantle convection. According to the stress patterns, the east and west China blocks and five seismic zones are identified. The tensional stresses exerted by the upwelling mantle convection flows under the crust of Tibet seem to be related to the Tibetan uplift. The compressional orogenic region from the southern tip of Lake Baikal, through Tien Shan, Hindu Kush and the Himalayas to northern Burma appears to be connected with the downwelling mantle convection flows. It is found that the directions of the subcrustal stresses under China are disposed perpendicularly to the major fault systems and seismic belts. The results of stress calculations show that the crust of north China should be in compression and that stresses within it should be sufficient to form the Shansi Graben and Linfen Basin Systems and fracture the lithosphere. This gives a possible explanation of why strong earthquakes occurred in north China which is an isolated narrow region of highest seismicity far from plate boundaries. The tensional stress fields, caused by the upwelling mantle convection flows, are found to be regions of structural kinship characterized by major concentrations of mineral and metal deposits in China.     

The modern geoid and ancient plate boundaries

The non-hydrostatic geoid is dominated by three large anomalies: an area of high gravity potential in the equatorial Pacific; another stretching from Greenland through Africa to the southwest Indian Ocean; and a semi-continuous low region passing from Hudson's Bay through Siberia to India and on to Antarctica. None of these three high-amplitude (greater than 60 m) and long-wavelength anomalies corresponds to present-day plate boundaries. However, if the modern geoid is plotted over the positions of continents and plate boundaries at 125 Ma B.P. (reconstructed relative to hotspots) a strong correlation emerges. The modern geoidal low corresponds in position to the areas of subduction surrounding the Pacific 125 Ma ago. The geoidal high now centered on Africa is entirely contained within ancient Pangaea, and the equatorial Pacific high overlies the location of the spreading centers preserved in the magnetic anomalies of the central Pacific. The most plausible cause of the large geoidal undulations is lower mantle convection only weakly coupled to plate motions. The correspondence between modern geoid and ancient plate boundaries implies either that the coupling was much more intimate in the past, or that there is a lag of at least 100 Ma in response of the lower mantle to upper mantle conditions.     

Geotectonic bipolarity, evidence from the pattern of active oceanic ridges bordering the Pacific and African plates

Active oceanic ridges are part of the global system of diverging plate boundaries encircling the Earth. They represent weak zones of the lithosphere. They are isostatically equilibrated. The system as a whole is considered to be well adapted to the present field of plate driving forces. The search for regularities in the pattern of active oceanic ridges may, therefore, provide valuable information as to the large-scale characteristics of structures and processes in the Earth’s mantle. Two large belts of active oceanic ridges are envisaged: (1) The semi-circular belt bordering the Pacific plate which extends from South of Tasmania to Northwest of Vancouver Island over a length of 20,000 km. It appears to encircle a center P1 in the central Pacific region. (2) The circum-African belt bordering the African plate which extends from the Azores to the Gulf of Aden over a length of 24,000 km. It appears to encircle a center A1 in central Africa. The attempt is made to determine the position of these centers. Extent and position of the ridge systems are described by 34 fixed points. Points R01–R20 mark the circum-African ridge system, points R21-R34 the Pacific ridge system. A least-squares adjustment is used to determine the optimum position of the centers P1 and A1. Center P1 of the Pacific ridge system is located at 169.8°W/2.6°S. Center A1 of the circum-African ridge system is located at 11.6°E/2.4°N. The location error of the centers is less than 2.8°. In view of the great extent of the ridge systems, and considering the fact that the location of P1 and A1 is based on independent data sets, the nearly antipodal and equatorial position of the centers is remarkable. The newly defined centers P1 and A1 are located close to the Pacific pole P, at 170°W/0°N, and the African pole A, at 10°E/0°N. Within the limits of error the center P1 coincides with pole P, the center A1 with pole A. Originally, these poles were introduced in order to describe a fundamental hemispherical symmetry which is apparent in the evolution of the Earth’s lithosphere during the last 180 Ma. The new results confirm the unique position of poles P and A in the global tectonic framework.