中国大陆及周边活动地块分布图

CN-plate-neighbor.dat | CN-block-L1.dat | CN-block-L1-deduced.dat | CN-block-L2.dat | CN-border-L1.dat

本文绘制了中国大陆及周边活动地块分布图。

  • CN-plate-neighbor.dat 为中国大陆周边板块边界数据
  • CN-block-L1.dat 为一级地块边界数据
  • CN-block-L1-deduced.dat 为一级地块推断边界数据
  • CN-block-L2.dat 为二级地块边界数据
  • CN-border-L1.dat 为国界线数据

数据由网友根据下图矢量化而成,并经本文作者做了一定的格式转换:

绘图效果如下:

中国大陆及周边活动地块分布图

#!/bin/bash

gmt begin BLOCK png
gmt gmtset MAP_ANNOT_OBLIQUE 6
gmt gmtset MAP_TICK_LENGTH_PRIMARY 0
gmt gmtset MAP_FRAME_TYPE plain

# ============底图
gmt basemap -R78/12/149/53+r -JB105/10/25/47/15 -Bx10 -By5 -BNEWS
gmt coast -Ggray95 -S83/216/238 -A5000 -Dh
gmt plot CN-border-L1.dat -W0.2p -Gwhite

# ============绘制板块边界
gmt plot CN-plate-neighbor.dat -W2.0p,2/138/210 -Sf0.5+t+l -G2/138/210
# ============绘制推断地块边界
gmt plot CN-block-L1-deduced.dat -W1.0p,2/138/210,-
# ============绘制二级地块边界
gmt plot CN-block-L2.dat -W1.0p,orange
# ============绘制一级地块边界
gmt plot CN-block-L1.dat -W1.0p,2/138/210

# ============活动地块标注
gmt text -F+f10p,37 << EOF
90 45 西域活动地块
122 46 东北亚活动地块
120 35 华北活动地块
110 27 华南活动地块
90 33 青藏活动地块
97 23 滇缅活动地块
EOF

# ============边界图例
gmt legend -DjTL+w1.5i+jTL+o0.2c/0.2c -F+gwhite+p0.5p --FONT_ANNOT_PRIMARY=8p,37 << EOF
S 0.3i f0.5+t+l 0.4i/0.3c 2/138/210 2.0p,2/138/210 0.7i 板块边界
S 0.3i - 0.50i 2/138/210 1.0p,2/138/210 0.7i  一级地块边界
S 0.3i - 0.44i - 1.0p,2/138/210,- 0.7i  推断地块边界
S 0.3i - 0.50i orange 1.0p,orange 0.7i 二级地块边界
EOF

# ============南海诸岛
gmt basemap -DjBR+w2c+t
gmt coast -R106/121/3/24 -JM2c -Ggray95 -S83/216/238 -A5000 -Dh 
gmt plot CN-border-L1.dat -R106/121/3/24 -JM2c -B0 -W0.2p -Gwhite
gmt end
gmt begin BLOCK png
gmt gmtset MAP_ANNOT_OBLIQUE 6
gmt gmtset MAP_TICK_LENGTH_PRIMARY 0
gmt gmtset MAP_FRAME_TYPE plain

REM ============底图
gmt basemap -R78/12/149/53r -JB105/10/25/47/15c -Bx10 -By5 -BNEWS
gmt coast -Ggray95 -S83/216/238 -A5000 -Dh
gmt plot CN-border-L1.dat -W0.2p -Gwhite

REM ============绘制板块边界
gmt plot CN-plate-neighbor.dat -W2.0p,2/138/210 -Sf0.5+t+l -G2/138/210
REM ============绘制推断地块边界
gmt plot CN-block-L1-deduced.dat -W1.0p,2/138/210,-
REM ============绘制二级地块边界
gmt plot CN-block-L2.dat -W1.0p,orange
REM ============绘制一级地块边界
gmt plot CN-block-L1.dat -W1.0p,2/138/210

REM ============活动地块标注
echo 90 45 西域活动地块 > blockname
echo 122 46 东北亚活动地块 >> blockname
echo 120 35 华北活动地块 >> blockname
echo 110 27 华南活动地块 >> blockname
echo 90 33 青藏活动地块 >> blockname
echo 97 23 滇缅活动地块 >> blockname
gmt text blockname -F+f10p,37

REM ============边界图例
echo S 0.3i f0.5+t+l 0.4i/0.3c 2/138/210 2.0p,2/138/210 0.7i 板块边界 > legend.dat
echo S 0.3i - 0.50i 2/138/210 1.0p,2/138/210 0.7i  一级地块边界 >> legend.dat
echo S 0.3i - 0.44i - 1.0p,2/138/210,- 0.7i  推断地块边界 >> legend.dat
echo S 0.3i - 0.50i orange 1.0p,orange 0.7i 二级地块边界 >> legend.dat
gmt legend legend.dat -DjTL+w1.5i+jTL+o0.1c/0.1c -F+gwhite+p0.5p --FONT_ANNOT_PRIMARY=8p,37

REM ============南海诸岛
gmt basemap -DjBR+w2c+t
gmt coast -R106/121/3/24 -JM2c -Ggray95 -S83/216/238  -A5000 -Dh
gmt plot CN-border-L1.dat -R106/121/3/24 -JM2c -B0 -W0.2p -Gwhite
del blockname legend
gmt end
pause
#!/bin/bash
gmt gmtset MAP_ANNOT_OBLIQUE 6
gmt gmtset MAP_TICK_LENGTH_PRIMARY 0
gmt gmtset MAP_FRAME_TYPE plain

R=78/12/149/53r
J=B105/10/25/47/15c
PS=BLOCK.ps

gmt psxy -R$R -J$J -T -K > $PS

# ============底图
gmt psbasemap -R$R -J$J -Bx10 -By5 -BNEWS -K -O >> $PS
gmt pscoast  -R$R -J$J -Ggray95 -S83/216/238 -A5000 -Dh -K -O>> $PS
gmt psxy CN-border-L1.dat -R$R -J$J -W0.2p -Gwhite -O -K >> $PS

# ============绘制板块边界
gmt psxy CN-plate-neighbor.dat  -R$R -J$J -W2.0p,2/138/210 -Sf0.5+t+l -G2/138/210 -O -K >> $PS
# ============绘制推断地块边界
gmt psxy CN-block-L1-deduced.dat -R$R -J$J -W1.0p,2/138/210,- -O -K >> $PS
# ============绘制二级地块边界
gmt psxy CN-block-L2.dat -R$R -J$J -W1.0p,orange -O -K >> $PS
# ============绘制一级地块边界
gmt psxy CN-block-L1.dat -R$R -J$J -W1.0p,2/138/210 -O -K >> $PS

# ============活动地块标注
gmt pstext -R$R -J$J -F+f10p,37 -O -K >> $PS << EOF
90 45 西域活动地块
122 46 东北亚活动地块
120 35 华北活动地块
110 27 华南活动地块
90 33 青藏活动地块
97 23 滇缅活动地块
EOF

# ============边界图例
gmt pslegend -R$R -J$J -DjTL+w1.5i+jTL+o0.2c/0.2c -F+gwhite+p0.5p  -O -K --FONT_ANNOT_PRIMARY=8p,37 >> $PS << EOF
S 0.3i f0.5+t+l 0.4i/0.3c 2/138/210 2.0p,2/138/210 0.7i 板块边界
S 0.3i - 0.50i 2/138/210 1.0p,2/138/210 0.7i  一级地块边界
S 0.3i - 0.44i - 1.0p,2/138/210,- 0.7i  推断地块边界
S 0.3i - 0.50i orange 1.0p,orange 0.7i 二级地块边界
EOF

# ============南海诸岛
gmt psbasemap -R$R -J$J -DjBR+w2c+t -K -O >> $PS
gmt pscoast -R106/121/3/24 -JM2c -Ggray95 -S83/216/238 -A5000 -Dh -K -O >> $PS
gmt psxy CN-border-L1.dat -R106/121/3/24 -JM2c -B0 -W0.2p -Gwhite -O -K >> $PS

gmt psxy -R$R -J$J -T -O >> $PS
rm gmt.*
gmt gmtset MAP_ANNOT_OBLIQUE 6
gmt gmtset MAP_TICK_LENGTH_PRIMARY 0
gmt gmtset MAP_FRAME_TYPE plain

set R=78/12/149/53r
set J=B105/10/25/47/15c
set PS=BLOCK.ps

gmt psxy -R%R% -J%J% -T -K > %PS%

REM ============底图
gmt psbasemap -R%R% -J%J% -Bx10 -By5 -BNEWS -K -O >> %PS%
gmt pscoast -R%R% -J%J% -Ggray95 -S83/216/238 -A5000 -Dh -K -O>> %PS%
gmt psxy CN-border-L1.dat -R%R% -J%J% -W0.2p -Gwhite -O -K >> %PS%

REM ============绘制板块边界
gmt psxy CN-plate-neighbor.dat -R%R% -J%J% -W2.0p,2/138/210 -Sf0.5+t+l -G2/138/210 -O -K >> %PS%
REM ============绘制推断地块边界
gmt psxy CN-block-L1-deduced.dat -R%R% -J%J% -W1.0p,2/138/210,- -O -K >> %PS%
REM ============绘制二级地块边界
gmt psxy CN-block-L2.dat -R%R% -J%J% -W1.0p,orange -O -K >> %PS%
REM ============绘制一级地块边界
gmt psxy CN-block-L1.dat -R%R% -J%J% -W1.0p,2/138/210 -O -K >> %PS%

REM ============活动地块标注
echo 90 45 西域活动地块 > blockname
echo 122 46 东北亚活动地块 >> blockname
echo 120 35 华北活动地块 >> blockname
echo 110 27 华南活动地块 >> blockname
echo 90 33 青藏活动地块 >> blockname
echo 97 23 滇缅活动地块 >> blockname
gmt pstext blockname -R%R% -J%J% -F+f10p,37 -O -K >> %PS%

REM ============边界图例
echo S 0.3i f0.5+t+l 0.4i/0.3c 2/138/210 2.0p,2/138/210 0.7i 板块边界 > legend
echo S 0.3i - 0.50i 2/138/210 1.0p,2/138/210 0.7i  一级地块边界 >> legend
echo S 0.3i - 0.44i - 1.0p,2/138/210,- 0.7i  推断地块边界 >> legend
echo S 0.3i - 0.50i orange 1.0p,orange 0.7i 二级地块边界 >> legend
gmt pslegend legend -R%R% -J%J% -DjTL+w1.5i+jTL+o0.1c/0.1c -F+gwhite+p0.5p -O -K --FONT_ANNOT_PRIMARY=8p,37 >> %PS%

REM ============南海诸岛
gmt psbasemap -R%R% -J%J% -DjBR+w2c+t -K -O >> %PS%
gmt pscoast -R106/121/3/24 -JM2c -Ggray95 -S83/216/238  -A5000 -Dh -K -O >> %PS%
gmt psxy CN-border-L1.dat -R106/121/3/24 -JM2c -B0 -W0.2p -Gwhite -O -K >> %PS%

gmt psxy -R%R% -J%J% -T -O >> %PS%
del gmt.* blockname legend
pause

参考文献

  1. 张培震,邓起东,张国民,等. 中国大陆的强震活动与活动地块[J]. 中国科学D辑, 2003, 33(z1):12-20.
  2. 张国民,马宏生,王辉,等. 中国大陆活动地块边界带与强震活动[J]. 地球物理学报,2005,(03):602-610.

更新历史

  • 2017-09-17:初稿;
  • 2017-09-28:更新了 CN-block-L1.datCN-block-L1-deduced.datCN-block-L2.dat,使一、二级地块边界数据不重合;
下一页
上一页