三维速度场是GNSS连续运行站监测地壳位移的一种主要表达方式，使用GMT可以方便的绘制速度场。本文以陆态网共享的GNSS数据为例，分享速度场的绘制方法。
数据下载
中国地震局GNSS数据产品服务平台：http://ftp.cgps.ac.cn/
现在需要申请账号，才可以使用FTP。如果没有账号，也可以使用小编已经整理好的数据（地址：https://www.jianguoyun.com/p/Dco7BIEQ2PCQBxjQ4esD）。
用到result.d2 和result.d4数据文件。
GMT模块
主要用到psvelo。
注意：该模块在不同的GMT版本中有较大更新，之前的祖传代码可能在新的GMT版本下无法使用。本例子使用GMT6版本。
水平速度场
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
| #!/bin/bash
gmt gmtset MAP_FRAME_WIDTH = 0.1c FORMAT_GEO_MAP = ddd:mm:ssF
gmt gmtset FONT_ANNOT_PRIMARY = 7p
gmt gmtset FONT_LABEL 7 MAP_LABEL_OFFSET 5p
gmt gmtset MAP_FRAME_TYPE plain MAP_TICK_LENGTH_PRIMARY -5p MAP_FRAME_PEN black
R=70/140/10/55
J=M5i
PS=xmd5.ps
gmt psbasemap -R$R -J$J -P -K -Baf -BnwSE >$PS
gmt pscoast -R$R -J$J -O -Dl -N1/0.1p,gray -W1/0.2p,gray -K -A1000 -I1/0.1p,lightblue>> $PS
awk '{print $1,$2,$5,$3,$6,$4,0}' ../../result.d4 | gmt psvelo -J -R -Se0.015c/0.95/0 -A0.15c+e+p0.75p -Gblue -W0.1p,blue -K -O >> $PS
gmt pstext -J -R -F+f7p+jML -K -O >> $PS << EOF
130 15 30\2611 mm/yr
EOF
gmt psvelo -J -R -Se0.015c/0.95/0 -A0.15c+e+p0.75p -W0.2p,blue -O >> $PS << EOF
130 14 30 0 1 1 0
EOF
gmt psconvert $PS -P -Tg -A
|
垂直方向速度场
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
| #!/bin/bash
R=70/140/10/55
J=M5i
PS=xmd6.ps
gmt gmtset MAP_FRAME_WIDTH = 0.1c FORMAT_GEO_MAP = ddd:mm:ssF
gmt gmtset FONT_ANNOT_PRIMARY = 7p
gmt gmtset FONT_LABEL 7 MAP_LABEL_OFFSET 5p
gmt gmtset MAP_FRAME_TYPE plain MAP_TICK_LENGTH_PRIMARY -5p MAP_FRAME_PEN black
gmt psbasemap -R$R -J$J -P -K -Baf -BnwSE >$PS
gmt pscoast -R$R -J$J -O -Dl -N1/0.1p,gray -W1/0.2p,gray -K -A1000 -I1/0.1p,lightblue >> $PS
awk '{print $2,$3,0,$4,0,$5,0}' ../../result.d2 | gmt psvelo -J -R -Se0.08c/0.95/0 -A0.15c+e+p0.75p -W0.1p,blue -K -O >> $PS
awk '{print $2,$3,0,$4,0.15,$5,0}' ../../result.d2 | gmt psvelo -J -R -Se0.08c/0.95/0 -A0.15c+e+p0.75p -Gred -W0.1p,red -K -O >> $PS
gmt pstext -J -R -F+f7p+jML -K -O >> $PS << EOF
120 15 20\2611 mm/yr
EOF
gmt psvelo -J -R -Se0.08c/0.95/0 -A0.15c+e+p0.75p -W0.2p,blue -O >> $PS << EOF
120 14 20 0 1 1 0
EOF
gmt psconvert $PS -P -Tg -A
|
中国大陆速度场特征
水平方向上喜马拉雅山正在向长春方向移动，一年的位移量超过3cm；沿海陆地正在向太平洋移动，速度也约为3cm/每年。在青藏高原和四川盆地、甘肃等地区有明显的位移切变，从而使得我国大陆的速度场呈现出螺旋特征。（此处蕴含1篇顶级SCI。）
垂直方向上，大部分位移不大，少数地方如天津出现特别大的沉降，也有少数地区出现较大的上升（武汉）。沿海区域的沉降特征是大概以山东为边界，北边上升，南边下沉，因此在分析海平面变化的时候需要注意加入验潮站沉降速度差异改正。（此处蕴含1篇顶级SCI。）
是不是忽然恍然大悟，关注本公众号，还可以搞几篇SCI呀。哈哈.