global_tectonics: 全球地质构造数据
global_tectonics 是一个全球地质构造数据集,提供 Shapefile、GMT、KML 格式数据的下载。
数据简介
global_tectonics 数据集中提供了四个数据文件:
boundaries.gmt: 板块边界数据
plates.gmt: 板块数据
global_gprv.gmt: 地质块体数据
oc_boundaries.gmt: 海陆边界线
更详细的数据说明请阅读数据的参考文献。
使用示例
boundaries.gmt
数据文件 boundaries.gmt
包含了 7 种类型的板块边界数据。每段数据的头段记录中包含了板块边界的名称、板块边界类型、相邻的两个板块名以及长度等信息。下面的示例使用 gmtconvert 模块的 -S 选项提取不同类型的板块边界并分别绘制:
#!/usr/bin/env bash
gmt begin boundaries
gmt basemap -JN15c -Rg -Ba60
gmt coast -G244/243/239 -S167/194/223
# 提取数据头段中不同类型的边界数据并绘图
gmt convert boundaries.gmt -S"type=spreading center" | gmt plot -W0.5p,orange -l"spreading center"
gmt convert boundaries.gmt -S"type=extension zone" | gmt plot -W0.5p,red -l"extension zone"
gmt convert boundaries.gmt -S"type=subduction zone" | gmt plot -W0.5p,black -l"subduction zone"
gmt convert boundaries.gmt -S"type=collision zone" | gmt plot -W0.5p,darkgreen -l"collision zone"
gmt convert boundaries.gmt -S"type=dextral transform" | gmt plot -W0.5p,blue -l"dextral transform"
gmt convert boundaries.gmt -S"type=sinistral transform" | gmt plot -W0.5p,pink -l"sinistral transform"
gmt convert boundaries.gmt -S"type=inferred" | gmt plot -W0.5p,gray -l"inferred"
# 设置图例属性
gmt legend -DJMR+o0.1c -F+p1p+gwhite
gmt end show
plates.gmt
数据文件 plates.gmt
包含了各个板块的多边形数据。每段数据的头段记录中包含了板块名字、属性、性质、面积等数据。使用 gmtconvert 模块可以提取出特定板块的数据进行绘图。下面的例子展示了如何根据头段记录提取大陆地壳板块、变形带板块,以及提取青藏高原板块、华北板块。
#!/usr/bin/env bash
gmt begin plates
gmt basemap -JN15c -Rg -Ba60
gmt coast -Gwhite -S167/194/223
# 从crust_type头段中提取出所有大陆地壳板块,并填充橘黄色,设置透明度80%
gmt convert plates.gmt -S"crust_type=continental" | gmt plot -Gorange@80 -l"continental crust plates"
# 从plate_type头段中提取出所有变形带板块,并填充红色,设置透明度50%
gmt convert plates.gmt -S"plate_type=deformation zone" | gmt plot -W0p -Gred@50 -L -l"deformation zone"
# 从poly_name头段中提取出青藏高原板块、华北板块,并分别填充蓝色和淡绿色
gmt convert plates.gmt -S"poly_name=Tibet" | gmt plot -W0.5p -L -Gblue -l"Tibet plate"
gmt convert plates.gmt -S"poly_name=North China" | gmt plot -W0.5p -Glightgreen -l"North China plate"
gmt end show
global_gprv.gmt
数据文件 global_gprv.gmt
的头段记录中包含了不同地质区域的名字、属性、地质历史事件、面积等数据。使用 gmtconvert 模块可以提取出特定地质区域的数据进行绘图。下面的例子展示了如何根据从头段记录 prov_type
提取克拉通、地盾、被动边缘,从 lastorogen
提取阿尔卑斯-喜马拉雅造山事件的区域,以及从 prov_name
中提取南海盆地和鄂尔多斯地块的数据。
#!/usr/bin/env bash
gmt begin global_gprv
gmt basemap -JD105/35/36/42/10c -R70/140/3/60 -Baf
gmt coast -Gwhite -S167/194/223
# 提取出所有克拉通、地盾、被动板块边缘的边界数据,并填充颜色,设置30%透明度
gmt convert global_gprv.gmt -S"prov_type=craton" | gmt plot -G168/171/210 -t30 -l"craton"
gmt convert global_gprv.gmt -S"prov_type=shield" | gmt plot -G217/205/228 -t30 -l"shield"
gmt convert global_gprv.gmt -S"prov_type=passive margin" | gmt plot -G196/233/236 -t30 -l"passive margin"
# 提取出所有和阿尔卑斯-喜马拉雅造山事件相关的地质区域边界数据,并填充红色,设置30%透明度
gmt convert global_gprv.gmt -S"lastorogen=Alpine-Himalayan" | gmt plot -Gred@30 -l"Alpine-Himalayan"
# 提取出南海盆地的边界数据,绘制出黑色轮廓;提取出鄂尔多斯地块的边界数据,绘制出橘黄色轮廓
gmt convert global_gprv.gmt -S"prov_name=South China Sea Basins" | gmt plot -W1.5p -l"South China Sea Basins"
gmt convert global_gprv.gmt -S"prov_name=Ordos Block" | gmt plot -W1.5p,orange -l"Ordos Block"
gmt end show
oc_boundaries.gmt
数据文件 oc_boundaries.gmt 包含了全球的海陆边界线。
#!/usr/bin/env bash
gmt begin oc_boundaries
gmt basemap -JN15c -Rg -Ba60
gmt coast -G244/243/239 -S167/194/223
gmt plot oc_boundaries.gmt -W0.5p,red
gmt end show
引用信息
Hasterok, D., Halpin, J. A., Collins, A. S., Hand, M., Kreemer, C., Gard, M. G., & Glorie, S. (2022). New Maps of Global Geological Provinces and Tectonic Plates. Earth-Science Reviews, 231, 104069. https://doi.org/10.1016/j.earscirev.2022.104069