一、研究背景与争议
亚洲季风系统深刻影响着亚洲大陆的降水格局。自晚中新世（约11 Ma前）以来，亚洲季风降水的变化机制一直是古气候学研究的热点，但学界存在争议。争议的核心在于如何解释C4植物（主要为热带、亚热带草本植物，采用C4光合作用途径）在过去地质时期的扩张。传统观点之一认为，C4植物的扩张指示了暖季降水(WSP)的增加，这可能与青藏高原的阶段性隆升有关，后者通过改变大气环流增强了亚洲季风。例如，巴基斯坦北部土壤碳酸盐碳同位素(δ13C)的显著转变被解读为C4植物在晚中新世的扩张，可能反映了季节性气候的增强和暖季降水(WSP)的增加。然而，也有反对意见，例如阿拉伯海西部、中国黄土高原(CLP)的多项研究，基于不同的代用指标（如叶蜡δD、土壤碳酸盐δ13C等），对降水变化的解读存在差异，甚至得出相反结论。
核心争议点：基于稳定碳同位素的代用指标进行季风降水重建时，C4植物的生态环境意义仍不明确，其扩张不一定单纯对应降水增加，还可能与气温、CO2浓度、季节性干旱增强、火频率增加等多重因素有关。
二、研究目标与方法
为解决上述争议，本研究摒弃了争议较大的稳定碳同位素指标，转而采用一种更为直接的植被记录方法——植硅体分析。
1. 科学目标：
- 重建东亚地区自晚中新世（约11 Ma前）以来的植被与季风气候演化历史，尤其关注上新世（约5.3至2.6 Ma前）这一关键时期。
- 量化分析C4草本植物的变化历史，并据此定量重建古降水量和古温度。
2. 研究材料与地点：
- 本研究选取位于中国关中地区的陕西蓝田作为关键研究点（北纬34.20°，东经109.24°）。
- 分析对象是从渭河盆地中采集的沉积物序列（覆盖灞河组、蓝田组、黄土-古土壤序列）中提取的植硅体化石。从133个样品中，筛选出38个保存完好（属于“尚好”到“很好”等级）的样品进行定量分析。
- 研究时间范围覆盖过去约11.0 Ma（自晚中新世开始）。
3. 研究方法论：
- 植硅体分析：植硅体是植物细胞中沉积的二氧化硅微体化石，其形态具有植物分类学特征。相比碳同位素指标，植硅体能直接区分主要的禾本科植物亚科（如羊茅亚科[Pooideae, 主要为C3植物]、黍亚科[Panicoideae, 主要为C4植物]、虎尾草亚科[Chloridoideae, 主要为C4植物]等），为古植被（C3/C4比例）和栖息地类型提供直接证据。
- 定量古气候重建：本研究结合了两套数据源。第一，建立了中国境内现代禾本科植物物种与气候条件（年均温MAT，年均降水MAP，最暖月降水WMMP，暖季[4-9月]降水WSP）的区域性物种-气候数据库。第二，对化石植硅体组合进行分析，判断优势的禾本科亚科。通过共存分析法和生态系统匹配法，将化石植硅体反映的植物群落构成与现代数据对比，进而定量重建过去的温度和降水。
关键数据处理与计算公式：文中通过对比森林指示植硅体（FI TOT）与禾草植硅体（GSSC）的比例估计了树木覆盖率：FI-t ratio = FI TOT/(FI TOT + GSSC)%。
文中使用了两种基于特定植硅体类型的比例来估算 C4植物生物量的方法：
- 最小值估算：基于黍亚科(Panicoideae)和虎尾草亚科(Chloridoideae)的植硅体比例，公式为 C4min (%) = (PAN + CHLOR)/(GSSC – OTHG) × (GSSC)/(FI TOT + GSSC)%。
- 最大值估算：基于所有PACMAD类群（包括黍亚科、虎尾草亚科、三芒草亚科、盔花亚科、小丽草亚科、芦竹亚科、扁芒草亚科）的植硅体比例，公式为 C4max (%) = (PACMAD TOT)/(GSSC – OTHG) × (GSSC)/(FI TOT + GSSC)%。
以上公式中，OTHG指其他未分类的禾草植硅体，被排除在计算之外以提高估算准确性。
三、主要研究结果
基于对蓝田剖面的详细植硅体记录和定量重建，研究得出以下主要发现：
1. 植被与草类群落的长期演化：
- 约11.0至4.2 Ma前（晚中新世至早上新世）：植被由C3和C4植物混合的草原主导，而非密闭森林。C4植物非常丰富（最小值估算平均~25%，最大值估算平均~75%），主要由适应温暖潮湿气候的PACMAD类（以黍亚科为主）占优势。
- 约4.2 Ma（约早上新世）的关键转折点：C4植物显著减少（最小值估算降至~7%，最大值降至~26%），而适应寒冷干燥气候的羊茅亚科(C3植物)逐渐在草类群落中占据主导地位（平均占比约67%）。同时，树木覆盖率（根据孢粉证据）也大幅下降。
- 约2.6 Ma至今（更新世及以后）：C3的羊茅亚科持续主导草类群落（平均占比高达93%），C4植物占比更低。这表明气候进一步变冷变干。
2. 古气候定量重建结果：
- 温度：年均温(MAT)呈现出清晰的阶段性下降趋势。在11.0至4.2 Ma期间较为温暖（11-15.3°C）；到4.2至2.6 Ma的上新世，MAT降至约9.7-11°C；进入2.6 Ma以来的更新世，MAT进一步降至约3.3-11°C。
- 降水：最关键的发现是季风降水的显著减少。最暖月降水(WMMP，亚洲季风季节的降水指示)的变化尤为明显。
- 11.0至4.2 Ma (温暖湿润期): MAP为 800–1673 mm；WMMP为 76–199 mm；WSP为 427–1197 mm。
- 4.2至2.6 Ma (冷却干燥期): MAP下降至 812–900 mm；WMMP下降至 76–130 mm；WSP下降至 417–540 mm。
- 2.6 Ma至今 (更冷干期): MAP下降至 441–900 mm；WMMP下降至 43–130 mm；WSP下降至 263–540 mm。
- 与其他记录的对比：本研究重建的温度和降水序列（尤其是约4.2 Ma的转折）与多项独立记录高度一致，包括黄土高原、临夏盆地、朝那等地的孢粉记录、哺乳动物化石丰度及物种组成变化（如喜湿润的Browser减少，喜干燥的Grazer增多）以及土壤碳酸盐δ13C值的正向偏移等。
四、与C4扩张相关争议的解释
本研究首次为东亚地区提供了关于C4草地的最早化石证据（约11.0 Ma前）。研究者承认，其植硅体数据（尤其是最小值估算）可能低估了某些环境下C4植物的生物量，但认为这种数据仍能可靠反映C4植物的相对变化趋势。研究认为，针对关键地点的碳同位素记录（如土壤碳酸盐δ13C）与植硅体记录之所以存在差异（例如，对某些层位植被的解释不同），可能有两个主要原因：
- 过去研究中对碳同位素信号（约-8‰至-9‰）的解读存在偏差，该数值可能主要反映了高水分压力下的纯C3植被信号，而植硅体证据则明确支持当时是C3/C4混合草原。
- 土壤碳酸盐的δ13C值不仅受植被类型和大气CO2影响，还对土壤呼吸通量（与初级生产力和年降水量密切相关）高度敏感。因此，δ13C的正向偏移可能更多地反映了降水减少和生产力的下降（导致土壤CO2通量降低），而非C4植物向C3植物转变的信号。这一论点得到了孢粉记录中同期树木覆盖率同步下降的有力支持。
五、核心结论与驱动机制讨论
论文的核心结论是：东亚上新世（约4.2 Ma前后）出现了显著的降温和季风降水减少事件。C4植物的减少、C3草类（羊茅亚科）的扩张以及年降水量、最暖月降水的定量数据都指向这一点。
研究者通过深入对比分析，支持并强化了亚洲晚新生代气候变化受全球冷却驱动的假说。
- 强相关性：本研究重建的羊茅亚科百分比（作为季风气候的定性代理指标）与全球温度和/或冰量记录之间显示出高达0.68的强相关性（在0.01水平上显著）。这表明亚洲季风变化与全球气候趋势紧密耦合。
- 物理机制：全球冷却通过影响热带辐合带(ITCZ)的位置来调控季风降水。降温会使ITCZ向赤道方向收缩，从而导致东亚季风降水带的位移和降水量减少。这一机制能够很好地解释季风强度与全球温度的正相关关系。
- 大气CO2水平的下降虽然是北半球冰期开启的诱因之一，但在上新世之前已低于触发冰期所需的阈值，因此CO2变化本身可能并非上新世中期气候变化事件的首要和直接驱动力，而冰盖反馈机制在放大全球冷却效应方面扮演了重要角色。
- 与全球海表温度不一致的时期：文中也指出，存在一段时期（晚中新世至早上新世，约11.0-4.2 Ma），东亚地区呈现温暖湿润的相对稳定状态，但这与全球海表温度记录的长期变冷趋势不完全一致。这表明，可能存在区域性的调节因素或其他驱动机制共同作用，但这些因素并未改变约4.2 Ma全球冷却对东亚季风的主导性影响。
最终观点：本文通过创新的植硅体记录，为研究上新世及更长时间尺度上东亚气候演变提供了独立的、高分辨率的直接证据。它有力地支持了在全球冰室-温室气候转型背景下，全球温度变化是东亚季风降水演变的关键强迫因子这一观点，深化了对区域气候对全球变化响应的理解。
参考文献：
Wang Hanlin, Lu Huayu*, Zhao Lin, Zhang Hongyan, Lei Fang & Wang Yichao, 2019. Asian monsoon rainfall variation during the Pliocene forced by global temperature change. Nature Communications, 10(1): 5272.
