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

    河流岸帶濕地沉積物重金屬分布對(duì)植被物種多樣性和底棲動(dòng)物群落特征的影響

    2021-04-07 03:19:36尤曉光劉靜玲
    生態(tài)學(xué)報(bào) 2021年5期
    關(guān)鍵詞:灤河沉積物河流

    馬 康,史 璇,尤曉光,3,劉靜玲,*

    1 北京師范大學(xué)環(huán)境學(xué)院水環(huán)境模擬國(guó)家重點(diǎn)實(shí)驗(yàn)室, 北京 100875 2 北京師范大學(xué)流域環(huán)境生態(tài)工程研發(fā)中心, 珠海 519087 3 交通運(yùn)輸部水運(yùn)科學(xué)研究院, 北京 100088

    河流生態(tài)系統(tǒng)中自然存在重金屬微量元素,主要來(lái)源于巖石侵蝕、水動(dòng)力遷移等地球化學(xué)作用。隨著社會(huì)經(jīng)濟(jì)發(fā)展,人類生產(chǎn)生活直接或間接向自然水體中排放大量重金屬,成為水體中重金屬的主要來(lái)源,對(duì)水生態(tài)系統(tǒng)造成不同程度生態(tài)風(fēng)險(xiǎn)。重金屬污染物具有高毒性、不可降解性和生物蓄積性等特點(diǎn),進(jìn)入河流后大量滯留在水體中,被水體中懸浮物吸附,沉積在河流沉積物中[1]。長(zhǎng)期沉積作用使河流沉積物成為河流重金屬主要受體,當(dāng)水環(huán)境條件改變,重金屬重新進(jìn)入水體,形成“二次污染”。因而,河流沉積物不僅是重金屬主要儲(chǔ)存庫(kù),也是重金屬潛在釋放源,即使污染源消除,受重金屬污染的沉積物也可以長(zhǎng)時(shí)間緩慢釋放重金屬,從而延長(zhǎng)重金屬污染時(shí)間,對(duì)河流生態(tài)系統(tǒng)健康具有長(zhǎng)期影響。

    河流沉積物不僅是河流重金屬污染物的重要貯存庫(kù)和釋放源,還是河流濕地生物如植物、大型無(wú)脊椎動(dòng)物的重要棲息地[2- 3]。沉積物釋放的重金屬通過(guò)水流遷移或富集作用進(jìn)入植物體內(nèi),對(duì)生物體產(chǎn)生毒性效應(yīng),使得抗逆性低的植被死亡,降低植物群落多樣性[4]。沉積物-水界面是底棲動(dòng)物主要生活場(chǎng)所,重金屬污染的細(xì)沉積物可被底棲動(dòng)物吞食[5],使無(wú)脊椎動(dòng)物生長(zhǎng)受到抑制,底棲動(dòng)物群落受到破壞[6- 7]。例如Costas等[8]研究發(fā)現(xiàn)沉積物重金屬含量與底棲生物群落結(jié)構(gòu)多樣性指標(biāo)顯著相關(guān),重金屬污染點(diǎn)位EPT昆蟲(chóng)(蜉蝣目Ephemeroptera、襀翅目Plecoptera和毛翅目Trichoptera)數(shù)目顯著高于參照點(diǎn)。類似的, Bere等[9]研究表明沉積物金屬改變了大型底棲無(wú)脊椎動(dòng)物群落物種組成,且重金屬污染程度較低時(shí)對(duì)底棲動(dòng)物群落組成的影響依舊顯著。

    灤河隸屬于海河流域,是京津冀城市群的重要水源地,也是北方重要的生物棲息地,其生態(tài)健康對(duì)于區(qū)域經(jīng)濟(jì)、社會(huì)以及生態(tài)等方面可持續(xù)發(fā)展具有重要意義。進(jìn)入21世紀(jì),隨著城市化進(jìn)程和工農(nóng)業(yè)發(fā)展,灤河水環(huán)境受到人類活動(dòng)嚴(yán)重影響,例如,流域徑流減少[10],部分河段富營(yíng)養(yǎng)化[11],水體重金屬含量超標(biāo)[12- 13],水和沉積物中存在化學(xué)品多環(huán)芳烴(PAHs)污染的潛在生態(tài)風(fēng)險(xiǎn)[14]。河流污染物對(duì)生物群落具有生物毒性和潛在生態(tài)風(fēng)險(xiǎn),因而,明確沉積物重金屬分布與濕地生物群落特征的關(guān)系是當(dāng)前河流生態(tài)管理和保護(hù)亟待解決的問(wèn)題。

    本研究以灤河岸帶濕地為研究對(duì)象,采用沉積物生物毒性效應(yīng)系數(shù)法和內(nèi)梅羅指數(shù)法評(píng)價(jià)沉積物重金屬污染生物毒性和生態(tài)風(fēng)險(xiǎn),通過(guò)計(jì)算植被物種多樣性指數(shù)和底棲動(dòng)物完整性指數(shù)探究灤河植被和底棲動(dòng)物群落特征,分析沉積物重金屬分布、植被和底棲群落特征的關(guān)系,旨在為灤河流域河流生態(tài)風(fēng)險(xiǎn)管理提供科學(xué)依據(jù)。

    1 材料與方法

    1.1 研究區(qū)概況

    灤河流域(115°34′—119°50′E,39°02′—42°43′N)地跨蒙古高原和華北平原,地勢(shì)由西北向東南傾斜,面積約4.47萬(wàn)km2。灤河發(fā)源于河北省豐寧滿族自治縣西北巴彥古圖爾山麓,上游屬內(nèi)蒙古高原區(qū),中游為華北山地丘陵區(qū),下游為冀北平原區(qū),由河北樂(lè)亭縣入渤海,水系呈羽狀分布,全長(zhǎng)約888 km。灤河流域年均降水量553.2 mm,夏季多冬季少,汛期為7—8月,枯水期為1—2月,受降雨量影響,徑流量年內(nèi)變化較大,平均徑流量47.9億m3,氣候?yàn)榈湫蜏貛Т箨懶约撅L(fēng)氣候[15]。

    1.2 樣點(diǎn)布設(shè)和數(shù)據(jù)獲取

    根據(jù)已有研究對(duì)灤河流域生態(tài)分區(qū),結(jié)合兩岸地形、氣候、植被、土壤類型等特點(diǎn),在河流上游高原區(qū)(L1—L4)、中游山區(qū)(L5—L10)、下游平原區(qū)(L11—L14),依次布設(shè)L1-大灘鎮(zhèn)(41°34′42″N, 115°59′18″E)、L2-白城子(41°19′05″N, 116°28′31″E)、L3-石人溝(41°21′23″N, 116°20′35″E)、L4-紅旗營(yíng)房(42°04′16″N, 116°42′37″E)、L5-外溝門(41°52′148″N, 115°33′32″E)、L6-郭家屯(41°34′53″N, 117°05′53″E)、L7-太平莊(41°13′28″N, 117°22′50″E)、L8-西溝(41°08′55″N, 117°29′52″E)、L9-張百灣(40°59′42″N, 117°29′30″E)、L10三道河(40°58′04″N, 117°42′52″E)、L11-遷西大橋(40°09′37″N, 118°18′43″E)、L12-馬蘭莊(40°06′53″N, 118°37′01″E)、L13-王家樓(39°34′37″N, 118°52′49″E)、L14-姜各莊(39°27′39″N, 119°08′12″E)共計(jì)14個(gè)樣地,于2014年7月從灤河源頭到河口進(jìn)行野外采樣和實(shí)地調(diào)查。

    圖1 灤河采樣點(diǎn)分布圖Fig.1 Locations of sampling sites of Luanhe River

    采用樣方法在灤河岸帶植被類型樣地進(jìn)行取樣,植被調(diào)查采樣面積參照草本層設(shè)計(jì),每個(gè)樣點(diǎn)處設(shè)置10個(gè)2 m×2 m樣方,總計(jì)140個(gè)樣方,記錄樣方植被覆蓋度、物種豐度、物種名稱等植被信息。底棲動(dòng)物樣品利用直徑1/16 m2采泥器在河岸帶區(qū)域采集,每個(gè)點(diǎn)位采集多于3個(gè)子樣品混合,樣品經(jīng)過(guò)0.595 mm篩網(wǎng)過(guò)濾,保存于75%乙醇,利用解剖鏡或顯微鏡分類鑒定底棲動(dòng)物分類單元(種或?qū)?,稱量?jī)糁?。沉積物樣品利用彼得遜采泥器在14個(gè)樣地采集0—10 cm深度表層沉積物,每個(gè)樣地在河岸帶采集3個(gè)樣品,同一樣地沉積物樣品混合后放入聚乙烯塑料袋中密封保存,在實(shí)驗(yàn)室風(fēng)干,剔除動(dòng)植物殘?bào)w等雜質(zhì)待測(cè)。沉積物樣品過(guò)100目尼龍篩,采用HF-HNO3-HClO4三酸混合,經(jīng)石墨消解儀(SH220N)消解后,利用島津火焰原子吸收分光光度計(jì)(GFA- 6880)測(cè)定沉積物重金屬Ni、Cr、Zn、Cu、Pb和Cd的含量[16]。實(shí)驗(yàn)中所用試劑均為優(yōu)級(jí)純,用水為超純水,分析過(guò)程中采用沉積物標(biāo)準(zhǔn)參考物質(zhì)(GSD- 9和GSD- 11)進(jìn)行質(zhì)量控制,各元素回收率在94%—105%之間,測(cè)量誤差均小于10%。

    1.3 數(shù)據(jù)統(tǒng)計(jì)與分析

    應(yīng)用EXCEL 2016軟件對(duì)所有數(shù)據(jù)進(jìn)行處理,采用SPSS 23.0軟件進(jìn)行統(tǒng)計(jì)學(xué)分析,采用Origin 2017軟件繪圖,進(jìn)行采用Canoco 5.0對(duì)生物群落指標(biāo)與沉積物環(huán)境因子進(jìn)行主成分分析(Principal Component Analysis,PCA)。

    1.3.1重金屬污染評(píng)價(jià)

    根據(jù)沉積物質(zhì)量基準(zhǔn)值(sediment quality guideline values, SQGs)[17]的生物毒性效應(yīng)范圍值計(jì)算效應(yīng)系數(shù)(effects range-median quotient, ERMQ)評(píng)價(jià)不同采樣點(diǎn)沉積物生物毒性危害等級(jí)[18],采用綜合潛在生態(tài)風(fēng)險(xiǎn)指數(shù)(Integrated potential ecological risk index,RI)對(duì)沉積物重金屬綜合生態(tài)風(fēng)險(xiǎn)進(jìn)行評(píng)價(jià)[19],其中ERMQ和RI評(píng)價(jià)等級(jí)劃分標(biāo)準(zhǔn)見(jiàn)表1,計(jì)算公式如下:

    表1 重金屬污染評(píng)價(jià)等級(jí)劃分

    1.3.2植被物種多樣性評(píng)價(jià)

    植被物種多樣性評(píng)價(jià)采用以下3種指標(biāo): 物種豐富度指數(shù)(S)、Shannon-Wiener指數(shù)(H)、Pielou均勻度指數(shù)(J),計(jì)算公式如下[22]:

    J=H′/lnS

    H′=-∑pilnpi

    其中,pi表示種i出現(xiàn)的頻度,即重要值,S為物種數(shù)(物種豐富度指數(shù)),N為個(gè)體總數(shù)。

    1.3.3底棲動(dòng)物完整性評(píng)價(jià)

    研究團(tuán)隊(duì)選取表征底棲動(dòng)物多樣性和豐度(Diversity and Abundance, D/A)、耐污敏感性(Sensitivity and Tolerance, S/T)及功能攝食類群(Functional Feeding Groups, FFGs)的31個(gè)候選指標(biāo),通過(guò)分布范圍分析、敏感性分析和冗余分析篩選顯著區(qū)分參照點(diǎn)和受損點(diǎn)的核心指標(biāo),核心指標(biāo)由蜉蝣目、襀翅目和毛翅目個(gè)體數(shù)之和百分比(EPT%)、物種耐污性綜合指數(shù)BI (Biotic Index)和收集者攝食類群百分比(Collect-Gatherers%)三個(gè)指標(biāo)組成,通過(guò)驗(yàn)證和標(biāo)準(zhǔn)化構(gòu)建底棲動(dòng)物完整性指數(shù)(B-IBI, Benthic Macroinvertebrate-based Index of Biotic Integrity),詳細(xì)方法見(jiàn)團(tuán)隊(duì)已發(fā)表的相關(guān)文獻(xiàn)[23- 24]。

    2 研究結(jié)果

    2.1 灤河表層沉積物重金屬分布特征

    灤河沉積物重金屬元素含量見(jiàn)表2,灤河沉積物重金屬元素Ni、Cr、Zn、Cu、Pb和Cd的平均值分別為33.66 mg/kg、37.64 mg/kg、36.06 mg/kg、12.79 mg/kg、6.76 mg/kg和0.09 mg/kg,其中僅Ni元素超過(guò)了背景值,說(shuō)明灤河水系無(wú)明顯重金屬污染。重金屬Ni、Cr、Zn、Cu、Pb和Cd的變異系數(shù)分別為0.42、0.43、0.78、1.39、0.65和0.60,說(shuō)明灤河水系沉積物重金屬空間分布具有較大差異,特別是Zn和Cu離散程度較大,分布極不均勻。查閱文獻(xiàn),灤河沉積物重金屬歷史值見(jiàn)表2,2010年重金屬Cr、Cu、Cd含量分別是背景值1.98、5.21,2.55倍,重金屬污染十分嚴(yán)重,2013年灤河沉積物重金屬含量顯著降低,本研究中沉積物Ni、Cr、Zn、Cu、Pb和Cd分別為背景值的1.09、0.55、0.46、0.59、0.31、0.93、0.02、0.09倍,沉積物質(zhì)量明顯改善。

    灤河沉積物重金屬效應(yīng)系數(shù)ERMQ和潛在生態(tài)風(fēng)險(xiǎn)指數(shù)RI結(jié)果如圖2所示。研究結(jié)果表明,灤河水系上游、中游和下游沉積物重金屬ERMQ范圍為0.09—0.12、0.09—0.16和0.07—0.28,在下游L14-姜各莊點(diǎn)位最低,在下游L12-馬蘭莊鎮(zhèn)點(diǎn)位最高。上、中、下游點(diǎn)位均位于低級(jí)和中低級(jí),但下游點(diǎn)位總體高于上中游點(diǎn)位,且波動(dòng)較大。綜合潛在風(fēng)險(xiǎn)指數(shù)上游、中游和下游范圍分別為31.13—47.75、23.35—66.11、22.49—112.37,在下游L14-姜各莊點(diǎn)位最低,在下游L12-馬蘭莊鎮(zhèn)點(diǎn)位最高,與ERMQ評(píng)價(jià)結(jié)果一致,所有點(diǎn)位風(fēng)險(xiǎn)等級(jí)處于低風(fēng)險(xiǎn)水平,下游重金屬潛在風(fēng)險(xiǎn)總體高于上中游。對(duì)不同河段沉積物重金屬效應(yīng)系數(shù)ERMQ和潛在生態(tài)風(fēng)險(xiǎn)指數(shù)RI進(jìn)行非參數(shù)檢驗(yàn)(Kruskal-Wallis法),檢驗(yàn)結(jié)果表明,灤河干流上中下游沉積物沉積物重金屬效應(yīng)系數(shù)ERMQ具有顯著性差異(P<0.05)。

    表2 灤河沉積物樣品重金屬元素含量(mg/kg)統(tǒng)計(jì)

    圖2 灤河沉積物重金屬效應(yīng)系數(shù)ERMQ和潛在風(fēng)險(xiǎn)指數(shù)RIFig.2 Effects range-median quotient (ERMQ) values and potential ecological risk indices (RI) of heavy metals in sediments

    2.2 灤河岸帶濕地植被物種分布特征

    野外調(diào)查共識(shí)別維管束植物共計(jì)50科93屬219種,其中本地物種273種,外來(lái)物種54種。植物群落以被子植物為主,主要科系依次為菊科(42種)、禾本科(32種)、蓼科(12種)、豆科(11種)、薔薇科(11種)、莎草科(9種)、十字花科(9種)、藜科(6種),分別占總物種數(shù)的19.2%、14.6%、5.5%、5.5%、5.5%、4.1%、4.1%、2.7%,9個(gè)科系植物占所有物種的70.0%。植物區(qū)系優(yōu)勢(shì)科明顯,灤河上游為高原草甸,主要可見(jiàn)“披堿草(ElymusdahuricusTurcz.)+委陵菜(PotentillachinensisSer.)”、“藨草(ScirpustriqueterLinn.)+節(jié)節(jié)菜(Rotalaindica(Willd.) Koehne)+地榆(SanguisorbaocinalisL.)”和“野艾蒿(ArtemisialavandulaefoliaDC.)+藜(ChenopodiumserotinumLinn.)”等為優(yōu)勢(shì)種的群落,草甸中伴生有水楊梅(GeumaleppicumJacq.)、地榆(SanguisorbaocinalisL.)、萹蓄(PolygonumaviculareL.)等,其余雜類草群落成分繁雜,優(yōu)勢(shì)種不明顯。中游為山區(qū)植被,主要可見(jiàn)以“委陵菜(PotentillachinensisSer.)+無(wú)芒雀麥(BromusinermisLayss. Linn.)”、“扁稈藨草(ScirpusplaniculmisFr. Schmidt)+老芒麥(ElymussibiricusLinn.)”和“澤芹(SiumsuaveWalt.)+大刺兒菜(Cephalanoplossetosum(willd.) kitam.)+地筍(LycopuslucidusTurcz.)”等為優(yōu)勢(shì)種的群落,草甸中伴生有冬葵(MalvacrispaLinn.)、地膚(Kochiascoparia(Linn.) Schrad.)、蒼耳(XanthiumsibiricumPatrin ex Widder.)等。下游為平原植被,主要可見(jiàn)“附地菜(Trigonotispeduncularis(Trev.) Benth. ex Baker et Moore)+水芹(ArtemisialavandulaefoliaDC.)”、“野艾蒿(ArtemisialavandulaefoliaDC.)+附地菜(Trigonotispeduncularis(Trev.) Benth. ex Baker et Moore)”和“白茅(Imperatacylindrica(Linn.) Beauv.)+蘆葦(Phragmitesaustralis(Cav.) Trin. ex Steud.)”等為優(yōu)勢(shì)種的群落,草甸中伴生有圓葉牽牛(Pharbitispurpurea(L.) Voigt)、菖蒲(AcoruscalamusL.)、蘿藦(Metaplexisjaponica(Thunb.) Makino)等。

    綜合物種數(shù)、個(gè)體數(shù)、分布特性等植被調(diào)查信息,計(jì)算植物群落物種多樣性指數(shù)(物種豐度指數(shù)S、Shannon-Wiener指數(shù)H和Pielou均勻度指數(shù)J)定量反映植物群落物種多樣性特征,結(jié)果如圖3所示。物種豐富度指數(shù)S范圍為5.75—12.33,中游L9-西溝點(diǎn)位最高,下游L15-姜各莊點(diǎn)位最低,上、中、下游平均值分別為8.65、8.37、7.20;Shannon-Wiener指數(shù)H范圍為2.25—3.30,中游L9-西溝點(diǎn)位最高,下游L13-馬蘭莊鎮(zhèn)點(diǎn)位最低,上、中、下游平均值分別為3.13、2.92、2.33;Pielou均勻度指數(shù)J范圍為21.75—40.91,中游L9-西溝點(diǎn)位最高,下游L14-王家樓村點(diǎn)位最低,上、中、下游平均值分別為30.53、30.57、23.47??傮w來(lái)看,上游和中游群落物種多樣性指數(shù)波動(dòng)較小,二者高于下游,說(shuō)明上游和中游植被種類最多,群落復(fù)雜程度最高,且群落分布較為均勻。對(duì)不同河段植物群落物種多樣性指數(shù)的非參數(shù)檢驗(yàn)(Kruskal-Wallis法)結(jié)果表明,灤河干流上中下游Shannon-Wiener指數(shù)H和Pielou均勻度指數(shù)J具有顯著性差異(P<0.05)。

    圖3 灤河植物群落物種多樣性指數(shù) Fig.3 Vegetation species diversity index in Luanhe River

    2.3 灤河底棲動(dòng)物群落分布特征

    灤河水系共采集底棲動(dòng)物105種,屬于節(jié)肢動(dòng)物、環(huán)節(jié)動(dòng)物和軟體動(dòng)物等三大類,其中節(jié)肢動(dòng)物76種,占比72.4%,環(huán)節(jié)動(dòng)物18種,占比17.1%,軟體動(dòng)物10種,占比9.5%。優(yōu)勢(shì)類群為節(jié)肢動(dòng)物門的昆蟲(chóng)綱,共33科71屬73種,占全部種類的69.5%。EPT昆蟲(chóng)對(duì)棲息地生境變化十分敏感,通常生活在清潔河流水環(huán)境中,常作為河流水系清潔的指示生物[29]。本研究中EPT昆蟲(chóng)總物種數(shù)為17種,占調(diào)查發(fā)現(xiàn)的水生昆蟲(chóng)物種數(shù)的16.2%,蜉蝣目7科10屬,襀翅目2科2屬,毛翅目4科4屬。

    灤河水系底棲動(dòng)物生物量和密度如圖4所示,生物量范圍為0.31—67.03 g/m2,最低點(diǎn)位為L(zhǎng)6-外溝門,最高點(diǎn)位為L(zhǎng)14-遷西大橋,上、中、下游平均值分別為4.13、1.35、33.31 g/m2。密度范圍為3.85—414.00 個(gè)/m2,最低點(diǎn)位為L(zhǎng)2-白城子,最高點(diǎn)位為L(zhǎng)12-馬蘭莊鎮(zhèn),上、中、下游平均值分別為132.46、99.83、190.00 g/m2。研究團(tuán)隊(duì)基于蜉蝣目、襀翅目和毛翅目個(gè)體數(shù)之和的百分比(EPT%)、物種耐污性綜合指數(shù)(Biotic Index, BI)和收集者攝食類群的百分比(Collect-Gatherers%)三個(gè)指標(biāo),建立底棲動(dòng)物完整性指數(shù)(B-IBI)表征底棲動(dòng)物群落的結(jié)構(gòu)功能特征,綜合評(píng)價(jià)河流生態(tài)系統(tǒng)健康。灤河干流水系14個(gè)點(diǎn)位B-IBI均值為0.41,變化范圍為0.15—0.79,最大值為L(zhǎng)6-外溝門,最小值為L(zhǎng)2-白城子,上中下游點(diǎn)位平均值分別為0.56、0.29、0.28,中游B-IBI指數(shù)高于上游和下游,對(duì)B-IRHI指數(shù)進(jìn)行等級(jí)劃分,優(yōu)級(jí)(0.75—1.0)、良(0.50—0.75)、中(0.25—0.50)、差(0—0.25)點(diǎn)位數(shù)分別占比7.14%、42.86%、28.57%、28.57%,說(shuō)明灤河干流底棲動(dòng)物完整性總體良好,但上游和下游較差?;贙ruskal-Wallis法的不同河段底棲動(dòng)物完整性指數(shù)非參數(shù)檢驗(yàn)結(jié)果為,下游B-IBI指數(shù)低于上游和中游,具有顯著性差異(P<0.05)。

    圖4 灤河底棲動(dòng)物密度、生物量和底棲動(dòng)物完整性指數(shù) (B-IBI) 值Fig.4 The macroinvertebrate density, biomass and benthic macroinvertebrate-based index of biotic integrity (B-IBI) in Luanhe River

    2.4 生物群落對(duì)沉積物污染因子的響應(yīng)

    對(duì)灤河水系沉積物重金屬、岸帶濕地植物和大型底棲動(dòng)物特征信息進(jìn)行主成分分析(圖5),PC1解釋了總方差的65.43%,PC2解釋了方差的22.99%。對(duì)第一排序軸作用最大的生境因子是生物毒性效應(yīng)系數(shù)(ERMQ),沉積物重金屬污染程度與物種多樣性指數(shù)呈負(fù)相關(guān),與底棲動(dòng)物完整性指數(shù)(B-IBI)呈負(fù)相關(guān)。重金屬Zn、Cu、Cr、Pb作為主要環(huán)境因子,顯著影響了灤河水系植物群落物種多樣性和底棲動(dòng)物群落。

    圖5 沉積物特征和植被、底棲動(dòng)物因子的主成分分析Fig.5 Principal component analysis (PCA) of macroinvertebrate species, vegetation communities, and sediment indicators

    3 討論

    河流沉積物重金屬污染具有長(zhǎng)期潛在的危害, 重金屬能對(duì)底棲生物和植被產(chǎn)生毒性效應(yīng),進(jìn)而可能通過(guò)遷移進(jìn)入食物鏈危害人類健康。本研究中,灤河干流岸帶濕地沉積物無(wú)明顯重金屬污染,說(shuō)明該流域水系較為清潔,這與已有研究結(jié)果基本一致[28]。查閱文獻(xiàn),分析灤河沉積物歷史研究發(fā)現(xiàn),2008—2013年表層沉積物重金屬含量顯著升高,污染較為嚴(yán)重,說(shuō)明該歷史時(shí)間段水體吸納了大量污染物,表層沉積物受到嚴(yán)重污染。灤河作為京津冀地區(qū)重要水源地,近年來(lái)管理部門投入大量資金進(jìn)行專項(xiàng)整治,重金屬污染水平逐漸降低,但部分河段水質(zhì)和沉積物重金屬含量仍舊超標(biāo),本研究中重金屬污染程度上中下游空間差異較大,污染程度隨河流流向呈上升趨勢(shì),下游河段沉積物重金屬污染相比于上中游河段仍處于較高污染水平。變異系數(shù)能夠表征不同樣點(diǎn)間重金屬含量離散和均勻程度[30],本研究中灤河岸帶濕地沉積物重金屬含量變異程度較大,一定程度上說(shuō)明人類活動(dòng)可能對(duì)灤河水系沉積物重金屬分布造成影響。

    河流重金屬來(lái)源與周邊城市分布和人類活動(dòng)密切相關(guān)[31],灤河上游和中游分別位于高原和山區(qū),生態(tài)環(huán)境較為原始,人口稀疏且工農(nóng)業(yè)不發(fā)達(dá),但下游為人口密集工業(yè)發(fā)達(dá)的平原區(qū),周邊有大量城市和農(nóng)田,工農(nóng)業(yè)生產(chǎn)可能是灤河重金屬的主要來(lái)源。環(huán)境中Zn 主要來(lái)源于燃煤和鉛鋅礦業(yè),Cu 主要來(lái)源于金屬冶煉等,且Zn、Cu、Cd、Pb 和Cr 普遍存在于農(nóng)藥化肥中,灤河下游流經(jīng)河北省唐山市,該地區(qū)人口密集,冶金工業(yè)發(fā)達(dá),且冶金過(guò)程中會(huì)燃燒大量煤炭,結(jié)合本研究結(jié)果分析可知,灤河重金屬主要來(lái)源為下游平原區(qū)工農(nóng)業(yè)生產(chǎn),特別是下游冶金業(yè)和農(nóng)藥化肥的施用。

    植被物種多樣性和底棲動(dòng)物完整性是評(píng)價(jià)生態(tài)系統(tǒng)健康和穩(wěn)定的關(guān)鍵指標(biāo),對(duì)生境健康具有一定指示作用。本研究中灤河上中游植物群落物種多樣性較高,下游較低,這與已有研究結(jié)果一致[32- 33]。河岸帶濕地植物群落與濕地生態(tài)健康密切相關(guān),特別是水體中污染物例如重金屬是植物的重要脅迫因子。重金屬遷移進(jìn)入植物體內(nèi)后能夠?qū)е律锬p傷,光合作用效率下降等毒性效應(yīng),進(jìn)而導(dǎo)致部分耐受性差的植物物種退化、死亡,群落物種多樣性降低。類似的,底棲動(dòng)物也是水環(huán)境質(zhì)量的重要指示生物,對(duì)生境變化十分敏感,當(dāng)其棲息環(huán)境受到擾動(dòng)時(shí),其內(nèi)部群落特征會(huì)相應(yīng)發(fā)生改變。本研究中,灤河下游沉積物重金屬污染程度高于上中游,但下游底棲動(dòng)物生物量和密度高于上中游河段,分析原因可能為底棲動(dòng)物群落是一個(gè)復(fù)雜動(dòng)態(tài)整體,隨著環(huán)境因子發(fā)生變化,其內(nèi)部群落特征也會(huì)相應(yīng)發(fā)生改變[34],僅考慮群落整體生物量和密度并不能準(zhǔn)確反映底棲動(dòng)物群落對(duì)環(huán)境因子的響應(yīng)。例如,李晉鵬等[35]分析水庫(kù)底棲動(dòng)物群落對(duì)重金屬污染的響應(yīng),發(fā)現(xiàn)昆蟲(chóng)綱搖蚊科密度隨沉積物中重金屬鉻含量升高而升高。李立強(qiáng)等[36]對(duì)洞庭湖大型底棲動(dòng)物研究發(fā)現(xiàn)寡毛類底棲動(dòng)物較為耐受重金屬污染,其密度與重金屬含量呈極顯著正相關(guān)。Wang等[2]研究表明了以搖蚊為主腹足類昆蟲(chóng)(以搖蚊為主)與沉積物鉻、鉛濃度呈極顯著正相關(guān)。以往研究說(shuō)明隨著河流沉積物污染程度增加,部分耐污種底棲動(dòng)物呈增加趨勢(shì),大型底棲動(dòng)物是河流生態(tài)系統(tǒng)健康的敏感生物,近年來(lái)許多研究者建立底棲動(dòng)物完整性指數(shù)(B-IBI)綜合反映底棲動(dòng)物群落結(jié)構(gòu)功能特征,評(píng)價(jià)河流生態(tài)系統(tǒng)健康[23- 24]。本研究中灤河下游底棲動(dòng)物完整性較差,說(shuō)明沉積物中重金屬污染改變了河流底棲動(dòng)物群落結(jié)構(gòu),耐污種的增加使得群落整體密度和生物量增加,這與已有研究結(jié)果基本一致[37]。

    河流岸帶濕地植物和底棲動(dòng)物群落與沉積物質(zhì)量密切相關(guān),例如沉積物粒徑、污染物含量、營(yíng)養(yǎng)鹽水平等特征的改變都可能引起生物群落變化[2]。本研究中沉積物重金屬污染水平與植物群落物種多樣性、底棲動(dòng)物完整性呈負(fù)相關(guān)關(guān)系,說(shuō)明重金屬污染降低了水環(huán)境質(zhì)量,對(duì)植被和底棲動(dòng)物及濕地棲息地完整性造成脅迫效應(yīng)。

    4 結(jié)論與建議

    (1)灤河干流沉積物總體清潔,但具有明顯空間差異性,河流下游沉積物重金屬生物毒性危害和潛在生態(tài)風(fēng)險(xiǎn)高于上游和中游,下游工農(nóng)業(yè)是灤河沉積物重金屬主要污染源。

    (2)生物群落調(diào)查共識(shí)別維管束植物219種,大型無(wú)脊椎底棲動(dòng)物105種,下游植物群落物種多樣性最低,上游和中游群落物種多樣性指數(shù)波動(dòng)較小,二者高于下游,說(shuō)明上游和中游植被種類最多,群落復(fù)雜程度最高,且群落分布較為均勻。灤河干流上中下游Shannon-Wiener指數(shù)H和Pielou均勻度指數(shù)J具有顯著性差異(P<0.05)。

    (3)底棲動(dòng)物完整性下游最差,說(shuō)明上中游河段生境較好,下游植被和底棲動(dòng)物受到環(huán)境因子干擾。灤河干流底棲動(dòng)物完整性總體良好,但上游和下游較差?;贙ruskal-Wallis法的不同河段底棲動(dòng)物完整性指數(shù)非參數(shù)檢驗(yàn)結(jié)果為,下游B-IBI指數(shù)低于上游和中游,具有顯著性差異(P<0.05)。

    (4)主成分分析結(jié)果顯示,對(duì)第一排序軸作用最大生境因子是生物毒性效應(yīng)系數(shù)(ERMQ),沉積物重金屬污染程度與物種多樣性指數(shù)呈負(fù)相關(guān),與底棲動(dòng)物完整性指數(shù)(B-IBI)呈負(fù)相關(guān),說(shuō)明沉積物污染水平是影響河流岸帶濕地生態(tài)健康的重要環(huán)境因子。

    河流岸帶沉積物是影響濕地環(huán)境質(zhì)量和生態(tài)健康重要基質(zhì),灤河流域管理者應(yīng)進(jìn)一步加強(qiáng)下游污染源治理,特別是控制工業(yè)點(diǎn)源污染排放和農(nóng)業(yè)面源污染,充分保障灤河水源地生態(tài)安全。河流岸帶濕地是重要生物棲息地,未來(lái)研究中應(yīng)重點(diǎn)關(guān)注河流沉積物典型污染物釋放及其在生物群落內(nèi)部的遷移轉(zhuǎn)化過(guò)程,進(jìn)而揭示河流岸帶濕地棲息地關(guān)鍵生態(tài)功能過(guò)程與污染物特征的響應(yīng)關(guān)系。建議流域管理部門將岸帶濕地植被多樣性和底棲動(dòng)物完整性指數(shù)B-IBI納入河流水文-水環(huán)境水生態(tài)一體化監(jiān)測(cè)與風(fēng)險(xiǎn)評(píng)價(jià)指標(biāo)體系中。

    參考文獻(xiàn)(References):

    [1] Almeida J A, Diniz Y S, Marques S F G, Faine L A, Ribas B O, Burneiko R C, Novelli E L B. The use of the oxidative stress responses as biomarkers in Nile tilapia (Oreochromisniloticus) exposed to in vivo cadmium contamination. Environment International, 2002, 27(8): 673- 679.

    [2] Wang X X, Su P, Lin Q D, Song J X, Sun H T, Cheng D D, Wang S Q, Peng J L, Fu J X. Distribution, assessment and coupling relationship of heavy metals and macroinvertebrates in sediments of the Weihe River Basin. Sustainable Cities and Society, 2019, 50: 101665.

    [3] Mirzaei R, Ghorbani H, Hafezi Moghaddas N, Martín J A R. Ecological risk of heavy metal hotspots in topsoils in the Province of Golestan, Iran. Journal of Geochemical Exploration, 2014, 147: 268- 276.

    [4] 季鵬飛, 許海, 詹旭, 朱廣偉, 鄒偉, 朱夢(mèng)圓, 康麗娟. 長(zhǎng)江中下游湖泊水體氮磷比時(shí)空變化特征及其影響因素. 環(huán)境科學(xué), 2020, 41(09): 4030-4041.

    [5] Dias V, Vasseur C, Bonzom J M. Exposure ofChironomusripariuslarvae to uranium: effects on survival, development time, growth, and mouthpart deformities. Chemosphere, 2008, 71(3): 574- 581.

    [6] Gale S A, King C K, Hyne R V. Chronic sublethal sediment toxicity testing using the estuarine amphipod,Melitaplumulosa(Zeidler): evaluation using metal‐spiked and field‐contaminated sediments. Environmental Toxicology and Chemistry, 2006, 25(7): 1887- 1898.

    [7] Dabney B L, Clements W H, Williamson J L, Ranville J F. Influence of metal contamination and sediment deposition on benthic invertebrate colonization at the North Fork Clear Creek Superfund Site, Colorado, USA. Environmental Science & Technology, 2018, 52(12): 7072- 7080.

    [8] Costas N, Pardo I, Méndez-Fernández L, Martínez-Madrid M, Rodríguez P. Sensitivity of macroinvertebrate indicator taxa to metal gradients in mining areas in Northern Spain. Ecological Indicators, 2018, 93: 207- 218.

    [9] Bere T, Dalu T, Mwedzi T. Detecting the impact of heavy metal contaminated sediment on benthic macroinvertebrate communities in tropical streams. Science of the Total Environment, 2016, 572: 147- 156.

    [10] 師忱, 袁士保, 史常青, 趙廷寧. 灤河流域氣候變化與人類活動(dòng)對(duì)徑流的影響. 水土保持學(xué)報(bào), 2018, 32(2): 264- 269.

    [11] 李永麗, 劉靜玲. 灤河流域不同時(shí)空水環(huán)境重金屬污染健康風(fēng)險(xiǎn)評(píng)價(jià). 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào), 2009, 28(6): 1177- 1184.

    [12] 解瑩, 李敘勇, 王慧亮, 李文贊. 灤河流域上游地區(qū)主要河流水污染特征及評(píng)價(jià). 環(huán)境科學(xué)學(xué)報(bào), 2012, 32(3): 645- 653.

    [13] 王旭, 單保慶, 郭伊荇, 趙鈺, 唐文忠. 灤河干流表層沉積物中營(yíng)養(yǎng)元素和重金屬含量分布特征及其風(fēng)險(xiǎn)評(píng)價(jià). 環(huán)境工程學(xué)報(bào), 2017, 11(10): 5429- 5439.

    [14] 曹治國(guó), 劉靜玲, 欒蕓, 李永麗, 韓守亮, 劉寶艷, 王少明. 灤河流域多環(huán)芳烴的污染特征、風(fēng)險(xiǎn)評(píng)價(jià)與來(lái)源辨析. 環(huán)境科學(xué)學(xué)報(bào), 2010, 30(2): 246- 253.

    [15] 董滿宇, 王磊鑫, 李潔敏, 江源, 徐霞. 1960—2017年灤河流域降水集中度與集中期時(shí)空變化特征. 北京師范大學(xué)學(xué)報(bào): 自然科學(xué)版, 2019, 55(4): 468- 475.

    [16] 孟博, 劉靜玲, 李毅, 史璇. 北京市涼水河表層沉積物不同粒徑重金屬形態(tài)分布特征及生態(tài)風(fēng)險(xiǎn). 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào), 2015, 34(5): 964- 972.

    [17] MacDonald D D, Ingersoll C G, Berger T A. Development and evaluation of consensus-based sediment quality guidelines for freshwater ecosystems. Archives of Environmental Contamination and Toxicology, 2000, 39(1): 20- 31.

    隨著團(tuán)場(chǎng)小城鎮(zhèn)建設(shè)步伐的加快,母親居住了幾十年的平房被拆除了,搬遷進(jìn)了團(tuán)里安排的廉租房。廉租房在二樓,面積不大,也就60多個(gè)平方,很多東西都被母親扔掉了,但是那臺(tái)老舊的飛人牌縫紉機(jī),卻被母親搬到樓房上,放在了陽(yáng)臺(tái)上光線最充足的地方。

    [18] Roach A C. Assessment of metals in sediments from Lake Macquarie, New South Wales, Australia, using normalisation models and sediment quality guidelines. Marine Environmental Research, 2005, 59(5): 453- 472.

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

    [20] Huang Y, Chen Q Q, Deng M H, Japenga J, Li T Q, Yang X E, He Z L. Heavy metal pollution and health risk assessment of agricultural soils in a typical peri-urban area in southeast China. Journal of Environmental Management, 2018, 207: 159- 168.

    [21] 郭西亞, 高敏, 張杰, 張海濤, 朱金格, 鄧建才. 陽(yáng)澄湖沉積物重金屬空間分布及生物毒害特征. 中國(guó)環(huán)境科學(xué), 2019, 39(2): 802- 811.

    [22] 常虹, 孫海蓮, 劉亞紅, 邱曉, 石磊, 溫超. 東烏珠穆沁草甸草原不同退化程度草地植物群落結(jié)構(gòu)與多樣性研究. 草地學(xué)報(bào), 2020, 28(1): 184- 192.

    [23] Shi X, Liu J L, You X G, Bao K, Meng B, Bin C. Evaluation of river habitat integrity based on benthic macroinvertebrate-based multi-metric model. Ecological Modelling, 2017, 353: 63- 76.

    [24] Shi X, Liu J L, You X G, Bao K, Meng B. Shared effects of hydromorphological and physico-chemical factors on benthic macroinvertebrate integrity for substrate types. Ecological Indicators, 2019, 105: 406- 414.

    [26] 郝紅, 高博, 王健康, 周懷東, 陸瑾, 殷淑華, 朱成. 灤河流域沉積物中重金屬分布特征及風(fēng)險(xiǎn)評(píng)價(jià). 巖礦測(cè)試, 2012, 31(6): 1000- 1005.

    [27] 尚林源, 孫然好, 王趙明, 汲玉河, 陳利頂. 海河流域北部地區(qū)河流沉積物重金屬的生態(tài)風(fēng)險(xiǎn)評(píng)價(jià). 環(huán)境科學(xué), 2012, 33(2): 606- 611.

    [28] 蘇虹程, 單保慶, 唐文忠, 王立碩, 王衛(wèi)民. 海河流域典型清潔水系表層沉積物中重金屬總體污染水平研究. 環(huán)境科學(xué)學(xué)報(bào), 2015, 35(9): 2860- 2866.

    [29] Lenat D R, Penrose D L. History of the EPT taxa richness metric. Bulletin of the North American Benthological Society, 1996, 13(2): 305- 307.

    [30] Wang S P, Wang Y H, Zhang R J, Wang W T, Xu D Q, Guo J, Li P Y, Yu K F. Historical levels of heavy metals reconstructed from sedimentary record in the Hejiang River, located in a typical mining region of Southern China. Science of the Total Environment, 2015, 532: 645- 654.

    [31] Zhang G L, Bai J H, Xiao R, Zhao Q Q, Jia J, Cui B S, Liu X H. Heavy metal fractions and ecological risk assessment in sediments from urban, rural and reclamation-affected rivers of the Pearl River Estuary, China. Chemosphere, 2017, 184: 278- 288.

    [32] 周緒申, 齊向華, 吳筱, 孫峰, 孫凱艷, 王乙震, 張俊, 羅陽(yáng). 灤河干流浮游植物多樣性及污染狀況. 南水北調(diào)與水利科技, 2015, 13(3): 448- 452, 462- 462.

    [33] 吳佳寧, 王剛, 路獻(xiàn)品, 郝彩蓮, 張誠(chéng), 秦天玲. 灤河流域浮游生物與底棲動(dòng)物分布特征調(diào)查研究. 環(huán)境保護(hù)科學(xué), 2014, 40(6): 1- 6.

    [34] Peng F J, Pan C G, Zhang N S, ter Braak C J F, Salvito D, Selck H, Ying G G, Van Den Brink P J. Benthic invertebrate and microbial biodiversity in sub-tropical urban rivers: Correlations with environmental variables and emerging chemicals. Science of the Total Environment, 2020, 709: 136281.

    [35] 李晉鵬, 趙愛(ài)東, 董世魁. 瀾滄江梯級(jí)水電開(kāi)發(fā)對(duì)漫灣庫(kù)區(qū)大型底棲動(dòng)物群落及重金屬沉積的影響. 生態(tài)環(huán)境學(xué)報(bào), 2019, 28(1): 117- 127.

    [36] 李利強(qiáng), 王丑明, 張屹, 黃代中, 田琪. 洞庭湖大型底棲動(dòng)物與表層沉積物重金屬研究. 生態(tài)環(huán)境學(xué)報(bào), 2016, 25(2): 286- 291.

    [37] 田鵬, 錢昶, 林佳寧, 代陽(yáng), 張遠(yuǎn), 丁森, 高欣. 灤河流域大型底棲動(dòng)物生物完整性指數(shù)健康評(píng)價(jià). 中國(guó)環(huán)境監(jiān)測(cè), 2019, 35(4): 50- 58.

    猜你喜歡
    灤河沉積物河流
    晚更新世以來(lái)南黃海陸架沉積物源分析
    河北省灤河“一河一策”方案編制與實(shí)施評(píng)估
    河北水利(2022年4期)2022-05-17 05:42:42
    渤海油田某FPSO污水艙沉積物的分散處理
    海洋石油(2021年3期)2021-11-05 07:43:12
    水體表層沉積物對(duì)磷的吸收及釋放研究進(jìn)展
    讓灤河流域水量豐、水質(zhì)好、生態(tài)美——河北省人大常委會(huì)通過(guò)關(guān)于加強(qiáng)灤河流域水資源保護(hù)和管理的決定
    美麗河北之灤河
    公民與法治(2020年7期)2020-05-11 02:14:10
    河流
    流放自己的河流
    啊,我的灤河燕山
    當(dāng)河流遇見(jiàn)海
    国产一区二区三区综合在线观看 | 国产精品国产av在线观看| 成人亚洲精品av一区二区| 97精品久久久久久久久久精品| 亚洲精品影视一区二区三区av| 久久久久久伊人网av| 久久精品国产自在天天线| 国产午夜福利久久久久久| 国产欧美亚洲国产| 日日摸夜夜添夜夜添av毛片| 美女被艹到高潮喷水动态| 国产乱人偷精品视频| 久久人人爽人人片av| 51国产日韩欧美| 精品人妻熟女av久视频| 插阴视频在线观看视频| 国产免费一区二区三区四区乱码| 久久人人爽av亚洲精品天堂 | 亚洲精品aⅴ在线观看| 人体艺术视频欧美日本| 你懂的网址亚洲精品在线观看| 成人亚洲精品一区在线观看 | 麻豆精品久久久久久蜜桃| 性色avwww在线观看| videossex国产| 狂野欧美激情性xxxx在线观看| 欧美少妇被猛烈插入视频| 国产亚洲午夜精品一区二区久久 | 欧美高清性xxxxhd video| 女的被弄到高潮叫床怎么办| 国产精品久久久久久精品电影小说 | 亚洲丝袜综合中文字幕| 老女人水多毛片| 在线亚洲精品国产二区图片欧美 | 国产综合懂色| 在线亚洲精品国产二区图片欧美 | 内地一区二区视频在线| 免费看日本二区| 国产69精品久久久久777片| 人妻系列 视频| 狂野欧美白嫩少妇大欣赏| 美女被艹到高潮喷水动态| 亚洲精品一区蜜桃| 亚洲精品一区蜜桃| tube8黄色片| 亚洲美女搞黄在线观看| 成年av动漫网址| 国产精品蜜桃在线观看| 性色av一级| 男人狂女人下面高潮的视频| 成年免费大片在线观看| 国产精品国产三级国产专区5o| 欧美xxxx性猛交bbbb| 久久久久国产精品人妻一区二区| 精品久久国产蜜桃| 啦啦啦在线观看免费高清www| 丰满人妻一区二区三区视频av| 99久久精品热视频| 搡女人真爽免费视频火全软件| 一区二区av电影网| 亚洲国产精品专区欧美| 亚洲av成人精品一二三区| 九九久久精品国产亚洲av麻豆| 草草在线视频免费看| 亚洲欧美成人综合另类久久久| 国产伦精品一区二区三区视频9| 成人黄色视频免费在线看| 日本色播在线视频| 亚洲国产av新网站| 亚洲av男天堂| www.色视频.com| 亚洲精品乱久久久久久| 亚洲av免费在线观看| 成人国产av品久久久| 亚洲精品乱码久久久久久按摩| 熟女电影av网| 免费观看在线日韩| 午夜福利网站1000一区二区三区| 亚洲精品aⅴ在线观看| 久久久亚洲精品成人影院| 亚洲自拍偷在线| 亚洲精品日韩av片在线观看| 欧美日韩在线观看h| 亚洲熟女精品中文字幕| 哪个播放器可以免费观看大片| 免费高清在线观看视频在线观看| 久久热精品热| 国内精品美女久久久久久| 麻豆久久精品国产亚洲av| 日韩av不卡免费在线播放| 亚洲精品日韩在线中文字幕| 亚洲精品,欧美精品| 哪个播放器可以免费观看大片| 看非洲黑人一级黄片| 男女下面进入的视频免费午夜| 精品人妻视频免费看| 免费高清在线观看视频在线观看| 亚洲熟女精品中文字幕| 国产伦精品一区二区三区视频9| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 欧美另类一区| 超碰av人人做人人爽久久| 久久精品久久精品一区二区三区| 王馨瑶露胸无遮挡在线观看| 亚洲精品国产av成人精品| 好男人在线观看高清免费视频| 狂野欧美白嫩少妇大欣赏| 久久精品久久久久久久性| 欧美 日韩 精品 国产| 久久精品国产亚洲av天美| 久久久a久久爽久久v久久| 精品视频人人做人人爽| 免费看不卡的av| 久久鲁丝午夜福利片| 网址你懂的国产日韩在线| 精品熟女少妇av免费看| 成人毛片a级毛片在线播放| 久久精品人妻少妇| 2021天堂中文幕一二区在线观| 精华霜和精华液先用哪个| 国产精品av视频在线免费观看| 我的女老师完整版在线观看| 国产成人福利小说| 午夜亚洲福利在线播放| 一区二区三区乱码不卡18| 久久久久网色| 国产一级毛片在线| 亚洲国产成人一精品久久久| av黄色大香蕉| 青春草视频在线免费观看| 免费看不卡的av| 色综合色国产| 身体一侧抽搐| 精品人妻视频免费看| 亚洲伊人久久精品综合| 亚洲婷婷狠狠爱综合网| 插阴视频在线观看视频| 免费黄色在线免费观看| 亚洲精品国产av蜜桃| 久热久热在线精品观看| 在现免费观看毛片| 亚洲va在线va天堂va国产| 丝袜脚勾引网站| 国产精品熟女久久久久浪| 午夜老司机福利剧场| 最近中文字幕高清免费大全6| 国产精品女同一区二区软件| 综合色av麻豆| 69av精品久久久久久| 婷婷色综合www| 国产黄色视频一区二区在线观看| 亚洲精品中文字幕在线视频 | 亚洲成色77777| 亚洲国产高清在线一区二区三| freevideosex欧美| 少妇人妻 视频| 欧美国产精品一级二级三级 | 麻豆成人av视频| 亚洲国产日韩一区二区| 亚洲经典国产精华液单| av网站免费在线观看视频| 亚洲自拍偷在线| 一级毛片 在线播放| 久久久久久国产a免费观看| 视频中文字幕在线观看| 18禁裸乳无遮挡免费网站照片| 国产 精品1| 免费看av在线观看网站| 国产 一区 欧美 日韩| 日韩电影二区| 天天一区二区日本电影三级| 精品国产露脸久久av麻豆| av在线蜜桃| 成人黄色视频免费在线看| 九色成人免费人妻av| 欧美激情在线99| 亚洲av免费高清在线观看| 日韩一区二区三区影片| 免费看日本二区| 青春草视频在线免费观看| 色综合色国产| 精品人妻偷拍中文字幕| 少妇丰满av| 一区二区av电影网| 啦啦啦在线观看免费高清www| 91在线精品国自产拍蜜月| 在线观看免费高清a一片| 国产精品嫩草影院av在线观看| 日本黄色片子视频| 成人一区二区视频在线观看| 国产成人freesex在线| av在线蜜桃| 男女下面进入的视频免费午夜| 又粗又硬又长又爽又黄的视频| av在线播放精品| 国产精品99久久99久久久不卡 | 草草在线视频免费看| 国产视频首页在线观看| 一区二区三区乱码不卡18| 国产69精品久久久久777片| 国产高清三级在线| 国产精品人妻久久久久久| 国产精品久久久久久久电影| 亚洲成人精品中文字幕电影| 国产精品久久久久久精品电影| 日韩成人av中文字幕在线观看| 最近手机中文字幕大全| 欧美丝袜亚洲另类| 精品久久久噜噜| 精品酒店卫生间| 亚洲精品视频女| 美女国产视频在线观看| 深爱激情五月婷婷| 能在线免费看毛片的网站| 又爽又黄无遮挡网站| 午夜精品一区二区三区免费看| 精品酒店卫生间| 午夜福利在线在线| 18禁裸乳无遮挡动漫免费视频 | 亚洲天堂国产精品一区在线| 欧美老熟妇乱子伦牲交| 插逼视频在线观看| 欧美bdsm另类| 精品一区在线观看国产| 高清欧美精品videossex| 热re99久久精品国产66热6| 国产成人a区在线观看| 大香蕉久久网| 日日啪夜夜撸| 国产国拍精品亚洲av在线观看| 嫩草影院新地址| 丝袜脚勾引网站| av在线老鸭窝| 欧美激情久久久久久爽电影| 免费看a级黄色片| 久久久国产一区二区| 婷婷色综合www| 别揉我奶头 嗯啊视频| 午夜免费男女啪啪视频观看| 亚洲av日韩在线播放| 免费av毛片视频| 亚洲色图av天堂| 欧美日韩精品成人综合77777| 色视频在线一区二区三区| 亚洲精品成人av观看孕妇| 天堂网av新在线| 全区人妻精品视频| 国产一区二区在线观看日韩| 国产视频首页在线观看| 一本一本综合久久| 老女人水多毛片| 精品酒店卫生间| 欧美日韩在线观看h| 亚洲欧洲日产国产| 国产亚洲av嫩草精品影院| 伊人久久国产一区二区| 熟女av电影| 亚洲精品乱码久久久久久按摩| 欧美精品国产亚洲| 男男h啪啪无遮挡| 国产精品久久久久久久久免| 亚洲精品一二三| 午夜福利高清视频| 我的女老师完整版在线观看| 天天躁日日操中文字幕| 免费在线观看成人毛片| av专区在线播放| 制服丝袜香蕉在线| 色吧在线观看| 男男h啪啪无遮挡| 在线亚洲精品国产二区图片欧美 | 91狼人影院| 精品人妻偷拍中文字幕| 中文精品一卡2卡3卡4更新| 亚洲欧美日韩卡通动漫| 99热这里只有是精品在线观看| 黄色欧美视频在线观看| 你懂的网址亚洲精品在线观看| 久久久成人免费电影| 国产免费一区二区三区四区乱码| 成年av动漫网址| 亚洲无线观看免费| 亚洲激情五月婷婷啪啪| 免费av毛片视频| 91精品伊人久久大香线蕉| 女人被狂操c到高潮| 人妻 亚洲 视频| 精品国产露脸久久av麻豆| 视频区图区小说| 校园人妻丝袜中文字幕| 又爽又黄无遮挡网站| av在线亚洲专区| 日产精品乱码卡一卡2卡三| 97在线视频观看| av女优亚洲男人天堂| 国产色婷婷99| 亚洲欧美精品自产自拍| 国产精品一及| 插阴视频在线观看视频| 久久久a久久爽久久v久久| 国产有黄有色有爽视频| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 中文乱码字字幕精品一区二区三区| 日韩视频在线欧美| 午夜福利网站1000一区二区三区| 精品午夜福利在线看| 国内精品宾馆在线| 欧美一级a爱片免费观看看| 在线a可以看的网站| 97人妻精品一区二区三区麻豆| 色视频在线一区二区三区| 亚洲成人一二三区av| 欧美日本视频| 在线观看av片永久免费下载| 日韩av免费高清视频| 亚洲精品日本国产第一区| 黄色欧美视频在线观看| 天堂中文最新版在线下载 | 男男h啪啪无遮挡| 国产久久久一区二区三区| 色网站视频免费| 国产精品国产三级国产av玫瑰| 春色校园在线视频观看| 国产成人a区在线观看| 午夜福利网站1000一区二区三区| 一边亲一边摸免费视频| 热99国产精品久久久久久7| 国产精品麻豆人妻色哟哟久久| 肉色欧美久久久久久久蜜桃 | 99热这里只有精品一区| 看免费成人av毛片| 最新中文字幕久久久久| 王馨瑶露胸无遮挡在线观看| 国产大屁股一区二区在线视频| 国产精品久久久久久久电影| 亚洲精品乱久久久久久| 最近的中文字幕免费完整| 日韩在线高清观看一区二区三区| 蜜臀久久99精品久久宅男| 舔av片在线| 欧美丝袜亚洲另类| 亚洲欧美日韩东京热| 又大又黄又爽视频免费| av国产精品久久久久影院| 狠狠精品人妻久久久久久综合| 免费av不卡在线播放| 亚洲一级一片aⅴ在线观看| 亚洲四区av| 欧美潮喷喷水| 亚洲真实伦在线观看| 嫩草影院新地址| 国产在线一区二区三区精| 国产熟女欧美一区二区| 精品久久久久久电影网| 天天躁日日操中文字幕| 亚洲精品国产av蜜桃| 欧美变态另类bdsm刘玥| 久久99精品国语久久久| 激情 狠狠 欧美| 久久人人爽人人片av| 国产又色又爽无遮挡免| 国产永久视频网站| 激情五月婷婷亚洲| 亚洲精品第二区| kizo精华| 亚洲真实伦在线观看| 日韩成人av中文字幕在线观看| 国产精品久久久久久久久免| 久久精品熟女亚洲av麻豆精品| 中文天堂在线官网| 色哟哟·www| 亚洲人与动物交配视频| 久久这里有精品视频免费| 国产乱人视频| 99九九线精品视频在线观看视频| 日本黄色片子视频| 亚洲国产日韩一区二区| 最新中文字幕久久久久| 天天躁日日操中文字幕| 国产高清三级在线| 成人毛片a级毛片在线播放| 日本猛色少妇xxxxx猛交久久| 亚洲欧美一区二区三区国产| 欧美激情国产日韩精品一区| 少妇人妻 视频| 一级毛片我不卡| 又黄又爽又刺激的免费视频.| 欧美三级亚洲精品| 日韩av免费高清视频| 人妻系列 视频| 伦精品一区二区三区| 少妇丰满av| 国产亚洲av嫩草精品影院| 高清毛片免费看| 我的老师免费观看完整版| 麻豆国产97在线/欧美| 国产精品久久久久久久久免| 九九在线视频观看精品| 一级毛片电影观看| 国产欧美日韩一区二区三区在线 | av天堂中文字幕网| 一边亲一边摸免费视频| 亚洲精品中文字幕在线视频 | 99热网站在线观看| 亚洲国产精品成人久久小说| 日日摸夜夜添夜夜添av毛片| 国产精品av视频在线免费观看| 男人和女人高潮做爰伦理| 人人妻人人爽人人添夜夜欢视频 | 国产高清国产精品国产三级 | 卡戴珊不雅视频在线播放| 国产 一区 欧美 日韩| 国内少妇人妻偷人精品xxx网站| 精品久久国产蜜桃| 极品教师在线视频| 国产午夜福利久久久久久| 欧美性感艳星| 亚洲国产精品国产精品| 我的老师免费观看完整版| 国产综合精华液| 99久久人妻综合| 午夜日本视频在线| 69av精品久久久久久| 秋霞伦理黄片| 在线观看一区二区三区| 亚洲精品国产av成人精品| 亚洲精品一二三| 久久久成人免费电影| 久久久久国产精品人妻一区二区| 不卡视频在线观看欧美| 国产精品国产三级国产av玫瑰| 日韩一区二区三区影片| 亚洲经典国产精华液单| 99久久中文字幕三级久久日本| 熟女人妻精品中文字幕| 免费看不卡的av| 亚洲欧美成人综合另类久久久| 国产精品一区www在线观看| 欧美激情久久久久久爽电影| 国产精品国产三级国产av玫瑰| 麻豆精品久久久久久蜜桃| 人人妻人人澡人人爽人人夜夜| 精品久久国产蜜桃| 特大巨黑吊av在线直播| 中文天堂在线官网| 一级毛片久久久久久久久女| 新久久久久国产一级毛片| 中文资源天堂在线| 亚州av有码| 超碰97精品在线观看| 欧美日韩综合久久久久久| 丰满乱子伦码专区| 直男gayav资源| 亚洲精品久久久久久婷婷小说| 毛片一级片免费看久久久久| 久久久精品免费免费高清| 亚洲精品日韩在线中文字幕| 亚洲国产精品国产精品| 日韩欧美一区视频在线观看 | 六月丁香七月| 在线播放无遮挡| 七月丁香在线播放| 男人狂女人下面高潮的视频| 日本免费在线观看一区| 欧美日韩视频高清一区二区三区二| 免费黄频网站在线观看国产| 国产精品成人在线| 亚洲av免费高清在线观看| 精品酒店卫生间| 亚洲内射少妇av| 亚洲国产精品国产精品| 亚洲成人精品中文字幕电影| 亚洲国产色片| 黄色一级大片看看| 听说在线观看完整版免费高清| 国产成人精品一,二区| 日韩精品有码人妻一区| 久久久久性生活片| 插逼视频在线观看| 亚洲aⅴ乱码一区二区在线播放| 国产欧美日韩精品一区二区| 亚洲国产精品专区欧美| 国产视频内射| 国产女主播在线喷水免费视频网站| 观看美女的网站| 精品久久久久久久人妻蜜臀av| 午夜福利高清视频| 91在线精品国自产拍蜜月| 久久精品熟女亚洲av麻豆精品| 亚洲av国产av综合av卡| 69人妻影院| 国产淫语在线视频| 交换朋友夫妻互换小说| 最近中文字幕高清免费大全6| 可以在线观看毛片的网站| 少妇人妻久久综合中文| 免费黄色在线免费观看| 三级国产精品欧美在线观看| 午夜福利视频精品| 国产精品不卡视频一区二区| 国产av不卡久久| 亚洲图色成人| eeuss影院久久| 少妇 在线观看| 在线看a的网站| 美女内射精品一级片tv| 亚洲性久久影院| 97在线视频观看| 国产色婷婷99| 少妇人妻久久综合中文| 国产视频内射| 男女啪啪激烈高潮av片| 国产亚洲精品久久久com| 成人免费观看视频高清| 亚洲最大成人av| 久久这里有精品视频免费| 国产精品久久久久久精品电影| 美女高潮的动态| 亚洲精品久久久久久婷婷小说| 男人爽女人下面视频在线观看| 国产熟女欧美一区二区| 天美传媒精品一区二区| 人人妻人人看人人澡| 大码成人一级视频| 亚洲精品国产av蜜桃| 97人妻精品一区二区三区麻豆| 久久人人爽人人片av| 国产高潮美女av| 插阴视频在线观看视频| 亚洲成色77777| 亚洲国产精品成人综合色| 免费av观看视频| 高清av免费在线| 国产色爽女视频免费观看| 国产伦精品一区二区三区视频9| 又爽又黄无遮挡网站| 自拍偷自拍亚洲精品老妇| 国产精品99久久99久久久不卡 | 我要看日韩黄色一级片| 赤兔流量卡办理| 男女国产视频网站| 天堂网av新在线| 国产成人91sexporn| 亚洲久久久久久中文字幕| 亚洲精品日韩av片在线观看| 亚洲成人久久爱视频| 纵有疾风起免费观看全集完整版| 亚洲va在线va天堂va国产| 人体艺术视频欧美日本| 久久6这里有精品| 欧美人与善性xxx| 免费大片18禁| 午夜视频国产福利| 色哟哟·www| 中文字幕久久专区| 亚洲自偷自拍三级| 在线观看人妻少妇| 国产精品秋霞免费鲁丝片| 亚洲国产成人一精品久久久| 国产午夜精品久久久久久一区二区三区| 免费黄色在线免费观看| 91久久精品电影网| 精品一区二区三区视频在线| 我的老师免费观看完整版| 99久国产av精品国产电影| 欧美xxⅹ黑人| 日韩免费高清中文字幕av| av福利片在线观看| 欧美激情在线99| 黄色配什么色好看| av国产精品久久久久影院| 99re6热这里在线精品视频| 嫩草影院精品99| 小蜜桃在线观看免费完整版高清| 国产男女内射视频| 国产欧美日韩精品一区二区| 啦啦啦在线观看免费高清www| 真实男女啪啪啪动态图| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 大又大粗又爽又黄少妇毛片口| av免费在线看不卡| 国产日韩欧美亚洲二区| 七月丁香在线播放| 高清毛片免费看| 亚洲伊人久久精品综合| 特大巨黑吊av在线直播| 男女下面进入的视频免费午夜| 午夜精品国产一区二区电影 | 成人漫画全彩无遮挡| 另类亚洲欧美激情| 在线观看国产h片| 国产午夜精品一二区理论片| 久久久精品欧美日韩精品| 欧美变态另类bdsm刘玥| 国产亚洲最大av| 波多野结衣巨乳人妻| 久久ye,这里只有精品| 国产精品人妻久久久久久| 少妇裸体淫交视频免费看高清| 伦精品一区二区三区| 少妇裸体淫交视频免费看高清| 国产成人一区二区在线| 少妇裸体淫交视频免费看高清| 亚洲欧美日韩东京热| 一二三四中文在线观看免费高清| 亚洲经典国产精华液单| 美女内射精品一级片tv| 日本-黄色视频高清免费观看| 国产亚洲av嫩草精品影院| 午夜福利网站1000一区二区三区| 国产精品久久久久久精品电影小说 | 国产一级毛片在线| 日本与韩国留学比较| 色视频www国产| 美女主播在线视频|