土壤Cd含量风险管理论文1材料与方法1.1土样采集为检测珠三角土壤元素变化，2022年建立了珠三角核心区的生态地球化学监控网络，共布设了98个监控点，以土壤为监控介质，以约160km2/点为监控点密度，土壤样点分布见图1。2022年对“监控网络”同点位表层土壤样品进行了采集。2022年和2022年，在监控点上分别采集了0~20cm的表层土壤样品。每个土壤样点均由4个副样混合组成，样品在常温下（20~24℃）风干，除去石子或其它杂物；土壤混合样过2mm的聚乙烯筛，然后在玛瑙研钵中研磨，过0.149mm筛[26]。为消除不同年份间的分析系统误差，2022年采集的“监控网络”表层土壤样品与2022年副样进行同批次分析。全部样品预处理后，在中国地质科学院廊坊物化探研究所进行分析测试。用电感耦合等离子体质谱仪(ICP-MS)测定了土壤重金属Cd的全量。分析过程中进行分析质量监控，分别插入8%的GSS-1、GSS-2、GSS-3和GSS-8等国家一级标准物质和5%的密码重复样，监控样的分析数据显示样品分析质量符合《生态地球化学评价样品分析技术要求（试行）》的要求[27]。1.2土壤重金属累积速率估算模型工业化程度越高，城市化发展速度越快，对土壤污染的“贡献”就越大，每年对土壤中污染物的叠加，不是等量的，对土壤的污染也不是等速的，而是加速的。因此，可认为土壤的污染过程可归纳为两个阶段：一是加速阶段，二是匀速阶段[29]。有研究表明：10a尺度土壤表层重金属元素积累符合线性关系[30~32]。根据研究区的工业化、城市化进程，产业结构以及人均GDP与环境的关系，2000年以后研究区土壤的污染过程进入匀速阶段[32]。因此，2000年以后土壤重金属的含量符合下式。1.3土壤重金属时空分布预测模型1）时空模式。对不同时间下变量Z叠加累积速率，预测不同时间下变量Z的空间分布。先将时间ti下位置xj的变量Z，叠加累积速率。2）空时模式。先对变量Z进行空间分布预测，再将各点位的变量Z叠加累积速率，预测变量Z在不同时间下的空间分布。2结果与讨论2.1土壤Cd和pH特征2022年和2022年土壤pH和Cd的描述性统计结果、两期数据的均值多重比较分别见表1、2。从表1可知，2022年和2022年土壤pH的平均值分别为5.65和6.02，属于微酸性土壤。两者均值多重比较结果表明（表2），pH的均值变化达到显著水平。从表2可知，土壤pH通过Levene检验，显著性为0.099，均值比较选用方差相等条件下的t检验，pH的检验结果拒绝H。假设，即从2022年至2022年研究区土壤pH发生了显著变化。总体上土壤pH有所升高，这与该区土地利用方式包含较多的建筑用地、居民用地、道路用地等，受石灰性填充物的影响较大有关。6a间土壤Cd的平均累积量不变（表1），但是两者均值多重比较结果表明（表2），Cd的均值变化达到显著水平。从表2可知，Cd没有通过Lev-ene检验，均值比较选用方差不相等条件下的t检验，检验结果拒绝H。假设，即从2022年至2022年该研究区土壤Cd发生了显著变化。2022年和2022年两期次土壤Cd平均含量（0.28mg/kg）是广东土壤背景值（0.06mg/kg）的5倍；与《土壤环境质量标准》（二级标准）相比较，2022年、2022年两期次土壤Cd平均含量与其相当，但是其最大值分别是《土壤环境质量标准》[33]（二级标准）的6倍和3.5倍，表明两期次土壤Cd含量在土壤中富集程度高，且有部分地区Cd含量超过了危害生态安全的阈值。变异系数在某种程度上定量刻画了数据的离散程度和变异性。在本研究中，2022年Cd的变异系数较2022年有降低，但变异水平仍属于强变异，暗示了土壤Cd主要来源于人为成因。以“土壤Cd”和“珠三角”为关键词，对CNKI数据库进行文献检索，得出2022年以来珠三角各类土壤中Cd含量（表3）。从表3可知，随着时间的推移，珠三角土壤Cd含量有累积趋势，佛山比较突出。2.2土壤Cd累积速率估算采用土壤重金属累积速率估算模型对研究区各地市土壤Cd进行累积速率估算（表4）。从表4可知，区域性土壤重金属Cd的年均累积速率KC2022-C2022为0.016mg/kg，对比研究区2022年与2022年珠三角工业企业周边蔬菜土壤Cd含量，土壤Cd的年均增长率KC2022-C2022为0.014mg/kg，为成都经济区土壤Cd累积速率的2.7倍[40],这可能与当地经济发展程度有一定联系。其中，东莞、...