用VLBI和SLR实测数据解算现时板块运动参数的方法

中国科学院上海天文台年刊
ＣＡＮｏ．１６，１９９５
１９９５年第１６期ＡＮＮＡＬＳＯＦＳＨＡＮＧＨＡＩＯＢＳＥＲＶＡＴＯＲＹＡＣＡＤＥＭＩＡＳＩＮＩ
用ＶＬＢＩ和ＳＬＲ实测数据解算现时
板块运动参数的方法
孙付平赵铭
主题词：甚长基线干涉测量（ＶＬＢＩ）—人卫激光测距（ＳＬＲ）—板块运动
提要
本文着重讨论了用空间大地测量（ＶＬＢＩ和ＳＬＲ）实测数据解算现时板块运动参数的方法，推导了有关公式。

根据作者的实践，阐述了应用中应注意的问题。

一、引言
板块构造学认为，地球表面岩石圈是由大小不等的２０多个板块构成。

由于地幔对流及其它因素的作用，这些板块之间存在着１ １０ｃｍ燉ｙｒ的相对运动。

利用海底扩张速率、转换断层方位角、地震滑动矢量等地质和地球物理观测资料，人们已建立了能定量反映板块间相对运动的地学板块运动模型，例如ＲＭ２［１］和ＮＵＶＥＬ １［２］。

根据板块运动模型可以估计出地球表面任意两点之间的相对运动。

但是，海底扩张速率是由海岭两侧年龄相同的磁异常条带的间距除以该磁异常条带的年龄而得到的，它反映的是最近几百万年里板块扩张的平均效果，而转换断层的演化则有更长的地质历史。

所以，地学板块运动模型估计的是最近几百万年内板块的平均运动。

甚长基线干涉测量（ＶＬＢＩ）、人卫激光测距（ＳＬＲ）等空间大地测量技术，目前已能以几个ｍｍ燉ｙｒ的精度测定板块间的运动。

这种实测的板块运动是由最近１０多年的观测资料解算的，与漫长的百万年时间尺度相比，完全可以认为是现时的板块运动。

现时的板块运动与地质历史上平均的板块运动相比较，将为我们研究板块运动的稳定性、板块的刚性等重大问题提供证据。

并且，要为地震等灾害预测工作服务，要维持一个高精度地球参考架，建立一个精确的现时板块运动模型也是很必要的。

地球物理研究表明，地面点的长期运动主要有三种：（１）大尺度板块运动：主要沿水平方向，量级是１－１０ｃｍ燉ｙｒ，是全球性的。

（２）局部构造形变：主要沿水平方向，量级是１－２０ｍｍ燉ｙｒ，是局部性的。

（３）冰期后地壳回弹：主要沿垂直方向，量级是１－１０ｍｍ燉ｙｒ，是区域性的。

实测的台站运动主要是以上三种运动的叠加。

因此，用实测数据求解现时板块运动参数时，必须尽可能地消除局部形变和冰期后地壳回弹的影响，这是应用中应注意的问题。

１９９４年４月１５日收到。

二、实测的板块运动观测量
１．ＶＬＢＩ
观测量［３］
测地ＶＬＢＩ是通过观测射电信号至两个ＶＬＢＩ接收机的时间延迟和延迟率来解算站间基线矢量的。

不同时段得到的同一ＶＬＢＩ
基线的各观测值投影到同一个初始参考架中，对这些观测量进行时间拟合即可得到该基线沿其长度、横向和垂向的三个时间变化率：爜 牏牐
、爴
牏牐
和爣 牏牐。

由定义可知［３］
，基线横向变化率爴 牏牐与初始参考架定向有关；而基线长度变化率爜 牏牐
只与尺度有关；基线垂向变化率爣 牏牐
是ＶＬＢＩ确定最差的一个分量，它主要受大气因素的影响。

所以用基线变化率求解板块运动参数时一般只采用基线的长度和横向变化率。

图１显示了ＶＬＢＩ基线变化的三个分量方向。

图１ＶＬＢＩ基线变化的三个分量Ｆｉｇ．１ＴｈｒｅｅｃｏｍｐｏｎｅｎｔｓｏｆＶＬＢＩ
ｂａｓｅｌｉｎｅｖａｒｉａｔｉｏｎ
其中，
牄
牏牐＝（牜 牐－牜 牏
）燉牜 牐－牜
牏，（１）牠
牏牐＝（牜
牐×牜
牏
）燉牜
牐×牜
牏，（２）
牎
牏牐＝牄
牏牐×牠
牏牐
（３）分别是基线沿长度、横向和垂向的单位矢量。

牠
牏牐的方向是从正面看由里向外。

为了更直观地分析各台站的运动，ＶＬＢＩ
数据处理中心还解算了各台站在地心参考架中的站
速度爼牨牏、爼牪牏和爼牫
牏，以及在站心参考架中的站速度爼牉牏、爼牕牏和爼牣
牏。

因站速度是绝对量，
解算时必须附加相应的约束条件。

例如把某一ＶＬＢＩ站（
如Ｗｅｓｔｆｏｒｄ
）的站速度及到另一站的方位变化固定为板块绝对运动模型（如ＮＮＲ－ＮＵＶＥＬ
１［４］
）的估计值，而其它站的站速度相对于它来确定。

实际上，ＶＬＢＩ
基线的长度、横向和垂向变化率与站心坐标速度之间还有如下关系：［５］
爜 牏牐爴 牏牐爣 熿燀燄燅牏牐＝ 牄 牏牐燈牕 牏 牄 牏牐燈牉 牏 牄 牏牐燈牣 牏牄 牏牐燈牕 牐牄 牏牐燈牉 牐牄 牏牐燈牣 牐 牠 牏牐燈牕 牏 牠 牏牐燈牉 牏 牠 牏牐燈牣 牏牠 牏牐燈牕 牐牠 牏牐燈牉 牐牠 牏牐燈牣 牐 牎 牏牐燈牕 牏 牎 牏牐燈牉 牏 牎 牏牐燈牣 牏牎 牏牐燈牕 牐牎 牏牐燈牉 牐牎 牏牐燈牣 熿燀燄燅牐爼牕
牏
爼牉牏爼牣牏爼牕
牐
爼牉牐
爼熿燀燄燅牣牐
（４）其中牕 、牉 和牣
分别为站心参考架中的当地北向、东向和天顶方向的单位矢量。

由（４）
式也可求解ＶＬＢＩ
站速度，因基线变化率是相对量，求解时同样需要上述约束条件。

测站地心坐标速度爼牨牏、爼牪牏和爼牫牏与站心坐标速度爼牕牏、爼牉牏和爼牣
牏之间也有如下关系：
爼牨牏爼牪牏爼熿燀燄燅牫牏＝ ｓｉｎ犧牏 ｃｏｓ犧牏ｓｉｎ 牏ｃｏｓ犧牏ｃｏｓ 牏ｃｏｓ犧牏 ｓｉｎ犧牏ｓｉｎ 牏ｓｉｎ犧牏ｃｏｓ 牏０ｃｏｓ 牏ｓｉｎ 熿燀燄
燅牏爼牉
牏爼牕
牏
爼熿燀燄燅牣
牏
（５）
其中 牏和犧牏是测站牏
的地心纬度和地心经度。

８
中国科学院上海天文台年刊
１９９５年
可见，ＶＬＢＩ提供三种板块运动观测量：（１）基线沿其长度和横向的变化率爜 牏牐和爴 牏牐。

它们都可能受到局部形变的影响，并且当基线较长时爜 牏牐
还可能受到冰期后地壳回弹的影响；（２）地心坐标速度爼牨牏、爼牪牏、爼牫牏。

可能受到局部形变和冰期后地壳回弹的影响；（３）站心坐标速度爼牕
牏
、爼牉牏、爼牣
牏。

其中水平分量可能受局部形变影响，
垂直分量则主要反映冰期后地壳回弹引起的变化。

２．ＳＬＲ观测量［６］
ＳＬＲ技术通过测量地面站至人造卫星的距离来解算地面站地心位置。

不同时段的观测会
得到不同的地心坐标。

由同一时段所得到的两点的站坐标可解算站间大地线距离，不同时段里所得距离的时间拟合即可得到站间大地线变化率爳 牏牐。

爳
牏牐
与参考架定向无关，只与尺度有关。

同样，提供某种约束条件（例如ＮＮＲ－ＮＵＶＥＬ
１），也可由不同时段的测站地心坐标求解测站的地心坐标速度爼牨牏、爼牪牏、爼牫牏和站心坐标速度爼牕牏、爼牉牏、爼牣
牏。

可见，ＳＬＲ也提供三种板块运动观测量：（１）站间大地线变化率爳 牏牐，因大地线与当地垂线
基本垂直，所以爳 牏牐对冰期后地壳回弹的垂直分量不敏感；（２）地心坐标速度；（３）站心坐标速度。

３．基线长度变化率爜 牏牐与大地线（弧长）变化率爳 牏牐
的关系因地球常作为球体来处理，大地线变化率近似为弧长变化率。

如图１所示，爜牏牐、爳牏牐为第牏站到第牐站的基线长和弧长，犜牏牐为两站间的地心角，
设爲是地球平均半经，则有爳牏牐＝爲燈犜牏牐（６）爜牏牐＝２
爲燈ｓｉｎ（犜牏牐２
）（７）
把上面两式对时间牠求导并进行代换，可得
爜 牏牐＝爳 牏牐燈ｃ
ｏｓ犜牏牐２
（８）因爜 牏牐和爳 牏牐
对板块间沿基线横切方向的运动不敏感，这种不足需要较多分布合理的基线长度或弧长变化率来弥补，所以实用中常需组合应用ＶＬＢＩ和ＳＬＲ两种数据。

如果两者尺度
一致，则可用（８）式进行两种观测量之间的变换。

值得注意的是，当站间地心角犜牏牐接近１
８０°时，数据中的误差会被显著地放大或缩小。

所以（８）式的应用应限定在一定的范围，如犜牏牐≤１
２０°。

三、用实测数据求解板块运动参数的数学模型
地球近似为椭球体，但在研究板块运动时，为方便起见常把地球近似为球体。

这种近似引起的误差是可以忽略的。

例如对于３毫秒的地心角（对应地面距离约１０ｃｍ）
来说，地球半经１０ｋｍ的变化引起的地面距离变化仅为０．１５ｍｍ。

对于几个ｃｍ燉ｙｒ的板块运动来说，地面点地心距变化的影响是很小的，应用中常取地球平均半径爲（爲＝６３７８１４０ｍ）
来计算。

板块的运动常采用旋转角速度矢量（也叫欧拉矢量）犓
来描述。

犓
的球坐标分量为牼、犎、犽，牼和犎是板块旋转极的经纬度，犽是旋转角速度。

犓
的直角坐标分量是犓牨、犓牪、犓牫，它们分别是犓
在牀、牁、牂轴上的投影。

板块运动有绝对和相对运动之分，
绝对运动是板块相对于某一参考架（
例如热点参考架）的运动，相对运动则是一个板块相对于另一板块的运动。

通常是求解板块之间的相对运动参数。

９
第１６期用ＶＬＢＩ和ＳＬＲ实测数据解算现时板块运动参数的方法
１．由基线变化率爜 牏牐和爴 牏牐
求解板块运动参数如图１所示，设牏站位于牑板块，牐站位于牓板块，两板块的欧拉矢量分别为犓
牑和犓
牓，
则有爜 牏牐
＝（犓 牓×牜 牐－犓 牑×牜 牏）燈牄 牏牐（９）
爴
牏牐＝（犓 牓×牜 牐－犓 牑×牜 牏）燈牠 牏牐
（１０）
设犓 牑牓＝犓 牓－犓
牑为牓
板块相对牑板块的欧拉矢量，则代入（９）和（１０）式变换后可得爜
牏牐＝（牜
牏×牜
牐）牜 牐－牜
牏
燈犓 牑牓（１１）
爴 牏牐＝犓 牑牓×牜 牐＋犓 牑×（牜 牐－牜 牏
）燈牠 牏牐（１２）可见爜 牏牐只与板块间相对运动有关，而爴 牏牐
与参考架有关。

由（１１）和（１２）式可以求解板块间的相对运动欧拉矢量和至少一个板块的绝对欧拉矢量。

２．由站间大地线变化率爳 牏牐
求解板块运动参数
把（１１）式代入（８）式即可得到爳 牏牐
与犓
牑牓的关系式：爳 牏牐＝爜 牏牐ｃｏｓ犜牏牐２＝（牜 牏×牜 牐）ｃｏｓ犜牏牐
２
燈牜 牐－牜
牏燈犓 牑牓
（１３）如果用球坐标表示，则有比较复杂的形式。

如图２所示，爮牏（犺牏、犧牏
）为牑板块上的测站，爮牐（犺牐、犧牐）为牓板块上的测站，两板块的相对欧拉矢量犓
牑牓
用球坐标表示为犎、牼、犽。

ＣＴＰ为协议地极，取爲为地球平均半径，则由图２可找出犓
牑牓与在△牠内弧长变化△爳牏牐的关系式为
［７］
△爳牏牐＝爲燈犽燈△牠燈ｓｉｎ犲燈ｓｉｎ犜（１４）
根据球面几何关系，（１４）
式又可变换为：爳 牏牐
＝犽燈爲ｓｉｎ爳牏牐
ｓｉｎ犎ｃｏｓ犺牏ｃｏｓ犺牐ｓｉｎ（犧牏 犧牐）＋ｃｏｓ犎ｃｏｓ犺牏ｓｉｎ犺牐ｓｉｎ（犧牐 牼） ｃｏｓ犎ｓｉｎ犺牏ｃｏｓ犺牐ｓｉｎ（犧牐［］ 牼）（１５）
图２Ω
ｋｌ与弧长变化△Ｓｉｊ
的关系Ｆｉｇ．２Ｒｅｌａｔｉｏｎｓｈｉｐｏｆａｒｃｌｅｎｇｔｈｖａｒｉａｔｉｏｎ△ＳｉｊａｎｄΩ
ｋｌ３．由站心坐标速度爼牕牏和爼牉
牏求解板块运动参数
设爼牕牏、爼牉牏、爼牣牏为牑板块上第牏站的站心坐标速度，因板块运动主要沿水平方向，理论上爼牣
牏
与板块运动无关则站速度水平分量牕牏和牉与板块绝对欧拉矢量
牑
牨牑牪牑牫
牑的关系式为０
１中国科学院上海天文台年刊
１９９５年
爼牕
牏
爼［］牉
牏＝爲ｓｉｎ犧牏爲ｃｏｓ犧牏
０ 爲ｓｉｎ犺牏ｃｏｓ犧牏 爲ｓｉｎ犺牏ｓｉｎ犧牏爲ｃｏｓ犺［
］
牏犓牨
牑犓牪牑
犓熿燀燄
燅牫
牑
（１６）
式中，爲为地球平均半径，犧牏、犺牏为第牏
站地心经纬度。

如图３所示，如果犓
牑用球坐标犎牑、牼牑、犽牑表示，则ｋ板块上爮牏站在时间△牠
内因板块运动产生的坐标变化△犺牏、△犧牏与犎牑、牼牑、犽牑的关系为［７］
ｓｉｎ（犺牏＋△犺牏）＝ｓｉｎ犎牑ｃｏｓ犲＋ｃｏｓ犎牑ｓｉｎ犲ｃｏｓ（犜－犽牑燈
△牠）ｃｏｓ（犧牏－牼牑＋△犧牏
）＝ｃｏｓ犎牑ｃｏｓ犲－ｓｉｎ犎牑ｃｏｓ（犜－犽牑燈△牠）ｓｉｎ（犜－犽牑烅烄烆燈
△牠）（１７）
式中犜是方位角，犲是旋转极到爮牏点的弧距。

对（１７）式线性化即可得到坐标变化△犺牏和△犧牏与板块运动参数犎牑、牼牑、犽牑的关系式：
△犺牏＝ｃｏｓ犎牑ｓｉｎ（犧牏－牼牑）燈犽牑
△牠△犧牏＝（ｓｉｎ犎牑－ｃｏｓ（犧牏－牼牑）ｔｇ犺牏ｃｏｓ犎牑）燈犽牑
烅
烄烆△牠（１８）
或者把（１８）
式变换为爼牕牏＝
爲燈△犺
牏△牠
＝爲燈犽牑ｃｏｓ犎牑ｓｉｎ（犧牏－牼牑）爼牉
牏＝爲燈△犧牏ｃｏｓ犺牏△牠＝爲燈犽牑（ｓｉｎ犎牑ｃｏｓ犺牏－ｃｏｓ（犧牏－牼牑）ｓｉｎ犺牏ｃｏｓ犎牑
烅烄烆）（１９）
若取犓牨牑＝犽牑ｃｏｓ犎牑ｃｏｓ犧牑、犓牪牑＝犽牑ｃｏｓ犎牑ｓｉｎ犧牑、犓牫
牑＝犽牑ｓｉｎ犎牑代入（
１６）式，则可证明（１６）式与图３△犺牏、△犧牏与犎牑、牼牑、犽牑的关系Ｆｉｇ．３Ｒｅｌａｔｉｏｎｓｈｉｐｏｆ犎牑、牼牑、犽牑
ａｎｄ△犺牏、△犧牏
（１９）式是等价的。

由（１６）或（１９）式解出的是板块的绝对欧拉矢量，它与定义站速度的参考架有关。

我们需要的是板块间的相对运动参数，因此还须把解得的绝对欧拉矢量换算为相对欧拉矢量（犓
牑牓＝犓
牓－犓
牑
）。

也可用所有板块上的站速度由（１６）式联合解算所有板块对的相对欧拉矢量和一个绝对欧拉矢量。

４．由地心坐标速度爼牨牏、爼牪牏、爼牫
牏求解板块运动参数设欧拉矢量为犓
牑的牑
板块上有牕个测站，则第牏站的站速度爼
牏（爼牨牏、爼牪牏、爼牫牏）与犓
牑（犓牨牑、犓牪牑、犓牫
牑）的关系为：爼牨牏爼牪牏爼熿燀燄燅牫牏＝０爲ｓｉｎ犺牏 爲ｃｏｓ犺牏ｓｉｎ犧牏 爲ｓｉｎ犺牏０爲ｃｏｓ犺牏ｃｏｓ犧牏爲ｃｏｓ犺牏ｓｉｎ犧牏 爲ｃｏｓ犺牏ｃｏｓ犧牏熿燀燄
燅０犓牨
牑犓牪
牑
犓熿燀燄燅牫牑
（２０）
根据前面导出的ＶＬＢＩ和ＳＬＲ实测数据与板块运动参数的关系式，
当观测量大于待解参数个数时，就可用最小二乘平差方法解算出板块运动参数。

１
１第１６期用ＶＬＢＩ和ＳＬＲ实测数据解算现时板块运动参数的方法
四、应用中需注意的问题
１．稳定台站的选择
在引言中我们提到，
实测的台站运动除观测误差外，主要是大尺度板块运动、局部构造形变和冰期后地壳回弹这三种长期运动的叠加。

我们感兴趣的是板块运动信息，但它却受到局部形变和冰期后地壳回弹的干扰。

特别是局部形变，它的影响有时很大。

对于局部形变的影响，可通过选择稳定的测站来消除。

板块的边界和内部一些地区虽然存在着形变，但板块的主体是稳定的。

位于板块稳定主体内的台站，其实测运动主要反映了板块的整体运动。

稳定的测站可由如下两种方法判断：（１）台站是否位于板块边界区或板块内明显的构造形变区；（２）同一板块内各台站间的基线长度（弧长）变化率是否接近于零，如采用站速度则要看各台站的站速度是否一致。

综合应用这两种方法，完全可以判定台站是否位于稳定区，选择位于稳定区的台站的观测量来求解板块运动参数，一般可以基本上消除局部形变的影响。

对于冰期后地壳回弹，因其主要沿垂直方向，对站间距离较长的基线长度变化率和地心坐标速度会有一定影响，而对弧长变化率和站心坐标速度的水平分量影响较小。

消除这种影响
也有两种方法：（１）用冰期后地壳回弹模型（如ＩＣＥ－３Ｇ）估计各站的运动，然后从实测数据中消除其影响。

但因目前的冰期后地壳回弹模型精度和分辨率还不太高，冰期后地壳回弹量级又比较小，所以这种方法不一定可行。

（２）采用对地壳垂直运动不敏感的观测量，如弧长变化率和站心坐标速度的水平分量，即采用（１３）或（１５）式、（１６）或（１９）式来求解板块运动参数，这在目前来说是可行的。

用基线变化率或弧长变化率求解板块运动参数时，要求基线两端测站都位于板块内稳定地区。

如果一端测站位于形变区，
观测的基线（或弧长）变化率就会受到形变的干扰，这个观测量就不能采用。

在目前观测的基线尚少，而位于形变区的台站又较多的情况下，完全满足要求的基线（或弧长）变化率观测量还比较少，所以用（１１）和（１２）式及（１３）或（１５）式求解板块运动参数的可靠性和精度都不太高。

与此相比，若用站心坐标速度（１６）或（１９）
式求解板块运动参数则有如下优点［５］：（１）因站速度是根据所有观测量导出的，它包含有更多的观测信息；（２）
站速度能更直观地反映出板块的整体运动及测站之间的相对形变，便于判定稳定台站；（３）站速度个数远小于基线（
或弧长）变化率个数，便于进行残差分析。

总之，无论是着眼于消除局部形变和冰期后地壳回弹的影响，或者是着眼于利用更多的观测信息，用站心坐标速度水平分量求解板块运动参数都具有优越性。

根据我们的实践，用（１６）
和（１９）式求解的结果基本一致，但因（１９）式是非线性形式，列误差方程时要进行近似处理，也比较繁琐，建议采用（１６）
式。

２．台站分布对解算结果的影响
就最小二乘估计方法而言，要精确估计板块运动参数，不仅要求观测量精度要高，还要求观测台站要尽可能多，并且尽可能均匀分布于整个板块稳定地区。

在实际应用中，观测台站的几何分布对估计结果的精度影响很大。

例如在太平洋和北美板块上，目前的台站分布比较均匀合理，所以求解的太平洋－北美板块相对运动欧拉矢量精度就比较高。

而欧亚板块上稳定的台站虽不算太少，但它们都位于西欧一小块地区，并且与北美几个台站的连接基线大体沿纬线方
２
１中国科学院上海天文台年刊
１９９５年
向，所以观测量对估计参数的约束不够强，致使估计结果误差较大。

３．ＶＬＢＩ和ＳＬＲ实测数据的联合应用
对于基线长度（弧长）变化率而言，ＶＬＢＩ与ＳＬＲ观测量之间主要表现为尺度差；
而对于站速度而言，ＶＬＢＩ与ＳＬＲ观测量之间却存在着尺度差和参考架定向之差。

因为目前各数据处理中心对ＶＬＢＩ和ＳＬＲ数据的处理都采用一致的光速牅和地球引力常数爢爩，并采用一致的相对论改正模型，所以两种观测量之间的尺度差完全可以忽略，主要区别是参考架定向之差（对站速度而言）。

因板块相对运动观测量爜 牏牐（ＶＬＢＩ）和爳 牏牐（ＳＬＲ）对板块间沿基线横切方向（牠 牏牐）
的相对运动不敏感，用它们求解板块运动参数时，要求基线在两板块上的几何分布尽可能覆盖所有板块
稳定区域，这样才能弥补爜 牏牐和爳 牏牐
包含板块运动信息量的不足。

在目前ＶＬＢＩ和ＳＬＲ观测站都不太多的情况下，
组合应用这两种观测量可加强对板块运动参数的约束。

当测站间地心角较小（犜牏牐≤１
２０°）时，ＶＬＢＩ和ＳＬＲ观测量可用（８）式相互换算，并用（１１）或（１３）式求解板块运动参数。

对于ＶＬＢＩ和ＳＬＲ站速度观测量，
必须统一在同一个参考架中才能组合应用。

统一的方法是：（１）由一些ＶＬＢＩ与ＳＬＲ并置站上的观测量解算两组站速度之间的平移和旋转参数，然后变换至同一参考架；（２）当并置站较少或不够时，可选用某个板块不动作为固定参考架，在这个板块的稳定地区应分布有较多的ＶＬＢＩ和ＳＬＲ站。

分别从所有ＶＬＢＩ（或ＳＬＲ
）台站的站速度中剔除位于板块内稳定地区的ＶＬＢＩ（或ＳＬＲ）站的整体运动部分，然后才能将两组站速度组合在一起使用。

这样做对估计板块相对运动参数没有影响。

参
考
文
献
［１］Ｍｉｎｓｔｅｒ，Ｊ．Ｂ．ａｎｄＪｏｒｄａｎ，Ｔ．Ｈ．，Ｊ．Ｇｅｏｐｈｙｓ．Ｒｅｓ．，８３（１９７８），５３３１－５３４８．［２］ＤｅＭｅｔｓ，Ｃ．ｅｔｃｒａｌ．，Ｇｅｏｐｈｙｓ．Ｊ．Ｉｎｔ
．，１０１（１９９０），４２５－４７５．［３］Ｃａｐｒｅｔｔｅ，Ｄ．Ｓ．ｅｔｃｒａｌ．，ＧｕｓｔａｌＤｙｎａｍｉｃｓＰｒｏｊｅｃｔＤａｔａＡｎａｌｙｓｉｓ１９９０（１９９１）．［４］Ａｒｇｕｓ，Ｄ．Ｆ．ａｎｄＧｏｒｄ，Ｒ．Ｇ．，Ｇｅｏｐｈｙｓ．Ｒｅｓ．Ｌｅｔｔ．，１８（１９９１），２０３９－２０４４．［５］Ｗａｒｄ，Ｓ．Ｎ．，Ｊ．Ｇｅｏｐｈｙｓ．Ｒｅｓ．，９５（１９９０），２１９６５－２１９８１．［６］Ｓｍｉｔｈ，Ｄ．Ｅ．ｅｔｃｒａｌ．，Ｊ．Ｇｅｏｐｈｙｓ．Ｒｅｓ．，９５（１９９０），２２０１３－２２０４１．［７］Ｄｒｅｗｅｓ，Ｈ．，Ｂｕｌｌ．Ｇｏｏｄ
．，５６（１９８２），７０－７９． ﹪﹦ ﹦ ﹪ ﹥ ﹨ ﹥﹦ ﹦ ﹫ ﹩ ﹢ ﹦ ﹫
﹢ ﹢ ﹦ ﹦ ﹨ ﹣﹫﹢ ﹥ ﹥﹢ ﹢
ＳｕｎＦｕｐｉｎｇ
ＺｈａｏＭｉｎｇ
（ＳｈａｎｇｈａｉＯｂｓｅｒｖａｔｏｒｙ，ＣｈｉｎｅｓｅＡｃａｄｅｍｙｏｆＳｃｉｅｎｃｅｓ
） ┎┌┄┇ ┈ＶｅｒｙＬａｎｇＢａｓｅｌｉｎｅＩｎｔｅｒｆｅｒｏｍｅｔｒｙ－ＳａｔｅｌｌｉｔｅＬａｓｅｒＲａｎｇｉｎｇ－Ｐｌａｔｅｍｏｔｉｏｎ
﹢ ┈┉┇ ┉
Ｉｎｔｈｉｓｐａｐｅｒ，ｔｈｅｍｅｔｈｏｄｓｆｏｒｔｈｅｄｅｔｅｒｍｉｎａｔｉｏｎｏｆｔｈｅｐｌａｔｅｍｏｔｉｏｎｐａｒａｍｅｔｅｒｓｆｒｏｍ
ＶＬＢＩａｎｄＳＬＲｄａｔａｈａｖｅｂｅｅｎｄｉｓｃｕｓｓｅｄ，ａｎｄｓｏｍｅｒｅｌａｖａｎｔｆｏｒｍｕｌａｓｈａｖｅｂｅｅｎｄｅｒｉｖｅｄｔｏｏ
．Ａｃｃｏｒｄｉｎｇｔｏｏｕｒｃｏｍｐｕｔａｔｉｏｎ，ｓｏｍｅｐｒｏｂｌｅｍｔｏｗｈｉｃｈｓｈｏｕｌｄｂｅｐａｉｄａｔｔｅｎｔｉｏｎｉｎｔｈｅｉｒａｐ ｐｌｉｃａｔｉｏｎｈａｖｅｂｅｅｎｄｅｓｃｒｉｂｅｄ
．３
１第１６期用ＶＬＢＩ和ＳＬＲ实测数据解算现时板块运动参数的方法。