某市农业土壤重金属的空间分布特征及来源解析doctypepdf.docx

某市农业土壤重金属的空间分布特征及来源解析doctypepdf.docx

《某市农业土壤重金属的空间分布特征及来源解析doctypepdf.docx》由会员分享，可在线阅读，更多相关《某市农业土壤重金属的空间分布特征及来源解析doctypepdf.docx（11页珍藏版）》请在冰豆网上搜索。

某市农业土壤重金属的空间分布特征及来源解析doctypepdf

第２９卷第１２期　２００８年１２月　环　境　科　学　ＥＮＶＩＲＯＮＭＥＮＴＡＬ　ＳＣＩＥＮＣＥ　Ｖｏ１．２９．Ｎｏ．１２　

Ｄｅｃ．，２００８　东莞市农业土壤重金属的空间分布特征及来源解析　蔡立梅　，马瑾　，周永章　，黄兰椿　，窦磊　，张澄博　，付善明　（１．中山大学地球科学系，广州５１０２７５；２．中山大学地球环境与地球资源研究中心，广州５１０２７５；３．广东省地质调查院，　广州　５１０８００）　摘要：

对东莞市农业土壤进行了系统采样分析，测定了表层土壤中ｃｕ、ｚｎ、Ｎｉ、ｃｒ、Ｐｂ、ｃｄ、Ａｓ和Ｈｇ的含量．利用ＧＩＳ技术和多元　

地统计结合的方法，研究了上述８种重金属元素的空间分布特征和来源．结果表明，ｃｕ、ｚｎ、　、Ｐｂ、ｃｄ和Ｈｇ的含量均超过广东　省土壤背景值，尤其是Ｈｇ，Ｃｄ和Ｐｂ含量分别达到０．２４、０．１２和６５．３８　ｍｇ／ｋｇ，远高于背景值，说明随着东莞市快速工业化，土壤　中重金属呈明显积聚的趋势．上述８种元素可以由３个主因子反映．因子１为“自然源因子”，主要包括元素ｃｕ、ｚｎ、Ｎｉ、ｃｒ和　Ａｓ，其高值区主要分布在研究区西部，主要受河流冲积物海相沉积物控制，低值区的空间分布受酸性火山喷出岩风化物和花岗　

岩风化物影响，因子１的空间分布主要受成土母质控制．因子２为“工业及交通源因子”，主要包括元素Ｐｂ和Ｈｇ，高值区主要　

分布在中部过渡区和西部平原区，高值区中心集中在中部过渡区的虎门镇和西部平原区的茶山镇，主要由该区工业污染所　

致．因子３为“工农业活动因子”，主要包括ｃｄ元素，其空间变异完全受人为随机因素控制．　

关键词：

农业土壤；重金属；多元地统计；空间分布；来源　中图分类号：

Ｘ５３文献标识码：

Ａ文章编号：

０２５０．３３０１（２００８）１２—３４９６．Ｏ７　Ｍｕｌｔｉｖａｒｉａｔｅ　Ｇｅｏｓｔａｔｉｓｔｉｃｓ　ａｎｄ　ＧＩＳ?

ｂａｓｅｄ　Ａｐｐｒｏａｃｈ　ｔｏ　Ｓｔｕｄｙ　ｔｈｅ　Ｓｐａｔｉａｌ　Ｄｉｓｔｒｉｂｕｔｉｏｎ　ａｎｄ　Ｓｏｕｒｃｅｓ　ｏｆ　Ｈｅａｖｙ　Ｍｅｔａｌｓ　ｉｎ　Ａｇｒｉｃｕｌｔｕｒａｌ　Ｓｏｉｌ　ｉｎ　ｔｈｅ　Ｐｅａｒｌ　Ｒｉｖｅｒ　Ｄｅｌｔａ，Ｃｈｉｎａ　ＣＡＩ　Ｌｉ—ｍｅｉ　，　＇　，ＭＡ　Ｊｉｎ　，　，ＺＨＯＵ　Ｙｏｎｇ．ｚｈａｎｇ　＇　，ＨＵＡＮＧ　Ｌａｎ．ｃｈｕｎ　，　，Ｄ０Ｕ　Ｌｅｉ　，　，ＺＨＡＮＧ　Ｃｈｅｎｇ－ｂｏ　’　，ＦＵ　Ｓｈａｎ一　?

　１，２　ｎｎｎｇ　（１．Ｄｅｐａｒｔｍｅｎｔ　ｏｆ　Ｅａｒｔｈ　Ｓｃｉｅｎｃｅ，Ｓｕｎ　Ｙａｔ—ｓｅｎ　Ｕｎｉｖｅｒｓｉｔｙ，Ｇｕａｎｇｚｈｏｕ　５１０２７５，Ｃｈｉｎａ；２．Ｃｅｎｔｅｒ　ｆｏｒ　Ｅａｒｔｈ　Ｅｎｖｉｒｏｎｍｅｎｔ　ａｎｄ　Ｒｅｓｏｕｒｃｅｓ，Ｓｕｎ　Ｙａｔ—ｓｅｎ　Ｕｎｉｖｅｒｓｉｔｙ，Ｇｕａｎｇｚｈｏｕ　５１０２７５，Ｃｈｉｎａ；３．Ｇｅｏｌｏｇｉｃａｌ　Ｉｎｖｅｓｔｉｇａｔｉｏｎ　Ｉｎｓｔｉｔｕｔｅ　ｏｆ　Ｇｕａｎｇｄｏｎｇ　Ｐｒｏｖｉｎｃｅ，Ｇｕａｎｇｚｈｏｕ　５１０８００，Ｃｈｉｎａ）　Ａｂｓｔｒａｃｔ：

Ｏｎｅ　ｈｕｎｄｒｅｄ　ａｎｄ　ｅｉｇｈｔｅｅｎ　ｓｕｒｆａｃｅ　ｓｏｉｌ　ｓａｍｐｌｅｓ　ｗｅｒｅ　ｃｏｌｌｅｃｔｅｄ　ｆｒｏｍ　ｔｈｅ　Ｄｏｎｇｇｕａｎ　Ｃｉｔ?

　　0１２期　蔡立梅等：

东莞市农业土壤重金属的空间分布特征及来源解析　万　．在快速工业化过程的同时，也给当地环境造　

成一定污染压力，尤其是重金属污染．本研究以多元　统计和地统计方法相结合为手段，对东莞市农地土　壤重金属的空间分布特征和来源进行了系统地分　

析，旨在为该区农业土壤污染综合防治提供参考，并　为快速工业化、城市化地区土壤资源的持续利用和　合理管理提供科学依据．　１材料与方法　１．１研究区概况　

东莞市地处珠江三角洲东部，位于广、深两大都　市之间，其地貌由网河平原和低山丘陵构成，２种截　

然不同的地貌和地理条件使东莞市分异演化出３类　

明显不同的地域（图１）．西北部为东江三角洲河网　

平原区，水资源丰富，主要成土母质为东江河流沉积　

物．其中位于三角洲上游东段的石排镇、横沥镇、企　石镇和桥头镇作为东深供水工程和东莞运河水源　圆　圈　Ｖｚ－］　东南丘陵区　镇　采样点　图１研究区分区及采样点分布　Ｆｉｇ．１　Ｚｏｎｅｓ　ａｎｄ　ｓａｍｐｌｅ　ｓｉｔｅｓ　ｏｆ　ｓｔｕｄｉｅｄ　ａｒｅａ　区，自１９９１年以来环境受到严格保护，此区称为水　乡上游保护区（工区），总面积２１３．０　ｋｍ２，占全市面　

积的８．６％．其它１１个镇总面积５４９．７　ｋｍｚ，占全市　总面积的２２．２％，该地域过境水资源丰富，乡镇企　业密集，可称之为西部平原区（Ⅱ区）．东南部为低山　丘陵地区（Ⅳ区），水资源短缺，包括９个镇，总面积　１　０１８．６　ｋｍ２，占全市面积的４１．２％，成土母质主要为　花岗岩、砂页岩；中部为山区与水乡的过渡地带，称　为中部过渡地带（Ⅲ区），包括常平镇、厚街镇和虎门　镇等８个镇（区），面积６９３．５　ｋｍ２，占全市面积的　２８．０％，成土母质为三角洲近代沉积物和花岗岩、红　色砂页岩　．　１．２样品采集　

根据工业布局、“三废”排放状况、灌溉水类型以　及农地面积进行布点（图１），共采集农地表层土壤　

样（０～２０　ｃｍ）１１８个，表土采用多点采样方式（１０～　

１５个点），等量混合后取１　ｋｇ混合样．　１．３样品处理与测试　土壤样品自然风干，剔除样品中的植物根系、有　

机残渣以及可见侵入体，用木棍碾碎并用玛瑙研钵　

研磨，过１００目尼龙筛．所有处理样品经标准方法进　行消化处理后，Ｈｇ总量采用冷原子吸收法测定，Ａｓ　

总量采用ＤＤＣＡｇ分光光度法测定，其它元素采用原　子吸收分光光度法测定．样品测试在广东省生态环　

境与土壤研究所分析测试中心完成，分析过程中所　

用的试剂均为优级纯．土壤样品分析过程中加入国　家标准物质样品ＧＳＳ一１作为未知样品的测定以进行　

分析质量控制，同时进行空白和试剂的实验¨　．　土壤样品基本性质的测定采用常规方法：

ｐＨ采　用电位法（水：

土＝２．５：

１），有机质含量采用重铬酸　钾外加热法测定，测试结果见表１．可见土壤偏酸　

性，有机质含量处于中等偏下水平．　表１东莞市农业土壤基本性质　Ｔａｂｌｅ　Ｉ　Ｂａｓｉｃ　ｐｒｏｐｅｒｔｉｅｓ　ｏｆ　ａ　ｃｕｌｔｕｒａｌ　ｓｏｉｌｓ　ｉｎ　Ｄｏｎｇｇｕａｎ　Ｃｉｔｙ　４　…，…　，　。

　２结果与讨论　

常规数理统计分析采用ＳＰＳＳ１３．０软件完成，地　………　统计分析采用Ｒ软件拟合半变异函数、建立拟合模　２．１土壤重金属的统计性描述　型．在ＡｒｃＧＩＳ９．１中?

　　环　境　科　学　２９卷　分另０为２１．８２、６６．１５、２０．５２、４３．０１、６５．３８、０．１２、　１２．７６和０．２４　ｍｇ／ｋｇ，均没有超过国家二级标准中的　

土壤环境质量限定值ｎ　（表３），但有部分样品的　

ｃｕ、Ｎｉ、ｃｄ和Ｈｇ含量超过了限定值，尤其是Ｈｇ含量　最高值已达到１．０１　ｍｇ／ｋｇ，远高于限定值０．３　ｍｇ／ｋｇ，　３０％土壤样品Ｈｇ含量超过了二级标准，应引起重　视．除ｃｒ和Ａｓ外，其余元素总量平均值都超过广东　省土壤背景值¨　，尤其是　、ｃｄ、Ｐｂ平均含量分别　达到０．２４、０．１２和６５．３８　ｍｇ／ｋｇ，远高于背景值，说　明随着东莞市快速工业化，土壤中重金属呈明显积　

聚的趋势．从变异系数看，ｃｄ和｝ｉｇ的变异系数相对　

较大，其余元素变异程度中等．　表２东莞市农业土壤重金属元素的统计性描述　Ｔａｂｌｅ　２　Ｄｅｓｃｒｉｐｔｉｖｅ　ｓｔａｔｉｓｔｉｃｓ　ｏｆ　ｔｈｅ　ｈｅａｖｙ　ｍｅｔａｌｓ　ｏｆ　ａｇｒｉｃｕｌｔｕｒａｌ　ｓｏｉｌｓ　ｉｎ　Ｄｏｎｇｇｕａｎ　Ｃｉｔｙ　表３土壤重金属评价标准　Ｔａｂｌｅ　３　Ｅｖａｌｕａｔｅｄ　ｓｔａｎｄａｒｄｓ　ｏｆ　ｈｅａｖｙ　ｍｅｔａｌ　ｃｏｎｃｅｎｔｒａｔｉｏｎｓ　ｉｎ　ｓｏｉｌ　ｅｎｖｉｒｏｎｍｅｎｔ　２．２土壤重金属的相关分析　

表４计算了土壤ｐＨ、有机质和土壤重金属元素　问的相关系数，从表４可以看出，ｃｕ、ｚｎ、Ｎｉ、Ｃｒ、Ｐｂ、　Ａｓ、Ｈｇ与有机质有较好的相关性，表明这些元素可　能来源于成土母质．Ｃｕ与ｚｎ、Ｎｉ、Ｃｒ、Ａｓ，ｚｎ与Ｎｉ、　ｃｒ、Ａｓ，Ｎｉ与ｃｒ、Ｐｂ、Ａｓ，ｃｒ与Ｐｂ，Ｐｂ与Ｈｇ之间存在　显著正相关．由此可初步推断，ｃｕ、ｚｎ、Ｎｉ、Ｃｒ和Ａｓ　

等５种元素可能同一来源，元素Ｐｂ与Ｈｇ同一来源，　

上述７种元素的来源与成土母质可能有关，而Ｃｄ元　素为其它来源．　表４东莞市农业土壤土壤ｐＨ、有机质和土壤重金属元素间相关系数　Ｔａｂｌｅ　４　Ｃｏｒｒｅｌａｔｉｏｎｓ　ｍａｔｒｉｘ　ｆｏｒ　ｔｈｅ　ｈｅａｖｙ　ｍｅｔａｌｓ，ｐＨ　ａｎｄ　ＳＯＭ　ｉｎ　ｔｈｅ　ａｇｒｉｃｕｌｔｕｒａｌ　ｓｏｉｌｓ　

ｐＨ　ＳＯＭ　Ｃｕ　Ｚｎ　Ｎｉ　ｐＨ　ＳＯＭ　Ｃｕ　

Ｚｎ　

Ｎｉ　

Ｃｒ　Ｐｂ　

Ｃｄ　

Ａｓ　Ｈｇ　１　

０．３８４　１　

０．５２１　０．４７７　

０．４６４　０．４３８　０．５２０　０．４８３　

Ｏ．４４３　０．１８６　

Ｏ．Ｏ９０　０．１２５　０．３９９　０．５２７　０．５ｏｏ　０．０４９　１　０．１２７　ｌ　

０．２０４　０．Ｏ６２　Ｉ）＊＊表示９９％置信区间显著，＊表示９５％置信区间显著　２．３土壤重金属的因子分析　

因子分析方法可用来判别土壤中重金属来　

源　”　．对８种元素进行因子分析，提取的３个因子解　

释了总方差的７５％，表５指示了８种元素所属类别．　

因子１主要包括ｃｕ、ｚｎ、Ｎｉ、Ｃｒ和Ａｓ等元素（表６），对　照讨论结果，这些元素含量与对应元素背景值相当，　１　Ｏ　０　Ｏ　

１　Ｏ　Ｏ　Ｏ　０　蜘唧　ｌ　０　Ｏ　０　Ｏ　Ｏ　跏　ｌ　Ｏ　Ｏ　Ｏ　Ｏ　Ｏ　Ｏ　瑚　栅　㈣　Ｏ　Ｏ　Ｏ　Ｏ　Ｏ　Ｏ　Ｏ　

一　一　一　一　一　一　。

螅幔恚穑欤澹蟆。

鳎澹颍濉?

冒_瑞_１２期　蔡立梅等：

东莞市农业土壤重金属的空间分布特征及来源解析　且这组元素与土壤有机质显著相关（表４），一般认为　Ｃｒ和Ｍ主要受地球化学成因影响…，　．Ｂ０ｍｖｋａ　在捷克东北部的研究也认为元素ｃｕ、ｚｎ、Ｎｉ和Ｃｒ主　要为地质来源．故将因子１作为“自然源因子”．　表５因子特征值　Ｔａｂｌｅ　５　Ｅｉｇｅｎｖａｌｕｅｓ　ｏｆ　ｆａｃｔｏｒｓ　元素Ｐｂ和№在因子２中占有较大的因子载荷　（表６）．从Ｐｂ和Ｈｇ在各区的分布可知，西部平原　区、中部过渡区和东南丘陵区Ｐｂ污染都较普遍，尤　以乡镇企业密集的西部平原区和中部过渡区最为严　

重，这说明土壤№、Ｐｂ的来源与该地区的工业有密　

切关系．在城市土壤中Ｐｂ常被作为机动车污染源的　

标识元素　，土壤Ｐｂ的来源与东莞市的公路网密　

集以及机动车辆的急剧增加有关．因子２可称为“工　业及交通源因子”．　因子３主要包括ｃｄ元素，ｃｄ一般可作为施用　

农药和化肥等农业活动的标识元素　。

　．据统计，自　１９９０年以来东莞市共施用化肥超过２００万ｔ（实物　量）、农药３．５万ｔｎ　，其施用量远高于全国平均用　量，我国农药和化肥的利用率普遍较低，约有７０％　流失到土壤、水和空气中　，长期的农业生产活动，　使ｃｄ在土壤中明显累积．此外，密集的工业活动也　

是ｃｄ的重要来源之一，刘庆等　研究浙江省慈溪　

市农田土壤重金属污染时，亦发现ｃｄ污染为工业　

“三废”排放．故因子３可称作“工农业活动因子”．　２．４因子的结构特征和空间分布特征　２．４．１因子的结构分析　半方差函数在原点处的数值称为块金常数　（Ｃ　），它由测定误差和小于最小采样尺度的非连续　性变异引起，属于随机变异；基台值（ｃ。

＋Ｃ　）通常　表示系统内的总变异，它是结构性变异和随机性变　

异之和．块金值Ｃ。

与基台值（Ｃ。

＋Ｃ　）之比是反映　

区域化变量空间异质性程度的重要指标，称为块金　

效应．该比值用于反映空间变异影响因素中区域因　

素（自然因素）和非区域因素（人为因素）的作用．当　［Ｃ。

／（Ｃ。

＋Ｃ　）］＜２５％时，表明变量的空间变异以　

结构性变异为主，变量具有强烈的空间相关性；当上　述比值２５％＜［ｃ０／（Ｇ０＋Ｃ　）］＜７５％时，变量为中　

等程度空间相关；［ｃ。

／（Ｃ。

＋Ｃ　）］＞７５％时，以随机　变异为主，变量的空间相关性很弱　．　利用Ｒ软件拟合因子的半方差函数理论模型　

及其相关参数见表７．根据表７，因子１和因子２的　

变异函数理论模型分别符合指数模型和高斯模型，　块金值／基台值分别为３３％和６７％，可见，因子１的　€筙x4_G_!

皩`4_

氢1皥v

冒4栶_３５００　环　境　科　学　２９卷　空间变异以土壤母质、地形、土壤类型等结构性变异　

为主，变量具有强烈的空间相关性；而因子２的空间　变异以随机性变异为主，空间变异主要受人为因素　

的影响，变量具有中等的空间相关性．从图２（ｃ）可　

看出，因子３的空间变异完全受人为随机因素的控　

制，在取样尺度内不具有空间相关性．因子的结构分　析也初步验证了因子分类的正确性．　０．８　０　６　凄ｏ．４　｛｝Ｌ　０．２　

ｌ?

Ｏ　

０．８　凄ｏ．６　。

一４　０　２　糕１　Ｏ　闰　酞　ｏ．５　图２因子的关变异函数　Ｆｉｇ．２　Ｓｅｍｉｖａｒｉｏｇｒａｍ　ｍａｐ　ｏｆ　ｆａｃｔｏｒｓ　２．４．２因子的空间分布特征及成因分析　对变量连续插值可更直观地查看其空间分布特　征，图３、４分别是利用ＡｒｃＧＩＳ地统计模块进行普通　Ｋｒｉｇｉｎｇ插值后，因子的插值结果．　

从图３可看出，因子１的高值区主要分布在研　

究区的西部，区域上主要属于西部平原区和中部过　

渡区．对照研究区的成土母质图（图５）可看出　，因　图３因子１的空间分布　Ｆｉｇ．３　Ｓｐａｃｉａｌ　ｄｉｓｔｒｉｂｕｔｉｏｎｓ　ｏｆ　ｆａｃｔｏｒ　１　［］＿２￣一０．６　

圈…０．６　０　５　

圜＿０＿５—０．２　一０．２～０　２　目０．２～０．７　

圈ｏ＿７￣５．０　图４因子２的空间分布　Ｆｉｇ．４　Ｓｐａｃｉａｌ　ｄｉｓｔｒｉｂｕｔｉｏｎｓ　ｏｆ　ｆａｃｔｏｒ　２　子１高值区与河流冲击物海相沉积物的空间分布十　分一致．而低值区的空问分布大致受酸性火山喷出　岩风化物和花岗岩风化物所控制．这说明因子１的　分布主要受成土母质的控制，与前面结构分析的结　论一致．　从图４可看出，因子２的高值区主要分布在中　

部过渡区和西部平原区，高值区中心集中在中部过　圈　¨　囹目圈　，　＃一　酬口园圆　。

┲仁欠从场?

区睹皦Y?

１２期　蔡立梅等：

东莞市农业土壤重金属的空间分布特征及来源解析　３５０１　目河流冲积物海相沉积物　口花岗岩风化物　曰变质岩风化物　目砂页岩风化物　口酸性火山喷出岩风化物　图５东莞市土壤母质　Ｆｉｇ．５　Ｓｏｉｌ　ｐａｒｅｎｔ　ｍａｔｅｒｉａｌｓ　ｉｎ　Ｄｏｎｇｇｕａｎ　Ｃｉｔｙ　渡区的虎门镇和西部平原区的茶山镇．主要是由于　

该区工业污染所致，西部平原区由于水源充足、交通　方便，一些资源消耗大、污染程度大的制造业在改革　开放初期开始进入并迅速发展，随着这些重污染产　业的不断发展，大量重金属污染物排放到环境中，从　而也成为东莞市重金属污染最严重的区域．据统　计¨　，２００４年底虎门镇共有各种类型企业１４０多　家，其中有５０多家电镀厂和３家电厂，茶山镇也分　布着东莞市２５家重点污染企业，且多是造纸厂和电　镀厂．这说明因子２的来源与该地区的工业有密切　

关系．　

土壤中重金属的空间分布特征受到土壤质地、　地形、大气沉降和人类活动等多种影响．王学军　

等　研究了多种尺度的土壤微量金属元素含量的　空间结构，从小尺度的北京东郊污灌区，中尺度的深　

圳地区到大尺度的内蒙古．其研究认为在中尺度和　

大尺度的情况下，成土母质成为影响土壤中微量金　属元素空间分布的首要因素．本研究认为在乡镇企　

业密集的东莞市，成土母质和高强度的人类活动共　同控制土壤中重金属元素的空间分布，其中ｃｕ、ｚｎ、　Ｎｉ、ｃｒ和Ａｓ等元素的空间分布主要受成土母质的控　制，而Ｐｈ、№和ｃｄ元素则主要受人类活动影响．　

３结论　（１）东莞市农业土壤中Ｃｕ、ｚｎ、Ｎｉ、Ｃｒ、Ｐｈ、ｃｄ、Ａｓ　和Ｈｇ的平均含量均低于国家二级标准．但有部分　

样品的ｃｕ、Ｎｉ、ｃｄ和Ｈｇ含量超过了二级标准，尤其　是Ｈｇ含量最高值远高于限定值，３０％土壤样品Ｈｇ　含量超过了二级标准，应引起重视．除ｃｒ和Ａｓ外，　

其余元素含量都超过广东省土壤背景值，尤其是　

№、ｃｄ和Ｐｂ含量远高于背景值，说明随着东莞市快　速工业化，土壤中重金属呈明显积聚的趋势．　

（２）东莞市农业土壤中８种重金属元素有着不　同的来源．相关分析、因子分析和空间分析表明８种　重金属元素可分成３类．因子１为“自然源因子”，主　要包括ｃｕ、ｚｎ、Ｎｉ、ｃｒ和Ａｓ等元素．因子２为“工业　

及交通源因子”，主要包括元素Ｐｂ和Ｈｇ．因子３为　“工农业活动因子”，主要包括Ｃｄ元素．　

（３）结构分析表明，因子１的空间变异以结构性　变异为主，空间变异影响因素主要为自然因素，具有　强烈的空间相关性．因子２的空问变异以随机性变　异为主，空间变异主要受人为因素的影响，具有中等　

的空间相关性．因子３的空间变异完全受人为随机　

因素的控制，在取样尺度内不具有空间相关性．　

（４）空间分析表明，因子１的高值区主要分布在　研究区的西部，区域上主要属于西部平原区和中部　过渡区，主要受河流冲击物海相沉积物控制．低值区　的空间分布大致受酸性火山喷出岩风化物和花岗岩　风化物控制．因子ｌ的空问分布主要受成土母质控　

制．因子２的高值区主要分布在中部过渡　和　部　

平原区，高值区中心集中在中部过渡区的虎¨镇和　

西部平原区的茶山镇，主要由于该区工业污染所致．　参考文献：

　［１］Ｌｅｅ　Ｃ　Ｓ，Ｌｉ　Ｘ，Ｓｈｉ　Ｗ，ｅｔ　ａ１．Ｍｅｔａｌ　ｃｏｎｔａｍｉ?

　　）3_

_-N?

３５０２　环　境　科　学　２９卷　ＭｃＧｒａｔｈ　Ｄ，Ｚｈａｎｇ　Ｃ，Ｃａｒｔｉｎ　０　Ｔ．Ｇｅｏｓｔａｔｉｓｔｉｃａｌ　ａｎａｌｙｓｅｓ　ａｎｄ　ｈａｓａｒｄ　ａｓｓｅｓｓｍｅｎｔ　ｏｎ　ｓｏｉｌ　ｌｅａｄ　ｉｎ　Ｓｉｌｖｅｒｍｉｎｅｓ　ａｒｅａ，Ｉｒｅｌａｎｄ［Ｊ］．　Ｅｎｖｉｒｏｎｍｅｎｔａｌ　Ｐｏｌｌｕｔｉｏｎ，２００４，１２７：

２３９—２４８．　Ｓａｂｙ　Ｎ，Ａｒｒｏｕａｙｓ　Ｄ，Ｂｏｕｌｏｎｎｅ　Ｌ，ｅｔ　ａ１　Ｇｅｏｓｔａｔｉｓｔｉｅａｌ　ａｓｓｅｓｓｍｅｕｔ　ｏｆ　Ｐｂ　ｉｎ　ｓｏｉｌ　ａｒｏｕｎｄ　Ｐａｒｉｓ，Ｆｒａｎｃｅ［Ｊ］．Ｓｃｉｅｎｃｅ　ｏｆ　ｔｈｅ　Ｔｏｔａｌ　Ｅｎｖｉｒｏｎｍｅｎｔ，２００６，３６７：

２１２－２２１．　Ｌｉｕ　Ｘ，Ｗｕ　Ｊ，Ｘｕ　Ｊ．Ｃｈａｒａｃｔｅｒｉｚｉｎｇ　ｔｈｅ　ｒｉｓｋ　ａｓｓｅｓｓｍｅｎｔ　ｏｆ　ｈｅａｖｙ　ｍｅｔａｌｓ　ａｎｄ　ｓａｍｐｌｉｎｇ　ｕｎｃｅｒｔａｉｎｔｙ　ａｎａｌｙｓｉｓ　ｉｎ　ｐａｄｄｙ　ｆｉｅｌｄ　ｂｙ　ｇｅｏｓｔａｔｉｓｔｉｃｓ　ａｎｄ　ＧＩＳ［Ｊ］．Ｅｎｖｉｒｏｎｍｅｎｔａｌ　Ｐｏｌｌｕｔｉｏｎ，２００６，１４１：

２５７—　２６４．　东莞市统计局．东莞市统计年鉴（２００５）［Ｍ］．北京：

中国统计　出版社，２００５．６８，８９．１１５．　陆发嘉．珠江三角洲土壤［Ｍ］．北京：

中国环境科学出版社，　１９８８．６４—６８．　刘凤枝．农业环境监测使用手册［Ｍ］．北京：

中国标准出版　社．２００１．　ＧＢ　１５６１８—１９９５．土壤环境质量标准［Ｓ］．　广东省环境监测中心站．广东省土壤环境背景值数据集［Ｒ］，　１９９０．　Ｍａｒｔｉｎ　Ｊ　Ａ　Ｒ，Ａｒｉａｓ　Ｍ　Ｌ，Ｃｏｒｂｉ　Ｊ　Ｍ　Ｇ．Ｈｅａｖｙ　ｍｅｔａｌｓ　ｃｏｎｔｅｎｔｓ　ｉｎ　ａｇｒｉｃｕｌｔｕｒａｌ　ｔｏｐｓｏｉｌｓ　ｉｎ　ｔｈｅ　Ｅｂｒｏ　ｂａｓｉｎ（Ｓｐａｉｎ）．Ａｐｐｌｉｃａｔｉｏｎ　ｏｆ　ｔｈｅ　ｍｕｌｔｉｖａｒｉａｔｅ　ｇｅｏｓｔａｔｉｓｔｉｅａｌ　ｍｅｔｈｏｄｓ　ｔｏ　ｓｔｕｄｙ　ｓｐａｔｉａｌ　ｖａｒｉａｔｉｏｎｓ［Ｊ］．　Ｅｎｖｉｒｏｎｍｅｎｔａｌ　Ｐｏｌｌｕｔｉｏｎ，２００６，１４４（３）：

１００１．１０２１．　Ｈａｎｅｓｃｈ　Ｍ，Ｓｃｈｏｌｇｅｒ　Ｒ，Ｄｅｋｋｅｒｓ　Ｍ　Ｊ．Ｔｈｅ　ａｐｐｌｉｃａｔｉｏｎ　ｏｆ　Ｆｕｚｚｙ　ｃｌｕｓｔｅｒ　ａｎａｌｙｓｉｓ　ａｎｄ　ｎｏｎ－ｌｉｎｅａｒ　ｍａｐｐｉｎｇ　ｔｏ　ａ　ｓｏｉｌ　ｄａｔａ　ｓｅｔ　ｆｏｒ　ｄｅｔｅｃｔｉｏｎ　ｏｆ　ｐｏｌｌｕｔｅｄ　ｓｉｔｅｓ［Ｊ］．Ｐｈｙｓｉｃａｌ　Ｃｈｅｍｉｃａｌ　Ｅａｒｔｈ　Ｓｃｉｅｎｃｅ，２００１，２６：

　８８５—８９】．　Ｂｏｒｕｖｋａ　Ｌ，Ｖａｃｅｋ　Ｏ，Ｊｅｈｌｉｃｋａ　Ｊ．Ｐｒｉｎｃｉｐａｌ　ｃｏｍｐｏｎｅｎｔ　ａｎａｌｙｓｉｓ　ａｓ　ａ　ｔｏｏｌ　ｔｏ　ｉｎｄｉｃａｔｅ　ｔｈｅ　ｏｒｉｇｉｎ　ｏｆ　ｐｏｔｅｎｔｉａｌｌｙ　ｔｏｘｉｃ　ｅｌｅｍｅｎｔｓ　ｉｎ　ｓｏｉｌｓ［Ｊ］．　Ｇｅｏｄｅｒｍａ，２００５，１２８：

２８９—３００．　

Ｇｒａｙ　Ｃ　Ｗ，ＭｃＬａｒｅｎ　Ｒ　Ｇ，Ｒｏｂｅｒｔｓ　Ａ　Ｈ　Ｃ，ｅｔ　ａ１．Ｔｈｅ　ｅｆｆｅｃｔ　ｏｆｌ　　?

冒_橡_