• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    基于GIS的農(nóng)業(yè)土壤重金屬風(fēng)險(xiǎn)評(píng)價(jià)研究*
    ——以南京市八卦洲為例

    2017-02-23 05:50:26宋金茜朱權(quán)姜小三趙海燕梁永紅羅永霞王強(qiáng)趙林麗
    土壤學(xué)報(bào) 2017年1期
    關(guān)鍵詞:農(nóng)業(yè)土壤八卦重金屬

    宋金茜朱 權(quán)姜小三?趙海燕梁永紅羅永霞王 強(qiáng)趙林麗

    (1 南京農(nóng)業(yè)大學(xué)資源與環(huán)境科學(xué)學(xué)院,南京 210095)

    (2 江蘇省耕地質(zhì)量與農(nóng)業(yè)環(huán)境保護(hù)站,南京 210036)

    基于GIS的農(nóng)業(yè)土壤重金屬風(fēng)險(xiǎn)評(píng)價(jià)研究*
    ——以南京市八卦洲為例

    宋金茜1朱 權(quán)1姜小三1?趙海燕1梁永紅2?羅永霞1王 強(qiáng)1趙林麗1

    (1 南京農(nóng)業(yè)大學(xué)資源與環(huán)境科學(xué)學(xué)院,南京 210095)

    (2 江蘇省耕地質(zhì)量與農(nóng)業(yè)環(huán)境保護(hù)站,南京 210036)

    農(nóng)產(chǎn)品產(chǎn)地環(huán)境安全受到廣泛關(guān)注,合理評(píng)價(jià)農(nóng)產(chǎn)品產(chǎn)地土壤重金屬的風(fēng)險(xiǎn)狀況,具有現(xiàn)實(shí)和長(zhǎng)遠(yuǎn)意義。以南京市八卦洲農(nóng)業(yè)土壤為例,選取對(duì)人體危害較大的As、Cd、Cr、Hg、Pb五種重金屬元素,對(duì)表層土壤進(jìn)行采樣分析,分別采用傳統(tǒng)標(biāo)準(zhǔn)方法、農(nóng)業(yè)部發(fā)布的《全國(guó)農(nóng)產(chǎn)品產(chǎn)地土壤重金屬安全評(píng)估技術(shù)規(guī)定》(產(chǎn)地評(píng)估標(biāo)準(zhǔn)方法)和改進(jìn)地累積綜合指數(shù)法對(duì)該地區(qū)農(nóng)業(yè)土壤重金屬風(fēng)險(xiǎn)狀況進(jìn)行評(píng)價(jià),運(yùn)用GIS地統(tǒng)計(jì)學(xué)方法分析重金屬風(fēng)險(xiǎn)等級(jí)的分布特征。結(jié)果表明:5種重金屬元素的平均含量分別為11.0、0.38、92.1、0.08和34.5 mg kg-1,除元素Hg以外,均超過(guò)南京市土壤環(huán)境背景值。產(chǎn)地評(píng)估標(biāo)準(zhǔn)方法表明八卦洲超過(guò)80%的耕地?zé)o重金屬風(fēng)險(xiǎn),風(fēng)險(xiǎn)區(qū)零散分布在八卦洲西南部、南部以及東部地區(qū),Cd和Pb是造成低風(fēng)險(xiǎn)的主要元素;傳統(tǒng)標(biāo)準(zhǔn)方法僅顯示了Cd元素的超標(biāo)特征;地累積指數(shù)法表明Cd、Cr、Pb三種元素受人為活動(dòng)影響較大,改進(jìn)地累積綜合指數(shù)法表明研究區(qū)總體生態(tài)狀況良好。傳統(tǒng)標(biāo)準(zhǔn)方法與產(chǎn)地評(píng)估標(biāo)準(zhǔn)方法能共同識(shí)別無(wú)風(fēng)險(xiǎn)的區(qū)域,而產(chǎn)地評(píng)估標(biāo)準(zhǔn)方法能更好地評(píng)價(jià)和劃分農(nóng)業(yè)土壤重金屬風(fēng)險(xiǎn)區(qū)域,建議與改進(jìn)地累積綜合指數(shù)法結(jié)合使用,共同評(píng)估該地區(qū)農(nóng)業(yè)土壤重金屬風(fēng)險(xiǎn)水平、空間分布以及生態(tài)風(fēng)險(xiǎn)特征。

    農(nóng)業(yè)土壤;重金屬;風(fēng)險(xiǎn)評(píng)價(jià);GIS

    土壤重金屬風(fēng)險(xiǎn)評(píng)價(jià)是土壤環(huán)境質(zhì)量評(píng)價(jià)的重要組成部分。土壤重金屬含量與母質(zhì)、土壤質(zhì)地、有機(jī)質(zhì)含量以及工業(yè)化、城市化和農(nóng)村集約化程度密切相關(guān)[1]。農(nóng)業(yè)土壤重金屬污染不僅會(huì)對(duì)土壤生態(tài)系統(tǒng)帶來(lái)不利影響,同時(shí)也會(huì)關(guān)系到食品安全和潛在的健康風(fēng)險(xiǎn)問(wèn)題[2]。開(kāi)展農(nóng)業(yè)土壤重金屬風(fēng)險(xiǎn)評(píng)價(jià)工作,可以了解農(nóng)業(yè)土壤重金屬的風(fēng)險(xiǎn)程度以及超標(biāo)元素類(lèi)型,為合理劃分農(nóng)產(chǎn)品產(chǎn)地環(huán)境土壤安全等級(jí)、有效利用土地資源、土壤修復(fù)工作等提供科學(xué)依據(jù)和指導(dǎo)建議。

    GIS技術(shù)應(yīng)用于土壤重金屬風(fēng)險(xiǎn)評(píng)價(jià)實(shí)現(xiàn)了從點(diǎn)狀數(shù)據(jù)到面狀信息的表達(dá)[3-4],其中地統(tǒng)計(jì)分析被廣泛應(yīng)用于土壤重金屬的空間變異研究[5-6]。在使用地統(tǒng)計(jì)學(xué)方法對(duì)評(píng)價(jià)指數(shù)進(jìn)行插值分析時(shí),針對(duì)不同的目標(biāo),對(duì)特定土地利用類(lèi)型的研究較為普遍,然而大多數(shù)研究未對(duì)相應(yīng)土地利用類(lèi)型的面積進(jìn)行提取分析,這可能會(huì)對(duì)插值評(píng)價(jià)結(jié)果中風(fēng)險(xiǎn)區(qū)域面積、范圍的判斷產(chǎn)生一定的影響。八卦洲是江北新區(qū)的一部分,以農(nóng)業(yè)和生態(tài)旅游業(yè)為主要發(fā)展特色,是重要的糧食蔬菜生產(chǎn)基地,對(duì)其農(nóng)業(yè)土壤進(jìn)行重金屬現(xiàn)狀風(fēng)險(xiǎn)評(píng)價(jià)十分必要。我國(guó)學(xué)者在參考不同的標(biāo)準(zhǔn)和方法的基礎(chǔ)上,對(duì)八卦洲不同區(qū)域土壤重金屬的濃度以及分布特征分別進(jìn)行了研究[7-11],雖取得了一些研究成果,但缺乏對(duì)評(píng)價(jià)方法以及評(píng)價(jià)標(biāo)準(zhǔn)的改進(jìn)和創(chuàng)新。而且由于采樣時(shí)期、區(qū)域、評(píng)價(jià)方法以及標(biāo)準(zhǔn)參數(shù)的差異,研究結(jié)果難以進(jìn)行比較,尚沒(méi)有一種公認(rèn)的可廣泛接受的模型或方法能全面評(píng)價(jià)農(nóng)業(yè)土壤重金屬的風(fēng)險(xiǎn)程度[12-13]。鑒于農(nóng)業(yè)部于2015年10月發(fā)布了《全國(guó)農(nóng)產(chǎn)品產(chǎn)地土壤重金屬安全評(píng)估技術(shù)規(guī)定》,其中修改了相關(guān)重金屬元素在農(nóng)產(chǎn)品產(chǎn)地土壤中的含量限值,并在單因子指數(shù)法評(píng)價(jià)的基礎(chǔ)上提出了最大單項(xiàng)指數(shù)法,用于采樣點(diǎn)多種重金屬元素的綜合評(píng)價(jià),因此,有必要采用該產(chǎn)地評(píng)估標(biāo)準(zhǔn)方法對(duì)八卦洲農(nóng)業(yè)土壤重金屬進(jìn)行風(fēng)險(xiǎn)評(píng)價(jià),并對(duì)原有評(píng)價(jià)方法進(jìn)行改進(jìn),與傳統(tǒng)標(biāo)準(zhǔn)方法評(píng)價(jià)結(jié)果進(jìn)行比較。

    本研究針對(duì)八卦洲農(nóng)業(yè)土壤中As、Cd、Cr、Hg、Pb五種重金屬元素,分別以南京市背景值、國(guó)家土壤環(huán)境二級(jí)標(biāo)準(zhǔn)值和農(nóng)業(yè)部發(fā)布的“農(nóng)產(chǎn)品產(chǎn)地土壤安全評(píng)估參比值”為基準(zhǔn),采用單因子與多因子結(jié)合的評(píng)價(jià)方法,對(duì)該地區(qū)農(nóng)業(yè)土壤重金屬進(jìn)行風(fēng)險(xiǎn)評(píng)價(jià),并通過(guò)ArcGIS10.2的地統(tǒng)計(jì)學(xué)方法對(duì)各綜合評(píng)價(jià)指數(shù)進(jìn)行插值分析,劃分風(fēng)險(xiǎn)等級(jí)分布區(qū),綜合評(píng)價(jià)八卦洲農(nóng)業(yè)土壤重金屬風(fēng)險(xiǎn)狀況,比較各評(píng)價(jià)方法與標(biāo)準(zhǔn)對(duì)評(píng)價(jià)結(jié)果的影響。

    1 材料與方法

    1.1 研究區(qū)概況

    八卦洲位于30°10′N(xiāo)、118°49′E,行政區(qū)域面積56 km2,處于南京市棲霞區(qū)西北部,與六合區(qū)、浦口區(qū)共同組成南京江北新區(qū),是長(zhǎng)江中僅次于崇明島和揚(yáng)中島的第三大島。八卦洲屬于亞熱帶季風(fēng)氣候區(qū),具有四季分明、溫暖濕潤(rùn)、雨量集中的特點(diǎn),降雨過(guò)程多發(fā)生在6—9月,其多年平均降雨量為1 000~1 100mm,多年平均氣溫為15.4℃,常年盛行東南風(fēng)和東北風(fēng)。洲上地勢(shì)低平,總體呈現(xiàn)西北略高、東南略低的格局。土壤類(lèi)型主要有水稻土和灰潮土,成土母質(zhì)主要為長(zhǎng)江沖積物,以種植特種蔬菜為主,是江蘇省最大的“野八鮮”生產(chǎn)基地,有“江中綠島”和“中國(guó)蘆蒿第一鄉(xiāng)”之稱(chēng)。洲兩岸匯聚了許多大型企業(yè),洲東南有南京新生圩港和南京新港開(kāi)發(fā)區(qū),洲北有南京鋼鐵集團(tuán)、南京熱電廠、華能南京電廠、南化集團(tuán)、揚(yáng)子石化等國(guó)家特大型企業(yè)。

    1.2 采樣點(diǎn)布設(shè)

    為了全面調(diào)查研究區(qū)土壤重金屬含量狀況和減少局部突變性信息的遺漏[14],采用ArcGIS10.2進(jìn)行網(wǎng)格布點(diǎn)。具體步驟:(1)創(chuàng)建大小為0.667 km2的網(wǎng)格(Create Fishnet):在八卦洲行政邊界圖上創(chuàng)建816.5m×816.5m的網(wǎng)格,并與該地土地利用圖合并(Union);(2)生成中心點(diǎn)(Feature To Point):將八卦洲的耕地選中為目標(biāo)區(qū)塊,在網(wǎng)格內(nèi)部生成中心點(diǎn);(3)篩選(Select by Attributes):每個(gè)網(wǎng)格內(nèi)僅保留一個(gè)中心點(diǎn),保證該點(diǎn)所代表的區(qū)塊面積大于0.334 km2;(4)校驗(yàn):使用該地區(qū)高精度的遙感圖逐點(diǎn)校正采樣點(diǎn)的位置,根據(jù)高清顯示的地物信息適當(dāng)刪除或移動(dòng)采樣點(diǎn)。最終生成66個(gè)采樣點(diǎn),對(duì)其轉(zhuǎn)換坐標(biāo)并添加經(jīng)緯度屬性,分幅制作采樣點(diǎn)位置圖以及相應(yīng)的經(jīng)緯度信息表。

    1.3 樣品采集與分析

    土壤樣品采集時(shí)間為2015年6月,將采樣點(diǎn)位置的相關(guān)圖表交由江蘇省耕地質(zhì)量與農(nóng)業(yè)環(huán)境保護(hù)站,由相關(guān)工作人員完成采樣,中心點(diǎn)由差分GPS定位并記錄,導(dǎo)入到ArcGIS10.2中根據(jù)經(jīng)緯度信息生成實(shí)際采樣點(diǎn)分布圖,見(jiàn)圖1。每個(gè)樣點(diǎn)采用梅花形采樣法采集5個(gè)以上子樣組成一個(gè)混合樣品,采樣深度為0~20cm,剔除土壤樣品與金屬采樣器接觸的部分后,放入聚乙烯塑料袋內(nèi)保存。將土壤樣品自然風(fēng)干,剔除樣品中植物根系、殘?jiān)翱梢?jiàn)侵入體,用木制工具碾碎,用瑪瑙研缽研磨過(guò)100目尼龍篩[15]。土樣的測(cè)定分析所用試劑、標(biāo)準(zhǔn)溶液配置以及儀器工作參數(shù)均遵循相關(guān)技術(shù)要求規(guī)范[16],分析數(shù)據(jù)由江蘇省耕地質(zhì)量與農(nóng)業(yè)環(huán)境保護(hù)站提供。As、Hg含量用氫化物發(fā)生原子熒光光譜法(HG-AFS)分析;Cd含量用石墨爐原子吸收光譜法(GF-AAS)分析;Cr、Pb含量采用電感耦合等離子體發(fā)射光譜法(ICP-AES)分析。樣品分析質(zhì)量水平由重復(fù)性檢驗(yàn)的方式監(jiān)控。

    圖1 八卦洲采樣點(diǎn)分布及地理位置圖Fig. 1 Distribution map of soil sampling points and geographical location map of the Baguazhou

    1.4 評(píng)價(jià)方法

    評(píng)價(jià)方法包括單因子評(píng)價(jià)和多因子綜合評(píng)價(jià),常用的單因子評(píng)價(jià)方法包括單因子指數(shù)法和地累積指數(shù)法等[17]。單因子指數(shù)法對(duì)土壤中的某一元素的風(fēng)險(xiǎn)程度進(jìn)行評(píng)價(jià),其計(jì)算公式為:

    式中,Pi為土壤中該元素的單因子指數(shù);Ci為元素的實(shí)測(cè)濃度;Si為評(píng)價(jià)標(biāo)準(zhǔn)(見(jiàn)表1)。農(nóng)業(yè)部發(fā)布的農(nóng)產(chǎn)品產(chǎn)地土壤重金屬安全評(píng)估參比值,在國(guó)家土壤環(huán)境質(zhì)量二級(jí)標(biāo)準(zhǔn)的基礎(chǔ)上進(jìn)行了修改和調(diào)整。由表1可以看出,農(nóng)業(yè)部標(biāo)準(zhǔn)對(duì)元素As、Pb、Cr按照水稻地、蔬菜地及其他農(nóng)產(chǎn)品產(chǎn)地進(jìn)行分類(lèi),有別于國(guó)家標(biāo)準(zhǔn)中按照水田和旱地進(jìn)行分類(lèi),并且對(duì)于元素Pb的限值進(jìn)行了較大的修改,其他元素的限值也進(jìn)行了不同幅度的調(diào)整。因此本文分別采用兩種評(píng)價(jià)標(biāo)準(zhǔn)進(jìn)行單因子指數(shù)的計(jì)算,以比較新老標(biāo)準(zhǔn)對(duì)評(píng)價(jià)結(jié)果的影響。

    在單因子指數(shù)的基礎(chǔ)上發(fā)展出來(lái)的多因子綜合評(píng)價(jià)方法包括最大單項(xiàng)指數(shù)法和內(nèi)梅羅綜合指數(shù)法等。其中最大單項(xiàng)指數(shù)法是農(nóng)業(yè)部環(huán)境檢測(cè)總站針對(duì)全國(guó)農(nóng)產(chǎn)品產(chǎn)地土壤規(guī)定的綜合評(píng)估方法,即按點(diǎn)位單項(xiàng)指數(shù)Pi的最大值Pimax來(lái)確定各個(gè)點(diǎn)位的安全等級(jí)。內(nèi)梅羅綜合指數(shù)法[17]同時(shí)考慮了污染物的平均值和最大值,其計(jì)算公式如下:

    地累積指數(shù)法由德國(guó)海德堡大學(xué)沉積物研究所的科學(xué)家Muller[19]在1969年提出,用于定量評(píng)價(jià)沉積物中的重金屬污染程度。該指數(shù)不僅反映了重金屬分布的自然變化特征,而且可以判別人為活動(dòng)對(duì)環(huán)境的影響[20]。其計(jì)算公式為:

    式中,Igeo為地累積指數(shù);Ci為元素的實(shí)測(cè)濃度;K是考慮各地巖石差異可能會(huì)引起背景值的變動(dòng)而取的系數(shù)(一般取值為K =1.5);Bi為該元素的地球化學(xué)背景值。

    在分別使用最大單項(xiàng)指數(shù)法和內(nèi)梅羅綜合指數(shù)法對(duì)研究區(qū)土壤重金屬進(jìn)行風(fēng)險(xiǎn)等級(jí)劃分時(shí),為了進(jìn)一步比較多因子評(píng)價(jià)方法的結(jié)果差異,本研究在地累積指數(shù)法的基礎(chǔ)上創(chuàng)建了改進(jìn)地累積綜合指數(shù)法。該方法考慮了重金屬的潛在生態(tài)毒性[21-22],對(duì)各元素的地累積指數(shù)進(jìn)行加權(quán)求和,以評(píng)價(jià)采樣區(qū)域由于人為因素對(duì)生態(tài)環(huán)境造成的綜合風(fēng)險(xiǎn)程度,計(jì)算公式為:

    表1 土壤重金屬元素相關(guān)標(biāo)準(zhǔn)[18]Table 1 Relevant standards for heavy metals in soil

    式中,I為改進(jìn)地累積綜合指數(shù);Igeo為各元素的地累積指數(shù);Ti為元素i的毒性因子權(quán)重;Tir為元素i的毒性系數(shù)[23](As 、Cd、 Cr、 Hg、 Pb的毒性系數(shù)分別為10、30、2、40、5);T為這五種重金屬毒性系數(shù)的和。表2為以上每種評(píng)價(jià)指數(shù)的分級(jí)標(biāo)準(zhǔn)。

    表2 評(píng)價(jià)指數(shù)分級(jí)標(biāo)準(zhǔn)Table 2 The classification standard of evaluation index

    1.5 統(tǒng)計(jì)分析與制圖

    利用Microsoft Excel 2013進(jìn)行數(shù)據(jù)的處理和統(tǒng)計(jì)分析,采用SPSS 20進(jìn)行重金屬元素的相關(guān)性分析。點(diǎn)位風(fēng)險(xiǎn)率只能揭示采樣點(diǎn)的重金屬風(fēng)險(xiǎn)特征,要對(duì)八卦洲農(nóng)業(yè)土壤進(jìn)行全面評(píng)價(jià)需借助軟件ArcGIS10.2中的地統(tǒng)計(jì)分析(Geostatistical Analyst)模塊進(jìn)行空間預(yù)測(cè)以及誤差驗(yàn)證。本研究側(cè)重對(duì)評(píng)價(jià)方法的比較以及風(fēng)險(xiǎn)等級(jí)的空間分布特征,因此分別對(duì)P、Pimax、I三種綜合評(píng)價(jià)指數(shù)進(jìn)行空間插值分析,選擇的空間預(yù)測(cè)模型為反距離權(quán)重(Inverse Distance Weighting,IDW)空間插值法。該方法是基于Tobler定理提出的一種空間確定性插值方法,其原理是通過(guò)計(jì)算未測(cè)量點(diǎn)附近各個(gè)點(diǎn)的測(cè)量值的加權(quán)平均來(lái)進(jìn)行插值,根據(jù)空間自相關(guān)性原理,在空間上越靠近的事物或現(xiàn)象越相似,則其在最近點(diǎn)處取得的權(quán)值最大[24]。反距離權(quán)重插值法僅考慮已知樣點(diǎn)與未知樣點(diǎn)的距離遠(yuǎn)近,而不考慮它們之間的空間分布及方位關(guān)系,適用于樣本數(shù)較少的空間插值。針對(duì)本研究目標(biāo),選擇反距離權(quán)重插值方法,插值底圖為八卦洲1∶2萬(wàn)耕地分布圖(2009年)[25],為了解插值的精度,對(duì)其進(jìn)行誤差驗(yàn)證。

    風(fēng)險(xiǎn)評(píng)價(jià)指數(shù)的空間分布制圖具體步驟為:①插值(IDW):創(chuàng)建訓(xùn)練子集和測(cè)試子集,樣本數(shù)分別為60和6,將訓(xùn)練子集進(jìn)行反距離權(quán)重插值;②驗(yàn)證(Validation):使用測(cè)試子集進(jìn)行驗(yàn)證,統(tǒng)計(jì)預(yù)測(cè)值與真實(shí)值之間的平均誤差(ME)和均方根誤差(RMSE),ME接近0,RMSE越小,則預(yù)測(cè)誤差越小,精度越高;③輸出(Export):將插值結(jié)果輸出為柵格單元大小為50 m×50 m的圖層,并用八卦洲的耕地分布圖對(duì)其進(jìn)行掩模提??;④重分類(lèi)(Reclassify):按照評(píng)價(jià)指數(shù)分級(jí)標(biāo)準(zhǔn)(表2)對(duì)風(fēng)險(xiǎn)程度進(jìn)行分級(jí),根據(jù)各風(fēng)險(xiǎn)等級(jí)對(duì)應(yīng)的像元個(gè)數(shù)計(jì)算該等級(jí)所占面積百分比,進(jìn)而制得八卦洲農(nóng)業(yè)土壤重金屬綜合風(fēng)險(xiǎn)等級(jí)分布圖。

    2 結(jié)果與討論

    2.1 土壤重金屬含量特征

    對(duì)土壤重金屬含量實(shí)測(cè)值進(jìn)行描述性統(tǒng)計(jì),結(jié)果見(jiàn)表3。5種重金屬元素的含量均值分別為11.0、0.38、92.1、0.08、34.5 mg kg-1,標(biāo)準(zhǔn)差分別為2.44、0.16、16.71、0.06、7.28 mg kg-1,除Hg以外,其余四種重金屬元素的含量均值和中值均超過(guò)南京地區(qū)土壤元素背景值。變異系數(shù)介于0.18~0.69之間,可以反映元素在土壤中的變異性和均勻性。變異程度大小為Hg>Cd>As>Pb>Cr,其中Hg和Cd的變異系數(shù)均達(dá)到0.4以上,分別屬于強(qiáng)變異和中等變異,空間變異相對(duì)顯著,易受人為活動(dòng)影響;As、Cr、Pb的變異系數(shù)均在0.2左右,屬于弱變異性,說(shuō)明這三種元素受外界影響比較小,空間變異相對(duì)不顯著,可能具有同源性。偏度是統(tǒng)計(jì)數(shù)據(jù)分布偏斜方向和程度的度量,偏度為零表示數(shù)值相對(duì)均勻地分布在平均值的兩側(cè),5種元素的偏度系數(shù)均為正,說(shuō)明數(shù)據(jù)為右偏離,可能是由于人為活動(dòng)等造成重金屬元素不斷富集的現(xiàn)象。峰度表征概率密度分布曲線在平均值處峰值高低的特征數(shù),峰度系數(shù)表明Cd、Hg的數(shù)據(jù)分布較為集中,As、Cr、Pb的數(shù)據(jù)分布較為分散。Pearson相關(guān)性分析表明(表4):As與Cr、Pb相互之間有極顯著相關(guān)性,進(jìn)一步證明了這三種元素的同源性。Cd、Cr、Hg和Pb相互之間相關(guān)性極顯著,容易有復(fù)合污染或同源性。

    表3 八卦洲土壤重金屬含量基本統(tǒng)計(jì)特征Table 3 Statistical characteristics of soil heavy metal contents in Baguazhou(n=66)

    表4 不同重金屬之間的相關(guān)性Table 4 Pearson correlation analysis of the 5 heavy metals in the soil(n=66)

    2.2 土壤重金屬風(fēng)險(xiǎn)評(píng)價(jià)指數(shù)

    將單因子指數(shù)與綜合評(píng)價(jià)指數(shù)分別進(jìn)行統(tǒng)計(jì)(表5),可以看出:根據(jù)國(guó)家土壤環(huán)境質(zhì)量二級(jí)標(biāo)準(zhǔn)得到單因子指數(shù)P1的平均值大小為Cd>As>Cr>Hg>Pb,均小于1,說(shuō)明土壤中五種重金屬元素的濃度均值均不超標(biāo),總體上符合農(nóng)業(yè)生產(chǎn)要求;內(nèi)梅羅綜合指數(shù)P的平均值為0.62,小于0.7,屬于安全生產(chǎn)級(jí)別。按照農(nóng)業(yè)部發(fā)布的《全國(guó)農(nóng)產(chǎn)品產(chǎn)地土壤重金屬安全評(píng)估技術(shù)規(guī)定》中的限值計(jì)算的單因子指數(shù)P2的平均值大小為Cd>As>Pb>Cr>Hg,均小于1,屬于無(wú)風(fēng)險(xiǎn)級(jí)別,與P1相比,該標(biāo)準(zhǔn)對(duì)元素Pb的規(guī)定更加嚴(yán)格;最大單項(xiàng)指數(shù)Pimax的均值為0.87,小于1,說(shuō)明采樣區(qū)整體土壤重金屬狀況良好。地累積指數(shù)Igeo的平均值大小為Cd>Cr>Pb>As>Hg,其中只有Cd和Cr的地累積指數(shù)均值大于0,分別為0.36和0.03,具有輕微風(fēng)險(xiǎn),改進(jìn)地累積綜合指數(shù)平均值為-0.52,小于0,說(shuō)明總體上人為活動(dòng)對(duì)重金屬的輸入未造成八卦洲農(nóng)業(yè)土壤生態(tài)環(huán)境的風(fēng)險(xiǎn)。

    根據(jù)評(píng)價(jià)指數(shù)分級(jí)標(biāo)準(zhǔn)(表2)計(jì)算有風(fēng)險(xiǎn)的樣本所占總體樣本的百分比,用點(diǎn)位風(fēng)險(xiǎn)率表示??傮w來(lái)看,大部分有風(fēng)險(xiǎn)的點(diǎn)位處于Ⅱ級(jí)低風(fēng)險(xiǎn)水平。其中,P1顯示分別有13個(gè)和1個(gè)樣本的Cd元素處于Ⅱ級(jí)和Ⅲ級(jí)水平,元素Cd的點(diǎn)位風(fēng)險(xiǎn)率約為21.2%;P2顯示元素Cd和Pb的點(diǎn)位風(fēng)險(xiǎn)率分別為19.7%和7.6%,除一個(gè)樣點(diǎn)外均處于Ⅱ級(jí)水平;Igeo顯示約80%、59%和32%的土壤樣點(diǎn)分別存在Cd、Cr和Pb元素輕微風(fēng)險(xiǎn)。內(nèi)梅羅綜合指數(shù)、最大單項(xiàng)指數(shù)以及改進(jìn)地累積綜合指數(shù)顯示分別有10.6%、23%和3%的土壤樣點(diǎn)處于風(fēng)險(xiǎn)水平。其中有極個(gè)別樣點(diǎn)達(dá)到Ⅲ級(jí)或Ⅳ級(jí),即中等—高風(fēng)險(xiǎn),需加密布點(diǎn)監(jiān)測(cè)。

    表5 土壤重金屬評(píng)價(jià)指數(shù)結(jié)果Table 5 Statistics of evaluation indices of soil heavy metals

    2.3 土壤重金屬風(fēng)險(xiǎn)指數(shù)分布特征

    將三種綜合評(píng)價(jià)方法得到的風(fēng)險(xiǎn)指數(shù)P、I、Pimax分別進(jìn)行空間插值(圖2)??梢钥闯觯趦?nèi)梅羅綜合指數(shù)P得到該區(qū)域有三個(gè)風(fēng)險(xiǎn)等級(jí):無(wú)風(fēng)險(xiǎn)的區(qū)域占八卦洲耕地的70.2%,符合安全生產(chǎn)級(jí)別;低風(fēng)險(xiǎn)區(qū)所占面積為27.5%,屬于安全生產(chǎn)的警戒線,這些地區(qū)的土壤重金屬風(fēng)險(xiǎn)需引起關(guān)注,加強(qiáng)監(jiān)測(cè)工作;高風(fēng)險(xiǎn)區(qū)所占面積僅為2.3%,主要分布在八卦洲南部和東部,分布較零散,主要受元素Cd的分布影響,這可能與八卦洲東南有較多大型工廠企業(yè)的分布有關(guān),需加緊污染源的排查和管控,開(kāi)展土壤修復(fù)工作。基于改進(jìn)地累積綜合指數(shù)I同樣得出三個(gè)風(fēng)險(xiǎn)等級(jí):八卦洲約98.7%的耕地處于Ⅰ級(jí),即考慮到各元素的生態(tài)毒性差異,大部分耕地的重金屬累積程度未對(duì)土壤生態(tài)環(huán)境造成危害;西南角極少地區(qū)有較為突出的重金屬累積以及生態(tài)風(fēng)險(xiǎn),可能與元素Cd、Pb以及Cr的復(fù)合污染有關(guān)?;谧畲髥雾?xiàng)指數(shù)Pimax得到該地區(qū)有四個(gè)風(fēng)險(xiǎn)等級(jí):Ⅰ級(jí)和Ⅱ級(jí)的區(qū)域面積分別占八卦洲耕地的81.5%和17.9%,其中Ⅱ級(jí),即低風(fēng)險(xiǎn)區(qū)主要分布在八卦洲南部和東部,Ⅲ級(jí)和Ⅳ級(jí)的區(qū)域面積共占0.6%,位于西南角地區(qū),主要是Cd和Pb兩種元素的超標(biāo)導(dǎo)致了最大單項(xiàng)指數(shù)的風(fēng)險(xiǎn)分布特征。對(duì)插值精度進(jìn)行驗(yàn)證得到,P、I、Pimax三項(xiàng)風(fēng)險(xiǎn)指數(shù)的預(yù)測(cè)值平均誤差ME分別為-0.05、-0.11、-0.02,即反距離權(quán)重插值法對(duì)三種評(píng)價(jià)指數(shù)的預(yù)測(cè)結(jié)果略微低于實(shí)測(cè)水平,均方根誤差RMSE分別為0.24、0.26、0.21,結(jié)果均在0.25左右,在允許的誤差范圍之內(nèi)。

    綜上所述,三種綜合評(píng)價(jià)方法均能識(shí)別八卦洲北部以及西部土壤安全區(qū)域,對(duì)于西南部重金屬風(fēng)險(xiǎn)較高的點(diǎn)位同樣能識(shí)別,后期需對(duì)該點(diǎn)加密監(jiān)測(cè),進(jìn)一步勘察。內(nèi)梅羅綜合指數(shù)法與最大單項(xiàng)指數(shù)法在采用不同的標(biāo)準(zhǔn)基礎(chǔ)上得出的風(fēng)險(xiǎn)分布規(guī)律大體相似,前者結(jié)合單因子指數(shù)的最大值與平均值進(jìn)行綜合風(fēng)險(xiǎn)評(píng)價(jià),平均值會(huì)削弱風(fēng)險(xiǎn)指數(shù),對(duì)于某些整體狀況良好,然而存在個(gè)別元素含量較高的區(qū)域判定可能有一定的影響,而后者通過(guò)最大值的方法可以很好地揭示有任何元素超標(biāo)的區(qū)域分布特征,對(duì)于土壤風(fēng)險(xiǎn)的排查以及重金屬元素的超標(biāo)統(tǒng)計(jì)具有重要意義。改進(jìn)地累積綜合指數(shù)法從人為活動(dòng)以及生態(tài)危害的角度評(píng)價(jià)土壤重金屬的風(fēng)險(xiǎn)特征,顯示了該地區(qū)的總體生態(tài)環(huán)境處于良好的狀態(tài)。

    圖2 內(nèi)梅羅綜合指數(shù)分級(jí)圖(a)、改進(jìn)地累積綜合指數(shù)分級(jí)圖(b)與最大單項(xiàng)指數(shù)分級(jí)圖(c)Fig. 2 Grading map of Nemerow comprehensive pollution index(a),Improved geo-accumulation comprehensive index(b)and Maximum single factor index(c)

    3 結(jié) 論

    八卦洲農(nóng)業(yè)土壤中元素Cd、Cr、Pb均有不同程度的超標(biāo),其中Cd元素的風(fēng)險(xiǎn)狀況應(yīng)引起重視。風(fēng)險(xiǎn)區(qū)分布以西南地區(qū)最嚴(yán)重,其次為南部和東部的部分地區(qū),這可能與八卦洲岸邊分布的大中型工廠企業(yè)有關(guān),具體風(fēng)險(xiǎn)來(lái)源有待進(jìn)一步研究。針對(duì)與人類(lèi)生活飲食緊密相關(guān)的糧食蔬菜產(chǎn)地,農(nóng)業(yè)部產(chǎn)地評(píng)估標(biāo)準(zhǔn)對(duì)土壤中部分元素的監(jiān)控更加嚴(yán)格,對(duì)農(nóng)業(yè)土壤普查和重金屬風(fēng)險(xiǎn)評(píng)價(jià)工作更加適宜。最大單項(xiàng)指數(shù)法簡(jiǎn)單直觀,針對(duì)土壤調(diào)查中超標(biāo)區(qū)域的判定以及風(fēng)險(xiǎn)級(jí)別的劃分有較好的適用性,改進(jìn)地累積綜合指數(shù)法可以反應(yīng)人為活動(dòng)對(duì)生態(tài)環(huán)境造成的影響,兩種方法可結(jié)合使用,全面評(píng)價(jià)農(nóng)產(chǎn)品產(chǎn)地土壤重金屬風(fēng)險(xiǎn)與生態(tài)環(huán)境安全,對(duì)農(nóng)業(yè)土壤的安全生產(chǎn)等級(jí)調(diào)查以及土壤修復(fù)工作具有十分重要的意義。

    致 謝感謝江蘇省耕地質(zhì)量與農(nóng)業(yè)環(huán)境保護(hù)站提供數(shù)據(jù)支撐,感謝國(guó)家地球系統(tǒng)數(shù)據(jù)共享平臺(tái)(http://nnu.geodata.cn)提供地理數(shù)據(jù)。

    [1]鐘曉蘭,周生路,李江濤,等. 長(zhǎng)江三角洲地區(qū)土壤重金屬污染的空間變異特征—以江蘇省太倉(cāng)市為例. 土壤學(xué)報(bào),2007,44(1):33—40

    Zhong X L,Zhou S L,Li J T,et al. Spatial variability of soil heavy metals contamination in the Yangtze River Delta—A case study of Taicang City in Jiangsu Province (In Chinese). Acta Pedologica Sinica,2007,44 (1):33—40

    [2]肖雪,趙南京,袁靜,等. 淮河流域某鎮(zhèn)農(nóng)業(yè)土壤重金屬含量特征及污染評(píng)價(jià). 光譜學(xué)與光譜分析,2014 (7):1785—1788

    Xiao X,Zhao N J,Yuan J,et al. Analysis of heavy metals distribution characteristics and pollution assessment in agricultural region soils of Huaihe basin (In Chinese). Spectroscopy and Spectral Analysis,2014(7):1785—1788

    [3]Burgos P,Madejón E,Pérez-de-Mora A,et al. Spatial variability of the chemical characteristics of a traceelement-contaminated soil before and after remediation. Geoderma,2006,130:157—175

    [4]Thomas R A,Bradley S. Spatial interpolation and image-integrative geostatistical prediction of mosquito vectors for arboviral surveillance. Geocarto International,2008,23(4):311—325

    [5]夏敏,趙炳梓,張佳寶. 基于GIS的黃淮海平原典型潮土區(qū)土壤重金屬積累研究. 土壤學(xué)報(bào),2013,50 (4):684—692

    Xia M,Zhao B Z,Zhang J B. GIS-based research on soil heavy metal accumulation in a fluvo-aquic soil area typical of the Huang-Huai-Hai Plain(In Chinese). Acta Pedologica Sinica,2013,50(4):684—692

    [6]林世滔,李琳,盧志紅,等. 基于GIS的江西省耕地土壤重金屬污染評(píng)價(jià)研究. 江西農(nóng)業(yè)大學(xué)學(xué)報(bào),2014 (5):1167—1172

    Lin S T,Li L,Lu Z H,et al. Evaluation of heavy metal pollution of farmland soil in Jiangxi Province based on GIS(In Chinese). Acta Agriculturae Universitatis Jiangxiensis,2014(5):1167—1172

    [7]吳昀昭. 南京城郊農(nóng)業(yè)土壤重金屬污染的遙感地球化學(xué)基礎(chǔ)研究. 南京:南京大學(xué),2005

    Wu Y Z. Heavy metal pollution in suburban soils of the Nanjing area—A feasibility study of remotesensing geochemistry(In Chinese). Nanjing:Nanjing University,2005

    [8]陳峰,尹春芹,蔣新,等. 基于GIS的南京市典型蔬菜基地土壤重金屬污染現(xiàn)狀與評(píng)價(jià). 中國(guó)環(huán)境監(jiān)測(cè),2008,24(2):40—45

    Chen F,Yin C Q,Jiang X,et al. Evaluation of heavy metal pollution in representative vegetable soils in Nanjing based on GIS(In Chinese). Environmental Monitoring in China,2008,24(2):40—45

    [9]周生路,廖富強(qiáng),吳紹華,等. 基于分等樣地的江蘇典型區(qū)農(nóng)用地土壤重金屬污染研究. 農(nóng)業(yè)工程學(xué)報(bào),2008,24(5):78—83

    Zhou S L,Liao F Q,Wu S H,et al. Farmland soil heavy metal pollution in typical areas of Jiangsu Province based on classification sample plots(In Chinese). Transactions of the Chinese Society of Agricultural Engineering,2008,24(5):78—83

    [10]馬宏瑞,張茜,季俊峰,等. 長(zhǎng)江南京段近岸沉積物中重金屬富集特征與形態(tài)分析. 生態(tài)環(huán)境學(xué)報(bào),2009,18(6):2061—2065

    Ma H R,Zhang Q,Ji J F,et al. Enrichment and speciation analysis of heavy metals in the sediments of Yangtze River(Nanjing section)(In Chinese). Ecology and Environment,2009,18(6):2061—2065

    [11]袁潤(rùn)杰,于高偉,邱曉蕾,等. 區(qū)域土壤重金屬空間差異及對(duì)蔬菜品質(zhì)安全影響的分析—以南京八卦洲、江心洲為例. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào),2015,34(8):1498—1507

    Yuan R J,Yu G W,Qiu X L,et al. Regionally spatial variation of soil heavy metals and their influences on vegetable quality:A case study of Baguazhou and Jiangxinzhou,Nanjing,China(In Chinese). Journal of Agro-Environment Science,2015,34(8):1498—1507

    [12]馬成玲,周健民,王火焰,等.農(nóng)田土壤重金屬污染評(píng)價(jià)方法研究——以長(zhǎng)江三角洲典型縣級(jí)市常熟市為例.生態(tài)與農(nóng)村環(huán)境學(xué)報(bào),2006,22(1):48—53

    Ma C L,Zhou J M,Wang H Y,et al. Methods for assessment of heavy metal pollution in cropland soils——A case study of Changshu(In Chinese). Journal of Ecology and Rural Environment,2006,22 (1):48—53

    [13]王顯煒,徐友寧,楊敏,等.國(guó)內(nèi)外礦山土壤重金屬污染風(fēng)險(xiǎn)評(píng)價(jià)方法綜述.中國(guó)礦業(yè),2009,18(10):54—56

    Wang X W,Xu Y N,Yang M,et al. Review on risk assessment methods for soil heavy metal contamination in mines at home and abroad(In Chinese). China Mining Magazine,2009,18(10):54—56

    [14]吳春發(fā),吳嘉平,駱永明,等. 冶煉廠周邊土壤重金屬污染范圍的界定與不確定性分析. 土壤學(xué)報(bào),2009,46(6):1006—1012

    Wu C F,Wu J P,Luo Y M,et al. Identification and uncertainty analysis of scope of heavy metal pollution in the vicinity of secondary copper smelteries(In Chinese). Acta Pedologica Sinica,2009,46(6):1006—1012

    [15]鮑士旦.土壤農(nóng)化分析.第3版.北京:中國(guó)農(nóng)業(yè)出版社,2000:20—22

    Bao S D. Soil and agricultural chemistry analysis(In Chinese). 3rd ed. Beijing:China Agriculture Press,2000:20—22

    [16]中國(guó)環(huán)境監(jiān)測(cè)總站.土壤元素的近代分析方法.北京:中國(guó)環(huán)境科學(xué)出版社,1992:39—262

    China National Environmental Monitoring Centre. The modern analysis method of soil elements(In Chinese). Beijing:China Environmental Science Press,1992:39—262

    [17]范拴喜,甘卓亭,李美娟,等. 土壤重金屬污染評(píng)價(jià)方法進(jìn)展. 中國(guó)農(nóng)學(xué)通報(bào),2010,26(17):310—315

    Fan S X,Gan Z T,Li M J,et al. Progress of assessment methods of heavy metal pollution in soil(In Chinese). Chinese Agricultural Science Bulletin,2010,26(17):310—315

    [18]國(guó)家環(huán)境保護(hù)總局. 土壤環(huán)境質(zhì)量標(biāo)準(zhǔn)(GB15618-1995). 北京:中國(guó)環(huán)境科學(xué)出版社,1995

    State Environmental Protection Administration of China. Environmental quality standard for soil(GB15618-1995)(In Chinese).Beijing:China Environmental Science Press,1995

    [19]Muller G. Index of geoaccumulation in sediment of the Rhine River. Geojournal,1969,2:108—118

    [20]滕彥國(guó),庹先國(guó),倪師軍,等. 應(yīng)用地質(zhì)累積指數(shù)評(píng)價(jià)沉積物中重金屬污染:選擇地球化學(xué)背景的影響. 環(huán)境科學(xué)與技術(shù),2002,25(2):7—9

    Teng Y G,Tuo X G,Ni S J,et al. Applying geoaccumulation index to assess heavy metal pollution in sediment:Influence of different geochemical background(In Chinese). Environmental Science & Technology,2002,25(2):7—9

    [21]Hakanson L. An ecological risk index for aquatic pollution control-A sedimentological approach(In Chinese). Water Research,1980,14(8):975—1001

    [22]郭平,謝忠雷,李軍,等. 長(zhǎng)春市土壤重金屬污染特征及其潛在生態(tài)風(fēng)險(xiǎn)評(píng)價(jià). 地理科學(xué),2005,25(1):108—112

    Guo P,Xie Z L,Li J,et al. Specificity of heavy metal pollution and the ecological hazard in urban soils of Changchun City(In Chinese). Scientia Geographica Sinica,2005,25(1):108—112

    [23]徐爭(zhēng)啟,倪師軍,庹先國(guó),等. 潛在生態(tài)危害指數(shù)法評(píng)價(jià)中重金屬毒性系數(shù)計(jì)算. 環(huán)境科學(xué)與技術(shù),2008,31(2):112—115

    Xu Z Q,Ni S J,Tuo X G,et al. Calculation of heavy metals' toxicity coefficient in the evaluation of potential ecological risk index(In Chinese). Environmental Science & Technology,2008,31(2):112—115

    [24]朱選,劉素霞.地理信息系統(tǒng)原理與技術(shù).上海:華東師范大學(xué)出版社,2006:171—173

    Zhu X,Liu S X. Principles and techniques of geographic information system(In Chinese). Shanghai:East China Normal University Press,2006:171—173

    [25]國(guó)家基礎(chǔ)條件平臺(tái):地球系統(tǒng)科學(xué)數(shù)據(jù)共享平臺(tái)—長(zhǎng)江三角洲科學(xué)數(shù)據(jù)中心. http://nnu.geodata.cn

    The National R&D Infrastructure and Facility Development Program of China:Data Sharing Infrastructure of Earth System Science .http://nnu.geodata.cn

    [26]夏增祿,李森照,李延芳,等. 土壤元素背景值及其研究方法. 北京:氣象出版社,1987:295

    Xia Z L,Li S Z,Li Y F,et al. Background values and the research method of soil elements(In Chinese). Beijing:China Meteorological Press,1987:295

    GIS-Based Heavy Metals Risk Assessment of Agricultural Soils —A Case Study of Baguazhou,Nanjing

    SONG Jinxi1ZHU Quan1JIANG Xiaosan1?ZHAO Haiyan1LIANG Yonghong2?LUO Yongxia1WANG Qiang1ZHAO Linli1
    (1 College of Resources and Environmental Sciences,Nanjing Agricultural University,Nanjing 210095,China)
    (2 Jiangsu Station of Farmland Quality and Agricultural Environment Protection,Nanjing 210036,China)

    【Objective】Safety of agricultural producing areas has become an extensive concern. Heavy metal pollution of agricultural soils not only brings about adverse effects on the soil ecosystem,but also is closely related to food safety and potential health hazard. How to assessment the risk of heavy metals pollution of agricultural soils is of some practical and far-reaching significance. Therefore,it is a very urgent task to investigate the status of heavy metals in agricultural soils and assess the risk with a variety of assessment methods. The findings in this study may help decision-makers select appropriate assessment methods for rational assessment and zoning of the agricultural producing area and remedy contaminated soils scientifically. 【Method】Taking the agricultural soil in Baguazhou Town of Nanjing as an example,this study focused on the five heavy metal elements,i.e. As,Cd,Cr,Hg and Pb,that are relatively more harmful than the others to the human body and analyzed the soil samples collected from the surface soil layer(0~20cm)in 2015. In October 2015,the Chinese Ministry of Agriculture promulgated the“National technical regulations for safety assessment of heavy metals in agricultural producing areas”(Standard Methods for Producing Areas),which has modified the criteria of some relevant elements and specified the use of a new method called the maximum single factor index method for assessing integrated risk of heavy metals in agricultural soils. It is necessary to adopt the specific standard and method to do the assessment. The study used the Traditional Standard Method(T-SM)and the Standard Method for Producing Areas(PA-SM)to assess the risk of heavy metals in the study area and compared the two methods for differences. Furthermore,an Improved Geoaccumulation Comprehensive Index Method(IG-M)was created to assess the ecological risk of heavy metals in the study area. With the aid of GIS geostatistics,spatial distribution of different heavy metals risk levels was characterized. The inverse distance weighting interpolation method was used in this study.【Result】Results show that the average concentration of As,Cd,Cr,Hg and Pb was11.0,0.38,92.1,0.08 and 34.5 mg kg-1,respectively. All of them,except for Hg,were beyond their background values in Nanjing. The elements of Cd,Cr and Pb exceeded the criteria for soil quality to a varying extent,especially Cd,which should be paid more attention to. The use of PA-SM found that over 80% of arable land in Baguazhou was riskfree and Cd and Pb were the main factors that posed a low risk. Risk areas scattered across the southwest,south and east parts of Baguazhou,which may be associated with the distribution of large and medium factories and plants,while the use of T-SM displayed that only Cd exceeded the standard. The use of PASM revealed that the risks of Cd,Cr and Pb were relatively high due to human activities,while the use of the IG-M disclosed that the study area as a whole was good in ecology.【Conclusion】Generally speaking,the agricultural environment of Baguazhou was good,but the risk area in the south needs intensive monitoring andadoption of pollution control measures. T-SM and PA-SM can work together to identify risk-free region,but PA-SM can better evaluate and divide risk areas. Therefore,it is recommended to use PA-SM in conjunction with IG-M to assess risk levels,spatial distribution and ecological risk of the heavy metals in the agricultural soil of the region.

    Agricultural soil;Heavy metals;Risk assessment;GIS

    X82

    A

    10.11766/trxb201603100033

    (責(zé)任編輯:檀滿(mǎn)枝)

    * 國(guó)家水體污染控制與治理科技重大專(zhuān)項(xiàng)(2012ZX07101-009)和江蘇省農(nóng)產(chǎn)品產(chǎn)地安全智慧管控系統(tǒng)項(xiàng)目(JSZC-G2015-039)共同資助 Supported by the National water pollution control and management technology major projects(2012ZX07101-009)and Intelligent management system for safety of agricultural producing area of Jiangsu Province,China(JSZC-G2015-039)

    ? 通訊作者 Corresponding author,E-mail:gis@njau.edu.cn;1597152683@qq.com

    宋金茜(1994—),女,河南新鄉(xiāng)人,碩士研究生,主要從事土壤資源環(huán)境遙感與信息系統(tǒng)研究。E-mail:412883774@qq.com

    2016-03-10;

    2016-07-30;優(yōu)先數(shù)字出版日期(www.cnki.net):2016-09-19

    猜你喜歡
    農(nóng)業(yè)土壤八卦重金屬
    伏羲創(chuàng)八卦
    幼兒100(2021年7期)2021-04-10 06:32:14
    重金屬對(duì)膨潤(rùn)土膨脹性的影響
    測(cè)定不同產(chǎn)地寬筋藤中5種重金屬
    中成藥(2018年8期)2018-08-29 01:28:16
    固原市部分農(nóng)業(yè)土壤無(wú)機(jī)磷形態(tài)及轉(zhuǎn)化
    綠色科技(2017年22期)2017-12-25 15:55:47
    農(nóng)業(yè)土壤中鎘污染現(xiàn)狀及污染途徑分析
    我國(guó)土壤污染現(xiàn)狀及污染防治制度概況
    ICP-AES、ICP-MS測(cè)定水中重金屬的對(duì)比研究
    世相八卦
    山東青年(2016年1期)2016-02-28 14:26:36
    世相八卦
    山東青年(2016年2期)2016-02-28 14:26:20
    再生水回灌中DOM對(duì)重金屬遷移與保留問(wèn)題研究
    国产成人精品久久久久久| 久久国产亚洲av麻豆专区| 伦精品一区二区三区| 91精品国产国语对白视频| 国产精品成人在线| 欧美一级a爱片免费观看看| 国产精品久久久久久精品电影小说| 久久精品久久精品一区二区三区| 日本av手机在线免费观看| 夫妻午夜视频| 日本vs欧美在线观看视频 | 在线观看美女被高潮喷水网站| 又大又黄又爽视频免费| 亚洲熟女精品中文字幕| 久久精品国产自在天天线| 亚洲真实伦在线观看| 国产成人精品福利久久| 成人二区视频| 国产精品久久久久成人av| 亚洲精品久久午夜乱码| 日韩大片免费观看网站| 日韩中字成人| 在线看a的网站| 精品人妻熟女毛片av久久网站| tube8黄色片| 熟妇人妻不卡中文字幕| 狂野欧美激情性bbbbbb| 欧美国产精品一级二级三级 | 最近手机中文字幕大全| 一边亲一边摸免费视频| 深夜a级毛片| 欧美亚洲 丝袜 人妻 在线| 久久鲁丝午夜福利片| 国产高清三级在线| 一级黄片播放器| 国国产精品蜜臀av免费| 国产精品福利在线免费观看| 黄色一级大片看看| 久久精品国产鲁丝片午夜精品| 18禁在线无遮挡免费观看视频| 国产一区亚洲一区在线观看| 亚洲av国产av综合av卡| 人人妻人人澡人人爽人人夜夜| 欧美97在线视频| 99热这里只有精品一区| 国产在视频线精品| 丰满迷人的少妇在线观看| 亚洲av.av天堂| 久久精品久久久久久噜噜老黄| 插逼视频在线观看| 成年人免费黄色播放视频 | 亚洲性久久影院| av专区在线播放| 自线自在国产av| 亚洲精品中文字幕在线视频 | 国产日韩欧美在线精品| 国产乱来视频区| 午夜老司机福利剧场| 国产综合精华液| av免费在线看不卡| 日产精品乱码卡一卡2卡三| 黑丝袜美女国产一区| 少妇被粗大的猛进出69影院 | a级一级毛片免费在线观看| 99热这里只有是精品在线观看| 日韩免费高清中文字幕av| 久久影院123| 久久久久视频综合| 婷婷色麻豆天堂久久| 欧美日韩综合久久久久久| 久久精品熟女亚洲av麻豆精品| 晚上一个人看的免费电影| 天堂俺去俺来也www色官网| 婷婷色综合www| 美女内射精品一级片tv| 久久99热这里只频精品6学生| 欧美日韩视频高清一区二区三区二| 国产精品一区二区在线不卡| 久久久精品94久久精品| 一边亲一边摸免费视频| 97超视频在线观看视频| 搡女人真爽免费视频火全软件| 色吧在线观看| 男女国产视频网站| 久久影院123| 有码 亚洲区| 免费不卡的大黄色大毛片视频在线观看| 七月丁香在线播放| 精品久久久久久久久亚洲| 街头女战士在线观看网站| 亚洲精品亚洲一区二区| 国产成人精品福利久久| 伦精品一区二区三区| 青春草视频在线免费观看| 精品视频人人做人人爽| 免费黄频网站在线观看国产| 国产精品国产三级国产专区5o| 免费黄色在线免费观看| 简卡轻食公司| 日韩亚洲欧美综合| 大片免费播放器 马上看| 久久青草综合色| 99九九线精品视频在线观看视频| 一本久久精品| 精品国产一区二区久久| 全区人妻精品视频| www.av在线官网国产| 在线精品无人区一区二区三| 成人午夜精彩视频在线观看| 97超碰精品成人国产| 国产亚洲av片在线观看秒播厂| 国产精品伦人一区二区| 乱码一卡2卡4卡精品| 99热这里只有精品一区| 久久久久久久国产电影| 免费看日本二区| 久久久久久人妻| 色视频www国产| 麻豆成人午夜福利视频| 一级毛片我不卡| 国模一区二区三区四区视频| 在线观看美女被高潮喷水网站| a级毛片免费高清观看在线播放| 大片电影免费在线观看免费| 亚洲内射少妇av| 国产国拍精品亚洲av在线观看| 亚洲精品乱久久久久久| 亚洲美女视频黄频| 爱豆传媒免费全集在线观看| 丰满人妻一区二区三区视频av| 极品人妻少妇av视频| 精品国产乱码久久久久久小说| 国产欧美另类精品又又久久亚洲欧美| 国产伦精品一区二区三区视频9| 69精品国产乱码久久久| 日本午夜av视频| 日本vs欧美在线观看视频 | 欧美日韩av久久| 中文乱码字字幕精品一区二区三区| 美女cb高潮喷水在线观看| 欧美日韩av久久| 99热6这里只有精品| 99热全是精品| 在线观看人妻少妇| 久久婷婷青草| 久久久久久久久久人人人人人人| 人体艺术视频欧美日本| 日韩一本色道免费dvd| 精品国产一区二区久久| 97在线视频观看| 能在线免费看毛片的网站| 韩国av在线不卡| 99热这里只有精品一区| av视频免费观看在线观看| 亚洲性久久影院| 人人妻人人澡人人爽人人夜夜| 中文天堂在线官网| 午夜av观看不卡| 十分钟在线观看高清视频www | 亚洲精品国产色婷婷电影| 嫩草影院入口| 精品国产国语对白av| 国产亚洲精品久久久com| 高清欧美精品videossex| a 毛片基地| 亚洲av综合色区一区| 免费看日本二区| 亚洲人成网站在线播| 97在线人人人人妻| 日韩欧美一区视频在线观看 | 国产精品免费大片| 97在线视频观看| 国产精品久久久久久av不卡| 黄色毛片三级朝国网站 | 人妻一区二区av| 亚洲国产日韩一区二区| 成人亚洲欧美一区二区av| 中文天堂在线官网| 久久久久精品久久久久真实原创| 免费久久久久久久精品成人欧美视频 | av福利片在线| 国产成人精品婷婷| 久久午夜综合久久蜜桃| 国产日韩一区二区三区精品不卡 | 久久99精品国语久久久| 亚洲怡红院男人天堂| 中国国产av一级| 能在线免费看毛片的网站| 女性被躁到高潮视频| 免费看av在线观看网站| 九九爱精品视频在线观看| 久久久久久久亚洲中文字幕| 亚洲丝袜综合中文字幕| 国产精品人妻久久久久久| 看十八女毛片水多多多| 国产成人精品无人区| 一区二区三区精品91| 日韩中文字幕视频在线看片| 国产精品人妻久久久影院| 日日摸夜夜添夜夜爱| 午夜福利影视在线免费观看| 亚洲四区av| 黄色怎么调成土黄色| av国产久精品久网站免费入址| 亚洲国产色片| 国产在视频线精品| 欧美日韩av久久| 麻豆精品久久久久久蜜桃| 精品亚洲乱码少妇综合久久| 乱码一卡2卡4卡精品| 欧美老熟妇乱子伦牲交| 人人妻人人澡人人爽人人夜夜| 成人毛片60女人毛片免费| 如日韩欧美国产精品一区二区三区 | 国产精品嫩草影院av在线观看| 亚洲人与动物交配视频| 欧美日韩亚洲高清精品| 乱码一卡2卡4卡精品| 日本黄色日本黄色录像| 成人毛片60女人毛片免费| 国产精品无大码| 高清av免费在线| 国产精品女同一区二区软件| 日韩免费高清中文字幕av| 久久精品夜色国产| 免费黄网站久久成人精品| 丁香六月天网| 99久久精品一区二区三区| 看免费成人av毛片| 久久久久久久久久久免费av| 婷婷色麻豆天堂久久| 男女边吃奶边做爰视频| 9色porny在线观看| 能在线免费看毛片的网站| 国产日韩欧美视频二区| 日韩一区二区视频免费看| 国产欧美亚洲国产| 亚洲精品,欧美精品| www.色视频.com| 9色porny在线观看| 国产黄色视频一区二区在线观看| 伊人久久国产一区二区| 亚洲欧美日韩东京热| 男的添女的下面高潮视频| 欧美日韩精品成人综合77777| 我的女老师完整版在线观看| 伦理电影免费视频| 99久国产av精品国产电影| 欧美3d第一页| 日韩视频在线欧美| 国产精品不卡视频一区二区| 在线看a的网站| 国产熟女午夜一区二区三区 | 免费高清在线观看视频在线观看| 国产成人精品婷婷| 丰满少妇做爰视频| 亚洲av不卡在线观看| 日韩电影二区| av天堂中文字幕网| 亚洲精品日韩av片在线观看| 亚洲自偷自拍三级| 亚洲va在线va天堂va国产| 99热这里只有精品一区| 人人妻人人爽人人添夜夜欢视频 | www.色视频.com| 久久午夜综合久久蜜桃| 亚洲精品国产av成人精品| 国产欧美日韩一区二区三区在线 | 成人美女网站在线观看视频| 国产一区二区三区综合在线观看 | 国产精品久久久久成人av| 啦啦啦啦在线视频资源| 亚洲欧美日韩东京热| 国产真实伦视频高清在线观看| 少妇的逼好多水| 欧美日韩国产mv在线观看视频| 老熟女久久久| 久久久久久久久久成人| 中文欧美无线码| 日韩av免费高清视频| 一级毛片我不卡| 欧美 日韩 精品 国产| 成人综合一区亚洲| 成人免费观看视频高清| 婷婷色麻豆天堂久久| 免费不卡的大黄色大毛片视频在线观看| 91精品伊人久久大香线蕉| 美女视频免费永久观看网站| 欧美变态另类bdsm刘玥| 国内少妇人妻偷人精品xxx网站| 亚洲一级一片aⅴ在线观看| 国产黄片视频在线免费观看| 水蜜桃什么品种好| 美女中出高潮动态图| 午夜福利在线观看免费完整高清在| 国产伦精品一区二区三区视频9| 精品少妇久久久久久888优播| 少妇人妻精品综合一区二区| 国产一区有黄有色的免费视频| 性高湖久久久久久久久免费观看| 欧美日韩在线观看h| 国产男女内射视频| 午夜福利视频精品| 国产欧美日韩一区二区三区在线 | 亚洲图色成人| 色94色欧美一区二区| 各种免费的搞黄视频| 麻豆成人午夜福利视频| 尾随美女入室| 丝袜在线中文字幕| 亚洲精品一二三| 国产成人精品久久久久久| 国产极品粉嫩免费观看在线 | 日产精品乱码卡一卡2卡三| 亚洲精品国产成人久久av| a 毛片基地| 中文资源天堂在线| 男女啪啪激烈高潮av片| 国产精品国产三级国产av玫瑰| 午夜福利视频精品| 如何舔出高潮| 精品酒店卫生间| 国产欧美日韩综合在线一区二区 | 在线 av 中文字幕| 涩涩av久久男人的天堂| 亚洲性久久影院| 大话2 男鬼变身卡| 一区二区av电影网| 观看美女的网站| 在现免费观看毛片| 亚洲不卡免费看| 久久精品熟女亚洲av麻豆精品| 亚洲国产日韩一区二区| 亚洲国产精品一区二区三区在线| 欧美日韩在线观看h| 男女边吃奶边做爰视频| 一区二区三区乱码不卡18| 国产一级毛片在线| av黄色大香蕉| 99国产精品免费福利视频| 两个人的视频大全免费| 26uuu在线亚洲综合色| 欧美性感艳星| 午夜激情福利司机影院| 日本欧美视频一区| 亚洲经典国产精华液单| 深夜a级毛片| 午夜久久久在线观看| 99久久综合免费| 中文字幕免费在线视频6| 人妻 亚洲 视频| 久久人妻熟女aⅴ| 少妇人妻 视频| 视频中文字幕在线观看| 国产成人aa在线观看| 午夜影院在线不卡| 亚洲经典国产精华液单| 国产精品国产三级专区第一集| 91久久精品国产一区二区三区| 欧美日韩国产mv在线观看视频| 婷婷色av中文字幕| 久久亚洲国产成人精品v| 亚洲激情五月婷婷啪啪| 黑人高潮一二区| 在线观看免费高清a一片| 热re99久久精品国产66热6| 日韩欧美一区视频在线观看 | 亚洲精品国产成人久久av| 中文字幕精品免费在线观看视频 | 欧美精品一区二区大全| 亚洲欧美成人综合另类久久久| 我要看日韩黄色一级片| 国产精品久久久久久精品电影小说| 国产深夜福利视频在线观看| 日韩,欧美,国产一区二区三区| 国产精品国产三级专区第一集| 亚洲av福利一区| 热re99久久国产66热| 一级爰片在线观看| av.在线天堂| 色哟哟·www| 日本黄色片子视频| 麻豆成人av视频| 久久精品久久久久久噜噜老黄| 久久青草综合色| 乱码一卡2卡4卡精品| 亚洲精品视频女| 国产高清有码在线观看视频| 久久久久久久久久久丰满| 亚洲电影在线观看av| 免费大片18禁| av在线老鸭窝| 久久免费观看电影| 日韩av免费高清视频| a级毛色黄片| 国产高清不卡午夜福利| 日本猛色少妇xxxxx猛交久久| 国产精品一区www在线观看| 国产在线视频一区二区| 国产精品成人在线| 国产成人免费无遮挡视频| 赤兔流量卡办理| 少妇人妻 视频| 亚洲欧美清纯卡通| 亚洲无线观看免费| 一级毛片 在线播放| 国产精品久久久久久精品电影小说| 国产成人一区二区在线| 一区二区av电影网| 男人和女人高潮做爰伦理| 国产探花极品一区二区| 亚洲av男天堂| 极品教师在线视频| freevideosex欧美| 能在线免费看毛片的网站| 97超碰精品成人国产| 最黄视频免费看| 免费黄色在线免费观看| 日本黄大片高清| 国产成人a∨麻豆精品| 久久这里有精品视频免费| 在线观看免费日韩欧美大片 | 一级毛片 在线播放| 精品人妻熟女av久视频| 日本黄色日本黄色录像| 久久久久视频综合| 亚洲精品成人av观看孕妇| 日韩熟女老妇一区二区性免费视频| 黑人高潮一二区| 亚洲av不卡在线观看| 在线观看av片永久免费下载| av卡一久久| 日韩一区二区视频免费看| 如日韩欧美国产精品一区二区三区 | 国产免费又黄又爽又色| 久久热精品热| 免费人妻精品一区二区三区视频| 日韩强制内射视频| 成人二区视频| 最近中文字幕2019免费版| 色吧在线观看| 精品人妻熟女av久视频| 在线精品无人区一区二区三| 国产伦精品一区二区三区视频9| 下体分泌物呈黄色| 欧美日本中文国产一区发布| 国产男人的电影天堂91| 三级国产精品片| 免费观看在线日韩| 校园人妻丝袜中文字幕| 91成人精品电影| 国产精品福利在线免费观看| 精品久久久精品久久久| 尾随美女入室| 一本色道久久久久久精品综合| 99热6这里只有精品| 国产黄片美女视频| 男女国产视频网站| 狠狠精品人妻久久久久久综合| 亚洲av二区三区四区| 欧美另类一区| 极品少妇高潮喷水抽搐| 日本av免费视频播放| 国产淫片久久久久久久久| 久久久久久久久久久免费av| 国产视频首页在线观看| 日本91视频免费播放| 一级,二级,三级黄色视频| 丰满少妇做爰视频| 亚洲在久久综合| 一级毛片电影观看| 精品亚洲成a人片在线观看| 乱码一卡2卡4卡精品| 又爽又黄a免费视频| 自拍偷自拍亚洲精品老妇| 人妻人人澡人人爽人人| 欧美激情国产日韩精品一区| 欧美少妇被猛烈插入视频| 精品人妻一区二区三区麻豆| 国产精品国产三级专区第一集| 最后的刺客免费高清国语| 国产成人一区二区在线| 国产欧美另类精品又又久久亚洲欧美| 男人狂女人下面高潮的视频| 日本黄大片高清| av在线播放精品| 99久久精品一区二区三区| 成人国产麻豆网| 少妇熟女欧美另类| 国产熟女欧美一区二区| 偷拍熟女少妇极品色| 少妇精品久久久久久久| 亚洲天堂av无毛| 精品一区二区三区视频在线| 伦理电影大哥的女人| 中国国产av一级| 久久婷婷青草| 亚洲国产欧美在线一区| 91成人精品电影| 99热国产这里只有精品6| 欧美xxxx性猛交bbbb| 插逼视频在线观看| 国产成人免费无遮挡视频| 尾随美女入室| 日本猛色少妇xxxxx猛交久久| 乱码一卡2卡4卡精品| 精品久久久久久久久亚洲| 麻豆乱淫一区二区| 国内少妇人妻偷人精品xxx网站| 性色avwww在线观看| 国产黄片美女视频| 久久久久久久精品精品| 赤兔流量卡办理| av又黄又爽大尺度在线免费看| 亚洲精品乱码久久久v下载方式| 国产永久视频网站| 日韩亚洲欧美综合| 欧美日韩一区二区视频在线观看视频在线| 久久精品久久精品一区二区三区| 夜夜爽夜夜爽视频| 大码成人一级视频| 日日啪夜夜爽| 午夜av观看不卡| 日韩亚洲欧美综合| .国产精品久久| 亚洲成人手机| 欧美日本中文国产一区发布| 成人影院久久| 少妇的逼水好多| 中文精品一卡2卡3卡4更新| 日日啪夜夜爽| 国产一区二区在线观看日韩| 99热这里只有是精品在线观看| 男人舔奶头视频| 亚洲欧美一区二区三区国产| 全区人妻精品视频| 亚洲va在线va天堂va国产| 中文字幕人妻丝袜制服| 伦理电影免费视频| 亚洲欧美成人精品一区二区| 男女边吃奶边做爰视频| 欧美97在线视频| 色吧在线观看| 久久精品国产a三级三级三级| 免费观看a级毛片全部| 在线免费观看不下载黄p国产| 日本-黄色视频高清免费观看| 最后的刺客免费高清国语| 国产高清不卡午夜福利| 91精品国产国语对白视频| 九草在线视频观看| 亚洲,一卡二卡三卡| 国产精品99久久久久久久久| 久久ye,这里只有精品| 80岁老熟妇乱子伦牲交| 国产亚洲午夜精品一区二区久久| 国产欧美日韩精品一区二区| 男女边吃奶边做爰视频| 免费人成在线观看视频色| 又粗又硬又长又爽又黄的视频| 极品少妇高潮喷水抽搐| 26uuu在线亚洲综合色| 晚上一个人看的免费电影| 久久99蜜桃精品久久| 水蜜桃什么品种好| av天堂中文字幕网| 一本久久精品| 国产免费福利视频在线观看| 国精品久久久久久国模美| 久久精品熟女亚洲av麻豆精品| 欧美日韩在线观看h| 国产无遮挡羞羞视频在线观看| 如日韩欧美国产精品一区二区三区 | 高清视频免费观看一区二区| 欧美性感艳星| 熟女av电影| 亚洲情色 制服丝袜| 日韩欧美一区视频在线观看 | 黑人巨大精品欧美一区二区蜜桃 | 亚洲精品中文字幕在线视频 | av又黄又爽大尺度在线免费看| 亚洲欧美一区二区三区国产| 麻豆成人av视频| 少妇人妻久久综合中文| 热re99久久精品国产66热6| 亚洲精品乱码久久久久久按摩| 亚洲欧美日韩卡通动漫| 国产成人午夜福利电影在线观看| 国产精品三级大全| 秋霞伦理黄片| 多毛熟女@视频| 日韩成人伦理影院| 久久综合国产亚洲精品| 国产伦在线观看视频一区| 日日爽夜夜爽网站| 国产黄片美女视频| 国产免费一区二区三区四区乱码| 欧美最新免费一区二区三区| 能在线免费看毛片的网站| 国产精品久久久久久久电影| 日韩免费高清中文字幕av| 国产av精品麻豆| 99久国产av精品国产电影| 美女大奶头黄色视频| 韩国av在线不卡| 国产午夜精品久久久久久一区二区三区| 国产无遮挡羞羞视频在线观看| 午夜av观看不卡| 老熟女久久久| 国产欧美另类精品又又久久亚洲欧美| 在线观看一区二区三区激情| 国产男女内射视频| 97超碰精品成人国产| 啦啦啦中文免费视频观看日本| 99久久精品一区二区三区|