0.1摄氏度——这是美国玉米带过去50年夏季升温的速度，在NOAA所有气候分区里几乎垫底。但正是这个数字，被气候学家Roy Spencer选中，用来论证气候模型"严重高估"变暖趋势。问题是：当你把视野从12个州扩大到全美9大区域，把夏季数据换成全年均值，把单条模型曲线换成144组模拟的分布——这幅图景还会一样吗？
这事最近有了新进展。几周前，法庭强制公开了一批邮件，其中Steve Koonin和Roy Spencer的对话耐人寻味。Koonin直截了当地质疑：玉米带的对比是不是"挑樱桃"（cherry-picking）？Spencer答应查证，但后续是否回复、结论如何，邮件里没再出现。公开渠道至今沉默。不过对话泄露了数据来源和方法，这给了我复现分析的机会——不是照着做，而是看看如果做得更规范，会得出什么。
数据从哪来：NOAA的现成货架
Spencer用的玉米带数据，来自NOAA气候分区系统（ClimDiv）。这是一个层级聚合产品：底层是单个气象站，往上汇总到"气候分区"、州、区域。官网提供9大区域、主要流域、农业区等多种切片。玉米带覆盖11个州——Spencer的图标注12州，可能是自制索引，这点差异对结论影响不大，但增加了复现难度。
我下载了代码261的区域加权观测数据，计算1973-2022年夏季（JJA）趋势：每十年0.1摄氏度，95%置信区间[-0.07, 0.27]。确实很低。但注意这个区间宽度：统计上几乎和零没区别。更关键的是，这是ClimDiv所有区域里几乎最低的夏季趋势。换成全年数据，区域间差异就小得多——噪声低了，但"对比效果"也弱了。
模型对比：被丢弃的信息
模型端，Spencer从ClimateExplorer提取CMIP6数据，范围约100-81°W、40-46°N。他选择了"每个模型的集合平均"——或者仅单次模拟。无论哪种，都丢掉了关键信息：同一模型的多次模拟（ensemble runs）能反映内部变率范围。
我提取了144组独立模拟（历史情景+SSP245），计算相同区域的JJA趋势。然后做两件事Spencer没做：第一，给观测值画上趋势的不确定性区间；第二，展示模型集合的散布，而非单点。
三幅被忽略的图景
第一，观测本身就不确定。 玉米带0.1°C/十年的趋势，置信区间横跨正负。简单说：这地方的夏季变暖信号，统计上还没从噪声里完全浮出来。选它来"代表"观测，就像在误差棒里精准挑了个左端点。
第二，模型不是单点，是分布。 144组模拟的趋势呈明显分布：部分高于观测，部分低于，大量落在观测置信区间内。只画集合平均，等于把一团云压成一根针，再指责针没戳中目标。
第三，敏感度分层后，"热模型"叙事变复杂。 按瞬态气候响应（TCR）给模型上色——IPCC AR6评估其可能范围为1.4-2.2°C，很可能范围1.2-2.4°C——会发现高敏感度模型确实倾向更高趋势，但低敏感度模型也并非"更准确"。观测落在分布的偏左尾，但远谈不上"所有模型都失败"。
方法选择的连锁反应
Spencer的操作序列，每一步都导向同一方向：选夏季而非全年（噪声大、区域差异大），选玉米带而非全美（趋势几乎最低），选集合平均而非全分布（隐藏模型内部变率）。单独看每项选择都有理由，但组合起来，就构成经典的"挑樱桃"——不是伪造数据，而是让数据说你想听的话。
Koonin在邮件里点出的正是这个。他的身份值得注意：前奥巴马政府能源部次长，近年以"气候温和派"著称，绝非环保激进派。连他都质疑的选择性，说明问题不在意识形态，而在方法论。
复现教会我什么
学术争议里，原始数据和方法的透明度比结论更重要。Spencer的数据来源（NOAA ClimDiv）和模型来源（ClimateExplorer）都是公开资源，这给了外人复现的可能。但"可复现"不等于"已复现"——邮件泄露前，没人知道他具体选了哪个区域代码、哪个模型输出选项。
更深层的问题是：公众讨论需要什么样的"简化"？气候模型和观测的对比，本质是高维问题（空间、时间、变量、集合成员）。任何可视化都是压缩，但压缩到只剩一根针和一个小点，就丧失了科学对话的基础。好的科普应该向上抽象（找更稳健的模式），而非向下钻取（找最极端的切片）。
还没说完的事
邮件泄露后，Spencer的公开沉默构成另一种信息。学术规范里，被质疑后查证、回应、修正或坚持，都是正常路径。但"答应看看然后没下文"，在公共话语场中留下悬置。对关注此事的读者来说，这比任何结论都更值得记住：科学是过程，过程的可见性就是可信度的一部分。
玉米带的夏季还会变暖吗？ClimDiv数据更新到2022年，未来几十年的答案藏在排放情景和内部变率的交织里。但关于"如何讨论变暖"这件事，Roy Spencer的案例提供了一个具体教材：数据是公共的，选择是私人的，而检验选择的唯一方式，是让更多人能复现那条没被选中的路径。
热门跟贴