2026.04.24
【プレスリリース】地球型ダイポール磁場の南北2つの磁極の安定解を発見―地球磁場反転の物理機構解明への大きな手がかり―
概要
よく知られているように、地球は「大きな磁石」となっており(ダイポール磁場※1が存在しており)、その磁力は、地球内部の外核にある液体状の鉄が熱対流することにより起きるダイナモ過程※2によって生み出されていると考えられています。そして、この地球の磁場は、数十万~千万年程度の不規則な間隔で、その極性が反転していることが、これまでの地質学的な研究で明らかになっています。しかし、その物理機構は未だ解明されていません。特に、磁場の極性(北向き・南向き)がどのように決定されるのかについても、十分には理解されていない状況にあります。
自然科学研究機構核融合科学研究所/総合研究大学院大学の長谷川裕記助教、大谷寛明准教授、佐藤哲也名誉教授の研究チームは、この極性決定機構に着目し、地球の外核と同じ球殻状のプラズマで生じる対流ダイナモに対して、3次元電磁流体シミュレーション※3を用いた詳細な研究を行いました。その結果、地球型ダイナモでは、磁場の極性(北向き・南向き)は対流の向きではなく、初期に存在するごく微小な磁気擾乱によってランダムに決定されることを初めて明らかにしました。しかも、与える磁気擾乱の微妙な違いによって、その極性が北向きの状態か南向きの状態のいずれかの状態に落ち着き、その状態が続くことが明らかになりました(磁極の双安定性)。したがって、地球磁場についても、その極性は40億年前の磁場発生時に存在していた微小な揺らぎによって偶然に選択された可能性があります。そして、その地球磁場の極性は、そのまま維持されるはずですが、現実には度々反転が起きています。このことから、地球磁場の反転は、今回の計算モデルでは考慮していない物理効果によって引き起こされている可能性が示唆されます。
この研究成果をまとめた論文がScientific Reportsに2026年3月3日に掲載されました。
研究背景
うっそうとした森を歩くとき、進むべき方向を知るために欠かせないコンパスは、地球に磁場が存在することを利用しています。方位磁針が北を指し示すことは2世紀頃までには知られており、17世紀には「地球は一つの巨大な磁石である」と考えられるようになりました。地球磁場の起源は、アインシュタインも重大な未解決問題であると考えていましたが、この地球のダイポール磁場(図1左上)は地球内部の外核にある液体状の鉄の熱対流運動にともなう電磁流体相互作用(ダイナモ過程)により発生していることを、核融合科学研究所のシミュレーショングループ(注)とカリフォルニア大学ロサンゼルス校(UCLA)のグループが、ほぼ同時期に、電磁流体シミュレーションによって証明しました。一方、これまでの地質学的な研究(古地磁気学)により、地球磁場は、数十万~千万年に一度程度の不規則な頻度で、その方向(極性)が逆転(反転)すること (図1左下)、さらに、その極性の反転には千年のオーダーの時間を要することが示されています。このようなダイポール磁場の生成や周期的反転に関しては、電磁流体シミュレーションの研究は数多くなされています。そして、地球型磁場の極性が不規則(非周期的)に何度も反転する現象については、核融合科学研究所のグループが、四半世紀前にその再現に成功しています[Li, Sato & Kageyama, Science (2002)](図1右)。ただ、残念なことに、磁極がどのような原因で不規則に反転を繰り返すのか、そのトリガー機構は全く不明で、前人未到の難題となっています。そこで、本研究では、反転現象の解明に近づくため、磁場極性がどのように決定されるのかという点に着目し、地球のような複雑さを取り除いた、単純な対流が存在する状況について、電磁流体シミュレーションを行いました。
図1:地球磁場の模式図(左上)、海洋底の岩石の残存磁場から得られた地球磁場反転の履歴(Ogg, Geologic Time Scale 2020, Elsevier (2020)のデータより作成)(左下)と、地球磁場反転シミュレーションの結果(右)。海嶺でマントル内部から湧き出した岩石は、冷え固まるときに地球磁場によって帯磁します。海嶺で生成された岩石は海嶺の両側へ移動していくため、岩石に残された磁場の極性を調べると、反転の履歴が縞模様のようになって現れます。シミュレーション結果は、四半世紀前の核融合科学研究所のシミュレーションで再現された不規則な反転現象の様子を示しています。
研究成果
核融合科学研究所/総合研究大学院大学の長谷川裕記助教、大谷寛明准教授、佐藤哲也名誉教授の研究チームは、地球の外核と同じ球殻状のプラズマで生じる対流ダイナモにおける磁場極性の決定機構を調べるため、インヤン格子※4を適用した3次元電磁流体シミュレーションコードを用いた詳細な研究を、核融合科学研究所のスーパーコンピュータ“プラズマシミュレータ「雷神」・「双星」”※5を使用して行いました。シミュレーションでは、図2に示したような、ほぼ定常なねじれた対流を与え、初期の磁場として50種類のランダムかつ微小な磁気擾乱を用いました。その結果、すべての計算で、ダイポール磁場が支配的な状態に到達し、その極性は、北向きと南向きの2つの状態がほぼ同確率で出現することが明らかになりました。対流の方向を反転させてもこの確率は変化せず、極性は背景の対流ではなく初期の微小な磁気擾乱によって決定されることが示されました。また、磁場の成長過程が、磁気拡散時間※6のタイムスケールで周期的に極性反転が繰り返される初期の時間帯と、それ以降の極性が定まり磁場が安定な状態に至る時間帯の2つに分けられることを見出しました。さらに、安定な状態に達したあと、再度、微小な磁気擾乱を加える計算を行ったところ、形成されたダイポール磁場に顕著な変化は見られませんでした。このことは、このダイポール磁場の安定状態は非常に強固で、小さな磁気擾乱では極性の反転は起きないことを示しています(図3)。
本研究は、球殻ダイナモにおいてダイポール磁場が2つの安定な平衡状態として存在することを示すとともに、地球磁場の極性が初期の微小揺らぎによって偶然に選択され、その後は強固に維持される可能性があることを示唆するものです。このことから、磁場の反転を引き起こすためには、通常の電磁流体理論の範疇には入っていないトリガー(例えば、ミクロなプラズマ不安定性などによる異常な磁気拡散率※6や粘性)によってこの安定状態が破られる必要があると、研究チームでは考えています。
研究成果の意義と今後の展開
これまでのシミュレーション研究では、磁場の反転は主に磁気拡散時間のタイムスケールで起こる周期的な現象として理解されてきました。一方で、実際の地球磁場の反転は、それよりもはるかに長い時間スケールで不規則に起こることが知られており、その物理機構は未解明のままです。核融合科学研究所のグループは、2000年代初頭に、このような長時間スケールの不規則な磁場反転の再現に成功しています(図1右)。本研究で明らかになった磁場の双安定性は、このような反転現象と一見矛盾するように見えますが、研究チームでは、この関係にこそ反転のトリガー機構を理解するための重要な手がかりがあると考えています。
今後は、この安定状態を破る要因を明らかにすることで、地球磁場反転の物理機構の解明につながることが期待されます。
電磁流体シミュレーションには、避けることができない差分(数値)誤差が存在します。流体のシミュレーションでは、連続的な空間を多くの格子点で分割して表現します(離散化と呼びます)が、この差分化(離散化)は局所的に数値的な磁場の散逸をともないます。研究チームでは、過去の不規則な反転現象のシミュレーション結果は、この差分誤差がもたらした結果ではないかと考えています。流体のシミュレーションにおいて、エネルギーなどが集中する場所(よどみ点)ではこの差分誤差が増大します。一方、現実のプラズマは連続体ではなく粒子から構成されており、“よどみ点”のような領域では、粒子集団が平衡状態からずれることにより、プラズマ不安定性が発生する可能性があります。研究チームでは、このようなミクロな効果が磁場反転のトリガーとなる可能性があると考えています。今後は、まず、高精度な電磁流体シミュレーションにより、従来モデルの範囲でどこまで再現できるかを精査するとともに、その限界を明らかにした上で、ミクロな物理過程を取り入れた新たなモデルの構築へと発展させていく予定です。
本研究は、日本シミュレーション学会と核融合科学研究所によって進められている「遠非平衡科学の確立に向けたシミュレーション研究プロジェクト~マクロ・ミクロ連結階層シミュレーションによる『未来を見る望遠鏡』構想~」の一部として実施しています。本プロジェクトは、災害をもたらす地震などの自然現象の予測や、それらの現象による被害の最小化なども目標のひとつとしています。一旦、地球磁場の反転が起き始めると、完全に反転が終了するまでの千年程度の間、地球磁場は非常に弱くなり、宇宙から飛来する高エネルギー宇宙線を防ぐことができなくなり、生態系や人類に多大な影響をおよぼす可能性がありますが、地球磁場反転機構の解明は、そのような悪影響を予測、最小化するためにも大変重要となります。
図2:シミュレーションで用いた対流の模式図とダイポール磁場の成長過程。上段の模式図において、水色の領域にプラズマがあり、白色矢印の箇所でプラズマに力を与え続けると、黒色矢印のような対流が生まれます。中段の図は、北向き極性と南向き極性、それぞれの場合のダイポール磁場の成長過程を示しています。下段の図では、初期の周期的な小振幅の反転が起きる時間帯(薄黄色の領域)、その後の極性が定まる安定な状態の時間帯(薄緑色の領域)を示しています。さらに、安定状態において微小な磁気擾乱を加えても、顕著な変化は見られないことがわかります。
図3:プラズマの流れのエネルギーと磁場の極性の関係を示した概念図。両者の関係は、W型のポテンシャル形状になっていると考えられ、どちらの谷に落ちるかは、初期の擾乱によって決まります。また、谷が深いために、少しの擾乱では、隣の谷に移動できないと考えられます。
(注)中日新聞(1996年9月8日朝刊)の第1面に掲載。
用語解説
論文情報
研究サポート
本研究は、大学共同利用機関法人自然科学研究機構「ネットワーク型研究加速事業」(01422301)の助成を受けたものです。また、本研究の数値計算は、核融合科学研究所の一般共同研究(プラズマシミュレータ共同研究)(NIFS24KISM005)の支援を受けて、プラズマシミュレータ「雷神」(NEC SX-Aurora TSUBASA)、並びに、同「双星」(NEC LXシリーズ)において行いました。
本件のお問い合わせ先
本件の広報及び取材の申し込みについて
管理部 総務企画課 対外協力係
電話: 0572-58-2016
E-mail: taigai-kakari(at)nifs.ac.jp
総合企画課 広報社会連携係
電話:046-858-1629
E-mail:kouhou1(at)ml.soken.ac.jp
(at)を@に変更して送信してください。