第３２卷 第７期２０１９年７月环 境 科 学 研 究Ｒｅｓｅａｒｃｈ ｏｆ Ｅｎｖｉｒｏｎｍｅｎｔａｌ ＳｃｉｅｎｃｅｓＶｏｌ.３２ꎬＮｏ.７Ｊｕｌｙꎬ２０１９收稿日期 ２０１８￣０６￣０４ 修订日期 ２０１９￣０２￣１９作者简介李彤１９９５￣ꎬ女ꎬ山东潍坊人ꎬｌｉｔｏｎｇ１９９５０８０２＠ １２６.ｃｏｍ.∗责任作者ꎬ吴泉源１９５９￣ꎬ男ꎬ山东龙口人ꎬ教授ꎬ博士ꎬ博导ꎬ主要从事地图学与地理信息系统研究ꎬｗｑｙ６４２０５８２＠ １６３.ｃｏｍ基金项目山东省自然科学基金项目Ｎｏ.ＺＲ２０１７ＢＤ０１１ꎻ中国博士后科学基金项目Ｎｏ.２０１７Ｍ６２２２５６Ｓｕｐｐｏｒｔｅｄ ｂｙ Ｎａｔｕｒａｌ Ｓｃｉｅｎｃｅ Ｆｏｕｎｄａｔｉｏｎ ｏｆ Ｓｈａｎｄｏｎｇ Ｐｒｏｖｉｎｃｅꎬ Ｃｈｉｎａ Ｎｏ.ＺＲ２０１７ＢＤ０１１ꎻ Ｃｈｉｎａ Ｐｏｓｔｄｏｃｔｏｒａｌ Ｓｃｉｅｎｃｅ ＦｏｕｎｄａｔｉｏｎＮｏ.２０１７Ｍ６２２２５６龙口市北部平原区不同土地利用类型对土壤重金属污染距离的影响李 彤ꎬ吴泉源∗ ꎬ姚 磊ꎬ曹见飞ꎬ刘玉红ꎬ宋仁超山东师范大学地理与环境学院ꎬ山东济南 ２５０３５８摘要为探究不同土地利用类型及其组成对土壤重金属污染距离的影响ꎬ以龙口市北部平原作为典型研究区域ꎬ系统采集５５个表层土壤样品并对ｗＰｂ、ｗＣｄ进行了测定ꎻ通过建立基于土壤样区的梯度缓冲区ꎬ利用Ｐｅａｒｓｏｎ相关性分析ꎬ研究不同土地利用类型及其组成对土壤重金属污染距离的影响ꎻ通过热点分析和空间插值方法ꎬ探究研究区内对Ｐｂ、Ｃｄ富集贡献较大的土地利用类型.结果表明①龙口市北部平原土壤中Ｐｂ、Ｃｄ受人类活动影响的区域差异明显ꎬｗＰｂ、ｗＣｄ平均值均超过了山东省土壤相应背景值ꎬ其中ꎬｗＰｂ超过背景值的样品比例达到９６􀆰 ３６％ꎬ研究区西部富集现象明显. ②土地利用类型面积占比对Ｐｂ、Ｃｄ的污染距离各异.其中ꎬ交通用地面积占比对Ｐｂ、Ｃｄ的污染距离大约在２ ０００、５００ ｍ处ꎬ建设用地面积占比对二者的污染距离分别在１ ０００、２００ ｍ处左右ꎬ而耕地及园地面积占比对二者的污染距离较小ꎬ均在５０ ｍ处左右. ③ ｗＰｂ、ｗＣｄ的空间分布差异明显ꎬ整体表现为西高东低. ④研究区内对Ｐｂ、Ｃｄ富集贡献最大的土地利用类型是耕地ꎬ其次是交通和建设用地.研究显示ꎬ不同土地利用类型对重金属污染距离的差异性是实际存在的ꎬ因此在进行重金属污染治理时ꎬ对不同的土地利用类型要根据其对重金属污染距离的差异加以区别.关键词土地利用类型ꎻ土壤ꎻ重金属ꎻ相关性中图分类号 Ｘ５３ 文章编号 １００１￣６９２９２０１９０７￣１２２４￣０７文献标志码 Ａ ＤＯＩ １０􀆰 １３１９８∕ｊ􀆰 ｉｓｓｎ􀆰 １００１￣６９２９􀆰 ２０１９􀆰 ０２􀆰 １５Ｅｆｆｅｃｔｓ ｏｆ Ｄｉｆｆｅｒｅｎｔ Ｌａｎｄ Ｕｓｅ Ｔｙｐｅｓ ｏｎ Ｐｏｌｌｕｔｉｏｎ Ｄｉｓｔａｎｃｅ ｏｆ Ｈｅａｖｙ Ｍｅｔａｌｓ ｉｎＮｏｒｔｈ Ｐｌａｉｎ ｏｆ Ｌｏｎｇｋｏｕ ＣｉｔｙＬＩ Ｔｏｎｇꎬ ＷＵ Ｑｕａｎｙｕａｎ∗ ꎬ ＹＡＯ Ｌｅｉꎬ ＣＡＯ Ｊｉａｎｆｅｉꎬ ＬＩＵ Ｙｕｈｏｎｇꎬ ＳＯＮＧ ＲｅｎｃｈａｏＣｏｌｌｅｇｅ ｏｆ Ｇｅｏｇｒａｐｈｙ ａｎｄ Ｅｎｖｉｒｏｎｍｅｎｔꎬ Ｓｈａｎｄｏｎｇ Ｎｏｒｍａｌ Ｕｎｉｖｅｒｓｉｔｙꎬ Ｊｉｎａｎ ２５０３５８ꎬ ＣｈｉｎａＡｂｓｔｒａｃｔ Ｉｎ ｏｒｄｅｒ ｔｏ ｅｘｐｌｏｒｅ ｔｈｅ ｅｆｆｅｃｔｓ ｏｆ ｄｉｆｆｅｒｅｎｃｅ ｏｆ ｌａｎｄ ｕｓｅ ｔｙｐｅ ａｎｄ ｉｔｓ ｃｏｍｐｏｓｉｔｉｏｎ ｏｎ ｓｏｉｌ ｈｅａｖｙ ｍｅｔａｌ ｐｏｌｌｕｔｉｏｎ ｄｉｓｔａｎｃｅꎬ ｔｈｅｎｏｒｔｈｅｒｎ ｐｌａｉｎ ｏｆ Ｌｏｎｇｋｏｕ Ｃｉｔｙ ｗａｓ ｔａｋｅｎ ａｓ ａ ｓｔｕｄｙ ａｒｅａ􀆰 Ｆｉｆｔｙ￣ｆｉｖｅ ｓｕｒｆａｃｅ ｓｏｉｌ ｓａｍｐｌｅｓ ｗｅｒｅ ｃｏｌｌｅｃｔｅｄ ａｎｄ Ｐｂ ａｎｄ Ｃｄ ｃｏｎｔｅｎｔｓ ｗｅｒｅｄｅｔｅｒｍｉｎｅｄ. Ｂｙ ｅｓｔａｂｌｉｓｈｉｎｇ ｇｒａｄｉｅｎｔ ｂｕｆｆｅｒ ｚｏｎｅ ｂａｓｅｄ ｏｎ ｓｏｉｌ ｓａｍｐｌｅ ａｒｅａꎬ ｔｈｅ ｉｎｆｌｕｅｎｃｅ ｏｆ ｌａｎｄ ｕｓｅ ｔｙｐｅ ａｎｄ ｉｔｓ ｃｏｍｐｏｓｉｔｉｏｎ ｏｎ ｓｏｉｌ ｈｅａｖｙｍｅｔａｌ ｐｏｌｌｕｔｉｏｎ ｕｎｄｅｒ ｄｉｆｆｅｒｅｎｔ ｇｒａｄｉｅｎｔｓ ｗａｓ ａｎａｌｙｚｅｄ. Ｔｈｒｏｕｇｈ ｔｈｅ ｈｏｔ ｓｐｏｔ ａｎａｌｙｓｉｓ ａｎｄ ｓｐａｔｉａｌ ｉｎｔｅｒｐｏｌａｔｉｏｎ ｍｅｔｈｏｄꎬ ｔｈｅ ｌａｎｄ ｕｓｅ ｔｙｐｅｓ ｉｎｔｈｅ ｓｔｕｄｙ ａｒｅａ ｔｈａｔ ｃｏｎｔｒｉｂｕｔｅ ｇｒｅａｔｌｙ ｔｏ Ｐｂ ａｎｄ Ｃｄ ｅｎｒｉｃｈｍｅｎｔ ｗａｓ ｅｘｐｌｏｒｅｄ. Ｔｈｅ ｒｅｓｕｌｔｓ ｓｈｏｗｅｄ ｔｈａｔ １ Ｐｂ ａｎｄ Ｃｄ ｉｎ ｔｈｅ ｓｏｉｌ ｏｆ ｔｈｅｓｔｕｄｙ ａｒｅａ ｗｅｒｅ ｓｉｇｎｉｆｉｃａｎｔｌｙ ａｆｆｅｃｔｅｄ ｂｙ ｈｕｍａｎ ａｃｔｉｖｉｔｉｅｓꎬ ａｎｄ ｔｈｅ ａｖｅｒａｇｅ ｖａｌｕｅ ｅｘｃｅｅｄｅｄ ｔｈｅ ｂａｃｋｇｒｏｕｎｄ ｖａｌｕｅ ｏｆ ｓｏｉｌ ｉｎ ＳｈａｎｄｏｎｇＰｒｏｖｉｎｃｅ. Ｔｈｅ ｐｒｏｐｏｒｔｉｏｎ ｏｆ Ｐｂ ｅｘｃｅｅｄｉｎｇ ｔｈｅ ｂａｃｋｇｒｏｕｎｄ ｖａｌｕｅ ｒｅａｃｈｅｄ ９６􀆰 ３６％ꎬ ａｎｄ ｔｈｅ ｗｅｓｔｅｒｎ ｅｎｒｉｃｈｍｅｎｔ ｐｈｅｎｏｍｅｎｏｎ ｉｓ ｏｂｖｉｏｕｓ.２ Ｉｎ ｄｉｆｆｅｒｅｎｔ ｇｒａｄｉｅｎｔｓꎬ ｔｈｅ ｐｅｒｃｅｎｔａｇｅ ｏｆ ｌａｎｄ ｕｓｅ ｔｙｐｅｓ ｈａｄ ｄｉｆｆｅｒｅｎｔ ｃｏｒｒｅｌａｔｉｏｎｓ ｗｉｔｈ Ｐｂ ａｎｄ Ｃｄ. Ｔｈｅ ｍａｘｉｍｕｍ ｃｏｒｒｅｌａｔｉｏｎ ｃｏｅｆｆｉｃｉｅｎｔｏｆ ｔｒａｆｆｉｃ ｌａｎｄ ｗｉｔｈ Ｐｂ ａｎｄ Ｃｄ ｒｅａｃｈｅｄ ２０００ ａｎｄ ５００ ｍ ｒｅｓｐｅｃｔｉｖｅｌｙꎬ ａｎｄ ｔｈｅ ｃｏｎｓｔｒｕｃｔｉｏｎ ｌａｎｄ ｗａｓ １０００ ａｎｄ ２００ ｍ. Ｔｈｅ ｃｏｒｒｅｌａｔｉｏｎｃｏｅｆｆｉｃｉｅｎｔ ｏｆ ｃｕｌｔｉｖａｔｅｄ ｌａｎｄ ａｎｄ ｇａｒｄｅｎ ｌａｎｄ ｗａｓ ｔｈｅ ｌａｒｇｅｓｔ ａｔ ｔｈｅ ｇｒａｄｉｅｎｔ ｏｆ ａｂｏｕｔ ５０ ｍ. ３ Ｔｈｅ ｓｐａｔｉａｌ ｄｉｓｔｒｉｂｕｔｉｏｎ ｄｉｆｆｅｒｅｎｃｅ ｏｆ Ｐｂａｎｄ Ｃｄ ｗａｓ ｏｂｖｉｏｕｓꎬ ａｎｄ ｔｈｅ ｏｖｅｒａｌｌ ｐｅｒｆｏｒｍａｎｃｅ ｗａｓ ｈｉｇｈ ｉｎ ｔｈｅ ｗｅｓｔ ａｎｄ ｌｏｗ ｉｎ ｔｈｅ ｅａｓｔ. ４ Ｔｈｅ ｃｕｌｔｉｖａｔｅｄ ｌａｎｄ ｉｎ ｔｈｅ ｓｔｕｄｙ ａｒｅａｃｏｎｔｒｉｂｕｔｅｄ ｔｈｅ ｍｏｓｔ ｔｏ Ｐｂ ａｎｄ Ｃｄ ｅｎｒｉｃｈｍｅｎｔꎬ ｆｏｌｌｏｗｅｄ ｂｙ ｔｒａｆｆｉｃ ａｎｄ ｃｏｎｓｔｒｕｃｔｉｏｎ ｌａｎｄ. Ｔｈｅ ｓｔｕｄｙ ｓｈｏｗｓ ｔｈａｔ ｄｉｆｆｅｒｅｎｔ ｌａｎｄ ｕｓｅ ｔｙｐｅｓ ｈａｖｅｄｉｆｆｅｒｅｎｔ ｈｅａｖｙ ｍｅｔａｌ ｐｏｌｌｕｔｉｏｎ ｄｉｓｔａｎｃｅｓ. Ｉｎ ｔｈｅ ｔｒｅａｔｍｅｎｔ ｏｆ ｈｅａｖｙ ｍｅｔａｌ ｐｏｌｌｕｔｉｏｎꎬ ｄｉｆｆｅｒｅｎｔ ｔｙｐｅｓ ｏｆ ｌａｎｄ ｕｓｅ ｓｈｏｕｌｄ ｂｅ ｄｉｓｔｉｎｇｕｉｓｈｅｄａｃｃｏｒｄｉｎｇ ｔｏ ｔｈｅ ｄｉｆｆｅｒｅｎｃｅ ｏｆ ｔｈｅ ｐｏｌｌｕｔｉｏｎ ｄｉｓｔａｎｃｅ.Ｋｅｙｗｏｒｄｓ ｌａｎｄ ｕｓｅ ｔｙｐｅｓꎻ ｓｏｉｌꎻ ｈｅａｖｙ ｍｅｔａｌꎻ ｃｏｒｒｅｌａｔｉｏｎ第７期李 彤等龙口市北部平原区不同土地利用类型对土壤重金属污染距离的影响 随着城镇化的发展ꎬ林地、草地等土地利用类型开始向耕地、建设用地等人为利用方面转变ꎬ交通运输、污水灌溉等人类活动直接或间接地提高了土壤重金属污染水平[１￣２].重金属污染直接影响到农作物及饮用水的安全性ꎬ进而威胁到人类的身体健康[３].目前ꎬ重金属污染已经成为世界性的环境问题ꎬ引起了学术界的广泛关注.影响土壤重金属分布和累积的因素可分为两类①自然因素ꎬ如成土母质和土壤理化性质ꎻ②人类活动ꎬ如交通、矿产开采、农药化肥、污水灌溉及工业废弃物等.重金属污染则主要是由于人类活动导致的ꎬ而土地利用类型是人类活动的集中体现.目前ꎬ对土壤重金属污染的研究主要侧重于重金属的来源解析[４￣５]、空间分布[６￣７]、迁移转化[８￣９]以及污染综合评价[１０￣１４]等方面ꎬ重点关注土壤重金属污染现状及其对环境的潜在风险ꎬ或仅局限于不同土地利用类型之间重金属含量的比较.基于不同梯度来探究土地利用类型及其结构组成与土壤重金属关系的研究较少.实际上ꎬ在不同的距离梯度上ꎬ土地利用类型与重金属含量的相关性是不同的ꎬ如王祖伟等[１５￣１６]发现ꎬ重金属含量与缓冲区内土地利用类型面积占比的相关性较强ꎬ且随梯度的增大而发生改变.可见ꎬ明确不同梯度上对重金属含量影响较大的土地利用类型以及某种土地利用类型在哪个梯度上对于缓解重金属污染更为重要ꎬ对土壤重金属管理有重要的现实意义.龙口市是２０１７年全国综合经济竞争力十强县市之一ꎬ作为山东省入选的两个县市之一ꎬ该地区尤其是北部平原区社会经济活动发达ꎬ近年来发展迅速ꎬ在工业、农业等方面均表现出巨大的发展潜力.在社会经济增长的同时ꎬ土壤重金属含量也在不断攀升ꎬ以致部分区域土壤Ｐｂ、Ｃｄ富集现象明显[１７]ꎻ且北部平原区地质背景单一ꎬ成土条件相似ꎬ重金属之间的差异主要是不同土地利用类型引起的ꎬ而不必考虑自然条件对土壤中重金属影响的差异.因此ꎬ该研究选择龙口市北部平原区作为典型研究区ꎬ基于实地采集的土壤样品ꎬ利用ＧＩＳ方法构建以采样点为中心的不同间距的缓冲区作为分析单元ꎬ分析不同梯度缓冲区内对Ｐｂ、Ｃｄ富集起主要影响的土地利用类型ꎬ以期为环境保护部门评价当地土壤环境质量治理和修复提供科学指导.１ 数据与方法１􀆰 １ 研究区概况研究区设在山东省烟台市的龙口市１２０１３′Ｅ１２０４４′Ｅ、３７２７′Ｎ３７４４′Ｎ北部沿海平原ꎬ南北长约１８ ｋｍꎬ东西宽约２２ ｋｍꎬ总面积约４００ ｋｍ２.从地质背景来看ꎬ研究区三面环海ꎬ成土条件相似ꎬ地质背景整体上具有一致性.成土母质以冲积和海积物等第四系沉积物为主ꎬ多形成潮土.南部发育有少量白垩系的火山熔岩和上元古界的变质岩ꎬ岩石风化多形成棕壤ꎬ东部则有少量喜马拉雅期形成的玄武岩.从气候条件来看ꎬ研究区属于温带季风性气候ꎬ年均气温约１２ ℃ꎬ年均降水量６００ ｍｍꎬ且主要集中在夏季.河流主要有泳汶河、黄水河等ꎬ多为季节性河流.地势平坦ꎬ适宜农业生产ꎬ盛产农作物主要有小麦、玉米、苹果、葡萄等.近年来ꎬ龙口市城镇化和工业化速度较快ꎬ特别是机械加工业的发展ꎬ且矿产资源较丰富ꎬ其中ꎬ北部平原区以煤矿为主.长期的煤炭开采ꎬ不仅出现大面积采空区塌陷问题ꎬ其固废堆放、废液排放也可能对周围土壤重金属环境造成威胁.此外ꎬ龙口市北部平原是胶东地区重要的粮食生产基地ꎬ污水灌溉历史悠久ꎬ且农作物是当地农村经济发展的支柱ꎬ农民为提高产量而盲目施肥ꎬ也可能使得该区土壤面临重金属污染的风险.因此ꎬ该研究选择龙口市北部平原为研究区ꎬ并根据龙口市国土资源局提供的高分影像数据２４小类ꎬ将其整合为六大类 耕地、园地、林草地、交通用地、水域和建设用地见图１ꎬ其中ꎬ由于林草地所占面积极少、水域中重金属含量不明ꎬ故这两类不考虑在内.１􀆰 ２ 土壤样品采集与重金属含量测定结合研究区土地利用类型现状ꎬ采用蛇形布点法在ＡｒｃＧＩＳ １０􀆰 ２软件中初步确定采样点的位置.其中ꎬ采样点间距基本控制在２ ｋｍ左右ꎬ并使其覆盖所有的土地利用类型.同时ꎬ在研究区内典型地物如电热厂、大型工厂等周围也设定采样点ꎬ使得采样点具有典型性ꎬ能够代表研究区内不同土地利用类型的重金属含量水平.在实际采样过程中ꎬ根据预设采样点的周围环境适当进行调整.通过了解研究区地块的平均宽度ꎬ使样方尽可能覆盖地块而更具有代表性ꎬ最终确定样方大小为２０ ｍ２０ ｍ.由于样方面积较小ꎬ且样方内地形变化小、土壤质地均匀ꎬ因此ꎬ在样方内采用梅花形布点取样法ꎬ选取样方４个顶点及中心点作为采样点坐标位置的土壤ꎬ清理土壤表层ꎬ垂直采集表层０２０ ｃｍ土壤并均匀混合ꎬ避免出现异常高值[１７].采用四分法取１ ｋｇ土壤装入干净的布袋并做好标记.记录采样点的地理位置及土地利用类型ꎬ土壤样品采样点的分布如图１所示.在实验室内将土壤样品于室温下自然风干、敲碎５２２１环 境 科 学 研 究第３２卷图１ 研究区土地利用现状及采样点分布Ｆｉｇ.１ Ｍａｐ ｏｆ ｌａｎｄ ｕｓｅ ｔｙｐｅ ａｎｄ ｓａｍｐｌｉｎｇ ｐｏｉｎｔｓ ｉｎ ｔｈｅ ｓｔｕｄｙ ａｒｅａ后ꎬ置于干净的平台上挑选出植物根系、砾石等杂物ꎬ在无污染的玛瑙研钵内研磨样品后ꎬ过１００目０􀆰 ２５ｍｍ尼龙筛.选用ＨＮＯ３￣ＨＣｌ￣ＨＣｌＯ４微波法消解土壤样品ꎬ消解结束后进行蒸发ꎬ在接近干燥时冷却至常温常压ꎬ用超纯水溶解土样ꎬ置于５０ ｍＬ容量瓶中定容. ｗＰｂ、ｗＣｄ采用石墨原子吸收分光光度法测定.采用国家一级土壤标准物质检验分析方法的精密度和准确度ꎬ以保证测试结果符合精度需要[４].１􀆰 ３ 研究方法１􀆰 ３􀆰 １ 污染距离分析研究[１５]表明ꎬ重金属含量与不同梯度下土地利用类型结构组成的相关性存在较大差异.为了更好地体现不同土地利用类型对土壤重金属污染距离的差异ꎬ以各土壤采样点为圆心ꎬ利用ＡｒｃＧＩＳ １０􀆰 ２软件的“缓冲区”工具分别建立５０、１００、１５０、２００、５００、１ ０００及２ ０００ ｍ间距的缓冲区.分别计算缓冲区内的土地利用类型面积占比ꎬ即基于研究区内土地利用类型划分ꎬ计算每一种类型的面积占缓冲区总面积的比例.以采样点ｗＰｂ、ｗＣｄ与各缓冲区内土地利用类型面积占比作为分析变量进行Ｐｅａｒｓｏｎ相关性分析ꎬ通过这些由近及远、由密到疏的缓冲区以及相关系数随梯度缓冲区半径增大的变化方向ꎬ得出不同土地利用类型面积占比对Ｐｂ、Ｃｄ影响最强时对应的污染距离.１􀆰 ３􀆰 ２ 热点分析与单因素方差分析基于热点分析ꎬ利用单因素方差分析方法ꎬ探究冷点、热点在关键污染距离内的土地利用类型面积之间的差异ꎬ识别研究区内对Ｐｂ、Ｃｄ富集贡献较大的土地利用类型.其中ꎬ方差分析对象选择显著性小于或等于０􀆰 ０５的热点和冷点.热点分析ｇｅｔｉｓ￣ｏｒｄ Ｇｉ∗ 是一种局部空间自相关分析方法.在全局空间集聚的前提下ꎬ通过计算Ｇｉ∗统计量来揭示集聚模式并反映高值或低值要素在空间上发生集聚的位置.显著的Ｇｉ∗正值表示高观测值的区域单元趋于空间集聚ꎬ即热点ｈｏｔｓｐｏｔꎻ显著的Ｇｉ∗负值则表示低观测值的区域单元趋于空间集聚ꎬ即冷点ｃｏｌｄｓｐｏｔ.１􀆰 ４ 数据处理在ＡｒｃＧＩＳ １０􀆰 ２中完成梯度缓冲区的建立、不同缓冲区内土地利用类型面积汇总、热点分析以及普通克里金插值.对数据进行平均值、标准差、变异系数等描述性统计ꎬＰｅａｒｓｏｎ相关性分析以及单因素方差分析等统计分析均在ＳＰＳＳ １７􀆰 ０中完成.２ 结果与分析２􀆰 １ 土壤重金属含量的描述性统计由表１可见ꎬ研究区土壤中ｗＰｂ、ｗＣｄ的平均值分别为３１􀆰 ９４、０􀆰 １２ ｍｇ∕ｋｇꎬ均超过山东省土壤相应背景值[１８]分别是背景值的１􀆰 ２６、１􀆰 ０７倍ꎬ但均未超过ＧＢ １５６１８ １９９５土壤环境质量标准[１９]二级标准限值.其中ꎬｗＰｂ超过山东省土壤背景值的样品比例达到９６􀆰 ３６％. Ｗｉｌｄｉｎｇ[２０]对变异系数ｃｏｅｆｆｉｃｉｅｎｔ ｏｆ ｖａｒｉａｔｉｏｎꎬＣＶ进行了分级ꎬＣＶ０􀆰 ３６为高度变异.该研究中ｗＰｂ的变异系数为１７􀆰 ３４％ꎬ属于中度变异ꎬ而ｗＣｄ的变异系数高达７９􀆰 ５５％ꎬ６２２１第７期李 彤等龙口市北部平原区不同土地利用类型对土壤重金属污染距离的影响 表１ 研究区土壤中ｗＰｂ、ｗＣｄ的描述性统计分析Ｔａｂｌｅ １ Ｄｅｓｃｒｉｐｔｉｖｅ ｓｔａｔｉｓｔｉｃａｌ ａｎａｌｙｓｉｓ ｏｆ ｗＰｂ ａｎｄ ｗＣｄ ｉｎ ｓｏｉｌｓ ｉｎ ｔｈｅ ｓｔｕｄｙ ａｒｅａ重金属ｗ∕ｍｇ∕ｋｇ最小值最大值平均值标准差山东省背景值ＧＢ １５６１８ １９９５二级标准限值变异系数∕％超过背景值的样品比例∕％Ｐｂ ２３􀆰 ２３ ５５􀆰 ３３ ３１􀆰 ９４ ５􀆰 ５４ ２５􀆰 ４０ ３００ １７􀆰 ３４ ９６􀆰 ３６Ｃｄ ０􀆰 ０３ ０􀆰 ５４ ０􀆰 １２ ０􀆰 ０９ ０􀆰 １０８ ０􀆰 ３０ ７９􀆰 ５５ ３８􀆰 １８属于高度变异.这与ＬＶ等[４]对山东日照土壤重金属的研究结果相似ꎬ说明Ｃｄ易受交通、建筑施工等人为活动的影响ꎬ并在空间分布上也存在较大差异.２􀆰 ２ 土地利用类型对土壤重金属的污染距离分析Ｐｅａｒｓｏｎ相关性分析结果见表２、３表明ꎬ除园地面积占比与ｗＰｂ、ｗＣｄ呈负相关外ꎬ其他土地利用类型面积占比均与ｗＰｂ、ｗＣｄ呈正相关ꎬ且交通用地和建设用地面积占比与ｗＰｂ、ｗＣｄ的相关系数及对应梯度均略高于其他土地利用类型.表２ 不同梯度缓冲区内土地利用类型面积占比与ｗＰｂ的相关性Ｔａｂｌｅ ２ Ｔｈｅ ｃｏｒｒｅｌａｔｉｏｎ ｃｏｅｆｆｉｃｉｅｎｔ ｏｆ ｌａｎｄｓｃａｐｅｐｅｒｃｅｎｔａｇｅ ｗｉｔｈ ｗＰｂ ｉｎ ｅａｃｈ ｂｕｆｆｅｒ ｚｏｎｅ土地利用类型梯度缓冲区∕ｍ５０ １００ １５０ ２００ ５００ １ ０００ ２ ０００耕地０􀆰 ２０ ０􀆰 ０９ ０􀆰 ０８ ０􀆰 ０９ ０􀆰 １０ ０􀆰 １１ ０􀆰 １６园地－０􀆰 １４ －０􀆰 １２ －０􀆰 １０ －０􀆰 ０８ －０􀆰 ０５ －０􀆰 １１ －０􀆰 １３交通用地０􀆰 ０２ ０􀆰 ０２ ０􀆰 ０６ ０􀆰 １０ ０􀆰 １３ ０􀆰 １４ ０􀆰 １７建设用地０􀆰 １０ ０􀆰 １１ ０􀆰 １１ ０􀆰 １５ ０􀆰 １５ ０􀆰 ２８ ０􀆰 ０５表３ 不同梯度缓冲区内土地利用类型面积占比与ｗＣｄ的相关性Ｔａｂｌｅ ３ Ｔｈｅ ｃｏｒｒｅｌａｔｉｏｎ ｃｏｅｆｆｉｃｉｅｎｔ ｏｆ ｌａｎｄｓｃａｐｅｐｅｒｃｅｎｔａｇｅ ｗｉｔｈ ｗＣｄ ｉｎ ｅａｃｈ ｂｕｆｆｅｒ ｚｏｎｅ土地利用类型梯度缓冲区∕ｍ５０ １００ １５０ ２００ ５００ １ ０００ ２ ０００耕地０􀆰 ２４ ０􀆰 １３ ０􀆰 １２ ０􀆰 １７ ０􀆰 ２３ ０􀆰 ２０ ０􀆰 ２６园地－０􀆰 １３ －０􀆰 １２ －０􀆰 １１ －０􀆰 ０９ －０􀆰 ０６ －０􀆰 １４ －０􀆰 ２１交通用地０􀆰 ０３ ０􀆰 ０６ ０􀆰 １２ ０􀆰 １６ ０􀆰 １９ ０􀆰 １５ ０􀆰 ０９建设用地０􀆰 １４ ０􀆰 １４ ０􀆰 １４ ０􀆰 １６ ０􀆰 １１ ０􀆰 １２ ０􀆰 ０７由表２可见ꎬ在０２ ０００ ｍ梯度范围内ꎬ交通用地和建设用地面积占比与ｗＰｂ的相关系数略高于其他土地利用类型ꎬ最大值分别出现在２ ０００和１ ０００ｍ梯度处ꎬ且在０１ ０００ ｍ梯度范围内均呈上升趋势.相对于交通用地及建设用地ꎬ耕地对Ｐｂ的污染范围较小ꎬ最大相关系数出现在５０ ｍ梯度处ꎬ超过５０ ｍ之后呈下降趋势ꎬ并在１５０ ｍ处出现拐点.园地对Ｐｂ的污染距离与耕地相似ꎬ在５０ ｍ处相关系数最小ꎬ５０５００ ｍ梯度内呈上升趋势ꎬ相关系数逐渐变大ꎬ然后在５００ ｍ处出现拐点ꎬ相关系数减小.ｗＣｄ与各土地利用类型面积占比的相关系数变化趋势见表３与ｗＰｂ基本相似ꎬ这与戴彬等[２１￣２２]得出的结论相同ꎬ他们发现ｗＰｂ与ｗＣｄ呈显著相关ꎬ且有共同的污染源.交通用地和建设用地对Ｃｄ的污染距离也很大ꎬ最大相关系数对应的污染距离分别趋向于５００、２００ ｍ.与ｗＰｂ相同ꎬ园地面积占比与ｗＣｄ呈负相关ꎬ极小值出现在５０ ｍ处.随着梯度缓冲区半径的增大ꎬ二者的相关性逐渐减小ꎬ并在５００ ｍ处出现拐点.耕地面积占比与ｗＣｄ的相关系数在５０２ ０００ ｍ梯度内出现的拐点较多ꎬ但整体呈下降趋势ꎬ在５０ ｍ处出现极大值.２􀆰 ３ 土地利用类型对土壤重金属的贡献研究２􀆰 ３􀆰 １ 土壤重金属含量的热点分析分别对研究区５５个采样点ｗＰｂ、ｗＣｄ的Ｇｉ∗进行计算ꎬ得出不同显著性水平下的热点和冷点ꎬ结果如图２所示.热点区域是重金属的高值聚集区ꎬ具有较高的污染风险ꎬ是土壤环境管理、修复的重点区域.由图２可见ꎬＰｂ、Ｃｄ的热点数量分别为２５、１８个ꎬ其中ꎬＰｂ的热点较多ꎬ且在０􀆰 ０１显著性水平下的热点为１０个ꎬ占比为４０％ꎻＣｄ的热点多分布在研究区西部ꎬ冷点集中分布在东北部.对研究区土壤中ｗＰｂ、ｗＣｄ进行普通克里金插值分析见图２ꎬ研究其空间分布情况.由图２可见ꎬ冷点、热点分布与重金属含量插值结果的空间分布趋势大体一致ꎬ热点及周边地区均为重金属富集明显的区域ꎬ而冷点区均为ｗＰｂ、ｗＣｄ较低的区域.２􀆰 ３􀆰 ２ 影响土壤重金属富集的关键土地利用类型分析各土地利用类型面积占比与重金属含量的最大相关系数对应的污染距离大多未超过１ ０００ ｍꎬ这与Ｂｉａｓｉｏｌｉ等[２３]得出的各种土地利用类型产生的污染物传输只能限制在大约１ ｋｍ范围内的结论一致ꎬ这均说明采样点重金属含量受到周围１ ０００ ｍ左右的７２２１环 境 科 学 研 究第３２卷图２ 研究区Ｐｂ、Ｃｄ的热点分析及空间分布特征Ｆｉｇ.２ Ｈｏｔ ｓｐｏｔ ａｎａｌｙｓｉｓ ａｎｄ ｓｐａｔｉａｌ ｄｉｓｔｒｉｂｕｔｉｏｎ ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ ｏｆ Ｐｂ ａｎｄ Ｃｄ土地利用类型的影响.因此ꎬ该研究选择半径为１ ０００ｍ的缓冲区进行Ｃｄ冷点和热点的单因素方差分析Ｐｂ的热点分析中未产生冷点ꎬ无法与热点进行对比分析.由图３可见ꎬ就平均值而言ꎬＣｄ热点区域的耕地、建设用地面积明显大于冷点区域ꎬ而冷点区域的园地面积相对较大.方差分析结果显示ꎬ耕地、园地面积在统计意义上有显著差异ꎬ而交通用地、建设用地面积没有统计意义上的差异.这说明在研究区内ꎬＣｄ冷点和热点区域在土地利用类型组成上的区别主要体现在耕地和园地上. Ｃｄ热点区域内耕地面积明显高于冷点区域ꎬ且园地面积与ｗＣｄ呈负相关关系见表３ꎬ这在一定程度上能够缓解土壤中重金属的富集.可见ꎬ土壤采样点成为Ｃｄ热点的主要原因是耕地这一土地利用类型对土壤重金属的污染ꎬ其次是交通用地和建设用地.注缓冲区半径为１ ０００ ｍ.图３ Ｃｄ冷点、热点区域的土地利用类型平均面积对比Ｆｉｇ.３ Ｃｏｍｐａｒｉｓｏｎ ｏｆ ｌａｎｄ ｕｓｅ ｔｙｐｅｓ ｉｎｃｏｌｄ ａｎｄ ｈｏｔ ｓｐｏｔｓ ｏｆ Ｃｄ３ 讨论通过Ｃｄ冷点和热点的单因素方差分析ꎬ得出热点出现的主要原因是耕地受到污染. ２０１６年烟台市统计年鉴显示ꎬ２０１５年龙口市耕地有效灌溉面积为１７ ３８０ ｈｍ２ꎬ而农作物播种总面积达２４ ７１０ ｈｍ２ꎬ水库及机电井供水量并不能满足耕种需求ꎬ且研究区内河流多为季节性河流ꎬ在耕种早期严重缺水ꎬ故引泳汶河、黄水河及各支流中的工业、生活污水进行灌溉.实际上ꎬ龙口市灌区农用地的污水灌溉历史已有３０ ａ之久.未经处理的工业和生活污水直接灌溉农田ꎬ给土壤环境造成很大压力.李春芳等[１７]也通过溯源分析发现ꎬ污水灌溉是龙口市北部平原区Ｐｂ、Ｃｄ富集的共同污染因素.可见ꎬ长期的污水灌溉是导致耕地土壤重金属富集的原因之一.当地农民为了提高种植产量而盲目施肥施药ꎬ造成重金属在土壤中的富集也是一个不可忽略的因素. ２０１１年烟台市统计年鉴数据显示ꎬ２０１０年龙口市农药使用量已经达到２ ５４２ ｔꎬ化肥使用实物量更是达到１２􀆰 ２９１０４ ｔ. ＬＶ等[２４]发现ꎬ耕地土壤中ｗＣｄ与ｗＴＰ呈显著相关.Ｃｄ常作为磷矿石中的杂质而存在于磷肥当中ꎬ通过磷肥的施用在土壤中累积.王飞等[２５]也发现ꎬ华北地区禽畜粪便有机肥中ｗＰｂ、ｗＣｄ普遍存在严重的超标现象.可见ꎬ耕地土壤重金属的富集与肥料施用有直接关系.另外ꎬ农用塑料薄膜使用不当以及农用车的普及也是导致耕地重金属富集的原因ꎬ如龙口市在２０１０年的农用塑料薄膜使用量达到１ １１５ ｔꎬ农用柴油量为１􀆰 ２１１０４ ｔ.但由于农业活动产生的粉尘、固体废物及有害气体等有限ꎬ大气活动对其传输８２２１第７期李 彤等龙口市北部平原区不同土地利用类型对土壤重金属污染距离的影响 过程的影响甚微ꎬ且土壤重金属具有弱移动性ꎬ使得受污染的耕地、园地等农用地对土壤重金属的污染距离较小.研究区西部与北部的交通以及工业生产、电厂等煤炭燃烧产生的“三态”废弃物也是热点成因分析的一个不可忽略的因素.研究区内工厂多位于西北部ꎬ通过查阅２０１６年烟台市统计年鉴ꎬ发现２０１５年龙口市废水排放量达３ ５５１􀆰 ６４１０４ ｔꎬ其中ꎬ工业废水量为１ ２０２􀆰 ６４１０４ ｔꎬ工业烟粉尘排放量达０􀆰 ４３１０４ ｔꎬ工业ＳＯ２排放量为２􀆰 ３３１０４ ｔꎬＮＯｘ排放量为２􀆰 ００１０４ ｔ.废水通过灌溉进入土壤ꎬ而这些气态重金属进入大气后ꎬ在大气中传输ꎬ最终通过干湿沉降进入土壤ꎬ所以这些通过大气沉降进入土壤的重金属移动距离相对较大ꎬ且不可忽略.王幼奇等[２６]也认为ꎬ通过沉降进入土壤的重金属必需得到重视.研究区西部的龙口城区及中部的黄城城区拥有相对密集的交通网络ꎬ随着龙口市私人汽车保有量及农用车的持续增长２０１５年家用汽车、助力车是２０１０年的２８３􀆰 ５０％、２４６􀆰 ６４％ꎬ汽车轮胎磨损及油泄漏等交通排放产生的有害气体及粉尘、烟尘也不可忽视[２７￣２８]ꎻ尤其是上下班期间的高峰期堵车严重ꎬ发动机运行时不完全燃烧[２９]ꎬ尾气排放量增加ꎬ其中的重金属随之在空气中传输ꎬ扩散到道路周围ꎬ由于其迁移受到风速、风向等的影响ꎬ移动距离相对较大.研究区内位于工业区的采样点主要在东海、图新等多个大型热电厂附近ꎬ重金属含量明显偏高. ＬＶ等[３０￣３１]也发现ꎬ煤炭燃烧是工厂附近土壤中Ｐｂ、Ｃｄ的主要来源之一.因此ꎬ交通及工业燃炭过程是重金属的重要来源ꎬ且重金属在该过程中的迁移能力较强ꎬ污染范围较大.虽然该研究中的采样点都具有典型性ꎬ可以代表研究区不同土地利用下的重金属污染水平ꎬ但由于采样点有限ꎬ可能会对结论造成一定的干扰ꎬ但研究结果证实了不同土地利用类型对土壤重金属污染距离的差异性是实际存在的ꎬ并有效揭示了其污染的梯度区间.４ 结论ａ研究区土壤中ｗＰｂ、ｗＣｄ的平均值均超过山东省土壤相应背景值ꎬ其中ꎬｗＰｂ超过背景值的样品比例达到９６􀆰 ３６％ꎬ但均未超过ＧＢ １５６１８１９９５二级标准限值.龙口市北部平原区的西部存在明显的重金属集聚现象ꎬ东部污染并不严重.ｂ土地利用类型面积占比对Ｐｂ、Ｃｄ的污染距离均呈规律性变化.各土地利用类型在相关系数达到最大时对应的污染距离从大到小依次为交通用地建筑用地耕地≈园地.这种明显的距离差异要求在进行重金属污染治理时ꎬ对不同的土地利用类型需要根据其对重金属污染距离的差异加以区别ꎬ这也是该研究的实际意义所在.ｃ对研究区内Ｐｂ、Ｃｄ富集贡献最大的土地利用类型是耕地ꎬ其次是交通用地和建设用地.因此ꎬ在进行该区土壤治理时ꎬ要重点关注耕地土壤的修复ꎬ如建立污水处理厂解决污水灌溉问题、采用超富集植物修复等.参考文献Ｒｅｆｅｒｅｎｃｅｓ[ １ ] ＧＵ Ｙ ＧꎬＧＡＯ Ｙ ＰꎬＬＩＮ Ｑ. Ｃｏｎｔａｍｉｎａｔｉｏｎꎬｂｉｏａｃｃｅｓｓｉｂｉｌｉｔｙ ａｎｄｈｕｍａｎ ｈｅａｌｔｈ ｒｉｓｋ ｏｆ ｈｅａｖｙ ｍｅｔａｌｓ ｉｎ ｅｘｐｏｓｅｄ￣ｌａｗｎ ｓｏｉｌｓ ｆｒｏｍ ２８ｕｒｂａｎ ｐａｒｋｓ ｉｎ ｓｏｕｔｈｅｒｎ Ｃｈｉｎａ′ｓ ｌａｒｇｅｓｔ ｃｉｔｙꎬ Ｇｕａｎｇｚｈｏｕ [ Ｊ].Ａｐｐｌｉｅｄ Ｇｅｏｃｈｅｍｉｓｔｒｙꎬ２０１６ꎬ６７５２￣５８.[ ２ ] ＫＡＹＥ Ｊ ＰꎬＧＲＯＦＦＭＡＮ Ｐ ＭꎬＧＲＩＭＭ Ｎ Ｂꎬｅｔ ａｌ.Ａ ｄｉｓｔｉｎｃｔ ｕｒｂａｎｂｉｏｇｅｏｃｈｅｍｉｓｔｒｙ [Ｊ]. Ｔｒｅｎｄｓ ｉｎ Ｅｃｏｌｏｇｙ ＆ Ｅｖｏｌｕｔｉｏｎꎬ２００６ꎬ２１４１９２￣１９９.[ ３ ] ＮＡＮＯＳ Ｎꎬ ＭＡＲＴＩＮ Ｊ Ａ Ｒ. Ｍｕｌｔｉｓｃａｌｅ ａｎａｌｙｓｉｓ ｏｆ ｈｅａｖｙ ｍｅｔａｌｃｏｎｔｅｎｔｓ ｉｎ ｓｏｉｌｓ ｓｐａｔｉａｌ ｖａｒｉａｂｉｌｉｔｙ ｉｎ ｔｈｅ Ｄｕｅｒｏ Ｒｉｖｅｒ ＢａｓｉｎＳｐａｉｎ[Ｊ].Ｇｅｏｄｅｒｍａꎬ２０１２ꎬ１８９∕１９０６５５４￣５６２.[ ４ ] ＬＶ ＪｉａｎｓｈｕꎬＺＨＡＮＧ ＺｕｌｕꎬＬＩＵ Ｙａｎｇꎬｅｔ ａｌ.Ｓｏｕｒｃｅｓ ｉｄｅｎｔｉｆｉｃａｔｉｏｎａｎｄ ｈａｚａｒｄｏｕｓ ｒｉｓｋ ｄｅｌｉｎｅａｔｉｏｎ ｏｆ ｈｅａｖｙ ｍｅｔａｌｓ ｃｏｎｔａｍｉｎａｔｉｏｎ ｉｎＲｉｚｈａｏ Ｃｉｔｙ[Ｊ].Ａｃｔａ Ｇｅｏｇｒａｐｈｉｃａ Ｓｉｎｉｃａꎬ２０１２ꎬ６７７１０９￣１２２.[ ５ ] 张一修ꎬ王济ꎬ秦樊鑫ꎬ等.贵阳市道路灰尘和土壤重金属来源识别比较[Ｊ].环境科学学报ꎬ２０１２ꎬ３２１２０４￣２１２.ＺＨＡＮＧ ＹｉｘｕꎬＷＡＮＧ ＪｉꎬＱＩＮ Ｆａｎｘｉｎꎬｅｔ ａｌ.Ｃｏｍｐａｒｉｓｏｎ ｏｆ ｓｏｕｒｃｅｉｄｅｎｔｉｆｉｃａｔｉｏｎ ｏｆ ｒｏａｄ ｄｕｓｔ ａｎｄ ｓｏｉｌ ｈｅａｖｙ ｍｅｔａｌｓ ｉｎ Ｇｕｉｙａｎｇ[Ｊ].Ａｃｔａ Ｓｃｉｅｎｔｉａｅ Ｃｉｒｃｕｍｓｔａｎｔｉａｅꎬ２０１２ꎬ３２１２０４￣２１２.[ ６ ] ＺＨＡＮ Ｙ ＦꎬＳＨＩ Ｘ ＺꎬＨＵＡＮＧ Ｂꎬｅｔ ａｌ.Ｓｐａｔｉａｌ ｄｉｓｔｒｉｂｕｔｉｏｎ ｏｆ ｈｅａｖｙｍｅｔａｌｓ ｉｎ ａｇｒｉｃｕｌｔｕｒａｌ ｓｏｉｌｓ ｏｆ ａｎ ｉｎｄｕｓｔｒｙ￣ｂａｓｅｄ ｐｅｒｉ￣ｕｒｂａｎ ａｒｅａ ｉｎＷｕｘｉꎬＣｈｉｎａ[Ｊ].Ｐｅｄｏｓｐｈｅｒｅꎬ２００７ꎬ１７１４４￣５１.[ ７ ] ＹＩ Ｙ ＪꎬＹＡＮＧ ＺꎬＺＨＡＮＧ Ｓ.Ｅｃｏｌｏｇｉｃａｌ ｒｉｓｋ ａｓｓｅｓｓｍｅｎｔ ｏｆ ｈｅａｖｙｍｅｔａｌｓ ｉｎ ｓｅｄｉｍｅｎｔ ａｎｄ ｈｕｍａｎ ｈｅａｌｔｈ ｒｉｓｋ ａｓｓｅｓｓｍｅｎｔ ｏｆ ｈｅａｖｙｍｅｔａｌｓ ｉｎ ｆｉｓｈｅｓ ｉｎ ｔｈｅ ｍｉｄｄｌｅ ａｎｄ ｌｏｗｅｒ ｒｅａｃｈｅｓ ｏｆ ｔｈｅ ＹａｎｇｔｚｅＲｉｖｅｒ Ｂａｓｉｎ[Ｊ].Ｅｎｖｉｒｏｎｍｅｎｔａｌ Ｐｏｌｌｕｔｉｏｎꎬ２０１１ꎬ１５９１０２５７５￣２５８５.[ ８ ] 上官宇先ꎬ秦晓鹏ꎬ赵冬安ꎬ等.利用大型土柱自然淋溶条件下研究土壤重金属的迁移及形态转化[Ｊ].环境科学研究ꎬ２０１５ꎬ２８７１０１５￣１０２４.ＳＨＡＮＧＧＵＡＮ Ｙｕｘｉａｎꎬ ＱＩＮ Ｘｉａｏｐｅｎｇꎬ ＺＨＡＯ Ｄｏｎｇ′ａｎꎬ ｅｔ ａｌ.Ｍｉｇｒａｔｉｏｎ ａｎｄ ｔｒａｎｓｆｏｒｍａｔｉｏｎ ｏｆ ｈｅａｖｙ ｍｅｔａｌｓ ｉｎ ｓｏｉｌｓ ｂｙ ｌｙｓｉｍｅｔｅｒｓｔｕｄｙ ｗｉｔｈ ｆｉｅｌｄ ｃｏｎｄｉｔｉｏｎ[Ｊ].Ｒｅｓｅａｒｃｈ ｏｆ Ｅｎｖｉｒｏｎｍｅｎｔａｌ Ｓｃｉｅｎｃｅｓꎬ２０１５ꎬ２８７１０１５￣ １０２４.[ ９ ] ＨＵＡＮＧ Ｈ Ｊꎬ ＹＵＡＮ Ｘ Ｚ. Ｔｈｅ ｍｉｇｒａｔｉｏｎ ａｎｄ ｔｒａｎｓｆｏｒｍａｔｉｏｎｂｅｈａｖｉｏｒｓ ｏｆ ｈｅａｖｙ ｍｅｔａｌｓ ｄｕｒｉｎｇ ｔｈｅ ｈｙｄｒｏｔｈｅｒｍａｌ ｔｒｅａｔｍｅｎｔ ｏｆｓｅｗａｇｅ ｓｌｕｄｇｅ[Ｊ].Ｂｉｏｒｅｓｏｕｒｃｅ Ｔｅｃｈｎｏｌｇｙꎬ２０１６ꎬ２００９９１￣９９８.[１０] ＭＡＡＳ ＳꎬＳＣＨＥＩＦＬＥＲ ＲꎬＢＥＮＳＬＡＭＡ Ｍꎬｅｔ ａｌ.Ｓｐａｔｉａｌ ｄｉｓｔｒｉｂｕｔｉｏｎｏｆ ｈｅａｖｙ ｍｅｔａｌ ｃｏｎｃｅｎｔｒａｔｉｏｎｓ ｉｎ ｕｒｂａｎꎬｓｕｂｕｒｂａｎ ａｎｄ ａｇｒｉｃｕｌｔｕｒａｌｓｏｉｌｓ ｉｎ ａ Ｍｅｄｉｔｅｒｒａｎｅａｎ ｃｉｔｙ ｏｆ Ａｌｇｅｒｉａ [ Ｊ ]. ＥｎｖｉｒｏｎｍｅｎｔａｌＰｏｌｌｕｔｉｏｎꎬ２０１０ꎬ１５８６２２９４￣２３０１.[１１] ＸＩＥ Ｙ ＦꎬＣＨＥＮ Ｔ ＢꎬＬＥＩ Ｍꎬｅｔ ａｌ.Ｓｐａｔｉａｌ ｄｉｓｔｒｉｂｕｔｉｏｎ ｏｆ ｓｏｉｌ ｈｅａｖｙ９２２１环 境 科 学 研 究第３２卷ｍｅｔａｌ ｐｏｌｌｕｔｉｏｎ ｅｓｔｉｍａｔｅｄ ｂｙ ｄｉｆｆｅｒｅｎｔ ｉｎｔｅｒｐｏｌａｔｉｏｎ ｍｅｔｈｏｄｓａｃｃｕｒａｃｙ ａｎｄ ｕｎｃｅｒｔａｉｎｔｙ ａｎａｌｙｓｉｓ[Ｊ].Ｃｈｅｍｏｓｐｈｅｒｅꎬ２０１１ꎬ８２３４６８￣４７６.[１２] ＮＯＴＴＥＮ Ｍ ＪꎬＯＯＳＴＨＯＥＫ Ａ ＪꎬＲＯＺＥＭＡ Ｊꎬｅｔ ａｌ.Ｈｅａｖｙ ｍｅｔａｌｃｏｎｃｅｎｔｒａｔｉｏｎｓ ｉｎ ａ ｓｏｉｌ￣ｐｌａｎｔ￣ｓｎａｉｌ ｆｏｏｄ ｃｈａｉｎ ａｌｏｎｇ ａ ｔｅｒｒｅｓｔｒｉａｌ ｓｏｉｌｐｏｌｌｕｔｉｏｎ ｇｒａｄｉｅｎｔ [ Ｊ]. Ｅｎｖｉｒｏｎｍｅｎｔａｌ Ｐｏｌｌｕｔｉｏｎꎬ２００５ꎬ１３８ １１７８￣１９０.[１３] ＫＥＳＨＡＶＡＲＺＩ Ｂꎬ ＴＡＺＡＲＶＩ Ｚꎬ ＲＡＪＡＢＺＡＤＥＨ Ｍ Ａꎬ ｅｔ ａｌ.Ｃｈｅｍｉｃａｌ ｓｐｅｃｉａｔｉｏｎꎬｈｕｍａｎ ｈｅａｌｔｈ ｒｉｓｋ ａｓｓｅｓｓｍｅｎｔ ａｎｄ ｐｏｌｌｕｔｉｏｎｌｅｖｅｌ ｏｆ ｓｅｌｅｃｔｅｄ ｈｅａｖｙ ｍｅｔａｌｓ ｉｎ ｕｒｂａｎ ｓｔｒｅｅｔ ｄｕｓｔ ｏｆ ＳｈｉｒａｚꎬＩｒａｎ[Ｊ].Ａｔｍｏｓｐｈｅｒｉｃ Ｅｎｖｉｒｏｎｍｅｎｔꎬ２０１５ꎬ１１９１￣１０.[１４] ＭＡＲＴＩＮ Ｊ Ａ ＲꎬＲＡＭＯＳ￣ＭＩＲＡＳ Ｊ ＪꎬＢＯＬＵＤＡ Ｒꎬｅｔ ａｌ.Ｓｐａｔｉａｌｒｅｌａｔｉｏｎｓ ｏｆ ｈｅａｖｙ ｍｅｔａｌｓ ｉｎ ａｒａｂｌｅ ａｎｄ ｇｒｅｅｎｈｏｕｓｅ ｓｏｉｌｓ ｏｆ ａＭｅｄｉｔｅｒｒａｎｅａｎ ｅｎｖｉｒｏｎｍｅｎｔ ｒｅｇｉｏｎ Ｓｐａｉｎ [Ｊ].Ｇｅｏｄｅｒｍａꎬ２０１３ꎬ２００∕２０１６１８０￣１８８.[１５] 王祖伟ꎬ王袆玮ꎬ侯迎迎ꎬ等.于桥水库水源地水体沉积物重金属空间分异与景观格局的关系[Ｊ].环境科学ꎬ２０１６ꎬ３７９３４２３￣３４２９.ＷＡＮＧ ＺｕｗｅｉꎬＷＡＮＧ ＹｉｗｅｉꎬＨＯＵ Ｙｉｎｇｙｉｎｇꎬｅｔ ａｌ. Ｒｅｌａｔｉｏｎｓｈｉｐｂｅｔｗｅｅｎ ｌａｎｄｓｃａｐｅ ｐａｔｔｅｒｎ ａｎｄ ｓｐａｔｉａｌ ｖａｒｉａｔｉｏｎ ｏｆ ｈｅａｖｙ ｍｅｔａｌｓ ｉｎａｑｕａｔｉｃ ｓｅｄｉｍｅｎｔｓ ｉｎ ｈｅａｄｗａｔｅｒｓ ａｒｅａ ｏｆ Ｙｕｑｉａｏ Ｒｅｓｅｒｖｏｉｒ [ Ｊ].Ｅｎｖｉｒｏｎｍｅｎｔａｌ Ｓｃｉｅｎｃｅꎬ２０１６ꎬ３７９３４２３￣３４２９.[１６] 李铖ꎬ李芳柏ꎬ吴志峰ꎬ等.景观格局对农业表层土壤重金属污染的影响[Ｊ].应用生态学报ꎬ２０１５ꎬ２６４１１３７￣１１４４.ＬＩ Ｃｈｅｎｇꎬ ＬＩ Ｆａｎｇｂａｉꎬ ＷＵ Ｚｈｉｆｅｎｇꎬ ｅｔ ａｌ. Ｉｍｐａｃｔｓ ｏｆ ｌａｎｄｓｃａｐｅｐａｔｔｅｒｎｓ ｏｎ ｈｅａｖｙ ｍｅｔａｌ ｃｏｎｔａｍｉｎａｔｉｏｎ ｏｆ ａｇ