Spectral interpolation in semi-implicit semi-Lagrangian methods for shallow water equations on the sphere

技術的概要：球面上の浅い水方程式に対する半陰的半ラグランジュ法におけるスペクトル補間

問題提起
半陰的半ラグランジュ（SISL）法は、クーラン・フリードリヒス・レウィ（CFL）条件で許容されるものよりも大きな時間刻みを可能にするため、運用上の全球気象・気候モデルの標準となっています。しかし、現在の運用スキームには重大な数値的不整合が存在します。予報変数の微分の空間離散化には、通常、スペクトル的に正確な手法（球面調和関数変換 [SHT] やダブルフーリエ球 [DFS] 法など）が採用されていますが、オフグリッドの離脱点における変数の評価に必要な補間ステップは、3 次または 5 次テンソル積ラグランジュ多項式などの低次手法に依存しています。

この不一致は、数値精度の相違を生み出します。低次補間は、計算コストの懸念（スペクトル補間に対する $O(n^4)$ に対して $O(n^2)$ の複雑さ）および単純な移流に対する理論的な誤差境界によってしばしば正当化されますが、浅い水方程式（SWE）のような複雑なシステムにおけるこのトレードオフを正当化する数値的証拠は文献において限られています。さらに、低次補間は過度の数値拡散や位相誤差を導入することが知られており、長期間の積分において精度や保存特性を劣化させる可能性があります。

手法
本研究は、球面上の SWE に対する SISL 法におけるスペクトル的に正確な補間の最初の数値的調査を提示します。著者らは、安定なフーリエ基底と 2 時間レベルの SETTLS スキームを利用する、吉村 [56] によって開発された DFS ベースの SISL モデルを基盤としています。

手法の中核的な革新は、空間微分に使用されるのと同じ DFS フーリエ基底に基づくスペクトル補間スキームへの、標準的なラグランジュ補間の置換です。計算コストの障壁に対処するため、著者らは非一様高速フーリエ変換（NUFFT）を採用しています。

• 補間メカニズム: 予報変数（特に半陰的システムにおける項 $R_v$ と $R_h$）は、部分正則性 DFS 基底を用いて展開されます。これらの展開は、非グリッドの離脱点において評価されます。
• 複雑性の低減: $O(n^2)$ の離脱点における二変量フーリエ級数の直接評価には $O(n^4)$ の演算が必要となります。級数を複素指数形式に再構成し、NUFFT（特に FINUFFT ライブラリ）を利用することで、計算複雑性は $O(n^2 \log^2(1/\epsilon) + n^2 \log n)$ に低減されます（ここで $\epsilon$ は精度許容誤差です）。これにより、補間ステップは残りのスペクトル SISL スキームと計算上一貫性を持つようになります。
• 実装: 著者らは、このスペクトル補間を 2 つの文脈で実装しました。
  1. DFS ベースの SISL モデル（吉村の安定なフーリエ基底を使用）。
  2. SHT ベースの SISL モデル（ガウス・ルジャンドル格子を使用）。ここでは、非等間隔の緯度点から二変量フーリエ補間関数を生成するために DFS 法が適応されました。

主要な貢献と結果
著者らは、標準的な SWE テストケース（ウィリアムソン TC1、TC2、TC5、TC6、およびガレフスキー・テストケース）のスイートを用いて、新しいスペクトル補間スキームを標準的な低次ラグランジュ補間（LAG）と比較評価しました。結果は、いくつかの明確な利点を示しています。

1. 精度と収束性:

   • TC1（コサインベルの移流）: スペクトル DFS 補間は、LAG と同じ空間収束率を維持しつつ、誤差を約 2 桁低く達成しました。長期間の積分（10 回転）において、DFS 法は LAG に比べて誤差を 3 桁削減しました。
   • TC2（地衡流）: LAG 補間は、おそらく支配的な空間誤差に起因して、期待される 2 次時間収束を示せませんでしたが、DFS 補間スキームは理論的な 2 次精度を達成しました。
   • TC5 および TC6（山岳流およびロスビー・ホーウィッツ波）: DFS 法は、特に高分解能および小時間刻みにおいて、高さ場において LAG よりも一貫して低い誤差を生み出しました。TC6 において、DFS 法は波の谷の形状をよりよく保持し、高次オイラー法による参照解により密接に一致しました。

2. 保存特性:

   • スペクトル補間は、質量とエネルギーの保存を著しく改善しました。TC1 において、120 日後の絶対質量誤差は、DFS で $1.21 \times 10^{-11}$ に対し、LAG では $4.34 \times 10^{-4}$ でした。
   • TC2 および TC5 において、DFS 法は時間経過とともに安定した質量およびエネルギー保存誤差を維持しましたが、LAG の誤差は増加するか、より顕著に振動する傾向がありました。

3. 数値拡散と小規模ダイナミクス:

   • ガレフスキー・テストケース: このテストは、非線形ジェットが乱流へと進化することを扱います。運動エネルギースペクトルのスペクトル分析により、LAG 補間が高波数エネルギーの急速な散逸を引き起こしたことが明らかになりました。対照的に、DFS 補間スキームはエネルギーカスケードを保持し、高波数スペクトルが乱流回転流で観測される理論的な $k^{-5/3}$ スケーリングに近づくことを可能にしました。LAG 補間は、この挙動を捉えることに失敗し、過度の数値減衰を示しました。

4. 計算効率:

   • ウォールクロック時間の比較により、低精度要件の場合、LAG 補間の方が高速であることが示されました。しかし、必要な精度が高まる（厳格な誤差許容誤差）につれて、DFS 補間法はより効率的になり、LAG 手法よりも短い時間で目標精度に到達します。これは、NUFFT 加速されたスペクトル補間が高精度応用に viable であることを確認しています。

重要性
本論文は、スペクトル的に正確な補間を SISL 法に組み込むことが、空間離散化と移流補間の間の数値的不整合を解決すると主張しています。この研究は、このアプローチが単なる理論的な改善ではなく、具体的な利点を提供することを示しています。

• 空間微分計算と補間ステップの間の「精度ギャップ」を排除します。
• 数値拡散を低減し、それによって長期間の積分にわたって微細な構造や波のような特徴を保持します。
• 質量とエネルギーの保存を強化し、長期的な気候シミュレーションにとって重要です。
• NUFFT 加速のおかげで、過大な計算コストを負うことなくこれらの改善を達成します。

著者らは、低次補間が実用的な妥協点であった一方で、新しいスペクトルアプローチは、大気ダイナミクスの高精度・長期的シミュレーションに対する優れた代替案を提供し、より複雑な 3 次元ダイナミカルコアへの拡張を動機づける、と結論付けています。