禹城地区土壤铅含量空间分布的指示克里格估值-地理学报.PDF

生态环境学报 2011, 20(12): 1912-1918 Ecology and Environmental Sciences E-mail: editor@ 禹城地区土壤铅含量空间分布的指示克里格估值 * 姚荣江，杨劲松 ，杨奇勇，黄标，孙维侠 中国科学院南京土壤研究所，江苏 南京 210008 摘要：以鲁西北禹城地区土壤重金属 Pb 含量为例，针对其存在特异值与偏态分布的特点，利用指示克里格法研究了土壤 Pb 含量小于特定阈值的条件概率分布，分析了该地区土壤 Pb 含量的空间分布特征，并对普通克里格和指示克里格的估值效果 进行交叉验证。结果表明：研究区土壤 Pb 含量总体较低，未超过国家标准的限值，但不同部位土壤 Pb 含量差异较大；普 通克里格法和指示克里格法形成的土壤 Pb 含量空间分布格局相似，但普通克里格法的平滑效果强于指示克里格法；研究区 域土壤环境质量总体良好，Pb 含量高值区主要集中在各乡镇驻地附近的农田；指示克里格对土壤 Pb 含量空间估值的精度优 于普通克里格，指示克里格法更适用于原始数据呈非正态分布且特异值对估值结果有较大影响的土壤属性。本研究结果为该 地区土壤环境质量评价和重金属污染防控提供科学依据。 关键词：禹城；铅；空间分布；指示克里格；条件概率 中图分类号：X53 文献标志码：A 文章编号：1674-5906 （2011 ）12-1912-07 随着城市化进程与工农业的迅速发展，土壤中 地的鲁西北禹城市为例，采用指示克里格法对该地 出现大量的重金属元素富集，这不仅会引起土壤甚 区土壤重金属 Pb 含量分布进行估值研究，揭示土 至耕地污染，还可通过食物链进入人体，危害人类 壤中 Pb 含量小于特定阈值的条件概率分布及其含 健康。因此，土壤中重金属的含量及其空间分布、 量空间分布特征，并对普通克里格和指示克里格的 环境质量、风险评价等研究都受到了人们的广泛关 估值效果进行交叉验证，旨在为该地区土壤环境质 注，再加上近年来 GIS 和地统计理论被引入这一研 量评价和重金属污染防控提供科学依据。 究中，推动了研究向前发展，目前这已成为土壤科 1 材料与方法 学的热点问题之一[1-2]。然而，人们在土壤重金属时 1.1 研究区概况 空变异研究中发现，由于受来源、形态及控制因素 禹城市地处山东省西北部，徒骇河中游，地理 等影响，区域土壤重金属往往存在特异值，这些特 位置介于 E 116°22′~ 116°45′，N 36°40′~ 37°12′，属 异值既非分析误差所致，也非采样方法等人为误差 暖温带半湿润季风气候，多年平均气温 13.1 ℃，降 所致，却往往使观测值偏离正态分布，影响变异函 雨量 616 mm ，蒸发量 1 379.5 mm，蒸降比 2.24 ， 数的稳健性[3-4]。目前解决途径主要有三种：一种是 降雨季节性分布不均，其中 6~9 月的夏季平均为 非线性变换（如平方根变换，对数变换，反正弦变 470.7 mm ，占全年的 76.4%。禹城市属典型的冲积 换等），但目前只有对数变换可以应用对数正态克 平原地质区，海拔高度介于 19.2~27.3 m，微地貌比 里格估值并转换回原尺度，其他变换没有相应的转 较复杂，现有耕地 5.05×104 ha ，占全市总面积的 换方法；另一个途径是剔除或用估计邻域和影响系 51.01% ，主要有潮土和潮化盐土两大土类，潮土占 数等方法处理特异值[5-6] ，但会改变实验数据的原始 总面积的 97.77%。历史上，禹城市为黄淮海