そこで本研究ではその余震活動の空間分布と海底地下構造との関係をより明確にす るために、2004年10月に震源域付近で海底地震計を用いた人工地震探査をおこなっ た。得られた記録から地震波速度構造を求め、余震分布やアスペリティ分布との比 較を行なうことが目的である。

観測では12台の海底地震計を用い、20lのエアガン2台を人工震源として屈折法地震 探査をおこなった。海底地震計は、本震付近を通り海溝と直交する測線上に設置さ れ、屈折法地震探査と同時に、50チャンネルのハイドロフォン・ストリーマーを用 いて反射波の観測もおこなった。しかし観測機器の不備のため反射波の記録を得る ことはできなかった。

回収された海底地震計のうち記録のとれていたもの11台分のデータを使用して波線 追跡法(ZELT ラベル)によるフォワード・モデリングをおこなった後、得られたモ デルを初期モデルとして2次元走時インバージョン(ZELT ラベル)をおこなった。ま た、浅部堆積層の速度構造を拘束するため、海洋研究開発機構(JAMSTEC)が2000年 に同一測線上でおこなったマルチチャンネル反射法地震探査の記録を解析に利用し た。

解析の結果深さ約15kmまでのP波速度構造を求めることができた。しかしスラブの 15kmよりも深い部分からの波を読み取ることができなかったため、本震付近でのプ レート境界の構造を得るに至らなかった。本研究で得られた速度構造モデルは、 2003年十勝沖地震の余震分布や周辺地域でおこなわれた過去の構造探査の結果と調 和的なものであった。