P波・S波の地球内部での伝播・反射・屈折をリアルタイムアニメーション。地球の層構造（地殻・マントル・外核・内核）とシャドーゾーンを可視化。
地震波が異なる速度の層の境界で屈折する様子は、光の屈折と同様にスネルの法則で記述されます。これが波の経路を決める基本です。
$$ \frac{\sin i}{v_1}= \frac{\sin r}{v_2}$$ここで、$i$は入射角、$r$は屈折角、$v_1$, $v_2$はそれぞれ境界の上下の層における地震波の速度です。速度が大きい層（$v_2 > v_1$）へ進むと、波は境界から遠ざかる方向（屈折角が大きくなる）に曲がります。この法則に従って、シミュレーターはP波の経路を計算しています。
P波とS波の速度は、媒質の弾性定数と密度によって決まります。固体中での速度は以下の式で与えられます。
$$ v_P = \sqrt{\frac{K + \frac{4}{3}G}{\rho}}, \quad v_S = \sqrt{\frac{G}{\rho}} $$$v_P$はP波速度、$v_S$はS波速度、$K$は体積弾性率、$G$は剛性率（せん断弾性率）、$\rho$は密度です。液体では剛性率$G=0$なので、$v_S=0$となりS波は伝わりません。これが、シミュレーターでS波が外核で消える物理的な理由です。
震源決定と地震早期警報：複数の観測点でP波とS波の到着時間差を測定し、シミュレーターと逆の計算を行うことで、震源の位置（緯度、経度、深さ）を特定します。また、最初に到達するP波の情報から、後から来る強いS波や表面波の到達前に警報を発するシステムの基礎理論となっています。
地球内部構造の探査（地震学）：世界中で観測された地震波の伝播時間と経路を、このシミュレーターのようなモデルと比較することで、地殻、マントル、外核、内核の厚さや密度、弾性特性を明らかにします。シャドーゾーンの存在は外核が液体であることを示す決定的な証拠でした。
地下資源探査（反射法地震探査）：人工的に発生させた地震波（主にP波）の地下層からの反射波を記録・解析します。シミュレーターで見た波の反射・屈折の原理を利用して、石油・天然ガスを含む地層や断層の位置・形状を画像化します。
建築物・地盤の耐震設計：特定の地点に到達する地震波の種類（P波、S波、表面波）やその特性（周期、振幅）を理解するために、波動伝播の基礎知識が不可欠です。地盤がS波をどのように増幅するかを評価し、建築物の設計に反映させます。
このシミュレーターを使い始めるときに、いくつか気をつけてほしいポイントがあるよ。まず、「波の速度は層ごとに固定」という前提を理解しておこう。実際の地球では、同じマントル内でも深さとともに圧力・温度・密度が変化し、速度は連続的に増加しているんだ。シミュレーターではそれを「層」で簡略化しているから、境界での急激な屈折が強調されて見える。実データではもっと滑らかな経路になることもあるんだ。
次に、「震源距離」と「観測点までの距離」を混同しないこと。シミュレーターの「震源距離」は、震源から見た地表上の点までの角度距離（例えば103度）だ。でも、実際に震源の深さが深いと、波が地表に出てくるまでの経路が変わる。例えば深さ600kmの地震だと、シャドーゾーンの範囲がシミュレーター（震源地表仮定）と変わってくる。常に「震源の深さ」を意識するクセをつけよう。
最後に、振幅（マグニチュード）の変化を過大解釈しないこと。スライダーで振幅を変えられるけど、これは単純化された表示だ。実際の地震波の振幅は、伝わる距離だけでなく、地盤の増幅効果や散乱、減衰によって大きく変わる。シミュレーターでS波の振幅を大きくしても、液体の外核では絶対に伝わらないという根本原理は変わらない。この「伝わるか伝わらないか」の質的な違いに注目するのが学習のコツだね。
このシミュレーターの背後にある波動伝播の考え方は、CAEエンジニアなら絶対に知っておくべき、様々な分野の基礎になってるんだ。まず真っ先に挙がるのは「非破壊検査」だ。例えば、構造物の内部きずを探る超音波探傷試験では、P波に似た縦波とS波に似た横波を使い分けて、欠陥の種類や位置を特定する。金属中の伝播も、異なる材質の境界で屈折・反射する様子は地震波とまったく同じ物理だ。
もう一つは「地盤工学」や「土木構造物の耐震設計」。ここでは、シミュレーターで見た「地盤の層構造」がそのまま応用される。地表付近の軟弱な層（地殻に相当）を伝わる波が、硬い岩盤（マントルに相当）との境界でどう屈折し、地表でどう増幅されるかを計算する。これが建物や橋の固有周期と共振しないかを評価する地盤応答解析の出発点になる。
さらに視野を広げると、「音響工学」や「水中ソナー」にも通じる。水中ではS波は伝わらない（液体だから）のでP波（音波）だけが伝播する。潜水艦のソナーが海中の物体を探知するとき、水温や塩分の違いによる音速の層構造（深海チャネル）を考慮する必要があるが、これもスネルの法則で記述される。地球物理学で培われた波動伝播のノウハウは、まさにマルチフィジックスシミュレーションの共通言語なんだ。
このシミュレーターに慣れてきたら、次は「逆問題」の考え方に挑戦してみよう。シミュレーターは地球の構造（速度モデル）が既知として波の経路を計算する（順問題）。しかし実務では、観測点で得られたP波、S波の到着時間データから、震源位置や地下構造を逆算する。この逆問題を解くための基本的な手法が「走時曲線」の解析だ。震源距離に対する到着時間をプロットした曲線の形から、速度構造を推定するんだ。シミュレーターでいろんなパラメータを変えながら、自分で簡易的な走時曲線をスケッチしてみると、理解がぐっと深まるよ。
数学的な背景を少し深掘りしたいなら、波動を記述する偏微分方程式に触れてみよう。地震波の伝播は、弾性体の運動方程式から導かれる波動方程式で表現される。一次元の単純なケースでは次の形だ： $$ \frac{\partial^2 u}{\partial t^2} = v^2 \frac{\partial^2 u}{\partial x^2} $$ ここで $u$ は変位、$v$ は波速だ。この式の解が進行波の形をとること、そして波速 $v$ が先ほどの $v_P$ や $v_S$ の式とどう結びつくかを調べるのが次のステップ。これが理解できれば、FEM（有限要素法）による地震波シミュレーションの論文も読み始められるはずだ。
最後に、具体的な次のトピックとして「表面波（レイリー波・ラブ波）」を薦める。シミュレーターで扱っているのは地球の内部を伝わる「実体波」（P波、S波）だけど、実際の地震で大きな揺れと長い周期の振動をもたらすのは、地表付近を伝わる表面波だ。この波の分散性（周期によって速度が変わる性質）を分析することで、地殻の厚さや浅部の構造を詳細に推定できる。内部構造から地表の現象へと視点を広げることで、地震工学への応用がよりリアルに感じられるようになるよ。