global_tectonics: 全球地质构造数据

撰写:

陈箫翰, 田冬冬

最近更新日期:

2023-10-04


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
../../_images/c3dbc38fe29161bcd14fd3162246efda.png

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
../../_images/3933076f1b38d27070b3f3b2696ac6ca.png

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
../../_images/c0d0c8c385f8eabbdc3d8a7856eae0f6.png

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
../../_images/5d2473d19f99c8e3614b022d1be8c612.png

引用信息

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