【適宜更新】最近作成した図たち

ここでは最近GMTで作成した図を紹介します。

2024/4/28 オクラホマ

図

スクリプト

1. 震源

   #!/bin/bash
   gmt begin oklahoma_hypo png
   gmt basemap -JM12 -R-103/-93/33/41 -Bafg -BWSNE+t"Oklahoma_2001-2023_USGS"
   gmt coast -Df -W0.25 -Glightgray -Slightcyan -LJBR+jTR+o0/1.5+w100+f -N2
   gmt makecpt -Cseis -T0/40/1 -Z
   awk -F"," '{print $3, $2, $4, $5*0.1}' query.csv | gmt plot -Sc -C -W0.25
   gmt colorbar -DJBR+jBL+o1.5/0+w-5/0.2 -Ba10f10+l"Depth (km)"
   gmt inset begin -DjBR+w2.5/1.4+o0.1c -F+gwhite+p0.25 -M0.01
       gmt set MAP_FRAME_TYPE = "plain"
       gmt basemap -JM2.5 -R-125/-65/24.6/50 -Bafg -BWSNE
       gmt coast -Df -W0.25 -Glightgray -Slightcyan -LJBR+jTR+o0/1.5+w100+f -N2
       gmt plot -W0.25,blue -L <<END
   -103 33
   -93 33
   -93 41
   -103 41
   END
   gmt inset end
   gmt basemap -JX12/-12 -R-103/-93/0/40 -Bxafg+l"longitude" -Byafg+l"depth" -BWSne -Y-15
   awk -F"," '{print $3, $4, $4, $5*0.1}'  query.csv | gmt plot -Sc -W0.01 -C
   gmt colorbar -DJBR+jBL+o1.5/0+w-5/0.2 -Ba10f10+l"Depth (km)"
   gmt end

2. N-T図

   #!/bin/bash
   input="query2.csv"
   # ${input}のうち2001から2023が含まれる行の数を数える
   gmt begin count png
   awk -F"," '{print $1}' ${input} | gmt histogram -JX12 -R2001/2023/0/4000 -Bxa1fg+a45+l"Year" -By1000a+l"Count" -BWSne+t"Oklahoma_2001-2023_USGS_M2.5+"  -Glightgray -W -D -T1
   gmt inset begin -DjTL+w5/2.8+o0.1c -F+gwhite+p0.25 -M0.01
       gmt set MAP_FRAME_TYPE = "plain"
       gmt basemap -JM5 -R-125/-65/24.6/50 -Bafg -BWSNE
       gmt coast -Df -W0.25 -Glightgray -Slightcyan -LJBR+jTR+o0/1.5+w100+f -N2
       gmt plot -W3.5,blue -L <<END
   -103 33
   -93 33
   -93 41
   -103 41
   END
   gmt inset end
   gmt end

   謝辞

   データ出典：USGS
   記して感謝申し上げます。USGSはデータ取得しやすい！便利！！

2024/4/17 豊後水道

図

スクリプト

1. 震源

   #!/bin/bash
   input="./h2023.csv"
   gmt begin bungo_suido png
   gmt basemap -JM12 -R132.2/132.6/33.0/33.4 -Bafg -BWSNE+t"2023/01/01-2023/12/31_JMA"
   gmt coast -Df -W0.25 -Glightgray -Slightcyan
   gmt makecpt -Cseis -T0/50/1 -Z
   awk -F"," '{print $11,$9,$13,$15*0.1}' ${input} | gmt plot -Sc -W0.25 -C
   gmt plot -Sa -W0.25 -C << END
   132.4 33.2 41 0.62
   END
   gmt colorbar -DJBR+jBL+o1.5/0+w-5/0.2 -Baf+l"depth"
   gmt basemap -JX12/-12 -R132.3/132.5/0/50 -Bxafg+l"longitude" -Byafg+l"depth" -BWSne -Y-15
   awk -F"," '{print $11,$13,$13,$15*0.1}' ${input} | gmt plot -Sc -W0.01 -C
   gmt plot -Sa -W0.25 -C -l"202404172314" << END
   132.4 41 41 0.65
   END
   
   gmt colorbar -DJBR+jBL+o2.5/0+w-5/0.2 -Baf+l"Depth"
   gmt end

2. 発震機構解

   #!/bin/bash
   # ローカルのファイル先を指定したりURLを入力してください。
   input="./mecha/mecha2011-2023.txt"
   rm -f mecha.txt mecha_aki.txt
   # データ収集（URLから取得する場合のみ）
   # wget $url -O $input
   ## 安芸・リチャーズ
   ## 1,2行目は削除してカンマ区切りにする
   sed -e 1,2d $input | sed -e 's/  */,/g' | awk -F"," '{print $4,$3,$5,$8,$9,$10,$6}' > mecha.txt
   ## 負の値の対応
   sed -e 's/;[0-9]*//g' mecha.txt | sed -e 's/-\([0-9]*\)-\([0-9]*\)/-\1/g' > mecha_aki.txt
   # GMT
   gmt begin bungo_suido_mecha png
   gmt basemap -JM12 -R132.2/132.6/33.0/33.4 -Bafg -BWSNE+t"2011-2023_F-net"
   gmt coast -Df -W0.25 -Glightgray -Slightcyan -LJBR+jTR+o0/1+c40+w10+f
   gmt makecpt -Cseis -T0/50/1 -Z
   gmt meca "mecha_aki.txt" -Sa0.5 -C
   gmt meca -Sa0.5 -C << END
   132.4 33.2 41 166.1 69.6 -114 6.2 0 0 202404172314
   END
   gmt colorbar -DJBR+jBL+o2/0+w-5/0.2 -Baf+l"Depth(km)"
   gmt end

注意事項

• 震源の位置はAQUAシステム メカニズム解カタログに基づくため今後変わる可能性があります。
• 2023年の気象庁一元化震源データは今後見直される可能性があります。

謝辞

以下のデータを使用させていただきました。ここに記して感謝申し上げます。

• 震源データとして、気象庁の「地震カタログ」を一部改変したもの。

気象庁｜震源データ

気象庁が提供するページです

www.data.jma.go.jp

• 発震機構解のデータとして、国立研究開発法人防災科学技術研究所のF-netデータ

防災科研機関リポジトリ

nied-repo.bosai.go.jp

• HinetPy

Dongdong Tian, 2020 , HinetPy: A Python package to request and process seismic waveform data from Hi-net. , doi: 10.5281/zenodo.3885779