📌 证据等级：A

地磁活跃期对呼吸道超额死亡的边际增量：气象条件配对年份设计

在气象条件几乎相同的两个年份里，地磁更活跃的那一年，呼吸道超额死亡是否更高？——这是把”天时”从”天气”里拆出来的关键一役。

正文

中医五运六气讲”天时”与”气候”共同决定疫病，但在现代统计里，”天时”（太阳活动与地磁扰动）和”气候”（温度、湿度）高度纠缠：地磁活跃的年份往往也是太阳活动的特定相位，而太阳活动又通过紫外线、气温间接影响流感流行。于是当研究者报告”地磁活跃期呼吸道死亡升高”时，批评者总可以反问一句：你拆得开吗？会不会只是气象条件的影子？

这个问题正是本系列此前一篇文章（#230，双通路拆分检验）留下的悬案：免疫通路（宿主昼夜节律、褪黑素）与气象通路（绝对湿度、温度）各有一批成熟证据，但从未在同一样本中把两条通路同时拆开。而”配对年份设计”（matched-year design）正是为这种质疑量身定做的减法实验。

配对年份设计的逻辑非常朴素：先不看地磁，只看气象。在二十年时间里，选出若干对”气象条件几乎相同”的年份——同样的冬季平均温度、同样的绝对湿度曲线、同样的流感流行强度。然后看每一对年份里，地磁更活跃的那一年，呼吸道超额死亡是否系统性更高。因为气象被”配对”掉了，剩下的差异只能归因于地磁。

2026年6月刚发表的一项全美研究，恰好把这个思路推进到了可执行的程度。Harfouch 等人用2000-2019年美国47个州的月度呼吸道感染（RTI）死亡数据，以地磁指数 Kp 的6个月移动平均作为暴露指标，用广义估计方程（GEE）泊松回归估计相对风险。结果发现：在1-12个月的所有暴露窗口中，6个月累积窗口与RTI死亡的关联最强；Kp 每升高一个四分位距，65-74岁、75-84岁、85岁以上三个年龄组的RTI死亡相对风险都显著升高（85岁以上女性 RR=1.16，95%CI 1.14-1.18；男性 RR=1.14，95%CI 1.12-1.16）。这是第一项在国家级尺度上把”地磁累积暴露”与”呼吸道感染死亡”直接挂钩的研究（PMID 42350784）。

为什么用6个月累积而不是当天Kp？这正是配对年份设计的灵魂所在。地磁对宿主的影响是慢性的、累积的——如果机制是昼夜节律紊乱和免疫调制，那么一次磁暴的当天效应很小，真正起作用的是连续几个月的”低强度持续扰动”。6个月窗口的峰值关联，与五运六气”客气加临、一气主六十日”的月尺度累积观，在时间结构上惊人地呼应：都是”一段时间的天时扰动”而非”一天的突发事件”。

但必须承认，PMID 42350784 是生态学研究，用的是全国月度汇总数据，温度、湿度只作为协变量进入模型，并没有做严格的”配对年份”。它的贡献是证明了”地磁→RTI死亡”关联在6个月窗口下稳健存在；而”边际增量”的直接检验——即气象完全匹配后地磁还能否贡献额外死亡——仍未有人做过。

这正是 #231 要立项的设计：气象条件配对年份 + 地磁活跃期事件研究法 + 6个月累积Kp暴露窗口 + 超额死亡估计。具体做法分四步：第一步，用绝对湿度（AH）和温度序列把2000-2025年划分成”气象相似年份组”（例如按冬季AH曲线的动态时间规整距离聚类）；第二步，在每一组内比较地磁活跃年与平静年的RTI超额死亡（超额死亡用 Serfling 类周期回归基线估算，PMID 19453440、9209852 给出了四种成熟方法）；第三步，用事件研究法（event study）检查地磁活跃期前后死亡率的时滞结构——如果免疫机制成立，峰值应出现在累积暴露之后数周到数月，而非当天；第四步，把6个月Kp累积窗口与其他窗口（1个月、12个月）做敏感性对比，确认窗口选择不是数据挖掘的产物。

这套设计的统计威力在于”配对”天然消除了气象混杂，而”超额死亡”天然吸收了年龄结构、慢性病负担等长期趋势。剩下的问题只有一个：地磁是否真的贡献了边际增量？

值得注意的是，太阳-流感关联本身是有争议的。Yeung 曾用1700-2000年的数据提出太阳黑子数可预测流感大流行（PMID 16806734），但 Towers 用更严格的统计方法重分析后认为，此前的阳性结论大多源于数据错误和任意分组，修正后没有显著关联（PMID 28847318）。这个对照非常重要：它说明”太阳-疾病”类关联极易被方法学瑕疵制造出来，配对年份设计恰恰是抵御这类虚假关联的防火墙——气象被固定后，地磁效应如果仍然出现，才配得上”边际增量”四个字。

证据与局限

已验证: 全美47州、2000-2019年数据显示，Kp指数6个月移动平均与RTI死亡显著相关，85岁以上人群RR最高（PMID 42350784）；韩国2001-2019年数据显示日Kp与总死亡、心血管死亡相关，且调整PM10后仍稳健（PMID 40768856）；263个美国城市2019年研究确认GMD与总死亡、心血管死亡相关（PMID 31511079）。

部分验证: 超额死亡估算方法成熟（Serfling周期回归、泊松回归、率差法四方法对比，PMID 19453440；121城市及时监测方法，PMID 9209852），但尚未有人将其与地磁暴露结合做”配对年份”检验；太阳-流感大流行关联存在正反两方证据（PMID 16806734 vs 28847318），结论依赖统计方法选择。

待验证: ①气象条件完全匹配后，地磁活跃期对RTI超额死亡的边际增量是否仍显著；②6个月累积Kp窗口是否优于其他窗口的机制解释（免疫累积 vs 气象代理）；③配对年份设计能否在北美以外（如东亚、欧洲）复制，检验纬度与地磁强度梯度。

Tag: 呼吸道感染