112 年 國立中央大學地球科學學系地球物理碩士班一般生《大氣動力學》

此題主要考察學生對動量方程式的理解，以及如何從一般形式推導到特定平衡條件下的形式，並應用於實際天氣現象的解釋。

(a) 梯度風方程的推導：
首先，我們需要將給定的水平動量方程式在流線座標下展開。流線座標系是以氣流的速度方向為基礎建立的座標系，通常包含切線方向（與速度同向）和法線方向（指向曲率中心）。然而，題目要求推導的是梯度風方程，梯度風是假設氣壓梯度力、科氏力與離心力（或虛擬離心力，若考慮旋轉座標系）平衡，但科氏力在切線方向的分量與速度變化率有關，而在法線方向的分量與速度大小有關。題目給的方程式是無科氏力項的「非慣性座標」下的動量方程式，但為了推導梯度風，我們通常會從考慮科氏力的情況出發。

若題目意指從給定方程式推導，則應理解 $\frac{DV}{Dt}$ 為總加速度，包含隨時間變化 ($\frac{\partial V}{\partial t}$) 和隨位置變化 ($V \cdot \nabla V$)。$\nabla_p \Phi$ 是壓力梯度力。

然而，標準的梯度風方程是從考慮科氏力開始推導的。假設我們從以下更完整的形式開始（在考慮科氏力的情況下）：
$\frac{\partial V}{\partial t} + (V \cdot \nabla) V + f k \times V = -\nabla_p \Phi$
其中 $f$ 為科氏參數。

在流線座標系中，速度 $V$ 沿著切線方向，其大小為 $|V|$。加速度 $(V \cdot \nabla) V$ 可以分解為切線方向的加速度 $\frac{\partial |V|}{\partial t} + |V|\frac{\partial V}{\partial s}$（其中 $s$ 是沿著流線的距離）以及法線方向的向心加速度 $\frac{|V|^2}{R}$，其中 $R$ 是流線的曲率半徑。

梯度風假設是穩態（steady state, $\frac{\partial V}{\partial t} = 0$）且無垂直運動（或垂直運動不影響水平動量平衡）。因此，水平動量方程式在流線座標系下，沿著法線方向（指向曲率中心）的分量為：
$\frac{|V|^2}{R} - f|V| = -\frac{\partial_n \Phi}{\partial n}$
其中 $\frac{\partial_n \Phi}{\partial n}$ 是壓力梯度力沿著法線方向的分量，通常與垂直於流線方向的壓力梯度有關。

如果我們從題目給定的方程式 $\frac{DV}{Dt} + fk \times V = -\nabla_p \Phi$ 出發，並且假設 $\frac{DV}{Dt}$ 包含了所有非慣性力的效應，包括曲率效應。題目要求的是「梯度風方程」，這通常是指在穩態下，氣壓梯度力、科氏力與離心力（或其等效項）達成平衡。

另一種可能的解釋是，題目給的 $\frac{DV}{Dt}$ 包含了 $(V \cdot \nabla)V$ 項，並且我們將其在流線座標下展開，考慮其法線方向的分量。
假設速度 $V$ 的方向為切線方向，法線方向指向曲率中心。
$V = |V| \hat{t}$
$(V \cdot \nabla) V$ 在法線方向的分量為 $\frac{|V|^2}{R}$。
則動量方程式在法線方向的分量為：
$\frac{|V|^2}{R} + f|V| = -\frac{\partial_n \Phi}{\partial n} \quad \text{(北半球，離心力與科氏力同向)}$
或
$\frac{|V|^2}{R} - f|V| = -\frac{\partial_n \Phi}{\partial n} \quad \text{(北半球，離心力與科氏力反向)}$
這取決於 $\frac{DV}{Dt}$ 的具體定義。

然而，標準的梯度風方程是從慣性座標系下的動量方程推導出來的。
$\frac{\partial V}{\partial t} + (V \cdot \nabla) V + f k \times V = -\nabla_p \Phi$
在穩態下 ($\frac{\partial V}{\partial t} = 0$)，並考慮氣流沿著近似圓形等壓線運動，其曲率為 $R$。
法線方向（指向圓心）的方程為：
$\frac{|V|^2}{R} - f|V| = -\frac{\partial \Phi}{\partial n}$
其中 $\frac{\partial \Phi}{\partial n}$ 是重力位在垂直於速度方向（指向低壓區）的梯度。
如果氣流是順時針（北半球反氣旋）或逆時針（北半球氣旋），則 $R$ 的符號和 $\frac{\partial \Phi}{\partial n}$ 的方向需要仔細定義。

通常，梯度風方程寫為：
$|V_g| = \frac{1}{f} \frac{\partial \Phi}{\partial n}$
這是對於平行於等壓線的風，即地轉風。

梯度風方程是考慮了氣流的曲率效應，使得風速與地轉風不同。
在北半球，對於氣旋（逆時針流），氣壓梯度力指向圓心，科氏力向外，離心力向外。
$\frac{|V_g|^2}{R} + f|V_g| = \frac{1}{\rho} \frac{\partial p}{\partial n}$
在 $p$ 座標下，$-\nabla_p \Phi$ 相當於 $-\frac{1}{\rho} \nabla_p p$。
如果我們將壓力梯度力寫為 $\frac{\partial \Phi}{\partial n}$，那麼在氣旋中，法線方向（指向圓心）的方程是：
$\frac{|V_g|^2}{R} = \frac{\partial \Phi}{\partial n} - f|V_g|$
這其中 $\frac{\partial \Phi}{\partial n}$ 是沿著法線方向的重力位梯度。

在反氣旋（順時針流）中，法線方向（指向圓心）的方程是：
$\frac{|V_g|^2}{R} = f|V_g| - \frac{\partial \Phi}{\partial n}$
或
$\frac{|V_g|^2}{R} = - \frac{\partial \Phi}{\partial n} - f|V_g|$
這取決於 $\frac{\partial \Phi}{\partial n}$ 的定義。

一般而言，梯度風方程可以寫成關於 $|V_g|$ 的二次方程：
$|V_g|^2 \pm f R |V_g| - R \frac{\partial \Phi}{\partial n} = 0$
其中 $\pm$ 號取決於氣旋或反氣旋。
在氣旋中，科氏力提供向外的力，但要平衡向內的壓力梯度力，所以需要額外的向心力 $\frac{|V_g|^2}{R}$。
$f|V_g| - \frac{|V_g|^2}{R} = \frac{\partial \Phi}{\partial n} \quad (\text{北半球氣旋})$
$-f|V_g| - \frac{|V_g|^2}{R} = \frac{\partial \Phi}{\partial n} \quad (\text{北半球反氣旋})$
其中 $\frac{\partial \Phi}{\partial n}$ 是垂直於速度方向，指向低壓或高壓的梯度。

如果題目給定的方程式 $\frac{DV}{Dt} + fk \times V = -\nabla_p \Phi$ 中的 $\frac{DV}{Dt}$ 已經包含了曲率效應，即 $\frac{DV}{Dt}$ 在法線方向的分量為 $\frac{|V|^2}{R}$。
那麼在穩態下，法線方向的分量為：
$\frac{|V|^2}{R} + fk \times V \cdot \hat{n} = -\nabla_p \Phi \cdot \hat{n}$
在北半球，如果 $V$ 是逆時針，$\hat{n}$ 指向圓心，那麼 $k \times V$ 指向圓心外，所以 $fk \times V \cdot \hat{n} = -f|V|$。
$\frac{|V|^2}{R} - f|V| = -\frac{\partial \Phi}{\partial n}$

本題主要考察學生對於大氣運動中黏滯效應的處理方式，特別是從分子黏滯力到渦動黏滯力的概念轉換，以及雷諾平均在紊流情況下的應用。

(a) 雷諾分解與平均的方法與意義：
方法：
大氣運動，特別是中小尺度的紊流（turbulence），其運動軌跡極其複雜且難以直接求解。為了在數值模式中處理這些效應，我們採用雷諾分解（Reynolds decomposition）與平均（Reynolds averaging）。
對於任何一個變量（例如速度分量 $u$），我們可以將其分解為一個時間平均值（或空間平均值，取決於尺度）和一個與平均值偏差的脈動值：
$u = \bar{u} + u'$
其中 $\bar{u}$ 是平均速度， $u'$ 是速度的脈動分量，其時間平均值為零 ($\overline{u'} = 0$)。
將此分解應用於動量方程式中的各項，例如速度 $u$ 和 $v$：
$u = \bar{u} + u'$
$v = \bar{v} + v'$
將這些代入動量方程式（例如，以 $x$ 方向動量為例）：
$\frac{\partial u}{\partial t} + u \frac{\partial u}{\partial x} + v \frac{\partial u}{\partial y} + w \frac{\partial u}{\partial z} + fv = -\frac{1}{\rho}\frac{\partial p}{\partial x} + (\text{其他力})$
進行時間平均後，由於平均運算子的線性性質，以及脈動量的平均為零，會產生一些新的平均項，特別是來自速度項乘積的平均：
$\overline{u \frac{\partial u}{\partial x}} = \overline{(\bar{u} + u') \frac{\partial (\bar{u} + u')}{\partial x}} = \overline{(\bar{u} + u') (\frac{\partial \bar{u}}{\partial x} + \frac{\partial u'}{\partial x})}$
$= \overline{\bar{u}\frac{\partial \bar{u}}{\partial x} + \bar{u}\frac{\partial u'}{\partial x} + u'\frac{\partial \bar{u}}{\partial x} + u'\frac{\partial u'}{\partial x}}$
$= \bar{u}\frac{\partial \bar{u}}{\partial x} + \overline{u'\frac{\partial u'}{\partial x}}$ （因為 $\overline{\bar{u}\frac{\partial u'}{\partial x}}=0$ 且 $\overline{u'\frac{\partial \bar{u}}{\partial x}}=\frac{\partial \bar{u}}{\partial x}\overline{u'} = 0$）
$= \bar{u}\frac{\partial \bar{u}}{\partial x} + \frac{\partial \overline{u'u'}}{\partial x}$ （假設 $\frac{\partial \bar{u}}{\partial x}$ 不隨時間變化，或者進行空間平均）
這樣就產生了雷諾應力項 $\frac{\partial \overline{u'u'}}{\partial x}$。

意義：

1. 處理紊流效應：大氣中的許多運動，特別是中小尺度運動，具有高度的紊流特性。分子黏滯力（如牛頓黏滯定律）在宏觀尺度下非常微弱，無法描述紊流所帶來的動量混合和耗散效應。雷諾分解與平均後產生的雷諾應力項（例如 $\overline{u'v'}$）可以被看作是額外的「渦動黏滯力」，它描述了脈動運動在平均流中的動量傳遞。
2. 簡化模型：直接模擬所有尺度的紊流運動（直接數值模擬, DNS）計算量極其龐大，在實際應用中不可行。雷諾平均方法將動量方程式中的複雜脈動項歸結為平均量和一些需要額外參數化的雷諾應力項，從而大大簡化了模型，使其可以在可接受的計算資源下運行。
3. 引入渦動黏滯係數：為了參數化雷諾應力，通常引入渦動黏滯係數（eddy viscosity, $K_m$ 或 $\nu_t$）。類似於牛頓黏滯定律，假設雷諾應力與平均速度梯度成正比，例如：
   $\overline{u'v'} = -K_m \frac{\partial \bar{v}}{\partial x}$ （這是一種簡化的形式，更準確應為 $\overline{u'v'} = -K_{xy} \frac{\partial \bar{u}}{\partial y}$ 或其他形式）
   通過引入渦動黏滯係數，將難以求解的雷諾應力轉化為與平均流場相關的參數，這些參數可以通過實驗觀測或更精細的模擬來確定。

(b) 雷諾平均動量方程式（中緯度綜觀尺度運動）：
對於中緯度綜觀尺度運動，我們通常假設：

1. 穩態近似：$\frac{\partial}{\partial t} \approx 0$ （或認為時間尺度遠大於脈動尺度）。
2. 地轉近似：科氏力是主要的水平力之一。
3. 忽略分子黏滯力：分子黏滯力遠小於渦動黏滯力。
4. 壓力座標：通常在壓力座標系下分析。

以 $x$ 方向動量方程式為例，在壓力座標系下的雷諾平均形式為：
$\frac{\partial \bar{u}}{\partial t} + \bar{u}\frac{\partial \bar{u}}{\partial x} + \bar{v}\frac{\partial \bar{u}}{\partial y} + \bar{\omega}\frac{\partial \bar{u}}{\partial p} - f\bar{v} = -\frac{\partial \Phi}{\partial x} + \frac{\partial \tau_{xp}}{\partial p}$
其中 $\bar{\omega} = \frac{dp}{dt}$ 是垂直速度，$\Phi$ 是重力位， $f$ 是科氏參數。
$\frac{\partial \tau_{xp}}{\partial p}$ 代表了垂直方向的渦動黏滯力通量梯度。
$\tau_{xp}$ 是壓力座標下的渦動黏滯應力，它包含了雷諾應力在垂直方向的效應。

如果我們考慮水平方向的渦動黏滯力，方程式會更複雜，通常會包含 $\frac{\partial \overline{u'u'}}{\partial x} + \frac{\partial \overline{u'v'}}{\partial y}$ 等項。
然而，在中緯度綜觀尺度，水平渦動黏滯力通常被認為比科氏力和壓力梯度力小，因此常被忽略或參數化。
如果我們主要考慮垂直渦動黏滯力的影響，則方程式如上所示。

更完整的形式，考慮水平和垂直渦動黏滯力（在 $x$ 方向）：
$\frac{\partial \bar{u}}{\partial t} + \bar{u}\frac{\partial \bar{u}}{\partial x} + \bar{v}\frac{\partial \bar{u}}{\partial y} + \bar{\omega}\frac{\partial \bar{u}}{\partial p} - f\bar{v} = -\frac{\partial \Phi}{\partial x} + \frac{\partial}{\partial x}(\overline{u'u'}) + \frac{\partial}{\partial y}(\overline{u'v'}) + \frac{\partial}{\partial p}(\overline{u'\omega'})$

核心觀念

本題考查由水平動量方程式推導「垂直相對渦度方程式」，並解釋各項對局地垂直渦度的作用。

定義水平風場為

$$\boldsymbol{V}_h=(u,v),$$

垂直相對渦度為

$$\zeta=\frac{\partial v}{\partial x}-\frac{\partial u}{\partial y},$$

行星渦度為

$$f=2\Omega\sin\phi,$$

因此絕對垂直渦度為

$$\eta=\zeta+f.$$

採用高度座標下的水平動量方程式：

$$\frac{Du}{Dt}-fv = -\frac{1}{\rho}\frac{\partial p}{\partial x}+F_x,$$ $$\frac{Dv}{Dt}+fu = -\frac{1}{\rho}\frac{\partial p}{\partial y}+F_y,$$

其中

$$\frac{D}{Dt} = \frac{\partial}{\partial t} +u\frac{\partial}{\partial x} +v\frac{\partial}{\partial y} +w\frac{\partial}{\partial z}$$

為三維物質導數，$\rho$ 為密度，$p$ 為氣壓，$F_x,F_y$ 為摩擦力或其他非保守力的單位質量分量。

---

解題方法：對水平動量方程式取旋度

將南北向動量方程式對 $x$ 微分，將東西向動量方程式對 $y$ 微分，再以前者減去後者：

$$\frac{\partial}{\partial x} \left( \frac{Dv}{Dt}+fu \right) - \frac{\partial}{\partial y} \left( \frac{Du}{Dt}-fv \right) = \frac{\partial}{\partial x} \left( -\frac{1}{\rho}\frac{\partial p}{\partial y}+F_y \right) - \frac{\partial}{\partial y} \left( -\frac{1}{\rho}\frac{\partial p}{\partial x}+F_x \right).$$

1. 物質導數項

由

$$\zeta=\frac{\partial v}{\partial x}-\frac{\partial u}{\partial y}$$

可得

$$\frac{\partial}{\partial x}\left(\frac{Dv}{Dt}\right) - \frac{\partial}{\partial y}\left(\frac{Du}{Dt}\right) = \frac{D\zeta}{Dt} + \zeta \left( \frac{\partial u}{\partial x} + \frac{\partial v}{\partial y} \right) + \frac{\partial w}{\partial x}\frac{\partial v}{\partial z} - \frac{\partial w}{\partial y}\frac{\partial u}{\partial z}.$$

整理後，水平散度項與垂直風切項會分別形成渦度的伸長與傾斜效應。

2. 科氏力項

科氏力部分為

$$\frac{\partial(fu)}{\partial x} + \frac{\partial(fv)}{\partial y}.$$

在中緯度近似下，$f=f(y)$，且

$$\frac{\partial f}{\partial y}=\beta, \qquad \frac{\partial f}{\partial x}=0.$$

因此

$$\frac{\partial(fu)}{\partial x} + \frac{\partial(fv)}{\partial y} = f\left( \frac{\partial u}{\partial x} + \frac{\partial v}{\partial y} \right) +\beta v.$$

其中 $\beta v$ 是南北移動造成的行星渦度改變。

3. 壓力梯度項

壓力梯度項的水平旋度為

$$\frac{\partial}{\partial x} \left( -\frac{1}{\rho}\frac{\partial p}{\partial y} \right) - \frac{\partial}{\partial y} \left( -\frac{1}{\rho}\frac{\partial p}{\partial x} \right).$$

展開後得到

$$\frac{1}{\rho^2} \left( \frac{\partial \rho}{\partial x}\frac{\partial p}{\partial y} - \frac{\partial \rho}{\partial y}\frac{\partial p}{\partial x} \right).$$

此項可寫成

$$\frac{1}{\rho^2} \left(\nabla\rho\times\nabla p\right)_z,$$

即斜壓效應。

4. 摩擦力項

摩擦力對垂直渦度的貢獻為

$$\frac{\partial F_y}{\partial x} - \frac{\partial F_x}{\partial y} = \left(\nabla\times\boldsymbol{F}\right)_z.$$

---

垂直渦度方程式

綜合上述各項，得到高度座標下的垂直相對渦度方程式：

$$\boxed{ \frac{D\zeta}{Dt} = -(\zeta+f) \left( \frac{\partial u}{\partial x} + \frac{\partial v}{\partial y} \right) -\beta v +\frac{\partial w}{\partial y}\frac{\partial u}{\partial z} -\frac{\partial w}{\partial x}\frac{\partial v}{\partial z} +\frac{1}{\rho^2} \left( \frac{\partial \rho}{\partial x}\frac{\partial p}{\partial y} - \frac{\partial \rho}{\partial y}\frac{\partial p}{\partial x} \right) +\left(\nabla\times\boldsymbol{F}\right)_z }$$

其中

$$\frac{D\zeta}{Dt} = \frac{\partial\zeta}{\partial t} +u\frac{\partial\zeta}{\partial x} +v\frac{\partial\zeta}{\partial y} +w\frac{\partial\zeta}{\partial z}.$$

若題目一採用的是無摩擦、正壓且忽略 $\beta$ 效應的簡化水平動量方程式，上式可簡化為

本題主要考察學生對位渦概念的理解，以及其在解釋大氣現象中的應用。

(a) Ertel 位渦與淺水位渦的定義與異同：

Ertel 位渦 (Potential Vorticity, PV)

• 物理定義：Ertel 位渦是絕熱、無摩擦、無黏滯的流體中，一個守恆的向量量。其定義為：
  $PV = (\zeta + f) \frac{\partial \theta}{\partial z}$
  在壓力座標系下，更常見的定義是：
  $PV = (\zeta + f) \frac{\partial \theta}{\partial p}$
  其中，$\zeta$ 是相對渦度， $f$ 是科氏參數，$\theta$ 是位溫，$\frac{\partial \theta}{\partial p}$ 是靜力穩定度（或稱為溫標梯度）。
  更一般的向量形式為 $PV = (\nabla \times V + f \mathbf{k}) \cdot \nabla \theta$。
• 物理意義：Ertel 位渦是描述流體運動的一個非常重要的守恆量。它結合了絕對渦度 ($\zeta + f$) 和溫標梯度 ($\frac{\partial \theta}{\partial p}$)。
  • 絕對渦度：代表了流體質點的旋轉程度。
  • 溫標梯度：代表了流體質點在垂直方向上的梯度。
    Ertel 位渦守恆意味著，如果一個流體單元（質點）的 PV 在初始時刻是一個定值，那麼它在運動過程中將保持這個值。這對於大尺度大氣運動的診斷和預報非常重要，因為它約束了流體的運動。
    PV 的守恆性可以用來理解很多天氣現象，例如 Rossby 波的傳播、大尺度氣旋的發展等。
    PV 的單位是 $K \cdot m^2 / (s \cdot kg)$ 或 $m^2 / (s \cdot kg)$ (在 SI 單位制下，如果 $\theta$ 歸一化)。

淺水位渦 (Shallow-Water Vorticity)

• 物理定義：淺水位系統是指一個厚度很薄（相對於水平尺度）的流體層。在這種系統中，水平運動是主要的，垂直運動可以忽略。
• 淺水位渦的定義：在淺水位系統中，垂直渦度 $\zeta = \frac{\partial v}{\partial x} - \frac{\partial u}{\partial y}$ 的守恆律是：
  $\frac{d\zeta}{dt} + \zeta (\nabla \cdot V) + (\nabla V \times \nabla h)_z = 0$
  其中 $h$ 是流體層的厚度（或水位高度），$\frac{d\zeta}{dt}$ 是總體導數，$\nabla \cdot V$ 是水平速度的散度。
  如果系統是不可壓縮的 ($\nabla \cdot V = 0$) 且沒有垂直運動，則：
  $\frac{d\zeta}{dt} = 0$
  或者，在考慮了 $\nabla \cdot V$ 的情況下，更準確的守恆量是 淺水位位渦 (Shallow-Water Potential Vorticity)：
  $PV_{SW} = \frac{\zeta + f}{h}$
  其守恆律為：
  $\frac{d}{dt} \left( \frac{\zeta + f}{h} \right) = 0$
  這意味著 $\frac{d\zeta}{dt} + f \frac{d}{dt}(\frac{1}{h}) = 0$。
  利用連續方程 $\frac{dh}{dt} + h (\nabla \cdot V) = 0$，即 $\frac{1}{h}\frac{dh}{dt} = -\nabla \cdot V$。
  $\frac{d}{dt}(\frac{1}{h}) = -\frac{1}{h^2}\frac{dh}{dt} = -\frac{1}{h}(-\nabla \cdot V) = \frac{\nabla \cdot V}{h}$。
  所以，$\frac{d}{dt} \left( \frac{\zeta + f}{h} \right) = \frac{1}{h}\frac{d\zeta}{dt} + (\zeta+f) \frac{\nabla \cdot V}{h} = \frac{1}{h}(\frac{d\zeta}{dt} + (\zeta+f)\nabla \cdot V) = 0$。
  這與上面的渦度守恆形式是等價的。
• 物理意義：淺水位渦描述了在一個薄流體層中，絕對渦度與厚度的比例是一個守恆量。當流體層變薄時（例如在氣旋發展中，空氣匯聚導致厚度減小），為了保持 PV 守恆，絕對渦度必須增加。這解釋了為何氣旋的中心渦度會增強。

異同：

• 相同點：
  • 兩者都是描述流體運動的守恆量（在特定條件下）。
  • 兩者都與渦度（絕對渦度）有關，並且渦度的變化是其核心內容。
  • 兩者都與流體的伸長/壓縮（或厚度變化）有關。

本題主要考察學生對斜壓不穩定性的理解，以及兩層準地轉模式在斜壓波發展中的應用。

(a) 斜壓不穩定與兩層準地轉模式：

何謂斜壓不穩定 (Baroclinic Instability)？

斜壓不穩定是指在中緯度大氣中，由氣壓梯度和溫度梯度（即斜壓性）所驅動的一種不穩定現象。在斜壓性大氣中，等壓面與等溫面（或等位溫面）不平行。這種不平行性導致了氣壓梯度力和溫度梯度力之間的差異，使得能量可以從平均流中轉移到擾動中，從而驅動擾動（如氣旋和反氣旋）的發展。

更具體地說，在斜壓大氣中，如果存在一個平均的西風，其速度隨高度增加（例如，一個溫帶氣旋的結構，低層冷空氣，高層暖空氣，或者低層暖舌，高層冷舌），那麼就可能發生斜壓不穩定。這種不穩定性會導致擾動（波）的發展，這些擾動具有傾斜的結構，並能夠將高層的暖空氣輸送到高緯度，將低層的冷空氣輸送到低緯度，從而平均大氣的溫度梯度。

兩層準地轉模式 (Two-Layer Quasi-Geostrophic Model)

兩層準地轉模式是一個簡化的數學模型，用於研究斜壓不穩定性。它假設大氣由兩層流體組成，中間有一光滑的界面，上下兩層的風和溫度 are assumed to be constant within each layer. 準地轉近似意味著我們忽略了非地轉效應（如加速度項），主要考慮了地轉平衡、散度、科氏力、壓力梯度力和渦度守恆。

在兩層準地轉模式中，斜壓不穩定性的發展與兩層之間的相對運動以及垂直風切變有關。擾動的發展通常是通過斜壓波（baroclinic waves）的形式出現的。

最容易發展的斜壓波波長：

在標準的兩層準地轉模式分析中，可以推導出斜壓波的頻率與波長之間的關係。
對於一個具有波数 $k$（對應波長 $\lambda = 2\pi/k$）的擾動，其增長率（imaginary part of frequency）可以表示出來。
在忽略了科氏參數 $f$ 的影響（或將其納入平均流），或者在考慮了標準科氏參數的情況下，最容易發展的斜壓波的波長 $\lambda_m$ 通常與模式的參數有關，例如大氣的平均尺度、溫度梯度、科氏參數等。

在一個典型的兩層準地轉模式分析中，最容易發展的斜壓波的波長 $\lambda_m$ 大約是：
$\lambda_m \approx \frac{2\pi L}{\sqrt{f^2/gH}}$
其中 $L$ 是大氣的特徵尺度（例如，Rossby 半徑），$f$ 是科氏參數，$g$ 是重力加速度，$H$ 是大氣的特徵厚度。

更具體的，對於一個具有平均西風 $U$ 的兩層模式，最容易發展的斜壓波的波長 $\lambda_m$ 與 $U$ 和 $f$ 有關。
在一種簡化的分析中，最容易發展的波長與 Rossby 半徑 $\approx \frac{U}{f}$ 有關。

一個常見的結果是，最容易發展的斜壓波的波長約為 3000-4000 公里 (km)。
這個波長對應於中緯度大氣中天氣系統（如溫帶氣旋）的典型尺度。
這個波長也可以表示為 $\lambda_m = \frac{2\pi \bar{U}}{f}$，其中 $\bar{U}$ 是平均風速，但這個公式更適用於 Rossby 波。

在兩層準地轉模式的標準分析中，最容易發展的斜壓波的波長 $\lambda_m$ 約為：
$\lambda_m \approx 2\pi \sqrt{\frac{gH}{f^2}}$
其中 $H$ 是大氣的有效厚度。如果取 $H \sim 10^3 km$，$f \sim 10^{-4} s^{-1}$，$g \sim 10 m/s^2$，則：
$\lambda_m \approx 2\pi \sqrt{\frac{10^3 \times 10^3 m \times 10 m/s^2}{(10^{-4} s^{-1})^2}} = 2\pi \sqrt{\frac{10^{10}}{10^{-8}}} m = 2\pi \sqrt{10^{18}} m = 2\pi \times 10^9 m = 20000 \pi km \approx 62800 km$。
這個結果不對。

我們需要更標準的兩層模式結果。
對於一個兩層模式，假設平均西風為 $U$。
斜壓不穩定性增長率 $\sigma$ 的最大值對應的波長 $\lambda_m$。
一個常見的結果是，最容易發展的斜壓波的波長 $\lambda_m$ 大約是：

本題主要考察學生對斜壓不穩定性的理解，以及兩層準地轉模式在斜壓波發展中的應用。

(a) 斜壓不穩定與兩層準地轉模式：

何謂斜壓不穩定 (Baroclinic Instability)？

斜壓不穩定是指在中緯度大氣中，由氣壓梯度和溫度梯度（即斜壓性）所驅動的一種不穩定現象。在斜壓性大氣中，等壓面與等溫面（或等位溫面）不平行。這種不平行性導致了氣壓梯度力與溫度梯度力之間的差異，使得能量可以從平均流中轉移到擾動中，從而驅動擾動（如氣旋和反氣旋）的發展。

更具體地說，在斜壓大氣中，如果存在一個平均的西風，其速度隨高度增加（例如，一個溫帶氣旋的結構，低層冷空氣，高層暖空氣，或者低層暖舌，高層冷舌），那麼就可能發生斜壓不穩定。這種不穩定性會導致擾動（波）的發展，這些擾動具有傾斜的結構，並能夠將高層的暖空氣輸送到高緯度，將低層的冷空氣輸送到低緯度，從而平均大氣的溫度梯度。

兩層準地轉模式 (Two-Layer Quasi-Geostrophic Model)

兩層準地轉模式是一個簡化的數學模型，用於研究斜壓不穩定性。它假設大氣由兩層流體組成，中間有一光滑的界面，上下兩層的風和溫度 are assumed to be constant within each layer. 準地轉近似意味著我們忽略了非地轉效應（如加速度項），主要考慮了地轉平衡、散度、科氏力、壓力梯度力以及渦度守恆。

在兩層準地轉模式中，斜壓不穩定性的發展與兩層之間的相對運動以及垂直風切變有關。擾動的發展通常是通過斜壓波（baroclinic waves）的形式出現的。

最容易發展的斜壓波波長：

在標準的兩層準地轉模式分析中，可以推導出斜壓波的頻率與波長之間的關係。
對於一個具有波数 $k$（對應波長 $\lambda = 2\pi/k$）的擾動，其增長率（imaginary part of frequency）可以表示出來。
在忽略了科氏參數 $f$ 的影響（或將其納入平均流），或者在考慮了標準科氏參數的情況下，最容易發展的斜壓波的波長 $\lambda_m$ 通常與模式的參數有關，例如大氣的平均尺度、溫度梯度、科氏參數等。

在一個典型的兩層準地轉模式分析中，最容易發展的斜壓波的波長 $\lambda_m$ 大約是：
$\lambda_m \approx \frac{2\pi L}{\sqrt{f^2/gH}}$
其中 $L$ 是大氣的特徵尺度（例如，Rossby 半徑），$f$ 是科氏參數，$g$ 是重力加速度，$H$ 是大氣的特徵厚度。

更具體的，對於一個具有平均西風 $U$ 的兩層模式，最容易發展的斜壓波的波長 $\lambda_m$ 與 $U$ 和 $f$ 有關。
在一種簡化的分析中，最容易發展的波長與 Rossby 半徑 $\approx \frac{U}{f}$ 有關。

本題主要考察學生對斜壓不穩定性的理解，以及兩層準地轉模式在斜壓波發展中的應用。

(a) 斜壓不穩定與兩層準地轉模式：

何謂斜壓不穩定 (Baroclinic Instability)？

斜壓不穩定是指在中緯度大氣中，由氣壓梯度和溫度梯度（即斜壓性）所驅動的一種不穩定現象。在斜壓性大氣中，等壓面與等溫面（或等位溫面）不平行。這種不平行性導致了氣壓梯度力與溫度梯度力之間的差異，使得能量可以從平均流中轉移到擾動中，從而驅動擾動（如氣旋和反氣旋）的發展。

更具體地說，在斜壓大氣中，如果存在一個平均的西風，其速度隨高度增加（例如，一個溫帶氣旋的結構，低層冷空氣，高層暖空氣，或者低層暖舌，高層冷舌），那麼就可能發生斜壓不穩定。這種不穩定性會導致擾動（波）的發展，這些擾動具有傾斜的結構，並能夠將高層的暖空氣輸送到高緯度，將低層的冷空氣輸送到低緯度，從而平均大氣的溫度梯度。

兩層準地轉模式 (Two-Layer Quasi-Geostrophic Model)

兩層準地轉模式是一個簡化的數學模型，用於研究斜壓不穩定性。它假設大氣由兩層流體組成，中間有一光滑的界面，上下兩層的風和溫度 are assumed to be constant within each layer. 準地轉近似意味著我們忽略了非地轉效應（如加速度項），主要考慮了地轉平衡、散度、科氏力、壓力梯度力以及渦度守恆。

在兩層準地轉模式中，斜壓不穩定性的發展與兩層之間的相對運動以及垂直風切變有關。擾動的發展通常是通過斜壓波（baroclinic waves）的形式出現的。

最容易發展的斜壓波波長：

在標準的兩層準地轉模式分析中，可以推導出斜壓波的頻率與波長之間的關係。
對於一個具有波数 $k$（對應波長 $\lambda = 2\pi/k$）的擾動，其增長率（imaginary part of frequency）可以表示出來。
在忽略了科氏參數 $f$ 的影響（或將其納入平均流），或者在考慮了標準科氏參數的情況下，最容易發展的斜壓波的波長 $\lambda_m$ 通常與模式的參數有關，例如大氣的平均尺度、溫度梯度、科氏參數等。

在一個典型的兩層準地轉模式分析中，最容易發展的斜壓波的波長 $\lambda_m$ 大約是：
$\lambda_m \approx \frac{2\pi L}{\sqrt{f^2/gH}}$
其中 $L$ 是大氣的特徵尺度（例如，Rossby 半徑），$f$ 是科氏參數，$g$ 是重力加速度，$H$ 是大氣的特徵厚度。

更具體的，對於一個具有平均西風 $U$ 的兩層模式，最容易發展的斜壓波的波長 $\lambda_m$ 與 $U$ 和 $f$ 有關。
在一種簡化的分析中，最容易發展的波長與 Rossby 半徑 $\approx \frac{U}{f}$ 有關。

本題主要考察學生對準地轉重力位趨勢方程的理解，以及如何利用地轉渦度平流和厚度平流來解釋 500 hPa 槽脊系統的演變。

首先，我們需要寫出準地轉重力位趨勢方程。這個方程描述了在準地轉近似下，某一等壓面（例如 500 hPa）上的重力位高度隨時間的變化率。它通常可以表示為：
$\frac{\partial \Phi}{\partial t} = -V_g \cdot \nabla \Phi - \omega \frac{\partial \Phi}{\partial p} + \text{其他項}$
在準地轉理論中，我們將重力位趨勢與渦度變化聯繫起來。
一個常用的準地轉重力位趨勢方程（在壓力座標下）是：
$\frac{\partial \Phi}{\partial t} = - \mathbf{V}_g \cdot \nabla_p \Phi - \omega \frac{\partial \Phi}{\partial p}$
其中 $\mathbf{V}_g$ 是地轉風，$\omega = dp/dt$ 是垂直速度，$\nabla_p \Phi$ 是壓力梯度力。
利用地轉關係 $\mathbf{V}_g = \frac{1}{f} \mathbf{k} \times \nabla_p \Phi$ 和 $\frac{\partial \Phi}{\partial p} = -\frac{RT}{p}$ (理想氣體)。

更為關鍵的是，重力位趨勢方程可以表示為渦度變化率的函數。
通過聯繫重力位和渦度，可以得到：
$\frac{\partial \zeta_g}{\partial t} = - \mathbf{V}_g \cdot \nabla_p (\zeta_g + f) - (\nabla \cdot \mathbf{V})_g$
其中 $\zeta_g$ 是地轉渦度，$\mathbf{V}_g$ 是地轉風。

一個更常用的形式，將重力位趨勢與渦度聯繫起來，是通過絕熱、無摩擦、不可壓縮的準地轉近似，可以推導出：
$\frac{\partial \Phi}{\partial t} \bigg|_{p=const} = \frac{1}{f_0} \left( - \mathbf{V}_g \cdot \nabla_p (\zeta_g + f) \right) - \frac{1}{f_0} \left( \omega \frac{\partial \Phi}{\partial p} \right)$
其中 $f_0$ 是平均的科氏參數。
我們知道 $\zeta_g + f$ 是絕對渦度。
$-\mathbf{V}_g \cdot \nabla_p (\zeta_g + f)$ 是絕對渦度的水平平流。
$-\omega \frac{\partial \Phi}{\partial p}$ 是垂直運動對重力位的影響。

另一種表達方式是，重力位趨勢與渦度的聯繫：
$\frac{\partial \Phi}{\partial t} \bigg|_{p=const} = - \frac{1}{f_0} \mathbf{V}_g \cdot \nabla_p (\zeta_g + f) + \text{垂直運動項}$
如果我們只考慮渦度平流對重力位趨勢的影響，則：
$\frac{\partial \Phi}{\partial t} \approx - \frac{1}{f_0} (\text{渦度平流})$
負號表示，如果渦度平流是正的（將高渦度區域的空氣平流到某處），則該處的重力位會下降（趨勢為負）。

然而，題目要求討論「地轉渦度平流」和「差異性的厚度平流」。
這通常來自於準地轉的權衡方程 (balance equation) 或重力位趨勢方程的另一種形式。
重力位趨勢方程可以寫成：
$\frac{\partial \Phi}{\partial t} = - \mathbf{V}_g \cdot \nabla_p \Phi - \omega \frac{\partial \Phi}{\partial p}$
利用地轉關係 $\nabla_p \Phi = f \mathbf{k} \times \mathbf{V}_g$，則 $-\mathbf{V}_g \cdot \nabla_p \Phi = -\mathbf{V}_g \cdot (f \mathbf{k} \times \mathbf{V}_g) = 0$ (因為 $\mathbf{V}_g$ 與 $\mathbf{k} \times \mathbf{V}_g$ 垂直)。
所以，
$\frac{\partial \Phi}{\partial t} = - \omega \frac{\partial \Phi}{\partial p}$
這看起來過於簡化。

正確的重力位趨勢方程（準地轉）可以表示為：
$\frac{\partial \Phi}{\partial t} = \frac{1}{f_0} \left( - \mathbf{V}_g \cdot \nabla_p (\zeta_g + f) \right) - \frac{1}{f_0} \left( \omega \frac{\partial \Phi}{\partial p} \right)$
此處 $\zeta_g$ 是地轉渦度，$\mathbf{V}_g$ 是地轉風，$\omega$ 是垂直速度。
$\frac{\partial \Phi}{\partial p}$ 是靜力穩定度項。
$\frac{\partial \Phi}{\partial t}$ 是重力位趨勢。

題目提到的「地轉渦度平流」是 $-\mathbf{V}_g \cdot \nabla_p (\zeta_g + f)$。
題目提到的「差異性的厚度平流」與垂直運動和溫度有關。
在準地轉近似下，厚度（例如 1000 hPa 和 500 hPa 之間的距離）與溫度成正比。
厚度變化率（例如 500-1000 hPa 的厚度）與溫度平流有關。
$\frac{\partial h_{500/1000}}{\partial t} \approx - \mathbf{V}_g \cdot \nabla (\bar{T})$。

讓我們採用一個更直接的表達形式，將重力位趨勢與渦度平流和厚度平流聯繫起來。
準地轉權衡方程 (Balance Equation) 和重力位趨勢方程結合，可以得到：
$\frac{\partial \Phi}{\partial t} = \frac{1}{f_0} \left( - \mathbf{V}_g \cdot \nabla_p (\zeta_g + f) \right) + \frac{1}{f_0} \left( \omega \frac{\partial \Phi}{\partial p} \right)$
其中：

1. 第一項：$\frac{1}{f_0} \left( - \mathbf{V}_g \cdot \nabla_p (\zeta_g + f) \right)$。

   • $-\mathbf{V}_g \cdot \nabla_p (\zeta_g + f)$ 是地轉渦度平流。
   • 它表示地轉風將絕對渦度平流到某個區域。
   • 如果某區域的地轉風將「高渦度」區域的空氣平流進來（例如，在氣旋的左側，風從高渦度區吹向低渦度區），則渦度平流為負，導致 $\frac{\partial \Phi}{\partial t}$ 為負，即重力位趨勢為負，表明該區域的重力位在下降。
   • 如果某區域的地轉風將「低渦度」區域的空氣平流進來（例如，在反氣旋的右側），則渦度平流為正，導致 $\frac{\partial \Phi}{\partial t}$ 為正，即重力位趨勢為正，表明該區域的重力位在上升。

2. 第二項：$\frac{1}{f_0} \left( \omega \frac{\partial \Phi}{\partial p} \right)$。

   • $\omega = dp/dt$ 是垂直速度。
   • $\frac{\partial \Phi}{\partial p}$ 是靜力穩定度。
   • $\omega \frac{\partial \Phi}{\partial p}$ 代表了垂直運動對重力位的影響。
   • 如果 $\omega > 0$ (上升運動，空氣壓力降低)，且 $\frac{\partial \Phi}{\partial p} < 0$ (通常情況下，高層壓力較低，重力位梯度為負)，則此項為負，導致 $\frac{\partial \Phi}{\partial t} < 0$。
   • 如果 $\omega < 0$ (下沉運動，空氣壓力升高)，則此項為正，導致 $\frac{\partial \Phi}{\partial t} > 0$。
   • 這個項與「厚度平流」有關。因為大氣的厚度（例如 500-1000 hPa 之間的距離）與平均溫度成正比。溫度平流會影響垂直運動，進而影響重力位趨勢。
   • 差異性的厚度平流：這通常指的是，在不同的高度層，水平溫度平流（或厚度平流）不同，導致垂直運動的產生。例如，暖空氣平流到某區域，會使該區域的厚度增加，導致上升運動；冷空氣平流到某區域，會使厚度減小，導致下沉運動。
   • 因此，$\omega \frac{\partial \Phi}{\partial p}$ 這一項，可以看作是從溫度平流/厚度平流的效應。

現在，我們將這些概念應用於 500 hPa 的槽(trough)和脊(ridge)系統。

500 hPa 槽 (Trough) 系統：
槽是高空等壓面上氣壓（或重力位）較低、等高線向低緯度（或低值區）凸出的區域。
在槽的前部（東側）：

• 地轉渦度平流：通常是將低渦度的空氣平流進來（或將高渦度的空氣平流出去）。例如，在一個發展中的氣旋系統中，槽的前部常常是相對低渦度的區域，或者風從低渦度區吹向高渦度區。如果我們考慮地轉風 $\mathbf{V}_g$ 吹向渦度增加的區域，則 $-\mathbf{V}_g \cdot \nabla (\zeta_g+f)$ 可能為負值（渦度平流為負）。
  • 如果渦度平流為負 (將低渦度平流進來)，則 $-\mathbf{V}_g \cdot \nabla (\zeta_g+f) < 0$。
  • 根據方程 $\frac{\partial \Phi}{\partial t} \approx \frac{1}{f_0} \left( - \mathbf{V}_g \cdot \nabla_p (\zeta_g + f) \right) + \dots$
  • 如果渦度平流為負，則 $\frac{\partial \Phi}{\partial t}$ 為負，重力位趨勢為負，表明 500 hPa 高度下降。
• 差異性的厚度平流（垂直運動）：在槽的前部，通常有暖平流（暖空氣從低緯度平流到高緯度），這導致 500-1000 hPa 的厚度增加，引起上升運動 ($\omega > 0$)。