地球辐射平衡与温室效应物理基础

当太阳的短波辐射穿越近1.5亿公里的太空抵达地球大气层顶端时，一个决定行星命运的物理过程随即启动。这些辐射中，约30%被云层、气溶胶及地表反射回太空，这一比例即行星反照率（α）；剩余70%被地表和大气吸收。这些被吸收的能量最终以长波红外辐射形式重新发射。在稳态气候系统中，入射的太阳辐射能量必须等于出射至太空的红外辐射能量，这一能量平衡方程决定了地球的平均温度。然而，这个看似简单的能量守恒关系背后，隐藏着使地球平均温度比无大气计算值高出约33摄氏度的关键机制——温室效应。理解现代气候变化，必须首先掌握这一最基础的行星能量平衡原理及其数学表达，进而深入温室气体如何通过分子层面的红外吸收与再发射，扰动这一精妙平衡。本章将从第一性原理出发，系统构建地球辐射平衡的理论框架，推导核心定律，引入辐射强迫概念，并揭示温室效应的物理本质。

行星能量平衡：从太阳常数到有效温度

地球气候系统的能量来源几乎完全依赖于太阳。在地球大气层顶端，垂直于太阳光线的单位面积上接收到的太阳辐射通量被称为太阳常数（S）。尽管名称中有“常数”二字，其值实际上存在微小变化（主要与地球轨道参数和太阳活动周期相关），目前公认的长期平均值约为 1361 ± 0.5 W/m²。这一数值是通过卫星精确测量得到的基准。

然而，地球是一个球体，接收太阳辐射的面积是圆形的横截面积（πR²），而辐射和重新发射发生的面积是整个球面（4πR²）。因此，全球平均的入射太阳辐射通量为 S/4。设行星反照率为 α，则被地球系统吸收的净太阳辐射通量（即吸收的短波辐射）为：

$$ F_{ absorbed } = \frac{S(1 - \alpha)}{4} $$

在稳态条件下，地球系统必须通过向太空发射长波红外辐射来平衡这一吸收的能量。根据Stefan-Boltzmann定律，一个黑体（完美吸收和发射辐射的物体）单位面积向所有方向发射的总辐射通量（辐射出射度）与其绝对温度的四次方成正比：

$$ F = \sigma T^4 $$

其中 σ 是 Stefan-Boltzmann 常数，其值为 5.670374419 × 10⁻⁸ W·m⁻²·K⁻⁴（2019年重新定义后的精确值）。如果地球是一个没有大气的均匀黑体，其表面温度 Tₑ 应满足吸收的太阳辐射等于发射的红外辐射：

$$ \frac{S(1 - \alpha)}{4} = \sigma T_e^4 $$

解出有效温度：

$$ T_e = \left( \frac{S(1 - \alpha)}{4\sigma} \right)^{1/4} $$

代入 S = 1361 W/m²，α = 0.3（地球当前平均反照率），计算得：

$$ T_e = \left( \frac{1361 \times (1 - 0.3)}{4 \times 5.67 \times 10^{-8}} \right)^{1/4} \approx \left( \frac{952.7}{2.268 \times 10^{-7}} \right)^{1/4} \approx (4.20 \times 10^9)^{1/4} \approx 255 \text{ K} $$

255 K（约 -18°C）是地球作为无大气黑体的理论有效辐射温度。然而，我们观测到的地球平均表面温度约为 288 K（15°C），高出约33 K。这33 K的温差正是自然温室效应的贡献。月球可作为对照：其反照率约0.12，无大气，计算有效温度约270 K，但实际表面温度昼夜极差极大（从约100 K到400 K），平均约220 K，低于计算值，部分因其热惯性低，部分因表面热辐射特性非理想黑体。地球的稳定高温差则明确指向大气层对出射红外辐射的阻挡作用。

Stefan-Boltzmann定律的物理起源与推导

Stefan-Boltzmann定律并非凭空而来，它根植于19世纪对热辐射本质的探索，最终由玻尔兹曼从热力学理论严格推导。其物理本质是：一个处于热平衡的腔体，其内壁单位面积单位时间发射的辐射能量与绝对温度的四次方成正比。推导可简述如下：

考虑一个温度为 T 的空腔，内壁为黑体。腔内辐射场具有各向同性。根据热力学第二定律，若将空腔分为两个温度分别为 T 和 T+dT 的部分，净能量流必须为零，这导致辐射能量密度 u 与 T⁴ 成正比。而辐射通量 F = (c/4) u（c为光速，1/4因子源于角度积分），故 F ∝ T⁴。比例常数 σ 由基本物理常数决定，其值可通过实验精确测定。

对于非黑体，需引入发射率 ε（0 ≤ ε ≤ 1），实际辐射通量为 F = εσT⁴。地球表面和大气并非完美黑体，尤其在特定红外波段，发射率可能显著低于1，这为温室气体选择性地吸收和发射辐射提供了物理基础。

辐射强迫：量化扰动的核心指标

当行星能量平衡因外部因素而改变时，我们使用辐射强迫（Radiative Forcing, RF）来量化这种扰动。根据政府间气候变化专门委员会（IPCC）的权威定义，辐射强迫是指“由于外部因素（如大气成分变化、太阳辐射变化）改变而使大气层顶（TOA）净（向下减向上）辐射通量发生的改变”，所有其他参数（如温度、云量）保持初始值不变。其单位是 W/m²，表示全球平均的辐射通量变化。

辐射强迫是连接外部驱动因子与气候系统初始响应的关键桥梁。正的辐射强迫（如CO₂增加）导致系统净得热，倾向于升温；负的强迫（如火山气溶胶增加反射）导致净失热，倾向于降温。IPCC第六次评估报告（AR6，2021）给出了1750年至2019年主要人为强迫因子的最佳估计值：

| 强迫因子 | 辐射强迫 (W/m²) | 主要机制 | |----------|----------------|----------| | 温室气体 (CO₂, CH₄, N₂O等) | +2.72 | 吸收地表红外辐射 | | 气溶胶（直接效应） | -0.45 | 散射太阳辐射 | | 气溶胶（云反照率效应） | -0.84 | 增加云反射率 | | 土地利用变化 | -0.03 | 改变地表反照率 | | 太阳活动变化 | +0.05 | 改变太阳辐照度 | | 总人为强迫 | +2.72 | — |

其中，CO₂贡献最大（约+2.16 W/m²），其次是甲烷（CH₄，约+0.54 W/m²）。值得注意的是，总强迫的不确定性主要来自气溶胶，其冷却效应部分抵消了温室气体的增暖效应。

关键洞察：辐射强迫是一个理论概念，它假设除该因子外气候系统其他一切保持不变（“固定一切”假设）。在实际气候响应中，温度变化会引发水汽、云、冰盖等反馈，这些反馈会放大或减小初始强迫，最终的气候变化幅度远大于辐射强迫单独作用的结果。这一区分是理解气候敏感性的基础。

温室效应的微观物理机制：分子振动与红外吸收

温室效应的核心在于大气中某些气体分子（温室气体）能够吸收地表发射的长波红外辐射。这并非所有气体都能做到，其能力根植于分子的电偶极矩变化。

分子振动模式与红外活性

多原子分子（如H₂O、CO₂、CH₄、N₂O）由原子通过化学键连接，可以发生对称伸缩、不对称伸缩、弯曲等振动模式。根据量子力学，这些振动具有特定的能级。当红外光子的能量（E = hν，h为普朗克常数，ν为频率）恰好等于分子两个振动能级之差时，分子会吸收该光子，从低能级跃迁到高能级。

关键条件：吸收必须导致分子电偶极矩的改变。对称分子如N₂、O₂（地球大气主要成分）在振动时电偶极矩不变，因此几乎不吸收红外辐射，它们不是温室气体。而CO₂（线性对称分子，但弯曲振动破坏对称性）、H₂O（永久电偶极矩）、CH₄（非对称振动）等，其振动模式会改变电偶极矩，因此是强红外吸收体。

吸收光谱与大气窗口

每种温室气体都有其特征吸收光谱，即吸收系数随波长（或波数）变化的曲线。下图（概念示意）展示了主要温室气体在红外波段（4-25 μm，即2500-400 cm⁻¹）的吸收特性：

```mermaid graph LR A[红外辐射频谱] --> B{大气吸收}; B --> C[H₂O: 宽吸收带<br>5-8 μm, >18 μm]; B --> D[CO₂: 强吸收带<br>4.3 μm, 15 μm]; B --> E[CH₄: 吸收带<br>3.3 μm, 7.7 μm]; B --> F[N₂O: 吸收带<br>4.5 μm, 7.8 μm, 17 μm]; C --> G[红外窗口<br>8-13 μm]; D --> G; E --> G; F --> G; G --> H[地表辐射可部分逃逸]; ```

红外窗口：在大气吸收光谱中，存在一些相对较弱的区域，最著名的是 8-13 μm（约1250-770 cm⁻¹） 波段。在此窗口，大气吸收较弱，地表和低层大气的热辐射可以较直接地逃逸至太空。这是地球散热的主要通道。温室气体浓度的增加，尤其是那些在窗口附近有吸收翼（wings）的气体（如CO₂在15 μm吸收带两侧的翼），会“堵塞”这一窗口，减少出射辐射，从而产生增暖效应。

吸收饱和与“ wings”效应

一个常见误解是：CO₂在15 μm附近已有很强吸收，增加浓度不会产生额外效果，因为吸收已“饱和”。这需要更精细的辐射传输理论来澄清。

在吸收带中心，由于吸收极强，来自地表的红外辐射确实在低层大气即被完全吸收，增加CO₂浓度对中心波段的吸收影响甚微。然而，在吸收带的两侧（即“wings”区域），吸收系数随浓度增加而线性增加，这些原本未饱和的波段因此被加强吸收。更重要的是，有效辐射高度（即大气变得透明、辐射得以逃逸至太空的高度）在吸收带中心位于低温的高层大气，而在吸收带翼部，有效辐射高度位于相对较暖的低层大气。当CO₂增加时，吸收带翼部增强，使得这些波段的有效辐射高度上升至更冷的高层，根据Stefan-Boltzmann定律（F = σT⁴），高层大气温度更低，发射的辐射通量更小，导致大气层顶净辐射强迫为正。因此，温室效应的增强主要发生在吸收光谱的翼部，而非中心。这一机制是定量计算辐射强迫的基础。

辐射传输简化模型：从单层到多层

为了直观理解温室效应，我们构建理想化的辐射传输模型。

单层大气模型（零维能量平衡模型）

假设： 1. 地表为黑体，温度 Tₛ。 2. 大气层为单层，位于地表之上，完全透明于太阳短波辐射（τₛ = 1），但对地表发射的长波辐射具有均匀的吸收率 λ（0 < λ < 1），且自身发射黑体辐射（发射率等于吸收率，由Kirchhoff定律）。 3. 大气层顶处于稳态，无其他能量源。 4. 所有辐射垂直传输。

设太阳入射为 S(1-α)/4。能量平衡方程如下：

大气层顶平衡：入射太阳辐射 = 地表发射穿过大气层顶的辐射 + 大气层发射至太空的辐射。 但大气层顶的向上辐射来自两部分：① 地表发射未被大气吸收的部分（1-λ）σTₛ⁴；② 大气层自身发射（向上部分）λσTₐ⁴（假设大气双向等量发射）。故： $$ \frac{S(1-\alpha)}{4} = (1-\lambda)\sigma T_s^4 + \lambda\sigma T_a^4 \quad (1) $$

大气层能量平衡：大气吸收的能量 = 它发射的能量。大气吸收来自地表发射的 λσTₛ⁴。它向上下各发射 λσTₐ⁴（总发射 2λσTₐ⁴）。故： $$ \lambda\sigma T_s^4 = 2\lambda\sigma T_a^4 \quad (2) $$

若 λ ≠ 0，由(2)得： $$ T_a^4 = \frac{T_s^4}{2} \quad \Rightarrow \quad T_a = 2^{-1/4} T_s \approx 0.8409 T_s $$

代入(1)： $$ \frac{S(1-\alpha)}{4} = (1-\lambda)\sigma T_s^4 + \lambda\sigma (T_s^4/2) = \sigma T_s^4 \left(1 - \frac{\lambda}{2}\right) $$

解得地表温度： $$ T_s = \left[ \frac{S(1-\alpha)}{4\sigma (1 - \lambda/2)} \right]^{1/4} $$

与无大气（λ=0）的有效温度 Tₑ = [S(1-α)/(4σ)]^{1/4} 相比： $$ T_s = \frac{T_e}{(1 - \lambda/2)^{1/4}} $$

当 λ > 0 时，分母小于1，Tₛ > Tₑ。例如 λ=0.78（近似大气总红外吸收率），则 (1-0.39)=0.61，0.61^{-1/4}≈1.15，Tₛ≈1.15×255 K≈293 K，接近观测值。该模型虽极度简化，但清晰展示了：大气对地表红外辐射的部分吸收，使得大气自身发射红外辐射，其中一部分返回地表（“大气逆辐射”），从而加热地表，使其温度高于无大气时的有效温度。

多层模型与饱和吸收

单层模型无法解释吸收带饱和问题。引入多层模型：设大气分为 N 层，每层有相同吸收率 λ（或更真实地，吸收系数随高度指数衰减）。可以证明，随着层数增加，地表温度升高，但每增加一层，温升幅度递减。当 λ 接近1（强吸收）时，即使增加浓度，只要吸收带中心已饱和，单层模型的结论（Tₛ 依赖于 λ）似乎会失效。然而，实际大气吸收并非在整个波段均匀饱和，且光谱积分效应使得“翼部”吸收始终存在。更严谨的模型需考虑光谱分辨的辐射传输方程：

$$ \frac{dI_\nu}{d\tau_\nu} = I_\nu - B_\nu(T) $$

其中 I_ν 是光谱辐射强度，τ_ν 是光学厚度（与吸收气体柱浓度成正比），B_ν(T) 是 Planck 黑体辐射函数。通过求解此方程并积分所有频率，可精确计算给定气体浓度变化下的辐射强迫。这正是 IPCC 报告和气候模型使用的标准方法。计算表明，当前CO₂浓度下，加倍（从 280 ppm 到 560 ppm）产生的全球平均辐射强迫约为 3.7 W/m²，常用近似公式为 ΔF = 5.35 ln(C/C₀) W/m²（适用于 C 在 200-1000 ppm 范围）。

光谱证据：卫星观测与吸收带识别

理论需由观测验证。卫星遥感，特别是红外光谱仪（如 NASA 的 AIRS、Aqua 卫星上的 CERES，以及更早的 Nimbus 系列），直接测量了地球大气层顶和大气内部的出射红外辐射光谱。

观测光谱清晰显示出： 1. 地表黑体辐射曲线：在无吸收的“窗口”波段（如 10-12 μm），观测到的辐射强度接近 300 K 黑体曲线，对应地表温度。 2. 吸收带特征：在 CO₂ 的 15 μm 带、H₂O 的 6.3 μm 和 >18 μm 带、O₃ 的 9.6 μm 带等位置，观测辐射远低于黑体曲线，表明这些波段被大气强烈吸收，有效辐射来源是较高层较冷的大气。 3. 时间变化趋势：对比 1970 年与 2020 年的卫星光谱，在 CO₂ 和 CH₄ 的吸收带中心，出射辐射强度略有下降（对应高层大气变冷？不，实际是辐射高度上升），而在吸收带翼部，下降更显著，这正是浓度增加导致辐射强迫的指纹。下图（基于 Myhre 等人研究的概念图）展示了 CO₂ 加倍前后理论出射光谱的变化：

```mermaid graph TD A[地表黑体辐射 ~300K] --> B[无CO₂大气: 光谱接近黑体曲线]; C[当前CO₂浓度] --> D[观测光谱: 15μm带深吸收]; D --> E[CO₂加倍后理论光谱]; E --> F[15μm带翼部进一步减弱<br>辐射强迫 ~3.7 W/m²]; F --> G[为恢复平衡, 地表与低层大气必须升温]; ```

这种光谱指纹的独特性和可检测性，是确认人为温室气体增加导致辐射强迫的决定性证据。它排除了太阳活动等外部因子作为主要驱动力的可能性，因为太阳变化会影响所有波段，而不会在特定气体吸收带产生特征性变化。

关键公式总结与能量平衡计算实践

本章核心公式构成了气候物理的基石：

1. Stefan-Boltzmann 定律：$$ F = \sigma T^4 $$ 描述黑体辐射出射度。 2. 行星有效温度：$$ T_e = \left( \frac{S(1-\alpha)}{4\sigma} \right)^{1/4} $$ 无大气时的平衡温度。 3. 辐射强迫（简化）：对于 CO₂，$$ \Delta F_{CO_2} \approx 5.35 \ln(C/C_0) \text{ (W/m}^2) $$。对于其他气体有类似公式或查表。 4. 单层大气模型地表温度：$$ T_s = \left[ \frac{S(1-\alpha)}{4\sigma (1 - \lambda/2)} \right]^{1/4} $$，λ 为大气红外吸收率。 5. 多层/光谱精确模型：需数值求解辐射传输方程：$$ \frac{dI_\nu}{d\tau_\nu} = I_\nu - B_\nu(T) $$，并垂直积分。

计算示例：工业革命以来CO₂的辐射强迫

已知：工业革命前（约1750年）CO₂浓度 C₀ = 280 ppm，2019年 C = 415 ppm。 $$ \Delta F = 5.35 \ln(415/280) = 5.35 \ln(1.482) = 5.35 \times 0.394 = 2.11 \text{ W/m}^2 $$ 这与 IPCC AR6 给出的 +2.16 W/m² 非常接近。若考虑所有温室气体（包括CH₄、N₂O、卤代烃等），总温室气体强迫约 +2.72 W/m²（2019年相对于1750年）。

计算示例：无温室效应的地球温度

取 S=1361 W/m², α=0.3, σ=5.67×10⁻⁸ W/m²K⁴。 $$ T_e = \left( \frac{1361 \times 0.7}{4 \times 5.67 \times 10^{-8}} \right)^{1/4} = \left( \frac{952.7}{2.268 \times 10^{-7}} \right)^{1/4} = (4.202 \times 10^9)^{1/4} $$ 计算：4.202e9 的四次方根。先平方根：√(4.202e9) ≈ 64800，再平方根：√64800 ≈ 254.6 K。 故 Tₑ ≈ 255 K (-18°C)。实际地表平均约 288 K (15°C)，自然温室效应贡献约 33 K。

###