世界の（過去10年間の）地震分布

　以前、日本の過去10年間の大地震（「最大震度６強以上」）についてプロットしたが、今回は世界の過去10年間の大地震についてプロットしてみた。

　日本の地震は気象庁のデータを利用したが、
　世界の地震はUSGS（アメリカ地質調査所）のデータを利用。
　USGSは日本の国土地理院のように地形図を作成しているところ。

目　次

マグニチュード６以上
マグニチュード５以上

マグニチュード６以上

　まずUSGS　Earthquake Hazards Program（earthquake.usgs.gov/）にアクセス。
　Latest Earthquakesで、直近24時間に起こったマグニチュード2.5以上（M2.5＋）の地震震源情報、地図表示。

　Search Earthquake Catalogで、規模、期間、地域を選択して表示。
　今回、

　規模　＝　マグニチュード６以上（M6＋）　・・・　CustomのMinimumを６
　期間　＝　Start　2011-01-01～End　2021-05-30
　地域　＝　World（世界）

を選択して検索　Search。

　検索結果　Search Resaltsは、1516データ。

　画面左のリストの一番下からダウンロード　DOWNLOADできる。
　データ形式はCSV、GeoJson、KML、QuakeML。

　Leaflet　＋　OpenStreetMapでプロットされているので、改めてプロットし直す必要はないが、せっかくなのでGeoJSONデータをダウンロードしてプロット。

2011年以降の世界の地震（「M6＋」）震源プロット　　　→
　マグニチュード７以上（M7＋）が赤丸、詳細はUSGSのページへ
　追）プレート境界も追加
　　　データのtimeはUNIX時間

　気象庁のマグニチュードとUSGSのマグニチュードはやや異なる。

マグニチュード５以上

　もっとたくさんプロットすべくマグニチュード4.5以上（M4.5＋）で検索したら極端に増えて78152データ。
　20000データが限度。
　マグニチュード５以上（M5＋）で検索したら限度内の18446データ。

　マグニチュード　Mと地震の発生頻度　Nの関係式は、
　log₁₀ N ＝ a － bM
　グーテンベルク－リヒターの関係式（Ｇ－Ｒ式）。印刷のグ氏とは別人。

　この式は以前こちらのページで触れたマグニチュード　Mと地震のエネルギー　Eの関係式
　1.5M ＋ 4.8 ＝ log₁₀ E　[J]
と同じ片対数の形式。傾きが違う。

　Mが増減　⇔　Eが指数関数的増減
　Mが増減　⇔　Nが指数関数的減増
という式。

　log₁₀ N’ ＝ a － bM’
　log₁₀ (N’／N) ＝ log₁₀ N’ － log₁₀ N ＝ b(M－M’)　・・・　aが消える
　N’／N ＝ 10^b(M－M’)
で、b値は0.9～1.0、とのこと。
　よってM－M’＝1、つまりMが１小さくなると発生頻度10^b倍（約８～10倍）、Mが２小さくなると10^2b倍（約63～100倍）。
　Mが1.5小さくなると10^1.5b倍（約22～約32倍）。