气候系统能量平衡与温室效应物理基础

气候系统的所有动态变化，本质上都可以追溯到能量的输入、储存与再分配过程。作为全球气候系统的核心驱动力，太阳辐射以电磁波形式抵达地球大气层顶端，其光谱分布、能量通量以及与大气的相互作用，构成了整个气候系统能量平衡的物理起点。理解这一过程的每一个细节，是把握气候变化物理机制、评估人为干预效应的基础。

太阳辐射的电磁波谱特性与大气衰减过程

太阳作为距离地球约1.5亿公里的G型主序星，其表面温度约为5778K，可近似为理想黑体，其辐射能量分布遵循普朗克黑体辐射定律。根据维恩位移定律，黑体辐射的峰值波长与温度成反比，即$\lambda_{max} = b/T$，其中维恩常数$b=2898\ \mu m\cdot K$，代入太阳表面温度可得峰值波长约为0.5μm，处于可见光波段（0.4-0.7μm），因此太阳辐射属于短波辐射范畴。在太阳光谱中，约43%的能量集中在可见光波段，49%在近红外波段（0.7-2.5μm），7%在紫外波段（0.2-0.4μm），剩余1%为远红外和射电波段。

在大气层顶端，垂直于太阳光线方向的单位面积接收到的太阳总辐射通量被称为太阳常数，当前卫星观测值为$S_0=1361\pm0.5\ W/m^2$，其长期变化幅度小于0.1%，是气候系统能量输入的核心参数。由于地球为球体，同一时间内只有半个球面接收太阳辐射，且太阳光线与地表的夹角随纬度和时间变化，因此全球平均的入射短波辐射需将太阳常数除以4，即$S_0/4\approx340\ W/m^2$，这一数值是气候系统能量收支计算的基准。

太阳短波辐射穿过大气层时，会与大气中的分子、气溶胶、云滴等发生相互作用，导致能量衰减，主要过程包括三类：第一是吸收，大气中的臭氧（O₃）主要吸收波长小于0.3μm的紫外辐射，水汽（H₂O）、二氧化碳（CO₂）等分子主要吸收近红外波段的辐射，平流层气溶胶和云滴也会吸收部分短波辐射；第二是散射，大气分子的瑞利散射对波长小于1μm的短波辐射作用显著，是天空呈现蓝色的原因，云滴和较大气溶胶的米氏散射则对所有波长的短波辐射都有较强散射作用；第三是反射，云层和对流层气溶胶的反射作用是最主要的短波辐射衰减机制，全球平均而言，约23%的入射太阳辐射被反射回太空，其中云层贡献了约17%，气溶胶贡献了约6%。

不同天气条件下大气的短波辐射透过率差异显著：晴空无云条件下，大气对短波辐射的透过率可达70%-80%，仅少量辐射被吸收和散射；厚云天气下，短波辐射透过率可降至30%以下，几乎一半以上的入射辐射被云层反射。这一特性使得云辐射反馈成为气候系统中最不确定的反馈过程之一，也是当前气候模式模拟差异的主要来源。

地表反照率的调控机制与反馈效应

地表反照率（Albedo）是描述地表反射短波辐射能力的核心参数，定义为地表反射的短波辐射通量与入射短波辐射通量的比值，即$\alpha = F_{reflected}/F_{incident}$，其取值范围为0（完全吸收）到1（完全反射）。不同下垫面的反照率差异极大：新积雪的反照率可达0.8-0.95，沙漠砂石的反照率为0.3-0.4，森林冠层的反照率为0.1-0.15，开阔海洋的反照率仅为0.06-0.1，城市建成区由于建筑、路面等高反射材料的使用，反照率约为0.15-0.25。全球平均地表反照率约为0.3，其中冰雪覆盖区域虽然面积占比不足10%，但对全球反照率的贡献超过30%。

反照率的变化会直接改变气候系统的能量收支，进而触发反馈效应。最典型的例子是北极海冰反照率反馈：北极夏季海冰的反照率约为0.6-0.7，而融化后的开阔海水反照率仅为0.1，当全球变暖导致北极海冰融化时，更多的太阳短波辐射被海水吸收，进一步加剧海冰融化，形成正反馈循环。量化来看，北极夏季海冰范围每减少100万平方公里，对应的反照率降低约0.03，导致的辐射强迫约为0.5W/m²，这一反馈强度可使北极地区的升温速率达到全球平均的2-3倍，即“北极放大效应”。

除海冰外，土地利用变化也会显著改变地表反照率：热带雨林砍伐后，反照率从0.12升高到0.18，会增加短波辐射反射，产生约-1.5W/m²的负辐射强迫，但同时会释放储存的碳，产生约+2W/m²的正辐射强迫，净效应仍为增温；中高纬度地区的森林砍伐则会因为反照率升高的负强迫超过碳释放的正强迫，产生微弱的冷却效应。此外，积雪覆盖的时间变化也会触发反照率反馈：全球变暖导致积雪期缩短、积雪量减少，反照率降低，进一步加剧变暖，这一效应在高山地区和北半球中高纬度尤为显著。

气候系统净辐射收支方程与自然稳态

气候系统的能量收支可以通过净辐射收支方程定量描述。在大气层顶，能量收入为被吸收的短波辐射，即$(1-\alpha)S_0/4$；能量支出为地球系统向太空发射的长波辐射。在自然稳态下，能量收入与支出相等，即： $$\frac{(1-\alpha)S_0}{4} = OLR$$ 其中$OLR$为出射长波辐射（Outgoing Longwave Radiation）。地球作为近黑体，其出射长波辐射服从斯特藩-玻尔兹曼定律，即$OLR = \sigma T_e^4$，其中$T_e$为地球的有效辐射温度，是假设地球为均匀黑体时，能产生相同出射长波辐射的温度，$\sigma=5.67\times10^{-8}\ W/(m^2\cdot K^4)$为斯特藩-玻尔兹曼常数。

将自然状态下的参数代入方程，可得$(1-0.3)\times340 = 5.67\times10^{-8}\times T_e^4$，解得$T_e\approx255K$，即-18℃，这就是没有温室效应时地球的平均温度。而实际观测到的全球地表平均温度约为288K，即

实际观测到的全球地表平均温度约为288K，即15℃，比无温室效应下的有效辐射温度255K高出约33K，这33K的温差正是自然温室效应的核心贡献——大气中的温室气体分子吸收了地表发射的长波辐射，并通过逆辐射将部分能量返还地表，如同给地球覆盖了一层透射短波、阻挡长波的“保温膜”。理解这一效应的物理本质，需要从黑体辐射的基本定律和分子的量子能级结构入手。

温室效应的物理本质：从黑体辐射到分子能级跃迁

斯特藩-玻尔兹曼定律与地球长波辐射特征

任何温度高于绝对零度的物体都会向外发射电磁辐射，其辐射通量服从斯特藩-玻尔兹曼定律：$M = \varepsilon \sigma T^4$，其中$M$为辐射出射度，$\varepsilon$为发射率（黑体$\varepsilon=1$，实际地表的发射率约为0.95-0.98，接近黑体），$\sigma=5.67\times10^{-8}\ W/(m^2\cdot K^4)$为斯特藩-玻尔兹曼常数，$T$为物体绝对温度。对于温度为255K的地球有效辐射层，其出射长波辐射的峰值波长可通过维恩位移定律计算：$\lambda_{max}=2898\ \mu m\cdot K / T\approx11.4\ \mu m$，处于大气窗口波段（8-12μm），这也是卫星观测地球出射长波辐射的核心波段。

地表吸收短波辐射后升温，以长波辐射形式向外释放能量，其辐射峰值波长约为10μm，与大气层顶的出射长波峰值接近。如果没有大气，地表发射的长波辐射将直接逃逸到太空，地表温度将稳定在255K；但实际大气中存在能够吸收长波辐射的温室气体，它们吸收地表长波辐射后升温，同时向上下两个方向再辐射长波：向下的部分被地表吸收，进一步加热地表，向上的部分最终逃逸到太空，这就是温室效应的核心物理过程。

分子振动-转动能级跃迁与温室气体的选择吸收特性

温室气体之所以能够吸收长波辐射，本质是其分子结构的量子能级与长波辐射的能量匹配。分子的总能量由四部分组成：电子能级（能量差1-10eV，对应紫外-可见光波段）、振动能级（能量差0.05-1eV，对应近红外-中红外波段）、转动能级（能量差0.0001-0.05eV，对应远红外-微波波段）和平动能级（连续分布，对应所有波段）。长波辐射的光子能量恰好与分子的振动和转动能级跃迁能量匹配，因此可以被温室气体分子吸收。

并非所有分子都能成为温室气体，其核心条件是：分子振动或转动过程中必须产生偶极矩的变化，才能与电磁波的电场分量耦合，实现能量的吸收。以大气中占比99%的氮气（N₂）和氧气（O₂）为例，二者均为同核双原子分子，无永久偶极矩，且振动过程中偶极矩始终保持为0，因此无法吸收长波辐射，不属于温室气体。而二氧化碳（CO₂）、甲烷（CH₄）、氧化亚氮（N₂O）等微量气体均为不对称分子，振动过程中会产生偶极矩变化，因此具有强长波辐射吸收能力：

• CO₂为线性对称分子，无永久偶极矩，但其不对称伸缩振动（4.3μm吸收峰）和弯曲振动（15μm吸收峰）过程中会产生偶极矩变化，其中15μm吸收峰恰好覆盖地球出射长波辐射的峰值波段，是CO₂增温效应的核心贡献波段；
• CH₄为正四面体结构，振动过程中偶极矩变化显著，核心吸收峰位于3.3μm和7.7μm，其中7.7μm波段处于地球出射长波辐射的高能量区域；
• N₂O为线性不对称分子，核心吸收峰位于4.5μm、7.8μm和17μm，17μm波段同样处于长波辐射的高能量区域。

<!-- 电子版侧边栏：温室气体的“窗口”利用策略 --> > 温室气体的增温效率不仅取决于吸收强度，还取决于其吸收波段是否处于大气窗口之外：大气窗口是大气吸收最弱的8-12μm波段，地表发射的长波辐射可以几乎无衰减地逃逸到太空。如果温室气体的吸收峰位于大气窗口内，哪怕浓度很低，也能显著减少出射长波辐射，产生强增温效应。这也是为什么氟利昂等人工合成温室气体虽然浓度极低，但GWP极高的原因：其吸收峰恰好位于大气窗口内。

辐射强迫的定义与计算

辐射强迫（Radiative Forcing, RF）是IPCC评估气候变化的核心指标，定义为：由于外部驱动因子（温室气体浓度变化、气溶胶、太阳活动、土地利用变化等）引起的、大气层顶净辐射通量的变化，单位为$W/m^2$。正的辐射强迫表示气候系统能量收入增加，将导致全球平均温度上升；负的辐射强迫则表示能量支出增加，将导致温度下降。辐射强迫的概念剔除了气候系统内部的快速调整过程（如气温变化导致的水汽含量变化），仅反映外部驱动因子对能量收支的直接干扰，因此具有明确的物理意义和可观测性。

对于混合均匀的长寿命温室气体（Well-Mixed Greenhouse Gases, WMGHGs，如CO₂、CH₄、N₂O），IPCC第六次评估报告（AR6）给出了经过验证的经验计算公式： 1. CO₂辐射强迫：$$\Delta F_{CO_2} = 5.35 \times \ln\left(\frac{C}{C_0}\right)$$ 其中$C$为当前CO₂浓度，$C_0$为工业革命前（1850年）的CO₂浓度（约280ppm），该公式在180-2000ppm浓度范围内误差小于1%。 2. CH₄辐射强迫：由于CH₄会与大气中的羟基自由基（OH）发生反应，同时其吸收波段与N₂O存在重叠，因此其辐射强迫计算需引入修正项： $$\Delta F_{CH_4} = 0.036 \times \left(\sqrt{M} - \sqrt{M_0}\right) - \left(f(M, N_2O) - f(M_0, N_{2O,0})\right)$$ 其中$M$为当前CH₄浓度，$M_0$为工业革命前浓度（约722ppb），$N_2O$为当前N₂O浓度，$N_{2O,0}$为工业革命前浓度（约271ppb），$f$为重叠吸收修正函数。 3. N₂O辐射强迫：同样需考虑与CH₄的重叠修正： $$\Delta F_{N_2O} = 0.12 \times \left(\sqrt{N} - \sqrt{N_0}\right) - \left(f(M, N) - f(M_0, N_0)\right)$$

IPCC AR6进一步提出了有效辐射强迫（ERF）的概念，用于描述考虑大气内部快速调整后的净辐射强迫：对于WMGHGs而言，其大气调整时间尺度仅为几天到几周，因此RF与ERF的差异可忽略不计；但对于气溶胶、土地利用变化等分布不均的驱动因子，ERF能更准确地反映其对气候系统的实际影响。

主要温室气体的辐射强迫效率与百年尺度全球增温潜势

辐射强迫效率

辐射强迫效率是指温室气体浓度每增加1倍时，产生的辐射强迫变化，单位为$W/(m^2\cdot \text{倍数})$。不同温室气体的辐射强迫效率差异极大：CO₂浓度翻倍（从280ppm到560ppm）的辐射强迫约为3.7$W/m^2$，是基准值；CH₄浓度翻倍的辐射强迫约为0.6$W/m^2$，N₂O浓度翻倍的辐射强迫约为0.9$W/m^2$；而部分人工合成温室气体如六氟化硫（SF₆）的辐射强迫效率是CO₂的23900倍，尽管其2023年浓度仅为11ppb，但辐射强迫已超过0.5$W/m^2$。

百年尺度全球增温潜势（GWP）

全球增温潜势（Global Warming Potential, GWP）是《联合国气候变化框架公约》和《巴黎协定》下核算温室气体排放的核心指标，定义为：某一温室气体在特定时间范围内（通常为100年）的累积辐射强迫，相对于等质量CO₂累积辐射强迫的比值，无量纲。其数学表达式为： $$GWP_i(T) = \frac{\int_0^T HRF_i(t) dt}{\int_0^T HRF_{CO_2}(t) dt}$$ 其中$HRF_i(t)$为第$i$种温室气体的瞬态辐射强迫函数，$T$为时间尺度。

GWP的取值完全取决于温室气体的大气寿命和辐射强迫效率：寿命越长，在时间尺度内的累积辐射强迫越高；辐射强迫效率越高，单位浓度的增温效应越强。IPCC AR6给出的主要温室气体100年尺度GWP如下表所示：

| 温室气体 | 大气寿命 | 100年GWP | 20年GWP | 辐射强迫效率（相对于CO₂） | |----------|----------|----------|---------|----------------------------| | CO₂ | 混合（数十年至数千年） | 1（基准） | 1（基准） | 1 | | CH₄ | 12.4年 | 27-30 | 82-87 | 0.6 | | N₂O | 109年 | 273 | 273 | 0.9 | | SF₆ | 3200年 | 25200 | 37700 | 23900 |

<!-- 纸书：表1.1 主要温室气体的辐射特性与增温潜势（IPCC AR6） --> <!-- 标签：GWP，辐射强迫效率，温室气体寿命 -->

GWP的局限性需要特别说明：首先，GWP仅反映辐射强迫的累积效应，未考虑不同气体的区域分布差异、对碳循环和其他气候系统的反馈效应；其次，GWP的取值高度依赖时间尺度，例如CH₄的20年GWP是100年GWP的3倍，这是因为CH₄大气寿命短，短期增温效应远强于长期；第三，GWP未考虑不同温室气体对气候系统的非辐射效应，例如CH₄会增加对流层臭氧浓度，产生额外的增温效应，这部分在传统GWP计算中未被纳入。

工业革命以来人为温室气体排放的能量平衡干扰量化

温室气体浓度的观测变化

自工业革命以来，人类活动通过化石能源燃烧、工业过程、土地利用变化、农业活动等途径，向大气中排放了大量温室气体。根据美国NOAA莫纳罗亚天文台（Mauna Loa Observatory）自1958年以来的连续观测，大气CO₂浓度从1958年的315ppm上升至2024年5月的425ppm，年均增长率约2.5ppm/年，过去10年的增长率已超过3ppm/年，为过去2000万年以来的最高水平。CH₄浓度从工业革命前的722ppb上升至2023年的约1930ppb，N₂O浓度从271ppb上升至约336ppb，均远超过去80万年的冰芯记录的自然波动范围（CO₂在180-300ppm之间波动，CH₄在300-800ppb之间波动，N₂O在200-300ppb之间波动）。

人为辐射强迫的量化计算

将观测到的温室气体浓度代入辐射强迫公式，可计算得到1850-2019年主要温室气体的辐射强迫贡献：

• CO₂的辐射强迫约为2.16$W/m^2$，占总人为温室气体强迫的约60%；
• CH₄的辐射强迫约为0.54$W/m^2$，占比约15%；
• N₂O的辐射强迫约为0.21$W/m^2$，占比约6%；
• 其他长寿命温室气体（如氢氟碳化物HFCs、全氟化碳PFCs等）的辐射强迫约为0.41$W/m^2$，占比约12%。

除温室气体外，人类活动排放的气溶胶、黑碳、土地利用变化等因子也会产生辐射强迫：硫酸盐气溶胶通过反射短波辐射产生约-1.3$W/m^2$的负强迫，是有记录以来最大的人为冷却因子；黑碳通过吸收短波辐射和降低冰雪反照率产生约+0.51$W/m^2$的正强迫；土地利用变化（如毁林、城市化）的直接反照率效应产生约-0.15$W/m^2$的负强迫，但间接的碳排放效应产生约+0.56$W/m^2$的正强迫，净效应为+0.41$W/m^2$。

综合所有驱动因子，1850-2019年全球净人为辐射强迫约为2.72±0.3$W/m^2$，其中温室气体的正强迫完全抵消了气溶胶等因子的负强迫，剩余的正强迫是当前全球变暖的核心驱动。

观测到的能量平衡响应

根据卫星观测，1979年以来大气层顶的出射长波辐射已减少约1.5$W/m^2$，恰好与同期温室气体增加导致的辐射强迫增量一致，这是温室效应增强的直接观测证据。同时，海洋作为气候系统最大的热库，已吸收了约90%的多余能量，1990-2020年海洋上层2000米的热含量每十年上升约0.9×10²²J，相当于每秒爆炸5颗广岛原子弹释放的能量。由于海洋的热惯性，当前全球平均温度仅响应了约70%的净辐射强迫增量，剩余30%的能量仍被海洋吸收，即使温室气体浓度立即稳定，温度仍会继续上升0.3-0.5℃，这一现象被称为“剩余变暖”。

零维能量平衡模型的计算实践案例

零维能量平衡模型（Zero-Dimensional Energy Balance Model, EBM）是理解气候系统能量平衡的最简化模型，其假设地球为均匀的黑体，无大气和海洋环流的热输送，反照率固定，仅考虑大气层顶的短波入射和长波出射。尽管模型极度简化，但其模拟的平衡温度变化与复杂气候模式的结果高度一致，是理解气候敏感度的核心工具。

模型基本方程

零维EBM的基本方程为： $$C \frac{dT}{dt} = \frac{(1-\alpha)S_0}{4} + \Delta F(t) - \varepsilon \sigma T^4 - \lambda (T - T_0)$$ 其中：

• $C$为系统的有效热容量，单位$J/(m^2\cdot K)$，对于包含上层海洋的地球系统，$C\app