Commutative Coarse-graining on the Sphere

技術要約：球面における可換な粗視化

問題提起
本論文は、球面上（$S^2$）のベクトル場における粗視化およびフィルタリングにおける根本的な不整合性に対処するものである。これは、ラージ・エディ・シミュレーション（LES）、マルチスケール動力学解析、ならびに地球物理学や天文学におけるデータ処理において極めて重要な操作である。標準的な球面上のフィルタリング手法は、通常、空間畳み込みに依存している。しかし、粗視化された方程式を導出するための基本的な整合性の要件として、フィルタリング演算子は空間微分演算子（勾配、発散、および回転）と可換でなければならない。

著者らは、一般的なフィルタリング手法がこの可換性を満たしていないことを示している。この失敗は、数値グリッドの離散化、座標系の選択、あるいは計量項によるアーティファクトではなく、球面の曲率に固有の性質である。具体的には以下の通りである：

1. 成分ごとの畳み込み（Componentwise Convolution）： 球面ベクトル成分（径方向、帯状、子午線方向）をスカラーとして扱い、それらを独立してフィルタリングすることは、幾何学的構造を無視することになる。これは、一定の速度を持つ純粋な帯状流がスケール不変であることを誤って示唆しており、球面の固有の勾配や渦度を無視している。
2. 周囲空間での畳み込み（Ambient Convolution）： 3次元ユークリッド空間（$\mathbb{R}^3$）の埋め込み空間内で、デカルト成分を平均化することで畳み込みを行う手法は、基底ベクトルの変化を考慮してはいるものの、一般的な流れに対して場の接線方向または径方向の性質を保持することには失敗する。これにより、「径方向と横方向の混合（radial-lateral mixing）」が生じ、フィルタリングされた場が本来の物理的特性（例：本来は非発散であるべきものが発散を持つようになる等）を失う原因となる。

手法および理論的枠組み
著者らは、先行研究[1]で提案された一般化された球面畳み込みの枠組みに基づいている。彼らの手法の核となるのは以下の通りである：

• 対称性解析： ユークリッド空間の畳み込みは並進等変性により微分と可換であるが、球面の畳み込みは回転等変性（$SO(3)$）に依存しなければならないことを認識している。
• スペクトル分解： エドモンズ（Edmonds）のベクトル球面調和関数基底を利用する。著者らは、標準的な「周囲空間（ambient）」の畳み込みが、この基底に対して、特定の成分において次数がシフトした固有値を持つ対角行列として作用することを示す。このシフトが、径方向、ポテンシャル、およびソレノイダル（渦なし）モードの混合を引き起こし、一般的な流れに対する微分演算子との可換性を崩壊させる。
• 一般化された畳み込みの構築： 可換性を回復するために、著者らは、エドモンズ基底におけるベクトル場の特定の成分に対して、スペクトル的にシフトされたカーネルを適用する一意な線形演算子である一般化された畳み込み（$G \overrightarrow{\otimes}$）を定義する。
  • この演算子は、すべてのベクトル球面調和モード（$Y_{n,j}, \Phi_{n,j}, \Psi_{n,j}$）に対して、固有値 $\dot{b}_G(n)$ を保持するように定義される。
  • これは、レジャンドル係数 $\dot{b}_G(n-1)$ および $\dot{b}_G(n+1)$ を持つシフトされたカーネル $G^{(1)}$ および $G^{(2)}$ を定義することで達成され、幾何学によって導入される次数のシフトを実質的に「打ち消す」ものである。
• 実用的な実装： 本論文では、この一般化された畳み込みが、ベクトル場をヘルムホルツのスカラー（径方向、ポテンシャル、およびソレノイダル）に分解し、それぞれにスカラー畳み込みを適用した後に再構成することによって、完全な球面調和変換を行わずに効率的に計算できることを示している。

主な貢献

1. 根本原因の特定： 本論文は、可換性の失敗が、座標の特異点やグリッド解像度に関連する問題ではなく、多様体 $S^2$ の幾何学的特性であることを厳密に確立している。
2. 演算子の唯一性： 著者らは、与えられたスカラーカーネルに対して、一般化された畳み込みが、球面上での径方向射影、発散、および回転と可換となる一意な線形演算子であることを主張している。これは、あらゆる一貫した線形粗視化手法（スペクトル切断、ウェーブレット変換、拡散、または機械学習による畳み込みを含む）が、この形式をとらなければならないことを意味する。
3. 誤差の定量化： 非可換フィルタによって導入される誤差の大きさを定量化し、それらが物理信号自体を上回る可能性があることを示している。

数値結果
著者らは、球面上における単層浅水方程式（SWE）を用いた不安定な帯状ジェットのシミュレーションを用いて、理論の検証を行っている。

• 発散テスト： 厳密に3次元発散フリー（$\nabla_3 \cdot \mathbf{v} = 0$）である流れに対し、一般化された畳み込みは、数値丸め誤差の範囲内（$O(10^{-22} - 10^{-20})$ s$^{-1}$）でこの性質を保持する。対照的に、周囲空間の畳み込みは、$O(10^{-2})$ s$^{-1}$ の偽の発散を生じさせ、これは実際の流れの横方向の発散と同等の大きさであり、18桁も大きい誤差である。
• 勾配テスト： 高さの勾配（$\nabla h$）をフィルタリングする場合、一般化された畳み込みは、相対誤差 $O(10^{-12})$ で球への接線性を維持する。周囲空間の畳み込みは、勾配の大きさの約20%に達する径方向成分を導入してしまう。成分ごとの畳み込みは、信号と同等の誤差（$O(1)$）をもたらす。
• 渦度テスト： ソレノイダル（発散フリー）な流れに対しては、周囲空間の畳み込みは回転（渦度）と可換である。これは、ソレノイダル基底関数（$\Psi_{n,j}$）が周囲空間の畳み込み演算子の固有関数であるためである。したがって、周囲空間の手法を用いて計算されたフィルタリング後の速度の渦度は、フィルタリングされたスカラー渦度と丸め誤差の範囲内で一致する。しかし、この可換性はソレノイダル場に限定された特殊なケースであり、一般的な（発散成分を持つ）流れには拡張されない。一般的な流れでは、周囲空間の畳み込みは接線性を保持できず、勾配との可換性にも失敗する。

意義と主張
本論文は、一般化された畳み込みが、球面上の一貫した粗視化方程式を導出するために必要な数学的枠組みを提供すると主張している。この枠組みがなければ、標準的なフィルタリングは、偽のソース（源）やシンク（吸収）を生み出し（発散誤差）、粗視化された動力学の物理を汚染する勾配誤差を導入することになる。

著者らは、この構成が以下のために不可欠であることを強調している：

• 全球的な流れに対する物理空間でのLES方程式の導出。
• 球面上のデータ同化および逆問題が、基礎となる物理的対称性を保持することを保証すること。
• 回転等変なニューラルネットワークなどの機械学習モデルを、ベクトル場やテンソル場へ正しく拡張し、幾何学的な不整合によって訓練範囲外で失敗しないようにすること。

本研究は、数十年にわたる地球物理学的モデリングの慣習に対する修正であり、惑星物理学、宇宙論、および慣性閉じ込め核融合など、球面幾何学が本質的な分野における将来の応用への指針となるものである。