气候系统的物理学导论 / Introduction to Climate System Physics

第 1 章：气候系统的物理学导论 / Introduction to Climate System Physics

1.1 引言：从能量收支出发理解气候系统

气候系统不是单一对象，而是一个由多个物理介质、热力学储层与生物地球化学循环耦合而成的层级结构。要理解为什么地球平均表面温度维持在 288 K 而不是按 Stefan-Boltzmann 定律预测的 255 K，要理解为什么二氧化碳浓度的微小变化会触发全球尺度的响应，必须从能量收支方程这一最底层的物理约束出发，并在此基础上建立"强迫—反馈"的形式化语言。

本章的核心任务，是为后续十四章铺设共同的物理基础。我们将依次回答四个问题：

1. 气候系统由哪些子系统构成？它们之间通过哪些通量耦合？ 2. 在没有扰动的理想状态下，全球能量平衡如何决定平衡温度？ 3. 当外部扰动或内部扰动发生时，温度响应如何用线性化的反馈参数 $\lambda$ 与 $\lambda_F$ 描述？ 4. 不同物理过程的时间尺度差异如何决定系统的"记忆"长度？

完成本章后，读者将掌握 IPCC AR6 物理气候系统评估所依赖的最核心概念框架，并具备自行计算全球能量收支、用 Python 复现 0-D 能量平衡模型的能力。

---

1.2 气候系统的五分量结构与耦合机制

1.2.1 子系统的物理定义

气候系统在 WMO 与 IPCC 的标准定义中由五个相互作用的分量组成：大气（Atmosphere）、海洋（Ocean）、冰冻圈（Cryosphere）、陆面（Land Surface）、生物圈（Biosphere）。每个分量都有独特的热容、响应时间与状态变量集合，不能被简化为单一介质。

• **大气**：氮、氧、氩、二氧化碳、水汽及痕量气体与气溶胶的混合气体层。垂直方向按温度廓线分为对流层、平流层、中间层、热层。对流层顶高度随纬度变化（赤道 ~17 km，极地 ~8 km），决定了水汽反馈与对流活动的有效层结。
• **海洋**：覆盖地球表面约 71% 的液态水储层，平均深度 ~3682 m，总热容约为大气的 1000 倍。其上层混合层（~50–200 m）响应季节—年际信号，深层海洋（>1000 m）则储存百年到千年尺度的热量与碳。
• **冰冻圈**：包括海冰、冰川、冰盖、积雪与冻土。其反射率（海冰约 0.5–0.7，新雪可达 0.9）和相变潜热使其成为最强的正反馈来源之一（冰—反照率反馈）。
• **陆面**：岩石圈上层、土壤、植被冠层与地下水。其反照率（沙漠 0.3–0.4，森林 0.1–0.2）、粗糙度与蒸散发能力显著影响地表能量分配。
• **生物圈**：陆地生态系统与海洋生态系统的总称，承载碳、氮、磷循环的生物驱动部分。陆地植被通过光合作用每年吸收约 120 GtC，海洋生物泵约 10 GtC。

1.2.2 耦合通量

五分量之间通过四类基本通量交换：能量（辐射、感热、潜热）、动量（风应力）、水分（蒸发—降水）、碳与其他生物地球化学要素。这些通量的大小、方向与时空分布决定了系统整体的耦合强度。

以能量通量为例，地表净辐射 $R_n$ 可分解为：

$$R_n = (1-\alpha_s) S_d + L^{\downarrow} - L^{\uparrow} - H - \lambda E - G$$

其中 $\alpha_s$ 为地表反照率，$S_d$ 为入射短波辐射，$L^{\downarrow}$ 与 $L^{\uparrow}$ 分别为长波下行与上行辐射通量，$H$ 为感热通量，$\lambda E$ 为潜热通量（$\lambda$ 为汽化潜热），$G$ 为地中热通量。该方程把大气顶（TOA）的辐射收支与近地层的湍流通量显式连接起来，是后续讨论陆面过程与边界层耦合的基础。

碳通量的耦合同样关键。大气 $\mathrm{CO}_2$ 浓度由陆地碳汇（约 $-2.5\,\mathrm{GtC/yr}$）、海洋碳汇（约 $-2.3\,\mathrm{GtC/yr}$）与人为排放（约 $+10\,\mathrm{GtC/yr}$）的代数和决定。当陆地或海洋碳汇响应气候变化而减弱时，大气 $\mathrm{CO}_2$ 的留存比例（airborne fraction, AF）上升，形成"碳—气候"耦合反馈。这是第 8 章的核心内容。

1.2.3 气候系统耦合结构图

下图展示了五分量之间主要的物质能量交换路径。节点采用 ASCII 标识以便在 Mermaid 中稳定渲染。

```mermaid flowchart TD A[Atmosphere] -->|LW radiation| A A -->|Precip, wind| O[Ocean] A -->|Precip, radiation| C[Cryosphere] A -->|Precip, radiation| L[Land Surface] A -->|CO2 exchange| B[Biosphere] O -->|SST, momentum| A O -->|Heat, carbon| C O -->|Carbon uptake| B C -->|Albedo, freshwater| A C -->|Sea level| L L -->|Albedo, ET| A L -->|Carbon, water| B B -->|CO2, latent heat| A B -->|Carbon| O B -->|Vegetation| L ```

这张图并非简单的"清单"，而是揭示了三个对后续章节至关重要的结构特征：第一，反照率与水汽分布主要由大气—冰冻圈—陆面三元耦合决定，因此辐射过程必然涉及多个子系统；第二，海洋是热量与碳的"主储层"，其缓慢的环流调节着整个系统的时间常数；第三，生物圈既是碳汇也是反馈源，其状态会因温度、降水与 $\mathrm{CO}_2$ 的变化而调整。

---

1.3 能量平衡方程与辐射平衡温度

1.3.1 顶层辐射平衡

定义太阳常数为 $S_0 \approx 1361\,\mathrm{W/m^2}$（TSI，在 1 AU 处垂直于太阳光束的辐照度）。地球是一个半径为 $R$ 的球体，截面积为 $\pi R^2$，而接收的辐射需要分布在 $4\pi R^2$ 的球面上。因此球面平均入射短波辐射（SW）为 $S_0/4 \approx 340\,\mathrm{W/m^2}$。

设行星反照率为 $\alpha_p \approx 0.30$，则吸收的短波辐射（ASR）为：

$$\mathrm{ASR} = \frac{S_0}{4}(1-\alpha_p) \approx 240\,\mathrm{W/m^2}$$

在没有大气的情况下，地球作为黑体向外辐射长波，其出射长波辐射（OLR）由 Stefan-Boltzmann 定律给出：

$$\mathrm{OLR} = \sigma T_e^4$$

其中 $\sigma = 5.670 \times 10^{-8}\,\mathrm{W\,m^{-2}\,K^{-4}}$。令 ASR = OLR，得到有效辐射温度（effective radiative temperature）：

$$T_e = \left(\frac{S_0(1-\alpha_p)}{4\sigma}\right)^{1/4} \approx 255\,\mathrm{K} = -18^{\circ}\mathrm{C}$$

这正是黑体行星在没有任何温室效应时的平衡温度。然而实际地表平均温度约为 $T_s \approx 288\,\mathrm{K}$，高出 33 K。这 33 K 的"温室温升"是大气中水汽、二氧化碳、甲烷、臭氧、氧化亚氮等痕量气体对长波辐射吸收与再发射的结果。它不是某种次要修正，而是气候系统存在的核心物理机制。

1.3.2 表面能量平衡与温室效应分解

更精确地，地表能量平衡方程可写为：

$$(1-\alpha_s)S^{\downarrow} + \mathrm{LW}^{\downarrow} = \mathrm{LW}^{\uparrow} + H + \lambda E + G$$

其中 $S^{\downarrow}$ 是到达地面的太阳辐射，$\mathrm{LW}^{\downarrow}$ 是大气向下的红外辐射，$\mathrm{LW}^{\uparrow} = \epsilon_s \sigma T_s^4$ 是地表上行长波辐射，$H$ 是感热，$\lambda E$ 是潜热，$G$ 是地中热通量。在多年平均稳态下，$G \approx 0$（地表储热长期为零），陆面 $H + \lambda E$ 与海面潜热主导再分配。

整个大气柱的净能量收支可以由顶层平衡条件推导：

$$\mathrm{ASR}_{TOA} = \mathrm{OLR}_{TOA}$$

而 OLR 主要由对流层上层（接近 200 hPa）的温度与水汽含量决定。这就是为什么"全球平均地表温度"实际上是辐射平衡温度与对流层温度递减率共同决定的产物：地表可以比辐射平衡温度更暖，因为它不是直接向太空辐射的黑体。

1.3.3 一个简单的 0-D 能量平衡模型

构建 0-D 能量平衡模型的关键假设是：行星反照率 $\alpha_p$ 与 OLR 可以相对表面温度 $T_s$ 线性化展开：

$$\alpha_p(T_s) \approx \alpha_0 + \alpha_T (T_s - T_0)$$

$$\mathrm{OLR}(T_s) \approx A + B T_s$$

参数 $B$ 是 Planck 反馈参数的负值，约为 $B \approx 3.3\,\mathrm{W\,m^{-2}\,K^{-1}}$（即 $T_s$ 上升 1 K，OLR 增加 3.3 W/m²）。参数 $A$ 隐含了温室气体的平均吸收强度。

系统对外部辐射强迫 $F$（单位 W/m²）响应的稳态温度变化为：

$$\Delta T_s = -\frac{F}{B} = \lambda F$$

其中 $\lambda = -1/B \approx 0.30\,\mathrm{K/(W/m^2)}$ 是气候敏感参数（climate sensitivity parameter）的简单 0-D 形式。该数值是 Planck 反馈单独作用下的上限；真实系统的 $\lambda$ 因水汽、冰雪、云等反馈机制叠加而变大（更强的响应）。

1.3.4 数值估算：吸收短波与射出长波

下面用具体数值说明 $340\,\mathrm{W/m^2}$ 的球面平均入射如何被分配：

• 大气顶入射短波（球面平均）：340 W/m²
• 反射回太空（云、地面、气溶胶合计）：≈ 100 W/m²（反射率 0.30）
• 大气顶吸收短波 ASR：≈ 240 W/m²
• 大气顶射出长波 OLR：≈ 240 W/m²（平衡态）
• 地表吸收短波（包含大气透射）：≈ 168 W/m²
• 地表吸收长波（大气逆辐射）：≈ 333 W/m²
• 地表发射长波：≈ 396 W/m²
• 感热 $H$（地表→大气）：≈ 24 W/m²
• 潜热 $\lambda E$（地表→大气）：≈ 81 W/m²

各分量加和后地表净能量为零，整个系统达到热力学稳态。注意 396 W/m² 的地表发射远大于 240 W/m² 的 OLR——多余的辐射被大气吸收并以逆辐射形式部分返回地表，这就是温室效应的本质。

---

1.4 强迫与反馈的形式化理论

1.4.1 辐射强迫的定义

IPCC AR6 采用了有效辐射强迫（Effective Radiative Forcing, ERF）的标准定义：它是在准平衡条件下，地气系统对施加的扰动所产生的 TOA 净辐射通量变化，并允许平流层温度调整但不允许地表温度与对流层调整。其形式化表达为：

$$F = \Delta N_{TOA}\big|_{T_s = T_0, q = q_0, T_{strat} = T_{strat}^{adj}}$$

其中 $N_{TOA}$ 是 TOA 净辐射通量（入减出），$q$ 是水汽分布，$T_{strat}^{adj}$ 是经过快速调整的平流层温度。

与传统的瞬时辐射强迫（Instantaneous Radiative Forcing, IRF）相比，ERF 更接近真实的能量扰动，因为它扣除了快速调整（如平流层冷却对 OLR 的反馈）。例如，$\mathrm{CO}_2$ 倍增的 ERF 约为 $3.93 \pm 0.47\,\mathrm{W/m^2}$（AR6 最佳估计），而其瞬时强迫约为 $3.7\,\mathrm{W/m^2}$——两者的差异来自平流层快速冷却。

1.4.2 主要强迫因子

| 强迫因子 | 1750–2019 累积 ERF (W/m²) | 不确定性范围 (90%) | |---|---|---| | $\mathrm{CO_2}$ | +2.16 | +1.90 ~ +2.41 | | $\mathrm{CH_4}$ | +0.54 | +0.43 ~ +0.65 | | $\mathrm{N_2O}$ | +0.21 | +0.18 ~ +0.24 | | 卤代烃 | +0.41 | +0.33 ~ +0.49 | | 臭氧（平流层 + 对流层） | +0.47 | +0.24 ~ +0.71 | | 气溶胶（总） | −1.1 | −1.7 ~ −0.4 | | 土地利用变化 | −0.20 | −0.30 ~ −0.10 | | 太阳活动 | +0.01 | −0.06 ~ +0.08 | | 人为净强迫 | +2.72 | +2.0 ~ +3.4 |

该表的关键信息是：人为温室气体总强迫约 $+3.3\,\mathrm{W/m^2}$ 被气溶胶冷却效应部分抵消，但净强迫仍显著为正，并且与观测到的全球能量不平衡（约 $+0.6\,\mathrm{W/m^2}$）自洽——多余的能量进入海洋，导致 0–2000 m 深度范围的海洋热含量上升约 $0.015\,\mathrm{K/decade}$。

1.4.3 反馈的形式化

设系统处于 $T_s^0$ 的平衡态，外部施加强迫 $F$ 后温度变为 $T_s^0 + \Delta T_s$。在扰动很小时，辐射响应可线性化为：

$$\Delta N_{TOA} = F + \sum_i \lambda_i \Delta T_s$$

其中 $\lambda_i$ 是第 $i$ 个反馈参数（feedback parameter），单位 $\mathrm{W\,m^{-2}\,K^{-1}}$。平衡时 $\Delta N_{TOA} = 0$，故：

$$\Delta T_s = -\frac{F}{\sum_i \lambda_i} \equiv \lambda F$$

这里 $\lambda = -1/\sum_i \lambda_i$ 是总气候敏感参数。各反馈参数之和为 $\sum \lambda_i$，其符号决定系统稳定性。

主要反馈参数（IPCC AR6 评估，单位 W/m²/K，符号约定：负值=稳定，正值=放大）：

• **Planck（普朗克）反馈**：$-3.3$。任何温度上升都会增加 OLR，是最基本的稳定机制。
• **水汽反馈**：$+1.8$。Clausius-Clapeyron 关系使饱和比湿随温度呈指数增长，增加温室效应。
• **温度递减率反馈**：$-0.6$ 至 $-1.2$。对流层上层增温幅度大于地表，部分抵消水汽效应。
• **地表反照率反馈**：$+0.3$。高纬海冰与积雪退缩降低行星反照率。
• **云反馈**：$+0.42$（AR6 最佳估计）。低云减少与高云增多的净效应，但仍是最大不确定源。

将所有反馈求和后，总反馈参数约为 $\sum \lambda_i \approx -1.3\,\mathrm{W\,m^{-2}\,K^{-1}}$，对应的平衡气候敏感度（ECS）为：

$$\mathrm{ECS} = -\frac{F_{2\times CO_2}}{\sum_i \lambda_i} \approx \frac{3.93}{1.3} \approx 3.0\,\mathrm{K}$$

AR6 评估的"最可能" ECS 范围为 2.5–4.0 K，最佳估计 3.0 K，概率分布 5–95% 为 2.0–5.0 K。

1.4.4 增益因子与放大倍数

一个有用的诊断量是增益因子（gain, $g$），定义为：

$$g = \frac{1}{1 - \sum_i \lambda_i / \lambda_{\text{Planck}}}$$

或者等价地：

$$g = \frac{\Delta T_s}{\Delta T_s^{\text{Planck only}}} = \frac{|\lambda_{\text{Planck}}|}{|\sum_i \lambda_i|}$$

若 $g > 1$，反馈总体上放大响应；若 $g < 1$，则抑制响应。AR6 最佳估计 $g \approx 2.5$，意味着如果没有 Planck 以外的其他反馈，温升仅为当前的 1/2.5。这定量展示了水汽、冰雪、云等反馈过程在气候敏感性中的放大作用。

1.4.5 一个工作示例：$\mathrm{CO_2}$ 加倍的瞬时响应

假设 $\mathrm{CO_2}$ 浓度从 280 ppm 增加到 560 ppm。其瞬时强迫由简化辐射传递关系估算：

$$F_{2\times} \approx 5.35 \ln(560/280) \approx 3.71\,\mathrm{W/m^2}$$

考虑平流层调整后，ERF 约为 $3.93\,\mathrm{W/m^2}$。若仅 Planck 反馈起作用，平衡温度变化为 $3.93 / 3.3 \approx 1.19\,\mathrm{K}$。但加入水汽、云、地表反照率等反馈后，最终温升放大为 3.0 K 左右，放大倍数 $g \approx 2.5$。这一放大幅度决定了为什么几十 ppm 的浓度变化能引起全球尺度的环境响应——并不是辐射直接作用的结果，而是反馈链式放大后的净效应。

---

1.5 多时间尺度耦合与 IPCC 参考框架

1.5.1 时间尺度谱

气候系统的每一类过程都有其特征响应时间，由热容与交换系数的比值决定。下表列出主要过程的