[S23P-07] 電子基準点リアルタイム解析システム(REGARD)による2024年日向灘の地震に伴う地震時地殻変動と断層モデルの即時推定
GNSS観測によって得られる変位データに基づき推定されるマグニチュードは、原理上、飽和しない。そのため、特に巨大地震の地震規模の即時推定において優位性を持つ。国土地理院では、2012年から東北大学大学院理学研究科との共同研究の下、国内約1,300点からなるGEONET(GNSS連続観測システム)のデータをリアルタイムに解析し、座標値を推定した上で、地震前後の座標差から有限断層モデルを自動推定する、電子基準点リアルタイム解析システム(REGARD)の開発を進めている。現行のREGARDでは日本海側の観測点を固定点としてリアルタイムキネマティック解析(RTK)を常時実施し、緊急地震速報をトリガーに2種類の有限断層モデル(矩形断層モデル・プレート境界におけるすべり分布モデル)を即時推定している。このうち矩形断層モデルは、最尤推定のほかMCMCを用いた推定を試験的に実施している。REGARDは2016年から本格運用されており、これまで2016年熊本地震や2024年能登半島地震をはじめとした多くの地震において地震発生後数分で断層モデルを推定した上で関係機関へ情報提供を行い、政府内の初動対応等に活用されてきた。2024年8月8日16時43分頃に発生した日向灘を震源とする地震では、推定された地震時地殻変動と有限断層モデル情報が、当日17時30分から臨時に開催された南海トラフ沿いの地震に関する評価検討会でも活用された。本発表ではこれらの推定結果について報告する。
REGARDによって、震源に近い宮崎県南部の電子基準点で東~南東向きの水平方向の地殻変動を推定した。特に震源に最も近い電子基準点「宮崎」では、約9 cmの東南東方向の水平変動が推定された。また、上下については水平と比べて精度は劣るが、震源近傍では沈降する地殻変動が推定された。これらはGEONET定常解析(R5)による後処理解析と比較しておおむね整合する。
矩形断層モデルの推定では、最尤推定アルゴリズムによって走向が北東―南西の逆断層すべりを持つモデルが地震発生後約2分で推定された。同推定による断層モデルの地震規模はMw6.91であり、これは国土地理院による暫定の後処理モデル(Mw7.10、www.gsi.go.jp/common/000260639.pdf)と整合的である。一方で、推定されたモデルパラメータは、REGARDでは長さ約40 km・すべり角約114度が推定されたのに対し、後処理モデルでは、傾斜角20度に固定した上で、長さ約25 km・すべり角約85度が推定され、両者には差異が認められる。これは、断層推定に用いた変位場のノイズレベルの違い、特に外れ値的な変位の存在、また推定時の断層幾何の拘束の違いが影響していると考えられる。
また、MCMCによる断層モデル(Mw7.15)は地震発生後1分で形状も含めてより後処理モデルとおおむね一致する結果が得られた。加えて、観測点配置や地殻変動場のノイズを反映した推定の不確実性が可視化されている。国土地理院では、今後も断層モデル推定とその入力データとなる座標値算出の高度化を進め、REGARDの精度とシステム安定性の向上に取り組んでいく。
REGARDによって、震源に近い宮崎県南部の電子基準点で東~南東向きの水平方向の地殻変動を推定した。特に震源に最も近い電子基準点「宮崎」では、約9 cmの東南東方向の水平変動が推定された。また、上下については水平と比べて精度は劣るが、震源近傍では沈降する地殻変動が推定された。これらはGEONET定常解析(R5)による後処理解析と比較しておおむね整合する。
矩形断層モデルの推定では、最尤推定アルゴリズムによって走向が北東―南西の逆断層すべりを持つモデルが地震発生後約2分で推定された。同推定による断層モデルの地震規模はMw6.91であり、これは国土地理院による暫定の後処理モデル(Mw7.10、www.gsi.go.jp/common/000260639.pdf)と整合的である。一方で、推定されたモデルパラメータは、REGARDでは長さ約40 km・すべり角約114度が推定されたのに対し、後処理モデルでは、傾斜角20度に固定した上で、長さ約25 km・すべり角約85度が推定され、両者には差異が認められる。これは、断層推定に用いた変位場のノイズレベルの違い、特に外れ値的な変位の存在、また推定時の断層幾何の拘束の違いが影響していると考えられる。
また、MCMCによる断層モデル(Mw7.15)は地震発生後1分で形状も含めてより後処理モデルとおおむね一致する結果が得られた。加えて、観測点配置や地殻変動場のノイズを反映した推定の不確実性が可視化されている。国土地理院では、今後も断層モデル推定とその入力データとなる座標値算出の高度化を進め、REGARDの精度とシステム安定性の向上に取り組んでいく。