《Systems》：Uncertainty-Aware Source Apportionment of Soil Heavy Metals in a Coal–Steel Industrial Zone Using Inter-Metal Graph-Informed Dirichlet Neural Posterior Estimation
编辑推荐：
工业化与采矿提升了土壤生态系统的重金属负荷，并对环境及人类健康造成风险。当相邻设施排放重叠的金属信号时，源解析变得困难，而传统受体模型仅返回贡献份额且未附不确定性陈述，进一步加剧了这一困难。研究人员针对土耳其宗古尔达克（Zonguldak）煤钢工业区，利用93
工业化与采矿提升了土壤生态系统的重金属负荷，并对环境及人类健康造成风险。当相邻设施排放重叠的金属信号时，源解析变得困难，而传统受体模型仅返回贡献份额且未附不确定性陈述，进一步加剧了这一困难。研究人员针对土耳其宗古尔达克（Zonguldak）煤钢工业区，利用93个地理参考表土站点与8种优先重金属，开发了一种概率化源解析框架。图形最小绝对收缩与选择算子（Graphical Lasso）估计以实测金属为节点的稀疏条件依赖图；在该图上的固定传播将每次观测扩展为24维关系表示；狄利克雷神经后验估计器将该表示映射至源贡献分布。狄利克雷参数化从构造上强制非负性与单位求和，并在每站提供95%可信区间；在留出模拟上的经验覆盖率为0.957（名义0.95）。全省平均后验均值贡献为：钢铁部门0.25、热电发电0.25、由成岩背景与未建模人为输入构成的复合残差类0.50。两个工业信号围绕Ere?li与?atala?z?集中，残差类在其他区域占主导。此处图编码金属间化学关系而非站点间地理关系；后验估计地图仅作事后可视化。该框架作为不确定性感知的方法学贡献提出，与既有受体模型的对比尚待开展。
研究背景方面，全球工业化、无节制采矿及集约农业使镉（Cd）、铅（Pb）、砷（As）、汞（Hg）等持久性有毒重金属进入土壤食物链，带来生物多样性和人体健康风险。准确溯源是风险管理前提。传统受体模型如正矩阵分解（PMF）、绝对主成分得分-多元线性回归（APCS-MLR）和化学质量平衡（CMB）存在三类短板：将样本视为独立而忽略空间依赖；基于线性假设难以刻画土壤异质非线性过程；仅给出点估计，贡献份额不强制非负与单位求和，不确定性只能靠事后重采样获得。标准机器学习把每站当作独立单元，抑制了金属间协方差所承载的多源叠加信息，且输出无内置可靠性度量，小样本下泛化性差。已有图神经网络与贝叶斯受体模型分别向预测和不确定量化推进，但未见将“金属间稀疏条件依赖”与“构造上约束于概率单纯形（probability simplex）的后验”结合的工作。宗古尔达克盆地同时布局Ere?li钢铁厂与?atala?z?燃煤电厂，二者经共享燃煤路径在Hg、As、Cu上指纹高度重叠，正是常规模型最易失效的场景，因此亟需一种附带校准不确定性陈述的源解析方法学。
研究人员开展的研究为：以土耳其宗古尔达克省93个地理参考表土站点（两轮采样，首轮36站已发表、次轮57站新增）的8种优先重金属（As、Cd、Cr、Cu、Hg、Ni、Pb、Zn，ICP-MS定量）为对象，构建“Graphical Lasso条件依赖图 + 固定图传播特征富集 + 模拟推理（SBI）下狄利克雷神经后验估计”的端到端概率管道。结论上，该框架在留出模拟取得0.957经验覆盖率（名义0.95），全省后验均值贡献为钢铁0.25、热电0.25、残差类0.50；工业信号聚于Ere?li与?atala?z?，残差类在余区占优；移除图或换用边际相关图仅使贡献偏移0.005–0.008，表明结论立于概率架构本身。论文发表于《Systems》，其系统科学贡献在于把8金属当作交互系统，将数据稀缺野外问题转为可摊销推理任务，并将显式不确定性嵌入环境受体建模。
关键技术方法概括：研究人员采用三部件耦合架构。其一，对log1p变换后观测用Graphical Lasso估计8金属稀疏精度矩阵，取偏相关为加权邻接，以条件独立剔除共享混杂带来的伪相关；该图固定不再反向传播。其二，依简单图卷积（SGC）原理做两次确定性邻域聚合，将8维原始签名扩为原始+一阶+二阶共24维图富集特征。其三，在物理混合模型下以对称狄利克雷先验生成36万对模拟观测-权重，噪声协方差取现场log1p值的Ledoit-Wolf收缩估计，训练两层128单元全连接网络（ReLU、Dropout 0.3、Softplus输出）近似后验为狄利克雷分布，用负对数似然优化，Adam学习率1e-3，批次2048，早停。样本队列即前述93站宗古尔达克表土。
研究结果部分，各小节结论如下。定量源解析与贡献份额：93站平均后验均值贡献钢铁0.25、热电0.25、残差0.50；残差类在0.6–0.9区间高频出现印证成岩背景主导，工业贡献呈左偏低值点源特征，与同区既往PMF识别的煤基与钢铁源身份一致（非数值对标）。地理参考分布与污染热点：后验均值制图显示热电影响?atala?z?沿海带及宗古尔达克中心城区，钢铁在Ere?li突出但范围较窄，残差在Zonguldak、Ere?li、?aycuma人口密集区偏高；图仅为事后可视化，模型不含站点坐标与距离，不学习空间自相关，但高贡献站伴随特征金属（Pb、Cd、Cu、Hg）高浓度，具地球化学一致性。不确定性分析与模型可靠性：工业核心39站95%可信区间均宽0.369，较背景54站0.312宽18.5%；背景站工业区间窄而残差区间最宽（0.381），工业核心站两工业源区间宽（0.384、0.375）源于燃煤指纹重叠致分割不定；模拟校准覆盖钢铁0.9678、热电0.9551、残差0.9470，均值0.9567，区间未膨胀即锐利且校准。参考剖面不确定性传播：自举重采样参考站使现场区间均宽升13.0%（钢铁扩18.8%），覆盖仍0.9550；留一参考站实验中位绝对偏移0.00637，主导源一致率98.5%，说明指纹由整体参考集支撑而非单站。消融与敏感：传播阶数0–3、邻接换Pearson/Spearman/无图、网络深度1–3、模拟量4.5万–36万、狄利克雷浓度0.6–2.0共18配置覆盖均0.954–0.960；图定义无关紧要（偏移0.005–0.008），仅先验浓度实质改变区间锐度与贡献（平均偏移0.0518），证明归因不依赖任意图选择。范围与局限：图连金属不连站点故未处理站间空间依赖；未与PMF/APCS-MLR/MCMC数值比对（现场无权重真值）；单时相表层土，不含pH/有机质，剖面垂向异质性未覆盖；参考剖面具盆地专属性，架构可迁移参数不可直迁。
讨论总结指出，该框架以架构贡献弥补传统受体“无内置不确定性”缺口：单纯形约束由狄利克雷构造保证而非事后强加，后验单次前向即得；工业重叠区不确定性反最高，点估计会掩盖此监管关键信息；图虽对归因数值无关，但提供偏相关层级的金属间物理化学解释（如Cr-Ni强正偏相关0.698示超基性母岩，Hg条件孤立示气相脱耦）。结论部分译文：本研究提出耦合稀疏图建模与贝叶斯推理的方法学框架，将土壤重金属源贡献恢复为概率分布。宗古尔达克盆地浓度异质性强，框架将其分解为附区间的源贡献；留出模拟经验覆盖0.957，区间未膨胀即校准且锐利；自举参考剖面使区间均宽增13.0%而覆盖保0.955，证明训练前剖面不确定性被传入后验而非被吸收；留一实验中位偏移0.00637。消融表明归因不依赖于传播阶数与图定义（换边际相关偏移0.005–0.008），仅狄利克雷先验浓度有实质影响，符合贝叶斯估计器应有表现。工业核心站区间比背景宽18.5%，因燃煤指纹重叠使两工业源分割不定，而背景站残差组成最不受约束；归因最不确定处恰为监管最需区间。框架只建金属间依赖不建站间空间依赖，未与既有受体模型数值对标，参考剖面具盆地专属性——这三点界定了后续引入站点邻接、同位素独立验证与外流域测试的明确方向。
