衛星の重力データで「重力の穴」を可視化してみた ― インド洋低ジオイド域、発見から75年越しの謎解きまで

海の上に、周囲より海面が最大100メートル前後も低い「窪み」が存在するのをご存知でしょうか。インド洋、スリランカ沖合を中心とした約300万平方キロメートルの海域。専門的には「インド洋低ジオイド域(Indian Ocean Geoid Low: IOGL)」、科学者の間では通称「重力の穴(Gravity Hole)」と呼ばれています。

今回はこの「重力の穴」について、NotebookLMのDeep Research機能でネット上の一次情報を横断調査し、発見の経緯と現在の科学的コンセンサスをまとめました。あわせて、GOCE・GRACEなど実際の重力観測衛星のデータを統合した「ジオイドモデル」を使い、この凹みを自分の手で可視化してみましたので、その方法と結果をご紹介します。

「重力の穴」とは何か

波や潮汐の影響を除いた平均海面（ジオイド）を計算すると、インド亜大陸の南、スリランカ沖合の海水は、地球全体の平均よりも最大で100メートル前後も地球の中心に近い側へ沈み込んでいます。原因は目に見える地形の凹みではなく、この領域の下にある地球内部（マントル）の密度が周囲より低く、重力が弱いことです。重力が弱い方向には海水が引き寄せられにくいため、結果としてジオイド面（重力が等しい面）がへこんで見えるわけです。

ただし勾配は非常になだらかで、水平方向に数百キロ単位でようやく1メートル沈む程度です。実際に船で通りかかっても、肉眼で坂道のような落ち込みが見えるわけではありません。オーストラリアの巨岩ウルル（エアーズロック）になぞらえて「地球のへそ」と呼ばれることもあります。

なぜ普通の衛星画像では見えないのか

「衛星画像を見れば凹みが写っているのでは」と思われるかもしれませんが、気象衛星の画像ではこの高低差を確認できません。例えば欧州の静止気象衛星Meteosatは雲・水蒸気・地表温度の観測に特化しており、重力やミリメートル単位の海面高度を測るセンサーは搭載していないため、可視光カメラ画像にこの凹みが写り込むことはありません。

「重力の穴」を可視化するには、重力そのものや海面高度を専用に測定する衛星データが必要です。具体的には、ESA（欧州宇宙機関）の重力観測衛星GOCEやNASAのGRACE、海洋高度計を搭載したJasonシリーズ・Sentinel-3・SWOTなどが挙げられます。

発見から解明まで、75年越しの謎解き

NotebookLMのDeep Researchで一次情報（Wikipedia、Forbes Africa、The Economic Times、European Geosciences Unionの記事等、54件のソースを調査・35件をインポート）にあたったところ、「重力の穴」の存在が知られてから原因が解明されるまで、実に75年もの歳月がかかっていたことが分かりました。

年出来事1923〜1938年オランダの地球物理学者フェリックス・ベニング・マイネスが、自作の三振り子重力計を潜水艦に搭載し世界各地で洋上重力測定を実施（総航程10万海里超）1948年マイネスが観測データを “Gravity Expeditions at Sea 1923–1938” (Vol.ジオイドの詳細な形状がより鮮明になる1979年F. J.GRACEと組み合わせ、より精細な重力場モデル（GOCO05c等）が作成される2023年5月インド科学大学院（IISc）のDebanjan Pal氏・Attreyee Ghosh氏らがGeophysical Research Lettersに発表。

「衛星で存在は分かっていたが、なぜそこにあるのかは長年謎のまま」という状態が2023年まで続いていたことになります。

2023年、ついに原因が特定される

2023年の研究（Pal & Ghosh, Geophysical Research Letters）は、過去1億4000万年分のプレート運動を再現する3次元マントル対流シミュレーションを行い、この重力異常の成因を突き止めました。要点は次の通りです。

• 約1億2000万年前、超大陸ゴンドワナからインドプレートが分離してユーラシア大陸へ北上する過程で、太古の「テチス海」の海洋プレート（スラブ）が次々とマントル深部（核とマントルの境界）へ沈み込んだ
• 沈み込んだ冷たく重いスラブが、東アフリカ下の深部にある巨大な低速度領域「アフリカン・ブロブ」（LLSVP: 超巨大マントル低速度層）を押しのけ、高温で密度の低いマントル物質（プルーム）を側方・上方へ押し出した
• このプルームがインド洋の下へ広がり、周囲より密度の低い（＝重力の弱い）領域を形成。地震波トモグラフィーでもこの領域の地震波速度が遅い（高温・低密度を示唆）ことが確認されている
• シミュレーション上、現在の「重力の穴」の形状は約2000万年前ごろに形成されたと推定される

研究者のAttreyee Ghosh氏は「地球上で最も顕著な重力異常の一つ」と述べており、数千万年前のプレート運動の痕跡が、現在も地球の重力場という形で残り続けていることになります。ただし、トモグラフィーモデルを密度へ変換する際に化学的な不均質性と温度差のどちらがどれだけ寄与しているかなど、細部は今も議論が続いている活発な研究領域です。

衛星の重力データから自分で可視化してみた

ここからは実践編です。GOCE・GRACEの生の観測データ（Level-1/2）から独自にジオイドを計算するには球面調和関数への大規模な最小二乗解析が必要で、研究機関レベルの処理になってしまいます。そこで今回は、GOCE・GRACE等19機の衛星による約11億6000万点の観測データを国際的な研究チーム（GOCOプロジェクト：グラーツ工科大学、ボン大学、ミュンヘン工科大学ほか）が既に解析・統合した重力場モデルGOCO06sを、配布元のICGEM（International Centre for Global Earth Models、GFZポツダム地球科学研究センター運営）から取得しました。認証不要・CC BY 4.0で誰でも使えます。

処理の流れは次の通りです。

1. ICGEMからGOCO06sの重力場係数ファイル（.gfc形式）をダウンロード
2. geoid_toolkit（Barthelmes (2009) の反復法による球面調和関数展開の実装）を使い、指定した緯度経度グリッド上でジオイド高（基準はWGS84楕円体からの偏差、単位はメートル）を計算。展開の次数はlmax=150
3. インド洋領域（緯度-60°〜30°、経度0°〜130°、解像度0.2°）と全球（解像度1.0°）の2種類のグリッドについて、海岸線・カラーバー・最深点ラベル付きのプレビュー画像とGeoTIFFを出力

Jason/Sentinel-3/SWOTなどの海面高度計データは海域の検証・補完に使われますが、配布元での利用者登録が必要なため、今回は認証不要でフルカバレッジのGOCE/GRACE系モデルのみで完結させています。

可視化結果

【図1】インド洋低ジオイド域（重力の穴）のジオイド高分布

インド亜大陸・スリランカのすぐ南から南インド洋にかけて、大きく青く窪んだ領域が一目で確認できます。計算上の最深点は北緯4.8°、東経79.2°（スリランカ南方沖）で、ジオイド高-106.1mでした。

【図2】全球ジオイドモデルにおける「重力の穴」の位置関係

視野を全球に広げると、インド洋の凹みが地球上で最も深い負の異常であることが分かります。対照的に北大西洋やニューギニア付近は正の異常（重力がやや強い＝ジオイドが盛り上がっている）になっており、地球の重力場が決して均一な球ではないことが視覚的に理解できます。

まとめ：数字で答え合わせをしてみる

今回自分で計算した最深点の値は-106.1m（北緯4.8°、東経79.2°）でした。これを先行研究と突き合わせてみると、

• 1979年のGEM10モデル：約-105m
• GOCE/GRACE統合モデルGOCO05cベースの先行研究：約-106m（3万km²、348フィート）
• 今回の計算（GOCO06s、lmax=150）：-106.1m（北緯4.8°、東経79.2°）

と、45年前の初期モデルから現在までほぼ同じ値に収束していることが確認できました。位置についても、学術文献で報告されている中心座標（おおよそ北緯3〜6°、東経77〜80°付近）とよく一致しています。手元のPCで、認証不要の公開データだけを使って、45年分の測地学の歴史が積み上げてきた数字を再現できたのは、素直に面白い体験でした。

一方で、「なぜそこが低いのか」という原因の解明には、衛星で存在が分かってから2023年まで、さらに数十年を要しています。数値としての「穴」を見つけることと、その成因を地球内部のダイナミクスから説明することの間には、まだ大きな距離があるのだと実感させられました。テチス海の消滅というはるか昔の出来事が、今なお海面の高さという形で観測できるというのは、地球のスケール感を考えさせられる話です。