王俊杰,黃陽,黃雪,3,左平,趙善道,鄒欣慶
(1.南京大學(xué) 地理與海洋科學(xué)學(xué)院, 江蘇 南京 210093;2.南京大學(xué) 海岸與海島開發(fā)教育部重點(diǎn)實(shí)驗(yàn)室,江蘇 南京 210093;3.北京師范大學(xué) 地理學(xué)與遙感科學(xué)學(xué)院,北京 100875)
江蘇鹽城國(guó)家級(jí)珍禽自然保護(hù)區(qū)又稱鹽城生物圈保護(hù)區(qū)(簡(jiǎn)稱鹽城保護(hù)區(qū)),是中國(guó)最大的濱海濕地自然保護(hù)區(qū),位于江蘇中部沿海,自北向南依次地跨響水、濱海、射陽、大豐、東臺(tái)五縣市的沿海灘涂,其地理范圍是32°20′N-34°37′N,119°29′E-121°16′E(李楊帆 等,2004)。其核心區(qū)是我國(guó)目前保留最完好的原生鹽沼濕地之一,自陸向海規(guī)律地分布有蘆葦灘、鹽蒿灘、互花米草灘、光灘(歐維新 等,2006;趙雪琴等,2010a;趙雪琴等,2010b)。但海岸帶生態(tài)系統(tǒng)作為開放的生態(tài)系統(tǒng),由于各類經(jīng)濟(jì)活動(dòng)強(qiáng)度的增大,導(dǎo)致其產(chǎn)生的各種工農(nóng)業(yè)污染物隨廢水排海、大氣顆粒沉降、地表徑流和降雨進(jìn)入近岸水體,并不可避免地隨近岸河流流入核心區(qū)(Frey et al,1978;王俊杰 等,2013;左平等,2010;黃向青等,2006;邱虎等,2010)。而重金屬污染在土壤中移動(dòng)性差,滯留時(shí)間長(zhǎng),不能被微生物降解(高園園等,2013),具有一定的毒性,通過食物鏈積累富集進(jìn)入魚類、鳥類體內(nèi),從而影響整個(gè)生態(tài)系統(tǒng)的健康(王俊杰等,2013;左平等,2010)。因此,本文基于鹽城自然保護(hù)區(qū)核心區(qū)外光灘的表層沉積物中的重金屬含量及其分布特征進(jìn)行分析,旨在為濱海保護(hù)區(qū)的開發(fā)利用和生態(tài)環(huán)境保護(hù)提供科學(xué)依據(jù)。
本研究于2011年9月23-25日自江蘇鹽城自然保護(hù)區(qū)核心區(qū)外圍光灘采集表層沉積物樣品,研究區(qū)域及采樣點(diǎn)位置見圖1。

圖1 研究區(qū)域及采樣點(diǎn)位置分布圖
從新洋港入海口的光灘地區(qū),自東南向西北整體沿平行于海岸線方向近似等間距取表層沉積物樣品25 個(gè)。其中樣品GT01-GT13 呈折線分布,折向?yàn)闁|北—西南,走向?yàn)橛蓶|南到西北;GT14-GT25依次按東南至西北的方向呈直線分布。由此將采樣的光灘區(qū)域分成5 個(gè)區(qū)段,前4 個(gè)區(qū)段均在折線范圍內(nèi):GT01-GT04 位于由西南到東北,即由陸向海采樣區(qū)的第1 段;第2 段是東北到西南方向的,即由海向陸的包括GT05-GT07 3 個(gè)樣點(diǎn);第3 段包含GT08-GT10,方向由西南到東北,即由陸向海的;GT11-GT13 位于第4 段,方向與第3 段相反。第5 段是平行于岸線的直線部分的GT14-GT25。現(xiàn)場(chǎng)用塑料自封袋密封,冷藏,帶回實(shí)驗(yàn)室。樣品分3 份,1 份用于粒度分析,1 份用于Cr、Ni、Cu、Zn、Cd、Pb 6 種重金屬分析,剩余一份用作實(shí)驗(yàn)室備份。
粒度測(cè)定:將25 個(gè)粒度實(shí)驗(yàn)樣品稱重,過16目篩,將篩出的≥1 mm 的顆粒稱重,計(jì)算所占樣品含量比例。過篩后的樣品預(yù)處理后使用Malvern Master Size 2000 型激光粒度儀進(jìn)行粒度分析。該儀器的測(cè)量范圍為0.02 ~ 2 000 μm,粒級(jí)分辨率為0.01Φ,重復(fù)測(cè)量的相對(duì)誤差<3%。使用矩法計(jì)算沉積物粒度參數(shù)(McManus,1988)。所有的樣品分析均在南京大學(xué)海岸與海島開發(fā)教育部重點(diǎn)實(shí)驗(yàn)室完成。
重金屬測(cè)定:樣品放入冰箱中冷凍,然后利用ALPHA-1-4 型凍干機(jī)低溫凍干,用瑪瑙研缽研磨后過孔徑100 μm 的尼龍篩。采用國(guó)標(biāo)濕消化法消解泥樣,用國(guó)家標(biāo)準(zhǔn)物質(zhì)進(jìn)行質(zhì)量控制以確保數(shù)據(jù)可靠性。用ICP-MS 元素分析儀測(cè)定樣品Cd、Cr、Cu、Pb、Zn、Ni 共6 個(gè)指標(biāo)(相對(duì)誤差<5%,回收率在95%~125%之間)。樣品分析在中國(guó)科學(xué)院南京地理與湖泊研究所國(guó)家重點(diǎn)實(shí)驗(yàn)室完成。
克里格差值分析:又稱空間局部插值法,是以變異函數(shù)理論和結(jié)構(gòu)分析為基礎(chǔ),在有限區(qū)域內(nèi)對(duì)區(qū)域化變量進(jìn)行無偏最優(yōu)估計(jì)的一種方法,是地統(tǒng)計(jì)學(xué)的主要內(nèi)容之一(王艷妮等,2009)。克里格法適用的條件:變異函數(shù)和相關(guān)分析的結(jié)果表明區(qū)域化的變量存在空間相關(guān)性。克里格方法與反距離權(quán)插值方法類似的是,兩者都通過對(duì)已知樣本點(diǎn)賦權(quán)重來求得未知樣點(diǎn)的值,可統(tǒng)一表示為:

式(1) 中,Z(x0)為未知樣點(diǎn)的值,Z(xi)為未知樣點(diǎn)周圍的已知樣本點(diǎn)的值,ωi為第i 個(gè)樣本點(diǎn)對(duì)未知樣點(diǎn)的權(quán)重。n 為已知樣本點(diǎn)的個(gè)數(shù)。
環(huán)境風(fēng)險(xiǎn)指數(shù)法:對(duì)污染環(huán)境進(jìn)行環(huán)境風(fēng)險(xiǎn)表征的方法,其規(guī)定了相應(yīng)的環(huán)境風(fēng)險(xiǎn)的劃分標(biāo)準(zhǔn),可以定量地度量重金屬污染的土壤或沉積物中樣品的環(huán)境風(fēng)險(xiǎn)程度大小(Rapant et al,2003)。計(jì)算公式如下:

式(2) 中,IERi為超過臨界限量的第i 種元素的環(huán)境風(fēng)險(xiǎn)指數(shù);ACi為第i 種元素的分析含量;RCi為第i 種元素的臨界限量;IER為待測(cè)樣品的環(huán)境風(fēng)險(xiǎn)。需要說明的是,如果ACi 環(huán)境風(fēng)險(xiǎn)指數(shù)法能定量反映重金屬污染風(fēng)險(xiǎn)程度的大小,能用數(shù)值來反映污染物對(duì)環(huán)境現(xiàn)狀的危害程度,Rapant 等應(yīng)用環(huán)境風(fēng)險(xiǎn)指數(shù)法對(duì)斯洛伐克共和國(guó)的環(huán)境進(jìn)行了風(fēng)險(xiǎn)分級(jí),分析了各種重金屬對(duì)環(huán)境污染的貢獻(xiàn)程度和對(duì)環(huán)境污染貢獻(xiàn)最大的重金屬元素,但這種方法不能反映出重金屬污染在這個(gè)時(shí)間和空間的變化特征(徐燕等,2008;范拴喜等,2010)。 Hakanson 潛在生態(tài)危害指數(shù)法:綜合考慮了重金屬的毒性、在土壤和沉積物中普遍的遷移轉(zhuǎn)化規(guī)律和評(píng)價(jià)區(qū)域?qū)χ亟饘傥廴镜拿舾行裕约爸亟饘賲^(qū)域背景值的差異,消除了區(qū)域差異影響,體現(xiàn)了生物有效性和相對(duì)貢獻(xiàn)及地理空間差異等特點(diǎn),是綜合反映重金屬對(duì)生態(tài)環(huán)境影響潛力的指標(biāo),適合于大區(qū)域范圍沉積物和土壤進(jìn)行評(píng)價(jià)比較。該方法是目前國(guó)內(nèi)外沉積物和土壤質(zhì)量評(píng)價(jià)中應(yīng)用最為廣泛的方法之一(Hakanson,1980;徐爭(zhēng)啟等,2008;徐燕等,2008;方淑波等,2012;李桂海等,2007;張麗旭等,2005;張少峰 等,2010;張淑娜等,2008)。其計(jì)算公式如下(Hakanson,1980): 式(4)、式(5) 中:Ei 表1 Hakanson 重金屬毒性系數(shù)表 自東南至西北平行于岸線的GT01-GT25 樣點(diǎn)表層沉積物的粒度特征如圖2 所示,其中實(shí)線(折線) 表示樣點(diǎn)的實(shí)際粒度的分布特征,虛線表示相應(yīng)折線的趨勢(shì)線。 圖2 平行海岸線由南向北表層沉積物粒度變化 研究區(qū)域內(nèi)沉積物種類是由粘土、粉砂和砂組成的;整體分布以粉砂為主,粘土次之,體現(xiàn)出鹽城保護(hù)區(qū)的淤泥質(zhì)海灘特征;而砂的含量則很少,大部分樣點(diǎn)均小于1%,僅在第2 段出現(xiàn)砂含量大于1%的樣點(diǎn),但也在10%以內(nèi)。由粉砂含量的趨勢(shì)線可見,第1 段至第5 段呈現(xiàn)粉砂含量減小,粘土含量增加,平均粒徑減小的趨勢(shì)。造成這種趨勢(shì)可能是因?yàn)楸眰?cè)的樣點(diǎn)GT25 離新洋港入海口較近,河流流速較緩,加上沿河植物的過濾作用,使得河流輸送的細(xì)粒物質(zhì)多沉積于此。同時(shí)此處離海距離較遠(yuǎn),水動(dòng)力條件較為單一,沉積物粒度受潮汐波浪作用較弱。西北GT25 附近的區(qū)域,由于水動(dòng)力環(huán)境穩(wěn)定,并且接收了大量來自河流的營(yíng)養(yǎng)物質(zhì),成為當(dāng)?shù)貤l件較好的泥螺養(yǎng)殖區(qū)。而隨著東南部的樣點(diǎn)離河流的距離增大,受到潮汐動(dòng)力的影響逐漸增大所致,因此平均粒徑呈現(xiàn)增大趨勢(shì)。而在GT07 處,砂和粉砂含量異常增高,表明此處水動(dòng)力較強(qiáng),能量較高,可能是波浪破碎的位置。 圖3-8 分別對(duì)應(yīng)自東南向西北的光灘地區(qū)表層沉積物中Cr、Cu、Ni、Zn、Cd 和Pb 6 種重金屬的分布情況,其中實(shí)線(折線) 代表各采樣點(diǎn)的實(shí)際重金屬含量,短虛線(曲線) 代表其趨勢(shì)線,長(zhǎng)虛線(直線) 代表該種重金屬在該采樣區(qū)域所對(duì)應(yīng)的背景值(江蘇省沿海海域及港河口污染調(diào)查監(jiān)測(cè)總結(jié)報(bào)告1978-1979)。由圖3 可見,所有樣點(diǎn)的Cr 含量均遠(yuǎn)遠(yuǎn)高于背景值(28.50 mg/kg),甚至達(dá)到其2~4 倍之多;由東南至西北呈現(xiàn)中部低、兩端高的趨勢(shì);最大值出現(xiàn)在第2 段的GT07 點(diǎn),最小值出現(xiàn)在第五段的GT21 點(diǎn)。其原因可能是Cr以海上運(yùn)輸為主要來源,通過潮汐波浪的作用向岸輸送物質(zhì)的過程中,一開始流速較大,粗粒很快沉積;隨著流速下降,離海距離增大,較細(xì)顆粒逐漸沉降下來。假設(shè)粗粒與細(xì)粒密度相同,且單位體積Cr 含量相同,則離海距離越遠(yuǎn),表層樣中Cr 含量越小,因此由第一段至第4 段呈現(xiàn)減小的趨勢(shì)。Ni、Zn、Cd、Pb 含量分布可分別見圖5、圖6、圖7 和圖8,其含量的整體分布趨勢(shì)較為一致,第1段至第2 段呈減小趨勢(shì),第2 段至第3 段呈增大趨勢(shì),第3 段至第4 段呈減小趨勢(shì),在第4 段達(dá)到最低值,第4 段至第5 段呈現(xiàn)增加趨勢(shì)。其中Ni 含量的分布地區(qū)差異最大,且所有點(diǎn)的含量均超過其背景值(20.00 mg/kg),最大值出現(xiàn)在第2 段的GT07 點(diǎn);Pb 最高值與最低值之間差異相對(duì)較小,且除了第4 段和第5 段的前半段以外,其他采樣點(diǎn)的Pb 含量均超過其背景值(17.63 mg/kg),最高點(diǎn)在第3 段GT08 點(diǎn),最低在第4 段GT12 點(diǎn);Zn 和Cd 的地區(qū)差異為最小,且所有采樣點(diǎn)的Zn、Cd含量均低于其各自背景值(Zn:78.10 mg/kg,Cd:0.72 mg/kg);Zn 最高值出現(xiàn)在第3 段的GT10 和第2 段的GT07 點(diǎn),Cd 最高值也在第2 段的GT07點(diǎn)。這4 種重金屬最高點(diǎn)出現(xiàn)在粒度分布中粘土含量最高處附近,是因?yàn)槠骄皆叫。w粒對(duì)重金屬的吸附作用越強(qiáng)。由圖4 可見,除了含量最高的第1 段GT04 點(diǎn)以外,其他點(diǎn)的Cu 含量均低于其背景值(15.99 mg/kg);以第4 段含量為最低,向南北含量呈增加趨勢(shì),第3 段和第4 段Cu 含量的變化較緩,第1 段和第5 段變化較為劇烈。 圖3 鹽城濱海灘涂東南至西北一線Cr 含量 圖4 鹽城濱海灘涂東南至西北一線Cu 含量分布 圖5 鹽城濱海灘涂東南至西北一線Ni 含量分布 圖6 鹽城濱海灘涂東南至西北一線Zn 含量分布 圖7 鹽城濱海灘涂東南至西北一線Cd 含量分布 圖8 鹽城濱海灘涂東南至西北一線Pb 含量分布 由插值結(jié)果可見(圖9-14),除了Cu(圖10)和Ni(圖11) 以外的其他4 種重金屬元素含量的最高值均集中在東南部的GT01-GT08 處(第1 段至第3 段),其中Cr、Cd 和Pb 在東南部區(qū)域分布均勻,呈現(xiàn)由南到北、由岸向海減小的趨勢(shì);而Zn 的含量分布較為不均,第3 段明顯高于南部的前兩段。Cu、Ni 則在西北和東南兩端含量均為較高,其中西北岸(第5 段后半段) 范圍更大且含量更高;Cu 的高值多以點(diǎn)狀分散分布,Ni 則以片狀分布。6 種重金屬元素空間分布的共同特點(diǎn)是低值都分布在第4 段和第5 段前半段。從實(shí)際地理位置分析可知,南部的GT01-GT08 點(diǎn)靠近漁民下海的通道,附近發(fā)現(xiàn)有大量生活生產(chǎn)垃圾,多為漁民從事漁業(yè)活動(dòng)后遺棄的勞作輔助工具等,這些垃圾經(jīng)過海水的腐蝕、分解產(chǎn)生大量的重金屬元素,導(dǎo)致該地區(qū)重金屬含量異常增高;北部的GT18-GT25點(diǎn)靠近新洋港入海口,獲得了較多河流輸送的重金屬元素,并且作為泥螺養(yǎng)殖的主要作業(yè)區(qū),由于漁民長(zhǎng)時(shí)間的勞作過程中帶入和產(chǎn)生了較多污染,導(dǎo)致該區(qū)域的重金屬含量普遍較高。 圖9 鹽城濱海灘涂Cr 含量空間分布(mg/kg) 圖10 鹽城濱海灘涂Cu 含量空間分布(mg/kg) 圖11 鹽城濱海灘涂Ni 含量空間分布(mg/kg) 圖12 鹽城濱海灘涂Zn 含量空間分布(mg/kg) 圖13 鹽城濱海灘涂Cd 含量空間分布(mg/kg) 圖14 鹽城濱海灘涂Pb 含量空間分布(mg/kg) 參照《中華人民共和國(guó)海洋沉積物質(zhì)量》(GB18668-2002),對(duì)鹽城濱海光灘6 種重金屬進(jìn)行比較發(fā)現(xiàn),Cu、Zn、Cd、Pb 4 種重金屬含量均低于國(guó)家I 類標(biāo)準(zhǔn);Cr 含量的均值高于國(guó)家I 類標(biāo)準(zhǔn),但低于Ⅱ類標(biāo)準(zhǔn),因此該沉積物屬國(guó)家Ⅱ類沉積物。Ni 含量近于其背景值,因沒有相關(guān)分類標(biāo)準(zhǔn),此處不做討論。 采用環(huán)境風(fēng)險(xiǎn)指數(shù)法對(duì)江蘇鹽城濱海灘涂表層沉積物中Cr、Ni、Cu、Pb 和綜合重金屬污染風(fēng)險(xiǎn)評(píng)價(jià)見圖15,其中Zn 和Cd 的環(huán)境風(fēng)險(xiǎn)指數(shù)全部為0,故未在圖中畫出。 圖15 鹽城濱海灘涂環(huán)境風(fēng)險(xiǎn)指數(shù)(IER) 由圖可見,對(duì)整個(gè)環(huán)境風(fēng)險(xiǎn)貢獻(xiàn)最大的元素是Cr,是鹽城保護(hù)區(qū)光灘區(qū)最主要的重金屬污染因素;Pb 和Ni 有不同程度的輕度污染;Zn、Cd 和Cu 的環(huán)境風(fēng)險(xiǎn)指數(shù)均在正常范圍內(nèi)。6 種重金屬元素的環(huán)境風(fēng)險(xiǎn)指數(shù)排序是Cr>Ni>Pb>Cu>Zn=Cd=0。總風(fēng)險(xiǎn)指數(shù)以第1 段和第2 段最大,第3 段、第4 段和第5 段前半段較小。 采用Hakanson 潛在生態(tài)危害指數(shù)法對(duì)江蘇鹽城保護(hù)區(qū)光灘表層沉積物中Cr、Ni、Cu、Zn、Cd、Pb 共6 種重金屬的綜合潛在生態(tài)危害評(píng)價(jià)結(jié)果見圖16。 圖16 鹽城濱海灘涂潛在生態(tài)危害指數(shù)(RI) 表2 潛在生態(tài)危害系數(shù)和危害指數(shù)與污染程度的關(guān)系 由表1 可知,毒性系數(shù)最大的重金屬元素為Cd,最小的是Zn。6 種重金屬在整個(gè)采樣區(qū)的潛在生態(tài)風(fēng)險(xiǎn)指數(shù)平均值的排序?yàn)?0>Cd>Cr>Ni>Pb>Cu>Zn,參考表2 潛在生態(tài)危害系數(shù)與污染程度的關(guān)系可知,遠(yuǎn)小于最低標(biāo)準(zhǔn)40,表明該區(qū)域的六種重金屬全部屬于非常輕微的生態(tài)危害;6 種重金屬的綜合潛在生態(tài)危害指數(shù)RI 值均在40 以下,遠(yuǎn)小于150,總體也屬非常輕微程度的生態(tài)危害。由圖16 可知,由東南到西北方向,RI 值整體呈現(xiàn)先增大、后減小、再增大的趨勢(shì),以第1 段和第2 段最大,第3 段、第4 段和第5 段前半段較小,與環(huán)境風(fēng)險(xiǎn)指數(shù)的變化規(guī)律大致相同,可能是由于東南部的第1、2 段靠近漁業(yè)活動(dòng)通道,受人類活動(dòng)干擾較大;北部的第5 段后半段由于近新洋港入海口,獲得較多來自陸地的重金屬元素,且受到該區(qū)泥螺養(yǎng)殖的經(jīng)濟(jì)活動(dòng)也有所影響。由此看出,人類生產(chǎn)活動(dòng)對(duì)重金屬分布產(chǎn)生了直接的影響,人類活動(dòng)較多的區(qū)域重金屬含量要高于人類活動(dòng)少的。 (1) 鹽城濱海濕地灘涂東南至西北表層沉積物中,Cr、Ni、Pb 的含量均高于背景值;所有樣品的Cr 含量都超過了背景值達(dá)1-4 倍,且高于國(guó)家I 類標(biāo)準(zhǔn),低于Ⅱ類標(biāo)準(zhǔn),屬Ⅱ類沉積物;Cu、Zn、Cd、Pb 4 種重金屬含量均低于國(guó)家I 類標(biāo)準(zhǔn),屬國(guó)家I 類沉積物。 (2) 從空間插值結(jié)果來看,Cr、Cd、Zn、Pb高值集中于東南部,表明東南部的漁業(yè)活動(dòng)及其產(chǎn)生的垃圾為其主要重金屬來源;Cu、Ni 高值在西北和東南兩處,表明其來源除南部的漁業(yè)活動(dòng)外,還有北部的陸源輸入。綜合表明,研究區(qū)重金屬的來源和富集主要受人類活動(dòng)控制。 (3) 采用環(huán)境風(fēng)險(xiǎn)指數(shù)法進(jìn)行分析發(fā)現(xiàn),Cr對(duì)整個(gè)區(qū)域的環(huán)境風(fēng)險(xiǎn)的影響最大,總風(fēng)險(xiǎn)指數(shù)以東南部第1、2 段最大,中部第3、4 段和第五段前半段最小;采用Hakanson 潛在生態(tài)風(fēng)險(xiǎn)指數(shù)法對(duì)數(shù)據(jù)進(jìn)行分析評(píng)價(jià),結(jié)果表明Cd 的潛在生態(tài)風(fēng)險(xiǎn)最大,RI 值由東南到西北方向整體呈現(xiàn)先增大、后減小、再增大的趨勢(shì),與綜合環(huán)境風(fēng)險(xiǎn)指數(shù)的變化規(guī)律大致一致,以新洋港入海口和漁業(yè)活動(dòng)最密集處為最高。 沉積物粒度特征、重金屬的元素特性、陸源輸入或海水輸運(yùn)、潮汐波浪動(dòng)力作用、生物化學(xué)作用等對(duì)鹽城保護(hù)區(qū)潮間帶光灘表層沉積物中所研究的6 種重金屬的分布和轉(zhuǎn)移均具有不同程度的影響,直接關(guān)系到重金屬在沉積物中的污染程度。通過新洋港河輸運(yùn)的各類陸上污染物已經(jīng)隨地表徑流進(jìn)入江蘇鹽城國(guó)家級(jí)自然保護(hù)區(qū)的核心區(qū),海上輸運(yùn)產(chǎn)生的污染物也隨著潮汐波浪作用廣泛影響鹽城濱海潮間帶光灘區(qū)域。近些年來蘇北沿海經(jīng)濟(jì)開發(fā)的范圍和強(qiáng)度逐年增大,耕地開墾、工業(yè)園區(qū)、城市化建設(shè)進(jìn)度不斷加快,三廢的產(chǎn)生量不斷增加,而處理能力相對(duì)滯后(左平,2012)。只有更嚴(yán)格地規(guī)范各類污染物的排放和治理,才能保障核心區(qū)生態(tài)系統(tǒng)的健康。由于生態(tài)系統(tǒng)物質(zhì)循環(huán)具有區(qū)域廣泛性的特點(diǎn),在嚴(yán)禁核心區(qū)內(nèi)進(jìn)行各種人類活動(dòng)的基礎(chǔ)上,還要對(duì)核心區(qū)外圍區(qū)域的各類經(jīng)濟(jì)活動(dòng),尤其是工業(yè)生產(chǎn)活動(dòng),進(jìn)一步加強(qiáng)管理,控制污染源,如限制化肥農(nóng)藥用量、限制灘涂工業(yè)污染項(xiàng)目,生產(chǎn)過程中增設(shè)監(jiān)督機(jī)制,實(shí)施達(dá)標(biāo)排放,避免直接排海。 Frey R W,Basan P B,1978.Coastal salt marshes.In Coastal sedimentary environments (pp.101-169) .Springer US. Hakanson L,1980.An ecological risk index for aquatic pollution control.A sedimentological approach.Water research,14 (8) :975-1001. McManus J, 1988. Grain size determination and interpretation. Techniques in sedimentology,408. Rapant S, Kordik J, 2003. An environmental risk assessment map of the Slovak Republic: application of data from geochemical atlases.Environmental Geology,44 (4) :400-407. 范拴喜,甘卓亭,李美娟,等,2010.土壤重金屬污染評(píng)價(jià)方法進(jìn)展.中國(guó)農(nóng)學(xué)通報(bào),26(17):310-315. 方淑波,賈曉波,安樹青,等,2012.鹽城海岸帶土壤重金屬潛在生態(tài)風(fēng)險(xiǎn)控制優(yōu)先格局.地理學(xué)報(bào),67(1) . 高園園,周啟星,2013.納米零價(jià)鐵在污染土壤修復(fù)中的應(yīng)用與展望.農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào),32(3):418-425. 黃向青,張順枝,霍振海,2006.深圳大鵬灣,珠江口海水有害重金屬分布特征.海洋湖沼通報(bào),(4):38-44. 江蘇省沿海海域及港河口污染調(diào)查監(jiān)測(cè)協(xié)作組,1979.江蘇省沿海海域及港河口污染調(diào)查監(jiān)測(cè)總結(jié)報(bào)告(1978-1979) . 李桂海,藍(lán)東兆,曹志敏,等,2007.廈門海域沉積物中的重金屬及其潛在生態(tài)風(fēng)險(xiǎn).海洋通報(bào),26(1):67-72. 李楊帆,朱曉東,鄒欣慶,等,2004.鹽城海岸濕地資源環(huán)境壓力與生態(tài)調(diào)控響應(yīng).自然資源學(xué)報(bào),19(6):754-760. 歐維新,楊桂山,高建華,2006.鹽城潮灘濕地對(duì)N,P 營(yíng)養(yǎng)物質(zhì)的截留效應(yīng)研究.濕地科學(xué),4(3):179-186. 邱虎,呂惠進(jìn),2010.江蘇鹽城濱海濕地現(xiàn)狀與保護(hù)對(duì)策研究.湖南農(nóng)業(yè)科學(xué):上半月,(11):58-61. 王俊杰,左平,黃陽,等,2013.鹽城新洋港表層沉積物重金屬分布特征與評(píng)價(jià).環(huán)境保護(hù)科學(xué),(5):57-62. 王艷妮,謝金梅,郭祥,2009.ArcGIS 中的地統(tǒng)計(jì)克里格插值法及其應(yīng)用.教育技術(shù)導(dǎo)刊,7(12):36-38. 徐燕,李淑芹,郭書海,等,2008.土壤重金屬污染評(píng)價(jià)方法的比較.安徽農(nóng)業(yè)科學(xué),36(11):4 615-4 617. 徐爭(zhēng)啟,倪師軍,庹先國(guó),等,2008.潛在生態(tài)危害指數(shù)法評(píng)價(jià)中重金屬毒性系數(shù)計(jì)算.環(huán)境科學(xué)與技術(shù),31(2):112-115. 張麗旭,任松,蔡健,2005.東海三個(gè)傾倒區(qū)表層沉積物重金屬富積特征及其潛在生態(tài)風(fēng)險(xiǎn)評(píng)價(jià).海洋通報(bào),24(2):92-96. 張少峰,林明裕,魏春雷,等,2010.廣西欽州灣沉積物重金屬污染現(xiàn)狀及潛在生態(tài)風(fēng)險(xiǎn)評(píng)價(jià).海洋通報(bào),29(4):450-454. 張淑娜,唐景春,2008.天津海域沉積物重金屬潛在生態(tài)風(fēng)險(xiǎn)危害評(píng)價(jià).海洋通報(bào),27(2):85-90. 趙雪琴,趙善道,左平,等,2010a.江蘇鹽城原生濕地表層沉積物中的重金屬污染評(píng)價(jià).環(huán)境保護(hù)科學(xué),36(1):64-68. 趙雪琴,趙善道,左平,等,2010b.江蘇鹽城原生鹽沼濕地表層水體Hg,As 分布特征及評(píng)價(jià).三峽環(huán)境與生態(tài),(1):17-20. 左平,李云,趙書河,等,2012.1976年以來江蘇鹽城濱海濕地景觀變化及驅(qū)動(dòng)力分析.海洋學(xué)報(bào),29(4):372-377.

2 結(jié)果與討論
2.1 粒度數(shù)據(jù)分析

2.2 重金屬含量分布特征






2.3 重金屬含量克里格插值分析






2.4 重金屬風(fēng)險(xiǎn)評(píng)價(jià)



3 結(jié)論