event、以下「L-SSE」という。)の解析を行いました。その結果、L-SSE域の上限および内部にすべり量の急変を見出すとともに、これまでの手法よりもすべり域内部の構造を詳細に推定することができました。急変位置は、地震発生帯の下限および深部低周波微動の上限とよく一致することを示しました。これは、従来なめらかな空間変化しか捉えられないと考えられてきた地殻変動観測データに対する解析で初めて見出された結果であり、今後の地殻変動解析研究の解像度向上にも大きく貢献するものです。適切に推定できることになります。様々な地震現象の発生メカニズムに関する統一的な理解にもつながります。この研究はJSPS科研費JP16H01564、JP25120005、JP25120009、JP25120011、JP16K17798、JP16K16108、JP25280090、JP15K20864、JSTさきがけJPMJPR1676からの助成を受けました。Reports」に7月21日付け(日本時間)で掲載予定です。フィリピン海プレートが沈み込んでいます。kmのプレート境界面では、約1年かけてゆっくりとすべるスロー地震(ゆっくり地震)(※2)の一種であり、約5~7年に一度、間欠的に発生するL-SSEと呼ばれる現象が観測されています。このL-SSE領域より浅部側は、1946年南海地震(M8.0)の西端や1968年日向灘地震(M7.5)の震源域(地震発生帯)となっており、将来の南海トラフ地震震源域の西限に関係する重要な場所です。slip、以下「ETS」という。)」と呼ばれるスロー地震が起きています。これら一連の現象は、全て同一のプレート境界面上で発生していると考えられています(図1)。プレート境界面上のL-SSEやETSの発生域では、通常の地震の発生頻度は高くありません。従って、浅部における大地震の発生メカニズムや準備状況を理解するためには、豊後水道におけるL-SSEの発生領域を正確に把握することが非常に重要です。類似した現象が観測されつつあります。このようなスロー地震現象のいくつかは、日本海溝や北米カスカディアなど、過去にマグニチュード8以上の海溝型地震を経験している沈み込み帯でも起きています。その中で豊後水道周辺は、比較的観測網が充実していることもあり、最も多様な現象が繰り返し確認されている地域です。豊後水道周辺で発生した3つのL-SSE(1997年、2003年および2010年)は、すでにすべり分布の解析が行われており、3回ともほぼ同じ中心から外側へ向かって緩やかに変化する分布が推定されています。すべり域は、深部低周波微動や地震発生帯まで連続的に広がっているようにも見えました。従来の手法では、地殻変動観測データの不十分さを補うためにすべり分布のなめらかさを前提として解析していることから、すべり域内部が本当になめらかな分布をしているかは、明らかではありませんでした。図1.豊後水道周辺の概略図。成果そこで本研究グループは、地殻変動観測データの不十分さを補い、知りたい情報を正確に抽出することが可能となるスパースモデリングと呼ばれる手法を適用し、L-SSEの詳細なすべり分布を推定しました。その結果、0.2m以上すべった領域の浅部側では、すべり量が0mまで急に減少しており、従来考えられていたよりも急な空間変化であることが示されました(図2)。0.1m以上すべった領域の内部(0.2m以上すべった領域の深部側)でも、すべり量が半分近くまで急に減少する部分があることが明らかになりました。上限とよく一致していました(図3・図4)。まで緩やかに広がっているのではなく、地震・L-SSE・深部低周波微動における発生領域の空間的なすみわけが、これまでよりも明瞭に見えるようになりました。本結果が解析手法や観測点分布に起因するものではないことを確認するため、従来のなめらかな分布および今回得られたすべりの急な変化を模擬データとして同様に解析したところ、本結果を支持する空間変化をほぼ再現することができました。スパースモデリングの深化と高次元データ駆動科学の創成」の公募研究課題「スパースモデリングによる地震発生予測のための地殻活動データからの情報抽出」として、地球科学と情報科学の異分野間の学融合共同研究によってなされたもので、分野によらない普遍的なデータ解析の方法論を探求するデータ駆動科学(※3)がもたらした画期的な成果といえます。図2.本研究で得られた2010年のL-SSEによるすべり分布。0.1m以上すべった領域を薄い赤、0.2m以上すべった領域を濃い赤で示している。○はすべりを推定したプレート境界面上の位置を地表に投影したもの。●は解析に使用した地殻変動観測点の位置。(b)図2(a)の灰色線ABに沿ったすべり量の分布(赤線)。従来の手法で得られたもの(灰線)。図3.本研究で得られた2010年L-SSEのすべり分布(赤線)と、他の地震学的知見との空間分布の比較。沈み込むプレート境界面での等温度線(点線)。図4.豊後水道周辺の概略図。発生領域のすみわけを見出すことができました。した詳細な空間分布をすべりの力学モデルにより再現することは、L-SSE領域の深部側で発生しているETSから浅部側で起きている地震発生帯に至るまで、異なる時間的特徴を持った様々なすべりの発生メカニズムの系統的な理解につながります。今後、本研究グループでは、引き続き豊後水道L-SSEの解析を進めるとともに、紀伊水道や東海地方の地殻変動観測データに対しても同手法を適用し、南海トラフ全域の地震発生サイクルシミュレーションを高度化していきます。余効すべりを解析し、同域の応力が集中する場所についても検討を進める予定です。これらを進めることにより、地震・スロー地震の発生メカニズムの解明に寄与することが期待されます。様々な背景を持つ研究者同士の密接な協働によって得られた成果です。近年、データ駆動科学と呼ばれる学融合的な新学術分野が注目を集めています。ここでは、通常は別々に研究することしかなかった様々な分野の自然科学者と情報科学者が共同研究を進めています。本研究では、医学や天文学で成果を挙げているスパースモデリングが、そのまま地球科学でも威力を発揮しました。このように、異分野の融合を推し進め、普遍的なデータ解析手法の開発と応用を実践するデータ駆動科学は、今後ますますその重要性を増していくことと考えられます。Email:k-somu@adm.k.u-tokyo.ac.jp用語解説(※1)スパースモデリング高次元データから本質的な低次元空間を抽出する数理的手法の総称。自然科学データから未知の科学的知見を抽出する革新的なデータ解析技術として、様々な自然科学分野で注目されつつある。スパースとは「少ない」や「疎」という意味。スパースモデリングでは、データへの適合性と同時に、拘束条件としてスパース性(変数の少なさなど)を取り入れることで、質・量の不足したデータからでも、知りたい情報を正確に抽出することができるため、医学や天文学など様々な観測・計測分野の解析に応用されつつある。元素含有量による津波堆積物の高精度識別手法の開発などが挙げられる。最近では、沖合の水圧計データから沿岸の津波高を瞬時に予測する手法の高精度化などにも応用されている。lasso)を適用した。スロースリップイベント(ゆっくりすべり)など、通常の地震よりもゆっくりとした断層すべりの総称。低周波微動は、通常の地震よりも低周波数成分の卓越した地震波を放出する現象。スロースリップイベントは、地震波を放出せずに、地下の断層がおよそ一日以上の期間をかけてゆっくりとすべる地殻変動現象である。そのうち、継続期間が数か月以上のものを長期的スロースリップイベントと呼ぶ。(※3)データ駆動科学データ駆動科学とは、観測・計測データを最大限に活用し科学的知見を導出するために、対象や分野によらない普遍的なデータ解析の方法論を探求する学融合的な分野である。最先端のデータ解析手法を共有することにより、異分野間の学問的交流が促進され、大きな相乗効果が期待される。