第1章 地震
３． 震源の分布（付録1）
世界の震源の分布を自分でつくってみる。なお、この方法は、京都大学理学研究科・地球惑星科学専攻・地球物理学分野 久家慶子氏のページ（2008年11月16日現在つながりません）で公開されていたものであり、氏の了解を得て引用している。なお、この方法のオリジナルは同大学環境学研究科の酒井 敏氏とのことである。また、ある程度のＥｘｃｅｌの知識も必要である。以下の手順でできると思うが、どうしてもできないときは一番下のＥｘｃｅｌファイルをダウンロードして実際に見ることもできる。
| （１） | データの入手とExcel形式で保存 |
| ａ．データの入手 | |
| ｂ．Excel形式で保存 | |
| （２） | 震源分布図 |
| ａ．グラフ化 | |
| ｂ.．ＸとＹの入れ替え | |
| ｃ．整形 | |
| ｄ．世界地図の挿入 |
ａ．データの入手
USGS（アメリカ地質調査所）データを公開しているので、それを入手する。下の順にアクセスする。もちろん、いきなり一番下にアクセスしてもよい。
http://earthquake.usgs.gov/ USGS 国立地震防災プログラム
http://earthquake.usgs.gov/research/index.php?areaID=13 地震データ（Seismic
Data）
http://neic.usgs.gov/neis/epic/ 地震カタログサーチ（Earthquake
Catalog Search）
http://neic.usgs.gov/neis/epic/epic_global.html グローバル（世界）
最後のページで、データ形式（1. Expanded File Format with Headers and Spacesを選択してテキストファイル形式のデータを得る）、データベース（USGS/NEIC (PDE) 1973 - Presentを選択）、さらにオプションの取得する日にちの範囲（今回は2000年1月1日〜2004年3月31日まで）、マグニチュード（4.5（最低）〜9．0（最大））、震源の深さ（0km〜100kmにする）を入力する。最後の震度はIntensityは記入する必要はない。一番下の[Submit Search]ボタンを押すと、要求したデータがダウンロードされブラウザに表示される。
|
|
| USGSのサイト | このようなデーが入手できる。 |
ブラウザの保存機能（[ファイル]→[名前を付けて保存]）で、テキストファイル形式（例えばshingenbunpu.txtという名を付ける）で保存する。これを Microsoft Excel で読み込む。読み込む際、まず読み込もうとするファイル形式を＜すべて（*.*）＞にしないとファイルが見えない。また読み込むとき、「スペースによって左または右に揃えられた固定長データ」（デフォルト）、取り込み開始行（ヘッダーはいらないので23行目から）を指定し、さらにデータ区切り箇所もマウスで指定する。区切る場所は、左からデータベース名（今回はPDE）、西暦年、月、日、世界標準時、緯度、経度、震源の深さ、マグニチュードである。区切り方はExcelに表示されるので、それに従う。データ形式はとりあえず[標準」で取り込む。18,000件以上のデータがある。最後の列（K列）は必要ないので削除しておく。
1行目に2行挿入して、全体の見出し、列の見出し（PDE、年、月、日、世界標準時、緯度、経度、震源の深さ、マグニチュード、マグニチュードの説明）を付ける。見やすいように、緯度（F列、+が北緯、-が南緯）、経度（G列、+が東経、-が西経）のセルの書式を［数値型]にして小数点以下は第1位まで、震源の深さ(km）も[数値型]にして小数点以下は表示しない、マグニチュードも[数値型]にして小数点以下は第1位までにしておく。今回は時刻は使用しないのでこのままでもよいし、気になるなら[時刻型]にしてもよい。
必要と思うなら、罫線なども加えて整形する。これを今度はExcel形式で保存する。
 |
 |
| 取り込み開始行をデータが始まる行である23にする。 | データの区切り位置を確認して、きちんと上のように区切る。 |
（２） 震源分布図
ａ．グラフ化
震源の緯度と経度を選択して、震源の分布図をつくる。Excelのグラフ作成機能（ウィザード）を利用する。
範囲を選択したら、ツールバーの「グラフ」を押してウィザードを起動して[散布図]を選んでグラフ化する。グラフのタイトルは≪世界の震源分布（2000年1月1日〜2004年3月31日 M4.5以上、深さ100km以下）≫と、グラフを作成する場所は「新しいシート」を選択する。こうしてつくったグラフは下を参照。
ｂ．XとYの入れ替え
ａ．でつくったグラフは縦軸が経度、横軸が緯度になっているので、これを入れ替える。グラフ上で右クリックしてメニューを出し、「系列」タグを押して、Ｘの値（Ｘ）とＹの値（Ｙ）を入れ替える。Ｘの値（Ｘ）の指定部分の“F”を“G”に、Yの値（Y)の指定部分の“G”を“F”にする。これで何となく世界地図らしいものが浮かんでくる。
 |
|
| 左のようになっているので、Xの値とYの値を入れ替える | |
ｃ．整形
まず、Ｘ軸の整形をする。Ｘ軸上で右クリックしてメニューを出し、「目盛り」タグで軸の最小値を-180、最大値を180にする。目盛り間隔はこのままでもよいが、とりあえず60、補助目盛り間隔は30とでもしておく。「表示形式」タグで、小数点以下の表示を0にする。
次にＹ軸を整形する。適当なＹ軸の線（例えば緯度0°の線）上で右クリックしてメニューを出し、「目盛り線の書式設定」→最小値を-90、最大値を90、目盛り間隔を30、補助目盛り間隔を15にする。Ｙ軸上で右クリックして「目盛り線の書式設定」を出し、「表示形式」タグで小数点以下の表示を0とする。
震源の位置を表わす紺色の点が大きいので小さくする。どれか一つを選んで右クリックして「データ系列の書式設定」メニューを出す。マーカーの大きさを2ポイントとする。ここで色も適当に変えてもよい。とりあえず赤にしてみる。
凡例が残っていたら削除する。 これでほとんど完成である。
ｃ．の結果。世界地図らしくなってくる。
ｄ．世界地図の挿入
別に用意しておく世界地図を埋め込む。グラフ上で右クリックして、「プロットエリアの書式設定」メニューを出す。右下の「塗りつぶし効果」を押して、「図」タグを選択→「図の選択」で図を選択する。ここではこの世界地図のgif画像（8KB）を利用する。
これで完成である。下図参照。
さらにシートに名を付けておくとよい。データが読み込まれたシートを「震源データ」、グラフを描いたシートを「世界の震源分布」としておく。
最終的なExcelファイルはこちら（約3MB）になる。
