葛朝霞等,气象学与气候学教程,中国水利水电出版社(第2版)
- 4.7 水平运动与垂直运动的关系
- 5.4 水分方程和可降水量
-
大气河流3D model: https://www.bilibili.com/video/BV12Y411b7XY/, https://svs.gsfc.nasa.gov/4960
-
王光谦院士答天河工程:https://www.qhu.edu.cn/jxky/51516.htm
君不见黄河之水天上来,奔流到海不复回。-- 李白
天上的水,从哪里来?
https://www.ndrc.gov.cn/fggz/fzzlgh/gjjzxgh/201604/P020191104623960850846.doc
- 欧拉:着眼于空间。温度$T$随时间$t$和空间位置(
$x, y, z$ )变化,$T = f(t, x, y, z)$,变化速率:
- 拉格朗日:着眼于质元。
$T$ 只随时间$t$变化,T = g(t),变量速率:$\frac{dT}{dt}$
1.$\nabla = \frac{\partial f}{\partial x} \boldsymbol{i} + \frac{\partial f}{\partial y} \boldsymbol{j} + \frac{\partial f}{\partial y} \boldsymbol{k}$,
$\nabla_h = \frac{\partial f}{\partial x} \boldsymbol{i} + \frac{\partial f}{\partial y} \boldsymbol{j}$ 2.拉格朗日描述与欧拉描述,https://www.zhihu.com/question/26129680
如果是密度$\rho$,公式也是同样的形式:
密度随时间的变化:可联想污染物的扩散,污染物密度的变化
$x$ 方向净流入量:
$$ \rho u \delta y \delta z \delta t - (\rho u + \frac{\partial \rho u}{\partial x} )\delta y \delta z \delta t =
- \frac{\partial \rho u}{\partial x} \delta y \delta z \delta t $$
同理可以得出$y$和$z$方向的流入量, 分别为:$- \frac{\partial \rho v}{\partial y} \delta x \delta z \delta t$、$- \frac{\partial \rho w}{\partial z} \delta x \delta y \delta t$,因此,
$$ \frac{\partial \rho}{\partial t} \delta x \delta y \delta z \delta t =
- ({\partial \rho u \over \partial x} + {\partial \rho v \over \partial y} + {\partial \rho w \over \partial z}) \delta x \delta y \delta z \delta t $$
$$ \begin{align*} \ \frac{\partial \rho}{\partial t} + {\partial \rho u \over \partial x} + {\partial \rho v \over \partial y} + {\partial \rho w \over \partial z} = 0 \
\frac{\partial \rho}{\partial t} + \nabla · (\rho V) = 0 \ \frac{\partial \rho}{\partial t} + V \nabla · \rho + \rho \nabla · V= 0 \
\end{align*} $$
:::info:Recall $$ \frac{d \rho}{dt} = {\partial \rho \over \partial t} + V \nabla · \rho $$ :::
1.连续方程的数学推导, https://www.bilibili.com/video/BV1d64y1Q7iw
拉格朗日:对于质元(研究对象恒定),只要体积不变(不可压缩),则$\rho$不变。
流体不可压缩,则其流动过程中密度不变(
地表和对流层顶,$w = 0$
-
$\mathbf{Div} > 0$ ,物质丢失,辐散;$\mathbf{a_w}$速度不断减小(如$-2m/s$到$-3m/s$)- 情况发生在地表(
$w = 0$ ),向上速度减小,因此向上走$w$为负,下沉运动; - 中间层$w$不确定,要看辐散的强度
- 情况发生在地表(
-
$\mathbf{Div} < 0$ ,物质积累,辐合;$\mathbf{a_w}$速度不断增加(如$-2m/s$到$1m/s$)- 情况发生在地表(
$w = 0$ ),向上速度变大,因此向上走$w$为正,上升运动; - 中间层$w$不确定,要看辐合的强度
- 情况发生在地表(
不可压缩气体:
-
z坐标中($t, x, y, z$),$w = \frac{dz}{dt}$
-
p坐标中($t, x, y, p$)
$$ {\partial u \over \partial x} + {\partial v \over \partial y} + {\partial \omega \over \partial p} = 0 $$
$\omega = \frac{dp}{dt}$ :p坐标中,垂直方向上的速度-
$\omega < 0$ ,向上运动; -
$\omega > 0$ ,向下运动。
-
- 向上运动:$w > 0$,$\omega < 0$,辐合,降水
- 向下运动:$w < 0$,$\omega > 0$,辐散,高温干旱
定义:单位面积上,整个大气柱中的水汽全部凝结,降落至地面的水量(单位$mm$)。
英文:Precipitable Water Vapor
, Vertical Integrated Water Vapor
- 单位面积(
$A=1m^2$ ),厚度为dz的湿空气块,其所具有的水量:
$$ \begin{align*} dW &= \rho_v A dz, (dp = - \rho g dz)\ &= -\frac{1}{g}\frac{\rho_v}{\rho} dp = -\frac{1}{g} q dp
\end{align*} $$
- 从地面($z=0, p=P_s$)到高空($z=z, p=P_t$),对其水汽含量进行积分:
上述公式没有解析解,只能以离散的形式进行求解:
把大气划分为$n$个区间($1$到$n+1$层)。
Recall:
$S = \int{f(x) dx}$ 的求解方法
- 从第$i$层到第$i+1$层的水汽含量:
- 因此整层的大气可将水量:
大气可降水量 -> 潜在的大气河流
Gimeno, et al., 2012, Reviews of Geophysics
定义:单位时间、单位面积,通过的水汽质量。
单位:$kg~m^{-2}s^{-1}$
- 通过的气体体积(V):(
$\Delta P × \Delta L) × v \Delta t$ - 通过的水汽气体质量(
$m_v$ ):$m_v = m q = q \rho (\Delta P × \Delta L × v \Delta t)$
由于是单位面积、单位时间,因此$m_v$需要除以$A$和$\Delta t$,因此水汽通量:
可以根据$q v \rho$进一步验证$f$的单位:$kg~m^{-2}s^{-1}$,:grinning:。
定义:单位时间、单位宽度,通过整层气柱的水汽质量,单位$kg~m^{-1}s^{-1}$。
单位:kg/(m·s)
与大气可降水量相比,整层水汽通量多乘了速度$v$。
- 水汽通量
- 水汽通量散度
$$ \begin{align*} d &= \nabla_h (q \rho v) \ \mathbf{Div} &= -\int_{p_s}^{p_t} \nabla_h(q \rho v) \frac{d p}{\rho g} \ &= -\frac{1}{g} \int_{p_s}^{p_t} \nabla_h(q v) d p
\end{align*} $$
- 水汽通量散度
-
$\mathbf{Div} < 0$ ,进多出少(降水),辐合,气柱的水汽增加+ : -
$\mathbf{Div} > 0$ ,进少出多(蒸发),辐散,气柱的水汽流失-
$\mathbf{Div}$ 可以类比径流。
Ma, 2022, GPC.
-
$100 kg ~m^{-2} s^{-1}$ 转化为$mm/d$
-
$1×10^{-6} g·cm^{-2}·s^{-1}·hPa^{-1}$ 转化为$mm/d$假设$850hPa$,可近似认为代表900hPa~800hPa的平均状态,即气压层厚度可以认为是100hPa。若该层的水汽通量散度为
$1×10^{-6} g·cm^{-2}·s^{-1}·hPa^{-1}$ ,转化为$mm/d$大概是多少?
然后乘气压层的厚度(如果认为dp = 100hPa),转化为mm约为$86.4mm/s$;
P = F/A, F = mg (力的单位N = kg m s-2,气压的单位Pa = N /m2 = kg m-1 s-2)
$dp/g$ 单位:$kg$,因此W的单位也是$kg/m^2$.
F的单位: $kg
m/s / m^2$ = $kgm^{-1} s^{-1}$
$\mathbf{Div}$ 的单位: $kgs^{-1} / m^2$ = $kgm^{-2} s^{-1}$
-
大气河流(丁一汇)
-
2011-08-01河南暴雨
-
2021-07-21河南暴雨
较难观测变量:$E$、$\Delta S$、$\Delta W$
- 地表水量平衡:
- 大气水汽平衡,降水水分丢失,蒸发水分补充
其中$RES_s$为storage的残差项,如人类用水活动(生活、工业、农业、水利工程),以及模型的模拟误差;$RES_w$为vapor中的残差项,
两式相加可以得到:
$$
- R - \mathbf{Div} = \Delta W + \Delta S $$
上式提供了验证$E$的另一个思路。How?
-
多年尺度,可以认为$\Delta S = 0$,$\Delta W = 0$,因此$\overline R = - \overline{\mathbf{Div}}$
可以根据径流资料,判断哪个再分析资料的$\mathbf{Div}$质量比较好;
-
根据水汽平衡,得出$E_{atm}$
-
$E_{atm}$ 作为观测值,评估$E$模拟好坏
-
Ma, F., Yuan, X., & Li, H. (2022). Characteristics and circulation patterns for wet and dry compound day-night heat waves in mid-eastern China. Global and Planetary Change, 213, 103839.
-
Roads, J., Kanamitsu, M., & Stewart, R. (2002). CSE Water and Energy Budgets in the NCEP–DOE Reanalysis II. Journal of Hydrometeorology, 3(3), 227–248.
-
Lorenz, C., & Kunstmann, H. (2012). The Hydrological Cycle in Three State-of-the-Art Reanalyses: Intercomparison and Performance Analysis. Journal of Hydrometeorology, 13(5), 1397–1420.
-
丁一汇, 柳艳菊, & 宋亚芳. (2020). 东亚夏季风水汽输送带及其对中国大暴雨与洪涝灾害的影响. 水科学进展, 31(5), 629–643.