《Sustainability》：Heavy Metal Contamination in Arid Coastal Saudi Arabia: An Integrated Pollution-Index and Multivariate Approach Toward Environmental Sustainability
编辑推荐：
摘要译文：土壤中潜在有毒元素（Potentially Toxic Elements, PTEs）的污染是一个全球性的环境问题，可严重影响生态系统、农业发展以及人类健康。为评估潜在有毒金属的污染程度，本研究在沙特阿拉伯延布（Yanbu）和胡莱斯（Khulais）
摘要译文：土壤中潜在有毒元素（Potentially Toxic Elements, PTEs）的污染是一个全球性的环境问题，可严重影响生态系统、农业发展以及人类健康。为评估潜在有毒金属的污染程度，本研究在沙特阿拉伯延布（Yanbu）和胡莱斯（Khulais）地区的不同地点采集了32份土壤样品，分析了八种PTEs（Co、Cr、Cu、Fe、Ni、Pb、V和Zn）。所用污染指数包括地累积指数（Geo-accumulation Index, Igeo）、富集因子（Enrichment Factor, EF）、污染因子（Contamination Factor, CF）和污染负荷指数（Pollution Load Index, PLI），并采用主成分分析（Principal Component Analysis, PCA）和层次聚类分析（Hierarchical Cluster Analysis, HCA）进行源识别。由于从土壤中吸收的金属可沿土壤—植物路径进入农作物，表征农业土壤中潜在有毒金属富集的程度及可能来源，是后续任何土地管理或食品安全评估的必要第一步，尤其是在一个大部分粮食依赖进口的国家。PCA和HCA的结果表明，阿拉伯地盾（Arabian Shield）相关矿物的基岩是这些潜在有毒元素可能的主要来源。各PTE的平均EF值按降序排列为：Co (0.93) > Ni (0.92) > Cu (0.80) > Zn (0.74) > V (0.65) > Cr (0.59) > Pb (0.27)。这表明研究区土壤未富集或仅轻微富集PTEs。根据地累积指数值，几乎所有PTE均被归类为实际无污染（Igeo < 0），例外是Ni，在胡莱斯的四个点位达到了无污染至中度污染等级（Igeo ≥ 0）。CF和PLI值表明土壤平均处于无污染水平。总之，综合污染指数和多元统计框架为该干旱沿海地区提供了一种实用的地球化学筛选工具；以这种方式监测土壤质量，对水和土地资源受限环境中的可持续土地管理具有积极贡献。
论文解读：干旱沿海沙特阿拉伯土壤重金属污染的综合评价与可持续性启示
一、研究背景与目的
潜在有毒元素（Potentially Toxic Elements, PTEs）在土壤中的污染是全球关注的环境问题，可对生态系统、农业发展和人体健康造成严重危害。尽管这些元素天然存在于土壤中，但采矿、冶炼、化石燃料燃烧、交通、农药、垃圾填埋、农业化肥等多种人为活动会显著提高其浓度。PTEs可通过土壤—作物—食物链途径发生生物累积和生物放大，进而威胁食品安全与人体健康。此外，高浓度的PTEs还会降低生态系统的功能容量，减少可耕地面积，降低作物产量，并污染农产品。对于沙特阿拉伯而言，这一问题尤为突出，因为该国约80%的粮食依赖进口，国内农业仅能满足一小部分消费需求，若本土农田遭受重金属污染，将进一步压缩本已有限的粮食生产能力，直接威胁国家粮食安全与可持续发展目标。
尽管已有研究在沙特阿拉伯其他地区运用污染指数结合多元统计方法评价土壤质量，但针对红海沿岸重要农业区——延布（Yanbu）与胡莱斯（Khulais）的类似研究尚属空白。此前的研究虽提供了基线信息，但既未提供完全可重复的分析，也未明确对比两个不同的农业区域以查明区域地球化学变异性。因此，本研究的核心目标在于：应用综合污染指数（Igeo、EF、CF、PLI）与多元统计方法（PCA和HCA），系统评估延布和胡莱斯地区土壤中Co、Cr、Cu、Fe、Ni、Pb、V、Zn八种PTEs的污染程度，并识别其潜在来源。本研究的创新之处不在于方法本身，而在于首次在该区域综合应用上述四类互补污染指数并与PCA/HCA结合，通过独立指数间的相互验证来交叉确认污染状态，将指数结果与多元源识别相关联以区分地质成因与人为贡献，并识别出需要重点监测的具体点位。
二、研究方法概述
研究人员在延布（15个样品，编号Y1–Y15）和胡莱斯（17个样品，编号Kh1–Kh17）两地共采集32份表层土样（深度≤10 cm），每个采样点由三个子样混合而成。样品经风干、研磨、过筛后，采用硝酸-盐酸（王水）消解，利用电感耦合等离子体—原子发射光谱法（Inductively Coupled Plasma–Atomic Emission Spectrometry, ICP–AES）在沙特阿拉伯吉达ALS地球化学实验室测定八种金属浓度。质量控制采用标准参考物质（EMOG-17和OREAS 45f）、方法空白和重复样检验，回收率为72%–104%，重复样相对百分差异（Relative Percent Difference, RPD）除Ni和Pb外均低于10%。数据分析使用IBM SPSS 23和Python 3.10，计算描述性统计、Pearson相关、Spearman秩相关，并通过Mann–Whitney U检验比较两地差异。PCA前进行Kaiser–Meyer–Olkin（KMO）和Bartlett球形检验确认适用性。污染指数包括地累积指数（Igeo）、富集因子（EF，以Fe为归一化元素）、污染因子（CF）和污染负荷指数（PLI），背景值采用Turekian和Wedepohl的页岩平均值。
三、主要研究结果
3.1 研究区PTEs点位变化
研究人员测定了32份土壤样品中八种PTEs的浓度。平均浓度由高到低依次为Fe（43,518.75 mg kg
?1）、V（78.34 mg kg
?1）、Ni（62.31 mg kg
?1）、Zn（61.44 mg kg
?1）、Cr（48.84 mg kg
?1）、Cu（35.34 mg kg
?1）、Co（16.84 mg kg
?1）、Pb（4.38 mg kg
?1）。相对标准偏差（RSD%）从11%（Zn）到59%（Ni）不等，表明不同金属的点位间变异性存在显著差异。所有金属的最高浓度均低于农业土壤的植物毒性阈值和作物安全吸收范围，表明研究区内没有任何单个采样点达到可能危害作物健康或食品安全的水平。
3.2 区域及全球对比
研究人员将本研究测得的PTEs浓度与其他干旱地区研究进行了对比。平均而言，本研究测得的多数金属浓度高于沙特阿拉伯Al Qassim中部地区、Al-Muzahmiya、埃及Qena、沙特阿拉伯Jazan、阿拉伯湾以及Turekian和Wedepohl的背景值，但低于埃及尼罗河三角洲Kafr El-Sheikh省的值。Fe浓度高于表2中除背景值外的所有对照。Pb浓度则低于表中所有其他地点。Ni、V、Zn均高于Huraymla、利雅得西北部、Al-Muzahmiya和Al Qassim中部地区的报道值。多数金属浓度接近自然地壳背景值，而非接近工业和农业影响区的升高值，表明研究区土壤总体上反映地质成因背景。
3.3 正态性检验
统计检验显示，Pb（偏度0.79）和V（偏度0.25）呈正偏态，其他金属呈负偏态。峰度值介于?2.01到0.52之间。除Cr和Zn外，其余PTEs均不符合正态分布（p < 0.05）。
3.4 Pearson相关分析
所有金属对均呈现显著正相关（p < 0.05，多数p < 0.01），其中Co–Cr（r = 0.94）、Co–Fe（r = 0.87）、Co–Ni（r = 0.99）、Co–V（r = 0.95）和Co–Zn（r = 0.55）相关性最强，暗示共同来源。Pb与Co、Cr、Cu、Ni呈负相关。Spearman秩相关也显示Co、Cr、Cu、Fe、Ni、V之间存在强正相关（系数0.79–0.96），指示共同的地质成因。Pb与其他金属呈稳定负相关（r = ?0.64至?0.78），Zn与多数金属呈中等正相关（r = 0.49–0.76）。
3.5 多元分析源识别
PCA之前，KMO = 0.79，Bartlett球形检验显著（χ2 = 522.43，p < 0.001），适合进行因子提取。Kaiser准则保留了一个主成分，解释了81.6%的总方差。这一单一优势成分与阿拉伯地盾火成岩和变质岩基底的风化作用作为金属共同来源的推断相符，也得到低EF值的支持。阿拉伯地盾基岩富含多种金属矿物，如铬铁矿、钒磁铁矿、闪锌矿、黄铜矿、方铅矿、磁铁矿、赤铁矿、镍黄铁矿和毒砂等。但需强调，这些矿物描述代表阿拉伯地盾的广泛地质背景，并非本研究样品中金属的实际赋存矿物相，后者需要XRD、SEM-EDS或顺序提取等专门的矿物学和地球化学分析。HCA（Q模式）将32个样品分为两个簇：C1含13个样品，C2含19个样品。除Y3和Y7外，所有延布样品均归入C1，而全部胡莱斯样品加上Y3和Y7归入C2。Y3和Y7的Fe浓度异常高（分别为46,700和53,900 mg kg
?1），导致其与胡莱斯主导的簇归并。C1簇中除Pb外所有PTEs浓度最低。R模式HCA将Pb与其他七种金属分离，与其在相关矩阵中的独特行为一致。
3.6 Mann–Whitney U检验
两地比较显示，所有八种金属在两地间均存在显著差异（p < 0.05）。胡莱斯地区的Co、Cu、Ni、V全部高于延布（U = 0），Cr也几乎完全分离（U = 0.5）。Fe因Y3和Y7两个高值样品而存在部分重叠（U = 14.0）。Pb则相反，延布浓度一致高于胡莱斯（U = 1.5）。Zn分离程度最低（U = 63.5，p = 0.02），与其多种来源（地质成因和农业活动）相符。
3.7 污染状况评价
Igeo结果显示，几乎所有PTE在几乎所有点位均为实际无污染（Igeo < 0），唯一例外是Ni在四个胡莱斯点位（Kh5、Kh7、Kh13、Kh14）达到无污染至中度污染等级（Igeo ≥ 0）。EF平均值为Co (0.93)、Ni (0.92)、Cu (0.80)、Zn (0.74)、V (0.65)、Cr (0.59)、Pb (0.27)，均低于2，表明研究区土壤基本未富集或仅轻微富集。CF平均值由高到低为Fe (0.93) > Ni (0.92) > Co (0.89) > Cu (0.79) > Zn (0.65) > V (0.60) > Cr (0.54) > Pb (0.22)。尽管平均CF小于1，但在大多数胡莱斯点位，Co、Cu、Fe、Ni分别达到了1 ≤ CF < 3的“中度污染”等级。PLI平均值为0.61，范围0.35–0.84，表明土壤总体无污染。综合四项指数，研究区绝大多数农田处于安全水平，但胡莱斯部分点位已出现需关注的迹象。
四、讨论与结论
研究讨论指出，四项污染指数不仅可作为描述性指标，还可作为规划者和决策者在批准农业用地之前使用的低成本筛选工具。大多数点位CF、PLI和Igeo值低于富集或污染阈值，EF值均低于2，指示主导来源为地质成因而非人为输入，这对延布和胡莱斯大部分农田是令人放心的。然而，胡莱斯部分点位Co、Cu、Fe、Ni的CF超过1，Ni的Igeo达到无污染—中度污染等级，这些点位值得更密切监测。将基于指数的决策规则纳入常规农业土地利用规划，可在作物种植前预先识别风险地块，而非等问题进入食物链后才应对。干旱土壤有机物含量低、缓冲额外压力的能力有限，因此保持接近自然背景的化学状态需要积极保护。研究同时指出局限性：王水消解提取的是伪总浓度而非全量；样本量有限，属筛查级评估；区域背景值尚未建立，采用了国际通用的页岩平均背景。结论部分强调，PCA和HCA结果支持阿拉伯地盾基底岩石是这些金属的合理主要来源；研究区土壤PTEs富集程度低或无富集；CF和PLI表明土壤平均无污染，但个别点位需关注。综合污染指数与多元框架为规划和政策制定者提供了一个实用、低成本的环境保护与土地管理筛选工具。该研究建议定期监测土壤质量，并增加采样密度以提高空间分辨率。该框架具有可复制性，其他干旱、粮食依赖进口的地区可借鉴以保护有限的耕地资源。将此类基于指数的筛选纳入常规农业和环境规划，有助于实现合理的土地管理实践，特别对自然恢复力有限的干旱土壤而言，早期检测污染对于长期环境管理至关重要。该研究成果发表在《Sustainability》期刊，为可持续土地管理和粮食安全提供了有价值的地球化学基线。
