著者

加藤 翔太 西田 究

雑誌

JpGU-AGU Joint Meeting 2020

巻号頁・発行日

2020-03-13

---

地震波干渉法は2観測点で観測された地震波形記録の相互相関関数を計算することにより、片方を仮想的な震源とし、もう片方を観測点とした場合の観測波形を推定する手法である(e.g. Snieder et al., 2013)。地震波干渉法解析では、地震波動場がランダムかつその強度分布が等方・均質であることを仮定する。ランダムな波動場として海洋波浪起源の脈動を解析に用いる場合には、周期5-20 sの帯域で表面波が卓越することが知られている。そのため、脈動を用いた地震波干渉法は地殻・上部マントルの3次元構造の推定に適している(e.g. Shapiro et al., 2005)。近年では、表面波だけではなく実体波の抽出が試みられている。その一例として、マントルの410/660 km不連続面からの反射P波の抽出が報告されている(Poli et al., 2012, Feng et al., 2013)。しかし、これらの反射P波を抽出した先行研究の対象地域は大陸に限られていた。本研究の目的は、防災科学技術研究所Hi-netの上下動記録に地震波干渉法を適用することにより深さ410/660 km不連続面からの反射P波を抽出し、日本列島下の不連続面をイメージングすることである。本研究では以下の手順で各観測点ペアに対する相互相関関数を計算した。用いた波形記録は防災科学技術研究所Hi-net観測点のうち西南日本に存在する240点の上下動記録(2007年-2018年)である。まず、Hi-netの上下動記録を2 Hzにダウンサンプリングした。その上で各観測点について翌日の観測波形との差を計算して元の観測波形の代わりに用いた(高木ほか、2019)。これは、Hi-netの機器ノイズ(Takagi et al., 2015)の相互相関関数への影響を抑えるためである。次に、得られた1日長の波形を1024 sの時間窓に分割し、周期5-10 sおよび10-20 sの平均2乗振幅によって時間窓を選択した。選択した時間窓について周波数領域で白色化を行い、周期1-10 sの成分について全観測点ペアの相互相関関数を計算した。まず4-th root vespagramを全観測点ペアに対する相互相関関数について計算した(Rost and Thomas 2002)。その結果、410 km不連続面の反射P波がオフセット距離0-300 kmで見られ、660 km不連続面の反射P波はオフセット距離50-100 kmで見られた。また、660 km不連続面の反射P波は410 km不連続面の反射P波に比べて弱いことがわかった。次に、得られた反射P波を不連続面の深度に変換するため、Common Middle Point (CMP)重合を行った(e.g. Stein and Wysession, 2003)。具体的には、オフセット距離が500 km以内の各観測点ペアについて反射点の位置でグループ分けを行い、各グループに対して不連続面が水平と仮定しCMP重合を行った。速度構造はJMA2001(上野ほか、2002)を用い、深度推定は410 km不連続面についてのみ行った。その結果、地域ごとに410 km不連続面深度の変動が見られ、特に東経134°-135°北緯33°-36°に反射点を持つグループでは不連続面が上昇している結果が得られた。これは、従来の地震波を用いた不連続面深度に関する研究(Tonegawa et al., 2005, Tono et al., 2005)と整合的である。本研究では地震波干渉法により西南日本の410 km不連続面深度の推定を行った。今後は用いる観測点を増やして対象領域を日本全国へと拡大するとともに、今回扱わなかった660 km不連続面の推定も行う予定である。謝辞:本研究では防災科学技術研究所のHi-netの上下動記録を用いました。記して感謝いたします。