[论文解读] 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

一句话总结

本文通过分析洛伦兹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^{*}} $ 缩放。