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

    基于小流域尺度的土壤重金屬分布與土地利用相關(guān)性研究
    ——以廈門市坂頭水庫(kù)流域?yàn)槔?/h1>
    2015-01-18 07:39:45李清良高進(jìn)波徐秋芳
    生態(tài)學(xué)報(bào) 2015年16期
    關(guān)鍵詞:生態(tài)研究

    李清良,吳 倩,高進(jìn)波,馬 軍,徐秋芳,俞 慎,*

    1 浙江農(nóng)林大學(xué),杭州 311300

    2 中國(guó)科學(xué)院城市環(huán)境研究所城市環(huán)境與健康重點(diǎn)實(shí)驗(yàn)室,廈門 361021

    3 中國(guó)科學(xué)院大學(xué),北京 100049

    基于小流域尺度的土壤重金屬分布與土地利用相關(guān)性研究
    ——以廈門市坂頭水庫(kù)流域?yàn)槔?/p>

    李清良1,2,吳 倩2,3,高進(jìn)波2,3,馬 軍2,徐秋芳1,*,俞 慎2,*

    1 浙江農(nóng)林大學(xué),杭州 311300

    2 中國(guó)科學(xué)院城市環(huán)境研究所城市環(huán)境與健康重點(diǎn)實(shí)驗(yàn)室,廈門 361021

    3 中國(guó)科學(xué)院大學(xué),北京 100049

    人類活動(dòng)改變流域集水區(qū)土地利用方式同時(shí)提高了土壤重金屬水平,并影響地表水重金屬負(fù)荷。以具有明確地理邊界的流域作為研究單元可能實(shí)現(xiàn)計(jì)量區(qū)別不同人類活動(dòng)或因其引起的土地利用方式對(duì)土壤重金屬空間分布及其環(huán)境風(fēng)險(xiǎn)的貢獻(xiàn)。但不同人類活動(dòng)及其改變的土地利用方式對(duì)土壤和地表水重金屬負(fù)荷的相對(duì)貢獻(xiàn)研究缺乏可操作性方法論。以我國(guó)典型的山塘水庫(kù)小流域—福建省廈門市坂頭水庫(kù)流域?yàn)槔?,結(jié)合遙感影像解譯、野外土壤重金屬水平調(diào)查和GIS地統(tǒng)計(jì)方法,針對(duì)小流域尺度土壤重金屬空間分布及其環(huán)境風(fēng)險(xiǎn)與土地利用方式相關(guān)性研究進(jìn)行方法論的探索。坂頭水庫(kù)流域覆蓋205km2,土地利用類型有林地、農(nóng)業(yè)用地、城鎮(zhèn)用地、綠地和水面等5大類。按照不同土地利用類型覆蓋面積為權(quán)重,以統(tǒng)計(jì)最小樣本數(shù)為準(zhǔn)則,在該流域內(nèi)共采集150個(gè)表層土壤樣品,針對(duì)6種具有人為源特征重金屬(Cu、Zn、Pb、Cr、Ni、Cd)土壤總量及富集水平進(jìn)行分析,以潛在生態(tài)風(fēng)險(xiǎn)指數(shù)評(píng)價(jià)流域土壤重金屬風(fēng)險(xiǎn),并結(jié)合土地利用類型分析了其空間特征。結(jié)果表明:與林地和綠地相比,城鎮(zhèn)用地和農(nóng)業(yè)用地表層土壤重重金屬水平顯著較高;土壤重金屬富集水平及其潛在生態(tài)風(fēng)險(xiǎn)依次為:城鎮(zhèn)用地>農(nóng)業(yè)用地>綠地>林地;城鎮(zhèn)化加重了流域土壤重金屬Cu、Zn和Cd的污染風(fēng)險(xiǎn);土壤Pb富集水平與土地利用方式無(wú)關(guān),可能為大氣沉降來(lái)源;所有土地利用方式土壤Cr和Ni總量低于區(qū)域土壤背景值,但城鎮(zhèn)用地的富集指數(shù)顯著高于林地?;谕恋乩梅绞降姆淳嚯x加權(quán)空間插值(LU-IDW)也清晰地揭示了土壤重金屬富集及其潛在生態(tài)風(fēng)險(xiǎn)與人類活動(dòng)密切相關(guān),尤其是城鎮(zhèn)化。以野外調(diào)查結(jié)合遙感衛(wèi)星圖像解譯和基于土地利用方式的空間插值的研究方法有效地揭示了研究流域土壤重金屬的空間分布特征及其潛在風(fēng)險(xiǎn)評(píng)價(jià)空間模式,為小流域尺度的環(huán)境質(zhì)量演變研究提供了方法論和案例。

    重金屬;表層土壤;土地利用類型;城鎮(zhèn)化過程;潛在生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)

    隨著全球人口的增長(zhǎng),土地利用方式由自然土地向農(nóng)業(yè)和城鎮(zhèn)人為利用方式改變,已經(jīng)成為全球變化重要組成之一[1]。人類面臨著如何平衡人類福祉用地需求和自然生態(tài)系統(tǒng)服務(wù)價(jià)值可持續(xù)性間的重大矛盾[2]。城鎮(zhèn)化在現(xiàn)代經(jīng)濟(jì)發(fā)展和人類福祉提升中起著重要的作用,但同時(shí)也改變了自然生態(tài)系統(tǒng)及其元素生物地球化學(xué)循環(huán),在不同程度上影響了人類居住環(huán)境[3- 4]。

    重金屬污染具有隱蔽性、不可降解性、生物富集性和弱移動(dòng)性等特點(diǎn)[5]。交通運(yùn)輸、工業(yè)生產(chǎn)、市政建設(shè)等人類活動(dòng)都向環(huán)境排放重金屬,經(jīng)大氣沉降或污水灌溉使土壤重金屬負(fù)荷大幅度增加,直接影響農(nóng)產(chǎn)品和飲用水質(zhì)量安全及人類健康[6- 9]。

    已報(bào)道針對(duì)土壤重金屬污染的研究包括來(lái)源解析[10- 11]、累積特征[12- 13]、空間分布特征[14- 15]、土壤環(huán)境質(zhì)量狀況與食品安全性[16]、土壤重金屬污染評(píng)價(jià)[17- 18]等各方面,主要集中于區(qū)域資源開發(fā)和農(nóng)田污灌等導(dǎo)致土壤重金屬污染及其帶來(lái)的環(huán)境和人體健康風(fēng)險(xiǎn),但對(duì)不同人類活動(dòng)及其引起不同土地利用變化對(duì)土壤重金屬負(fù)荷貢獻(xiàn)的研究一直停留在定性表征水平。以具有明顯地理邊界的流域?yàn)檠芯繂卧?,其相?duì)封閉性可實(shí)現(xiàn)宏微觀相結(jié)合的環(huán)境計(jì)量學(xué)研究,即從流域尺度實(shí)現(xiàn)土壤重金屬空間分布與土地利用方式間的計(jì)量研究。但相關(guān)研究較為少見,尤其針對(duì)城鎮(zhèn)化對(duì)流域土壤和地表水的重金屬負(fù)荷貢獻(xiàn)研究鮮見報(bào)道。以流域?yàn)閱卧暮晡⒂^相結(jié)合研究方法論可在空間尺度上實(shí)現(xiàn)重金屬分布及其環(huán)境風(fēng)險(xiǎn)評(píng)價(jià)插值繪圖,為流域數(shù)字管理及其污染物排放減量提供科學(xué)依據(jù)及對(duì)策。

    本文以福建省廈門市坂頭水庫(kù)流域?yàn)閷?duì)象,在空間水平研究不同土地利用方式下研究表層土壤中具人為源特征重金屬元素(Cu、Zn、Pb、Ni、Cd、Cr)負(fù)荷及其空間分布特征,并應(yīng)用土壤重金屬富集指數(shù)評(píng)價(jià)其潛在生態(tài)風(fēng)險(xiǎn),以期為我國(guó)城鎮(zhèn)化流域土壤環(huán)境質(zhì)量評(píng)價(jià)及土地利用規(guī)劃提供科學(xué)依據(jù)和指導(dǎo)。本文還應(yīng)用本研究團(tuán)隊(duì)改良的基于土地利用方式的融合空間插值方法進(jìn)行重金屬元素分布及其潛在生態(tài)風(fēng)險(xiǎn)的空間分析。

    1 研究方法

    1.1 研究小流域概況

    研究小流域位于福建省廈門市西北部坂頭水庫(kù)流域,地理位置為24°39′12″ N、118°0′58″ E,流域面積約為205 km2,屬亞熱帶季風(fēng)性氣候。流域土壤以紅壤為主。根據(jù)流域氣候特征可以分為濕季和干季,濕季為每年的4月至10月,干季為每年的11月至3月,全年平均降雨量約1200 mm。根據(jù)2010年DEM數(shù)據(jù),該流域可分為5個(gè)次級(jí)流域,即水源地水庫(kù)次級(jí)流域(Ⅰ)、農(nóng)業(yè)城鎮(zhèn)復(fù)合次級(jí)流域(Ⅱ)、農(nóng)林城鎮(zhèn)復(fù)合次級(jí)流域(Ⅲ)、城鎮(zhèn)開發(fā)區(qū)次級(jí)流域(Ⅳ)和海灣次級(jí)流域(Ⅴ);經(jīng)解譯2010年LandSat ETM+影像和2012年實(shí)地校正,該流域土地利用類型主要為林地、農(nóng)業(yè)用地、城鎮(zhèn)用地、綠地和水面(圖1)。

    圖1 研究流域地理位置及土地利用方式和采樣點(diǎn)分布圖 Fig.1 Location of the studied watershed and map of land use type and sampling points

    1.2 表層土壤樣品采集與預(yù)處理

    土壤采樣方案:以次級(jí)流域?yàn)閱卧鶕?jù)土地利用類型及面積設(shè)置土壤采樣點(diǎn)數(shù)量,為統(tǒng)計(jì)需要每種土地利用類型至少保證3個(gè)表層土壤采樣點(diǎn)。因此,在次級(jí)流域I和II內(nèi)使用500m×500 m的采樣格網(wǎng),次級(jí)流域IV和V內(nèi)使用400m×400m的采樣格網(wǎng),次級(jí)流域(Ⅲ)內(nèi)使用200m×200m的采樣格網(wǎng),共設(shè)150個(gè)采樣點(diǎn),其中:林地25個(gè),農(nóng)業(yè)用地49個(gè),城鎮(zhèn)用地63個(gè),綠地13個(gè)(圖1)。2013年3月6日—4月7日根據(jù)預(yù)設(shè)布點(diǎn)進(jìn)行野外采樣,記錄GPS坐標(biāo);實(shí)際采樣過程中對(duì)于落在河流、道路等地采樣點(diǎn)調(diào)整到附近樣地。每個(gè)采樣點(diǎn)均按照10m×10m正方形4個(gè)頂點(diǎn)和中心點(diǎn)分別采樣,各取表層(0—10cm)土壤約1kg,現(xiàn)場(chǎng)均勻混合后用四分法從中選取1kg樣品作為代表該采樣點(diǎn)的混合樣品。土壤經(jīng)登記編號(hào)后,用聚氯乙烯塑料袋封裝,帶回實(shí)驗(yàn)室內(nèi)自然風(fēng)干,剔除植物殘?bào)w和石塊,磨碎、過100目(0.15mm)尼龍網(wǎng)篩,混勻、保存?zhèn)溆?。為防止樣品污染,在采樣和保存過程中避免與重金屬器皿直接接觸。

    1.3 土壤重金屬總量分析

    處理后土壤樣品重金屬總量分析采用強(qiáng)酸消解和ICP-OES分析測(cè)定[19]。每批處理都用長(zhǎng)江流域水系沉積物標(biāo)準(zhǔn)物(GSD- 13)做回收率實(shí)驗(yàn),并以混合酸溶液做空白實(shí)驗(yàn)。重金屬Cu、Zn、Pb、Cr、Cd和Ni的平均回收率分別為(104±8)%、(81±5)%、(97±4)%、(102±9)%、(98±6)%、(105±9)%。

    實(shí)驗(yàn)中所有器皿均在10%硝酸溶液中酸化24h以上,所用的試劑均為優(yōu)級(jí)純。

    1.4 潛在生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)方法

    流域表層土壤重金屬潛在生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)采用瑞典學(xué)者Hakanson 1980年提出的潛在生態(tài)風(fēng)險(xiǎn)指數(shù)法(PERI)[20- 21],其計(jì)算公式如下:

    (1)

    (2)

    (3)

    1.5 空間數(shù)據(jù)差值和統(tǒng)計(jì)

    空間數(shù)據(jù)插值:應(yīng)用本研究團(tuán)隊(duì)改良的基于土地利用類型融合空間反距離加權(quán)差值法(LU-IDW),即在ArcGIS(9.3)平臺(tái)上對(duì)4種陸地土地利用類型斑塊內(nèi)分別應(yīng)用反距離加權(quán)插值成圖,利用流域土地利用矢量圖融合得到流域表層土壤重金屬總量和風(fēng)險(xiǎn)指數(shù)連續(xù)空間分布圖。該方法是對(duì)不同土地利用類型進(jìn)行分別插值,然后融合成圖,區(qū)別于應(yīng)用模型一次插值成圖的傳統(tǒng)方法。

    表層土壤重金屬總量、富集指數(shù)、潛在生態(tài)風(fēng)險(xiǎn)因子及其指數(shù)在不同土地利用方式間的差異進(jìn)行了方差分析(ANOVA,P<0.05,SPSS(18.0)系統(tǒng))。

    2 結(jié)果

    2.1 流域不同土地利用方式表層土壤重金屬總量及其空間分布特征

    在調(diào)查的6種重金屬中,土地利用方式僅對(duì)表層土壤Cr和Ni總量有顯著的影響,而對(duì)土壤Cd、Cu、Pb和Zn總量影響不顯著(圖2,P<0.05)。農(nóng)業(yè)用地和城鎮(zhèn)用地表層土壤Cr總量相當(dāng),平均值分別為14.7和14.6mg/kg,均顯著大于林地表層土壤(8.5mg/kg)。表層土壤Ni總量則以城鎮(zhèn)用地最高(平均值0.82mg/kg),但僅顯著高于林地表層土壤(0.37mg/kg)。城鎮(zhèn)用地表層土壤Cd、Cu、Pb和Zn總量雖然由于數(shù)值變化區(qū)間太大,與其他土地利用方式差異統(tǒng)計(jì)不顯著,但平均值均高于其他用地方式(圖2)。

    圖2 研究流域不同土地利用類型下表層土壤Cu、Ni 、Zn、Pb、Cr和Cd總量Fig.2 Total contents of Cu, Ni, Zn, Pb, Cr and Cd in topsoil of different land use types in the studied watershed

    空間插值結(jié)果與統(tǒng)計(jì)結(jié)果一致,即城鎮(zhèn)用地和農(nóng)業(yè)用地使流域表層土壤重金屬富集,6種重金屬元素雖呈現(xiàn)出不同的空間分布特征,但趨勢(shì)一致(圖3)。

    2.2 流域表層土壤重金屬潛在生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)

    圖3 研究流域表層土壤重金屬總量(Cu、Ni、Zn、Cr、Cd、Pb)空間分布特征Fig.3 Spatial variation of total content of heavy metals (Cu, Ni, Zn, Cr, Cd and Pb) in topsoil of the studied watershed

    表1 研究流域不同土地利用類型表層土壤重金屬富集指數(shù)統(tǒng)計(jì)值

    3 討論

    本研究小流域表層土壤重金屬總量空間分布特征表明,人類活動(dòng)引起的土地利用方式改變提高了土壤重金屬水平。與林地相比較,隨著人類活動(dòng)強(qiáng)度逐漸遞增,綠地、農(nóng)業(yè)用地和城鎮(zhèn)用地方式表層土壤重金屬總量相應(yīng)提高(圖2)。Biplob Das等[26]亦證實(shí)了土地利用類型決定土壤重金屬水平,在其研究的淡水流域,土壤重金屬總量隨湖泊<林地<農(nóng)業(yè)用地<城鎮(zhèn)用地<礦區(qū)用地等土地利用方式逐步提高。Li等[19]進(jìn)一步發(fā)現(xiàn)公園表層土壤重金屬總量隨郊區(qū)到城區(qū)的城鎮(zhèn)化水平梯度遞增。但是,Van Landeghem等[25]在美國(guó)德克薩斯州流域的研究則認(rèn)為土地利用方式不影響流域表層土壤重金屬水平。

    表2 研究流域不同土地利用類型表層土壤重金屬潛在生態(tài)風(fēng)險(xiǎn)因子和潛在生態(tài)風(fēng)險(xiǎn)指數(shù)(RI)

    *平均值±標(biāo)準(zhǔn)差; 不同字母表示差異水平P<0.05,NS表示不顯著

    圖4 研究流域表層土壤重金屬潛在生態(tài)風(fēng)險(xiǎn)指數(shù)空間分布 Fig.4 Spatial distribution of the potential ecological risk index of heavy metals in topsoil of the studied watershed

    流域生態(tài)系統(tǒng)中重金屬濃度的提高直接或間接地威脅著生物受體和生態(tài)系統(tǒng)健康[31]。表層土壤中重金屬元素除了其母質(zhì)風(fēng)化釋放一部分外,主要來(lái)自于各種人類活動(dòng)[32]。比如,Cu、Zn和Pb是交通污染源特征重金屬元素[33],汽車輪胎磨損和尾氣排放均含有Zn和Cd的釋放[34]。農(nóng)業(yè)施用磷肥、含磷復(fù)合肥以及以城鎮(zhèn)垃圾、污泥為原料有機(jī)肥中均混雜有Cu、Zn、Pb、Cr、Cd和Ni等重金屬[35],因此,這幾種重金屬被稱為具有人為排放特征的污染重金屬元素,即人為源重金屬元素。本研究流域的表層土壤重金屬含量空間分布圖印證了這些人為輸入途徑,重金屬Cu、Zn、Pb、Cr、Cd和Ni的高值區(qū)均分布在流域南部城鎮(zhèn)用地、中部部分農(nóng)業(yè)用地以及主要交通線(圖3)。本研究結(jié)果表明研究流域尚處于城鎮(zhèn)化初期,流域內(nèi)廢棄物處理設(shè)施尚不完善,使得鄰近區(qū)域土壤重金屬富集水平相對(duì)較高,將進(jìn)一步導(dǎo)致流域水體重金屬水平提高。該研究流域已作為城鎮(zhèn)化環(huán)境與生態(tài)響應(yīng)的長(zhǎng)期定位研究站,將進(jìn)一步在時(shí)間尺度進(jìn)行驗(yàn)證。同時(shí),本研究以野外調(diào)查結(jié)合遙感衛(wèi)星圖像解譯和基于土地利用方式的空間插值的研究方法有效地揭示了研究流域土壤重金屬的空間分布特征及其潛在生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)空間模式。該方法論存在進(jìn)一步改善的需要,但其多學(xué)科結(jié)合將有助于流域尺度的環(huán)境科學(xué)研究。

    4 結(jié)論

    本文研究結(jié)果表明:

    (1)流域土地利用方式影響著表層土壤重金屬總量、富集水平和潛在生態(tài)風(fēng)險(xiǎn),人類活動(dòng)強(qiáng)度較大的城鎮(zhèn)用地和農(nóng)業(yè)用地方式對(duì)表層土壤重金屬負(fù)荷影響最大;

    (2)研究流域土壤重金屬富集水平相對(duì)較低,以Cu、Zn、Pb、Cd富集為主,Cr和Ni總量低于區(qū)域背景值,但其土地利用方式差異顯著。流域土壤重金屬的潛在生態(tài)風(fēng)險(xiǎn)為輕度水平,除城鎮(zhèn)用地和農(nóng)業(yè)用地的個(gè)別重金屬外,6種重金屬總量、富集水平和潛在生態(tài)風(fēng)險(xiǎn)因子均隨城鎮(zhèn)用地 > 農(nóng)業(yè)用地 > 綠地 > 林地遞減;

    (3)以野外調(diào)查結(jié)合遙感衛(wèi)星圖像解譯和基于土地利用方式的空間插值的研究方法有效地揭示了小流域尺度土壤重金屬的空間分布特征及其潛在生態(tài)風(fēng)險(xiǎn)的空間模式。

    [1] Martinuzzi S, Januchowski-Hartley S R, Pracheil B M, McIntyre P B, Plantinga A J, Lewis D J, Radeloff V C. Threats and opportunities for freshwater conservation under future land use change scenarios in the United States. Global Change Biology, 2014, 20(1): 113- 124.

    [2] Foley J A, DeFries R, Asner G P, Barford C, Bonan G, Carpenter S R, Chapin F S, Coe M T, Daily G C, Gibbs H K, Helkowski J H, Holloway T, Howard E A, Kucharik C J, Monfreda C, Patz J A, Prentice I C, Ramankutty N, Snyder P K. Global consequences of land use. Science, 2005, 309(5734): 570- 574.

    [3] Kaye J P, Groffman P M, Grimm N B, Baker L A, Pouyat R V. A distinct urban biogeochemistry?. Trends in Ecology & Evolution, 2006, 21(4): 192- 199.

    [4] Grimm N B, Faeth S H, Golubiewski N E, Redman C L, Wu J G, Bai X M, Briggs J M. Global change and the ecology of cities. Science, 2008, 319(5864): 756- 760.

    [5] 唐發(fā)靜, 祖艷群. 土壤重金屬空間變異的研究方法. 云南農(nóng)業(yè)大學(xué)學(xué)報(bào), 2008, 23(4): 558- 561.

    [6] Li X D, Lee S L, Wong S C, Shi W Z, Thornton I. The study of metal contamination in urban soils of Hong Kong using a GIS-based approach. Environmental Pollution, 2004, 129(1): 113- 124.

    [7] Zhang C S. Using multivariate analyses and GIS to identify pollutants and their spatial patterns in urban soils in Galway, Ireland. Environmental Pollution, 2006, 142(3): 501- 511.

    [8] 石寧寧, 丁艷鋒, 趙秀峰, 王強(qiáng)盛. 某農(nóng)藥工業(yè)園區(qū)周邊土壤重金屬含量與風(fēng)險(xiǎn)評(píng)價(jià). 應(yīng)用生態(tài)學(xué)報(bào), 2010, 21(7): 1835- 1843.

    [9] Mcgratha D, Zhang C S, Cartona O T. Geostatistical analyses and hazard assessment on soil lead in Silvermines area, Ireland. Environmental Pollution, 2004, 127(2): 239- 248.

    [10] Qishlaqi A, Moore F. Statistical analysis of accumulation and sources of heavy metals occurrence in agricultural soils of Khoshk River Banks, Shiraz, Iran. American-Eurasian Journal of Agricultural &Environmental Sciences, 2007, 2(5): 565- 573.

    [11] 吳紹華, 周生路, 楊得志, 廖富強(qiáng), 張紅富, 任奎. 宜興市近郊土壤重金屬來(lái)源與空間分布研究. 科學(xué)通報(bào), 2008, 53(增刊): 162- 170.

    [12] 陳同斌, 鄭袁明, 陳煌, 吳泓濤, 周建利, 羅金發(fā), 鄭國(guó)砥. 北京市不同土地利用類型的土壤砷含量特征. 地理研究, 2005, 24(2): 229- 235.

    [13] 郭朝暉, 肖細(xì)元, 陳同斌, 廖曉勇, 宋杰, 武斌. 湘江中下游農(nóng)田土壤和蔬菜的重金屬污染. 地理學(xué)報(bào), 2008, 63(1): 3- 11.

    [14] Zhao Y F, Shi X Z, Huang B, Yu D S, Wang H J, Sun W X, Oboern I, Blomback K. Spatial distribution of heavy metals in agricultural soils of an industry-based peri-urban area in Wuxi, China. Pedosphere, 2007, 17(1): 44- 51.

    [15] Wei C Y, Wang C, Yang L S. Characterizing spatial distribution and sources of heavy metals in the soils from mining-smelting activities in Shuikoushan, Hunan Province, China. Journal of Environmental Sciences, 2009, 21(9): 1230- 1236.

    [16] Notten M J M, Oosthoek A J P, Rozema J, Aerts R. Heavy metal concentrations in a soil-plant-snail food chain along a terrestrial soil pollution gradient. Environmental Pollution, 2005, 138(1): 178- 190.

    [17] Chen T B, Zheng Y M, Lei M, Huang Z C, Wu H T, Chen H, Fan K K, Yu K, Wu X, Tian Q Z. Assessment of heavy metal pollution in surface soils of urban parks in Beijing, China. Chemosphere, 2005, 60(4): 542- 551.

    [18] 陳同斌, 宋波, 鄭袁明, 黃澤春, 鄭國(guó)砥, 李艷霞, 雷梅, 廖曉勇. 北京市蔬菜和菜地土壤砷含量及其健康風(fēng)險(xiǎn)分析. 地理學(xué)報(bào), 2006, 61(3): 297- 310.

    [19] Li H B, Yu S, Li G L, Liu Y, Yu G B, Deng H, Wu S C, Wong M H. Urbanization increased metal levels in lake surface sediment and catchment topsoil of waterscape parks. Science of the Total Environment, 2012, 432: 202- 209.

    [20] Hakanson L. An ecological risk index for aquatic pollution control-a sedimentological approach. Water Research, 1980, 14(8): 975- 1001.

    [21] 于瑞蓮, 胡恭任, 戚紅璐, 陳福料. 泉州市不同功能區(qū)大氣將塵重金屬污染及生態(tài)風(fēng)險(xiǎn)評(píng)價(jià). 環(huán)境化學(xué), 2010, 29(6): 1086- 1090.

    [22] 陳振金, 陳春秀, 劉用清, 吳瑜端, 楊孫楷, 盧昌義. 福建省土壤環(huán)境背景值研究. 環(huán)境科學(xué), 1992, 13(4): 70- 75.

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

    [24] 胡二邦. 環(huán)境風(fēng)險(xiǎn)評(píng)價(jià)實(shí)用技術(shù)和方法[M]. 北京: 中國(guó)環(huán)境科學(xué)出版社, 2000.

    [25] VanLandeghem M M, Meyer M D, Cox S B, Sharma B, Patino R. Spatial and temporal patterns of surface water quality and ichthyotoxicity in urban and rural river basins in Texas. Water Research, 2012, 46(20): 6638- 6651.

    [26] Das B, Nordin R, Mazumder A. Watershed land use as a determinant of metal concentrations in freshwater systems. Environmental Geochemistry and Health, 2009, 31(6): 595- 607.

    [27] 孟昭虹, 周嘉, 鄭元福. 哈爾濱市城市土壤重金屬生態(tài)風(fēng)險(xiǎn)評(píng)價(jià). 水土保持研究, 2009, 16(2): 152- 159.

    [28] Yang Y G, Campbel C D, Clark L, Cameron C M, Paterson E. Microbial indicators of heavy metal contamination in urban and rural soils. Chemosphere, 2006, 63(11): 1942- 1952.

    [29] 武永鋒, 劉叢強(qiáng), 涂成龍. 貴陽(yáng)市土壤重金屬污染及其潛在生態(tài)風(fēng)險(xiǎn)評(píng)價(jià). 礦物巖石地球化學(xué)通報(bào), 2007, 26(3): 254- 257.

    [30] 姚榮江, 楊勁松, 孟慶峰, 張春銀, 吉榮龍. 蘇北沿海某灘涂區(qū)土壤重金屬含量及其污染評(píng)價(jià). 環(huán)境科學(xué)研究, 2012, 25(5): 512- 518.

    [31] Fairbrother A, Wenstel R, Sappington K, Wood W. Framework for metals risk assessment. Ecotoxicology and Environmental Safety, 2007, 68(2): 145- 227.

    [32] 孟慶峰, 楊勁松, 姚榮江, 李冬順, 張春銀, 吉榮龍. 灘涂圍墾區(qū)土壤重金屬調(diào)查及生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)——以鹽城市弶港鎮(zhèn)為例. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào), 2011, 30(11): 2249- 2257.

    [33] 肖鵬飛, 李法云, 付寶榮, 王效舉. 土壤重金屬污染及其植物修復(fù)研究. 遼寧大學(xué)學(xué)報(bào):自然科學(xué)版, 2004, 31(3): 279- 283.

    [34] 李波, 林玉鎖, 張孝飛, 徐亦鋼, 俞飛. 滬寧高速公路兩側(cè)土壤和小麥重金屬污染狀況. 農(nóng)村生態(tài)環(huán)境, 2005, 21(3): 50- 53, 70- 70.

    [35] 王軍, 陳振樓, 王初, 葉明武, 沈靜, 聶智凌. 上海崇明島蔬菜地土壤重金屬含量與生態(tài)風(fēng)險(xiǎn)預(yù)警評(píng)估. 環(huán)境科學(xué), 2007, 28(3): 647- 653.

    Distribution of heavy metals in topsoils affected by land use patterns at a small watershed scale: a case study in the Bantou Reservoir watershed in Xiamen,China

    LI Qingliang1,2, WU Qian2,3, GAO Jinbo2,3, MA Jun2, XU Qiufang1,*, YU Shen2,*

    1ZhejiangA&FUniversity,Hangzhou311300,China2KeyLaboratoryofUrbanEnvironmentandHealth,InstituteofUrbanEnvironment,ChineseAcademyofSciences,Xiamen361021,China3UniversityofChineseAcademyofSciences,Beijing100049,China

    Anthropogenic activities increase levels of heavy metals in catchment soils while changing land use in a watershed and further impairing the receiving water. Using a watershed as a study unit might help to determine the contributions of various anthropogenic activities, or their driven land use changes to spatial loadings and environmental risks of heavy metals in catchment soils because of its close geographic boundary. Few studies have been performed at a watershed scale to understand spatial relationships between anthropogenic land use changes and heavy metal enrichment in catchment soils, and have combined this with an overall potential ecological risk of heavy metal enrichment. This study integrated satellite image interpretation, field surveys of soil heavy metal level, and geo-statistics to explore relationships of soil heavy metal loads and their potential ecological risk in an urbanizing watershed in Xiamen, China. The Bantou Reservoir watershed covers 205 km2and consists of forested land, agricultural land, urban land, green land, and water surface. Interpreted by the digital elevation map (DEM), the studied watershed is comprised of five subwatersheds. In total, 150 soil-sampling points were assigned in the watershed based on coverage of each land use type with a minimum of three samples in each subwatershed to satisfy statistical requirements. A composite sample of two top- 10 cm soil samples was collected in a 10 × 10 m square at each sampling point. Six heavy metals, Cu, Zn, Pb, Cr, Ni, and Cd were analyzed in these topsoil samples. The enrichment factor and potential ecological risk factor of each heavy metal for each land use type, and total risk index of six heavy metals for each land use type were calculated. Results indicated that urban and agricultural land use significantly increased heavy metal loads in topsoils in comparison with forest and green land use in the watershed. The enrichment factor of heavy metals calculated on the basis of the regional background values showed that urban land use has higher values than agricultural land use, followed by green and forest land use. Accordingly, the potential ecological risk factor of each heavy metal followed the same pattern. Urban land use increased the enrichment of Cu, Zn, and Cd in topsoil of the watershed in comparison with other land uses. Lead enrichment factor was not related to land use type, suggesting that Pb was derived from atmospheric deposition. Regardless of land use types in the watershed, total Cr and Ni contents in topsoil were below the regional background values, although significant increases were observed within urban land in comparison with forest land. The land use based-inverse distance weighting (LU-IDW) spatial interpolation revealed that heavy metal loads and potential ecological risks are clearly related to anthropogenic activities, especially to urbanization. This study demonstrated that the methodology we applied to understand relationships between land use change with urbanization and heavy metal enrichment and even their ecological risks at a watershed scale is feasible. The integration of traditional field survey and satellite image interpretation techniques can serve as a powerful methodology for watershed-scale studies on environmental changes.

    heavy metal; topsoil; land use pattern; urbanization; potential ecological risk assessment

    國(guó)家自然科學(xué)基金面上項(xiàng)目(41371474);中國(guó)科學(xué)院知識(shí)創(chuàng)新工程重要方向項(xiàng)目(KZCX-YW-JC402);創(chuàng)新團(tuán)隊(duì)國(guó)際合作伙伴計(jì)劃(KZCX-YWT08)

    2014- 02- 25;

    2014- 06- 13

    10.5846/stxb201402250327

    *通訊作者Corresponding author.E-mail: xuqiufang@zafu.edu.cn;syu@iue.ac.cn

    李清良,吳倩,高進(jìn)波,馬軍,徐秋芳,俞慎.基于小流域尺度的土壤重金屬分布與土地利用相關(guān)性研究——以廈門市坂頭水庫(kù)流域?yàn)槔?生態(tài)學(xué)報(bào),2015,35(16):5486- 5494.

    Li Q L, Wu Q, Gao J B, Ma J, Xu Q F, Yu S.Distribution of heavy metals in topsoils affected by land use patterns at a small watershed scale: a case study in the Bantou Reservoir watershed in Xiamen, China.Acta Ecologica Sinica,2015,35(16):5486- 5494.

    猜你喜歡
    生態(tài)研究
    FMS與YBT相關(guān)性的實(shí)證研究
    “生態(tài)養(yǎng)生”娛晚年
    2020年國(guó)內(nèi)翻譯研究述評(píng)
    遼代千人邑研究述論
    住進(jìn)呆萌生態(tài)房
    生態(tài)之旅
    生態(tài)之旅
    生態(tài)之旅
    視錯(cuò)覺在平面設(shè)計(jì)中的應(yīng)用與研究
    科技傳播(2019年22期)2020-01-14 03:06:54
    EMA伺服控制系統(tǒng)研究

    99视频精品全部免费 在线| 五月玫瑰六月丁香| 国产伦精品一区二区三区视频9| 一级,二级,三级黄色视频| 亚洲国产精品一区二区三区在线| 在线观看免费视频网站a站| 综合色丁香网| 91在线精品国自产拍蜜月| 亚洲精品国产色婷婷电影| 91久久精品国产一区二区三区| 免费人妻精品一区二区三区视频| 内射极品少妇av片p| 中文字幕亚洲精品专区| 成人二区视频| 建设人人有责人人尽责人人享有的| 亚洲av成人精品一区久久| 不卡视频在线观看欧美| 99视频精品全部免费 在线| 在线观看免费日韩欧美大片 | 蜜桃久久精品国产亚洲av| 乱人伦中国视频| 成人综合一区亚洲| 亚洲av二区三区四区| 国产成人精品久久久久久| 成人毛片60女人毛片免费| 看十八女毛片水多多多| 久久女婷五月综合色啪小说| 2022亚洲国产成人精品| 九九久久精品国产亚洲av麻豆| 久久久国产一区二区| xxx大片免费视频| 搡女人真爽免费视频火全软件| 国产在线一区二区三区精| 99久久精品国产国产毛片| 精品一区二区免费观看| 如日韩欧美国产精品一区二区三区 | 人人妻人人澡人人看| 亚洲av成人精品一二三区| 亚洲精品色激情综合| 国国产精品蜜臀av免费| 国产亚洲午夜精品一区二区久久| 91久久精品国产一区二区三区| 男男h啪啪无遮挡| 少妇 在线观看| 免费黄频网站在线观看国产| a级片在线免费高清观看视频| 久久久久久久久久久免费av| .国产精品久久| 国产亚洲5aaaaa淫片| 亚洲av成人精品一二三区| 日韩制服骚丝袜av| 日韩伦理黄色片| 九九久久精品国产亚洲av麻豆| 欧美日韩精品成人综合77777| 看非洲黑人一级黄片| 看非洲黑人一级黄片| 亚洲欧美成人综合另类久久久| 国产一区亚洲一区在线观看| xxx大片免费视频| 午夜免费男女啪啪视频观看| 国产亚洲精品久久久com| 国产黄色视频一区二区在线观看| 久久人人爽人人片av| 人妻系列 视频| 国产欧美日韩精品一区二区| 亚洲成人av在线免费| 日本av免费视频播放| a级毛片免费高清观看在线播放| 99久久综合免费| 水蜜桃什么品种好| 亚洲图色成人| 成人午夜精彩视频在线观看| 九九在线视频观看精品| 国内少妇人妻偷人精品xxx网站| 如日韩欧美国产精品一区二区三区 | 狂野欧美白嫩少妇大欣赏| 中文在线观看免费www的网站| 精品亚洲成a人片在线观看| 国产亚洲av片在线观看秒播厂| 极品少妇高潮喷水抽搐| 97在线人人人人妻| 日本与韩国留学比较| 菩萨蛮人人尽说江南好唐韦庄| 国产精品三级大全| 精品久久久精品久久久| 人妻系列 视频| 国产深夜福利视频在线观看| 久久久国产欧美日韩av| 大又大粗又爽又黄少妇毛片口| 精品国产一区二区三区久久久樱花| 观看av在线不卡| 国产高清国产精品国产三级| 国产精品99久久99久久久不卡 | 精品少妇久久久久久888优播| 精品人妻熟女av久视频| 亚洲国产最新在线播放| 香蕉精品网在线| 欧美+日韩+精品| av在线观看视频网站免费| a级毛片免费高清观看在线播放| 国内精品宾馆在线| 91精品国产九色| 午夜影院在线不卡| 色5月婷婷丁香| 女人久久www免费人成看片| 男女免费视频国产| 熟女电影av网| 韩国高清视频一区二区三区| 亚洲va在线va天堂va国产| 色视频在线一区二区三区| 一边亲一边摸免费视频| 精品少妇内射三级| 女的被弄到高潮叫床怎么办| 国产欧美日韩精品一区二区| 大陆偷拍与自拍| 成年美女黄网站色视频大全免费 | 久久久精品免费免费高清| 色94色欧美一区二区| 交换朋友夫妻互换小说| 男人舔奶头视频| 国产亚洲5aaaaa淫片| 日韩精品有码人妻一区| 日本vs欧美在线观看视频 | 日韩欧美精品免费久久| 乱系列少妇在线播放| 丝袜在线中文字幕| 午夜免费观看性视频| 欧美日本中文国产一区发布| 久久久午夜欧美精品| 亚洲av.av天堂| 国产黄片视频在线免费观看| 曰老女人黄片| 三上悠亚av全集在线观看 | 内射极品少妇av片p| 亚洲av二区三区四区| 我的老师免费观看完整版| 国产一区二区在线观看av| 国产视频首页在线观看| 秋霞在线观看毛片| 国产成人午夜福利电影在线观看| 亚洲精品日韩在线中文字幕| 黑人猛操日本美女一级片| 久久久久网色| 日韩精品免费视频一区二区三区 | 日韩欧美精品免费久久| 中文在线观看免费www的网站| 插逼视频在线观看| 一区二区av电影网| 午夜福利视频精品| 国产精品国产三级国产专区5o| 香蕉精品网在线| 国产深夜福利视频在线观看| 十八禁网站网址无遮挡 | 王馨瑶露胸无遮挡在线观看| 嘟嘟电影网在线观看| 91精品一卡2卡3卡4卡| 久久久国产精品麻豆| 日韩中字成人| 国产精品一区二区在线不卡| 精品人妻偷拍中文字幕| 国产精品一区二区三区四区免费观看| 我的老师免费观看完整版| 一级毛片 在线播放| 久久人人爽av亚洲精品天堂| 成年人午夜在线观看视频| 午夜影院在线不卡| 三级国产精品欧美在线观看| 日韩三级伦理在线观看| 性色avwww在线观看| 青春草国产在线视频| 久久久久久久久久人人人人人人| 中文资源天堂在线| 久久久久久久久久久久大奶| 一本色道久久久久久精品综合| 人人澡人人妻人| 国产成人aa在线观看| 免费少妇av软件| 少妇高潮的动态图| 日产精品乱码卡一卡2卡三| 亚洲美女黄色视频免费看| 中文字幕制服av| 国产伦理片在线播放av一区| 熟女av电影| 十八禁网站网址无遮挡 | 久久久久久久精品精品| 女性被躁到高潮视频| 高清在线视频一区二区三区| 在线播放无遮挡| 欧美一级a爱片免费观看看| 国产精品嫩草影院av在线观看| 日韩一区二区三区影片| 男女无遮挡免费网站观看| 99视频精品全部免费 在线| 免费播放大片免费观看视频在线观看| 青青草视频在线视频观看| 国产黄色免费在线视频| 成年av动漫网址| 伊人久久国产一区二区| 日韩,欧美,国产一区二区三区| av网站免费在线观看视频| 国产精品免费大片| 亚洲精品国产色婷婷电影| 午夜日本视频在线| 亚洲怡红院男人天堂| 免费观看无遮挡的男女| 妹子高潮喷水视频| 91精品一卡2卡3卡4卡| 少妇人妻一区二区三区视频| 高清不卡的av网站| 香蕉精品网在线| 久久毛片免费看一区二区三区| 久久99一区二区三区| 少妇人妻一区二区三区视频| 国产日韩欧美在线精品| 在线精品无人区一区二区三| 午夜av观看不卡| 最近中文字幕高清免费大全6| 99热6这里只有精品| 久久久精品免费免费高清| 日产精品乱码卡一卡2卡三| av国产精品久久久久影院| 国产爽快片一区二区三区| 精品酒店卫生间| 欧美bdsm另类| 中国国产av一级| 色视频在线一区二区三区| 日韩不卡一区二区三区视频在线| 亚洲精品国产成人久久av| 青青草视频在线视频观看| 一本久久精品| 中文字幕亚洲精品专区| 在线精品无人区一区二区三| 亚洲精品久久久久久婷婷小说| 亚洲精品久久午夜乱码| 亚洲不卡免费看| 亚洲成色77777| av免费在线看不卡| 国产中年淑女户外野战色| 亚洲国产日韩一区二区| 免费不卡的大黄色大毛片视频在线观看| 成年人午夜在线观看视频| 中文资源天堂在线| 日本wwww免费看| 午夜福利影视在线免费观看| 国产又色又爽无遮挡免| 国产一区二区在线观看日韩| 一级毛片我不卡| 免费人成在线观看视频色| 欧美日韩在线观看h| 99九九在线精品视频 | 亚洲精品久久久久久婷婷小说| 久久韩国三级中文字幕| 精品亚洲成国产av| 在线观看av片永久免费下载| 黑人猛操日本美女一级片| 99热这里只有精品一区| 日韩成人伦理影院| 2022亚洲国产成人精品| 18禁在线播放成人免费| 麻豆成人午夜福利视频| 欧美激情国产日韩精品一区| 99久久精品热视频| 欧美精品亚洲一区二区| 少妇高潮的动态图| 大码成人一级视频| 下体分泌物呈黄色| 亚洲国产日韩一区二区| 国产成人精品婷婷| 国产精品国产三级国产专区5o| 国产探花极品一区二区| 草草在线视频免费看| 午夜老司机福利剧场| videos熟女内射| 日本-黄色视频高清免费观看| 啦啦啦在线观看免费高清www| 日韩欧美精品免费久久| 国产亚洲av片在线观看秒播厂| 国产精品99久久久久久久久| 亚洲精品自拍成人| 99久久中文字幕三级久久日本| 亚洲成人手机| 纯流量卡能插随身wifi吗| 亚洲性久久影院| 久久99蜜桃精品久久| 天美传媒精品一区二区| 久久精品国产亚洲网站| 国产探花极品一区二区| 少妇人妻 视频| 校园人妻丝袜中文字幕| 99热这里只有是精品在线观看| 伦精品一区二区三区| 少妇人妻一区二区三区视频| av又黄又爽大尺度在线免费看| 中文欧美无线码| 国产亚洲欧美精品永久| 最近手机中文字幕大全| 久久久国产欧美日韩av| 在线看a的网站| 成人漫画全彩无遮挡| 免费黄网站久久成人精品| 欧美最新免费一区二区三区| 久久久久视频综合| 51国产日韩欧美| 大又大粗又爽又黄少妇毛片口| 婷婷色综合www| 丁香六月天网| 男的添女的下面高潮视频| 观看美女的网站| 午夜精品国产一区二区电影| 九九在线视频观看精品| 18禁裸乳无遮挡动漫免费视频| 99久久精品国产国产毛片| 中文字幕人妻熟人妻熟丝袜美| 午夜福利在线观看免费完整高清在| 国产视频首页在线观看| 日日啪夜夜爽| 亚洲精品aⅴ在线观看| 久久综合国产亚洲精品| 久久精品国产亚洲av天美| 日本免费在线观看一区| 黑人猛操日本美女一级片| 国产精品蜜桃在线观看| 精品国产露脸久久av麻豆| 一级毛片 在线播放| 黄色视频在线播放观看不卡| 色5月婷婷丁香| 99久久精品热视频| 狂野欧美激情性bbbbbb| 少妇高潮的动态图| 亚洲内射少妇av| 香蕉精品网在线| 欧美最新免费一区二区三区| 26uuu在线亚洲综合色| 看十八女毛片水多多多| 国产极品粉嫩免费观看在线 | 天堂中文最新版在线下载| 中文精品一卡2卡3卡4更新| 麻豆乱淫一区二区| 伦理电影大哥的女人| 亚洲精品乱久久久久久| 亚洲欧美一区二区三区国产| 十八禁网站网址无遮挡 | 久久精品国产亚洲网站| 中国美白少妇内射xxxbb| 搡女人真爽免费视频火全软件| 国产精品久久久久成人av| 国产成人午夜福利电影在线观看| 99热国产这里只有精品6| 精品久久久精品久久久| 国产午夜精品一二区理论片| 国产精品福利在线免费观看| 国语对白做爰xxxⅹ性视频网站| 国产又色又爽无遮挡免| 久久久久久久久大av| 麻豆乱淫一区二区| 黄片无遮挡物在线观看| 国产69精品久久久久777片| 日韩人妻高清精品专区| 少妇熟女欧美另类| 国产精品无大码| 99热这里只有是精品在线观看| 汤姆久久久久久久影院中文字幕| 热re99久久精品国产66热6| av播播在线观看一区| 国产精品一二三区在线看| 免费人妻精品一区二区三区视频| 亚洲精品456在线播放app| 如日韩欧美国产精品一区二区三区 | 老司机影院毛片| 爱豆传媒免费全集在线观看| 欧美日韩视频高清一区二区三区二| 九色成人免费人妻av| 韩国高清视频一区二区三区| 精品国产一区二区三区久久久樱花| 搡老乐熟女国产| 熟妇人妻不卡中文字幕| 一级,二级,三级黄色视频| 高清欧美精品videossex| 国产成人freesex在线| 丰满人妻一区二区三区视频av| 青春草国产在线视频| 亚洲成色77777| 亚洲欧洲日产国产| 午夜精品国产一区二区电影| 全区人妻精品视频| 2018国产大陆天天弄谢| 麻豆乱淫一区二区| 在线天堂最新版资源| 日韩强制内射视频| 18禁动态无遮挡网站| 亚洲第一区二区三区不卡| 青春草亚洲视频在线观看| av线在线观看网站| 自拍偷自拍亚洲精品老妇| 免费观看的影片在线观看| 人人妻人人爽人人添夜夜欢视频 | 国产高清国产精品国产三级| 乱人伦中国视频| 国产极品天堂在线| 免费观看性生交大片5| av在线老鸭窝| 国产熟女午夜一区二区三区 | 免费黄色在线免费观看| 国产熟女欧美一区二区| 国产综合精华液| 国产淫片久久久久久久久| 十八禁高潮呻吟视频 | 亚洲国产精品成人久久小说| 五月开心婷婷网| 日本av手机在线免费观看| 欧美97在线视频| 黄片无遮挡物在线观看| 日韩欧美精品免费久久| 久久精品熟女亚洲av麻豆精品| 自拍偷自拍亚洲精品老妇| 美女福利国产在线| 国产在线男女| 最近最新中文字幕免费大全7| 国产无遮挡羞羞视频在线观看| 免费观看在线日韩| 青青草视频在线视频观看| 午夜免费观看性视频| 色吧在线观看| 国产精品久久久久久久久免| 我的老师免费观看完整版| av国产精品久久久久影院| a级一级毛片免费在线观看| 中国三级夫妇交换| 成年av动漫网址| 一区在线观看完整版| 最近2019中文字幕mv第一页| 国产精品.久久久| 亚洲欧洲精品一区二区精品久久久 | 成年av动漫网址| 99re6热这里在线精品视频| 欧美精品国产亚洲| 亚洲国产成人一精品久久久| 久久久久国产网址| 我的女老师完整版在线观看| 99热6这里只有精品| 人妻制服诱惑在线中文字幕| 久久久久视频综合| 男人和女人高潮做爰伦理| 中文字幕久久专区| 国产av码专区亚洲av| 日韩强制内射视频| 免费少妇av软件| 国产一区二区三区av在线| 国产老妇伦熟女老妇高清| 人妻一区二区av| a级毛片免费高清观看在线播放| 乱人伦中国视频| 最新的欧美精品一区二区| 国产片特级美女逼逼视频| 少妇高潮的动态图| 亚洲欧洲日产国产| 乱码一卡2卡4卡精品| 日韩精品有码人妻一区| 在线免费观看不下载黄p国产| av有码第一页| xxx大片免费视频| 国产精品伦人一区二区| 婷婷色麻豆天堂久久| 国产黄片视频在线免费观看| 午夜福利,免费看| 2018国产大陆天天弄谢| 国产精品麻豆人妻色哟哟久久| 蜜臀久久99精品久久宅男| 人妻少妇偷人精品九色| 日本免费在线观看一区| 亚洲婷婷狠狠爱综合网| 国产在线免费精品| 亚洲欧洲精品一区二区精品久久久 | 国产黄色免费在线视频| 久久久久国产精品人妻一区二区| 亚洲精品亚洲一区二区| 一级黄片播放器| 免费高清在线观看视频在线观看| 天堂俺去俺来也www色官网| 久久99精品国语久久久| 欧美3d第一页| 久久精品国产自在天天线| 久久精品久久久久久久性| 久久99蜜桃精品久久| 国产一级毛片在线| 欧美精品一区二区免费开放| 日韩 亚洲 欧美在线| 国产精品一区二区性色av| 欧美 亚洲 国产 日韩一| 成人亚洲精品一区在线观看| 人人澡人人妻人| a级一级毛片免费在线观看| 久久久久视频综合| 久久精品国产亚洲av天美| 少妇丰满av| 丝袜在线中文字幕| 国产永久视频网站| 夜夜骑夜夜射夜夜干| 久久久久久久久久人人人人人人| 在线观看免费高清a一片| 黄色日韩在线| 欧美日韩av久久| 五月伊人婷婷丁香| 中文字幕久久专区| 最近中文字幕高清免费大全6| 一级毛片黄色毛片免费观看视频| 国产乱来视频区| 日韩制服骚丝袜av| 免费看不卡的av| 亚洲婷婷狠狠爱综合网| 亚洲精品亚洲一区二区| 乱人伦中国视频| 欧美少妇被猛烈插入视频| av有码第一页| 国产美女午夜福利| 精品国产一区二区久久| 欧美变态另类bdsm刘玥| 91精品国产国语对白视频| 日韩欧美 国产精品| 欧美日韩亚洲高清精品| 日韩av免费高清视频| 人妻制服诱惑在线中文字幕| 麻豆成人午夜福利视频| av福利片在线| 久久这里有精品视频免费| 黄色毛片三级朝国网站 | 国产日韩欧美视频二区| 久久久久久久久大av| 亚洲欧美成人精品一区二区| 日韩欧美精品免费久久| 亚洲精品,欧美精品| 中文在线观看免费www的网站| 我要看日韩黄色一级片| 在线天堂最新版资源| 日本黄色片子视频| 亚洲自偷自拍三级| 成年人免费黄色播放视频 | 久久久久久久久久人人人人人人| 久久人人爽人人片av| 久久久久久伊人网av| 久久亚洲国产成人精品v| 国语对白做爰xxxⅹ性视频网站| av在线观看视频网站免费| 在线观看人妻少妇| 欧美日韩视频高清一区二区三区二| 免费观看在线日韩| 久热久热在线精品观看| 大片电影免费在线观看免费| 伦理电影大哥的女人| 亚洲欧美一区二区三区国产| 极品人妻少妇av视频| 亚洲精品一二三| 婷婷色综合www| 国产 一区精品| 国产一区有黄有色的免费视频| 校园人妻丝袜中文字幕| 欧美精品一区二区大全| 黑人巨大精品欧美一区二区蜜桃 | 桃花免费在线播放| 欧美激情国产日韩精品一区| 成年人午夜在线观看视频| 最新中文字幕久久久久| 久久久欧美国产精品| 人妻人人澡人人爽人人| 少妇人妻一区二区三区视频| 久久精品熟女亚洲av麻豆精品| 久久韩国三级中文字幕| 日日啪夜夜撸| 欧美精品高潮呻吟av久久| 黄色怎么调成土黄色| 午夜福利网站1000一区二区三区| 内地一区二区视频在线| 狂野欧美激情性xxxx在线观看| 欧美区成人在线视频| 国产无遮挡羞羞视频在线观看| 国产成人精品婷婷| 色94色欧美一区二区| 人妻系列 视频| 日韩中文字幕视频在线看片| 人妻制服诱惑在线中文字幕| 最近2019中文字幕mv第一页| 亚洲av男天堂| 久久精品国产亚洲av涩爱| 国产精品一区www在线观看| 精品一区二区免费观看| 国产成人一区二区在线| 男人爽女人下面视频在线观看| 成人国产麻豆网| 免费看不卡的av| 亚洲四区av| 日韩精品有码人妻一区| 久久亚洲国产成人精品v| 国产精品嫩草影院av在线观看| 又大又黄又爽视频免费| 97超碰精品成人国产| 99久国产av精品国产电影| 3wmmmm亚洲av在线观看| 国产成人午夜福利电影在线观看| 欧美 亚洲 国产 日韩一| 少妇丰满av| a级片在线免费高清观看视频| 日本欧美国产在线视频| 日韩不卡一区二区三区视频在线| 亚洲成人av在线免费| 国产视频内射| 少妇的逼好多水| 国产精品一区二区在线观看99| 另类精品久久| 国国产精品蜜臀av免费| 免费观看av网站的网址| 自线自在国产av| 免费大片18禁| 美女大奶头黄色视频| 免费黄频网站在线观看国产| 91久久精品国产一区二区成人|