《Toxics》：Exceedance Probability and Toxicity-Value Sensitivity in Health-Risk Assessment of a Pb–Zn Slag Reprocessing Site Under the 2025 IRIS Arsenic Values
编辑推荐：
筛查级风险评估报告致癌风险（carcinogenic risk, CR）和危害指数（hazard index, HI）的单点估计，既不传达超标概率，也不传达其对修订毒性值的敏感性。研究人员在哈萨克斯坦Shymkent一处活动性铅–锌渣再处理场地，通过用蒙特卡洛
筛查级风险评估报告致癌风险（carcinogenic risk, CR）和危害指数（hazard index, HI）的单点估计，既不传达超标概率，也不传达其对修订毒性值的敏感性。研究人员在哈萨克斯坦Shymkent一处活动性铅–锌渣再处理场地，通过用蒙特卡洛模拟（Monte Carlo simulation）重新分析已发表数据集（24个土壤样品），针对成年户外工人，在源评估所用毒性值和2025年1月综合风险信息系统（IRIS）无机砷值下对两者进行了考察。在原始毒性值下，超过CR = 10?4的概率为70%（95% bootstrap区间59–80%），超过HI = 1的概率为82%。将2025年值应用于相同浓度，得到高21倍的平均风险以及99%和97%的超标概率。因此，筛查结论在两套毒性值下均成立，而定量估计则不然。砷浓度主导输出方差。由于摄入和皮肤接触承担工人超过99%的风险，风险降低需要阻断这些途径的措施：堆场周边工程覆盖、限制进入和职业卫生。植物稳定化（phytostabilization）对现场工人未模拟出获益；其价值在于限制场外粉尘，而后者未被建模。砷生物可及性（bioaccessibility）检测和更密集采样是优先事项。
研究背景与问题：历史Pb–Zn渣再处理在循环经济中处于矛盾位置。二次金属回收可减少永久处置和原生开采需求，但也会扰动已固结废体，在封闭、抑尘和径流控制不足时重新活化遗留污染物。哈萨克斯坦南部Shymkent前铅厂渣堆经约七十年火法冶金运行且当前仍在再处理，是区域严重污染场地之一。场地管理决策常以确定性评估为基础，报告致癌风险（carcinogenic risk, CR）和危害指数（hazard index, HI）单点估计并与阈值比较，但不能说明在合理暴露范围内超标概率，也不能说明各输入对输出方差的贡献。概率方法可提供超标概率与敏感性排序，但在后苏联中亚冶金遗留场地应用有限。另一个较少审查的不确定性是毒性值修订。2025年1月美国环境保护署（US EPA）综合风险信息系统（IRIS）修订无机砷评估，将口服斜率因子由1.5升至32 mg·kg
?1·d
?1，并降低口服参考剂量至6.0 × 10
?5 mg·kg
?1·d
?1，吸入单位风险未修订。修订前完成的砷主导场地评估未重算，其管理结论是否幸存尚不明确。因此，研究人员在Shymkent活动Pb–Zn渣再处理场地，利用已发表24个土壤样品数据，针对成年户外工人，在源评估毒性值和2025 IRIS砷值下开展概率重分析和毒性值敏感性分析，检验已发表分区与干预优先顺序是否成立。其意义在于提出可转移框架，用于判断管理结论能否在输入变异和毒性值修订后幸存，尤其适用于单一强偏态污染物主导风险谱的评估。
关键技术与方法：研究人员未开展新采样，而是重分析Pernebayev等已发表地球化学数据集。样本队列来源为该场地24条样带土壤样品和两个远端参考样品，于2026年4月24日采集，采样深度0–50 cm，沿四条样带距堆边0、1、5、10、50和100 m设置，参考样位于300 m和1.5 km；八种金属为Pb、Zn、Cd、Cu、As、Hg、Ni和Cr。方法包括最大似然拟合对数正态分布、Shapiro–Wilk检验、对不符合者采用经验重采样加对数正态核；用US EPA风险评价指南（RAGS）成人户外工人三途径暴露模型（摄入、吸入、皮肤接触）和蒙特卡洛模拟100,000次迭代，在源评估毒性值与2025 IRIS值两套参数下计算；用非参数bootstrap 1000次量化采样不确定性；用Spearman秩相关进行输入敏感性排序；并开展5、10、20年管理情景分析。
研究结果：
3.1 浓度分布：通过分布拟合与检验发现，八种金属中六种符合对数正态，Zn和Hg不符合并用经验重采样；As的几何标准差（geometric standard deviation, GSD）为3.94，As浓度范围为95.2至10,900 mg·kg
?1，跨两个数量级以上。Pb、Zn、Cd、Cu和Ni强相关，As与除Hg外金属弱相关，因此独立采样假设对结论影响小。替代分布模型使超标概率变化约4个百分点，但未改变结论。
3.2 CR和HI分布：通过蒙特卡洛输出分布分析发现，两套毒性值下CR和HI均强右偏。源评估值下CR四分位距为7.91 × 10
?5至6.38 × 10
?4，99百分位达8.20 × 10
?3，均值7.38 × 10
?4高于中位数2.26 × 10
?4；HI均值8.03、中位数3.12。2025 IRIS值下极上尾模拟值可超过1，最大7.34，已超出线性低剂量外推范围，应视为尾部权重指示。
3.3 超标概率：通过对模拟分布计算超标概率发现，源评估值下P(CR > 10
?4)为69.9%，P(CR > 10
?3)约六分之一；2025 IRIS值下分别为99.3%和84.0%。P(HI > 1)由82.1%升至97.3%，P(HI > 10)由19.1%升至57.4%。Bootstrap显示源评估值下P(CR > 10
?4)为59–80%，2025值下为98.4–99.9%；所有重复中均超过50%。相对口服生物有效性（relative oral bioavailability, RBA）降低对源评估值结果影响更大。
3.4 已发表点估计位置：通过将已发表点估计置于模拟分布中发现，源评估值下23.1%模拟暴露超过已发表站点均值CR，2.8%超过最坏情况；2025值下分别为89.0%和52.7%。2025值下平均CR为源评估值的21.3倍，平均HI为4.7倍。毒性值修订对估计的改变远大于确定性到概率方法的改变，后者约15–30%。
3.5 污染物贡献和敏感性：通过贡献分解和Spearman秩相关发现，As在两种参数下均占平均CR的99.8%，Cr占0.2%，Cd和Ni合计低于0.05%；As占源评估值下平均HI的91.7%，Pb占6.9%，Cd占1.0%。As浓度主导CR，ρ = +0.867，超过下一输入三倍；HI中As的ρ分别为+0.852和+0.904。暴露持续时间对HI无影响；吸入率影响可忽略，源评估值下吸入占砷CR的0.17%，2025值下为0.01%，摄入和皮肤接触分别占71.4%和28.4%。
3.6 区域风险分布：通过分区计算发现，Zone 1（0–10 m，n = 16）和Zone 2（10–100 m，n = 8）中，Zone 2平均CR为Zone 1的1.81倍，平均HI高1.72倍，但Zone 1超标概率更高（77.5%对60.3%）。Zone 1污染均匀，Zone 2空间斑块化且上尾更重，由西部样带50和100 m次高值驱动。留一法显示排除两个西部样带样本之一可反转排序；分层bootstrap的Zone 2:Zone 1平均CR比95%区间为0.16–9.9，区域比较未被数据解决。
3.7 情景预测：通过5、10、20年S0、S1、S2情景分析发现，S1即Zone 1工程围堵使平均CR降至无干预的约0.716、0.680、0.662（源评估值下），2025值下类似，约减少三分之一；S2即加Zone 2植物稳定化与S1差异小于0.3%。植物稳定化粉尘减少效率在0–100%变化对结论无实质影响。Zone 1围堵效果强烈依赖受体在两区间的时间分配，但筛查结论不依赖该分配。
讨论部分总结：讨论区分观测结果、模型推导结果和管理含义。4.1指出点估计不传达超标概率：确定性估计在本数据集近似分布均值，但概率分析说明69.9%的合理暴露超过10
?4，建议强右偏分布场地的筛查评估报告超标概率。4.2指出2025 IRIS修订使平均CR提高21.2倍、HI提高4.7倍，定性结论不变而定量估计不稳健；修订前后评估不可比，应报告毒性值年份并重算。4.3指出As主导精度，应优先加密As采样而非仅做职业暴露调查；区域上Shymkent全国随机评估未包括As，而本研究表明As而非Pb主导风险。4.4指出距源距离不是风险良好代理，已发表三层次分区未获风险数据支持；Zone 1均匀污染、Zone 2斑块严重，均值和超标概率应同时报告，并偏向主导风向加密采样。4.5指出受体决定干预有效性：Zone 1围堵对职业工人有效但取决于停留时间，植物稳定化对现场工人无额外建模获益，其价值在限制场外粉尘和未建模受体；仅职业受体评估会系统性偏好屏障措施，需扩展至居民。4.6总结限制：单次采样、总浓度、时间代表性、Hg挥发未考虑、分析质控缺失、采样深度、毒值路线外推、区域面积和时间分配、样本量仅24、RBA假设、颗粒物排放因子默认、情景参数来自文献以及单一受体。
研究结论翻译：对来自一个活动性Pb–Zn渣再处理场地的已发表数据集进行概率重分析表明，2025年1月IRIS对无机砷的修订对定量风险估计的改变远大于任何建模选择：对于相同实测浓度，估计的成年户外工人平均致癌风险高21倍，估计超过US EPA风险管理范围上限（10
?4）的概率从约70%（95% bootstrap区间59–80%）增加到超过98%。因此，源评估的定性筛查结论在两套毒性值下均成立，而修订前后所做的估计不可比；该结论对抽样不确定性和所测试的替代浓度模型具有稳健性，但仍以暴露模型的占用时间（occupancy）、生物有效性（bioavailability）和颗粒物排放（particulate-emission）假设为条件。相对于已发表的确定性评估，概率分析将该数据集的中心估计仅改变约15%，但增加了点估计无法传达的内容：超标概率以及已发表值在分布中的位置。主要不确定性来源是砷浓度分布，其主导输出方差并基于24个样品。这一小空间数据集使区域间风险排序未能解决，且仅限于单一职业受体意味着结果不扩展到居民或儿童。
