109 年 國立臺灣大學大氣系碩士班甲組《大氣動力學》

本題旨在考驗學生對大氣熱力學基本概念的理解，特別是熱力風方程式的推導及其物理意義。熱力風方程式描述了在水平溫度梯度存在的情況下，大氣中垂直風速的變化，是理解斜壓大氣結構的關鍵。

【推導過程】
我們從流體靜力平衡和水平運動方程開始。
假設流體靜力平衡：
$\frac{\partial p}{\partial z} = -\rho g$
其中 $p$ 為壓力，$z$ 為高度，$\rho$ 為密度，$g$ 為重力加速度。

水平運動方程（忽略科氏力以外的力，如摩擦力，並考慮壓力梯度力）：
$\frac{\partial u}{\partial t} + u \frac{\partial u}{\partial x} + v \frac{\partial u}{\partial y} + w \frac{\partial u}{\partial z} - f v = -\frac{1}{\rho} \frac{\partial p}{\partial x}$
$\frac{\partial v}{\partial t} + u \frac{\partial v}{\partial x} + v \frac{\partial v}{\partial y} + w \frac{\partial v}{\partial z} + f u = -\frac{1}{\rho} \frac{\partial p}{\partial y}$
其中 $u, v, w$ 分別為 $x, y, z$ 方向的速度分量，$f$ 為科氏參數。

在準地轉近似下，我們忽略時間變化項和水平平流項，且認為水平運動主要由壓力梯度力和科氏力平衡，即地轉平衡：
$-f v = -\frac{1}{\rho} \frac{\partial p}{\partial x} \quad \implies \quad v = \frac{1}{\rho f} \frac{\partial p}{\partial x}$
$f u = -\frac{1}{\rho} \frac{\partial p}{\partial y} \quad \implies \quad u = -\frac{1}{\rho f} \frac{\partial p}{\partial y}$

現在我們考慮垂直方向的速度變化，即 $w$ 的變化。從連續性方程（不可壓縮流體）或更一般的質量守恆方程，結合理想氣體定律，可以推導出與溫度相關的項。然而，熱力風方程式的推導通常是從熱力學第一定律出發，結合流體靜力平衡來獲得。

考慮等壓面上的運動。在等壓面上，壓力梯度力與重力垂直分量平衡。如果存在溫度梯度，則等壓面會傾斜。
從流體靜力平衡方程，我們有：
$\frac{\partial p}{\partial z} = -\rho g$
考慮理想氣體定律：$\rho = \frac{p}{R T}$，其中 $R$ 是乾空氣的氣體常數，$T$ 是絕對溫度。
$\frac{\partial p}{\partial z} = -\frac{p}{R T} g \quad \implies \quad \frac{\partial (\ln p)}{\partial z} = -\frac{g}{R T}$
$\frac{\partial T}{\partial z} = -\Gamma_d = -\frac{g}{c_p}$
其中 $\Gamma_d$ 是乾絕熱直減率，$c_p$ 是定壓比熱。

現在我們考慮壓力場和溫度場的關係。
$\frac{\partial}{\partial x} \left( \frac{\partial p}{\partial z} \right) = \frac{\partial}{\partial z} \left( \frac{\partial p}{\partial x} \right)$
$\frac{\partial}{\partial x} (-\rho g) = \frac{\partial}{\partial z} (\rho f v)$
$-\frac{\partial \rho}{\partial x} g = f \left( v \frac{\partial \rho}{\partial z} + \rho \frac{\partial v}{\partial z} \right)$
使用理想氣體定律 $\rho = \frac{p}{RT}$：
$-\frac{\partial}{\partial x} \left( \frac{p}{RT} \right) g = f \left( \frac{p}{RT} \frac{\partial \rho}{\partial z} + \rho \frac{\partial v}{\partial z} \right)$
$-\frac{g}{RT} \frac{\partial p}{\partial x} + \frac{gp}{RT^2} \frac{\partial T}{\partial x} = f \left( \frac{p}{RT} \frac{\partial \rho}{\partial z} + \rho \frac{\partial v}{\partial z} \right)$
這個推導比較複雜。更直接的方法是從等壓面和等溫面的傾斜關係出發。

考慮等壓面 $p = \text{const}$。其法向量與重力方向（垂直向下）的夾角 $\alpha$ 滿足：
$\tan \alpha = \frac{1}{g} \frac{\partial p}{\partial h}$
其中 $h$ 是水平方向。

或者，從熱力學第一定律和運動方程結合。
考慮絕熱過程，溫度變化率為：
$\frac{dT}{dt} = \frac{w R T}{c_p p} \Gamma_d$
這也不是直接推導熱力風。

標準推導方法：
從流體靜力平衡和水平運動方程出發。
$\frac{\partial p}{\partial z} = -\rho g$
$\frac{\partial u}{\partial t} + \mathbf{v} \cdot \nabla u - f v = -\frac{1}{\rho} \frac{\partial p}{\partial x}$
$\frac{\partial v}{\partial t} + \mathbf{v} \cdot \nabla v + f u = -\frac{1}{\rho} \frac{\partial p}{\partial y}$
在準地轉近似下，我們主要關注的是風速隨高度的變化。
考慮垂直方向的速度分量 $w$。
從流體靜力平衡，我們有：
$\frac{\partial}{\partial x} \left( \frac{\partial p}{\partial z} \right) = \frac{\partial}{\partial z} \left( \frac{\partial p}{\partial x} \right)$
$\frac{\partial}{\partial x} (-\rho g) = \frac{\partial}{\partial z} (\rho f v)$
$-g \frac{\partial \rho}{\partial x} = f \left( v \frac{\partial \rho}{\partial z} + \rho \frac{\partial v}{\partial z} \right)$
這仍然很複雜。

簡潔推導：
考慮等壓面。在等壓面上，壓力梯度力與重力垂直分量平衡。
$\nabla_p p = 0$
$\nabla_p = \nabla - \frac{\nabla z}{\partial z / \partial p} \frac{\partial}{\partial p}$
$\nabla_p = \nabla - \mathbf{k} \times \nabla \psi$
其中 $\psi$ 是地衡勢。

或者，考慮等壓面方程：
$\frac{\partial p}{\partial z} = -\rho g$

本題考查的是大氣動力學中的核心概念——準地轉近似（Quasi-Geostrophic Approximation）。準地轉近似是分析大尺度天氣系統（如鋒面、氣旋、反氣旋）時常用的簡化方法，它基於幾個關鍵的無量綱參數和物理假設。

【準地轉近似的基礎】
準地轉理論的核心是假設大氣運動是「近似」地轉平衡的，並且其偏差（即非地轉成分）是小的，但正是這些小的偏差驅動了天氣系統的演變，特別是垂直運動和位勢的變化。

【關鍵的無量綱參數】
準地轉近似的有效性可以通過幾個無量綱參數來判斷：

1. 羅 स्वातंत्र्य數 (Rossby Number, Ro)：
   $Ro = \frac{U}{f L}$
   其中 $U$ 是特徵水平速度，$f$ 是科氏參數，$L$ 是特徵水平尺度。
   對於大尺度天氣系統（如行星尺度波動、中尺度氣旋），其水平尺度 $L$ 很大，而特徵速度 $U$ 相對較小，因此 $Ro \ll 1$。這意味著慣性力（$U^2/L$）遠小於科氏力（$fU$）。準地轉理論就是基於 $Ro \ll 1$ 的假設。

2. 馬赫數 (Mach Number, Ma)：
   $Ma = \frac{U}{c}$
   其中 $c$ 是聲速。由於大氣中的風速遠小於聲速，Ma 總是 $\ll 1$。這保證了我們可以忽略可壓縮性效應，並使用流體靜力平衡。

3. 弗勞德數 (Froude Number, Fr) 或 絕熱數 (Richardson Number, Ri) 的變種，與垂直結構相關。
   在準地轉理論中，一個重要的參數是準地轉參數 $\beta$ （注意這裡的 $\beta$ 不同於 Rossby 波的 $\beta$，而是代表地球自轉參數 $f$ 隨緯度變化的率 $df/dy$）：
   $\beta = \frac{df}{dy}$
   另一個相關的無量綱參數是斜壓參數 (Prandtl Number, Pr)，或者更常提及的是絕熱參數 (Brunt-Väisälä frequency, $N$) 和 大氣尺度參數 $L$ 的組合。
   在準地轉理論中，我們常常引入一個與垂直尺度和水平尺度相關的參數，例如：
   $\epsilon = \frac{U}{NL}$
   或者，在考慮斜壓性時，會用到斜壓性參數 (baroclinicity parameter)，例如：
   $\mu = \frac{L}{H}$
   其中 $H$ 是大氣的特徵垂直尺度（如比例高度）。

   在準地轉理論中，我們關注的是非地轉速度 $\mathbf{v}' = \mathbf{v} - \mathbf{v}_g$ 和 非地轉垂直速度 $w'$。準地轉理論的關鍵在於，雖然慣性力相對於科氏力很小，但非地轉速度可以通過渦度方程和熱力風關係來驅動天氣系統的演變，特別是通過垂直運動（$w$）來影響溫度的重新分佈，進而影響位勢的變化。

【準地轉近似的主要特性】

1. 近似地轉平衡：水平運動主要由壓力梯度力和科氏力平衡。

本題旨在考驗學生對大氣動力學中渦度概念的理解，以及如何推導描述渦度隨時間演變的渦度方程，並解釋其物理意義。渦度方程是分析天氣系統（如氣旋、反氣旋）發展和演變的關鍵工具。

【渦度方程的推導】
我們從水平運動方程開始。以 $x$ 方向為例：
$\frac{\partial u}{\partial t} + u \frac{\partial u}{\partial x} + v \frac{\partial u}{\partial y} + w \frac{\partial u}{\partial z} - f v = -\frac{1}{\rho} \frac{\partial p}{\partial x}$
以 $y$ 方向為例：
$\frac{\partial v}{\partial t} + u \frac{\partial v}{\partial x} + v \frac{\partial v}{\partial y} + w \frac{\partial v}{\partial z} + f u = -\frac{1}{\rho} \frac{\partial p}{\partial y}$
其中 $u, v, w$ 分別是 $x, y, z$ 方向的速度分量，$f$ 是科氏參數，$\rho$ 是密度，$p$ 是壓力。

我們感興趣的是垂直渦度 $\zeta = \frac{\partial v}{\partial x} - \frac{\partial u}{\partial y}$。
對 $x$ 方向方程關於 $y$ 求偏導，對 $y$ 方向方程關於 $x$ 求偏導，然後相減。

對 $x$ 方向方程關於 $y$ 求偏導：
$\frac{\partial}{\partial y} \left( \frac{\partial u}{\partial t} \right) + \frac{\partial}{\partial y} \left( u \frac{\partial u}{\partial x} \right) + \frac{\partial}{\partial y} \left( v \frac{\partial u}{\partial y} \right) + \frac{\partial}{\partial y} \left( w \frac{\partial u}{\partial z} \right) - f \frac{\partial v}{\partial y} = -\frac{\partial}{\partial y} \left( \frac{1}{\rho} \frac{\partial p}{\partial x} \right)$
$\frac{\partial^2 u}{\partial t \partial y} + u \frac{\partial^2 u}{\partial y \partial x} + \frac{\partial u}{\partial y} \frac{\partial u}{\partial y} + v \frac{\partial^2 u}{\partial y^2} + w \frac{\partial^2 u}{\partial y \partial z} + \frac{\partial w}{\partial y} \frac{\partial u}{\partial z} - f \frac{\partial v}{\partial y} = -\frac{1}{\rho} \frac{\partial^2 p}{\partial y \partial x} - \frac{\partial \rho^{-1}}{\partial y} \frac{\partial p}{\partial x}$

對 $y$ 方向方程關於 $x$ 求偏導：
$\frac{\partial}{\partial x} \left( \frac{\partial v}{\partial t} \right) + \frac{\partial}{\partial x} \left( u \frac{\partial v}{\partial x} \right) + \frac{\partial}{\partial x} \left( v \frac{\partial v}{\partial y} \right) + \frac{\partial}{\partial x} \left( w \frac{\partial v}{\partial z} \right) + f \frac{\partial u}{\partial x} = -\frac{\partial}{\partial x} \left( \frac{1}{\rho} \frac{\partial p}{\partial y} \right)$
$\frac{\partial^2 v}{\partial t \partial x} + u \frac{\partial^2 v}{\partial x^2} + \frac{\partial u}{\partial x} \frac{\partial v}{\partial x} + v \frac{\partial^2 v}{\partial x \partial y} + w \frac{\partial^2 v}{\partial x \partial z} + \frac{\partial w}{\partial x} \frac{\partial v}{\partial z} + f \frac{\partial u}{\partial x} = -\frac{1}{\rho} \frac{\partial^2 p}{\partial x \partial y} - \frac{\partial \rho^{-1}}{\partial x} \frac{\partial p}{\partial y}$

兩式相減，並利用 $\frac{\partial^2 u}{\partial t \partial y} = \frac{\partial^2 v}{\partial t \partial x}$ 和 $\frac{\partial^2 p}{\partial y \partial x} = \frac{\partial^2 p}{\partial x \partial y}$。
$\left( u \frac{\partial^2 u}{\partial y \partial x} + \frac{\partial u}{\partial y} \frac{\partial u}{\partial y} + v \frac{\partial^2 u}{\partial y^2} + w \frac{\partial^2 u}{\partial y \partial z} + \frac{\partial w}{\partial y} \frac{\partial u}{\partial z} - f \frac{\partial v}{\partial y} \right) - \left( u \frac{\partial^2 v}{\partial x^2} + \frac{\partial u}{\partial x} \frac{\partial v}{\partial x} + v \frac{\partial^2 v}{\partial x \partial y} + w \frac{\partial^2 v}{\partial x \partial z} + \frac{\partial w}{\partial x} \frac{\partial v}{\partial z} + f \frac{\partial u}{\partial x} \right)$
$= \left( -\frac{1}{\rho} \frac{\partial^2 p}{\partial y \partial x} - \frac{\partial \rho^{-1}}{\partial y} \frac{\partial p}{\partial x} \right) - \left( -\frac{1}{\rho} \frac{\partial^2 p}{\partial x \partial y} - \frac{\partial \rho^{-1}}{\partial x} \frac{\partial p}{\partial y} \right)$

整理左邊：
$u \left( \frac{\partial^2 u}{\partial y \partial x} - \frac{\partial^2 v}{\partial x^2} \right) + \frac{\partial u}{\partial y} \frac{\partial u}{\partial y} + v \left( \frac{\partial^2 u}{\partial y^2} - \frac{\partial^2 v}{\partial x \partial y} \right) + w \left( \frac{\partial^2 u}{\partial y \partial z} - \frac{\partial^2 v}{\partial x \partial z} \right) + \frac{\partial w}{\partial y} \frac{\partial u}{\partial z} - \frac{\partial w}{\partial x} \frac{\partial v}{\partial z} - f \left( \frac{\partial v}{\partial y} + \frac{\partial u}{\partial x} \right)$
注意到 $\frac{\partial v}{\partial y} + \frac{\partial u}{\partial x} = \nabla \cdot \mathbf{v}$，這是水平散度。
令 $\zeta = \frac{\partial v}{\partial x} - \frac{\partial u}{\partial y}$，則 $\frac{\partial \zeta}{\partial t} = \frac{\partial^2 v}{\partial t \partial x} - \frac{\partial^2 u}{\partial t \partial y}$。
我們需要將上面的式子重新組合。

一個更系統的方法是從矢量形式的運動方程出發，取其旋度。
運動方程（忽略黏滯和非保守力）：
$\frac{d \mathbf{v}}{d t} = - \frac{1}{\rho} \nabla p - g \mathbf{k}$
$\frac{\partial \mathbf{v}}{\partial t} + (\mathbf{v} \cdot \nabla) \mathbf{v} = - \frac{1}{\rho} \nabla p - g \mathbf{k}$
取旋度 $\nabla \times$:
$\nabla \times \frac{\partial \mathbf{v}}{\partial t} + \nabla \times [(\mathbf{v} \cdot \nabla) \mathbf{v}] = \nabla \times \left( - \frac{1}{\rho} \nabla p \right) - \nabla \times (g \mathbf{k})$
左邊：
$\frac{\partial}{\partial t} (\nabla \times \mathbf{v}) + \nabla \times [(\mathbf{v} \cdot \nabla) \mathbf{v}] = \frac{\partial \boldsymbol{\zeta}}{\partial t} + \nabla \times [(\mathbf{v} \cdot \nabla) \mathbf{v}]$
其中 $\boldsymbol{\zeta} = \nabla \times \mathbf{v}$ 是渦度矢量。
右邊：
$\nabla \times \left( - \frac{1}{\rho} \nabla p \right) = - \nabla \left( \frac{1}{\rho} \right) \times \nabla p - \frac{1}{\rho} (\nabla \times \nabla p)$
由於 $\nabla \times \nabla p = 0$，所以
$\nabla \times \left( - \frac{1}{\rho} \nabla p \right) = - \nabla \left( \frac{1}{\rho} \right) \times \nabla p = \frac{1}{\rho^2} \nabla \rho \times \nabla p$
這個項稱為散度渦度生成項 (divergence vorticity generation term) 或 非保守力渦度生成項 (non-conservative vorticity generation term)。如果密度和壓力場是同向或反向的，則此項為零。

$\nabla \times (g \mathbf{k}) = 0$。

所以，
$\frac{\partial \boldsymbol{\zeta}}{\partial t} + \nabla \times [(\mathbf{v} \cdot \nabla) \mathbf{v}] = \frac{1}{\rho^2} \nabla \rho \times \nabla p$

現在處理 $\nabla \times [(\mathbf{v} \cdot \nabla) \mathbf{v}]$。
展開 $(\mathbf{v} \cdot \nabla) \mathbf{v} = (u \frac{\partial}{\partial x} + v \frac{\partial}{\partial y} + w \frac{\partial}{\partial z}) (u \mathbf{i} + v \mathbf{j} + w \mathbf{k})$。
$\nabla \times [(\mathbf{v} \cdot \nabla) \mathbf{v}] = \nabla \times (\mathbf{v} \times (\nabla \times \mathbf{v}) + (\nabla \cdot \mathbf{v}) \mathbf{v})$。
這是一個複雜的矢量恆等式。
$\nabla \times [(\mathbf{v} \cdot \nabla) \mathbf{v}] = (\mathbf{v} \cdot \nabla) \boldsymbol{\zeta} - (\boldsymbol{\zeta} \cdot \nabla) \mathbf{v} + (\nabla \cdot \mathbf{v}) \boldsymbol{\zeta} - (\boldsymbol{\zeta} \cdot \nabla) \mathbf{v}$。
或者，使用 $\nabla \times (\mathbf{a} \times \mathbf{b}) = \mathbf{a} (\nabla \cdot \mathbf{b}) - \mathbf{b} (\nabla \cdot \mathbf{a}) + \mathbf{b} (\nabla \cdot \mathbf{a}) - \mathbf{a} (\nabla \cdot \mathbf{b})$。
$\nabla \times [(\mathbf{v} \cdot \nabla) \mathbf{v}] = \nabla \times (\mathbf{v} \times \boldsymbol{\omega})$，其中 $\boldsymbol{\omega}$ 是渦度。
恆等式：$\nabla \times (\mathbf{a} \times \mathbf{b}) = (\mathbf{b} \cdot \nabla) \mathbf{a} - (\mathbf{a} \cdot \nabla) \mathbf{b} + \mathbf{a}(\nabla \cdot \mathbf{b}) - \mathbf{b}(\nabla \cdot \mathbf{a})$。
令 $\mathbf{a} = \mathbf{v}, \mathbf{b} = \boldsymbol{\zeta}$。
$\nabla \times (\mathbf{v} \times \boldsymbol{\zeta}) = (\boldsymbol{\zeta} \cdot \nabla) \mathbf{v} - (\mathbf{v} \cdot \nabla) \boldsymbol{\zeta} + \mathbf{v}(\nabla \cdot \boldsymbol{\zeta}) - \boldsymbol{\zeta}(\nabla \cdot \mathbf{v})$。
由於 $\nabla \cdot \boldsymbol{\zeta} = 0$（無散度渦度），
$\nabla \times (\mathbf{v} \times \boldsymbol{\zeta}) = (\boldsymbol{\zeta} \cdot \nabla) \mathbf{v} - (\mathbf{v} \cdot \nabla) \boldsymbol{\zeta} - \boldsymbol{\zeta}(\nabla \cdot \mathbf{v})$。
所以 $\nabla \times [(\mathbf{v} \cdot \nabla) \mathbf{v}] = -[(\boldsymbol{\zeta} \cdot \nabla) \mathbf{v} - (\mathbf{v} \cdot \nabla) \boldsymbol{\zeta} - \boldsymbol{\zeta}(\nabla \cdot \mathbf{v})]$。

回到垂直渦度 $\zeta = \frac{\partial v}{\partial x} - \frac{\partial u}{\partial y}$ 的時間變化率。
從上面計算的兩式相減，並考慮左邊的對稱性：
$\frac{\partial \zeta}{\partial t} + \dots = \dots$
一個更簡潔的推導是利用連續性方程和流體靜力平衡。
考慮絕熱、無摩擦、不可壓縮流體。
運動方程：
$\frac{du}{dt} - fv = -\frac{1}{\rho} \frac{\partial p}{\partial x}$
$\frac{dv}{dt} + fu = -\frac{1}{\rho} \frac{\partial p}{\partial y}$
連續性方程：
$\frac{d\rho}{dt} + \rho \nabla \cdot \mathbf{v} = 0$
對於不可壓縮流體，$\rho = \text{const}$，所以 $\nabla \cdot \mathbf{v} = \frac{\partial u}{\partial x} + \frac{\partial v}{\partial y} + \frac{\partial w}{\partial z} = 0$。
若考慮可壓縮流體，則 $\nabla \cdot \mathbf{v} = -\frac{1}{\rho} \frac{d\rho}{dt}$。

我們對運動方程取 $z$ 方向的旋度。
$\frac{\partial}{\partial t} (\nabla \times \mathbf{v}) + \nabla \times [(\mathbf{v} \cdot \nabla) \mathbf{v}] = \nabla \times \left( - \frac{1}{\rho} \nabla p \right) - \nabla \times (f \mathbf{k} \times \mathbf{v})$
$\frac{\partial \zeta}{\partial t} + \nabla \times [(\mathbf{v} \cdot \nabla) \mathbf{v}] = \frac{1}{\rho^2} \nabla \rho \times \nabla p - f \nabla \times (\mathbf{k} \times \mathbf{v})$
其中 $\mathbf{k} \times \mathbf{v} = -v \mathbf{i} + u \mathbf{j}$。
$\nabla \times (\mathbf{k} \times \mathbf{v}) = \nabla \times (-v \mathbf{i} + u \mathbf{j}) = \frac{\partial u}{\partial x} - (-\frac{\partial v}{\partial y}) = \frac{\partial u}{\partial x} + \frac{\partial v}{\partial y} = \nabla_h \cdot \mathbf{v}$，其中 $\nabla_h$ 是水平散度算子。

展開 $\nabla \times [(\mathbf{v} \cdot \nabla) \mathbf{v}]$：
$\nabla \times [(\mathbf{v} \cdot \nabla) \mathbf{v}] = \nabla \times [ u \frac{\partial \mathbf{v}}{\partial x} + v \frac{\partial \mathbf{v}}{\partial y} + w \frac{\partial \mathbf{v}}{\partial z} ]$
$= \frac{\partial}{\partial x} (u \mathbf{v}) \times \mathbf{i} + \dots$
這部分展開比較繁瑣。
一個重要的恆等式是 $\nabla \times [(\mathbf{v} \cdot \nabla) \mathbf{v}] = (\mathbf{v} \cdot \nabla) \boldsymbol{\zeta} - (\boldsymbol{\zeta} \cdot \nabla) \mathbf{v} + \boldsymbol{\zeta}(\nabla \cdot \mathbf{v}) - \mathbf{v}(\nabla \cdot \boldsymbol{\zeta})$。
由於 $\nabla \cdot \boldsymbol{\zeta} = 0$，且對於不可壓縮流體 $\nabla \cdot \mathbf{v} = 0$。
所以 $\nabla \times [(\mathbf{v} \cdot \nabla) \mathbf{v}] = (\mathbf{v} \cdot \nabla) \boldsymbol{\zeta} - (\boldsymbol{\zeta} \cdot \nabla) \mathbf{v}$。

因此，
$\frac{\partial \boldsymbol{\zeta}}{\partial t} + (\mathbf{v} \cdot \nabla) \boldsymbol{\zeta} - (\boldsymbol{\zeta} \cdot \nabla) \mathbf{v} = \frac{1}{\rho^2} \nabla \rho \times \nabla p - f (\nabla_h \cdot \mathbf{v})$
$\frac{d \boldsymbol{\zeta}}{d t} - (\boldsymbol{\zeta} \cdot \nabla) \mathbf{v} = \frac{1}{\rho^2} \nabla \rho \times \nabla p - f (\nabla_h \cdot \mathbf{v})$

我們感興趣的是垂直渦度 $\zeta = \zeta_z = \frac{\partial v}{\partial x} - \frac{\partial u}{\partial y}$。
所以我們需要考慮 $\frac{d \zeta_z}{d t}$。
$\frac{d \zeta_z}{d t} = \frac{\partial \zeta_z}{\partial t} + (\mathbf{v} \cdot \nabla) \zeta_z$

本題考查的是行星邊界層（Planetary Boundary Layer, PBL）中的幾個重要概念：Ekman spiral、Ekman pumping 和 spin-down effect。這些概念對於理解地表與大氣之間的能量和動量交換，以及邊界層內風場的結構至關重要。

【Ekman Spiral (艾克曼螺旋)】
概念： Ekman spiral 描述了在考慮了科氏力、壓力梯度力（通常假設為地衡風）以及摩擦力的情況下，風速矢量隨高度在行星邊界層內如何變化，形成一個螺旋狀的結構。

推導與原理：
考慮邊界層內某一層次的流體質點，其運動受到以下力的作用：

1. 壓力梯度力 (Pressure Gradient Force, PGF)：通常假設為提供一個恆定的地衡風 $\mathbf{v}_g = (u_g, v_g)$。
2. 科氏力 (Coriolis Force)：$-f \mathbf{k} \times \mathbf{v}$，其中 $\mathbf{v} = (u, v)$ 是該高度的風速。
3. 摩擦力 (Friction Force)：通常假設為與風速成正比（分子摩擦）或與風速梯度成正比。最簡單的分子摩擦形式是 $F_f = -k \mathbf{v}$，其中 $k$ 是一個摩擦係數。更真實的邊界層模型會考慮湍流交換，其有效摩擦係數會隨高度變化。

在穩態、水平無輻散且摩擦力為 $-k\mathbf{v}$ 的情況下，運動方程為：
$\frac{d \mathbf{v}}{d t} = 0$
$-f \mathbf{k} \times \mathbf{v} - k \mathbf{v} = - \frac{1}{\rho} \nabla p$
假設 $\nabla p$ 提供一個地衡風 $\mathbf{v}_g$。
$-f \mathbf{k} \times \mathbf{v} - k \mathbf{v} = - f \mathbf{k} \times \mathbf{v}_g$
寫成速度分量：
$-f v - k u = -f v_g$
$f u - k v = f u_g$
解這個聯立方程組，得到：
$u = u_g - \frac{f v_g + k u_g}{k^2 + f^2} k = u_g - \frac{k}{k^2 + f^2} (k u_g + f v_g)$
$v = v_g + \frac{f u_g - k v_g}{k^2 + f^2} k = v_g + \frac{k}{k^2 + f^2} (k v_g - f u_g)$
這個形式比較複雜。

更常見的 Ekman 層模型：
假設摩擦力是與渦度方程中的摩擦項相關，且在邊界層內，摩擦力可以近似為一個與風速成正比的力 $F_f = K \nabla^2 \mathbf{v}$，其中 $K$ 是渦黏係數。
在邊界層內，我們假設地衡風 $\mathbf{v}_g$ 是常數，但摩擦力隨高度變化。
運動方程（穩態，忽略水平平流）：
$-f \mathbf{k} \times \mathbf{v} + \frac{\partial}{\partial z} \left( K \frac{\partial \mathbf{v}}{\partial z} \right) = - \frac{1}{\rho} \nabla p$
假設 $\nabla p$ 提供一個恆定的地衡風 $\mathbf{v}_g$。
$-f \mathbf{k} \times \mathbf{v} + \frac{\partial}{\partial z} \left( K \frac{\partial \mathbf{v}}{\partial z} \right) = - f \mathbf{k} \times \mathbf{v}_g$
令 $\mathbf{v} = \mathbf{v}_g + \mathbf{v}'$，其中 $\mathbf{v}'$ 是非地轉風。
$-f \mathbf{k} \times (\mathbf{v}_g + \mathbf{v}') + \frac{\partial}{\partial z} \left( K \frac{\partial (\mathbf{v}_g + \mathbf{v}')}{\partial z} \right) = - f \mathbf{k} \times \mathbf{v}_g$
$-f \mathbf{k} \times \mathbf{v}' + \frac{\partial}{\partial z} \left( K \frac{\partial \mathbf{v}'}{\partial z} \right) = 0$
因為 $\mathbf{v}_g$ 是常數，所以 $\frac{\partial \mathbf{v}_g}{\partial z} = 0$。
這是一個關於 $\mathbf{v}'$ 的微分方程。

Ekman 螺旋的特徵：

• 地表風速：在地表附近，摩擦力最大，風速最小，且風向與地衡風方向成一個角度（在北半球，地表風向偏離地衡風向右約 15-30 度）。
• 隨高度變化：隨著高度增加，摩擦力減小，科氏力變得相對更重要。風速逐漸趨近於地衡風速，同時風向也逐漸向地衡風方向旋轉。
• 螺旋結構：將不同高度的風速矢量連接起來，就形成了一個向右（北半球）或向左（南半球）傾斜的螺旋結構。
• Ekman 層厚度：通常定義為摩擦力效應顯著的高度範圍，大約為 $H_{Ekman} \sim \sqrt{K/f}$，其中 $K$ 是渦黏係數。

物理意義：
Ekman spiral 表明，由於地表摩擦力的作用，邊界層內的風速和風向會隨高度發生顯著變化，並且與自由大氣中的地衡風不同。這種變化是理解風場結構、輻合輻散以及與海洋相互作用的基礎。

【Ekman Pumping (艾克曼抽吸/泵動)】
概念： Ekman pumping 是指在行星邊界層頂部，由於 Ekman spiral 的存在，導致的垂直運動（向上或向下）。它描述了由邊界層內的水平輻合或輻散引起的，通過 Ekman 層將動量傳遞到自由大氣的過程。

原理：
在 Ekman 層內，風速矢量隨高度旋轉。如果邊界層頂部的風速（通常接近地衡風）與邊界層底部的平均風速（受到摩擦力顯著影響）之間存在水平輻合或輻散，就會在 Ekman 層內產生淨的垂直運動。

考慮一個水平尺度 $L$ 的區域。如果該區域存在水平輻散，例如一個反氣旋的中心，則邊界層底部的平均風場會向外輻散。為了補償這種輻散，必須有從自由大氣向下的垂直運動，將空氣「泵」入邊界層。反之，在一個氣旋的中心，邊界層底部存在水平輻合，則必須有向上的垂直運動，將空氣「泵」出邊界層。

更嚴格地說，Ekman pumping 描述了邊界層頂部的垂直速度 $w_{top}$ 與邊界層頂部風場的水平散度 $\nabla_h \cdot \mathbf{v}_{top}$ 之間的關係。
在 Ekman 層內，平均風速 $\overline{\mathbf{v}}$ 和邊界層頂部風速 $\mathbf{v}_{top}$ （通常近似為地衡風 $\mathbf{v}_g$）之間存在差異。
如果我們對 Ekman 層內的運動方程進行積分，並考慮邊界條件（地表無滑移，頂部連續性），可以推導出：

核心觀念

本題考查斜壓不穩度的兩個重點：

1. 能量轉換：基本流的平均可用位能如何轉換為擾動可用位能，再轉換為擾動動能。
2. 不穩定波長：為何過短與過長的波動增長率較小，以及何種尺度最容易發展。

題目中的「$0'$」依大氣動力學慣例解讀為位溫擾動 $\theta'$，第二個熱通量則解讀為南北向擾動速度與位溫的相關項 $\overline{v'\theta'}$。

---

(a) 能量轉換與不穩定波長

1. 斜壓基本流的結構

斜壓大氣具有水平溫度梯度，因此等壓面上的位溫不均勻：

$$\frac{\partial \overline{\theta}}{\partial y}<0$$

表示北半球中緯度地區通常向北溫度降低。

由熱成風關係，若基本流為西風 $\overline{u}$，則：

$$\frac{\partial \overline{u}}{\partial z} = -\frac{g}{f\theta_0} \frac{\partial \overline{\theta}}{\partial y}$$

在北半球 $f>0$ 且 $\partial\overline{\theta}/\partial y<0$，因此：

$$\frac{\partial \overline{u}}{\partial z}>0$$

亦即西風隨高度增強，形成垂直風切。

這種水平溫度梯度與垂直風切所儲存的能量，是斜壓不穩度的主要能量來源。

---

2. 平均可用位能轉換為擾動可用位能

斜壓波中的暖空氣向極地方向輸送，冷空氣向赤道方向輸送，因此北半球通常有：

$$\overline{v'\theta'}>0$$

其中：

• $v'>0$：向北運動；
• $\theta'>0$：暖位溫擾動；
• $v'\theta'>0$：暖空氣向北、冷空氣向南。

此熱量輸送會減弱平均南北溫度梯度，使平均溫度場的可用位能降低，並將能量供給擾動。可概念性表示為：

$$\text{平均可用位能} \longrightarrow \text{擾動可用位能}$$

這是斜壓不穩度能夠成長的第一階段。

---

3. 擾動可用位能轉換為擾動動能

擾動中的暖空氣上升、冷空氣下沉時，通常有：

$$\overline{w'\theta'}>0$$

其中：

• $w'>0$ 且 $\theta'>0$：暖空氣上升；
• $w'<0$ 且 $\theta'<0$：冷空氣下沉。

垂直浮力通量為：

$$\overline{w'b'} = \frac{g}{\theta_0}\overline{w'\theta'}$$

因此 $\overline{w'\theta'}>0$ 表示浮力對垂直運動作正功，使擾動動能增加：

$$\text{擾動可用位能} \longrightarrow \text{擾動動能}$$

整體能量鏈為：

$$\text{平均可用位能} \longrightarrow \text{擾動可用位能} \longrightarrow \text{擾動動能}$$

擾動動能增加後，波幅便會增長，形成斜壓不穩度。

---

4. 不穩定波長的特性

斜壓波並非所有波長都有相同的增長率。

波長過短

波長很短時，水平溫度擾動的尺度小，垂直結構受到穩定層結的強烈限制。垂直運動需要克服浮力恢復力，且擾動容易受到三維結構與穩定度的抑制，因此增長率降低。

波長過長

波長很長時，雖然擾動可以大尺度輸送熱量，但水平梯度與波數太小，波動的旋轉與垂直耦合效率降低，因此增長時間很長，增長率也趨近於零。

中等尺度

介於兩者之間的波長最容易將：

• 水平溫度梯度；
• 垂直風切；
• 水平輻合與輻散；
• 垂直運動；

有效結合，因此具有最大增長率。實際中緯度天氣系統的斜壓波多為數千公里尺度，正是最有利於不穩定發展的範圍。

以經典 Eady 模型為例，最大增長率約為：

$$\sigma_{\max} \approx 0.31 \frac{f}{N} \left| \frac{\partial \overline{u}}{\partial z} \right|$$

最容易成長的波數約滿足：

$$k\frac{NH}{f}\approx 1.6$$

因此最適波長約為：

$$\lambda_{\max} \approx \frac{2\pi NH}{1.6f} \approx 3.9\frac{NH}{f}$$

這說明斜壓不穩度具有一個最有利的中等波長，而非波長越長或越短越容易發展。

---

(b) $y-z$ 剖面、傾斜等位溫面與運動向量

1. 等位溫面的傾斜

北半球中緯度地區有：

$$\frac{\partial\overline{\theta}}{\partial y}<0, \qquad \frac{\partial\overline{\theta}}{\partial z}>0$$

等位溫面滿足 $d\theta=0$，所以：

$$\frac{d z}{d y} = - \frac{\partial\overline{\theta}/\partial y} {\partial\overline{\theta}/\partial z} >0$$

因此在 $y-z$ 剖面上，等位溫面向北方逐漸升高：

本題考查學生對條件不穩定性 (Conditional Instability) 的理解，包括其定義、數學條件以及與大氣邊界層空氣性質的關聯。

(a) 條件不穩度、$\theta_e$ 和 $\theta$ 的垂直分布圖

核心概念：

• 假定溫度直減率 (Adiabatic Lapse Rate, ALR)：
  • 乾絕熱直減率 ($\Gamma_d$)：空氣塊在未飽和狀態下上升時，每升高 1 公里，溫度下降約 $9.8 \ ^{\circ}\text{C/km}$。
  • 濕絕熱直減率 ($\Gamma_s$)：空氣塊在飽和狀態下上升時，由於水汽凝結釋放潛熱，其溫度下降率小於 $\Gamma_d$，通常約為 $4-7 \ ^{\circ}\text{C/km}$，且隨溫度升高而減小。
• 位溫 ($\theta$)：表示空氣塊絕熱膨脹或壓縮到標準參考壓力（通常為 1000 hPa）時的溫度。它是一個衡量大氣穩定性的保守量，在絕熱過程中保持不變。
• 濕位溫 ($\theta_e$)：表示空氣塊絕熱地移動到 1000 hPa，並且其中所有水汽都凝結成液態水並被移除後，其溫度。它是一個衡量大氣中總潛熱的保守量，在絕熱和飽和過程中保持不變。

垂直分布圖的繪製與解釋：

我們通常在大氣垂直探測圖（如Skew-T log-p 圖）上分析這些參數。這裡我們繪製一個簡化的垂直剖面圖，橫軸為溫度，縱軸為高度（或壓力）。

1. $\theta$ (位溫) 的垂直分布：

   • 在穩定大氣中，位溫 $\theta$ 隨高度 增加。這意味著上層空氣的位溫高於下層空氣，所以空氣塊上升時會遇到位溫更高的環境，受到向下的浮力，系統穩定。
   • 在 中性大氣 中，位溫隨高度 不變。
   • 在 不穩定大氣 中，位溫隨高度 減小。

2. $\theta_e$ (濕位溫) 的垂直分布：

   • $\theta_e$ 的行為類似於 $\theta$，但它包含了水汽的潛熱效應。
   • 在 穩定大氣 中，濕位溫 $\theta_e$ 隨高度 增加。
   • 在 中性大氣 中，濕位溫隨高度 不變。
   • 在 不穩定大氣 中，濕位溫隨高度 減小。

3. 條件不穩定性 (Conditional Instability)：
   條件不穩定性是指大氣的穩定性取決於空氣塊是否飽和。

   • 未飽和狀態下：大氣是 穩定的。這意味著未飽和的空氣塊上升時，會比周圍環境冷（其溫度下降率為 $\Gamma_d$），受到向下的浮力。
   • 飽和狀態下：大氣是 不穩定 的。這意味著飽和的空氣塊上升時，會比周圍環境暖（其溫度下降率為 $\Gamma_s < \Gamma_d$），受到向上的浮力。

   如何體現在圖上？
   在垂直分布圖上，條件不穩定性對應於：

   • 未飽和空氣塊的穩定性：未飽和空氣塊的溫度（沿乾絕熱線 $\theta$）的垂直梯度比環境溫度（或環境 $\theta$）的垂直梯度要小。也就是說，未飽和空氣塊上升時，其溫度下降率 $\Gamma_d$ 小於環境溫度下降率。
   • 飽和空氣塊的不穩定性：飽和空氣塊的溫度（沿濕絕熱線 $\theta_e$）的垂直梯度比環境溫度（或環境 $\theta$）的垂直梯度要大。也就是說，飽和空氣塊上升時，其溫度下降率 $\Gamma_s$ 大於環境溫度下降率。

   繪製圖示：
   在高度 $z$ 處，繪製：

   • 環境溫度廓線：通常表示為一條隨高度遞減的曲線。
   • 位溫 ($\theta$) 廓線：在條件不穩定大氣中，未飽和層的 $\theta$ 廓線應該隨高度 增加（表示穩定）。
   • 濕位溫 ($\theta_e$) 廓線：在條件不穩定大氣中，飽和層的 $\theta_e$ 廓線應該隨高度 減小（表示不穩定）。
   • 注意：通常 $\theta_e \ge \theta$。在未飽和層，空氣塊沿著 $\theta$ 線移動；在飽和層，空氣塊沿著 $\theta_e$ 線移動。

   簡化圖示說明：
   假設在某高度 $z_1$ 以下，空氣是未飽和的；在 $z_1$ 以上，空氣是飽和的。

   • 在 $z_1$ 以下，環境溫度廓線和 $\theta$ 廓線都隨高度增加，且 $\theta$ 廓線的斜率（表示環境的穩定性）應大於乾絕熱直減率 $\Gamma_d$（即環境 $\theta$ 隨高度增加的速率大於 $\Gamma_d$）。
   • 在 $z_1$ 以上，環境溫度廓線和 $\theta_e$ 廓線都隨高度減小，但 $\theta_e$ 廓線的斜率（表示環境的濕穩定性）應小於濕絕熱直減率 $\Gamma_s$（即環境 $\theta_e$ 隨高度減小的速率小於 $\Gamma_s$）。

   總結圖示特徵：

   • 未飽和層：環境溫度廓線比乾絕熱線 $\Gamma_d$ 更陡峭（$\theta$ 隨高度增加）。
   • 飽和層：環境溫度廓線比濕絕熱線 $\Gamma_s$ 更平緩（$\theta_e$ 隨高度減小）。
   • 這意味著，未飽和空氣塊上升時，其溫度比環境溫度下降得快，所以穩定；飽和空氣塊上升時，其溫度比環境溫度下降得慢，所以不穩定。

(b) 條件不穩度數學條件與邊界層空氣性質

數學條件：
大氣處於條件不穩定狀態，是指：

1. 未飽和層的穩定性：
   $\frac{\partial \theta}{\partial z} > 0$
   或者，用溫度和高度表示，環境溫度直減率 $\Gamma_{env} < \Gamma_d$。

本題旨在考驗學生對大尺度環流、熱力學過程（加熱與冷卻）以及動力學過程（斜壓不穩定性）之間相互作用的理解，特別是它們如何共同塑造大氣的三胞環流和維持熱力風平衡。

核心概念：

1. 三胞環流 (Hadley Cell, Ferrel Cell, Polar Cell)：這是地球大尺度環流的基本模型，描述了赤道與極區之間能量傳輸的過程。
   • 哈德里環流 (Hadley Cell)：在熱帶地區，赤道附近空氣受熱上升，向高緯度輸送，在高空冷卻下沉，然後在地表向赤道輻合。
   • 費雷環流 (Ferrel Cell)：在中緯度地區，由哈德里環流和極地環流之間相互作用產生，是一個次級環流，空氣在地表向極區輻合，在高空向赤道輻合。
   • 極地環流 (Polar Cell)：在極區，極地空氣冷卻下沉，向赤道輻合，然後在高空向極區輸送。
2. 熱力風平衡 (Thermal Wind Balance)：描述了溫度梯度如何影響風的垂直切變。在水平溫度梯度存在的情況下，風速隨高度的變化與溫度梯度方向和大小有關。
3. 斜壓不穩定性 (Baroclinic Instability)：是中緯度地區天氣系統（氣旋、反氣旋）發展的主要機制，源於溫度梯度和垂直風切變的相互作用，將位能轉化為動能。
4. 熱帶加熱與極區冷卻：這是地球大氣能量收支的主要驅動力，赤道地區接收的太陽輻射多於散失的輻射，而極區則相反。這種能量差異驅動了大氣的環流。

討論：

1. 熱帶加熱與極區冷卻驅動大尺度環流：

   • 赤道地區的強烈加熱：使得赤道空氣變暖，密度減小，在其上方產生低壓，空氣上升。上升的暖濕空氣在高空向兩極輸送。
   • 極區的強烈冷卻：使得極地空氣變冷，密度增大，在其下方產生高壓，空氣下沉。下沉的冷空氣在地表向赤道輸送。
   • 這種由東西向溫度差異（赤道暖、極區冷）驅動的 南北向的溫度梯度 是產生大尺度環流（尤其是哈德里環流）和維持熱力風平衡的根本原因。

2. 三胞環流與熱力風平衡的關係：

   • 哈德里環流：在熱帶地區，雖然溫度梯度較小，但上升和下沉運動也產生一定的環流。其高空向兩極的輸送和低空向赤道的輸送，構成了一個相對 準地轉 的環流。
   • 極地環流：極區的冷卻導致空氣下沉，並在地表向低緯度輻合。在高空，冷空氣向極區輸送。
   • 費雷環流（中緯度）：這是 斜壓不穩定性最活躍 的區域。
     • 在中緯度地區，存在顯著的 從南向北的溫度梯度（南方暖、北方冷）。
     • 這個溫度梯度與地轉風的關係由 熱力風方程式 描述：$\frac{\partial u_g}{\partial z} = -\frac{g}{f T} \frac{\partial T}{\partial y}$。在中緯度，$\frac{\partial T}{\partial y} < 0$（北方冷），$f>0$（北半球），因此 $\frac{\partial u_g}{\partial z} > 0$，表示 西風隨高度增加。
     • 這種隨高度增加的西風（即存在垂直風切變）是斜壓不穩定性的必要條件。

本題考查學生對渦度方程、輻合/輻散與風場變化的關係，特別是關於西風的變化。

(a) 輻合增加西風的解釋

核心概念：

• 輻合 (Convergence)：指空氣在水平方向向某一點或區域匯聚。
• 渦度 (Vorticity)：描述流體旋轉的程度，垂直渦度 $\zeta = \frac{\partial v}{\partial x} - \frac{\partial u}{\partial y}$。
• 西風 (Zonal Wind)：即 $u$ 分量的風，在東亞地區通常指從西向東的風。
• 天氣示意圖 $uv > 0$：這表示 $u$ 和 $v$ 同號。在北半球，如果 $u>0$（西風）且 $v>0$（南風），則 $uv>0$。如果 $u<0$（東風）且 $v<0$（北風），則 $uv>0$。題目中提到「西風」，通常指 $u>0$，所以 $uv>0$ 可能意味著 $u>0$ 和 $v>0$（西風和南風的組合）。然而，題目中的 $uv>0$ 條件本身並不直接指向西風的增加，更可能是在討論某種特定的運動模式。我們將重點放在「輻合可以增加西風」這一論述上。

解釋：
要解釋「單純的輻合可以增加西風」，我們需要藉助渦度方程。
渦度方程（簡化形式，忽略摩擦和傾斜項）：
$\frac{\partial \zeta}{\partial t} = - \zeta (\nabla \cdot \mathbf{v}) + \frac{1}{\rho^2} (\nabla \rho \times \nabla p)_z - f (\nabla_h \cdot \mathbf{v})$
其中 $\zeta$ 是垂直渦度，$\nabla \cdot \mathbf{v}$ 是散度，$\nabla_h \cdot \mathbf{v} = \frac{\partial u}{\partial x} + \frac{\partial v}{\partial y}$ 是水平散度。
輻合意味著 $\nabla_h \cdot \mathbf{v} < 0$。

如何輻合增加西風 ($u$)？
這通常與 斜壓性 和 科氏力效應 有關。

1. 科氏力效應 (Coriolis Force Effect)：
   在北半球 ($f>0$)，輻合 ($\nabla_h \cdot \mathbf{v} < 0$) 會導致渦度增加（根據渦度方程的 $-f (\nabla_h \cdot \mathbf{v})$ 項）。
   如果一個區域原本是西風 ($u>0$)，並且存在南風 ($v>0$)，則 $uv>0$。
   輻合會增加該區域的渦度。如果該渦度增加是 氣旋式 的（即正渦度增加），這通常與 低壓系統的發展 有關。
   低壓系統的發展往往伴隨著 暖空氣在西側（高空）和東側（低空）的輸送。
   在中緯度地區，西風隨高度增加（熱力風關係）。一個發展中的低壓系統，其壓力梯度力會改變原有的地衡風平衡。

2. 斜壓性與溫度梯度：
   在中緯度地區，通常存在從南向北的溫度梯度（南方暖、北方冷）。根據熱力風關係：
   $\frac{\partial u_g}{\partial z} = -\frac{g}{f T} \frac{\partial T}{\partial y}$
   如果 $\frac{\partial T}{\partial y} < 0$ (北方冷)，則 $\frac{\partial u_g}{\partial z} > 0$，即西風隨高度增加。

   輻合如何影響西風？

   • 情況一：暖區輻合。如果輻合發生在一個暖區，並且伴隨著上升運動，這有利於氣旋的發展。氣旋的發展通常會加劇 西風的垂直切變，使得高空的西風增強。
   • 情況二：冷區輻合。如果輻合發生在冷區，可能與冷鋒後的西北氣流有關。
   • 更普遍的解釋：輻合本身並不直接「增加」西風。輻合是 低壓系統發展 的一個標誌。低壓系統的發展（尤其是氣旋）在中緯度地區，通常伴隨著 暖空氣被輸送到高空，冷空氣被輸送到低空。由於中緯度西風隨高度增加，暖空氣在高空被輸送，這就意味著高空的西風會增強。

   從流體動力學角度：
   考慮一個簡單的例子：一個純粹的輻合區域，沒有溫度梯度，只有摩擦力。這種情況下，輻合會逐漸減慢水平速度。
   但是，在大氣中，輻合通常發生在 低壓系統 中，而低壓系統的發展又與 斜壓不穩定性 和 溫度梯度 緊密相關。
   關鍵在於：輻合通常是低壓系統發展的結果，而低壓系統的發展（特別是氣旋）在中緯度地區，會加強高空的西風。

本題考查學生對行星尺度波動（Rossby waves 和 Kelvin waves）的理解，包括其產生機制、傳播方式以及能量傳輸特性。

(a) Rossby Waves (羅斯貝波) 的傳送機制與能量頻散關係

核心概念：

• Rossby 波：是一種行星尺度（長波長）的地球大氣或海洋中的波動，其產生和維持主要與地球自轉參數 $f$ 隨緯度變化的率 $\beta = df/dy$ 有關。

• 傳送機制：Rossby 波的傳播是通過 絕對渦度守恆 的概念來理解的。當一個空氣塊（或水塊）在南北方向移動時，由於 $f$ 的變化，其絕對渦度會發生改變。為了維持絕對渦度守恆（在忽略其他效應的情況下），空氣塊必須通過水平運動來補償這種變化。

  • 在北半球，緯度越高，$f$ 越大。
  • 如果一個空氣塊從低緯度（$f$ 小）向高緯度（$f$ 大）移動，為了保持絕對渦度 $\zeta_a = \zeta + f$ 不變，其相對渦度 $\zeta$ 必須減小。這通常通過 向西的運動 來實現（例如，氣旋的相對渦度減小）。
  • 反之，如果空氣塊從高緯度向低緯度移動，其相對渦度必須增大，通常通過 向東的運動。
  • 傳播方向：這種通過相對渦度調整來適應 $f$ 變化的機制，使得 Rossby 波的 總體傳播方向是向西的（相對於地球參考系）。這也是為什麼我們在天氣圖上看到的中緯度高空槽脊結構，其移動速度相對較慢，並且整體上有向西傳播的趨勢。

• 能量頻散 (Energy Dispersion)：

  • 頻散 (Dispersion)：指不同波長的波在介質中傳播速度不同的現象。
  • Rossby 波是頻散的：這意味著，不同波長（或波数 $k$）的 Rossby 波具有不同的 相位速度 ($c_p$) 和 群速度 ($c_g$)。
  • 相位速度 ($c_p$)：波形狀傳播的速度。對於 Rossby 波，其相位速度為：
    $c_p = U - \frac{\beta}{k^2}$
    其中 $U$ 是平均緯向風速，$k$ 是水平波数，$k^2 = k_x^2 + k_y^2$。
    這表明，長波長（小 $k$）的 Rossby 波比短波長（大 $k$）的 Rossby 波傳播得更快。
  • 群速度 ($c_g$)：能量和信息傳播的速度。對於 Rossby 波，群速度通常也向西傳播，但其速度與相位速度不同，且與波数有關。
  • 能量頻散的影響：由於不同波長的 Rossby 波以不同的速度傳播，一個由多種波長組成的初始擾動，會隨著時間推移而散開。這使得 Rossby 波的結構和傳播行為變得複雜。
  • 向西傳播的極限：當平均緯向風速 $U$ 減小或變為東風時，相位速度 $c_p$ 可能變為負值，表示波向西傳播。當 $U - \frac{\beta}{k^2} < 0$ 時，波向西傳播。如果 $U > \frac{\beta}{k^2}$，則波向東傳播。
  • 頻散關係 (Dispersion Relation)：$c_p$ 與 $k$ 的關係稱為頻散關係。對於 Rossby 波，頻散關係表明，長波長（小 $k$）的波比短波長（大 $k$）的波具有更大的向西（或向東）傳播速度。

總結：
Rossby 波的傳播機制主要是通過空氣塊適應緯度變化下的科氏參數變化，導致相對渦度的調整，其總體傳播方向是向西的。它們是頻散的，不同波長（波数）的 Rossby 波具有不同的傳播速度，長波長波的傳播速度通常比短波長波更快。

(b) Kelvin Waves (開爾文波) 的傳送機制與傳送方向