准地转方程

地转运动是指由科里奥利力和水平压力梯度力之间的精确平衡产生的风，准地转 (Quasi-geostrophic, QG) 运动是指科里奥利力和压力梯度力几乎平衡的流动，但并不能排除惯性的影响。

## 起源

大气和海洋流动发生在水平长度尺度上，远大于垂直长度尺度，因此可以使用浅水方程来描述。罗斯贝数是一个无量纲数，它与科里奥利力的强度相比，表征惯性强度。准地转方程是小罗斯贝数极限下浅水方程的近似，因此惯性力比科里奥利力和压力小一个数量级。如果罗斯贝数等于 0，则准地转方程成为精确的地转方程。

准地转方程首先由儒勒·查尼提出。

## 单层 QG 方程的推导

在笛卡尔坐标中，地转风的分量是

${f_0} {v_g} = {\partial \Phi \over \partial x}$ (1a)

${f_0} {u_g} = - {\partial \Phi \over \partial y}$ (1b)

其中${\Phi}$是位势高度。

地转涡度

$${\zeta_g} = {\hat{\mathbf{k}} \cdot \nabla \times \mathbf{V_g}}$$

因此可以用位势高度表示为

${\zeta_g} = {{\partial v_g \over \partial x} - {\partial u_g \over \partial y} = {1 \over f_0} \left({ {\partial^2 \Phi \over \partial x^2} + {\partial^2 \Phi \over \partial y^2}}\right) = {1 \over f_0}{\nabla^2 \Phi}}$ (2)

式（2）可用于从已知位势高度场${\Phi (x,y)}$ 找到${\zeta_g (x,y)}$。也可以通过反转拉普拉斯算子从已知分布${\zeta_g}$来确定${\Phi}$。

准地转涡度方程可由下式得到${x}$和${y}$准地转动量方程的分量，然后可以从水平动量方程导出：

${D\mathbf{V} \over Dt} + f \hat{\mathbf{k}} \times \mathbf{V} = - \nabla \Phi$ (3)

式（3）中的物质导数定义为

${{D \over Dt} = {\left({\partial \over \partial t}\right)_p} + {\left({\mathbf{V} \cdot \nabla}\right)_p} + {\omega {\partial \over \partial p}}}$ (4)

其中${\omega = {Dp \over Dt}}$是运动后的压力变化。

水平速度${\mathbf{V}}$可以分为地转部分${\mathbf{V_g}}$和非地转部分${\mathbf{V_a}}$

${\mathbf{V} = \mathbf{V_g} + \mathbf{V_a}}$ (5)

准地转近似的两个重要假设是

- ${\mathbf{V_g} \gg \mathbf{V_a} }$，或者，更准确地说${{|\mathbf{V_a}| \over |\mathbf{V_g}|}} \sim O(\text{罗 斯 贝 数 })$ .

- β平面近似：${f = f_0 + \beta y}$, ${\frac{\beta y}{f_0} \sim O(\text{罗 斯 贝 数 })}$

第二个假设证明，在地转近似中，让科里奥利参数具有恒定值${f_0}$是合理的，并通过${f_0 + \beta y}$近似其在科里奥利力项中的变化。但是，由于运动后的加速度（在（1）中作为科里奥利力和压力梯度力之间的差值给出）取决于实际风与地转风的偏离，因此不允许简单地替换科里奥利力这一项中的地转速度。（3）式中的加速度可以重写为

${f \hat{\mathbf{k}} \times \mathbf{V} + \nabla \Phi} = {(f_0 + \beta y)\hat{\mathbf{k}} \times (\mathbf{V_g} + \mathbf{V_a}) - f_0 \hat{\mathbf{k}} \times \mathbf{V_g}} = {f_0 \hat{\mathbf{k}} \times \mathbf{V_a} + \beta y \hat{\mathbf{k}} \times \mathbf{V_g} }$ (6)

因此，近似水平动量方程具有形式

${D_g \mathbf{V_g} \over Dt} = {-f_0 \hat{\mathbf{k}} \times \mathbf{V_a} - \beta y \hat{\mathbf{k}} \times \mathbf{V_g}}$ (7)

用其分量表达方程（7）。

${{D_g u_g \over Dt} - {f_0 v_a} - {\beta y v_g} = 0}$ (8a)

${{D_g v_g \over Dt} + {f_0 u_a} + {\beta y u_g} = 0}$ (8b)

我们进行${{\partial (8b) \over \partial x} - {\partial (8a) \over \partial y}}$，并注意到地转风是无辐散的（即 ${\nabla \cdot \mathbf{V} = 0}$)，可得涡量方程为

${{D_g \zeta_g \over Dt} = - f_0 \left ({{\partial u_a \over \partial x}+{\partial v_a \over \partial y}} \right) - \beta v_g }$ (9)

因为${f}$只取决于${y}$ （亦即${{D_g f \over Dt} = \mathbf{V_g} \cdot \nabla f = \beta v_g}$) 并且地转风的散度基于连续性方程可以写成含${\omega}$的形式：

$${{\partial u_a \over \partial x}+{\partial v_a \over \partial y}+{\partial \omega \over \partial p}=0}$$

因此式 (9) 可以化为

${{\partial \zeta_g \over \partial t} = {-\mathbf{V_g} \cdot \nabla ({\zeta_g + f})} - {f_0 {\partial \omega \over \partial p}} }$ (10)

### 引入位势倾向

定义位势倾向${\chi = {\partial \Phi \over \partial t}}$，并且注意到偏微分可能被反转，等式（10）可以重写为${\chi}$的关系式

${{1 \over f_0}{\nabla^2 \chi} = {-\mathbf{V_g} \cdot \nabla \left({{1 \over f_0}{\nabla^2 \Phi} + f} \right)} + {f_0 {\partial \omega \over \partial p}}}$ (11)

等式（11）的右侧取决于变量${\Phi}$和${\omega}$ .依赖于这两个变量的类比方程可以从热力学能量方程导出

${{{ \left ({{\partial \over \partial t} + {\mathbf{V_g} \cdot \nabla}} \right) \left({-\partial \Phi \over \partial p} \right)}-\sigma \omega}={kJ \over p}}$ (12)

其中${\sigma = {-R T_0 \over p}{d \log \Theta_0 \over dp}}$和${\Theta_0}$是对应于基态温度的位温。在对流层中部，${\Theta_0}$ ≈ ${2.5 \times 10^{-6} \mathrm{m}{^2}\mathrm{Pa}^{-2}\mathrm{s}^{-2}}$ .

将 (12) 乘以${f_0 \over \sigma}$并对${p}$微分，结合${\chi}$的定义可得

${{{\partial \over \partial p} \left ({{f_0 \over \sigma}{\partial \chi \over \partial p}} \right )}=-{{\partial \over \partial p}\left({{f_0 \over \sigma}{\mathbf{V_g} \cdot \nabla}{\partial \Phi \over \partial p}}\right)}-{{f_0}{\partial \omega \over \partial p}}-{{f_0}{\partial \over \partial p}\left({kJ \over \sigma p}\right)}}$ (13)

简单起见，${J}$设为 0，消除 (11) 和 (13) 中的${\omega}$得出

${{ \left({\nabla^2 + {{\partial \over \partial p} \left({{f_0^2 \over \sigma}{\partial \over \partial p}}\right)}}\right){\chi}}=-{{f_0}{\mathbf{V_g} \cdot \nabla}\left({{{1 \over f_0}{\nabla^2 \Phi}}+f}\right)}-{{\partial \over \partial p}\left({{-}{f_0^2 \over \sigma}{\mathbf{V_g} \cdot \nabla}\left({\partial \Phi \over \partial p}\right)}\right)}}$ (14)

方程（14）通常被称为位势倾向方程。它将局部位势趋势（项 A）与涡度平流分布（项 B）和厚度平流（项 C）联系起来。

### 使用准地转位涡度的相同恒等式

使用微分的链式法则，C 项可以写为

${-{{\mathbf{V_g} \cdot \nabla}{\partial \over \partial p}\left({{f_0^2 \over \sigma}{\partial \Phi \over \partial p}}\right)}-{{f_0^2 \over \sigma}{\partial \mathbf{V_g} \over \partial p}{\cdot \nabla}{\partial \Phi \over \partial p}}}$ (15)

但基于热成风关系。

${{f_0{\partial \mathbf{V_g} \over \partial p}}={\hat{\mathbf{k}} \times \nabla \left({\partial \Phi \over \partial p} \right)}}$ .

换句话说， ${\partial \mathbf{V_g} \over \partial p}$垂直于${\nabla ({\partial \Phi \over \partial p})}$，式（15）中的第二项消失。

第一项可以与式（14）中的项 B 组合，当除以${f_0}$可以用守恒方程的形式表示

${{\left({{\partial \over \partial t}+{\mathbf{V_g} \cdot \nabla}}\right)q}={D_g q \over Dt}=0}$ (16)

其中${q}$是由下式定义的准地转位涡

${q = {{{1 \over f_0}{\nabla^2 \Phi}}+{f}+{{\partial \over \partial p}\left({{f_0 \over \sigma}{\partial \Phi \over \partial p}}\right)}}}$ (17)

方程(17)的三项从左到右分别是地转相对涡度、行星涡度和伸展涡度。

## 推论

当一个气团在大气中移动时，它的相对涡量、行星涡量和拉伸涡量可能会发生变化，但式（17）表明，随着地转运动，三者之和必须是守恒的。

式 (17) 可用于用已知高度场${\Phi}$找到${q}$。或者，它也可以用于预测给定初始分布的位势场的演变${\Phi}$和合适的边界条件通过使用反演过程。

更重要的是，准地转系统将五变量原始方程简化为一个方程系统，其中所有变量如${u_g}$, ${v_g}$和${T}$可以从位涡${q}$或高度场${\Phi}$导出。

另外，因为${\zeta_g}$和${\mathbf{V_g}}$都被定义为${\Phi(x,y,p,t)}$，涡量方程可用于诊断垂直运动，前提是两者的场${\Phi}$和${\partial \Phi \over \partial t}$是已知的。