✍️ 陈箫翰 , 田冬冬 • 📅 2025-07-30

# global_tectonics: 全球地质构造数据

global_tectonics 是一个全球地质构造数据集，提供 Shapefile、GMT、KML 格式数据的下载。

- 官方网站: https://github.com/dhasterok/global_tectonics

## 数据简介与GMT格式下载

global_tectonics 数据集中提供了四个数据文件：

- 板块边界数据 : `boundaries.gmt`

- 板块数据 : `plates.gmt`

- 地质块体数据 : `global_gprv.gmt`

- 海陆边界线 : `oc_boundaries.gmt`

更详细的数据说明请阅读数据的参考文献。

## 使用示例

### boundaries.gmt

数据文件 `boundaries.gmt` 包含了 7 种类型的板块边界数据。每段数据的头段记录中 包含了板块边界的名称、板块边界类型、相邻的两个板块名以及长度等信息。 下面的示例使用 convert 模块的 -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` 包含了各个板块的多边形数据。每段数据的头段记录中包含了板块名字、属性、性质、面积等数据。 使用 convert 模块可以提取出特定板块的数据进行绘图。 下面的例子展示了如何根据头段记录提取大陆地壳板块、变形带板块，以及提取青藏高原板块、 华北板块。

    #!/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` 的头段记录中包含了不同地质区域的名字、属性、地质历史事件、 面积等数据。使用 convert 模块可以提取出特定地质区域的数据进行绘图。 下面的例子展示了如何根据从头段记录 `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