基于共同维度融合伽马与X射线光谱的土壤侵蚀表征

《Soil Science Society of America Journal》：Data fusion of gamma and x-ray spectra via Common Dimensions for soil erosion characterization

　　

本研究探讨了能量色散X射线荧光（EDXRF）和伽马射线光谱法（GRS）原始光谱在表征土壤再分布方面的协同应用。研究人员采用了一种名为共同维度（ComDim）的多块无监督机器学习技术，对从农业大田采集的31个土壤样本（0-3厘米深度）的核数据块和原子数据块进行共享潜模式识别。ComDim分析揭示了三个共同维度（CDs），解释了68.34%的累积方差。CD1主要由EDXRF数据块驱动（显著度≈90%），有效地区分了采样位置不同的样本；而CD3则主要由GRS数据块主导（显著度≈95%），并根据采集深度区分了样本。使用斯皮尔曼等级相关系数（ρ）进行的定量验证表明，CD得分与分析变量之间存在显著关联，尤其是CD3与铍-7（⁷Be）活度之间，这为土壤再分布评估提供了统计基础。这些发现表明，与单个数据块分析或传统的数据拼接方法相比，ComDim多块框架能更深入地揭示土壤动态信息，同时最大限度地减少预处理假设并保留原始光谱信息。

**基于共同维度融合伽马与X射线光谱的土壤侵蚀表征研究解读**

**一、 研究背景与意义**

土壤侵蚀是导致农业土壤退化的主要因素，直接影响表层土壤有机碳和必需养分的空间分布，从根本上改变其物理和化学性质，威胁全球超过15亿人口的生计。准确表征土壤组成和推断土壤再分布过程对于可持续土地管理至关重要。光谱学和核技术是评估土壤动态与健康的两种互补方法。能量色散X射线荧光（EDXRF）技术可用于测定土壤元素组成，而伽马射线光谱法（GRS）则用于量化土壤中的放射性核素浓度。其中，铍-7（⁷Be）作为一种天然存在的宇生放射性核素，半衰期短（53.3天），通过湿沉降沉积在土壤表层并迅速被土壤颗粒吸附，其垂直分布通常局限于土壤剖面上部30毫米，因此非常适合于评估短期土壤再分布。

尽管这两种技术能力强大，但各自存在局限。基于放射性核素的方法依赖于初始沉积空间均匀性、代表性未受扰动的参考库存等假设，且受探测器性能和核素半衰期限制。而EDXRF等地化学测量虽能提供详细的元素组成信息，却不能直接定量估算土壤再分布的幅度或时间。以往研究通常将这两种技术平行使用但独立解释，这可能限制了识别与土壤组成和再分布过程相关的联合模式。数据融合技术，特别是中水平融合（特征级融合），通过从不同数据源提取相关特征并将其拼接成单一矩阵，有望提供比单一技术更准确的知识和更好的推断。然而，将EDXRF与GRS在多块分析框架内进行整合的研究相对有限。

共同维度（ComDim）分析是一种基于中水平融合方法的多块无监督机器学习工具，它能够同时表示多个数据集，识别跨放射性核素和元素组成的共享模式。因此，本研究旨在利用ComDim融合从农业大田采集的土壤样本的EDXRF和GRS数据，并将此方法的性能与更传统的单块探索性分析进行比较。这项研究发表于《Soil Science Society of America Journal》，其重要意义在于探索了一种新的数据整合途径，以增强对土壤侵蚀和沉积物结合养分迁移的监测能力，为土地管理者识别侵蚀风险区域和指导可持续实践提供更有效的工具。

**二、 关键技术方法**

为开展研究，研究人员在巴西巴拉那州北部一个农业大田的缓坡上，于九个采样点采集了0-1厘米、1-2厘米和2-3厘米三个深度的土壤样本，并在附近平坦非侵蚀区选取了一个参考点，额外采集了5厘米深度的样本（用于获取低于检测限的⁷Be浓度），共计31个样本。土壤类型主要为“典型富铁红壤”（Rhodic Ferralsol）和“红壤性富铁红黏土”（Rhodic Nitisol），黏土含量高。

关键技术方法包括：1. **GRS测量**：使用铊激活碘化钠（NaI[Tl]）闪烁探测器测量样本的自然辐射，活时间为86，400秒，以获得统计上稳健的结果，光谱范围切割为238至2800 keV。2. **EDXRF测量**：使用手持式仪器在15 kV条件下进行测量，以利于较轻元素的荧光激发，光谱能量范围为0.98–15.00 keV。3. **数据融合与ComDim分析**：将两种技术切割并归一化后的原始光谱作为输入数据块，运用ComDim算法提取共享的共同维度。分析使用Julia编程语言和Jchemo.jl软件包进行，通过保留原始光谱信息并最小化预处理假设，旨在探索性地识别与土壤再分布过程相关的联合光谱模式。

**三、 研究结果**

**4.1 XRF和GRS光谱及浓度分布结果**
平均XRF光谱显示存在Al、Si、P、S、K、Ca、Ti、Mn、Fe、Cu和Zn等元素，其中Fe和Ti的峰强度较高，这与研究区域红壤和黏土中铁铝氧化物（如赤铁矿、针铁矿、三水铝石）为主的特征一致。平均GRS光谱识别出了²³⁸U衰变链（²¹⁴Pb和²¹⁴Bi）、²³²Th衰变链（²⁰⁸Tl、²¹²Pb和²²⁸Ac）、⁴⁰K以及⁷Be的特征发射。垂直分布剖面显示，在0-3厘米深度范围内，元素浓度相对稳定，而⁷Be活度则随深度增加呈指数下降，5厘米参考样本的⁷Be活度低于检测限，符合其作为短期土壤再分布示踪剂的预期行为。

**4.2 ComDim分析结果**
ComDim分析计算了前10个共同维度，其中前三个CDs解释了约68%的总方差，因此被保留用于进一步解释。显著度分析表明，CD1主要由XRF数据块驱动（λ≈90%），而CD3主要由GRS数据块主导（λ≈95%）。斯皮尔曼相关性分析证实，CD1得分与XRF获得的元素变量（特别是Fe、Ti、K、Ca和P）呈强显著正相关（ρ > 0.65），表明该维度主要与土壤地球化学变异相关。CD3得分则与⁷Be活度呈强显著正相关（ρ = 0.79），表明该维度捕获了与放射性核素衍生的土壤再分布过程相关的变异。

在CD1×CD3得分空间中，样本分布显示CD1能根据采样位置区分样本，这与XRF数据块在该维度占主导地位一致，反映了元素组成的空间变异。CD3则能根据采样深度区分样本，较浅深度的样本通常具有较高的得分，这与GRS数据块在CD3占主导地位及其与⁷Be的强相关性一致，表明CD3反映了与近期表层过程（影响放射性核素再分布）相关的变异。

**四、 讨论与结论**

**讨论部分总结：**
ComDim揭示的结构表明，其捕获的变异主要分配于元素组成（CD1）和放射性核素活度（CD3）之间，反映了XRF和GRS数据集的互补性。CD1的高显著度表明土壤矿物基质对共享结构有强烈影响，在该研究区域的富铁红壤和黏土中，这反映了氧化物空间分布的差异。CD3与⁷Be活度的关系为了解表层过程提供了额外见解，CD3得分可解释为原始土壤表面保存情况的代理指标，表明多块框架对作为放射性核素和地球化学标记物（如Fe和Ti）载体的细土壤颗粒的再分布很敏感。

样本在ComDim空间中的分布也反映了这种分离。CD1的变异主要与土壤组成差异有关，而沿CD3的分离则捕获了与放射性核素再分布相关的深度依赖性变化。同一地点不同深度采集的样本在CD1上离散有限，表明整个采样剖面元素组成相对稳定。然而，CD3上的变异遵循⁷Be预期的深度依赖性行为，强化了其对近地表过程的敏感性。需要指出的是，ComDim得分空间中样本的接近程度反映了其地球化学和放射性特征的相似性，而非直接指示土壤稳定性。相对于参考条件（代表未受扰动基线）的解释更为重要，远离参考点的样本可能反映了侵蚀、搬运或沉积的影响。

尽管空间模式复杂，但多元结构与独立推导的再分布速率之间的关系表明，即使在低扰动条件下，ComDim对细微的地球化学变化也很敏感。本研究样本量较小（n=31）且聚焦于特定土壤类型，未来研究需要在更大、更多样化的数据集上验证该方法的可转移性及其对土壤水分等环境因素的敏感性。尽管存在这些限制，目前的结果证明了多块数据融合方法在通过研究核与原子光谱数据的协同作用来增强土壤再分布和沉积物结合养分迁移监测方面的潜力。

**研究结论翻译：**
ComDim被证明是一种强大的多块算法，可有效评估来自互补光谱技术的变量的潜在影响。EDXRF和GRS原始光谱的融合能够识别出区分采样位置和深度的联合模式。具体而言，CD1（与EDXRF强相关）通过地球化学特征区分样本，而CD3（与GRS强相关）则成功地根据采样深度描述了样本。

通过斯皮尔曼相关性进行的统计验证证实，共同维度得分可作为⁷Be活度和元素浓度的可靠代理指标，从而在多块输出与已建立的土壤再分布示踪剂之间建立了定量联系。然而，当前研究受限于小样本量（n=31）和对特定土壤类型（富铁红壤和黏土）的关注。需要利用更大、更多样化的数据集进行进一步研究，以评估该方法的可转移性及其对环境因素（如土壤水分）的敏感性。

虽然这些结果表明原始光谱特征能有效捕获与再分布相关的变异，但利用这些特征将区域分类为受扰动或未受扰动仍不清楚。这种整合方法通过研究核与原子光谱数据之间的协同作用，增强了土壤再分布和沉积物结合养分迁移的监测。