Fractional-Order Viscoelastic Modeling of Seismic P-, S-, and Rayleigh-Wave Attenuation: Marmara Seismotectonic Context, Fault Geometry, and a Reproducible Strong-Motion Framework

地球の表面深くでは、地面は単純な固体の岩石の塊としては振る舞いません。地震が発生すると、エネルギーは波として外側へと伝わっていきますが、これらの波が地殻を移動するにつれて、エネルギーを失い、その形を変えていきます。減衰として知られるこのプロセスこそが、震源から離れた場所での揺れの感じ方が、破壊の直近での感じ方と異なる理由です。何十年もの間、科学者たちは標準的な物理学の規則を用いて、このエネルギーの損失を記述しようとしてきました。そこでは、地球の地殻を、完全に弾性のあるスプリング、あるいは単純で粘り気のある流体として扱ってきました。しかし、現実の世界は、これら二つの極端なモデルが許容するものよりも、はるかに複雑です。地面は岩石、流体、そして亀裂が混ざり合った混沌とした混合物であり、標準的な物理学では捉えきれない方法で、過去の動きを「記憶」しているのです。ここで、分数微積分学と呼ばれる数学の一分野が有用となります。物質の反応を記述するために整数を用いる代わりに、分数微積分学は、その中間的な値を用いることを可能にし、より少ない仮定で岩石の「記憶」をモデル化する方法を提供します。

イスタンブールのすぐ南に位置する、激しい地質活動が見られるマルマラ海に焦点を当てた新しい研究において、タヤン・デミルとアタカン・コチイギットの研究者らは、地震波がどのように減衰していくかを記述するための、より精密な数学的枠組みを構築しました。マルマラ海は、テクトニクスプレートが擦れ合う巨大な地殻の割れ目である北アナトリア断層の上に位置しています。この断層は主要都市の真下を走っているため、この特定の領域において地震波がどのように伝わり、弱まっていくのかを正確に理解することは、より安全な建物を設計する上で極めて重要です。研究者たちは、岩石の複雑で履歴依存の挙動を考慮に入れた分数モデルが、伝統的なモデルよりも地震波の減衰をより良く説明できるかどうかを検証することを目指しました。彼らは次にいつ大きな地震が起こるかを予測しようとしたのではなく、事象が記録された後に発生する揺れを、いかに正確に再構成するかという点に焦点を当てました。

チームはまず、地質学的な舞台の設定から始めました。彼らは、テキルダーダ、中央盆地、クンブルガズといった、海の下を通る断層の特定のセグメントを調査し、1999年のイズミット地震や、より小規模な2019年のシリブリ地震といった、この地域における地震の歴史を調べました。彼らが開発している数学が、断層帯の実際の形状や構造に基づいている必要があるため、この地質学的背景は不可欠でした。次に、彼らは岩石の挙動を記述する新しい方法を導入しました。地球の地殻を、瞬時に跳ね返る単純なスプリングとしてではなく、どのように引き伸ばされたか、あるいは押し潰されたかという記憶を時間をかけて保持していくような材料として想像してみてください。研究者たちは、この記憶効果を記述するために、カプト（Caputo）微分と呼ばれる数学的ツールを使用しました。このツールにより、岩石の物理的特性と、波が伝播する際のエネルギー損失との関係を結びつけるモデルを作成することができました。

この新しいモデルを用いて、研究者らは3種類の異なる地震波——最も速く最初に到達するP波、地面を横方向に揺らすS波、そして地表を転がるように進むレイリー波——がどのように振る舞うべきかを導き出しました。彼らは、自分たちのモデルが常にエネルギーの損失をもたらすことを数学的に証明しました。これは、地震に関するあらゆる現実的な記述において必要な条件です。もしモデルが、波が伝わるにつれてエネルギーを得ることを示唆した場合、それは物理的に不可能です。彼らの方程式は、分数的なアプローチが、幅広い周波数にわたって現実の岩石の挙動と一致する、滑らかで連続的なエネルギー損失の曲線を生み出すことを示しました。これは、複雑になりすぎることを避けずに異なる周波数のデータに適合させることに苦慮することが多い、従来のモデルに対する重要な改善です。

この理論が実際に機能するかどうかをテストするため、チームはコンピュータ生成データを用いた制御実験を行いました。彼らは、波がどのように減衰すべきかという正確なルールを知っている合成地震シナリオを作成しました。そして、その波のパターンのみに基づいて、それらのルールを特定するように彼らの分数モデルに求めました。結果は驚くべきものでした。分数モデルは、非常に高い精度で正しいパラメータを復元することに成功し、整数を用いた伝統的なモデルと比較して、予測誤差を87パーセント近く減少させました。伝統的なモデルでは、研究者が数値を無理に適合させなければならず、結果として実際の波の挙動との一致は不十分なものとなりました。分数的な「中間」の記憶効果を許容することで、分数モデルは波の損失の物理をより自然に捉えることができたのです。

この研究は、この枠組みを実世界のデータにどのように適用できるかについても検討しました。研究者らは、2019年のシリブリ地震を最初の実世界テストの最適な候補として特定しました。彼らは、エンジニアや地震学者がこの事象からの実際の記録をダウンロードし、彼らの新しい方程式を使用して揺れの強度を再構成するための明確な計画を概説しました。これにより、建築家が耐震建築を設計するために用いる数値である、最大加速度や最大速度といった地面の運動の精密な測定値を算出することが可能になります。著者らは、このツールは既に起こったことを分析するためのものであり、未来を予測するためのものではないと慎重に述べています。これは次の大きな破壊がいつ起こるかを教えてくれるものではありませんが、一度発生した際の揺れの激しさを理解するための、より鋭いレンズを提供してくれるのです。

厳密な数学的証明と明確な地質図、そして実データ分析への計画を組み合わせることで、本論文は抽象的な理論から実践的な応用へと議論を進めています。これは、分数的なアプローチが単なる数学的な好奇心の対象ではなく、マルマラ海における地震波の複雑な現実を記述するための、実行可能でより正確な方法であることを示しています。この研究は、断層帯の層状で混沌とした地質学と、エンジニアリングの安全に必要な、クリーンで予測可能な方程式との間の溝を埋めるものです。数千件の実際の地震記録に対する完全な較正は将来の課題ですが、本研究はその基礎となる道筋を整え、その手法が理論およびシミュレーションにおいて有効であることを証明し、イスタンブールの足元にある実際の地面に対してテストされる準備が整っていることを示しました。