气候系统的能量平衡

气候系统的能量平衡

引言

能量守恒是自然界最基本的物理定律之一，它断言在一个孤立系统中，能量既不会凭空产生，也不会凭空消失，只能从一种形式转化为另一种形式，或者从一个物体转移到另一个物体。这一看似简单的原理，却是理解地球气候系统运作机制的核心钥匙。气候系统的所有动态过程——从太阳辐射的入射与反射，到大气的运动与海洋的环流，从云的形成与消散，到冰川的进退与海平面的升降——本质上都是能量在不同形式之间转换与在不同空间之间传输的过程。因此，要深入理解气候变化的物理机制，首先必须建立起对气候系统能量平衡的清晰认识。

地球气候系统是一个由大气圈、水圈、冰冻圈、岩石圈和生物圈五个子系统组成的复杂耦合系统。这个系统接收来自太阳的短波辐射能量，同时向太空发射长波辐射能量。在地质时间尺度上，这种入射辐射与出射辐射之间的平衡关系决定了地球的长期平均温度；在季节乃至年际时间尺度上，这种平衡的扰动则驱动着天气的变化与气候的波动。当这种能量平衡被外力打破时——无论是由于太阳活动的变化、地球轨道参数的改变，还是大气成分的改变——气候系统就会做出响应，这种响应可能是温度的升降、降水模式的改变，或者更为剧烈的气候突变。

本章将从能量守恒定律出发，系统阐述地球气候系统的能量收支原理。首先介绍太阳常数与入射辐射的空间分布特征，分析地球轨道参数对天文辐射的影响。其次详细讲解地气系统的辐射平衡方程，包括短波辐射的入射、反射与吸收，以及长波辐射的发射与散逸。通过兰利图解法说明辐射平衡温度的计算过程。进一步探讨非辐射能量传输过程，特别是水汽相变潜热传输和感热通量输送。最后讨论能量平衡方程的不稳定性与气候敏感性的概念引入，为后续章节奠定坚实的物理基础。这些内容将帮助读者建立理解气候变化所必需的能量学框架，为深入探讨后续章节中更为复杂的动力学过程做好铺垫。

太阳常数与入射辐射的空间分布

太阳常数的基本概念

太阳常数是指在日地平均距离处，垂直于太阳辐射方向单位面积上单位时间内接收到的太阳辐射总能量，通常记作$S_0$。根据美国国家航空航天局（NASA）和世界气象组织（WMO）推荐的数值，太阳常数约为$1361 \ \mathrm{W/m^2}$。这一数值是通过数十年的卫星观测和地面校准不断精进的，目前其测量精度已达到约$\pm 0.5 \ \mathrm{W/m^2}$的量级。值得注意的是，太阳常数并非真正恒定不变，而是存在约$0.1\%$的周期性变化，这种变化与太阳黑子活动周期（约11年）密切相关。

理解太阳常数的物理含义需要考虑太阳的辐射特性。太阳是一个接近理想的黑体辐射源，其有效辐射温度约为5778 K。根据斯特藩-玻尔兹曼定律，黑体的辐射通量密度与其温度的四次方成正比：

$$E = \sigma T^4$$

其中$\sigma$为斯特藩-玻尔兹曼常数，约为$5.67 \times 10^{-8} \ \mathrm{W/(m^2 \cdot K^4)}$。将太阳有效温度代入，可得太阳表面的辐射通量密度约为$6.3 \times 10^7 \ \mathrm{W/m^2}$。太阳的总辐射功率约为$3.85 \times 10^{26} \ \mathrm{W}$，这些能量以球面波的形式向太空辐射。在日地平均距离（$1.496 \times 10^{11} \ \mathrm{m}$）处，由于辐射能量分散到更大的球面上，辐射通量密度降低到太阳常数的水平。

太阳常数$S_0$乘以地球的截获面积$\pi R_E^2$（其中$R_E$为地球平均半径，约6371 km）得到地球接收太阳辐射的总功率，约为$1.74 \times 10^{17} \ \mathrm{W}$。这一数值看起来巨大，但将其与地球表面积$4\pi R_E^2$相除，可得地球平均接受的太阳辐射通量密度约为$\frac{S_0}{4} \approx 340 \ \mathrm{W/m^2}$。这个因子$1/4$的出现是因为地球是一个球体，太阳光线只照射到朝向太阳的半个地球表面，而地球各点接收的辐射能量与其相对于太阳的方向有关。

入射辐射的纬度分布特征

地球球形几何的直接影响是：不同纬度接收的太阳辐射能量存在显著差异。在春分和秋分时节，太阳直射点在赤道，日照时数在各纬度相等，入射辐射的纬度分布由太阳高度角决定。对于北纬$\phi$度的某一点，春分日正午的太阳高度角为$90° - \phi$，到达单位水平面的太阳辐射通量密度为：

$$Q = S_0 \cdot \sin(90° - \phi) = S_0 \cdot \cos\phi$$

这一简单公式揭示了入射辐射随纬度增加而递减的规律：赤道处（$\phi = 0°$）接收最大辐射，而极地（$\phi = 90°$）则为零。当然，实际观测中需要考虑大气层的削弱作用（光学质量）、日地距离的季节变化以及日照时数的影响，但纬度分布的基本格局不变。

将上述分析推广到全球，可以得到全球平均的入射短波辐射分布。对于任意纬度$\phi$，假设日地距离变化忽略不计，春分日（或秋分日）单位纬向带接收的太阳辐射功率为：

$$dP = S_0 \cos\phi \cdot (R_E \cos\phi \cdot d\lambda) \cdot (R_E d\phi) = S_0 R_E^2 \cos^2\phi \cdot d\phi \cdot d\lambda$$

其中$d\lambda$为经向宽度。将此表达式对经度积分（从0到$2\pi$），得到单位纬向带接收的总功率：

$$P(\phi) = 2\pi S_0 R_E^2 \cos^2\phi$$

进一步对纬度从$-90°$到$90°$积分，可验证总功率为$\pi R_E^2 S_0$，即地球截获的总太阳辐射功率。单位纬向带的辐射功率随纬度的分布可用下式描述：

$$\bar{Q}(\phi) = \frac{S_0}{4}(1 + 0.477 \cos 2\phi)$$

这一公式表明，赤道地区接收的辐射约为全球平均的141%，而极地则不足全球平均的一半。这种辐射能量的纬度不均匀分布，是驱动大气环流和海洋环流的根本能量来源，也是形成气候带分异的基本原因。

日地距离的变化与天文辐射的季节变异

地球围绕太阳公转的轨道是一个椭圆，太阳位于椭圆的一个焦点上。根据开普勒第二定律，地球在近日点附近运动速度较快，在远日点附近运动速度较慢。地球轨道的偏心率约为0.0167，这意味着日地距离在一年中会有约3.3%的变化。目前地球的近日点（约$1.471 \times 10^{11} \ \mathrm{m}$）出现在每年1月初，此时北半球正值冬季；远日点（约$1.521 \times 10^{11} \ \mathrm{m}$）出现在7月初，此时北半球正值夏季。这种配置的净效应是：北半球的冬季略暖（南半球的夏季略冷），夏季略凉（南半球的冬季略暖），在一定程度上缓和了北半球季节变化的幅度。

日地距离的变化直接影响到达地球的太阳辐射总量。太阳辐射通量密度与日地距离的平方成反比：

$$S = \frac{S_0}{(1-e)^2}$$

其中$e$为相对日地距离偏差。当地球处于近日点时，$S$约为$1.034 S_0$；处于远日点时，$S$约为$0.966 S_0$。这种约7%的差异虽然不大，但对于年际和更长尺度的气候波动仍有一定影响。

然而，日地距离的季节变化对气候的影响远不及另一个因素——地轴倾斜角（黄赤交角）显著。地轴相对于轨道平面的倾斜角约为23.44°，且倾斜方向（指向）在一年中基本保持不变（忽略岁差运动的缓慢变化）。这意味着，在夏至日，北极倾向太阳，北半球获得较长的日照时间和较高的太阳高度角；在冬至日，北极背离太阳，北半球获得较短的日照和较低的太阳高度角。正是这种倾斜导致了季节的更替和天文辐射的显著年变化。

天文辐射（不考虑大气影响的太阳辐射）可以用日地距离、太阳赤纬（日光直射点的纬度）和日照时数来精确计算。对于纬度$\phi$的某一点，天文辐射日总量（日照时间内接受的辐射能量）可表示为：

$$Q_{daily} = \frac{S}{\pi} \left(\frac{\pi}{180}\right) d \left[h_0 \sin\phi \sin\delta + \cos\phi \cos\delta \sin h_0\right]$$

其中$S$为考虑日地距离的太阳常数修正值，$d$为日地距离（天文单位），$\delta$为太阳赤纬，$h_0$为半日弧长（日出到正午的时角），且$\cos h_0 = -\tan\phi \tan\delta$。

这一公式清楚地表明，天文辐射随纬度、季节和日地距离而变化。对于北纬45°的地点，夏至日的天文辐射日总量约为冬至日的3.5倍；而对于北极点，夏至日24小时均有日照（极昼），天文辐射日总量远大于赤道。

地球轨道参数对天文辐射的影响

轨道参数的基本定义

地球围绕太阳的公转轨道并非完美椭圆，而是受到其他行星引力扰动的影响，轨道参数在长时间尺度上会发生周期性变化。气候学中最为关注的轨道参数有三个：偏心率、倾角（地轴倾斜角）和岁差。

偏心率$e$描述椭圆轨道偏离圆形的程度，其定义为：

$$e = \frac{a-b}{a}$$

其中$a$为轨道半长轴，$b$为轨道半短轴。偏心率的变化范围约为0.005至0.058，当前值约为0.0167。当偏心率增大时，地球轨道的椭圆程度加剧，近日点与远日点的日地距离差异增大，这会影响不同季节接受的太阳辐射强度。

倾角（或称黄赤交角）$\varepsilon$是地轴相对于轨道平面的倾斜角度，当前值约为23.44°。倾角的变化范围约为22.1°至24.5°，周期约为41000年。倾角决定了太阳直射点移动的范围，也决定了极地地区极昼极夜的范围。当倾角增大时，高纬度地区的季节变化加剧：夏季更热、冬季更冷，年均辐射总量变化不大，但年较差增大。

岁差$\psi$描述地轴指向在空间中的缓慢旋转运动，周期约为25700年（完整的岁差周期）。岁差运动意味着，北半球夏季处于近日点还是远日点的配置会逐渐改变。当前，北半球夏季接近远日点，南半球夏季接近近日点，因此北半球季节变化相对缓和，南半球季节变化相对剧烈。在约11000年后，这种配置将完全颠倒，北半球将经历更为剧烈的季节变化。

米兰科维奇理论与冰期旋回

20世纪初期，塞尔维亚地球物理学家米卢廷·米兰科维奇（Milutin Milanković）提出了一个革命性的假说：地球轨道参数的周期性变化通过改变北半球高纬度地区夏季的天文辐射分布，引发了第四纪冰期旋回。这一理论后来被称为米兰科维奇理论，虽然在提出之初遭遇质疑，但在20世纪70年代深海沉积物和冰芯记录得到验证后，已成为理解第四纪气候变化的基石。

米兰科维奇理论的核心是冰期发生的机制。在北半球夏季，如果高纬度地区接受的太阳辐射足够低，使得夏季不足以完全融化上一冬天的积雪，积雪就会积累形成冰盖。白色冰盖的高反照率会反射更多太阳辐射，进一步降低地表吸收的辐射能量，形成正反馈机制。这种雪-冰-反照率反馈使得小的天文辐射变化被放大，最终导致冰盖的扩张或收缩。

定量研究表明，夏至日北纬65°的天文辐射变化与第四纪冰期旋回高度吻合。该纬度夏季的天文辐射变化主要受两个因素控制：岁差运动改变近日点相对于夏季的位置，倾角变化改变夏季的太阳高度角。当岁差使得北半球夏季处于远日点且倾角较小时，北纬65°的夏季天文辐射显著降低，有利于冰盖保存；当配置相反时，冰盖趋于消融。

地质记录显示，第四纪期间存在约10万年的冰期-间冰期旋回，这与偏心率的周期吻合。偏心率虽然本身对天文辐射的影响很小（其变化引起的辐射变化仅约0.1%），但它调制了岁差效应的影响强度：当偏心率较大时，近日点与远日点的辐射差异增大，岁差效应被放大；当偏心率接近零（轨道接近圆形）时，岁差效应被压制。因此，虽然倾角（约41000年周期）和岁差（约26000年周期）在冰芯记录中也有体现，但偏心率的10万年周期主导了冰期旋回的总体节律。

轨道参数变化对当代气候的潜在影响

米兰科维奇理论主要解释的是地质时间尺度（万年至十万年）的冰期旋回，但其思想对于理解当代和未来气候变化也具有启示意义。虽然当代轨道参数变化的速度极其缓慢（每世纪倾角变化约0.0002°），其对气候的直接影响远不及人类活动排放的温室气体，但轨道参数变化在理解自然气候变率中仍占重要地位。

过去100万年的气候记录表明，在约40万年前开始的最后一个冰期旋回中，间冰期（暖期）持续时间较长且温度较高，这与轨道参数的特殊配置有关。此外，在理解全新世（过去11700年）气候变化时，岁差的影响也不可忽视。研究表明，全新世早期北半球夏季太阳辐射较高，这与早期人类文明的农业发展可能存在关联。

在现代气候模拟中，地球轨道参数的变化通常作为边界条件输入气候模式。虽然其直接影响很小，但通过海洋热量传输、冰盖-反照率反馈和温室气体浓度的间接联系，轨道参数变化可能对百年尺度气候变率有微弱但可检测的影响。对于长期气候预测（千年至万年级别），轨道参数变化是必须考虑的因素，这也是理解未来冰期何时会再次降临的关键。

地气系统的辐射平衡方程

短波辐射的入射、反射与吸收

到达大气层顶的太阳辐射是气候系统的外部能量来源。如前所述，在日地平均距离处，大气层顶垂直于太阳方向的太阳辐射通量密度为太阳常数$S_0 \approx 1361 \ \mathrm{W/m^2}$。然而，并非所有这些辐射都能到达地球表面。当太阳辐射穿过大气层时，会发生散射、反射和吸收等过程，一部分能量被直接反射回太空，一部分被大气吸收，余下的才能到达地表。

描述大气对太阳短波辐射影响的常用参数是行星反照率$\alpha$，定义为被地球-大气系统反射回太空的短波辐射与入射短波辐射之比。当前全球平均行星反照率约为0.30，意味着约30%的入射太阳辐射被反射回去，余下约70%被地球-大气系统吸收。全球平均的短波辐射吸收为：

$$Q_{abs,SW} = \frac{S_0}{4}(1-\alpha) \approx \frac{1361}{4} \times 0.70 \approx 238 \ \mathrm{W/m^2}$$

这一数值约等于地球-大气系统向太空发射的长波辐射通量密度，表明在气候系统长期平均状态下，入射短波辐射与出射长波辐射基本平衡，全球能量收支接近于零。

行星反照率的空间分布高度不均匀。陆地的平均反照率约为0.15-0.20，海洋的反照率更低（约0.05-0.15），但冰雪覆盖地区的反照率可高达0.60-0.90。云的平均反照率约为0.50，是影响行星反照率最重要的因素之一。大气中的气溶胶（如硫酸盐、火山灰、灰尘等）也会通过散射和吸收影响大气层的反照率。典型的晴空大气反照率约为0.15，云天可达0.40-0.50，全球平均约0.23。

到达地面的短波辐射（也称总辐射或入射辐射）由两部分组成：直接辐射和散射辐射。前者是太阳光束直接照射到地面的部分，后者是经过大气散射后从各方向到达地面的部分。总辐射与大气层顶入射辐射之比取决于太阳高度角、大气透明度和云量等因素，可用比尔-朗伯定律描述：

$$Q = Q_0 \cdot T \cdot \sin h_\odot$$

其中$Q_0$为大气层顶的入射辐射，$h_\odot$为太阳高度角，$T$为大气透射率（考虑云和气溶胶影响后的大气透明度）。晴空无云时，$T$可高达0.75-0.85；厚云层覆盖时，$T$可低至0.1以下。

长波辐射的发射与散逸

地球表面和大气层向太空发射的能量主要以热红外长波辐射的形式进行。根据斯特藩-玻尔兹曼定律，任何温度高于绝对零度的物体都会发射热辐射，其辐射通量密度与温度的四次方成正比。地球表面温度约为288 K（15°C），据此计算的地表黑体辐射通量密度约为：

$$E_{surf} = \sigma T_{surf}^4 \approx 5.67 \times 10^{-8} \times 288^4 \approx 390 \ \mathrm{W/m^2}$$

然而，实际地表发射的长波辐射并不能全部到达大气层顶，因为大气中的水汽、二氧化碳、甲烷等温室气体会对长波辐射产生吸收和散射。大气层吸收的长波辐射会被再次发射，其中向太空方向发射的部分会散逸到宇宙空间，向地面方向发射的部分则形成大气逆辐射（大气向地面的长波辐射），这是温室效应的物理基础。

根据基尔霍夫辐射定律，在热平衡条件下，物体对某波长辐射的吸收率等于其发射率。因此，大气中温室气体的吸收率越高，其发射率也越高。当温室气体浓度增加时，大气层的吸收和发射能力增强，这会改变大气层向太空和向地面的长波辐射分配，导致温室效应的增强。

描述长波辐射传输的物理过程比短波辐射更为复杂，因为大气中存在多种温室气体和水汽，它们在不同波长范围内具有不同的吸收特性。长波辐射在大气中的传输过程可用辐射传输方程来描述：

$$\frac{dI_\lambda}{ds} = -\kappa_\lambda \rho I_\lambda + \kappa_\lambda \rho I_{\lambda,B}(T)$$

其中$I_\lambda$为波长$\lambda$的辐射强度，$\kappa_\lambda$为质量吸收系数，$\rho$为大气密度，$I_{\lambda,B}(T)$为温度$T$下的黑体辐射强度（普朗克函数）。方程右边第一项描述辐射被吸收导致的衰减，第二项描述热发射导致的增强。

在实际气候模式中，通常不逐波长求解辐射传输方程，而是将长波辐射谱分成若干波段，使用宽带近似方法（如k分布法）来提高计算效率。大气科学和气候学中常用的简化方法是分别考虑晴空大气和有云大气两种情况，计算大气层顶的出射长波辐射（Outgoing Longwave Radiation, OLR），这是卫星遥感监测地球辐射平衡的核心参数。

辐射平衡方程的完整形式

综合以上分析，地气系统的辐射平衡方程可表述为：大气层顶的短波入射辐射与长波出射辐射之差等于零（长期平均），或等于系统能量存储的变化率（短期变化）。用数学公式表示：

$$\frac{1}{4}S_0(1-\alpha) = OLR$$

其中$OLR$（Outgoing Longwave Radiation）为大气层顶出射长波辐射通量密度，约等于$\sigma T_e^4$，$T_e$为有效辐射温度。

这一方程描述的是全球平均值的情况。在空间分布上，不同纬度的辐射平衡存在显著差异。低纬度地区接收的太阳短波辐射大于全球平均值，而高纬度地区则小于平均值。根据前面的分析，北纬15°处接收的天文辐射约为全球平均的1.2倍，而北纬60°处仅约为0.7倍。然而，实际观测显示大气层顶的长波出射辐射随纬度的变化幅度要小得多，赤道与极地之间的OLR差异仅约20%。这种“热带过剩、高纬亏损”的短波辐射分布与相对均匀的长波辐射出射分布之间的矛盾，需要通过大气和海洋的热量输