第一章 阿伦尼乌斯的阁楼:一次改变世界的计算 / Chapter 1 Arrhenius's Attic: A Calculation That Changed the World
第一章 阿伦尼乌斯的阁楼:一次改变世界的计算
Chapter 1 Arrhenius's Attic: A Calculation That Changed the World
1895 年深秋,斯德哥尔摩。一间被雪压得吱呀作响的阁楼里,34 岁的斯万特·阿伦尼乌斯(Svante Arrhenius)正伏在一张宽大的橡木桌上,面前摊着瑞典地质调查所寄来的岩石样本清单、几十张写满积分号的手稿、以及一支被他咬得几乎秃了的铅笔。他要做的事情,用今天的话说叫"数值模拟"——但 1895 年没有电子计算机,没有 FORTRAN,甚至没有像样的手摇计算器。他要用手,把大气中二氧化碳浓度从 300 ppm 翻到 600 ppm 之后地球会暖多少这件事,一步步算出来。
整整一年之后,他在《Philosophical Magazine》上发表了 41 页的论文《On the Influence of Carbonic Acid in the Air upon the Temperature of the Ground》(论大气中碳酸对地面温度的影响)。结论是一个数字:5–6 °C。
这个数字本身重要,但更重要的,是它第一次把"温室效应"从哲学猜想变成了可量化的科学问题。在阿伦尼乌斯之前,所有人都知道"大气像毯子";在阿伦尼乌斯之后,人们开始争论"这条毯子究竟有多厚"。本章要回答的核心问题是:温室效应这个概念,是如何从模糊猜想发展为可量化科学问题的? 我们将沿 1824—1859—1896 这条主线,把接力棒一棒一棒地捡起来。
---
1.1 傅里叶的玻璃罩:1824 年的一次"错误" / Fourier's Glass Dome: A "Mistake" in 1824
故事必须从一位数学家讲起。让·巴普蒂斯特·约瑟夫·傅里叶(Jean-Baptiste Joseph Fourier)一生最出名的成就是热传导方程,但他在地球科学上的一个"错误"反而改变了一切。
1824 年,傅里叶向法国科学院宣读了一篇短文《Remarques générales sur les températures du globe terrestre et des espaces planétaires》(关于地球及行星空间温度的总体评注)。他注意到一个奇怪的事实:地球表面平均温度大约是 15 °C(288 K),但如果把地球当成一个简单的黑体(blackbody),只考虑它从太阳那里接收的能量再以热辐射的形式发出去,温度应该只有 −18 °C(255 K)。两者之间有 33 K 的差值。傅里叶写道:
> *「The temperature of the ground can be augmented by the interposition of the atmosphere, which acts like the glass of a hothouse, because the heat easily penetrates the air in the form of light, but has more difficulty in escaping in the form of obscure heat.」* > > (地面温度可以被大气所"插入"的大气所抬高,大气的作用就像温室的玻璃,因为热量容易以光的形式穿透空气,却较难以"不可见热"的形式逃逸。)
注意,傅里叶在这里其实犯了一个物理上的小错误——玻璃温室保温的真正机制不是阻挡红外辐射,而是阻止对流(convection),这是 19 世纪末瑞利勋爵(Lord Rayleigh)才指出来的。但是,这个"错误"的比喻——大气像玻璃罩——太有生命力了,它变成了"温室效应"这个名词的来源。
傅里叶的工作为什么重要?因为他做了三件事: 1. 提出了问题:地球为什么比它"应该"更暖? 2. 给出了框架:大气对辐射有选择性透过率(selective transmissivity); 3. 留下了空白:他没有量化那 33 K 里有多少来自大气、多少来自海洋和陆地。
整整三十多年,没有人去补这个空白。
---
1.2 丁达尔的分光仪:把猜想钉死在实验台上 / Tyndall's Spectrophotometer: Nailing the Hypothesis to the Lab Bench
1859 年,爱尔兰物理学家约翰·丁达尔(John Tyndall)在一间伦敦皇家学院的地下室里,架起了一台他自己改进的分光仪(spectrophotometer)。他要做的事情非常专一:测量不同气体对红外辐射的吸收能力。
丁达尔的实验装置可以简化为四块:① 一个发热的铜立方体作为红外光源;② 一个充满待测气体的黄铜管;③ 一个灵敏的温差电堆(thermopile)作为红外探测器;④ 一个盖着冰水的遮光箱来隔绝背景辐射。他把管子里的空气依次替换为氧气、氮气、二氧化碳、甲烷、水蒸气,记录温差电堆的读数变化。
结果让所有人都意外了。氧气和氮气几乎不吸收红外辐射——它们是"透明"的。而二氧化碳和水蒸气(后来的甲烷、一氧化二氮)则强烈吸收特定波长的红外。丁达尔在 1861 年的论文中兴奋地写道:
> *「The atom links itself to the ether, and the two vibrate together; but the vibrations of the atom and of the ether must be isochronous before they can combine. It is, in fact, the sympathy of vibrations that determines the absorption.」* > > (原子与以太耦合,二者共同振动;但原子与以太的振动必须同步才能结合。事实上,正是振动的"同频共振"决定了吸收的发生。)
这"振动的同频共振"这句话,是 19 世纪实验物理学最美的一段表述。丁达尔虽然没有量子力学,但他凭直觉抓住了本质:分子吸收红外辐射,是因为入射光子的频率与分子的某些振动模式的本征频率(natural frequency)相匹配。
从今天的角度看,丁达尔的发现打开了三个关键认识: 1. 大气是"有色"的——它对可见光(来自太阳)几乎透明,但对红外(来自地表)有选择性地吸收。 2. 温室气体是少数派——CO₂、CH₄、N₂O、H₂O 才是"主角",N₂ 和 O₂ 反而是"背景"。 3. 物理量化的方向——下一步要算的不是"温度差多少度",而是"哪条吸收线被多强地打开"。
但是丁达尔也只走到了这里。他没有建立一个数学模型,把"分子吸收强度"和"全球温度"连接起来。
---
1.3 阿伦尼乌斯的近似:他能做的和不能做的 / Arrhenius's Approximation: What He Could and Couldn't Do
接下来是本书最关键的一次接力。1896 年,阿伦尼乌斯发表了他的长文。一个常常被忽略的事实是:他不是一个人在战斗。他借用了他的同事、瑞典地质学家阿维德·霍格姆(Arvid Högbom)的数据——霍格姆在研究火山与石灰岩循环时估算过地质时期大气 CO₂ 浓度的变化范围,从 1 倍到 3 倍、5 倍、10 倍不等。阿伦尼乌斯把这些地质情景代入他的辐射-对流平衡方程(radiative-convective equilibrium equation),计算每一个情景下的全球平均温度变化。
他做的核心数学,本质上是一个一维的辐射收支模型。让我们把它拆开来看,因为这个推导的思路,今天我们还在用。
1.3.1 入射能量
太阳常数(solar constant)在 1900 年前后的测量值大约是 $S_0 \approx 1360 \text{ W m}^{-2}$,这与今天的精确值 $1361 \text{ W m}^{-2}$ 几乎一致。地球单位面积接受的平均能量是:
$$F_{\text{in}} = \frac{S_0 (1 - \alpha)}{4}$$
其中 $\alpha \approx 0.30$ 是行星反照率(planetary albedo),即被反射回太空的太阳光比例。$4$ 这个分母来自几何投影——地球是个球,但只有朝向太阳的截面能截获阳光,所以全球平均要除以 4。代入数字:
$$F_{\text{in}} = \frac{1360 \times 0.7}{4} = 238 \text{ W m}^{-2}$$
1.3.2 出射能量(无大气情况)
如果没有大气,地表就像一颗 255 K 的黑体,按照斯特藩-玻尔兹曼定律(Stefan-Boltzmann law)向外辐射:
$$F_{\text{out}} = \sigma T^4$$
其中 $\sigma = 5.67 \times 10^{-8} \text{ W m}^{-2} \text{ K}^{-4}$。平衡时 $F_{\text{out}} = F_{\text{in}}$,解出:
$$T = \left(\frac{238}{5.67 \times 10^{-8}}\right)^{1/4} \approx 255 \text{ K} = -18\,°\text{C}$$
这就是傅里叶看到的那个"应有温度"。但地表实际是 288 K(15 °C),多了 33 K。
1.3.3 加入吸收层
阿伦尼乌斯在他的模型里,假设大气由若干个等温吸收层(isothermal absorbing layers)组成,每一层对红外辐射的吸收率(absorptivity)记为 $a_i$。他需要解的不是单一温度,而是一组辐射传递方程(radiative transfer equation):
$$F_{\text{top}} = \sigma T_s^4 (1 - a_1) + \sigma T_1^4 a_1 (1 - a_2) + \cdots$$
其中 $T_s$ 是地表温度,$T_i$ 是第 $i$ 层大气的温度。每多一层吸收,逃逸到太空的红外就少一份,于是地表必须升高温度才能重新建立平衡。
关键来了:吸收率 $a_i$ 本身依赖于气体的浓度和分子的吸收截面(absorption cross-section)。阿伦尼乌斯用了一个现在叫做比尔-朗伯定律(Beer-Lambert law)的近似:
$$a_i \approx 1 - e^{-\kappa_i N_i}$$
其中 $\kappa_i$ 是吸收系数(absorption coefficient),$N_i$ 是该层气体的柱密度(column density)。当浓度翻倍,$N_i$ 翻倍,$\kappa_i N_i$ 翻倍,吸收率非线性地增大。
1.3.4 阿伦尼乌斯的关键假设
在用今天的标准看,阿伦尼乌斯的模型充满了"可以原谅的粗糙":
- **忽略了水蒸气反馈**——实际上水蒸气是最强的温室气体,CO₂ 升高会通过升温间接增加大气水蒸气含量,从而**正反馈**放大升温。**真实的温升远大于纯 CO₂ 效应**。
- **把对流处理为固定通量**——他用了**平衡对流调整**(convective adjustment)的假设,但没有参数化对流强度对温度梯度的依赖。
- **没有云**——他假设天空万里无云,而实际上云既是温室(吸收红外)又是遮阳(反射太阳),这是 20 世纪气候模型最棘手的部分。
- **纬向与季节都做了平均**——他给的是全球年平均,但真实地球两极升温比赤道大得多(**极地放大效应**,polar amplification)。
正是这些"不能做的",让他的数字 5–6 °C 看起来比今天 IPCC 估计的 2.5–4 °C(CO₂ 翻倍、考虑所有反馈的平衡气候敏感度,equilibrium climate sensitivity,ECS)偏高了约 1.5–3 倍。但要在 1895 年把水蒸气反馈、云、极地放大都算清楚,那不是人类大脑能在一年里手工完成的。
---
1.4 数字的回声:5–6 °C 是怎么传到今天的 / Echoes of a Number: How 5–6 °C Travelled to Today
阿伦尼乌斯的论文发表后,几乎没有引起媒体关注。倒是他的同事、也是他的学术对手,克努特·奥格斯特罗姆(Knut Ångström)1901 年发表了一篇反驳文章,质疑他高估了 CO₂ 的吸收效应。奥格斯特罗姆做了一个聪明的实验:用一对发射红外线的铜板(一个密封有 CO₂,一个没有),发现两者辐射几乎相同。他据此宣称"CO₂ 浓度的变化对气候没有显著影响"。
这个反驳其实犯了一个错误:他用的红外源温度太高,辐射峰在短波(近红外),恰好落在 CO₂ 吸收的窗口外。今天我们用维恩位移定律(Wien's displacement law)一眼就能看出来:
$$\lambda_{\text{max}} = \frac{b}{T}, \quad b = 2.898 \times 10^{-3} \text{ m K}$$
地表温度 288 K 对应的辐射峰在 10 微米附近,正是 CO₂ 在 15 微米的吸收带旁边。所以地表向太空发出的红外,会强烈被 CO₂ 吸收;而奥格斯特罗姆那个高温铜板(~700 K)辐射峰在 4 微米,恰好"躲过"了 CO₂ 的吸收带。这是一个历史级别的"实验设计错误"——但它在公众心理上,让温室效应沉默了将近半个世纪。
直到 1938 年,英国工程师乔治·卡伦达(Guy Callendar)重新拾起阿伦尼乌斯的接力棒。他用 1900–1935 年全球气温记录和化石燃料排放记录,第一次给出了 CO₂ 浓度上升与全球升温相关的统计证据。他算出的"翻倍 CO₂ 升温"大约是 2 °C——比阿伦尼乌斯保守得多,却与现代中位估计非常接近。卡伦达之后是汉斯·苏斯(Hans Suess)和罗杰·雷维尔(Roger Revelle)的同位素测量,查尔斯·基林(Charles Keeling)的 Mauna Loa CO₂ 曲线(1958 年至今的"基林曲线",Keeling Curve),以及真锅淑郎(Syukuro Manabe)和理查德·韦瑟尔德(Richard Wetherald)1967 年用一维辐射-对流模型给出的、与阿伦尼乌斯原理上几乎相同但结果精细得多的估算。
把这条线索拉直,我们看到的是:
```mermaid flowchart TD F["Fourier 1824<br/>Atmospheric barrier hypothesis<br/>(qualitative)"] T["Tyndall 1859<br/>Infrared absorption measured<br/>(experimental)"] A["Arrhenius 1896<br/>First hand calculation<br/>(quantitative, 5-6 C)"] C["Callendar 1938<br/>Statistical evidence<br/>(observational)"] M["Manabe & Wetherald 1967<br/>1D radiative-convective model<br/>(computational)"] K["Keeling 1958-present<br/>CO2 time series<br/>(monitoring)"] I["IPCC AR6 2021<br/>ECS best estimate 3 C<br/>(consensus)"]
F --> T T --> A A --> C C --> M K --> M M --> I A -. gap of ~70 years .-> M ```
注意这里的"接力"特征:每一步都不是"取代"前一步,而是叠加新的工具。阿伦尼乌斯用傅里叶的概念、丁达尔的实验数据;真锅用阿伦尼乌斯的方程、阿洛伊斯(Arrhenius)的物理,但加上了计算机、水蒸气反馈和云的处理。
把这条历史脉络对照今天的 IPCC AR6 评估:平衡气候敏感度(ECS,CO₂ 翻倍时全球平均地表温度的最终变化)的最可能值是 3 °C,很可能范围(likely range)是 2.5–4 °C。阿伦尼乌斯的 5–6 °C 偏高了 1.5–3 倍,但数量级正确——在一个没有计算机、没有量子力学、没有卫星、没有风洞的时代。
---
1.5 微观与宏观的桥梁:分子振动为何能撬动行星 / The Micro-Macro Bridge: Why Molecular Vibrations Move Planets
到这里,读者会自然地追问一个非常物理的问题:CO₂ 在大气中只有 420 ppm(百万分之 420),按体积算只占 0.042%。这么稀薄的一点点气体,怎么能左右一颗行星 12,700 公里直径的能量平衡?
答案分三层。
第一层:吸收截面的不对称。 太阳表面温度约 5800 K,按维恩位移定律,它辐射的峰值波长在 0.5 微米(可见光)。地表 288 K 的辐射峰值在 10 微米(远红外)。N₂、O₂ 这些双原子分子没有电偶极矩(electric dipole moment),它们的振动模式不能与电磁波耦合,所以对可见光和红外都"沉默"。但 CO₂ 是一个线性三原子分子(O=C=O),它有四种振动模式(vibrational modes):对称伸缩、非对称伸缩、两个简并的弯曲。后三种具有电偶极矩的变化,因此能强烈吸收红外——特别是 4.3 微米(伸缩带)和 15 微米(弯曲带)这两个基带(fundamental bands)。
第二层:黑体辐射函数的覆盖。 真正决定温室效应强弱的不是单一吸收线,而是"整个红外光谱里有多少能量落进了吸收带"。普朗克黑体辐射函数告诉我们,288 K 黑体在 13–17 微米区间集中了地表向上辐射能量的相当大份额——而 CO₂ 的 15 微米吸收带恰好"罩"住了这个区间。少量气体 × 高浓度能量 × 强吸收截面 = 显著升温。
第三层:吸收的"饱和"与"翼区"。 这是更深的一层。CO₂ 浓度从 0 增加到几百 ppm,15