《Environmental Pollution》：Integration of geostatistical and hot spot analysis to identify drivers of radionuclide spatial heterogeneity in the high-mountain lake basin topsoil
编辑推荐：
本研究评估了亚美尼亚塞凡湖流域（高加索地区最大的高山淡水水体）表土中放射性核素空间异质性的驱动因素，采用集成统计与空间分析方法。共分析170个土壤样品的Ra-226、Th-232、K-40、Cs-137及总β活度，同时测定原位剂量。统计分析显示天然放射性核素之
本研究评估了亚美尼亚塞凡湖流域（高加索地区最大的高山淡水水体）表土中放射性核素空间异质性的驱动因素，采用集成统计与空间分析方法。共分析170个土壤样品的Ra-226、Th-232、K-40、Cs-137及总β活度，同时测定原位剂量。统计分析显示天然放射性核素之间存在强相互关联（采用Spearman检验），Th-232和总β活度呈双峰分布模式，表明存在不同的源区。土地利用分析确认其对K-40和Cs-137分布有显著影响，Cs-137在扰动低的土壤中较高，耕地中最低，而K-40在森林土壤中明显升高。采用地统计学方法（包括经验贝叶斯克里金(Empirical Bayesian Kriging, EBK)和Getis-Ord Gi? (Gi?)热点分析）模拟这些参数的空间异质性。分析识别出流域西南部火山岩主导区存在Ra-226、Th-232和镭等效活度(radium equivalent activity, RaEq)的统计显著热点，东部湖相衍生土壤中为冷点。K-40在西北部农作区有额外富集，提示长期农业实践的可能贡献。Cs-137在极小扰动土壤中表现出有限但明确的热点聚集，反映沉积后沉降过程。辐射风险评估显示Ra-226和Th-232是西南部定居区的主要风险贡献因子。研究表明，集成统计推断、空间建模和放射性核素地球化学可实现基于过程的土壤放射性解释，并为复杂山地环境的污染评估、放射生态靶向监测和管控优先级排序提供可迁移框架。
论文解读：《Environmental Pollution》发表的“Integration of geostatistical and hot spot analysis to identify drivers of radionuclide spatial heterogeneity in the high-mountain lake basin topsoil”
研究背景方面，土壤是天然和人为引入放射性核素的主要储库，其定量评估对确立辐射安全标准、保护环境及理解放射性核素行为的空间与环境过程具有重要意义。天然放射性核素（K-40、U-238衰变系包括Ra-226、Th-232衰变系）的空间变异性主要受母岩地质来源和矿物成分控制；人工核素Cs-137作为大气核试验与事故的遗留物，主要滞留于表层土壤，其再分布受土地利用、侵蚀、沉积、扰动等地表过程强烈影响；而土地利用与管理（耕作、混拌、施肥等）会重塑天然主控的空间格局。塞凡湖流域（亚美尼亚Gegharkunik州）作为高加索地区最大的高山淡水系统，生态敏感且社会经济关键，区内地质构造复杂（富火山岩建造，如玄武岩、凝灰岩(tuff)、安山岩等），并有集约化农业（马铃薯、畜牧等）和采矿活动（如Sotk金矿，属黄铁矿、毒砂等硫化物多金属矿床），存在天然放射性本底偏高与技术增强天然放射性物质(TENORM)的复合风险。然而目前该流域缺乏空间集成且经统计验证的土壤放射性研究，难以区分地质本底、地表过程与人为活动的各自贡献，也无法科学支撑监测分区与风险管控，因此研究人员开展了此项研究。
为开展研究，研究人员用到的主要关键技术方法包括：在塞凡湖所在的Gegharkunik州基于约5 km网格布设采样，采集170个表土（0–10 cm）样品并混合为复合样，样品经干燥、过筛（<2 mm）后，采用高纯锗(HPGe)γ谱仪测定Ra-226、Th-232、K-40、Cs-137活度（用LabSOCS无源效率校准，用IAEA-447、IAEA-465标准物质质控）；用iMatic?无气自动α/β计数系统按ISO 18589-6测定总β活度；用Automess 6150AD剂量率仪在距地1 m处测原位γ剂量率；基于活度计算镭等效活度(RaEq)、空气吸收剂量率(ODRA)、年有效剂量当量(AEDE)、超额终身癌症风险(ELCR)，并将活度转为等效元素浓度（eK、eRa、eTh）及比值（eTh/eRa、eTh/eK、eRa/eK）以支持地球化学解释；统计上用R语言进行Shapiro-Wilk正态性检验、Spearman秩相关与聚类热图、Kruskal-Wallis检验及Dunn事后检验（Bonferroni校正）以分析土地利用（耕地、森林、草甸、牧场）的影响；空间上在ArcMap 10.6中用经验贝叶斯克里金(EBK)插值生成连续分布面并输出预测标准误，采用固定距离带5000 m的Getis-Ord Gi?统计识别热点（高值聚集）与冷点（低值聚集）。
研究结果部分，3.1节统计特征与相关：研究人员通过描述统计、直方图、箱线图及Shapiro-Wilk检验发现所有变量原始数据均显著偏离正态分布，Ra-226极端正偏（偏度3.7，峰度19.5）且对数变换亦不正态，指示受离散天然热点控制；Th-232和总β呈双峰分布，暗示至少两个地球化学源区（火山岩 vs 湖相沉积）；Cs-137强非正态、均值高于中位数，右尾长，反映不均匀分布与少数高值点；K-40与原位剂量率相对对称。Spearman相关热图与层次聚类显示Ra-226与Th-232最强正相关（ρ=0.89），Ra-226与K-40次之（ρ=0.69），三者与原位剂量率聚为一支，代表统一天然辐射与剂量集群；Cs-137与天然核素关联较弱，主要受地表过程控制；海拔相关虽显著但弱于核素间相关，地形仅为次级修饰因子。
3.2节土地利用对核素浓度的影响：Kruskal-Wallis检验表明土地利用显著影响Cs-137（p<0.01）和K-40（p<0.01），对Ra-226、Th-232无显著影响。Dunn事后检验显示Cs-137在牧场显著高于耕地和草甸，因牧场扰动少利于表层滞留，耕地因长期耕作将污染层混入深层而稀释表面活度；K-40在森林土壤显著偏高，与森林多位于未扰动的火山质高地（天然富K矿物）及有机质细粒保K有关，但也叠加了土地利用信号；Ra-226、Th-232不受近期农业管理影响，说明它们主要由岩性/矿物成分主导，其偏高不代表人为污染而是天然本底富集。
3.3节EBK空间分布：EBK插值得Ra-226、Th-232、RaEq、总β的高值区一致分布于流域西南部/南部火山岩区，低值区在东部/东北湖相沉积、砂质土区，证实岩性为主要驱动，双峰统计对应“火山富集域”与“沉积低背景域”两类环境；K-40除对应天然背景外，在西北部农作区有独立高值区，提示地质背景叠加长期农业管理（如钾肥施用、马铃薯栽培）的混合信号；Cs-137呈碎片化斑块，高值多在高位、低扰动区（利于沉降滞留），低值在耕作低地（耕作混拌稀释），与土地利用结论吻合；Sotk金矿区未见天然核素升高，因赋存于蛇绿岩（通常贫 primordial 核素）且仅为采掘点，矿石外运至阿拉拉特金提取厂处理，无原地尾矿/加工残留广泛影响周边土壤；任何采石/加工活动的影响可能是局地的，需针对性调查。
3.4节热点分析：Getis-Ord Gi?统计显示Ra-226、Th-232、RaEq的热点（Gi_Bin≥2）集中在西南火山岩区（分别占样本7.6%、5.3%、5.3%），冷点（Gi_Bin≤?2）主要在湖以东沉积区（分别为5.9%、7.1%、7.1%），验证了EBK的梯度格局且具有统计显著性；K-40热点除了部分与火山区重叠外，在西北部农作区形成独立热点（占4.7%），冷点多在东部（7.6%），反映地质与农业复合控制；总β热点/冷点模式近似天然核素（西部热、东部冷），说明U-Th系主导了总β空间分异；Cs-137仅有2.4%热点（高位低扰区）、4.7%冷点（中部农耕平原），显著占比最低（共7.1%），分布受沉积后过程控制，未形成大规模强聚集。研究人员由此区分出流域内两大天然源域：西南火山岩富放射性核素域与东部湖相沉积低背景域；K-40额外叠加西北部农业信号；Cs-137为遗留沉降示踪。
3.5节源控制判别（地质 vs 人为）：综合统计、土地利用、EBK、Gi?及等效元素浓度表明，Ra-226、Th-232、RaEq热点与火山岩重合，强相关、eTh/eRa均值4.79（高于上陆壳UCC的3.89），反映富Th火山母岩加上风化中Th（较U系）更低迁移性导致表土相对Th富集，属岩性主控；约7%样eRa、12%样eTh超出UCC值，但属局部而非全流域污染；K-40为混合信号，西北部耕区eTh/eK、eRa/eK较低，指示相对K富集可能来自长期钾肥与耕作累积；Cs-137完全受沉降后地表过程（扰动、耕作、侵蚀）重塑，是历史人为污染的示踪剂；采石活动在当下尺度未系统对应热点，需靶向调查；Sotk矿区无天然核素异常，因为岩性贫 primordial 核素且加工在外地。集成EBK与Gi?能有效分辨区域梯度、局部异常与统计显著簇，从而区分地质驱动域与人为/地表过程影响。
3.6节剂量与辐射风险：原位剂量率、AEDE、ELCR空间格局与Ra-226、Th-232、RaEq一致，高值在南/西南部火山岩区；大部分区域AEDE低于全球均值0.07 mSv/y，但南部部分人口集中区AEDE>0.07 mSv/y，ELCR在两处南部超全球均值2.9×10-4，Ra-226与Th-232是主要贡献因子；K-40、Cs-137热点对年剂量与ELCR塑造有限。尽管若干采石场邻近南部聚居区，本研究尺度未发现系统叠加，但仍建议将这些区作为优先靶向监测区。研究说明天然地质放射性、集约土地利用、遗留沉降信号可共存于敏感高山湖流域，Ra-226/Th-232偏高在重叠居民点、采石、频繁扰动区才具环境意义；Cs-137主要指示历史污染留存与再分布。
讨论部分总结：研究人员在讨论中指出，Ra-226和Th-232强相关、共现热点且明确关联西南火山岩区，证实岩性是高天然放射性的主要控制因素；东部湖相衍生土壤一致为冷点，说明母质决定低天然本底；等效元素浓度与比值显示局部相对Th富集与火山源及风化中Th更低迁移性有关。K-40受混合控制，西北耕区独立富集与农业管理（历史钾肥、混拌、细粒保蓄）有关，更低eTh/eK、eRa/eK支持相对K富集。Cs-137碎片化分布、高位低扰区高值、耕地低值，证明现今格局受土地利用重塑，是历史沉降示踪剂。辐射风险绝大部分低于全球阈值，但南部部分聚居区AEDE与ELCR超标，需优先监测。研究贡献在于提供了一个可迁移框架，集成空间统计与放射性核素地球化学以识别土壤放射性驱动过程，改进污染评估、监测优先级与复杂山地放射生态管理。
结论部分翻译：本研究应用集成统计与地统计学框架以检验塞凡湖流域表土放射性核素异质性主要由地质背景控制还是被人为地表过程修改。经验贝叶斯克里金(EBK)与Getis-Ord Gi?互补使用实现了高分辨率空间插值与热点/冷点的显著分类，有助于区分地质与人为控制。Ra-226和Th-232显示强相关、共现热点、明确关联南部及西南部火山岩区，证实岩性为天然放射性主要控制因素；东部湖相与沉积衍生土壤一致形成冷点，表明母质关联的较低天然本底。等效元素浓度与比值显示与火山源有关的相对局部Th富集及选择性风化保留。K-40还显示混合控制：部分地质成因之外，西北部耕作区明确富集与土地利用显著影响，提示农业管理（可能为历史肥料投入与土壤混拌）的附加作用；耕区较低eTh/eK和eRa/eK进一步支持相对K富集。Cs-137呈碎片化空间分布，高位极小扰区较高、耕地因耕作混拌而偏低，表明现今格局受土地利用控制；因而Cs-137可作为山地土壤扰动历史与历史污染持续存在的环境示踪剂。辐射风险评估显示大部分区低于全球安全阈值，但局部南部和西部区超出AEDE与ELCR参考值，强调需在采石场与聚居区附近靶向监测。本研究的主要贡献是提供了一个可迁移框架，集成空间统计与放射性核素地球化学以识别土壤放射性控制过程，改进污染评估、监测优先级排序与复杂山地环境的放射生态管理。
