NNR constraint in ITRF2000 and the new global plate motion model NNR-ITRF2000VEL

There are 54 sites employed by ITRF2000 for ITRF2000 orientation. The deficiencies are obvious. First, these sites cannot well represent the rotation rate of the earth crust because there is no selected site in five out of fourteen tectonic plates and three of fourteen plates only have one site each. Second, the total angular momentum of the crust is non-vanishing in ITRF2000, even though it is declared that No Net Rotation (NNR) with respect to NNR-NUVEL1A is imposed on ITRF2000 construction according to the documentations of ITRF2000. So the NNR condition in conventional terrestrial reference system (CTRS) realization cannot be satisfied in ITRF2000. In this paper, the criteria of site selection for estimating the Euler vectors are suggested; the Tisserand system constraint equation in ITRF construction is derived; and as a product, the global plate motions can be obtained from the ITRF2000 construction.     

关于全球板块运动模型ITRF2000VEL若干问题探讨

利用最新的国际地球自转服务（IERS）发布的国际地球参考架ITRF2000速度场，建立了一个独立于任何其它模型、安全基于现代空间大地测量实测结果的现今全球板块运动模型ITRF2000VEL，较百万年平均地质模型NNR－NUVELIA更能真正反映全球板块现今运动特征，相比ITRF96,晚接近于百万年地质模型NNR－NUVELIA，而且ITRF2000参考架在定向、原点和尺度的定义较ITRF96和ITRF97参考架有了改进，但ITRF200VEL模型存在诸多问题：总角动量|L|=0.127，即不为零，与协议参考架（CTRF）不符，存在整体旋转，并建立无旋转NNR－ITRF2000VEL模型，台站不均匀分布全球板块、部分板块界线不明确以及有的板块不满足刚性特征等，这些对ITRF2000的高精度应用和长期维持、地球自转参数的长期变化都将产生一定的影响。     

关于ITRF96参考架整体旋转性的探讨

利用国际地球自转服务IERS发布的国际地球参考架ITRF96的速度场 ,建立了一个基于现代空间大地测量实测结果的全球板块运动模型ITRF96VEL ,该模型与NNR -NUVEL1A模型有着较好的一致性 .利用该模型 ,本文对ITRF96参考架相对无整体旋转参考架是否存在一个整体旋转的问题进行了初步的探讨 ,结果表明ITRF96参考架相对无整体旋转参考架有一个整体旋转 ,旋转角速度为 0 1 61°/Ma ,旋转极指向南纬 50 5° ,东经 65 5°,这与要求ITRF96无整体旋转的定义不相符 ,这个整体旋转将会对国际地球参考架的高精度应用和长期维持 ,特别是对地球自转参数长期变化的研究产生一定的影响     

关于ITRF96参考架整体旋转性的探讨

利用国际地球自转服务IERS发布的国际地球参考架ITRF96的速度场 ,建立了一个基于现代空间大地测量实测结果的全球板块运动模型ITRF96VEL ,该模型与NNR -NUVEL1A模型有着较好的一致性 .利用该模型 ,本文对ITRF96参考架相对无整体旋转参考架是否存在一个整体旋转的问题进行了初步的探讨 ,结果表明ITRF96参考架相对无整体旋转参考架有一个整体旋转 ,旋转角速度为 0 1 61°/Ma ,旋转极指向南纬 50 5° ,东经 65 5°,这与要求ITRF96无整体旋转的定义不相符 ,这个整体旋转将会对国际地球参考架的高精度应用和长期维持 ,特别是对地球自转参数长期变化的研究产生一定的影响     

Global tectonics and the plate motion obtained from the ITRF97 station velocity vectors

~~Global tectonics and the plate motion obtained from the ITRF97 station velocity vectors@马宗晋 @任金卫 @张进~~     

Glacial isostasy and plate motion

The influence of glacial-isostatic adjustment (GIA) on the motion of tectonic plates is usually neglected. Employing a recently developed numerical approach, we examine the effect of glacial loading on the motion of the Earth’s tectonic plates where we consider an elastic lithosphere of laterally variable strength and the plates losely connected by low viscous zones. The aim of this paper is to elucidate the physical processes which control the GIA-induced horizontal motion and to assess the impact of finite plate-boundary zones. We show that the present-day motion of tectonic plates induced by GIA is at, or above, the order of accuracy of the plate motions determined by very precise GPS observations. Therefore, its contribution should be considered when interpreting the mechanism controlling plate motion.     

利用中国地壳运动观测网络研究中国大陆相对于ITRF97板块模型形变

在建立全球ITRF97板块运动模型的基础上,利用"中国地壳运动观测网络"79个GPS基本站的数据,建立我国新的地壳运动方向图和块体运动模型.通过与NNR-NUVEL1A地质模型比较认为,ITRF97板块运动模型反映了现今十几年跨度的地壳运动,在研究我国现今几年到十几年时间跨度的地壳形变时,地壳运动背景场应采用基于ITRF97实测速度场建立欧亚板块运动模型.     

基于热点参考系的板块绝对运动模型

本文在最新一代板块相对运动模型（MORVEL模型）的基础上,通过最小二乘反演观测热点方向数据建立板块绝对运动模型.经过对前人用于约束板块绝对运动模型的所有热点数据与MORVEL模型的一致性作系统的统计学检验,发现存在离群数据.相应地,提出两种系统剔除离群数据的新方法：逐一剔除法和全局搜索法.结合对总残差的卡方检验和对残差频率分布的正态检验,逐步筛选最优模型,最终得到一个基于热点参考系的板块绝对运动新模型（T87模型）.该模型能够合理拟合全球分布的87个热点方向数据,而模型预测的板块绝对运动速率比观测热点火山的迁移速率系统偏小,最大偏差达到4 cm·a^-1.这样的偏差可由观测热点速率的系统误差或地幔返回流造成的热点间运动解释.不论是哪种解释成立,本文结果表明热点方向数据能够独立、有效地定义全球热点参考系.这样定义的热点参考系可方便地应用于板块绝对运动、地幔对流及真极移的研究.

     

Reconciling strong slab pull and weak plate bending: The plate motion constraint on the strength of mantle slabs

Although subducting slabs undergo a bending deformation that resists tectonic plate motions, the magnitude of this resistance is not known because of poor constraints on slab strength. However, because slab bending slows the relatively rapid motions of oceanic plates, observed plate motions constrain the importance of bending. We estimated the slab pull force and the bending resistance globally for 207 subduction zone transects using new measurements of the bending curvature determined from slab seismicity. Predicting plate motions using a global mantle flow model, we constrain the viscosity of the bending slab to be at most ~ 300 times more viscous than the upper mantle; stronger slabs are intolerably slowed by the bending deformation. Weaker slabs, however, cannot transmit a pull force sufficient to explain rapid trenchward plate motions unless slabs stretch faster than seismically observed rates of ~ 10^− 15 s^− 1. The constrained bending viscosity (~ 2 × 10²³ Pa s) is larger than previous estimates that yielded similar or larger bending resistance (here ~ 25% of forces). This apparent discrepancy occurs because slabs bend more gently than previously thought, with an average radius of curvature of 390 km that permits subduction of strong slabs. This gentle bending may ultimately permit plate tectonics on Earth.     

弹性板块运动模型研究进展

传统板块构造理论认为板块是一个刚体,实际上板块是可变形的.板块内部几年到几十年时间尺度的变形主要是弹性变形,因此应当用弹性模型描述板块运动.推导了板块的弹性运动方程,由空间大地测量新的观测成果建立了菲律宾海、太平洋和澳大利亚板块的弹性运动模型.发现三个板块内部都存在明显的水平形变.板内应变场的空间变化有明显的规律：板块边界附近的应变率最大,从边界向内部逐渐减小;在板块扩散边界附近,主张应变率大于主压应变率,主张应变轴基本上与边界的扩张方向一致;在俯冲边界附近,主压应变率大于主张应变率,主压应变轴基本上与板块的俯冲方向一致;在走滑兼有俯冲性质的边界附近,最大剪应变的方向与边界断裂的走向基本一致.由GPS观测得到的主压应变轴与由震源机制解得到的主压应力轴方向具有很好的一致性.板内的应力－应变场基本上遵循广义胡克定律.     

采用剪切波分裂方向约束的板块绝对运动模型研究

板块绝对运动（即岩石圈板块相对于深部地幔的运动）导致软流圈深度存在较强的剪切作用,为软流圈呈现较强的地震波速各向异性提供了一种物理解释.相应地,软流圈地震各向异性的实测数据为反演板块绝对运动提供了一种定量的约束.本文利用前人发表的由474个剪切波分裂数据组成的全球软流圈地震各向异性方向数据集,结合板块相对运动模型MORVEL,通过加权最小二乘法反演板块绝对运动.计算结果表明,由实测数据约束的板块绝对运动模型不能将观测数据拟合到原始数据测量误差要求的统计水平上,反映出以地震各向异性方向指示板块绝对运动方向存在原始数据测量误差之外的系统性误差.该误差或可归因于板块绝对运动控制软流圈地震各向异性的物理机制复杂性以及小尺度地幔流动的区域复杂性.因此,地震各向异性方向只能从统计平均的意义上约束板块绝对运动方向,两者间的差别除了测量误差外还应包含复杂物理成因的模型误差.在原始数据测量误差的基础上增加20°的模型误差并且剔除一个离群数据后,本文得到了最优拟合剪切波分裂数据集的SKS473模型,并与由热点数据反演得到的板块绝对运动模型进行对比.由于约束板块绝对运动的地震各向异性数据和热点数据都存在误差较大、地理分布不均的局限性,结合两类数据的联合反演或可成为未来建立更高精度板块绝对运动模型的有效途径.

     

板块运动与地震各向异性及应力场的相关性统计分析

利用收集到的各种来源共计7 959组的地震各向异性观测数据和21 750组应力场数据,结合板块绝对运动模型计算给出的各板块的运动规律,分别统计分析了板块运动与地震各向异性及应力场的相关性,并对板块运动对地震各向异性及应力场特征产生的影响进行了分析． 统计结果表明,阿拉伯、 加勒比、 胡安德富卡、 北美、 纳兹卡、 太平洋和南美板块上地震各向异性与板块运动均具有较好的相关性,而非洲、 南极洲、 澳大利亚、 欧亚、 印度和菲律宾板块上二者的相关性则相对较差. 讨论分析发现,板块运动拖动软流圈流动、 橄榄岩晶格优选方位、 化石各向异性和地幔流动或岩石圈流动等因素均在一定程度上控制并影响着地震各向异性与板块运动的一致性． 而板块基底拖曳力、 洋脊推力、 浮力作用和碰撞及俯冲作用等多种因素共同制约了板块运动与应力场的相关性,使得非洲、 可可斯、 欧亚、 胡安德富卡、 北美、 纳兹卡、 菲律宾和南美板块上二者的相关性较好,其它板块上其相关性则较差． 对于俯冲带地区,由于俯冲机制的复杂性和软流圈、 岩石圈地幔流动方向的不确定性,其板块运动与地震各向异性及应力场的相关性图像表现复杂,需要结合具体的俯冲带构造进行近一步研究．     

Determination of Euler parameters of Philippine Sea plate and the inferences

Euler vectors of 12 plates, including Philippine Sea plate (PH), relative to a randomly fixed Pacific plate(PA) were determined by inverting the 1122 data from NUVEL-1 global plate motion model, earthquake slip vectors along Philippine Sea plate boundary, and GPS observed velocities. Euler vectors of Philippine Sea plate relative to adjacent plates are also gained. Our results are well consistent with observed data and can satisfy the geological and geophysical constraints along the Caroline(CR)-PH and PA-CR boundaries. Deformation of Philippine Sea plate is also discussed by using the plate motion Euler parameters.     

Vertical coherency function model of spatial ground motion

Many studies have focused on horizontal ground motion, resulting in many coherency functions for horizontal ground motion while neglecting related problems arising from vertical ground motion. However, seismic events have demonstrated that the vertical components of ground motion sometimes govern the ultimate failure of structures. In this paper, a vertical coherency function model of spatial ground motion is proposed based on the Hao model and SMART 1 array records, and the validity of the model is demonstrated. The vertical coherency function model of spatial ground motion is also compared with the horizontal coherency function model, indicating that neither model exhibits isotropic characteristics. The value of the vertical coherency function has little correlation with that of the horizontal coherency function. However, the coherence of the vertical ground motion between a pair of stations decreases with their projection distance and the frequency of the ground motion. When the projection distance in the wave direction is greater than 800 meters, the coherency between the two points can be neglected.     

一种新的全球对流层天顶延迟模型GZTD

对流层延迟是GNSS导航定位主要误差源之一,主要受气象参数(如总气压、温度和水汽压等)的影响,具有变化随机性强的特点.本文利用 GGOS Atmosphere提供的2002—2009年全球天顶对流层延迟格网时间序列研究了全球对流层天顶延迟的时空变化特征.并以此为基础对全球天顶对流层延迟(Zenith Troposphere Delay, ZTD)进行建模,提出了一种基于球谐函数的全球非气象参数对流层天顶延迟改正模型——GZTD模型.实验对比结果表明考虑ZTD经纬向变化的GZTD模型内符合精度全球统计结果(bias:0.2 cm,RMS:3.7 cm)优于只考虑ZTD纬向变化的UNB3m (bias:3.4 cm,RMS:6.0 cm)、UNB4 (bias:4.7 cm,RMS:7.4 cm)、UNB3 (bias:4.0 cm,RMS:7.0 cm)和EGNOS (bias:4.5 cm,RMS:6.9 cm)等模型.使用全球385个IGS站进行外符合检验,统计结果表明GZTD模型(bias:-0.02 cm,RMS:4.24 cm)同样优于其它模型.GZTD模型具有改正效果良好、使用简单、所需参数少等优点.     

A stochastic model of the Fourier phase of strong ground motion

ＡｓｔｏｃｈａｓｔｉｃｍｏｄｅｌｏｆｔｈｅＦｏｕｒｉｅｒｐｈａｓｅｏｆｓｔｒｏｎ　ｇｒｏｕｎｄｍｏｔｉｏｎＺｈｅｎ－ＰｅｎｇＬＩＡＯ；（廖振鹏）ａｎｄＸｉｎｇＪＩＮ（金星）（ＩｎｓｔｉｔｕｔｅｏｆＥｎｇｉｎｅｅｒｉｎｇＭｅｃｈａｎｉｃｓ，ＳｔａｔｅＳ...     

Seismic ground motion relationships in southern China based on stochastic finite-fault model

The characteristics of seismic ground motions in southern China are difficult to determine statistically due to a lack of strong ground motion data. In this study, a stochastic finite-fault ground motion model was adopted to simulate the seismic ground motions at bedrock for southern China, based on parameters derived from small and medium earthquakes that have occurred in the region. From these, the response spectra was estimated. A set of ground motion attenuation relationships was then developed based on simulated peak ground motions and response spectral parameters through regression, which would be applicable for use in engineering practice. Through comparisons, it was demonstrated that the proposed ground motion relationships are generally consistent with those obtained from other reported ground motion attenuation models for southern China.     

Seismic ground motion relationship in southern China based on stochastic finite-fault model

The characteristics of seismic ground motions in southern China are difficult to determine statistically due to a lack of strong ground motion data. In this study, a stochastic finite-fault ground motion model was adopted to simulate the seismic ground motions at bedrock for southern China, based on parameters derived from small and medium earthquakes that have occurred in the region. From these, the response spectra was estimated. A set of ground motion attenuation relations hipswas then developed based on simulated peak ground motions and response spectral parameters through regression, which would be applicable for use in engineering practice. Through comparisons, it was demonstrated that the proposed ground motion relationships are generally consistent with those obtained from other reported ground motion attenuation models for southern China.     

上地幔俯冲板块的动力学过程:数值模拟

大洋板块俯冲到地幔转换带,进而可形成不同的形态：板块可以停滞在660 km不连续面,抑或穿过地幔转换带进入下地幔.这些不同的俯冲模式可进一步影响到海沟的运动.为更好地理解上地幔中俯冲板片的变形行为以及俯冲过程与海沟运动之间的关系,本文通过建立一系列高精度二维热-力学自由俯冲的数值模型,揭示了俯冲板块在上地幔中的变形方式及其与地幔转换带之间的相互作用过程.模拟结果显示,在俯冲板块与地幔转换带的相互作用过程中,其动力学过程可以分为以海沟后撤主导、海沟前进主导以及稳定型海沟等三种主要动力学类型.对于年龄较老,厚度较大的俯冲板块容易形成海沟后撤型俯冲,俯冲板块停滞在660 km不连续面.相反,年龄较小,塑性强度较小的板块容易形成海沟前进型俯冲,俯冲板块穿越660 km不连续面.

     

A simple ground‐motion prediction model for cumulative absolute velocity and model validation

Cumulative absolute velocity (CAV) is an important ground motion intensity measure used in seismic hazard analysis. Based on the Next Generation Attenuation strong motion database, a simple ground‐motion prediction equation is proposed for the geometric mean of as‐recorded horizontal components of CAVs using mixed regression analysis. The proposed model employs only four parameters and has a simple functional form. Validation tests are conducted to compare the proposed model with the recently developed Campbell–Bozorgnia (CB10) model using subsets of the strong motion database, as well as several recent earthquakes that are not used in developing the model. It is found that the predictive capability of the proposed model is comparable with the CB10 model, which employs a complex functional form and more parameters. The study also corroborates previous findings that CAV has higher predictability than other intensity measures such as the peak ground acceleration. The high predictability of CAV warrants the use of the proposed simple model as an alternative in seismic hazard analysis. Copyright © 2012 John Wiley & Sons, Ltd.