global_tectonics: 全球地质构造数据

global_tectonics 是一个较新的全球地质构造数据集,提供 shp, gmt, KML 格式数据的下载。

使用示例

数据文件 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 的头段记录中包含了每个板块的名字、属性、性质、面积等数据。使用 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 的头段记录中包含了不同地质区域的名字、属性、地质历史事件、面积等数据。使用 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

引用信息

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