地下の速度構造

茨城県南西部の地下の速度構造を地震波トモグラフィ法を用いて推定しました．
図１．茨城県南西部の鉛直断面．

(a) 気象庁カタログの震源データ． (b) S波速度の分布． (c) Vp/Vs比の分布． (d) 発震機構のP軸（橙）とT軸（青）の分布． (b)～(d)のドットはDDトモグラフィ法で再決定した震源分布． (b)～(d)の赤線は本研究で推定したプレート境界，破線は厚さ7 kmを仮定して推定したスラブモホ面．

図１から，次のことがわかります．

• 震源再決定によって，震源分布が締まった．
• 赤実線に沿う地震は，P軸・T軸分布から低角逆断層型であることがわかる．
  つまり赤実線はプレート境界を示すことになる．
• 厚さ7 km程度の低Vsと高Vp/Vs域が赤実線の下に分布している．これはスラブ地殻の厚さと整合する．

また，ここで示した茨城県南西部以外の地域でも

• 厚さ7 km程度の低Vsと高Vp/Vs域が分布している．
• それらは人工地震探査や別の解析手法によって得られた結果とも整合する．

これらの情報から，厚さ7 km程度の低Vsと高Vp/Vs域の上面がプレート境界であると考えられます． 関東から九州南部の広範囲に対して地震波トモグラフィ法を適用し，地下構造の推定およびプレート形状の推定を行いました．

推定したプレート形状

図２．フィリピン海プレートおよび太平洋プレート上面のコンター．

赤および紫線はフィリピン海プレート，橙破線は太平洋プレートの上面の形状を示します． 矢印は陸のプレート（ユーラシア，アムール，およびオホーツクプレート）に対するフィリピン海プレートおよび太平洋プレートの沈み込みベクトルです． 太平洋プレートは日本列島の東沖の海溝（千島海溝，日本海溝，伊豆・小笠原海溝）から沈み込み， 一方，フィリピン海プレートは日本列島の南沖のトラフ・海溝（相模トラフ，駿河トラフ，南海トラフ，琉球海溝）から沈み込んでいます．

関東では陸のプレートと太平洋プレートの間にフィリピン海プレートが押し入るように沈み込んでいるため，とても複雑な場となっています． 関東の青破線はフィリピン海プレートの北東端を示します． 太平洋プレートが既に居座っているため，ここまでしかフィリピン海プレートは沈み込んでいないと考えられています．

太平洋プレートに比べ，フィリピン海プレートの形状はぐにゃぐにゃしていることが見て取れます． これは主にプレートの年齢の違いによるものでしょう．太平洋プレートに比べて若く温かいフィリピン海プレートは変形しやすいのです． 関東の青領域はフィリピン海スラブと太平洋スラブの接触域を示します． 太平洋プレートが東から押しているため，接触域を介して， 変形しやすいフィリピン海プレートが東西に皺の寄った形状となっているのではないかと考えています．

フィリピン海プレート，太平洋プレートに関するページもご参照ください．

まとめ・今後の課題

地震波トモグラフィ法を用いて速度構造や震源位置を推定し，厚さ7 km程度の低Vsと高Vp/Vs域の情報からプレート形状を推定しました． 地震は主にプレート運動に起因します．地震を考える上で，地下の速度構造やプレート形状の情報は必要不可欠です．

ここでの解析には自然地震を用いているため，真の発震時間がわかりません．そのため結果にはどうしても誤差が多く含まれます． 一方，人工地震を用いれば発震時間は明確なため伝播時間に関する誤差は少ないですが， 発震場所が限られるため全国の地下構造を面的に推定することが困難です．また莫大な費用がかかります． このように一長一短ありますが，計測機器の進歩，データの蓄積や解析手法の開発などによって， 地下の様子は徐々に明らかになるでしょう．