文章目录
大气和天气系统的数学推导
本课程旨在通过物理原理和数学模型，解释全球重要天气系统（如急流、气旋和大尺度波动）的空间结构和演变。这些天气系统对气候变化、极端天气事件（如强降水、热浪和干旱）有重要影响。课程强调理解和应用理论以解决问题，并使用数学方程作为基础工具。此外，还包括实际实验环节，以帮助学员将理论应用于真实案例。课程也涉及数值模型的近似解，这些解需要结合物理理论来简化系统，从而预测气候和天气模式的演变。
在这一篇博客中，我们将深入探讨与大气和天气系统相关的数学公式推导。这些公式基于流体力学、热力学和大气动力学原理，并通过具体的推导过程，帮助我们更好地理解大气中的旋转运动和地转平衡。
1.1 Thermaldynamics and vertical structure
1.1.1 Ideal gas law
我们假设大气可以近似为理想气体，单位体积理想气体定律为：
p = ρ R T \colorbox{yellow}{$p = \rho RT$} p=ρRT
其中：( p ) 是气压， ( ρ \rho ρ) 是密度，( R ) 是理想气体常数， ( T ) 是温度。
1.1.2 hydrostatic balance with gravitational acceleration
在大气中，垂直压力梯度力和重力达到平衡，称为 静力学平衡。平衡方程为：
1
ρ
∂
p
∂
z
=
−
g
\colorbox{yellow}{$\frac{1}{\rho} \frac{\partial p}{\partial z} = -g$}
ρ1∂z∂p=−g
这里，( z ) 表示高度，( g ) 为重力加速度。
为了求解压力随高度的变化，我们假设温度随高度恒定不变，此时可以推导出压力的指数衰减关系：
p = p 0 e − z / H \colorbox{yellow}{$p = p_0 e^{-z/H}$} p=p0e−z/H
其中 ( H = R T g H = \frac{RT}{g} H=gRT ) 是 压强尺度高度（pressure scale-height）。
1.1.3 位温与布伦特-瓦依萨拉频率（Brunt-Väisälä frequency)
位温（Potential Temperature）（
θ
\theta
θ）用于描述大气中空气团的热力学状态。它表示一个空气团在绝热过程中如果被移动到某一标准压力（通常是1000 hPa）时的温度。假设大气中的空气团是绝热上升或下降的：
θ
=
T
(
p
0
p
)
R
/
c
p
\colorbox{yellow}{$\theta = T \left(\frac{p_0}{p}\right)^{R/c_p}$}
θ=T(pp0)R/cp
在稳定大气中，绝热温度递减率(adiabatic lapse rate) 可以定义为：
∂
T
∂
z
=
−
g
c
p
=
−
Γ
a
\colorbox{yellow}{$\frac{\partial T}{\partial z} = -\frac{g}{c_p}=-\Gamma_a$}
∂z∂T=−cpg=−Γa
- 干绝热 − 10 K k m − 1 -10K km^{-1} −10Kkm−1
- stably stratified: Γ < Γ a \Gamma < \Gamma_a Γ<Γa; Γ \Gamma Γ是实际的温度梯度
当空气团在这种稳定的环境中被扰动时，它会围绕其原始高度来回振荡，这种振荡的频率称为布伦特-魏萨拉频率（Brunt-Väisälä frequency） ( N ) ，表示为：
N
2
=
g
T
(
Γ
a
−
Γ
)
N^2 = \frac{g}{T} (\Gamma_a - \Gamma)
N2=Tg(Γa−Γ)
另一种表示是用位温:
N
2
=
g
θ
∂
θ
∂
z
\colorbox{yellow}{$N^2 = \frac{g}{\theta} \frac{\partial \theta}{\partial z}$}
N2=θg∂z∂θ
当 ( N 2 > 0 N^2 > 0 N2>0 ) 时，大气是稳定的。
- 在热带和温带地区的对流层中，N 的典型值约为( 1 0 − 2 s − 1 10^{-2} s^{-1} 10−2s−1)，这对应于一个大约为 Γ = 6 K / k m \Gamma=6K/km Γ=6K/km的温度梯度。
- N是描述天气系统的一个重要参数，通常也被称为“层结”（stratification），它表示大气层的稳定程度。
1.2 Horizontal structure
在全球大气系统中，温度和位温（potential temperature）随高度和纬度的变化决定了大气的稳定性和风场的分布。本节通过长期气候数据和再分析数据，探讨了大气的 温度、位温结构以及风场的特征 \colorbox{#E0FFFF}{温度、位温结构以及风场的特征} 温度、位温结构以及风场的特征。
1.2.1 大气温度(T)和位温( θ \theta θ)的水平结构
根据全球气候分析数据，我们可以观察到温度和位温 θ \theta θ随纬度和高度的变化特征：
-
温度随高度的变化：温度在对流层中随高度下降，特别是在热带地区的对流层顶附近（约100 hPa 或10 kPa），温度降至最低值，通常低于200K。而在整个大气中，位温 θ \theta θ随高度增加，显示出大气的稳定分层特征。 ∂ T ∂ z < 0 \frac{\partial T}{\partial z} < 0 ∂z∂T<0，位温等值线在赤道附近向下弯曲（说明赤道地区的位温较高），而温度等值线则向上拱起。
-
等位温面（Isentropic surfaces）：在对流层中，等位温面倾斜，从热带向极地逐渐上升。这意味着空气团如果在绝热条件下从低纬度移动到高纬度，它将在等位温面上升高，这解释了空气团的浮力变化和垂直运动特征。
一些重要的点包括：
- 位温面在对流层顶以上非常接近，表明平流层的静力稳定性比对流层大约高两倍。
- 热带对流层顶常接近一个恒定的位温值，通常约为380K。
2. 经向风场与大气环流结构
图2展示了大气的经向风场结构，重点关注的是冬季（DJF季节）的平均风速和位温分布：
-
中纬度的西风急流：北半球冬季在中纬度地区可以观察到显著的西风急流。图中显示出北半球的急流强度较大，对应着对流层顶高度的突然变化，这被称为“副热带急流”（subtropical jet）。
-
300K的等位温面：在热带地区，300K的等位温面接近地表（大约对应27°C的温度），但随着纬度增加，这个等位温面向对流层顶倾斜。这被称为“分隔等位温面”（dividing isentrope），尽管其存在的原因尚不完全明确，但它表示了低纬度和高纬度之间的重要分界。
-
热力平衡与急流结构：在中纬度地区，随着高度的增加，温度梯度和风速都变强。这种变化保持了热力学平衡（thermal wind balance），即随着纬度上温度的降低，水平温度梯度增加，进而导致更强的西风急流。
此外，另一张高分辨率的观测图展示了苏格兰上空的急流结构，其最大风速可达90 m/s（约200 mph），显示了急流结构的剧烈变化以及热力平衡的维持情况。
1.3 Large scale waves in mid-latitudes
1.4 Large scale waves in the tropics
1.5 Cyclones
1.6 Dominance of rotational flow
一个关键点是，这些风的运动主要受旋转流动（rotational flow）控制，且几乎是无辐散的（non-divergent）。
在大气与海洋流体动力学模块中，我们引入了描述流动旋转成分的量——涡度（vorticity）。通常，任何水平流动可以分解为两部分（称为 Helmholtz 分解）：
V ⃗ = V ⃗ non-divergent + V ⃗ irrotational \vec{V} = \vec{V}_{\text{non-divergent}} + \vec{V}_{\text{irrotational}} V=Vnon-divergent+Virrotational
V ⃗ = k ^ × ∇ h ψ + ∇ h χ \vec{V} = \hat{k} \times \nabla_h \psi + \nabla_h \chi V=k^×∇hψ+∇hχ
其中， ψ \psi ψ 称为流函数（streamfunction）， χ \chi χ 为速度势（velocity potential）， k ^ \hat{k} k^是垂直方向的单位向量， ∇ h \nabla_h ∇h 表示水平梯度算符，在笛卡尔坐标系 ( x , y ) (x, y) (x,y)中展开为：
∇ h = i ^ ∂ ∂ x + j ^ ∂ ∂ y \nabla_h = \hat{i} \frac{\partial}{\partial x} + \hat{j} \frac{\partial}{\partial y} ∇h=i^∂x∂+j^∂y∂
旋转流动部分沿着流函数的等值线方向运动，其速度等于流函数梯度的大小。在笛卡尔坐标系下，速度分量可以表示为：
u = − ∂ ψ ∂ y + ∂ χ ∂ x , v = ∂ ψ ∂ x + ∂ χ ∂ y u = -\frac{\partial \psi}{\partial y} + \frac{\partial \chi}{\partial x}, \quad v = \frac{\partial \psi}{\partial x} + \frac{\partial \chi}{\partial y} u=−∂y∂ψ+∂x∂χ,v=∂x∂ψ+∂y∂χ
涡度是流体中某点的局部自旋量，定义为速度的旋度（curl）。其垂直分量可以表示为：
ξ = k ^ ⋅ ( ∇ × u ) = ∂ v ∂ x − ∂ u ∂ y = ∂ 2 ψ ∂ x 2 + ∂ 2 ψ ∂ y 2 = ∇ h 2 ψ \colorbox{yellow}{$\xi = \hat{k} \cdot (\nabla \times \mathbf{u}) = \frac{\partial v}{\partial x} - \frac{\partial u}{\partial y} = \frac{\partial^2 \psi}{\partial x^2} + \frac{\partial^2 \psi}{\partial y^2} = \nabla_h^2 \psi$} ξ=k^⋅(∇×u)=∂x∂v−∂y∂u=∂x2∂2ψ+∂y2∂2ψ=∇h2ψ
可以看出，流函数看起来像是涡度的平滑版本。
流函数和势函数、散度；旋度。。。（待写）
1.7 Planetary rotation and absolute vorticity
在大尺度大气运动中旋转运动占主导的原因是地球自转和大气的强分层。地球围绕南北极之间的轴旋转，自转速率矢量 Ω \Omega Ω 从地心指向北极，其大小等于地球的自转速率。
绝对涡度（absolute vorticity）是行星涡度（planetary vorticity）和相对于地球表面的流动涡度之和：
ζ
=
2
Ω
⃗
⏟
p
l
a
n
e
t
a
r
y
+
∇
⃗
×
u
⃗
⏟
r
e
l
a
t
i
v
e
\zeta=\underbrace{2\vec{\Omega}}_{planetary}+\underbrace{\vec\nabla\times\vec{u}}_{relative}
ζ=planetary
2Ω+relative
∇×u
绝对涡度的垂直分量为：
ζ = 2 Ω sin ϕ + ξ \zeta = 2\Omega \sin\phi + \xi ζ=2Ωsinϕ+ξ
其中， ϕ \phi ϕ 是纬度， f = 2 Ω sin ϕ f = 2\Omega \sin\phi f=2Ωsinϕ是科氏参数（Coriolis parameter）， ξ \xi ξ是相对涡度（relative vorticity）。需要注意的是，根据数学约定，正的涡度与绕着自转轴逆时针旋转（当从上方看时）相关。因此，如果相对涡度的垂直分量在地球表面为正，则从上方看，流体必须绕该点逆时针运动。
如果相对涡度 ξ \xi ξ 与科氏参数 f f f 具有相同符号，则称这种流动为“气旋式”（cyclonic）。这意味着如果气旋逆时针旋转（在北半球），那么相对涡度在北半球为正，而在南半球则为负。
1.8 Geostrophic balance（待完善）
1.8.1 地转风平衡的背景
为了理解为什么在大尺度上旋转流动是主导的，我们需要回到相对于地球旋转表面的水平动量方程(牛顿第二定律)：
D V ⃗ D t + f k ⃗ × V ⃗ = − 1 ρ ∇ h p \colorbox{yellow}{$\frac{D\vec{V}}{Dt} + f\vec{k} \times \vec{V} = -\frac{1}{\rho}\nabla_h p$} DtDV+fk×V=−ρ1∇hp
其中， V \mathbf{V} V 是速度矢量， ρ \rho ρ 是密度， p p p 是气压， f f f 是科氏参数（Coriolis parameter）。
在整个课程中，我们使用局部平面近似（local planar approximation），这样可以使用局部的笛卡尔坐标系 ( x , y , z ) (x, y, z) (x,y,z)，而不需要处理球面几何的复杂性。不过， 这种近似仅在典型水平长度尺度 L 远小于地球半径 a 时有效 。 \colorbox{#E0FFFF}{这种近似仅在典型水平长度尺度 $L$远小于地球半径 $a$时有效}。 这种近似仅在典型水平长度尺度 L远小于地球半径 a时有效。
罗斯贝数（Rossby Number, (Ro)）通过对上述方程的左侧项进行无量纲化而得到：
R o = i n e r t i a l r o t a t i o n = V f 0 L Ro =\frac{inertial}{rotation}= \frac{V}{f_0 L} Ro=rotationinertial=f0LV
其中，(V) 是特征的水平风速。罗斯贝数越小，行星自转对所考虑运动类型的影响越大。具有小罗斯贝数的运动将经历地转调整（geostrophic adjustment）。当某个初始时间点出现压力异常（例如压力梯度指向较低压力的方向）时，流体会沿压力梯度相反的方向加速，但由于科氏力的存在，最终趋向于地转平衡。
1.8.2 地转风的定义
在小罗斯贝数情况下（即 V f 0 L \frac{V}{f_0 L} f0LV很小），局部加速度项相对较小，大气处于地转平衡（geostrophic balance），此时压力梯度力和科氏力相等，保持地转风在当前状态下平衡：
f k ^ × V g ⃗ = − 1 ρ ∇ h p \colorbox{yellow}{$f\hat{k} \times \vec{V_g} = -\frac{1}{\rho}\nabla_h p$} fk^×Vg=−ρ1∇hp
假设压力扰动 (p’) 在水平上变化比科氏参数 (f) 或密度 (\rho) 快得多，我们可以定义地转流函数（geostrophic streamfunction） ψ g \psi_g ψg：
ψ g = p ′ ρ 0 f 0 \colorbox{yellow}{$\psi_g = \frac{p'}{\rho_0 f_0}$} ψg=ρ0f0p′
地转流函数 ψ g \psi_g ψg 的推导(可跳过)
表达地转风的速度分量
首先，将地转平衡方程展开成分量形式。我们用速度分量 V g ⃗ = ( u g , v g ) \vec{V_g} = (u_g, v_g) Vg=(ug,vg) 表示地转风，那么：
f ( − v g u g ) = − 1 ρ ( ∂ p ∂ x ∂ p ∂ y ) f \begin{pmatrix} -v_g \\ u_g \end{pmatrix} = -\frac{1}{\rho} \begin{pmatrix} \frac{\partial p}{\partial x} \\ \frac{\partial p}{\partial y} \end{pmatrix} f(−vgug)=−ρ1(∂x∂p∂y∂p)
展开得到两个分量方程：
f v g = 1 ρ ∂ p ∂ x 和 f u g = − 1 ρ ∂ p ∂ y f v_g = \frac{1}{\rho} \frac{\partial p}{\partial x} \quad \text{和} \quad f u_g = -\frac{1}{\rho} \frac{\partial p}{\partial y} fvg=ρ1∂x∂p和fug=−ρ1∂y∂p3. 定义地转流函数
为了将压力梯度和风速与一个标量函数（即地转流函数 (\psi_g)）关联，我们做如下假设：
u g = − ∂ ψ g ∂ y , v g = ∂ ψ g ∂ x u_g = -\frac{\partial \psi_g}{\partial y}, \quad v_g = \frac{\partial \psi_g}{\partial x} ug=−∂y∂ψg,vg=∂x∂ψg4. 地转流函数的关系式
我们将以上速度分量的表达式代入地转平衡方程：
f ∂ ψ g ∂ x = 1 ρ ∂ p ∂ x 和 − f ∂ ψ g ∂ y = 1 ρ ∂ p ∂ y f \frac{\partial \psi_g}{\partial x} = \frac{1}{\rho} \frac{\partial p}{\partial x} \quad \text{和} \quad -f \frac{\partial \psi_g}{\partial y} = \frac{1}{\rho} \frac{\partial p}{\partial y} f∂x∂ψg=ρ1∂x∂p和−f∂y∂ψg=ρ1∂y∂p
这两个方程可以重写为：
∂ ψ g ∂ x = 1 f ρ ∂ p ∂ x 和 ∂ ψ g ∂ y = − 1 f ρ ∂ p ∂ y \frac{\partial \psi_g}{\partial x} = \frac{1}{f \rho} \frac{\partial p}{\partial x} \quad \text{和} \quad \frac{\partial \psi_g}{\partial y} = -\frac{1}{f \rho} \frac{\partial p}{\partial y} ∂x∂ψg=fρ1∂x∂p和∂y∂ψg=−fρ1∂y∂p5. 将地转流函数与气压扰动关联
为了进一步简化这个关系，我们假设密度 (\rho) 在局部范围内近似为常数，即 (\rho \approx \rho_0)。这样，我们可以定义一个新的量——气压扰动 (p’)，代表相对于一个基准值（例如平均气压）的变化：
p ′ = p − p 0 p' = p - p_0 p′=p−p0
在这种情况下，我们可以定义地转流函数 ψ g \psi_g ψg 为：
ψ g = p ′ ρ 0 f 0 \psi_g = \frac{p'}{\rho_0 f_0} ψg=ρ0f0p′
其中 f 0 f_0 f0 是一个常数（通常是参考纬度的科氏参数）。
这个定义将地转流函数 ψ g \psi_g ψg与气压扰动直接关联起来，从而简化了地转平衡的数学表达，并且使得我们可以利用 ψ g \psi_g ψg来描述风速的分量：
u g = − ∂ ψ g ∂ y , v g = ∂ ψ g ∂ x u_g = -\frac{\partial \psi_g}{\partial y}, \quad v_g = \frac{\partial \psi_g}{\partial x} ug=−∂y∂ψg,vg=∂x∂ψg
1.2.3 常数f地转平衡与涡度
假设 (f_0) 是一个常数，这被称为“常数f地转平衡”（constant-f geostrophy）。需要注意的是，由于科氏参数 f f f 随纬度变化，且加速度项并不为零（而是与罗斯贝数同阶）， ψ g \psi_g ψg 并不等于完整的流函数。
根据这个定义，地转风是无辐散的（non-divergent）：
∇ ⋅ V g ⃗ = ∂ u g ∂ x + ∂ v g ∂ y = − ∂ 2 ψ g ∂ x ∂ y + ∂ 2 ψ g ∂ y ∂ x = 0 \nabla \cdot \vec{V_g} = \frac{\partial u_g}{\partial x} + \frac{\partial v_g}{\partial y} = -\frac{\partial^2 \psi_g}{\partial x \partial y} + \frac{\partial^2 \psi_g}{\partial y \partial x} = 0 ∇⋅Vg=∂x∂ug+∂y∂vg=−∂x∂y∂2ψg+∂y∂x∂2ψg=0
此外，我们可以定义地转相对涡度（geostrophic relative vorticity）向量 ξ g \xi_g ξg：
ξ g = ∇ × V g ⃗ = ( − ∂ v g ∂ z , ∂ u g ∂ z , ∂ v g ∂ x − ∂ u g ∂ y ) \xi_g = \nabla \times \vec{V_g} = \left(-\frac{\partial v_g}{\partial z}, \frac{\partial u_g}{\partial z}, \frac{\partial v_g}{\partial x} - \frac{\partial u_g}{\partial y}\right) ξg=∇×Vg=(−∂z∂vg,∂z∂ug,∂x∂vg−∂y∂ug)
由于地转风没有垂直分量，因此涡度在x和y方向上只有一个分量。其垂直分量可以表示为：
ξ g = ∂ v g ∂ x − ∂ u g ∂ y = ∇ h 2 ψ g \colorbox{yellow}{$\xi_g = \frac{\partial v_g}{\partial x} - \frac{\partial u_g}{\partial y} = \nabla_h^2 \psi_g$} ξg=∂x∂vg−∂y∂ug=∇h2ψg
扫一扫
1937
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
