[論文レビュー] Initialization and the quasi-geostrophic slow manifold

Stephen M. Cox, A. J. Roberts|ArXiv.org|Mar 7, 2003

Quantum chaos and dynamical systems参考文献 31被引用数 4

ひとこと要約

本論文は、LorenzのL86モデルにおける準地転流的慢い多様体の分析を通じて、大気モデルにおける初期化を調査し、最も遅い不変多様体（SIM）と真の慢い多様体が指数的に小さい項を除いて同一であることを示した。重力波の分散が避けられず、初期重力波活動の二乗に反比例することを示し、正規形が波の分離と初期化手順における本質的予測分散を明らかにした。

ABSTRACT

Atmospheric dynamics span a range of time-scales. The projection of measured data to a slow manifold, ${\cal M}$, removes fast gravity waves from the initial state for numerical simulations of the atmosphere. We explore further the slow manifold for a simple atmospheric model introduced by Lorenz and anticipate that our results will relevant to the vastly more detailed dynamics of atmospheres and oceans. Within the dynamics of the Lorenz model, we make clear the relation between a slow manifold $\cal M$ and the ``slowest invariant manifold'' ({\SIM}), which was constructed by Lorenz in order to avoid the divergence of approximation schemes for $\cal M$. These manifolds are shown to be identical to within exponentially small terms, and so the {\SIM} in fact shares the asymptotic nature of $\cal M$. We also investigate the issue of balancing initial data in order to remove gravity waves. This is a question of how to compute an ``initialized'' point on $\cal M$ whose subsequent evolution matches that from the measured initial data that in general lie off $\cal M$. We propose a choice based on the intuitive idea that the initialization procedure should not significantly alter the forecast. Numerical results demonstrate the utility of our initialization scheme. The normal form for Lorenz' atmospheric model shows clearly how to separate the dynamics of the different atmospheric waves. However, its construction demonstrates that {\em any} initialization procedure must eventually alter the forecast--the time-scale of the divergence between the initialized and the uninitialized solutions is inevitable and is inversely proportional to the square of initial level of gravity-wave activity.

研究の動機と目的

• 低次元的大気モデルにおける慢い多様体と最も遅い不変多様体（SIM）の関係を明確化すること。
• 測定された初期条件を慢い多様体に射影することで、非物理的重力波を抑えることによる大気データの初期化の課題に対処すること。
• 予測歪みを最小限に抑える物理的根拠に基づく初期化スキームを構築し、初期重力波成分に起因する避けられない分散を考慮すること。
• L86モデルの正規形を分析することで、ロスビー波と重力波のダイナミクスの分離と初期化の本質的制限を明らかにすること。

提案手法

• 漸近解析と正規形理論を用いて、L86モデルの慢い多様体を構築し、速い運動と遅い運動を分離する。
• 慢い多様体とSIMを比較し、指数的に小さい項を除いて同一であることを示し、従来の近似スキームにおける重要な曖昧さを解消した。
• 特にヤコビ楕円関数を含む正規形変換に生じる積分を評価するために、等高線積分と留数計算を適用した。
• 発散級数を用いた慢い多様体の漸近展開を導出し、複素周波数における特異点によって収束半径が制限されることを示した。
• 正規形を用いて初期データを慢い多様体に射影するため、予測のずれを最小化する新しい初期化手法を提案した。
• 初期重力波振幅の二乗に反比例する時間スケールで、初期化済みと未初期化の解の間の分散が生じることを示した。

実験結果

リサーチクエスチョン

• RQ1 L86モデルにおいて、最も遅い不変多様体（SIM）と真の慢い多様体は数学的に同一であるか。もしそうであれば、どの程度の精度で同一であるか？
• RQ2 初期大気データを、予測歪みを最小限に抑えるように慢い多様体に射影する方法は何か？
• RQ3 非線形相互作用が初期化済みと未初期化の解の間の分散を引き起こす役割を果たすメカニズムは何か。また、これは初期重力波エネルギーにどのように依存するか？
• RQ4 正規形と漸近展開は、L86モデルにおいてロスビー波と重力波のダイナミクスの分離をどのように明らかにするか？
• RQ5 初期化精度の根本的限界は何か。初期波活動の観点から、どのように定量化されるか？

主な発見

• 最も遅い不変多様体（SIM）と真の慢い多様体は、指数的に小さい項を除いて同一であり、SIMが慢い多様体の有効な近似であることが確認された。
• 初期化済みと未初期化の解の間の分散の時間スケールは、初期重力波振幅の二乗に反比例しており、初期化精度の根本的限界を示している。
• L86モデルの正規形は、ロスビー波と重力波のダイナミクスの明確な分離を示しており、重力波振幅がロスビー波振幅の関数として表現されている。
• 慢い多様体の漸近級数は、複素周波数における特異点によって収束半径が制限されるため発散するが、正規形により物理的に意味のある近似が可能である。
• 提案された初期化スキームは、予測の変更を最小限に抑え、数値的に有効であることが示された。分散は非線形波相互作用に起因して避けられず生じる。
• 正規形積分の留数解析により、慢い多様体への主な寄与は複素平面における二重極から来るものであり、その留数は $ \frac{1}{mW^{*}} e^{-K^\prime / W^{*}} $ に比例することが判明した。