动力气象-大气运动与尺度分析（2024.9 - 2025.1）

前言

本课程使用教材为《动力气象学》（气象出版社）文章根据课堂笔记和教材内容整理订正编写，文章内容如有错误，请指正。
动力气象与天气学原理不割裂，动力气象在天气学原理的基础上，增加了更多数学物理计算和推导的过程，本章部分内容在《天气学原理》已有涉及，可以与关于天原的文章结合分析。

关于热成风、p坐标系、地转风、地转偏差的内容详见：
天气学原理-大气运动基本特征

---

动力气象学研究问题基本假设

1. 连续介质假设：不考虑大气的分子结构，大气是连续流体
2. 理想气体假设和可压缩性质： 大气是理想气体(没有凝结)，是可压缩的(散度不等于0)

研究对象及其性质

中高纬度多尺度的大气运动，目前动力气象主要理论集中在大尺度运动(天气尺度 Synoptic scale)的大气

大尺度大气的运动特征：准定常、准水平、准静态、准地转、无辐散

1. 考虑地球自转，即地转偏向力，大气准地转
2. 考虑地球重力作用，大气垂直运动速度很小，大气准水平运动，同时大气重力与垂直气压梯度力近似平衡，即准静力平衡
3. 大气是层结流体，大气密度随高度变化，层结分为稳定与不稳定层结
4. 大气中含有水分，对流凝结会释放潜热
5. 大气下边界不均匀，海陆地形分布影响大气，同时边界层大气运动具有湍流性

大气运动基本方程组

定义一些地球基本常数
$$\begin{array}{r}\Omega =7.292\times 10^{-5}{s}^{-1}\\ r=6.371\times 10^{6}m\\ p=1.01\times 10^{5}pa\\ ...\end{array}\text{\hspace{1em}(地转角速度)(地球半径)(标准大气压)}$$

任意矢量对时间的导数

矢量求导，定义 $\vec{A}$ 为任意矢量

$$\frac{d_a\vec{A}}{dt}=\frac{dA}{dt}\vec{i}+\frac{dA}{dt}\vec{j}+\frac{dA}{dt}\vec{k}+A_x\frac{d_a\vec{i}}{dt}+A_y\frac{d_a\vec{j}}{dt}+A_z\frac{d_a\vec{k}}{dt}$$

画图推导，容易得到

$$\frac{d_a\vec{i}}{dt}=\vec{\Omega }\times \vec{i}$$

其中 $\vec{\Omega }$ 为地球旋转角速度，其他方向同理

可得到

$$\frac{d_a\vec{A}}{dt}=\frac{d\vec{A}}{dt}+\vec{\Omega }\times \vec{A}$$

可以视为两个坐标系矢量对时间求导的转换算子

借助此算子可以计算速度和加速度公式

$$\frac{d_a\overrightarrow{V_a}}{dt}=\frac{d\vec{V}}{dt}+2\vec{\Omega }\times \vec{V}+\vec{\Omega }\times (\vec{\Omega }\times \vec{R})$$

绝对加速度 = 相对加速度 + 科氏加速度 + 向心加速度

运动方程

将上述加速度公式代入到牛顿第二定律可得到

$$\frac{d\vec{V}}{dt}=\sum \vec{F}+{\Omega }^2\vec{R}-2\vec{\Omega }\times \vec{V}$$

相对加速度 = 真实力 + 惯性离心力 + 地转偏向力

真实力：质量力(万有引力)与表面力(气压梯度力、摩擦力)

$$\sum \vec{F}=-\frac{GM}{r^2}\frac{\vec{r}}{r}-\frac{1}{\rho }\nabla p+\overrightarrow{F_f}$$

非真实力：惯性离心力、地转偏向力、重力(万有引力与惯性离心力的合力)

最终合并得到基本的矢量表达的运动方程

$$\frac{d\vec{V}}{dt}=-\frac{1}{\rho }\nabla p-2\vec{\Omega }\times \vec{V}+\vec{g}+\overrightarrow{F_f}$$
右式后四项分别为气压梯度力项、地转偏向力项、重力项和摩擦力项。

质量守恒与连续方程

Lagrange观点方程

$$\frac{d\rho }{dt}+\rho \nabla \cdot \vec{V}=0$$

Eular 观点方程

$$\frac{\partial \rho }{\partial t}+\nabla \cdot (\rho \vec{V})=0$$

其中 $\nabla \cdot \vec{V}$ 为速度散度，对应单位体积体胀速度，即某一气块膨胀和收缩的程度；

$\nabla \cdot \rho \vec{V}$ 为质量散度，对应单位体积流体质量通量，即某一空间内流体质量流入流出之间的差异。

热力学能量方程与状态方程

干空气状态方程，动力气象学一般研究干空气，不考虑水汽

$$p=\rho RT⟺\alpha p=RT$$

其中 $\alpha =\frac{1}{\rho }$ 为比容

根据热力学第一定律，得到热力学能量方程

$$c_v\frac{dT}{dt}+p\frac{d\alpha }{dt}=\frac{dQ}{dt}$$

其中 $c_v=718J/(kg\cdot K)$ 为干空气定容比热

气块吸收热量，内能增加并且对外做功

状态方程对时间求导，代入上述热力学能量方程，可得到热力学能量方程的另一种表达方式

$$c_p\frac{dT}{dt}+R\frac{dp}{dt}=\frac{dQ}{dt}$$

其中 $c_p=1005J/(kg\cdot K)$ 为干空气定压比热

位温及其表示的热力学方程

大气干绝热过程与泊松方程
$$\frac{T_0}{T}=(\frac{p_0}{p}{)}^{R/c_p}$$

气块干绝热过程温度递减率
$${\gamma }_d=-\frac{dT}{dz}\approx 9.8K/km$$
环境温度递减率，对流层平均而言温度递减率大约6.5K/km
$$\gamma =-\frac{\partial T}{\partial z}\approx 6.5K/km$$

位温：气压为 $p$，温度为 $T$ 的气块绝热移动到海平面 $p_s=1000hpa$ 时对应的温度
$$\theta =T(\frac{p_s}{p}{)}^{R/c_p}$$

取对数，对时间求导，可以得到热力学能量方程，也可以证明干绝热过程中位温守恒

$$\frac{d\ln \theta }{dt}=\frac{1}{c_{p}T}\frac{dQ}{dt}$$

运动方程、连续方程、状态方程、热力学能量方程联立称为大气运动方程组。

尺度分析

尺度：反映物理量的量值，例如速度 $u$
$$u=U{u}^{\ast }$$

其中 $U$ 为特征值有量纲，$u^{\ast }$ 为无量纲量

动力学尺度

动力学尺度主要可以分为空间尺度、时间尺度和运动尺度

空间尺度分为水平尺度 $L$ 和垂直厚度尺度 $D$

运动尺度分为水平速度尺度 (特征值为 $U$ ) 和垂直速度尺度(特征值为 $W$ )，观测表明，速度场的时空变动尺度可以达到本身的量级

时间尺度：系统演变经历一个阶段所需要的特征时间，对于大尺度系统 $\tau =L/U$，对于中小尺度系统 $\tau >L/U$

大气运动各个尺度特征值如下：

运动类型	$L(m)$	$D(m)$	$U(m/s)$	$\tau (s)$
大尺度	$10^6$	$10^4$	$10$	$10^5$
中尺度	$10^5$	$10^4$	$10$	$10^5$
小尺度	$10^4$	$10^3\sim 10^4$	$10$	$10^4$

热力学尺度

热力学变量：包括气压、密度、位温、温度等，热力学变量分解为表征其基本状态的基本热力学变量(仅与 $z$ 有关)和扰动量

$$p=\overline{p(z)}+p^{\prime }(x,y,z,t)\sim P+\Delta P$$

由于其时空变动值相对于其本身很小，扰动量远小于其本身的量级，即 $\frac{\Delta p}{p}<<1$

基本态的垂直厚度尺度(标高)：$H=R{T}^{\ast }/g\sim 10^{4}m$，一般是指随着高度增加，气压减少为起始高度气压的 $1/e$ 时对应的高度，注意，它隶属于热力学变量的尺度，不同于空间尺度中的垂直厚度尺度 $D$，并且 $D\le H$

热力学基本变量随高度的变化可以达到其本身的量级

$$\frac{d\bar{p}}{dz}=-\bar{\rho }g\sim \pi g\sim \frac{P}{H}$$

其中 $\pi$ 为密度 $\rho$ 的特征值

扰动热力学变量的时空变动值可以达到扰动量本身的量级，例如对于气压

$$\frac{\partial p}{\partial x}=\frac{\partial p^{\prime }}{\partial x}\sim \frac{\Delta P}{L}$$

水平运动方程尺度分析

将水平运动方程展开到两个方向，并忽略摩擦力，气压梯度力中的密度项可以如下化简(Taylor展开法)

$$-\frac{1}{\rho }=-\frac{1}{\bar{\rho }+{\rho }^{\prime }}\approx -\frac{1}{\bar{\rho }}(1-\frac{{\rho }^{\prime }}{\bar{\rho }})⟹-\frac{1}{\rho }\sim \frac{1}{\bar{\rho }}$$

代入水平运动方程得到三项

$$\begin{array}{rl}\frac{\Delta p}{\pi L}& \\ f_{0}U& \\ U^2/L& \end{array}\text{\hspace{1em}(特征气压梯度力)(特征地转偏向力)(特征速度平流)}$$
可以得到
$$\frac{\Delta p}{\pi L}\sim f_{0}U+\frac{U^2}{L}$$

定义 Rossby数：水平惯性力与地转偏向力的比值

$$Ro=\frac{U^2/L}{f_{0}U}=\frac{U}{f_{0}L}$$

• 当 $Ro\ge 1$ 则 $\frac{\Delta p}{f_{0}L}\sim \frac{U^2}{L}$
• 当 $Ro<<1$ 则 $\frac{\Delta p}{f_{0}L}\sim f_{0}U$ ，大气运动加速度很小，满足准地转，为大尺度大气运动的特征

定义Kibel数：惯性运动时间尺度和运动时间尺度的比值
$$\epsilon =\frac{U/\tau }{fU}=1/f\tau$$
其中 $1/f$ 可以看作惯性运动的特征时间尺度

当 $\epsilon <<1$ 时，表示慢过程，准定常，反之为快过程非定常

连续方程尺度分析

对连续方程做一定的转换

$$\frac{d\ln \rho }{dt}+\frac{\partial u}{\partial x}+\frac{\partial v}{\partial y}+\frac{\partial w}{\partial z}=0$$

将密度项展开为基本态和扰动项

$$\frac{d\ln \rho }{dt}=\frac{d\ln \bar{\rho }}{dt}+\frac{d}{dt}(\frac{{\rho }^{\prime }}{\bar{\rho }})$$

并再局地展开和忽略小项，原方程可以得到

$$\frac{\partial u}{\partial x}+\frac{\partial v}{\partial y}\sim \frac{\partial w}{\partial z}$$

因为根据观测已知

$$\frac{\partial u}{\partial x}+\frac{\partial v}{\partial y}\le \frac{U}{L}$$

得到
$$\frac{W}{D}\le \frac{U}{L}⟹W\le \frac{UD}{L}$$

可以说明大尺度大气运动准水平的原因，同时也能给出 $W$ 的上限

零级近似和一级近似

将方程每一项进行尺度分析，之后代入大尺度运动的特征值，进行计算

运动方程简化

零级近似：保留特征值最大项

$$\{\begin{array}{ll}-\frac{1}{\rho }\frac{\partial p}{\partial x}+fv=0& \\ -\frac{1}{\rho }\frac{\partial p}{\partial y}-fu=0& \\ -\frac{1}{\rho }\frac{\partial p}{\partial z}-g=0& \text{(静力平衡)}\end{array}$$

零级近似是准定常的，不能作为预报方程，可以作为诊断方程，揭示了中纬度大气的基本特征

一级近似：保留特征值最大项和相比于最大项小一个量级的项
$$\{\begin{array}{l}\frac{\partial u}{\partial t}+u\frac{\partial u}{\partial x}+v\frac{\partial u}{\partial y}=-\frac{1}{\rho }\frac{\partial p}{\partial x}+fv\\ \frac{\partial v}{\partial t}+u\frac{\partial v}{\partial x}+v\frac{\partial v}{\partial y}=-\frac{1}{\rho }\frac{\partial p}{\partial y}-fu\\ -\frac{1}{\rho }\frac{\partial p}{\partial z}-g=0\end{array}$$

连续方程简化

零级近似
$$\frac{\partial u}{\partial x}+\frac{\partial v}{\partial y}=0$$

反映中纬度大气具有准水平无辐散的特点

一级近似：保留原来的形式

热力学方程简化

零级近似：

$$\frac{\partial T}{\partial t}+u\frac{\partial T}{\partial x}+v\frac{\partial T}{\partial y}+w({\gamma }_d-\gamma )=0$$
此简化方程表明，在天气尺度中，温度局地变化主要由温度平流和温度垂直绝热变化引起

一级近似即为方程本身，没有变化

地转参数简化

地转参数随着纬度非线性变化，对地转参数在y方向上进行Taylor展开
$$f=f_0+\beta y+o(y)$$
其中 $\beta =\frac{2{\Omega }_0\cos {\varphi }_0}{a}=\frac{df}{dy}{|}_{{\varphi }_0}$

$f$ 平面近似：纬度变化不大的情况下，近似认为地转参数 $f$ 为常数

$\beta$ 平面近似：纬度变化不大的情况下，近似认为地转参数 $f$ 的变化为线性，即 $\beta =\frac{\partial f}{\partial y}$ 为常数。