《地球系统中溶解性有机质的分子特性：一项荟萃分析》研究报告
作者及发表信息
本研究由Han Lei、Hu Ang*、Mzuka Hellen Lucas、Chen Xingting、Shen Ji、Wang Jianjun共同完成。第一作者及主要研究人员来自中国科学院南京地理与湖泊研究所湖泊与环境国家重点实验室（南京 211135），合作作者来自南京大学地理与海洋科学学院（南京 210023）。通讯作者为胡昂（A. Hu）。该研究发表于学术期刊 Journal of Earth Science（《地球科学杂志》）2025年12月第36卷第6期（第2674–2688页）。
学术背景
本研究聚焦于环境科学和生物地球化学领域的一个核心主题：溶解性有机质（Dissolved Organic Matter， DOM）。DOM是全球最大的活性碳库，在水圈和陆圈的碳循环、生态系统功能以及众多物理、化学和生物过程中扮演着关键角色。尽管其重要性不言而喻，但科学界对于DOM分子层面的全球性理解仍存在不足。这些分子层面的属性，例如分子量、化学计量比（如H/C、O/C比）、不饱和度、芳香性以及碳的氧化状态，如何在地球不同系统（水体与陆地）中变化、彼此之间有何关联，以及受何种地理和环境因子驱动，仍是亟待解决的科学问题。
以往的研究多集中于单一或少数生态系统（如冰川、海洋、河流、湖泊或土壤）内的DOM特性，缺乏跨系统的综合比较和全球尺度的模式分析。近年来，基于傅里叶变换离子回旋共振质谱（Fourier Transform Ion Cyclotron Resonance Mass Spectrometry， FT-ICR MS） 的强大分析能力，科学家们发展了一系列DOM分子指数来评估其生物可利用性和化学特性，例如：质量数（Mass）、双键当量（Double Bond Equivalent， DBE）、修正芳香性指数（Modified Aromaticity Index， AImod）、碳的名义氧化态（Nominal Oxidation State of Carbon， NOSC）、氢碳比（H/C）和氧碳比（O/C）。这些指数分别反映了DOM分子的大小（质量）、生物难降解性（DBE、AImod、H/C）和氧化状态（NOSC、O/C）。然而，一个涵盖多种DOM分子属性、跨越全球各类生态系统的综合分析仍然缺失。
因此，本研究旨在通过一项大规模的荟萃分析（Meta-analysis） ，系统性地揭示DOM关键分子属性在全球水体与陆地系统中的分布规律、内在联系、纬度格局及其环境驱动因子。研究目标具体包括：（1）比较水体系统和陆地系统中DOM分子属性的差异；（2）探索不同分子属性之间的相关性；（3）识别主导全球DOM性质变异的关键属性；（4）描绘DOM属性随纬度的变化模式；（5）评估环境因子（如气候、pH）对DOM分子特性的影响。最终，此研究将为全面理解DOM在全球碳循环中的作用及其动态变化提供关键见解。
详细研究流程
本研究是一项基于文献数据整合的荟萃分析，其工作流程严谨而系统，主要包括以下几个步骤：
第一步：数据收集与筛选。 研究团队通过Google Scholar和Web of Science核心数据库，使用“有机质”、“FT-ICR MS”和“van Krevelen”等关键词，系统地检索了截至2022年6月发表的相关文献。数据筛选遵循严格标准： 1. 研究的必须是天然生态系统中的DOM。 2. 数据必须基于负离子模式电喷雾电离（negative-ion ESI）的FT-ICR MS分析获得，以确保数据来源的一致性，减少不同仪器带来的不确定性。 3. 能够从原始数据或文献中提取出加权平均的分子水平参数，包括：质量数（m/z）、DBE、AImod、NOSC、H/C和O/C。这些参数通过将单个分子的信息乘以其相对丰度后求和，再除以总丰度计算得出。 4. 原始FT-ICR MS数据需可获取或已公开报告，以便核对或计算。
经过筛选，研究最终纳入了204篇文献中来自2707个站点的数据。这些数据覆盖了海洋、河流、湖泊、水库、工程水体（如废水）、湿地、泥炭地和土壤等多种生境，形成了一个具有广泛地理和环境代表性的全球数据集。此外，研究还从WorldClim数据库为每个样本点提取了经纬度、高程以及19项生物气候变量（如年平均温度BIO1、年降水量BIO12等），并收集了pH等环境变量，用于后续的驱动因子分析。特别值得注意的是，本研究在先前仅关注H/C和O/C比的研究基础上，首次大规模整合并分析了质量、DBE、AImod和NOSC这四个关键分子指数。
第二步：统计分析方法。 为从海量数据中提炼出普适性规律，并控制不同研究和生境带来的随机效应，研究采用了多种先进的统计模型： 1. 差异性检验：使用Wilcoxon秩和检验来比较水体系统与陆地系统之间DOM属性的整体差异，以及不同生境之间的差异。 2. 相关性分析：采用Spearman秩相关分析来探究六个DOM分子指数之间的关联强度。 3. 主成分分析（Principal Component Analysis， PCA）：对六个DOM指数构成的矩阵进行PCA，以简化数据维度，识别主导DOM性质变异的主要方向（如生物可利用性维度、氧化状态维度），并可视化不同系统和生境的DOM属性分布。 4. 纬度模式分析：使用线性或二次模型拟合DOM属性与绝对纬度之间的关系，并通过Akaike信息准则（AIC） 选择最佳拟合模型。 5. 环境驱动因子分析： * 变量筛选：为避免多重共线性，首先使用方差膨胀因子（VIF） 将初始的36个潜在驱动因子缩减为15个（VIF），包括纬度、高程、pH及部分关键气候变量。 * 重要性排序：采用随机森林（Random Forest）模型（设置500棵树）评估各环境变量对DOM属性的相对重要性。 * 量化影响：运用线性混合效应模型（Linear Mixed-Effects Models），将“研究来源”作为随机因子，以控制样本间的非独立性，定量评估单个环境变量对DOM属性的影响大小和统计显著性。模型解释力通过条件R²来表示。
主要研究结果
1. DOM属性在系统和生境间的差异： 分析显示，水体系统中的DOM在质量数、DBE和AImod上均显著高于陆地系统，意味着水体DOM整体具有更大的分子量、更高的不饱和度和芳香性，可能包含更多难降解的分子组分。其中，河流沉积物中的这三个属性值最高（平均质量数447， DBE 12， AImod 0.46）。在氧化状态方面，废水（NOSC = 0.226 ± 0.06）和海水（NOSC = 0.133 ± 0.06） 的平均NOSC值最高，表明其DOM氧化程度更高。而在陆地系统中，泥炭地的DOM具有显著更低的质量数、AImod和NOSC，而农田土壤的DBE和AImod值相对最高，反映了不同土地利用方式对DOM分子特性的显著塑造作用。
2. DOM属性间的内在关联： 研究发现了DOM分子属性之间存在系统性的、可预测的相关关系，且这种关系因系统而异。 * 属性类别内：在生物可利用性属性内部，DBE和AImod与H/C比呈高度负相关（分子越不饱和、越芳香，氢饱和度越低）。在氧化状态属性内部，NOSC与O/C比在大多数情况下呈正相关（碳氧化态越高，氧含量通常也越高），但在海水中出现了微弱的负相关。 * 属性类别间：在陆地系统中，质量数与氧化状态（NOSC和O/C）呈强烈正相关，且质量数与生物可利用性属性（如与H/C负相关）也关联紧密。这表明在陆地环境中，更大分子的DOM往往氧化程度更高，氢饱和度更低。相比之下，在水体系统中，这种质量-氧化态的强关联较弱。此外，在陆地系统中，NOSC与生物可利用性属性（DBE、AImod、H/C）的关系更强。
3. 全球DOM性质变异的主导因素： PCA分析清晰地揭示了DOM属性变异的主要维度。第一主成分（PC1，解释约45%变异）的载荷主要体现在H/C（负）、DBE和AImod（正）上，可解释为“生物可利用性维度”。第二主成分（PC2，解释约31%变异）的载荷主要体现在NOSC和O/C（正）上，可解释为“氧化状态维度”。而质量数对两个主成分的贡献都较小。这一结果明确表明，在全球尺度上，H/C和AImod（即生物可利用性相关属性）是导致DOM性质变异的最主要因素。PCA图还显示，水体系统和陆地系统的DOM属性在PC1-PC2空间中有明显的分离，表明两者在分子组成上存在本质差异。
4. DOM属性的纬度格局： 研究首次在全球尺度上揭示了DOM分子属性具有显著的、可预测的纬度分布模式，尤其在水体系统中。 * DBE和AImod在水体中呈现显著的“驼峰型（hump-shaped）”纬度格局，其峰值出现在大约绝对纬度30°–50°的中纬度地区。这意味着中纬度水体的DOM具有最高的不饱和度和芳香性。 * NOSC在水体中则呈现单调递增的格局，即向高纬度地区，DOM的氧化程度逐渐升高。 * 这种纬度模式在海洋、河水和水库水等具体生境中同样显著。例如，在河水中，DBE、AImod和NOSC都表现出强烈的驼峰型模式。 * 在陆地系统中，仅AImod表现出随纬度升高而单调下降的趋势。
5. 环境驱动因子： 随机森林和线性混合模型分析共同证实，DOM分子属性的空间变异受到地理、气候和环境因子的显著驱动。 * 水体系统：年平均温度（BIO1） 是对DOM综合属性（PCOA1）影响最大的因子，其次是地理变量和pH。具体到单个属性，DBE、AImod和NOSC主要受地理变量（如纬度）控制，而质量数主要受年平均温度影响。此外，气候极端性（如等温性、最湿月降水量）对DOM的生物可利用性属性有强烈影响。 * 陆地系统：DOM属性主要受年平均温度、降水季节性、pH和温度年较差等因子驱动。其中，生物可利用性属性（DBE、AImod）主要受降水相关气候变量影响，而氧化状态（NOSC）主要受温度相关变量影响。
研究结论与意义
本研究通过一项全面的全球荟萃分析，系统地揭示了DOM关键分子属性在地球不同生态系统中的分布规律、内在联系及其环境驱动机制。
核心结论包括： 1. 系统差异：水体DOM比陆地DOM具有更高的分子量、不饱和度和芳香性，氧化态则在废水和海水中最高。 2. 内在关联：DOM属性间存在系统性关联，且在陆地系统中，分子大小与氧化态、氧化态与生物难降解性之间的联系更强。 3. 变异主导：H/C和AImod（生物可利用性相关属性）是全球DOM性质变异的最主要贡献者。 4. 纬度格局：在水体中，DOM的不饱和度和芳香性在中纬度（30°–50°）达到峰值，而氧化态则向高纬度单调递增。 5. 关键驱动：气候（尤其是年平均温度）和环境因子（如pH）是塑造DOM分子特性的关键驱动力，气候极端事件的影响也不容忽视。
科学价值： 1. 整合与升华：该研究首次将质量、DBE、AImod、NOSC等多个关键分子指数与经典的H/C、O/C比相结合，在全局尺度上提供了对DOM分子特性的全景式认识，填补了该领域的知识空白。 2. 揭示新规律：发现了中纬度水体DOM不饱和度和芳香性最高的“驼峰型”纬度格局，这是一个此前未被广泛报道的重要空间模式，为理解全球碳循环的区域差异提供了新视角。 3. 阐明驱动机制：定量解析了气候均值与极端事件、地理因素及pH等环境变量对DOM不同属性的差异化影响，增强了我们预测未来气候变化（如变暖、干旱）如何改变DOM性质及碳循环过程的能力。 4. 架起桥梁：通过关联分子属性（微观）与宏观环境梯度，研究深化了对DOM生物地球化学循环内在机制的理解，为建立更准确的地球系统模型提供了重要的参数化和验证依据。
应用价值：研究成果对于评估全球碳收支、预测气候变化下的生态系统反馈、以及管理水环境和土壤碳库具有指导意义。例如，理解不同纬度/生境DOM的分子顽固性，有助于改进全球碳循环模型中DOM降解速率的参数化。
研究亮点