关于板块

地球是由多个刚体板块构成的，在Bernese中认可的板块模式是NNR–NUVEL–1A。
关于NUVEL-1A的相关文献如下：

DeMets. C., R. G. Gordon, D. F. Argus, and S. Stein, Effect of recent revisions to the geomagnetic reversal time scale on estimate of current plate motions, Geophys. Res. Lett., vol. 21, no. 20, 2191-2194, 1994.
(Revised version of NUVEL-1, Euler poles are same)
[Note 1] Model parameters for Philippine Sea Plate are based on Seno et al. (JGR, 1993)
[Note 2] Velocities calculated based on NUVEL-1 and NUVEL-1A may be 4.5% and 2% faster than those measured by space geodetic methods by using VLBI/SLR (Gordon, Nature 1993)

这个NUVEL-1A板块划分图如下：
谷歌的kml文件下载地址：https://download.csdn.net/download/zzh_my/12171984

NNR–NUVEL–1A板块模型

当然还有其他板块划分，例如MORVEL板块模型

MORVEL板块模型

---

计算地球上任何一点的相对和绝对板块运动方向和速度，公式：
参考：http://ofgs.aori.u-tokyo.ac.jp/~okino/platecalc_new.html

rotation vector
MωF= MωP - FωP = (ωx, ωy, ωz)
rotation pole
plat = atan(ωz / sqrt (ωx*ωx+ωy*ωy))
plon = atan(ωy/ωx)
rotation rate
rate = sqrt(ωx*ωx+ωy*ωy*+ωz*ωz)
velocity at a point (tlat, tlon)
Vew = radius * cos (tlat)*rate*(sin(plat)-cos(plat)*tan(tlat)*cos(plon-tlon))
Vns = radius * rate * cos(plat)*sin(tlon-plon)
Velocity = sqrt (Vew*Vew + Vns*Vns)
Direction = atan(Vew/Vns)

例如在纬度45，经度126，板块运动速率：

plate velocity : 90 [mm/yr]=[km/my]
direction: 105.7 [deg. from North]
north component of velocity : -24 [mm/yr]=[km/my]
east component of velocity : 87 [mm/yr]=[km/my]
Parameters used for calculation
Plate Model: NUVEL-1A
Moving Plate: eu
Fixed Plate: pa
rotation vector omega: 0.000529,-0.007235,0.013123
Latitude of Euler pole: 61.0664653390412 [deg.]
Longitude of Euler pole: -85.818175722285 [deg.]
Angular velocity: 0.859127791068401 [deg./m.y.]
Latitude inputted: 45 [deg.]
Longitude inputted: 126 [deg.]

学习参考