廣東省揭陽市土壤鉻的空間結(jié)構(gòu)及分布特征研究①

王涵植，溫漢輝，蔡立梅, 4*，羅 杰，穆桂珍，王秋爽，蔣慧豪，王 碩

廣東省揭陽市土壤鉻的空間結(jié)構(gòu)及分布特征研究①

王涵植1, 2，溫漢輝3，蔡立梅1, 2, 4*，羅 杰1,2，穆桂珍1,2，王秋爽1, 2，蔣慧豪1, 2，王 碩1, 2

(1長江大學(xué)油氣資源與勘探技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室，武漢 430100；2長江大學(xué)資源與環(huán)境學(xué)院，武漢 430100；3 廣東省有色金屬地質(zhì)局940隊(duì)，廣東清遠(yuǎn) 511500；4中國科學(xué)院廣州地球化學(xué)研究所礦物學(xué)與成礦學(xué)重點(diǎn)實(shí)驗(yàn)室，廣州 510640)

開展對(duì)廣東省揭陽市土壤Cr的調(diào)查研究，系統(tǒng)采集了表層(0 ~ 20 cm)土壤樣1 330個(gè)和深層(> 150 cm)土壤樣331個(gè)，測定表層和深層土壤Cr含量，并據(jù)此研究該地區(qū)Cr的空間結(jié)構(gòu)及分布特征。利用富集因子法分析了Cr的污染狀況，結(jié)合Geoda、GS+ 等軟件分析了Cr的結(jié)構(gòu)特征、Moran指數(shù)與空間集聚效應(yīng)，采用GIS空間分析技術(shù)和方差分析對(duì)Cr的空間含量分布以及影響因素進(jìn)行了研究。結(jié)果表明，揭陽市土壤Cr不存在明顯的污染狀況，表層土壤平均Cr含量為22.39 mg/kg，低于揭陽市土壤背景值27.88 mg/kg，遠(yuǎn)低于廣東省土壤背景值50.50 mg/kg。富集因子表明研究區(qū)土壤Cr富集程度多為無富集和輕微富集，分別占到75.11% 和23.46%；Moran指數(shù)值為0.52，表明研究區(qū)有正相關(guān)集聚，“熱點(diǎn)”主要分布在榕城區(qū)及揭東區(qū)的東南部、揭西縣的中南部、惠來縣的北部以及西南部，“冰點(diǎn)”主要分布在揭西縣的西南部以及惠來縣南邊沿海地區(qū)，這與Cr含量的空間分布特征基本一致。土壤Cr主要受到成土母質(zhì)的影響，成土母質(zhì)為粉砂巖和頁巖中的Cr含量較高，沙性母質(zhì)中的Cr含量較低；不同土壤類型中水稻土Cr含量較高，風(fēng)沙土含量較低；不同土地利用方式下的Cr含量差異不明顯，Cr含量并未受到人為活動(dòng)的明顯影響。

土壤鉻；空間結(jié)構(gòu)；分布特征；影響因素；揭陽市

土壤是Cr進(jìn)入生物地球化學(xué)營養(yǎng)鏈的重要媒介，Cr可以通過農(nóng)作物吸收，經(jīng)食物鏈進(jìn)入人體，對(duì)人體造成直接或潛在危害[1]。土壤中Cr含量異常引起農(nóng)作物生長不良、生態(tài)變異和產(chǎn)量下降[2]，在土壤環(huán)境中一般以三價(jià)和六價(jià)兩種形態(tài)存在，其危害程度因化合價(jià)的差異而明顯不同，后者有明顯的致癌致畸作用，同時(shí)Cr也是動(dòng)物和人體所必需微量元素，缺Cr會(huì)出現(xiàn)糖耐量受損，甚至導(dǎo)致糖尿病、動(dòng)脈粥樣硬化和高血壓等疾病[3-4]。有研究表明，土壤中Cr的濃度和空間變化主要受到成土母質(zhì)的影響[5-6]，但并不是唯一的因素，冶煉廠和金屬加工(例如電鍍)以及農(nóng)業(yè)活動(dòng)中施用磷肥等人類活動(dòng)會(huì)導(dǎo)致周圍環(huán)境的Cr濃度增高[7- 8]。

地統(tǒng)計(jì)學(xué)在研究土壤中元素空間分布規(guī)律上得到了廣泛的使用。如武婕等[9]利用地統(tǒng)計(jì)學(xué)對(duì)山東省土壤有機(jī)質(zhì)和微量元素進(jìn)行空間分析時(shí)，發(fā)現(xiàn)在微量元素中, 除B不受影響外, Fe、Mn、Cu和Zn與土壤類型、耕層質(zhì)地、坡度、土地利用類型和地貌類型密切相關(guān)；何博等[10]利用地統(tǒng)計(jì)學(xué)對(duì)典型城市化地區(qū)土壤8種重金屬進(jìn)行空間插值并探討其生態(tài)風(fēng)險(xiǎn)分布，發(fā)現(xiàn)重金屬超標(biāo)地區(qū)集中在研究區(qū)北部。局部Moran’s I是一種在鑒定土壤重金屬熱點(diǎn)、空間自相關(guān)分析方面的有用工具，對(duì)探索重金屬的遷移過程以及影響因素有一定的幫助。揭陽是廣東省發(fā)展速度最快的城市之一，成為粵東地區(qū)新的經(jīng)濟(jì)發(fā)展極，但經(jīng)濟(jì)發(fā)展的同時(shí)，重化、五金、印染紡織等高污染企業(yè)所產(chǎn)生的“三廢”經(jīng)過沉降、風(fēng)化和淋溶進(jìn)入土壤[11]，具有一定的環(huán)境風(fēng)險(xiǎn)，因此，對(duì)土壤中Cr含量的空間結(jié)構(gòu)和分布特征的研究已成為當(dāng)代土壤污染的研究熱點(diǎn)之一。本文以揭陽市為研究區(qū)域，利用地統(tǒng)計(jì)學(xué)方法以及空間自相關(guān)分析技術(shù)，系統(tǒng)地研究該地區(qū)土壤Cr的結(jié)構(gòu)特點(diǎn)和空間分布規(guī)律，為揭示該地區(qū)土壤Cr的富集特征、空間分布特征及影響因素，進(jìn)一步評(píng)價(jià)該區(qū)土壤環(huán)境質(zhì)量提供理論支持。

1 材料與方法

1.1 研究區(qū)域概況

研究區(qū)位于廣東省東南部，陸地面積5 240.5 km2。境內(nèi)屬亞熱帶季風(fēng)性濕潤氣候，光照充足，降水充沛。年均氣溫 21.4 °C，年平均降雨量1 723 mm。研究區(qū)成土母質(zhì)主要由花崗巖、粉砂巖和第四紀(jì)沉積物組成，且以花崗巖及第四紀(jì)沉積物為主。土壤類型以赤紅壤和水稻土為主，北部丘陵和南部沿海地區(qū)分布有少量黃壤和風(fēng)沙土。土地利用方式可分為林地及未利用地、建筑用地和農(nóng)耕用地，農(nóng)作物以水稻為主，是廣東省重要的農(nóng)業(yè)大市。

1.2 樣品采集與分析

表層樣品按照1 km × 1 km的網(wǎng)格單元進(jìn)行布設(shè)采樣，采樣深度為0 ~ 20 cm，1個(gè)樣品由主采樣點(diǎn)周圍100 m范圍內(nèi)3 ~ 5處多點(diǎn)采集組合。深層樣品按照2 km × 2 km的網(wǎng)格單元進(jìn)行布設(shè)采樣，低山丘陵土層較薄地區(qū)，適當(dāng)放稀，但保證采樣大格(4 km × 4 km)有樣點(diǎn)分布，采樣深度為150 ~ 200 cm，人工填土地區(qū)加大采樣深度或移動(dòng)點(diǎn)位。共采集表層土壤組合樣1 330個(gè)，深層母質(zhì)土壤組合樣331個(gè)(圖1)。在農(nóng)業(yè)地區(qū)，采樣點(diǎn)布設(shè)在農(nóng)田、菜地等土層較厚地帶，在城鎮(zhèn)地區(qū)，采樣位置選擇在公園、林地等堆積歷史較長的區(qū)域。避開有明顯點(diǎn)狀污染地段，表層土壤采樣點(diǎn)離主干公路、鐵路100 m以外，且采樣時(shí)間避開施肥期。土壤樣品自然風(fēng)干，除去雜物，研磨后過20目尼龍篩，裝入樣品袋后送樣測試分析。

分析過程嚴(yán)格按照《多目標(biāo)區(qū)域地球化學(xué)調(diào)查規(guī)范》(DD2005-1)的要求進(jìn)行，土壤Cr含量采用X射線熒光光譜儀(S4 PIONEER)，經(jīng)粉末壓片(稱樣4 g、硼酸鑲邊墊底)直接測定，采用離子選擇性電極法(ISE)測定土壤pH，有機(jī)質(zhì)含量采用重鉻酸鉀外加熱法測定。分析過程中所用的試劑均為優(yōu)級(jí)純，所用的水均為超純水(亞沸水)。試驗(yàn)過程中每50個(gè)樣品加入1個(gè)空白樣、4個(gè)土壤標(biāo)準(zhǔn)樣品(GSS-1)和兩個(gè)重復(fù)樣進(jìn)行質(zhì)量控制，以確保精度和準(zhǔn)確性。

1.3 富集因子法

富集因子法是評(píng)價(jià)人類活動(dòng)對(duì)土壤元素富集影響的重要參數(shù)[13]，其中參比元素應(yīng)選用在地殼中穩(wěn)定存在且含量相對(duì)較高的元素[14-16]，本文采用Ti作為參比元素。富集因子的計(jì)算公式為：

式中：CCr 表示 Cr 的實(shí)測濃度(mg/kg)；Cref表示參比元素的實(shí)測濃度(mg/kg)；BCr表示Cr的背景濃度(mg/kg)；Bref表示參比元素的背景濃度(mg/kg)。在此運(yùn)用 Sutherland[17]提出的分類方法將污染元素按不同等級(jí)劃分：EF≤1，無富集；1＜EF≤2，輕微富集；2＜EF≤5, 中度富集；5＜EF≤20，顯著富集；20＜EF≤40，強(qiáng)烈富集；EF＞40，極強(qiáng)富集。

1.4 局部莫蘭指數(shù)法

局部Moran’s I是一種常用的檢驗(yàn)變量局部空間相關(guān)性的重要參數(shù)，Moran’s I可用如下公式進(jìn)行計(jì)算：

1.5 數(shù)據(jù)處理

描述性統(tǒng)計(jì)、Box-Cox轉(zhuǎn)換、Pearson相關(guān)分析和回歸分析等借助SPSS 20.0、Minitab軟件進(jìn)行統(tǒng)計(jì)分析；半變異函數(shù)圖、局部Moran’s I、含量分布圖用GS+、Geoda、ArcGIS 10.2軟件繪制。

2 結(jié)果與討論

2.1 Cr 含量的統(tǒng)計(jì)性描述

揭陽市表層土壤Cr含量的描述性統(tǒng)計(jì)結(jié)果如表1所示。Cr含量變化范圍1.40 ~ 301.90 mg/kg，由于表層土壤Cr含量符合對(duì)數(shù)正態(tài)分布，可用幾何均值表示其平均值，為22.40 mg/kg。揭陽市面積僅占廣東省面積2.92%，用廣東省土壤背景值來衡量揭陽土壤Cr含量情況或許存在偏差，故利用深層土壤有效數(shù)據(jù)計(jì)算揭陽市土壤背景值。深層土壤Cr含量呈偏度分布，經(jīng)Box-Cox轉(zhuǎn)換，Kolmogorov-Smirnov檢驗(yàn)數(shù)據(jù)后符合正態(tài)分布(值=0.20 > 0.05)，其最優(yōu)化轉(zhuǎn)換系數(shù)λ = 0.30。深層Cr含量平均值由Box-Cox的逆運(yùn)算計(jì)算得出，最終算出深層Cr均值為27.88 mg/kg，此值可作為研究區(qū)土壤Cr背景值。由此看出，表層土壤Cr含量22.40 mg/kg較低于揭陽市土壤背景值27.88 mg/kg，遠(yuǎn)低于廣東省土壤背景值50.50 mg/kg。由此可推斷揭陽市土壤不存在明顯的Cr污染問題。Cr的變異系數(shù)為88.7%，屬于中等程度變異[12]，說明Cr含量可能受到一些外源因素的影響。

表1 研究區(qū)土壤Cr含量統(tǒng)計(jì)描述(mg/kg)

2.2 Cr元素的富集程度

以揭陽市土壤Cr 元素背景值為參比值，以Ti 為參比元素，對(duì)所有表層土壤樣品的Cr 元素富集因子(EF)進(jìn)行計(jì)算，結(jié)果見表2。揭陽市表層土壤Cr 富集因子均值為0.81，整體來看處于無污染水平。按照提出的分類方法分類，研究區(qū)土壤Cr 積累程度多為無富集和輕微富集，分別占到75.11% 和23.46%，中度富集占1.20%，極少數(shù)地區(qū)存在顯著富集(0.23%)，無強(qiáng)烈富集和極強(qiáng)富集。

表2 Cr元素所占富集程度百分比

2.3 Cr 含量的空間結(jié)構(gòu)特征

數(shù)據(jù)的非正態(tài)分布會(huì)造成實(shí)驗(yàn)半變異函數(shù)畸變，增大估計(jì)誤差，半變異函數(shù)的波動(dòng)甚至?xí)谏w其固有的結(jié)構(gòu)，因此需轉(zhuǎn)換數(shù)據(jù)消除比例效應(yīng)[18]。本研究中，表層土壤Cr含量經(jīng)對(duì)數(shù)轉(zhuǎn)換后滿足正態(tài)分布，從而消除了比例效應(yīng)。

根據(jù)異性半變異函數(shù)的計(jì)算劃分角度組，容限為±22.5°。研究計(jì)算了揭陽市土壤中Cr的實(shí)驗(yàn)半變異函數(shù)，并根據(jù)其特點(diǎn)進(jìn)行了理論模型的擬合，模型擬合參數(shù)見表3。

表3 土壤Cr理論半變異函數(shù)模型擬合參數(shù)

2檢驗(yàn)達(dá)到顯著(＜0.01)，證明理論變異函數(shù)對(duì)于實(shí)驗(yàn)變異函數(shù)的擬合效果較好。研究區(qū)土壤Cr含量具有很好的空間結(jié)構(gòu)性，表現(xiàn)出明顯的塊金效應(yīng)，塊金值0為0.018，在基臺(tái)值中所占的比例約為25.7%，這部分變異通常被解釋為隨機(jī)效應(yīng)。Cambardella等[19]認(rèn)為0/0+1值可以反映區(qū)域化變量的空間相關(guān)程度，比值小于0.25代表空間相關(guān)性很強(qiáng)，大于0.75代表空間相關(guān)性較弱，因此研究區(qū)土壤Cr含量具有較強(qiáng)的空間相關(guān)性。Cr 的半變異函數(shù)隨著步長的增加而上升，逐漸達(dá)到了基臺(tái)值，變程為225.0 km(圖2)，相關(guān)性范圍基本覆蓋整個(gè)研究區(qū)域。

圖2 研究區(qū)土壤Cr含量的半變異函數(shù)

Fig. 2 Semi-variogram of soil Cr content in studied area

Moran指數(shù)值為0.52 (圖3A)，表明研究區(qū)有正相關(guān)集聚。Moran散點(diǎn)圖表明研究區(qū)表層土壤Cr的集聚類型共有4種，每個(gè)象限代表著一種集聚類型，第一象限的樣點(diǎn)為高-高集聚，即表示這些樣點(diǎn)的Cr含量相對(duì)較高，且周圍均為較高值，該象限的樣點(diǎn)通常被叫做“熱點(diǎn)”；第三象限為低-低集聚，即該象限的樣點(diǎn)Cr含量相對(duì)較低，且周圍均為較低值，該象限的樣點(diǎn)通常被叫做“冰點(diǎn)”。這兩個(gè)象限的共同點(diǎn)均為呈現(xiàn)正向的空間相關(guān)性。Lisa聚類圖(圖3B)中顯示了“冰點(diǎn)”與“熱點(diǎn)”的分布特征，研究區(qū)“熱點(diǎn)”主要分布在榕城區(qū)及揭東區(qū)的東南部、揭西縣的中南部、惠來縣的北部以及西南部；“冰點(diǎn)”主要分布在揭西縣的西南部以及惠來縣南邊沿海地區(qū)。

2.4 Cr的空間分布特征

為了解揭陽市土壤Cr的空間分布特征，利用ArcGIS軟件進(jìn)行普通Kriging插值得到Cr的含量分布圖(圖4)。由圖4可知，Cr含量的空間分布主要受到成土母質(zhì)和地形的控制。表層土壤Cr高值集中在榕城區(qū)和惠來縣的北部以及西南部；土壤Cr低值在揭西縣的西南部以及惠來縣南面沿海等地區(qū)。榕城區(qū)主要位于榕江流域的中段，河流分布密集，成土母質(zhì)主要為粉砂巖與榕江河流沉積物，Cr易隨榕江河流沉積物等進(jìn)行遷移至此，揭陽市地勢自西向東傾斜，西北部多山地丘陵，地勢較高，南部為大面積的沿海沙灘地，不利于Cr的富集，Cr含量相對(duì)偏低，這與前面Lisa聚類圖表現(xiàn)一致。近些年研究表明，在重金屬污染的土壤中，Cr的空間結(jié)構(gòu)和分布主要受成土母質(zhì)的控制，如張慧等[21]通過對(duì)松嫩平原的肇源縣表層土壤重金屬Cr空間分布格局的研究發(fā)現(xiàn)，該地區(qū)土壤中的Cr主要受成土母質(zhì)的控制；呂建樹和何華春[22]對(duì)江蘇沿海帶土壤重金屬分布研究發(fā)現(xiàn)，該地區(qū)土壤中的Cr為自然來源元素，受到成土母質(zhì)的控制。本研究結(jié)果都與此相符。

圖3 研究區(qū)土壤Cr的Moran散點(diǎn)圖(A)和Lisa聚類圖(B)

圖4 表層土壤Cr含量空間分布

2.5 表層土壤Cr 水平空間分布的影響因素

表層土壤重金屬水平空間分布主要受自然因素、人為干擾等因素的影響。其中在自然因素中成土母質(zhì)起主導(dǎo)作用，而土壤類型是成土過程作用的結(jié)果，對(duì)重金屬的分布可能有一定的影響，同時(shí)人類活動(dòng)范圍導(dǎo)致土地性質(zhì)的改變對(duì)局部重金屬分布也會(huì)產(chǎn)生一定的作用。

2.5.1 成土母質(zhì)的影響 成土母質(zhì)是表層土壤主要的物質(zhì)基礎(chǔ)，其類型差異會(huì)直接影響表層重金屬含量[23]。本文利用ANOVA法統(tǒng)計(jì)不同成土母質(zhì)下表層土壤Cr含量(表4)，結(jié)果表明不同母質(zhì)間土壤Cr含量存在顯著性差異。各母質(zhì)中表層Cr平均含量由高到低依次排序?yàn)椋悍凵皫r＞第四紀(jì)沉積物＞頁巖＞凝灰?guī)r＞花崗巖。粉砂巖、頁巖通常是在水動(dòng)力條件較安靜、沉積速度緩慢環(huán)境下形成，分布于湖泊、河流三角洲等地帶，揭陽市中部、南部和東南部都是廣闊肥沃的榕江沖積平原和濱海沉積平原，有利于Cr元素的遷移及在河流沉積物形成的土壤中富集。

表4 不同母質(zhì)類型、土壤類型及土地利用方式下土壤Cr含量(mg/kg)

注：表中小寫字母不同表示同一影響因素不同類型間表層土壤Cr含量差異顯著(< 0.05)。

2.5.2 土壤類型的影響 研究區(qū)土壤類型主要為赤紅壤和水稻土，在北部丘陵和南部沿海區(qū)域分布有少量的黃壤和風(fēng)沙土。用ANOVA法按土壤類型對(duì)表層土壤進(jìn)行統(tǒng)計(jì)(表4)，結(jié)果表明黃壤和風(fēng)沙土之間以及水稻土和黃壤之間Cr含量存在顯著性差異。平均含量由高到低依次排序?yàn)椋核就?赤紅壤>黃壤>風(fēng)沙土。水稻土為水成或半水成土，具水耕淀積層[24]，有利于含Cr物質(zhì)的富集。六價(jià)鉻受pH的影響很大，pH堿性條件下能穩(wěn)定存在，赤紅壤的pH條件下六價(jià)Cr比例一般很低，絕大部分為三價(jià)Cr，而三價(jià)Cr極易被土壤膠體吸附和形成沉淀，以高嶺土、氧化鐵和氧化鋁為主的赤紅壤吸附三價(jià)Cr的能力較強(qiáng)[25]，而以伊利石為主的黃壤吸附三價(jià)Cr的能力較弱。風(fēng)沙土是發(fā)育于風(fēng)成沙性母質(zhì)的土壤，砂性母質(zhì)的流動(dòng)性較強(qiáng)，不利于Cr元素的富集。

2.5.3 土地利用方式的影響 人類通過改變土地利用方式，影響土壤性質(zhì)和土壤生產(chǎn)力，從而影響土壤元素含量和土壤環(huán)境[26]。本文依據(jù)研究區(qū)土地利用的特點(diǎn)，將土地利用類別分為農(nóng)耕用地、建筑用地和林地及未利用地3種方式。ANOVA統(tǒng)計(jì)分析表明這3種土地利用方式下表層土壤Cr含量并沒有顯著性差異且平均含量都小于揭陽市土壤背景值(27.88 mg/kg)。因此可以推斷研究區(qū)土壤Cr含量并未受到明顯的人類活動(dòng)影響。

3 結(jié)論

1)揭陽市表層土壤平均Cr含量為22.39 mg/kg，低于揭陽市土壤背景值27.88 mg/kg，富集因子法分析表明研究區(qū)土壤Cr富集程度多為無富集和輕微富集，分別占到75.11% 和23.46%，整體來看無明顯污染問題。

2)用地統(tǒng)計(jì)學(xué)方法對(duì)揭陽市土壤中Cr進(jìn)行空間分析，Cr含量高值區(qū)分布在成土母質(zhì)以粉砂巖為主的榕城區(qū)中部以及惠城區(qū)北部，低值分布在西北部山地丘陵以及東南部的沿海沙灘地地帶；Moran指數(shù)表明研究區(qū)存在正相關(guān)集聚地區(qū)，Lisa聚類圖可以看出研究區(qū)“熱點(diǎn)”主要分布在榕城區(qū)及揭東區(qū)的東南部、揭西縣的中南部、惠來縣的北部以及西南部，“冰點(diǎn)”主要分布在揭西縣的西南部以及惠來縣南邊沿海地區(qū)，這與土壤Cr的空間分布特征基本一致。

3)方差統(tǒng)計(jì)分析了不同成土母質(zhì)、土壤類型、土地利用方式下土壤Cr含量的影響差異。研究區(qū)土壤Cr主要受到成土母質(zhì)的影響，粉砂巖和頁巖中的Cr含量較高，同時(shí)受到榕江、練江等流域的河流動(dòng)力影響而易于遷移；不同土壤類型中水稻土中Cr含量較高，風(fēng)沙土中Cr含量較低，可能是因?yàn)樗就辆哂兴恋韺右子贑r的富集，而風(fēng)沙土流動(dòng)性較強(qiáng)，Cr難以富集；不同土地利用方式下Cr含量差異不明顯，表明Cr含量并未受到人為活動(dòng)的明顯影響。

[1] 李勇, 周永章, 張澄博, 等. 基于局部Moran's I和GIS的珠江三角洲肝癌高發(fā)區(qū)蔬菜土壤中Ni、Cr的空間熱點(diǎn)分析[J]. 環(huán)境科學(xué), 2010, 31(6): 1617–1623.

[2] Salvador-Blanes S, Cornu S, Bourennane H, et al. Controls of the spatial variability of Cr concentration in topsoils of a central French landscape[J]. Geoderma, 2006, 132(1/2): 143–157.

[3] Brahushi F, D?rfler U, Schroll R, et al. Stimulation of reductive dechlorination of hexachlorobenzene in soil by inducing the native microbial activity[J]. Chemosphere, 2004, 55(11): 1477–1484.

[4] Pohl H R, McClure P R, Fay M, et al. Public health assessment of hexachlorobenzene[J]. Chemosphere, 2001, 43(4–7): 903–908.

[5] Chen M, Ma L Q, Harris W G. Baseline concentrations of 15 trace elements in Florida surface soils[J]. Journal of Environmental Quality, 1999, 28(4): 1173–1181.

[6] Facchinelli A, Sacchi E, Mallen L. Multivariate statistical and GIS-based approach to identify heavy metal sources in soils[J]. Environmental Pollution, 2001, 114(3): 313– 324.

[7] Lindstr?m M. Urban land use influences on heavy metal fluxes and surface sediment concentrations of small lakes[J]. Water, Air, and Soil Pollution, 2001, 126(3/4): 363–383.

[8] 蔡立梅, 馬瑾, 周永章, 等. 珠江三角洲典型區(qū)農(nóng)業(yè)土壤鎳的空間結(jié)構(gòu)及分布特征[J]. 中山大學(xué)學(xué)報(bào)(自然科學(xué)版), 2008, 47(4): 100–104.

[9] 武婕, 李玉環(huán), 李增兵, 等. 南四湖區(qū)農(nóng)田土壤有機(jī)質(zhì)和微量元素空間分布特征及影響因素[J]. 生態(tài)學(xué)報(bào), 2014, 34(6): 1596–1605.

[10] 何博, 趙慧, 王鐵宇, 等. 典型城市化區(qū)域土壤重金屬污染的空間特征與風(fēng)險(xiǎn)評(píng)價(jià)[J]. 環(huán)境科學(xué), 2019, 40(6): 2869–2876.

[11] Liu Y G, Zhang H Z, Zeng G M, et al. Heavy metal accumulation in plants on Mn mine tailings[J]. Pedosphere, 2006, 16(1): 131–136.

[12] 雷志棟, 楊詩秀, 許志榮, 等. 土壤特性空間變異性初步研究[J]. 水利學(xué)報(bào), 1985, 9(9): 10–21.

[13] 李娟娟, 馬金濤, 楚秀娟, 等. 應(yīng)用地積累指數(shù)法和富集因子法對(duì)銅礦區(qū)土壤重金屬污染的安全評(píng)價(jià)[J]. 中國安全科學(xué)學(xué)報(bào), 2006, 16(12): 135–139.

[14] Hasan A B, Kabir S, Selim Reza A H M, et al. Enrichment factor and geo-accumulation index of trace metals in sediments of the ship breaking area of Sitakund Upazilla (Bhatiary-Kumira), Chittagong, Bangladesh[J]. Journal of Geochemical Exploration, 2013, 125: 130–137.

[15] Fang G C, Wu Y S, Chang S Y, et al. Size distributions of ambient air particles and enrichment factor analyses of metallic elements at Taichung Harbor near the Taiwan Strait[J]. Atmospheric Research, 2006, 81(4): 320–333.

[16] Haris H, Aris A Z. The geoaccumulation index and enrichment factor of mercury in mangrove sediment of Port Klang, Selangor, Malaysia[J]. Arabian Journal of Geosciences, 2013, 6(11): 4119–4128.

[17] Sutherland R A. Bed sediment-associated trace metals in an urban stream, Oahu, Hawaii[J]. Environmental Geology, 2000, 39(6): 611–627.

[18] 張朝生, 章申, 何建邦. 長江水系沉積物重金屬含量空間分布特征研究——空間自相關(guān)與分形方法[J]. 地理學(xué)報(bào), 1998, 53(1): 87.

[19] Cambardella C A, Moorman T B, Parkin T B, et al. Field-scale variability of soil properties in central Iowa soils[J]. Soil Science Society of America Journal, 1994, 58(5): 1501–1511.

[20] 李秋萍, 李長建, 肖小勇, 等. 中國農(nóng)業(yè)碳排放的空間效應(yīng)研究[J]. 干旱區(qū)資源與環(huán)境, 2015, 29(4): 30–35.

[21] 張慧, 馬鑫鵬, 史曉磊, 等. 區(qū)域土壤Cd和Cr空間分布的影響因素研究[J]. 土壤, 2018, 50(5): 989–998.

[22] 呂建樹, 何華春. 江蘇海岸帶土壤重金屬來源解析及空間分布[J]. 環(huán)境科學(xué), 2018, 39(6): 2853–2864.

[23] 張慧, 鄭志志, 楊歡, 等. 基于多元統(tǒng)計(jì)和地統(tǒng)計(jì)的肇源縣表層土壤重金屬來源辨析[J]. 土壤, 2017, 49(4): 819–827.

[24] 趙其國. 水稻土的類型特征及其管理[J]. 土壤, 1992, 24(4): 169–175.

[25] 周琴, 姜軍, 徐仁扣. Cu(Ⅱ)、Pb(Ⅱ)和Cd(Ⅱ)在紅壤膠體和非膠體顆粒上吸附的比較[J]. 土壤學(xué)報(bào), 2018, 55(1): 131–138.

[26] 吳文斌, 楊鵬, 唐華俊, 等. 土地利用對(duì)土壤性質(zhì)影響的區(qū)域差異研究[J]. 中國農(nóng)業(yè)科學(xué), 2007, 40(8): 1697– 1702.

Spatial Structure and Distribution of Soil Chromium in Jieyang City, Guangdong Province

WANG Hanzhi1,2, WEN Hanhui3, CAI Limei1, 2, 4*, LUO Jie1,2, MU Guizhen1,2, WANG Qiushuang1, 2, JIANG Huihao1, 2, WANG Shuo1, 2

(1 Key Laboratory of Exploration Technologies for Oil and Gas Resources, Ministry of Education, Yangtze University, Wuhan 430100, China; 2 College of Resources and Environment, Yangtze University, Wuhan 430100, China; 3 Non-ferrous Metals Geological Bureau 940 Team, Guangdong Province, Qingyuan, Guangdong 511500, China; 4 Key Laboratory of Mineralogy and Metallogeny, Guangzhou Institute of Geochemistry, Chinese Academy of Sciences, Guangzhou 510640, China)

A systematic investigation of soil chromium (Cr) in Jieyang City of Guangdong Province was carried out, 1 330 surface soil samples and 331 deep soil samples were collected, the content, spatial structure and distribution Cr were measured and analyzed, Cr pollution status was analyzed using enrichment factor method, structural characteristics, Moran index and spatial agglomeration effect of Cr were analyzed by software such as Geoda and GS+, and spatial distribution and influence factors of Cr were analyzed by GIS and variance analysis. The results showed that there was no obvious soil Cr pollution in the studied area, the average Cr content of surface soil was 22.39 mg/kg, which was lower thanthe background value of 27.88 mg/kg in Jieyang City, and much lower than the background value of 50.50 mg/kg in Guangdong Province. Cr accumulation in the studied area was mostly non-accumulation and slight accumulation, accounting for 75.11% and 23.46%, respectively. Moran index was 0.52, indicating that the studied area had positive correlation aggregations. The “hot spots” were mainly distributed in the southeastern part of Rongcheng and Jiedong districts, the south-central part of Jiexi County, the north and southwest of Huilai County; the “ice spots” were mainly distributed in the southwest of Jiexi County and the coastal area south of Huilai County, which were basically consistent with the spatial distribution of Cr content. Soil Cr was mainly affected by the parent material, higher Cr in siltstone and shale but lower in sandy parent material. Cr content in paddy soil was higher than other soil types, and lower in sandy soil, but no obvious difference between different land use types, and human activities affected Cr content insignificantly.

Soil Cr; Spatial structure; Distribution characteristics; Influencing factor; Jieyang City

X825

A

10.13758/j.cnki.tr.2020.05.010

王涵植, 溫漢輝, 蔡立梅, 等. 廣東省揭陽市土壤鉻的空間結(jié)構(gòu)及分布特征研究. 土壤, 2020, 52(5): 941–947.

國家自然科學(xué)基金項(xiàng)目(41203061)、教育部油氣資源勘探技術(shù)重點(diǎn)實(shí)驗(yàn)室開放基金項(xiàng)目(K2018-19)、湖北省自然科學(xué)基金項(xiàng)目(2015CFB603)、湖北省教育廳科學(xué)技術(shù)研究重點(diǎn)項(xiàng)目(D20161301)、有機(jī)地球化學(xué)國家重點(diǎn)實(shí)驗(yàn)室開放基金項(xiàng)目(OGL-201408)和長江大學(xué)大學(xué)生創(chuàng)新創(chuàng)業(yè)訓(xùn)練項(xiàng)目(2016006,2017223)資助。

王涵植(1995—)，男，湖北鄂州人，碩士研究生，主要研究方向?yàn)榄h(huán)境地球化學(xué)。E-mail: cjdxwhz@163.com