地球辐射平衡与温室效应 / Earth's Radiative Balance and the Greenhouse Effect

地球气候系统的稳定性建立在全球辐射能量收支的精确平衡之上。太阳辐射作为几乎唯一的外部能量来源，以短波形式入射地球，而地球则通过发射长波红外辐射将能量返还太空。这一平衡的任何系统性偏差——无论是自然变异还是人为干扰——都会导致地球系统的能量积累或亏损，进而驱动全球温度变化、冰盖消融和海洋环流重组。理解这一平衡的物理机制，是量化人类活动对气候影响、预测未来气候 trajectories 并制定有效减缓策略的科学基石。本章将系统阐述地球辐射平衡的核心物理原理，从太阳辐射的光谱特征到温室气体的辐射传输过程，构建起从基础定律到定量预测的完整逻辑链条。

1.1 太阳辐射输入：光谱特性与行星尺度分配

太阳作为近似黑体，其辐射光谱遵循普朗克定律，峰值波长位于可见光波段（约0.5 μm），总辐射功率由斯蒂芬-玻尔兹曼定律描述。在地球轨道位置，太阳常数（Solar Constant） $S_0$ 定义为大气顶垂直于太阳光束的单位面积接收的辐射通量，其最新校准值约为1361 W/m²（Kopp & Lean, 2011）。由于地球球面几何，全球平均入射太阳辐射为 $S_0/4 \approx 340$ W/m²。

太阳辐射在大气顶的分布受日地距离、太阳活动周期（约0.1%变异）和轨道参数（米兰科维奇周期）调制，但后两者影响的时间尺度远超工业世纪，本章聚焦于稳定状态下的辐射收支。关键的是，太阳辐射并非全部被地球系统吸收：行星反照率（$\alpha$）描述了被云、气溶胶、冰雪和地表反射回太空的太阳辐射比例。根据CERES卫星观测（2000-2020年），全球平均行星反照率约为0.298（Loeb et al., 2020），意味着约101 W/m²的入射太阳辐射被反射，净吸收太阳辐射为 $F_{\text{sun}} = (S_0/4)(1-\alpha) \approx 239$ W/m²。这一数值必须与地球系统向太空发射的长波辐射精确平衡，以维持零辐射强迫状态。

1.2 地球长波辐射与黑体辐射定律

地球表面及大气在吸收太阳辐射后，以长波红外辐射（波长4-100 μm）形式冷却。若地球是一个均匀的绝对黑体（吸收率=发射率=1），根据斯蒂芬-玻尔兹曼定律，其发射的向上长波辐射通量应为 $F_{\text{out}} = \sigma T_e^4$，其中 $\sigma = 5.67 \times 10^{-8}$ W·m⁻²·K⁻⁴ 为斯蒂芬-玻尔兹曼常数，$T_e$ 为有效辐射温度。在零辐射强迫平衡下，$F_{\text{out}} = F_{\text{sun}}$，因此：

$$T_e = \left( \frac{F_{\text{sun}}}{\sigma} \right)^{1/4} = \left( \frac{(S_0/4)(1-\alpha)}{\sigma} \right)^{1/4}$$

代入数值计算得 $T_e \approx 255$ K（-18°C）。然而，全球平均地表温度 $T_s$ 约为288 K（15°C），大气低层温度也显著高于255 K。这一33 K的温差即温室效应的直接体现：大气对地表长波辐射的部分吸收和再发射，使得大气顶的等效辐射高度升高至更冷的大气层，为维持能量平衡，地表必须升温以增加其发射通量。

维恩位移定律 $\lambda_{\text{max}} T = b$（$b \approx 2898$ μm·K）进一步阐明了光谱差异：太阳辐射峰值在0.5 μm，而255 K黑体峰值在11.4 μm，288 K地表峰值在10.1 μm。大气中某些气体（如H₂O、CO₂、CH₄）在8-13 μm波段有强吸收带，恰好覆盖了地表辐射峰值，这是温室效应光谱选择性的物理基础。

1.3 辐射传输方程（RTE）的推导与求解

辐射传输方程是描述电磁辐射在介质中传播时强度变化的微分方程，其一般形式为：

$$\frac{dI_\nu}{ds} = -\kappa_\nu I_\nu + j_\nu$$

其中 $I_\nu$ 是光谱辐射强度（W·m⁻²·sr⁻¹·Hz⁻¹），$s$ 是传播路径，$\kappa_\nu$ 是质量吸收系数（m²/kg），$j_\nu$ 是质量发射系数（W·sr⁻¹·kg⁻¹·Hz⁻¹）。对于热平衡介质，基尔霍夫定律给出 $j_\nu = \kappa_\nu B_\nu(T)$，$B_\nu(T)$ 为普朗克黑体辐射谱强度：

$$B_\nu(T) = \frac{2h\nu^3}{c^2} \frac{1}{e^{h\nu/kT} - 1}$$

在平面平行大气假设下（垂直方向z，向上为正），忽略散射（对长波辐射近似合理），RTE简化为：

$$\mu \frac{dI_\nu(z)}{dz} = -\kappa_\nu(z) I_\nu(z) + \kappa_\nu(z) B_\nu(T(z))$$

其中 $\mu = \cos\theta$ 为天顶角余弦。定义光学厚度 $\tau_\nu(z) = \int_z^\infty \kappa_\nu(z') dz'$，方程可积分求解。对于上行辐射（$\mu>0$），边界条件为大气顶（$\tau_\nu=0$）无入射辐射，地表（$\tau_\nu=\tau_{\nu,0}$）处 $I_\nu(0) = \epsilon_\nu B_\nu(T_s)$（$\epsilon_\nu$ 为地表比辐射率，通常≈1）。解为：

$$I_\nu^+(0) = \int_0^{\tau_{\nu,0}} B_\nu(T(\tau_\nu)) e^{-\tau_\nu/\mu} \frac{d\tau_\nu}{\mu}$$

大气顶净通量需对立体角和频率积分：

$$F_{\text{TOA}} = 2\pi \int_0^\infty \int_0^1 I_\nu^+(0) \mu d\mu d\nu - F_{\text{sun}}(1-\alpha)$$

实际求解需处理连续光谱，常用方法有：

• **逐线计算（Line-by-Line, LBL）**：使用HITRAN等数据库，对每条吸收线精确积分，计算成本极高。
• **带模式（Band Model）**：将光谱离散化为若干“带”，每带内参数平均，如 correlated-k 方法。
• **简化两流近似**：将辐射场近似为向上和向下两个流，适用于快速估算。

1.4 温室效应的定量描述：从单层到连续大气

单层大气模型是理解温室效应的经典框架。假设大气为单层，完全吸收地表长波辐射（在吸收带内），且自身以黑体辐射 $B(T_a)$ 向上下发射。地表能量平衡：

$$\frac{S_0}{4}(1-\alpha) + \epsilon \sigma T_a^4 = \sigma T_s^4$$

大气能量平衡（假设大气无太阳辐射吸收）：

$$\epsilon \sigma T_s^4 = 2 \epsilon \sigma T_a^4$$

其中 $\epsilon$ 为大气在长波波段的吸收率（0<$\epsilon$<1）。联立解得：

$$T_s = \left( \frac{1}{1-\epsilon/2} \cdot \frac{(S_0/4)(1-\alpha)}{\sigma} \right)^{1/4} = T_e \left(1-\frac{\epsilon}{2}\right)^{-1/4}$$

当 $\epsilon=1$（完全温室），$T_s = 2^{1/4} T_e \approx 303$ K，但实际大气非完全吸收（有“大气窗”），且吸收率随气体浓度变化，故 $T_s$ 介于 $T_e$ 与303 K之间。此模型虽粗糙，但清晰展示了温室气体增加 $\epsilon$ 将升高 $T_s$。

连续大气模型更接近真实。考虑大气垂直分层，每层有温度 $T(z)$ 和吸收系数 $\kappa_\nu(z)$。温室效应增强的本质是：增加CO₂等气体浓度，提高大气在吸收带的光学厚度，使有效辐射高度上升。由于温度随高度递减（对流层平均递减率6.5 K/km），辐射高度升高意味着该处温度更低，为保持大气顶出射辐射通量不变（短期近似），地表必须升温以补偿。这一机制可通过辐射强迫概念量化。

1.5 辐射强迫与气候敏感度

辐射强迫（Radiative Forcing, RF） 定义为：在给定气候状态下，某因子（如CO₂浓度变化）改变后，大气顶净（太阳减去长波）辐射通量的变化，所有其他参数（如温度、水汽、云）保持不变的条件下（IPCC, 2021）。单位：W/m²。正强迫导致系统能量积累。

对于CO₂，RF可通过简化公式估算（Myhre et al., 1998）：

$$\Delta F_{\text{CO}_2} = 5.35 \ln\left( \frac{C}{C_0} \right)$$

其中 $C$ 和 $C_0$ 分别为扰动前后CO₂浓度（ppm）。例如，从280 ppm（工业前）到415 ppm（2020年），$\Delta F \approx 5.35 \ln(415/280) \approx 2.0$ W/m²。该公式基于LBLRTM计算，适用于1-10倍变化，误差约±10%。

其他温室气体（CH₄、N₂O、卤代烃）的RF计算需考虑吸收带重叠，通常使用叠加公式或更复杂的参数化。气溶胶的RF为负值（反射太阳辐射），但不确定性大。

气候敏感度描述全球平均地表温度对辐射强迫的响应。平衡气候敏感度（Equilibrium Climate Sensitivity, ECS） 定义为：大气CO₂浓度加倍（$\Delta F_{2\times\text{CO}_2} \approx 3.7$ W/m²）后，全球平均地表温度达到新平衡时的变化量 $\Delta T_{2\times}$。根据物理原理和模型，ECS likely range为2.0-4.5°C（IPCC AR6），中位数约3.0°C。瞬时气候响应（TCR）则指CO₂以每年1%速率加倍时，20年时的温度变化，通常小于ECS（约1.5-2.5°C），因海洋热惯性未完全响应。

关系可简化为：$\Delta T_s = \lambda \Delta F$，其中 $\lambda$ 为气候反馈参数（W/m²/°C），$\lambda = 1/\text{ECS} \times \Delta F_{2\times}$。负反馈（如长波辐射增加）使 $\lambda>0$，但水汽、冰反照率、云等反馈放大初始变暖，使ECS远大于无反馈响应（约1.2°C）。

1.6 辐射传输代码示例：简化一维模型

以下Python代码实现一个两流近似的一维辐射传输模型，用于计算CO₂浓度变化引起的大气顶辐射强迫。模型假设：

• 大气为单层，温度 $T_a$ 由能量平衡确定。
• 地表温度 $T_s$ 输入（或迭代求解）。
• 大气在长波波段有吸收率 $\epsilon$，与CO₂浓度相关。
• 太阳辐射仅考虑净吸收，忽略大气吸收。

```python import numpy as np import matplotlib.pyplot as plt

常数

S0 = 1361.0 # 太阳常数 (W/m²) alpha = 0.298 # 行星反照率 sigma = 5.67e-8 # Stefan-Boltzmann 常数 (W/m²/K⁴) F_sun = (S0/4) * (1 - alpha) # 净吸收太阳辐射 (W/m²)

def epsilon_from_co2(co2_ppm, eps0=0.78, co2_ref=415): """ 简化：吸收率与CO₂浓度对数成正比，参考值 eps0 对应 co2_ref。 实际关系更复杂，此处为教学简化。 """ return eps0 + 0.1 * np.log(co2 / co2_ref) # 系数0.1为示意

def radiative_forcing(co2_ppm, Ts=288.0): """ 计算给定CO₂浓度和地表温度下的大气顶净辐射通量。 返回: F_toa (W/m²) """ eps = epsilon_from_co2(co2_ppm)

单层大气能量平衡: F_sun + eps * sigma * Ta**4 = sigma * Ts**4

大气平衡: eps * sigma * Ts**4 = 2 * eps * sigma * Ta**4 => Ta = Ts / 2**0.25

Ta = Ts / (2**0.25)

大气顶向上长波辐射: 大气向上发射 + 透过大气的表面辐射

F_out = eps * sigma * Ta4 + (1 - eps) * sigma * Ts4 F_toa = F_sun - F_out return F_toa

计算辐射强迫：从工业前280 ppm到当前415 ppm

co2_hist = 280 co2_now = 415 F_hist = radiative_forcing(co2_hist) F_now = radiative_forcing(co2_now) RF = F_now - F_hist print(f"工业前({co2_hist} ppm)大气顶净辐射: {F_hist:.2f} W/m²") print(f"当前({co2_now} ppm)大气顶净辐射: {F_now:.2f} W/m²") print(f"CO₂引起的辐射强迫: {RF:.2f} W/m²")

计算不同CO₂浓度下的强迫，绘制曲线

co2_range = np.linspace(200, 1000, 100) RF_range = [radiative_forcing(c) - radiative_forcing(280) for c in co2_range]

plt.figure(figsize=(8,5)) plt.plot(co2_range, RF_range, 'b-') plt.axvline(x=415, color='r', linestyle='--', label='当前水平 (415 ppm)') plt.axhline(y=3.7, color='g', linestyle=':', label='CO₂加倍强迫 (~3.7 W/m²)') plt.xlabel('CO₂浓度 (ppm)') plt.ylabel('辐射强迫 (W/m²)') plt.title('CO₂浓度变化与辐射强迫（简化模型）') plt.legend() plt.grid(True) plt.show() ```

此代码输出当前CO₂增加约2.0 W/m²（与公式5.35 ln(415/280)=2.0一致），并显示强迫与对数关系。注意：实际RF计算需考虑光谱重叠、非瞬时混合等，此处仅为教学示意。更真实的模型需使用LBLRTM或带模式代码。

1.7 卫星观测验证：CERES辐射不平衡测量

卫星观测为辐射平衡理论提供了直接验证。NASA的CERES（Clouds and the Earth's Radiant Energy System）仪器自2000年起持续测量大气顶的短波（0.3-5 μm）和长波（5-200 μm）辐射通量。关键发现包括：

1. 平均辐射不平衡：2000-2020年间，CERES测得全球平均净辐射流入为 $+0.9 \pm 0.15$ W/m²（Loeb et al., 2021），正值表示系统能量积累，与海洋热含量上升（约0.5-1.0 W/m²）一致，证实正辐射强迫存在。

2. 光谱验证：CERES的长波波段观测显示，在CO₂吸收带（如15 μm）出射辐射减少，而大气窗（8-12 μm）变化较小，与LBL模型预测的CO₂增加导致的辐射冷却高度下降吻合。

3. 强迫分解：通过“能量约束”方法，将观测到的净辐射变化分解为各强迫因子（温室气体、气溶胶、云、反照率变化）的贡献。例如，温室气体增加贡献约+2.5 W/m²（1990s至2010s），气溶胶变化约-0.5 W/m²，云变化约-0.2 W/m²，总和与观测净强迫匹配（Smith et al., 2021）。

4. 空间模式：辐射强迫的空间分布与温室气体增加的理论模式一致：热带海洋地区净辐射增加最显著，极地放大效应在长波辐射中可见，但受云反馈调制。

这些观测不仅验证了辐射传输理论的正确性，还揭示了云和气溶胶反馈的复杂性，为气候模型提供关键约束。

1.8 模型局限与前沿挑战

尽管辐射平衡框架坚实，但实际应用面临多重挑战：

• **光谱精度**：LBL计算需数小时处理单天光谱，业务模型依赖带模式，吸收线参数（如HITRAN）仍有不确定性，尤其对微量气体。
• **垂直分辨率**：大气垂直结构（特别是温度廓线）对辐射通量敏感，但气候模型垂直层有限（通常30-100层），可能平滑关键逆温层。
• **地表异质性**：地表温度、反照率、比辐射率空间变化大，全球平均引入误差。
• **云-辐射相互作用**：云同时影响短波反射和长波吸收，其微物理（液水路径、粒子大小）对辐射强迫高度敏感，但云参数化是气候模型最大不确定性来源。
• **瞬变响应**：辐射强迫定义假设其他参数不变，但实际水汽、云等会快速调整（“快速调整”），产生“有效辐射强迫”（ERF），更接近真实气候响应。例如，CO₂增加直接导致平流层冷却，对流层水汽增加，这些调整需在完整大气-海洋耦合模型中计算。

当前研究前沿包括：利用高光谱卫星（如CLARREO）提升观测精度；发展机器学习加速LBL计算；通过“辐射 kernels”方法量化反馈贡献；以及将辐射传输模块与云微物理更紧密耦合。

1.9 本章总结

本章从基础物理定律出发，系统构建了地球辐射平衡的理论框架。太阳辐射输入经行星反照率调制后，净吸收约239 W/m²；地球通过长波辐射出射，若为黑体则有效温度255 K，但实际地表温度288 K，差值源于大气对长波辐射的吸收再发射——即温室效应。辐射传输方程定量描述了这一过程，其求解依赖于光谱吸收特性。温室气体浓度增加通过提高大气光学厚度，抬升有效辐射高度，产生正辐射强迫（如CO₂加倍约3.7 W/m²）。气候敏感度（ECS）量化了强迫向温度的转换，当前最佳估计为2.0-4.5°C，其中反馈过程（水汽、冰、云）起放大作用。简化代码示例演示了强迫的对数依赖性，而CERES卫星观测直接验证了辐射不平衡的存在及其光谱特征，为理论提供坚实证据。

然而，辐射平衡仅是气候系统的起点。下一章将深入探讨大气和海洋