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

    基于5種大型底棲動(dòng)物評(píng)價(jià)指數(shù)的河流生態(tài)健康評(píng)價(jià)

    2017-02-09 01:13:02毛建忠黎征武王旭濤鄧培雁
    水資源保護(hù) 2017年1期
    關(guān)鍵詞:東江樣點(diǎn)河流

    盛 蕭,毛建忠,曹 然,黎征武,王旭濤,鄧培雁

    (1.華南師范大學(xué)化學(xué)與環(huán)境學(xué)院,廣東 廣州 510006;2.云南省水文水資源局,云南 昆明 650106; 3.珠江流域水環(huán)境監(jiān)測(cè)中心,廣東 廣州 510611)

    ?

    基于5種大型底棲動(dòng)物評(píng)價(jià)指數(shù)的河流生態(tài)健康評(píng)價(jià)

    盛 蕭1,毛建忠2,曹 然1,黎征武1,王旭濤3,鄧培雁1

    (1.華南師范大學(xué)化學(xué)與環(huán)境學(xué)院,廣東 廣州 510006;2.云南省水文水資源局,云南 昆明 650106; 3.珠江流域水環(huán)境監(jiān)測(cè)中心,廣東 廣州 510611)

    大型底棲動(dòng)物;生物評(píng)價(jià)指數(shù);河流生態(tài)健康評(píng)價(jià);東江

    河流作為環(huán)境介質(zhì)和營(yíng)養(yǎng)物質(zhì)的載體,為水生生物提供棲息場(chǎng)所[1]。當(dāng)河流受到干擾時(shí),會(huì)使河流生態(tài)系統(tǒng)的結(jié)構(gòu)和功能發(fā)生變化,最終使河流健康受損。藻類、浮游動(dòng)物、魚類及大型底棲動(dòng)物等最常被用來監(jiān)測(cè)和評(píng)價(jià)河流健康狀況[2]。大型底棲動(dòng)物因種類繁多、分布廣泛、易于辨認(rèn)、能綜合反映物理、化學(xué)、生物等多種生態(tài)脅迫對(duì)水環(huán)境的累積效應(yīng)等特點(diǎn)[3],在河流健康評(píng)估中廣泛應(yīng)用,并取得了良好的效果。

    水質(zhì)生物學(xué)評(píng)價(jià)是指通過對(duì)水體中水生生物的調(diào)查或直接檢測(cè)來評(píng)價(jià)水體的生物學(xué)質(zhì)量[4]。德國(guó)科學(xué)家Kolkwitz于1902年首次提出并應(yīng)用生物評(píng)價(jià)指示河流有機(jī)污染后,許多學(xué)者逐漸認(rèn)識(shí)到生物評(píng)價(jià)在水質(zhì)的評(píng)價(jià)與管理中的不可替代性,并開始進(jìn)行大量的研究[5]。其中,基于大型底棲動(dòng)物評(píng)價(jià)指數(shù)的發(fā)展較為迅速,應(yīng)用較多,這些指數(shù)主要包括:①單一生物指數(shù)[6],依據(jù)指示物種存在的數(shù)量來定性判斷,如Goodnight-Whitley指數(shù)、Chandler計(jì)分系統(tǒng)、family biotic index(FBI)指數(shù)等;②生物多樣性指數(shù)[7],通過定量采樣,運(yùn)用統(tǒng)計(jì)學(xué)分析方法來建立,如Shannon-Wiener多樣性指數(shù)、Margalef多樣性指數(shù)等;③生物完整性指數(shù)[8],又稱多度量生物指數(shù),由Karr[8]于1981年提出,通過多種生物指數(shù)綜合反映水體的生物學(xué)狀況,從而評(píng)估河流水質(zhì)。

    Goodnight指數(shù)是最早的基于大型底棲無(wú)脊椎動(dòng)物的指數(shù),Wright等于1933年通過計(jì)算寡毛類的密度來反映水體污染程度,并在此基礎(chǔ)上發(fā)展了Goodnight指數(shù)和Goodnight-Whitley修正指數(shù)(GBI指數(shù))[5]。average score per taxon(ASPT)指數(shù)是在biological monitoring working party(BMWP)指數(shù)的基礎(chǔ)上發(fā)展而來,考慮到不同類群的耐污值可能會(huì)受總物種數(shù)的影響,Armitage等[9]提出了ASPT指數(shù),以降低總類群中出現(xiàn)的偶見種對(duì)評(píng)價(jià)結(jié)果的影響,該指數(shù)在應(yīng)用過程中僅要求將底棲動(dòng)物鑒定到科,既減少了工作量,又減少了鑒定錯(cuò)誤所帶來的誤差。FBI指數(shù)由美國(guó)學(xué)者Hilsenhoff[10]于1988年提出,為降低物種鑒定的難度和節(jié)省時(shí)間,實(shí)現(xiàn)河流健康的快速評(píng)價(jià),他在其所建立的Hilsenhoff biological index(HBI)指數(shù)的基礎(chǔ)上,提出了FBI指數(shù),有效地推動(dòng)了biological index(BI)指數(shù)在美國(guó)的廣泛應(yīng)用。

    生物多樣性指數(shù)是目前常用作水質(zhì)生物學(xué)評(píng)價(jià)的一種指數(shù),主要是通過群落中各物種的組成狀況來反映水體污染對(duì)群落造成的影響,從而劃分水質(zhì)級(jí)別[11],由于其在評(píng)價(jià)過程中存在的局限性,一般不單獨(dú)用作水質(zhì)評(píng)價(jià)。主要包括Shannon-Wiener多樣性指數(shù)、Simpson多樣性指數(shù)、Pielous均勻度指數(shù)和Margalef豐富度指數(shù)4種指數(shù),而在國(guó)際上使用較多的是Shannon-Wiener多樣性指數(shù)。

    benthic index of biotic integrity(B-IBI)指數(shù)[12]是在一系列參數(shù)中選擇能夠反映水環(huán)境中大型底棲動(dòng)物的種類組成、營(yíng)養(yǎng)關(guān)系、豐度和生長(zhǎng)狀況的生物參數(shù),相對(duì)于BI指數(shù),能夠傳達(dá)較多的河流生態(tài)完整性的信息。該方法在國(guó)外已得到較為廣泛的應(yīng)用。在國(guó)內(nèi),楊蓮芳等[13]于1992年首次將B-IBI指數(shù)引進(jìn),隨后,王備新等[14-15]應(yīng)用B-IBI指數(shù)分別評(píng)價(jià)了安徽黃山溪流和香溪河的河流健康狀況。

    東江是珠江的3大水系之一,水資源豐富,是廣東省諸多城市的主要生產(chǎn)和生活水源,同時(shí)擔(dān)負(fù)著香港地區(qū)的用水供需任務(wù)。然而近年來,隨著東江流域周邊城市的經(jīng)濟(jì)快速發(fā)展,各種工業(yè)在東江沿岸相繼落戶,對(duì)東江的水環(huán)境造成巨大威脅[16]。因此,本研究選用不同類型的大型底棲動(dòng)物評(píng)價(jià)指數(shù)對(duì)東江進(jìn)行評(píng)價(jià),并進(jìn)行對(duì)比,根據(jù)評(píng)價(jià)結(jié)果綜合分析不同生物評(píng)價(jià)指數(shù)在東江河流生態(tài)健康評(píng)估中的適用性,以期為東江流域綜合管理策略的制定與實(shí)踐提供理論基礎(chǔ)和科學(xué)依據(jù)。

    1 研究方法

    1.1 研究區(qū)域概況及采樣點(diǎn)分布

    東江發(fā)源于江西省尋烏縣,向西南流入廣東境內(nèi),經(jīng)龍川、惠州、東莞等縣市流入珠江,在獅子洋出虎門入海,干流全長(zhǎng)562 km,流域面積達(dá)35 340 km2。東江流域?qū)僦心蟻啛釒Ъ撅L(fēng)氣候,多年平均氣溫為20.4°C,多年平均降雨量1 750 mm[17]。全流域共設(shè)置24個(gè)采樣點(diǎn),包括干流及其所有一級(jí)支流,采樣點(diǎn)分布見圖1。

    干流樣點(diǎn):1—南農(nóng)村,2—車田水,6—柏埔河,17—寶華寺,20—潼湖農(nóng)場(chǎng)。支流樣點(diǎn):魚潭江 3—三多村;浰江 4—林寨鎮(zhèn),5—東水鎮(zhèn);康禾河 7—若壩大橋;秋香江 8—瓦溪鎮(zhèn),9—鳳安鎮(zhèn);石壩河 10—石壩橋;公莊河 15—公莊鎮(zhèn);赤水河 16—柏塘鎮(zhèn);新豐江 14—新豐縣福水陂;增江 11—正果大橋,12—黃竹瀝村,13—隔水橋,22—雁塔大橋,24—廟潭大橋;淡水河 18—三和鎮(zhèn);石馬河 19—東莞雍景花園;沙河 21—龍華橋;西福河 23—神崗橋(I)圖1 東江采樣點(diǎn)分布示意圖

    1.2 樣品的采集、處理與分析

    采樣安排在枯水期,分多組,于3日內(nèi)完成采樣。利用直徑30 cm、孔徑500 μm的尼龍紗D形網(wǎng)采集大型底棲動(dòng)物樣品。每個(gè)采樣點(diǎn)在100 m長(zhǎng)的范圍內(nèi)(水深小于1.5 m),采集長(zhǎng)度為3~10 m,面積為0.9~3.0 m2。每個(gè)樣點(diǎn)采集3個(gè)平行樣,以減少誤差。將D形網(wǎng)所采底泥倒入水桶,用水柔和攪動(dòng)桶內(nèi)底質(zhì),用孔徑為300 μm的篩網(wǎng)過濾,重復(fù)多次,直至目測(cè)無(wú)大型底棲動(dòng)物漂浮,再將底質(zhì)倒入白瓷盤中進(jìn)行仔細(xì)挑揀。隨后裝入廣口瓶中,用70%乙醇固定保存,帶回實(shí)驗(yàn)室作進(jìn)一步分揀。在體視顯微鏡下鑒定、分類,將每個(gè)樣點(diǎn)采集的大型底棲動(dòng)物按不同種類準(zhǔn)確地統(tǒng)計(jì)個(gè)體數(shù)[18]。種類鑒定大部分至科或?qū)?少數(shù)鑒定到種[19]。

    1.3 水樣采集與理化指標(biāo)測(cè)定

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

    本文選用了大型底棲動(dòng)物單一生物指數(shù)中的GBI、ASPT和FBI指數(shù),生物多樣性指數(shù)Shannon-Wiener指數(shù),生物完整性指數(shù)B-IBI指數(shù)來比較不同指數(shù)對(duì)東江健康評(píng)估的區(qū)別。GBI、ASPT、FBI和Shannon-Wiener多樣性指數(shù)計(jì)算方法[21]如下:

    (1)

    (2)

    (3)

    (4)

    式中:IGBI為GBI指數(shù)值;IASTP為ASTP指數(shù)值;IFBI為FBI指數(shù)值;H′為Shannon-Wiener多樣性指數(shù)值;N為樣品中大型底棲動(dòng)物個(gè)體總數(shù);Noil為樣品中寡毛類個(gè)體總數(shù);ti是每個(gè)類群的敏感值;n為總的科級(jí)分類單元數(shù);gi是大型底棲動(dòng)物第i個(gè)科級(jí)分類單元的個(gè)體數(shù);mi是科級(jí)分類單元i的耐污值;Ni為第i個(gè)物種的個(gè)體數(shù);S為總種類數(shù)。

    該區(qū)域B-IBI指數(shù)值(IB-IBI)的計(jì)算參見相關(guān)研究[22]。其中,參照點(diǎn)和受損點(diǎn)主要根據(jù)水質(zhì)等級(jí)、底質(zhì)類型、棲境復(fù)雜性、土地利用狀況等來確定,候選參數(shù)主要是依據(jù)東江采樣的實(shí)際情況選取了30個(gè)分屬4種類型(包括群落結(jié)構(gòu)與功能、功能攝食類群、干擾耐受和棲境)的生物參數(shù)。然后通過分布范圍檢驗(yàn)、敏感性分析及相關(guān)性分析對(duì)參數(shù)進(jìn)行逐步篩選,最終選取了6個(gè)參數(shù)(包括總分類單元、EPT分類單元、EPT相對(duì)豐度、前3位優(yōu)勢(shì)分類單元相對(duì)豐度、刮食者相對(duì)豐度和集食者相對(duì)豐度)組成B-IBI指數(shù)。5種類型的大型底棲動(dòng)物生物評(píng)價(jià)指數(shù)健康等級(jí)的最終劃分標(biāo)準(zhǔn)[23]見表2。

    表2 生物評(píng)價(jià)指數(shù)等級(jí)劃分標(biāo)準(zhǔn)

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

    大型底棲動(dòng)物的數(shù)據(jù)處理利用軟件EXCEL2010進(jìn)行,不同生物評(píng)價(jià)指數(shù)的正態(tài)分布檢驗(yàn)、相關(guān)性分析及與環(huán)境因子間的線性回歸分析均利用軟件SPSS 19.0完成。

    2 結(jié)果與分析

    2.1 不同生物評(píng)價(jià)指數(shù)對(duì)東江的水質(zhì)評(píng)估

    利用5種生物評(píng)價(jià)指數(shù)綜合分析東江水質(zhì)狀況,評(píng)價(jià)結(jié)果見表3。結(jié)果表明,對(duì)于24個(gè)樣點(diǎn)中的具體樣點(diǎn),根據(jù)不同生物評(píng)價(jià)指數(shù)所評(píng)價(jià)出的水質(zhì)等級(jí)不完全相同。其中,部分樣點(diǎn),5種生物指數(shù)的評(píng)價(jià)結(jié)果相近,如干流中的樣點(diǎn)6、支流中的樣點(diǎn)9等;然而有不少樣點(diǎn),5種指數(shù)的評(píng)價(jià)等級(jí)相差甚遠(yuǎn),以支流中的樣點(diǎn)15為例,FBI和B-IBI指數(shù)都將其評(píng)為水質(zhì)“較差”,而GBI、ASPT和Shannon-Wiener多樣性指數(shù)則分別將其評(píng)價(jià)為“健康”、“極差”和“一般”,由此可見,僅依靠單一的生物評(píng)價(jià)指數(shù),難以準(zhǔn)確地表征和反映東江流域內(nèi)具體監(jiān)測(cè)樣點(diǎn)的健康狀況。

    綜合5種指數(shù)的評(píng)價(jià)結(jié)果,從整體上來看,在東江流域,干流水質(zhì)上游相比下游較好,上游樣點(diǎn)除樣點(diǎn)2(車田水)的ASPT指數(shù)和B-IBI指數(shù)評(píng)價(jià)為水質(zhì)“較差”外,其余樣點(diǎn)評(píng)價(jià)結(jié)果都好于“一般”等級(jí),調(diào)查顯示,樣點(diǎn)2周邊的人類活動(dòng)較強(qiáng),這可能是造成其水質(zhì)差的一個(gè)因素。支流中,水質(zhì)相對(duì)較好的為魚潭江、浰江、康禾河、秋香江、增江上游河段;水質(zhì)一般的為新豐江、增江中游段、西福河;而水質(zhì)較差為沙河、公莊河;水質(zhì)極差即受到嚴(yán)重污染的為淡水河、石馬河、增江下游河段??傮w而言,東江流域上游河段的水質(zhì)優(yōu)于下游,河流水質(zhì)整體處于中等狀態(tài)。

    表3 不同生物評(píng)價(jià)指數(shù)對(duì)東江水質(zhì)的評(píng)價(jià)結(jié)果

    圖2 不同生物評(píng)價(jià)指數(shù)結(jié)果分布

    圖2為東江流域不同生物評(píng)價(jià)指數(shù)結(jié)果分布情況,由圖2可見,5種生物評(píng)價(jià)指數(shù)的健康評(píng)價(jià)等級(jí)在總體上存在較大差異。利用GBI指數(shù)評(píng)價(jià)河流健康等級(jí)的結(jié)果相對(duì)偏高,全流域中達(dá)到“健康”的樣點(diǎn)數(shù)占總樣點(diǎn)數(shù)的54.2%,其次為“一般”,占16.7%;ASPT指數(shù)評(píng)價(jià)的“健康”樣點(diǎn)占總樣點(diǎn)數(shù)的33.3%,其次為“極差”,占29.2%,ASPT指數(shù)表現(xiàn)的趨勢(shì)與GBI指數(shù)較為相似,唯一的區(qū)別是,ASPT指數(shù)在最后表現(xiàn)為上升,GBI指數(shù)則呈現(xiàn)下降;FBI指數(shù)和B-IBI指數(shù)的曲線趨勢(shì)則表現(xiàn)極其相似,均呈現(xiàn)為先下降,后上升,再下降的趨勢(shì),達(dá)到“一般”和“較差”的樣點(diǎn)數(shù)均占總樣點(diǎn)數(shù)的29.2%,其次為“健康”,分別占20.8%和16.7%。與以上4種方法評(píng)價(jià)結(jié)果不盡相同的是Shannon-Wiener多樣性指數(shù),其在總體上呈現(xiàn)為倒U型,其中,水質(zhì)等級(jí)達(dá)到“一般”的樣點(diǎn)數(shù)為最多,占58.3%,無(wú)達(dá)到“健康”的樣點(diǎn),其他評(píng)價(jià)等級(jí)的樣點(diǎn)數(shù)百分比幾乎持平(12.5%~16.7%)??傮w而言,B-IBI指數(shù)的評(píng)價(jià)結(jié)果在各評(píng)價(jià)等級(jí)間的比例最為平均(12.5%~29.2%),而其他指數(shù)在各等級(jí)間的比例則有較大的起伏。

    2.2 不同生物評(píng)價(jià)指數(shù)的相關(guān)性

    表4 不同生物評(píng)價(jià)指數(shù)的K-S檢驗(yàn)結(jié)果

    表5 不同生物評(píng)價(jià)指數(shù)間的Pearson相關(guān)性分析結(jié)果

    圖3 不同生物評(píng)價(jià)指數(shù)與環(huán)境因子的回歸分析

    2.3 不同生物評(píng)價(jià)指數(shù)與環(huán)境因子的關(guān)系

    3 討 論

    3.1 不同生物評(píng)價(jià)指數(shù)的比較

    對(duì)于干流,5種指數(shù)的評(píng)價(jià)結(jié)果差異不大,除靠近源頭的樣點(diǎn)2外,干流中上游河段的水質(zhì)總體較好,下游河段較差。這與王兆印等[24]的研究結(jié)果一致,其從東江上游到河口選取12個(gè)樣點(diǎn)進(jìn)行底棲動(dòng)物調(diào)查,發(fā)現(xiàn)東江流域底棲動(dòng)物生物多樣性在上游和中游保持在較高水平,而到下游則迅速降低到0。

    而對(duì)于支流,則呈現(xiàn)不同程度的差異。GBI指數(shù)的評(píng)價(jià)結(jié)果表明,東江的大多數(shù)支流,其下游河段的水質(zhì)相對(duì)上游較好,這與其他指數(shù)評(píng)價(jià)結(jié)果相反。以增江段為例,樣點(diǎn)13屬于增江上游,其水質(zhì)GBI評(píng)價(jià)結(jié)果為“較差”,而在對(duì)該區(qū)域的實(shí)際調(diào)查時(shí),此樣點(diǎn)周邊的棲境較為復(fù)雜,植被多樣性較高,水體透明度高,由此說明,僅將寡毛類個(gè)體數(shù)存在的比例用來評(píng)價(jià)河流水質(zhì)狀況,其精確性有所限制。王博等[25]在對(duì)東江干流的43個(gè)點(diǎn)位進(jìn)行底棲動(dòng)物調(diào)查,并結(jié)合多種指數(shù)對(duì)東江水質(zhì)進(jìn)行評(píng)價(jià)時(shí),發(fā)現(xiàn)GBI指數(shù)的評(píng)價(jià)結(jié)果較適合東江的實(shí)際情況,但從本文的研究結(jié)果來看,該指數(shù)并不完全適用于東江支流的評(píng)價(jià),也即不能完全反映東江水質(zhì)的整體實(shí)際情況。對(duì)于樣點(diǎn)3、5和18,ASPT指數(shù)的評(píng)價(jià)結(jié)果與其他4種指數(shù)相差較大,表現(xiàn)為相反的情況;對(duì)于樣點(diǎn)12,Shannon-Wiener指數(shù)與ASPT指數(shù)評(píng)價(jià)結(jié)果一致,而其他3種指數(shù)評(píng)價(jià)一致,造成這一系列差異的原因可能是多種因素的綜合影響。

    總體而言,5種指數(shù)對(duì)東江支流的評(píng)價(jià)結(jié)果在總體上一致:東江上游支流的水質(zhì)相對(duì)下游較好。其中,除GBI外,其他指數(shù)的評(píng)價(jià)結(jié)果表明增江上游水質(zhì)較下游好。何琦等[26]對(duì)東江支流增江進(jìn)行底棲硅藻的采樣調(diào)查,發(fā)現(xiàn)敏感種在增江的上游出現(xiàn)頻率較高,表明增江中上游河段水質(zhì)較清潔,而位于下游河段的2采樣點(diǎn),以耐污種占優(yōu)勢(shì),表明增江下游河段污染嚴(yán)重。由此可見,需結(jié)合多種生物指數(shù)對(duì)東江河流水質(zhì)進(jìn)行綜合評(píng)價(jià),才能提高評(píng)價(jià)結(jié)果的可靠性。

    從對(duì)東江河流樣點(diǎn)的評(píng)估等級(jí)結(jié)果可知,5種生物指數(shù)也呈現(xiàn)一定的差異。造成所有這一系列差異的主要原因可能有:①不同的評(píng)價(jià)指數(shù)存在自身的局限性,如Patrick[27]認(rèn)為多樣性指數(shù)作出的水質(zhì)評(píng)判結(jié)果相當(dāng)模糊,因?yàn)椴煌乃w污染類型對(duì)物種多樣性的影響程度不一致。房英春等[28]認(rèn)為多樣性指數(shù)評(píng)價(jià)中沒考慮生物體本身耐污性的差異,尤其是耐污種的敏感問題,可能會(huì)造成比實(shí)際值偏高的水質(zhì)評(píng)價(jià)結(jié)果。②底棲動(dòng)物的敏感值和耐污值的確定不明確,我國(guó)在對(duì)基于耐污值和敏感值的指數(shù)計(jì)算時(shí)主要是參照國(guó)外和國(guó)內(nèi)的相關(guān)研究,盡管目前王備新等[29]建立了一套適用我國(guó)華東地區(qū)底棲動(dòng)物的耐污值,但是否適用于其他區(qū)域尚待考證。③不同生物評(píng)價(jià)標(biāo)準(zhǔn)中等級(jí)范圍的劃分對(duì)結(jié)果也會(huì)造成一定的影響。如B-IBI指數(shù)最終的評(píng)價(jià)標(biāo)準(zhǔn)具有一定的主觀性,研究者在對(duì)該指數(shù)的構(gòu)建過程中,如何確定參照點(diǎn)和候選參數(shù),以及最終使用何種方法計(jì)算綜合參數(shù)都可能會(huì)對(duì)評(píng)價(jià)結(jié)果產(chǎn)生影響。

    3.2 不同生物評(píng)價(jià)指數(shù)在東江的適用性

    河流的健康狀況因地理位置、地貌特征、棲境、生物狀況、污染類型等多種因素的綜合作用存在差異,因此在選取不同生物評(píng)價(jià)指數(shù)評(píng)價(jià)河流健康狀況時(shí)應(yīng)慎重考慮生物評(píng)價(jià)指數(shù)在該河流中的適用性。

    BOD和DO均是與耗氧相關(guān)的水化指標(biāo)。BOD一般可用來表征有機(jī)污染的情況,多樣性指數(shù)可用來反映水體有機(jī)污染對(duì)群落結(jié)構(gòu)造成的影響[21],當(dāng)水體受到有機(jī)污染、無(wú)機(jī)還原性物質(zhì)污染時(shí),DO的質(zhì)量濃度就會(huì)降低[32]。在本研究中,Shannon-Wiener多樣性指數(shù)與BOD呈現(xiàn)較顯著的負(fù)相關(guān)關(guān)系;B-IBI指數(shù)與DO呈現(xiàn)出明顯的正相關(guān)關(guān)系,隨著DO質(zhì)量濃度的增加,B-IBI指數(shù)呈現(xiàn)上升的趨勢(shì)。由此可見,Shannon-Wiener多樣性指數(shù)和B-IBI指數(shù)對(duì)水體中的有機(jī)污染有較好的指示作用。

    電導(dǎo)率表征水體中溶解的總離子量,可綜合反映流域內(nèi)土地利用對(duì)河流生態(tài)系統(tǒng)的影響,通常電導(dǎo)率越高,人類活動(dòng)對(duì)河流的干擾作用越大[33]。渠曉東等[22]在運(yùn)用標(biāo)準(zhǔn)化方法篩選參照點(diǎn)構(gòu)建B-IBI指數(shù)時(shí),發(fā)現(xiàn)B-IBI指數(shù)與電導(dǎo)率存在顯著的對(duì)數(shù)曲線關(guān)系,隨著電導(dǎo)率的升高,B-IBI值逐漸下降,電導(dǎo)率高于1 000 μS/cm后,B-IBI指數(shù)維持在一個(gè)較低的水平。本研究也有一致的結(jié)果,B-IBI指數(shù)與電導(dǎo)率存在顯著的負(fù)相關(guān)關(guān)系,不同的是B-IBI值在電導(dǎo)率達(dá)到600 μS/cm后才緩慢下降,這可能與在B-IBI指數(shù)構(gòu)建時(shí)研究者所選用的方法不同而產(chǎn)生的差異相關(guān)。同時(shí),本研究中發(fā)現(xiàn)Shannon-Wiener多樣性指數(shù)與電導(dǎo)率也有顯著的負(fù)相關(guān)關(guān)系,因此,可認(rèn)為這2種指數(shù)能夠反映土地利用對(duì)東江流域水質(zhì)的影響。

    4 結(jié) 論

    選取GBI、ASPT、FBI、Shannon-Wiener多樣性和B-IBI指數(shù)對(duì)東江流域的主要干支流進(jìn)行綜合評(píng)價(jià),結(jié)果表明:

    a. 東江流域的整體健康狀態(tài)處于中等水平,河流上游水質(zhì)多達(dá)到“良好”,中游水質(zhì)多為“一般”,下游水質(zhì)則為“較差”或“極差”。

    b. 5種生物評(píng)價(jià)指數(shù)大部分都具有顯著的相關(guān)性,其中,FBI指數(shù)與其他4種指數(shù)都達(dá)到了顯著水平。但由于多種因素的綜合影響,不同指數(shù)對(duì)各樣點(diǎn)的評(píng)價(jià)等級(jí)存在一定差異。

    [ 1 ] 張晶, 董哲仁, 孫東亞, 等. 基于主導(dǎo)生態(tài)功能分區(qū)的河流健康評(píng)價(jià)全指標(biāo)體系[J]. 水利學(xué)報(bào), 2010, 41(8): 883-892. (ZHANG Jing, DONG Zheren, SUN Dongya, et al. Complete river health assessment index system based on eco-regional method according to dominant ecological functions[J]. Journal of Hydraulic Engineering, 2010, 41(8): 883-892. (in Chinese))

    [ 2 ] 戴紀(jì)翠, 倪晉仁. 底棲動(dòng)物在水生態(tài)系統(tǒng)健康評(píng)價(jià)中的作用分析[J].生態(tài)環(huán)境, 2008, 17(6): 2107-2111. (DAI Jicui, NI Jinren. Roles of benthos in the aquatic ecosystem health assessment[J]. Ecology and Ecosystem, 2008, 17(6): 2107-2111.(in Chinese))

    [ 3 ] AURA C M, RABUBU P O, HERRMANN J. A preliminary macroinvertebrate index of biotic integrity foe bioassessment of the Kipkaren and Sosiani Rivers, Nzoia River Basin, Kenya[J]. Lake and Reservoirs: Research and Management, 2010, 15(2): 119-128.

    [ 4 ] 王備新. 大型底棲無(wú)脊椎動(dòng)物水質(zhì)生物評(píng)價(jià)研究[D]. 南京: 南京農(nóng)業(yè)大學(xué), 2003.

    [ 5 ] HERING D, MOOG O, SANDIN L, et al. Overview and application of the AQEM assessment system[J]. Hydrobiologia, 2004, 516(1): 1-20.

    [ 6 ] 李鎮(zhèn), 張巖, 袁建平, 等. 大型底棲無(wú)脊椎動(dòng)物在河流健康評(píng)價(jià)中的發(fā)展趨勢(shì)[J]. 南水北調(diào)與水利科技, 2011, 9(4): 96-101. (LI Zhen, ZHANG Yan, YUAN Jianping, et al. Review of the role of benthic macro-invertebrates as indicators on river health assessment[J]. South-to-North Water Diversion and Water Science & Technology, 2011, 9(4): 96-101.(in Chinese))

    [ 7 ] NAQINEZHADA A, HAMZEH B, ATTARB F. Vegetation-Environment relationships in the Alder Wood Communities of Caspian Lowlands, N. Iran (toward an ecological classification) Flora [J]. Flora-Morphology, Distribution, Functional Ecology of Plants, 2008, 203(7): 567-577.

    [ 8 ] KARR J R. Assessment of biotic integrity using fish communities[J]. Fisheries, 1981, 6(6): 21-27.

    [ 9 ] ARMITAGE P D, MOSS D, WRIGNT J F, et al. The performance of a new biological water quality score system based on macroinvertebrates over a wide range of unpolluted running-water sites[J]. Water Research, 1983, 17(3): 333-347.

    [10] HILSENHOFF W L. Rapid field assessment of organic pollution with a family-level biotic index[J]. Journal of the North American Benthological Society, 1988, 7(1): 65-68.

    [11] 吳東浩, 王備新, 張?jiān)? 等. 底棲動(dòng)物生物指數(shù)水質(zhì)評(píng)價(jià)進(jìn)展及在中國(guó)的應(yīng)用前景[J]. 南京農(nóng)業(yè)大學(xué)學(xué)報(bào), 2011, 34(2): 129-134.(WU Donghao, WANG Beixin, ZHANG Yong, et al. Advances in the use of biotic index for water quality bioassessment with benthic macroinvertebrate and its perspective in China[J]. Journal of Nanjing Agricultural University, 2011, 34(2): 129-134.(in Chinese))

    [12] KERANS B L, KARR J R. A benthic index of biotic integrity (B-IBI) for rivers of the Tennessee Valley[J].Ecological Application, 1994, 11(4):768-785.

    [13] 楊蓮芳, 李佑文, 戚道光, 等. 九華河水生昆蟲群落結(jié)構(gòu)和水質(zhì)生物評(píng)價(jià)[J]. 生態(tài)學(xué)報(bào), 1992, 12(1): 8-15.(YANG Lianfang, LI Youwen, QI Daoguang, et al. Community structure of aquatic insect and biomonitoring of water quality in Jiuhuahe river[J]. Acta Ecologica Sinica, 1992, 12(1): 8-15.(in Chinese))

    [14] 王備新, 楊蓮芳, 胡本進(jìn), 等. 應(yīng)用底棲動(dòng)物完整性指數(shù)B-IBI評(píng)價(jià)溪流健康[J]. 生態(tài)學(xué)報(bào), 2005, 25(6): 1481-1490.(WANG Beixin, YANG Lianfang, HU Benjin, et al. A preliminary study on the assessment of stream ecosystem health in south of Anhui Province using benthic-index of biotic integrity[J]. Acta Ecologica Sinica, 2005, 25(6): 1481-1490.(in Chinese))

    [15] 渠曉東. 香溪河大型底棲動(dòng)物時(shí)空動(dòng)態(tài)、生物完整性及小水電站的影響研究[D]. 北京: 中國(guó)科學(xué)院研究院, 2006.

    [16] 鄧培雁, 張婉, 王旭濤, 等. 水質(zhì)對(duì)東江流域附生硅藻群落的影響研究[J]. 生態(tài)學(xué)報(bào), 2015, 35(6): 1-13. (DENG Peiyan, ZHANG Wan, WANG Xutao, et al. The effects of water quality on epilithic diatoms communities of Dongjiang river basin[J]. Acta Ecological Sinica, 2015, 35(6): 1-13.(in Chinese))

    [17] 鄧培雁, 雷遠(yuǎn)達(dá), 劉威, 等. 七項(xiàng)河流附著硅藻指數(shù)在東江的適用性評(píng)估[J]. 生態(tài)學(xué)報(bào), 2012, 32(16): 5014-5024. (DENG Peiyan, LEI Yuanda, LIU Wei, et al. Exploration of benthic diatom indices to evaluate water quality in rivers in the Dongjiang Basin[J]. Acta Ecological Sinica, 2012, 32(16): 5014-5024.(in Chinese))

    [18] 孟偉, 張遠(yuǎn), 梁渠東, 等. 河流生態(tài)調(diào)查技術(shù)方法[M]. 北京: 科學(xué)出版社, 2011: 112-119.

    [19] TACHET H, RICHOUX P, BOURNAUD M,等. 淡水無(wú)脊椎動(dòng)物系統(tǒng)分類、生物及生態(tài)學(xué)[M]. 劉威, 王旭濤, 黃少峰, 譯. 北京: 中國(guó)水利水電出版社, 2015.

    [20] 國(guó)家環(huán)境保護(hù)總局, 《水和廢水監(jiān)測(cè)分析方法》編委會(huì). 水和廢水監(jiān)測(cè)分析方法[M]. 4版. 北京: 中國(guó)環(huán)境科學(xué)出版社, 2006.

    [21] LI Li, ZHENG Binghui, LIU Lusan. Biomonitoring and bioindicators used for river ecosystems: definitions, approaches and trends[J]. Procedia Environmental Sciences, 2010,2(1): 1510-1524.

    [22] 渠曉東, 劉志剛, 張遠(yuǎn). 標(biāo)準(zhǔn)化方法篩選參照點(diǎn)構(gòu)建大型底棲動(dòng)物生物完整性指數(shù)[J]. 生態(tài)學(xué)報(bào), 2012, 32(15): 4661-4672.(QU Xiaodong, LIU Zhigang, ZHANG Yuan. Discussion on the standardized method of reference sites selection for establishing the benthic-index of biotic integrity[J]. Acta Ecologica Sinica, 2012, 32(15): 4661-4672.(in Chinese))

    [23] 耿世偉, 渠曉東, 張遠(yuǎn), 等. 大型底棲動(dòng)物生物評(píng)價(jià)指數(shù)比較與應(yīng)用[J]. 環(huán)境科學(xué), 2012, 33(7): 2281-2287.(GENG Shiwei, QU Xiaodong, ZHANG Yuan, et al. Comparison and application of biological indices of Macroinvertebrates in river health assessment[J]. Environmental Science, 2012, 33(7): 2281-2287.(in Chinese))

    [24] 王兆印, 程?hào)|升, 段學(xué)花, 等. 東江河流生態(tài)評(píng)價(jià)及其修復(fù)方略[J]. 水利學(xué)報(bào), 2007, 38(10): 1228-1235.(WANG Zhaoyin, CHENG Dongsheng, DUAN Xuehua, et al. Assessment of ecological system in East River and corresponding ecological restoration strategies[J]. Journal of Hydraulic Engineering, 2007, 38(10): 1228-1235.(in Chinese))

    [25] 王博, 劉全儒, 周云龍, 等. 東江干流底棲動(dòng)物群落結(jié)構(gòu)與水質(zhì)生物學(xué)評(píng)價(jià)[J]. 水生態(tài)學(xué)雜志, 2011, 32(5): 43-49.(WANG Bo, LIU Quanru, ZHOU Yunlong, et al. The community structure of zoobenthos and bioassessment of water quality of the Dongjiang River[J]. Journal of Hydroecology, 2011, 32(5): 43-49. (in Chinese))

    [26] 何琦, 韋桂峰, 劉靜, 等. 增江底棲硅藻種類多樣性及分布特征[J]. 熱帶亞熱帶植物學(xué)報(bào), 2011, 19(3): 245-253.(HE Qi, WEI Guifeng, LIU Jing, et al. Species diversity and spatial distribution of benthic diatom in Zeng River, Guangdong Province, southern China[J]. Journal of Tropical and Subtropical Botany, 2011, 19(3): 245-253.(in Chinese))

    [27] PATRICK R. Use of algae, especially diatoms, in the assessment of water quality[M]. Philadelphia: American Society for Testing and Materials, 1973: 76-95.

    [28] 房英春, 田春, 肖友紅. 蘇子河浮游藻類多樣性與水質(zhì)的生物評(píng)價(jià)[J]. 水利漁業(yè), 2007, 27(1): 57-58.(FANG Yingchun, TIAN Chun, XIAO Youhong. The diversity of planktonic algae and bioassessment of water quality in Suzi River[J]. Reservoir Fisheries, 2007, 27(1): 57-58. (in Chinese))

    [29] 王備新, 楊蓮芳. 我國(guó)東部底棲無(wú)脊椎動(dòng)物主要分類單元耐污值[J]. 生態(tài)學(xué)報(bào), 2004, 24(12): 2768-2775.(WANG Beixin, YANG Lianfang. A study on tolerance values of benthic macroinvertebrate taxa in eastern China[J]. Acta Ecologica Sinica, 2004, 24(12): 2768-2775.(in Chinese))

    [30] 金相燦, 劉鴻亮, 屠清瑛, 等. 中國(guó)湖泊富營(yíng)養(yǎng)化[M]. 北京: 中國(guó)環(huán)境科學(xué)出版社, 1990.

    [31] BOURASSA N, MORIN A. Relationships between size structure of invertebrate assemblages and trophy and substrate composition in streams[J]. Journal of the North American Benthological Society, 1995, 14(3): 393-403.

    [32] 奚旦立, 孫裕生. 環(huán)境監(jiān)測(cè)[M]. 4版.北京: 高等教育出版社, 2010.

    [33] WELCH E B, JACOBY J M, MAY C W. River Ecology and Management[M]. New York: Springer-Verlag, 1998.

    River health assessment based on five biological indices for macroinvertebrates

    SHENG Xiao1, MAO Jianzhong2, CAO Ran1, LI Zhengwu1, WANG Xutao3, DEGN Peiyan1

    (1.SchoolofChemistryandEnvironment,SouthChinaNormalUniversity,Guangzhou510006,China; 2.HydrologyandWaterResourcesBureauofYunnanProvince,Kunming650106,China; 3.PearlRiverWaterEnvironmentMonitoringCenter,Guangzhou510611,China)

    Based on the monitoring data for macroinvertebrates from 24 sampling sites, five indices including the Goodnight-Whitley modified index (GBI), the average score per taxon (ASPT) index, the family biotic index (FBI), the Shannon-Wiener index, and the benthic index of biotic integrity (B-IBI) were used to assess the health status of the main stream and tributaries of the Dongjiang River Basin. The applicability of the indices to the basin was studied. The results show that the overall health status of the Dongjiang River Basin was at a middle level but varied in different parts of the basin: the upper basin had a healthy status, the central part had a moderately heathy status, and the lower part had an unhealthy or extremely unhealthy status. Analysis of the responses of different indices to various types of human activities shows that all the indices, excluding the ASPT index, can reflect the influence of human activities to various degrees. Of the indices, the GBI index was a good indicator of phosphor contamination, the FBI index was a good indicator of nitrite nitrogen contamination, the Shannon-Wiener index could indicate organic pollution well while it had a significantly negative correlation with conductivity, and the B-IBI index had a strongly positive correlation with dissolve oxygen, being a good indicator to ammonia contamination. In conclusion, the comprehensive assessment of the Dongjiang River Baisn with different biological indices can not only reflect the whole condition of the basin reliably, but also indicate the effect of different types of human activities on the river ecosystem of the basin.

    macroinvertebrate; biological index; river health assessment; Dongjiang River

    10.3880/j.issn.1004-6933.2017.01.015

    廣東省水利科技創(chuàng)新項(xiàng)目(2014-16)

    盛蕭(1990—),女,碩士研究生,研究方向?yàn)榄h(huán)境生態(tài)學(xué)。E-mail:winterlilies@163.com

    鄧培雁,教授。E-mail:dpy213@126.com

    Q958.12

    A

    1004-6933(2017)01-0075-08

    2016-04-18 編輯:王 芳)

    猜你喜歡
    東江樣點(diǎn)河流
    小麥條銹病田間為害損失的初步分析
    湖北植保(2022年4期)2022-08-23 10:51:52
    基于空間模擬退火算法的最優(yōu)土壤采樣尺度選擇研究①
    土壤(2021年1期)2021-03-23 07:29:06
    河流
    流放自己的河流
    基于分融策略的土壤采樣設(shè)計(jì)方法*
    萬(wàn)綠水 東江請(qǐng)?jiān)瓌?chuàng)歌曲|全國(guó)征集
    嶺南音樂(2016年4期)2016-11-21 09:25:24
    東江本地早快速投產(chǎn)配套技術(shù)研究
    浙江柑橘(2016年4期)2016-03-11 20:12:58
    當(dāng)河流遇見海
    湖南東江水庫(kù)移民農(nóng)業(yè)事業(yè)單位體制改革探討
    基準(zhǔn)地價(jià)評(píng)估中異常樣點(diǎn)剔除方法研究
    亚洲一区二区三区色噜噜| 女生性感内裤真人,穿戴方法视频| 国产精品一区二区三区四区久久 | 日本撒尿小便嘘嘘汇集6| 少妇的丰满在线观看| svipshipincom国产片| √禁漫天堂资源中文www| 国产高清videossex| 久久久久久久久免费视频了| 免费在线观看视频国产中文字幕亚洲| 一区二区三区高清视频在线| 波多野结衣高清无吗| 男女下面进入的视频免费午夜 | www.精华液| 99精品在免费线老司机午夜| 黄网站色视频无遮挡免费观看| 久久久国产欧美日韩av| 亚洲激情在线av| 亚洲九九香蕉| 制服诱惑二区| www.熟女人妻精品国产| 精品国产国语对白av| 不卡一级毛片| 日本 av在线| 成在线人永久免费视频| 午夜精品在线福利| 精品一区二区三区四区五区乱码| 91老司机精品| 日韩精品中文字幕看吧| 日韩欧美 国产精品| 国产av在哪里看| 欧美午夜高清在线| 欧美性猛交黑人性爽| a级毛片在线看网站| 日本撒尿小便嘘嘘汇集6| 国产aⅴ精品一区二区三区波| 亚洲自偷自拍图片 自拍| 在线观看www视频免费| АⅤ资源中文在线天堂| 久久午夜亚洲精品久久| 日本三级黄在线观看| 久久久久久久久中文| 日韩中文字幕欧美一区二区| 在线观看一区二区三区| 久久香蕉国产精品| tocl精华| 一本精品99久久精品77| 婷婷亚洲欧美| 国产不卡一卡二| videosex国产| 听说在线观看完整版免费高清| 久久性视频一级片| 很黄的视频免费| 精品午夜福利视频在线观看一区| 中文资源天堂在线| 老司机靠b影院| 婷婷六月久久综合丁香| 欧美性猛交黑人性爽| 两性夫妻黄色片| 丰满人妻熟妇乱又伦精品不卡| 女人爽到高潮嗷嗷叫在线视频| 悠悠久久av| 一区二区日韩欧美中文字幕| 亚洲精品中文字幕一二三四区| 真人一进一出gif抽搐免费| 丰满人妻熟妇乱又伦精品不卡| 欧美日韩瑟瑟在线播放| 久久精品人妻少妇| 久久久久久久午夜电影| 亚洲色图 男人天堂 中文字幕| 高清在线国产一区| 一级a爱视频在线免费观看| 午夜影院日韩av| 深夜精品福利| 成年人黄色毛片网站| 99久久综合精品五月天人人| 久久精品国产综合久久久| 国产精品99久久99久久久不卡| 男人舔女人的私密视频| 亚洲无线在线观看| 制服诱惑二区| 可以在线观看毛片的网站| 国产午夜福利久久久久久| 精品人妻1区二区| 天天躁夜夜躁狠狠躁躁| 欧美国产精品va在线观看不卡| 在线观看免费视频日本深夜| 日韩欧美一区视频在线观看| 男女视频在线观看网站免费 | 亚洲在线自拍视频| 一边摸一边做爽爽视频免费| 欧美黑人欧美精品刺激| 99久久无色码亚洲精品果冻| 国产精品永久免费网站| 国产在线精品亚洲第一网站| 高清在线国产一区| 18禁黄网站禁片午夜丰满| 日韩av在线大香蕉| 人成视频在线观看免费观看| 欧美丝袜亚洲另类 | 99久久国产精品久久久| 国产aⅴ精品一区二区三区波| 免费在线观看视频国产中文字幕亚洲| 亚洲成人精品中文字幕电影| 亚洲精品中文字幕一二三四区| 老司机午夜十八禁免费视频| 麻豆一二三区av精品| 色哟哟哟哟哟哟| 久久99热这里只有精品18| 亚洲aⅴ乱码一区二区在线播放 | 色老头精品视频在线观看| 亚洲第一电影网av| 国产伦人伦偷精品视频| 久久精品成人免费网站| 亚洲国产欧美一区二区综合| 国产一区在线观看成人免费| 黄频高清免费视频| 丰满的人妻完整版| 2021天堂中文幕一二区在线观 | 成人三级黄色视频| 亚洲成人精品中文字幕电影| 亚洲国产看品久久| 97碰自拍视频| 在线免费观看的www视频| 成人永久免费在线观看视频| 成人18禁在线播放| 成人亚洲精品av一区二区| 国产成人啪精品午夜网站| 国产精品 欧美亚洲| 国产片内射在线| 啦啦啦韩国在线观看视频| 人妻丰满熟妇av一区二区三区| 国产激情偷乱视频一区二区| 欧美大码av| 国产又黄又爽又无遮挡在线| 欧洲精品卡2卡3卡4卡5卡区| 色尼玛亚洲综合影院| 欧美日韩亚洲国产一区二区在线观看| 亚洲第一欧美日韩一区二区三区| 欧美激情久久久久久爽电影| 美女午夜性视频免费| 亚洲七黄色美女视频| 婷婷丁香在线五月| 亚洲午夜精品一区,二区,三区| 一a级毛片在线观看| 日本免费a在线| 久久精品国产99精品国产亚洲性色| 香蕉久久夜色| 男女做爰动态图高潮gif福利片| 男女视频在线观看网站免费 | 高潮久久久久久久久久久不卡| 黄网站色视频无遮挡免费观看| 中文字幕人妻丝袜一区二区| 人人妻人人澡人人看| 精品国产乱子伦一区二区三区| 国产精品亚洲一级av第二区| 琪琪午夜伦伦电影理论片6080| 一进一出抽搐动态| 成人永久免费在线观看视频| 亚洲激情在线av| 一个人观看的视频www高清免费观看 | 很黄的视频免费| 热re99久久国产66热| 黑人巨大精品欧美一区二区mp4| 两个人免费观看高清视频| 国产午夜精品久久久久久| 久久天躁狠狠躁夜夜2o2o| 一级作爱视频免费观看| 中文在线观看免费www的网站 | 在线观看www视频免费| 男女下面进入的视频免费午夜 | xxx96com| 久久久久久人人人人人| 亚洲av美国av| 黄片大片在线免费观看| 欧美国产日韩亚洲一区| 亚洲一区二区三区不卡视频| 久久久国产欧美日韩av| 久久中文字幕人妻熟女| 午夜老司机福利片| 日韩大尺度精品在线看网址| 在线观看舔阴道视频| 一边摸一边做爽爽视频免费| 亚洲精品粉嫩美女一区| 麻豆国产av国片精品| 91大片在线观看| 成人特级黄色片久久久久久久| 久久天堂一区二区三区四区| 久久性视频一级片| 在线国产一区二区在线| 午夜福利欧美成人| 黄片大片在线免费观看| 一级毛片女人18水好多| 婷婷精品国产亚洲av在线| 露出奶头的视频| aaaaa片日本免费| 搡老熟女国产l中国老女人| 久久国产乱子伦精品免费另类| 亚洲av第一区精品v没综合| 精品久久蜜臀av无| 又紧又爽又黄一区二区| 日韩欧美一区视频在线观看| www.www免费av| 一级毛片精品| 无遮挡黄片免费观看| 欧美zozozo另类| 精品久久蜜臀av无| 国产精品免费视频内射| 女性生殖器流出的白浆| 一进一出抽搐动态| 无限看片的www在线观看| 午夜福利一区二区在线看| 亚洲aⅴ乱码一区二区在线播放 | 少妇被粗大的猛进出69影院| 男人舔女人下体高潮全视频| 中文字幕人成人乱码亚洲影| 日本 av在线| 亚洲精品色激情综合| 精品人妻1区二区| 欧美精品亚洲一区二区| 国产爱豆传媒在线观看 | 国产主播在线观看一区二区| 日韩成人在线观看一区二区三区| 日韩欧美 国产精品| 亚洲国产精品成人综合色| 成在线人永久免费视频| 色综合亚洲欧美另类图片| 成年版毛片免费区| 免费观看精品视频网站| 两个人看的免费小视频| 两个人看的免费小视频| 亚洲精品一卡2卡三卡4卡5卡| 18禁观看日本| av在线播放免费不卡| 国产精品香港三级国产av潘金莲| 波多野结衣av一区二区av| 欧美av亚洲av综合av国产av| 欧美在线黄色| 久久精品国产亚洲av香蕉五月| 欧美日本亚洲视频在线播放| 一夜夜www| 中文字幕精品亚洲无线码一区 | 欧美在线黄色| 麻豆av在线久日| 特大巨黑吊av在线直播 | 亚洲精品久久国产高清桃花| 亚洲av中文字字幕乱码综合 | 中亚洲国语对白在线视频| 国产精品永久免费网站| 麻豆成人av在线观看| 久久99热这里只有精品18| 日本一区二区免费在线视频| 欧美乱码精品一区二区三区| 亚洲国产精品sss在线观看| cao死你这个sao货| 日韩大尺度精品在线看网址| 午夜福利18| 久久亚洲真实| 精品国内亚洲2022精品成人| 国产av一区二区精品久久| 一级片免费观看大全| 可以免费在线观看a视频的电影网站| 亚洲精品一卡2卡三卡4卡5卡| 香蕉av资源在线| 欧美成人一区二区免费高清观看 | 大香蕉久久成人网| 精品国产国语对白av| 午夜免费观看网址| av超薄肉色丝袜交足视频| 99国产精品99久久久久| 国产av在哪里看| 这个男人来自地球电影免费观看| 欧美乱色亚洲激情| 免费av毛片视频| 久久天躁狠狠躁夜夜2o2o| 亚洲av成人一区二区三| 精品国产美女av久久久久小说| 亚洲人成电影免费在线| 亚洲国产日韩欧美精品在线观看 | 黄色视频不卡| 真人一进一出gif抽搐免费| 亚洲国产日韩欧美精品在线观看 | 亚洲人成伊人成综合网2020| 成人国产综合亚洲| 岛国视频午夜一区免费看| 桃色一区二区三区在线观看| 欧美激情久久久久久爽电影| 久久国产乱子伦精品免费另类| 亚洲国产欧美一区二区综合| 午夜日韩欧美国产| 免费av毛片视频| 欧美成狂野欧美在线观看| 精品高清国产在线一区| 久久精品国产亚洲av香蕉五月| 禁无遮挡网站| 又大又爽又粗| 亚洲国产欧美网| 91麻豆精品激情在线观看国产| 亚洲成av人片免费观看| 亚洲欧美日韩高清在线视频| 在线观看一区二区三区| 嫁个100分男人电影在线观看| 成人18禁高潮啪啪吃奶动态图| 日韩欧美国产一区二区入口| 精品免费久久久久久久清纯| 国产亚洲欧美精品永久| 国产亚洲av高清不卡| 香蕉av资源在线| www国产在线视频色| 女同久久另类99精品国产91| 丁香欧美五月| 一本一本综合久久| 亚洲精品av麻豆狂野| 不卡一级毛片| 女性生殖器流出的白浆| 精品欧美一区二区三区在线| av中文乱码字幕在线| 中亚洲国语对白在线视频| 久久午夜综合久久蜜桃| 丁香欧美五月| 一区二区三区精品91| 美女国产高潮福利片在线看| АⅤ资源中文在线天堂| 中出人妻视频一区二区| 亚洲成人久久性| 极品教师在线免费播放| 欧美日韩亚洲国产一区二区在线观看| 久久婷婷成人综合色麻豆| 亚洲av电影在线进入| 欧美国产日韩亚洲一区| 精品久久久久久,| 满18在线观看网站| 免费观看精品视频网站| 亚洲精品一区av在线观看| 一级作爱视频免费观看| 欧美国产精品va在线观看不卡| 国产精品爽爽va在线观看网站 | 亚洲天堂国产精品一区在线| av天堂在线播放| 免费搜索国产男女视频| 国产又黄又爽又无遮挡在线| 国产欧美日韩一区二区三| 午夜激情福利司机影院| 精品一区二区三区视频在线观看免费| 亚洲人成77777在线视频| 精品日产1卡2卡| 精品电影一区二区在线| 亚洲七黄色美女视频| 在线观看舔阴道视频| 国产91精品成人一区二区三区| 俺也久久电影网| 欧美日韩精品网址| or卡值多少钱| 亚洲在线自拍视频| 亚洲片人在线观看| 狂野欧美激情性xxxx| 亚洲人成网站在线播放欧美日韩| 国产午夜精品久久久久久| 亚洲天堂国产精品一区在线| 久久久国产欧美日韩av| 亚洲国产欧美一区二区综合| 亚洲专区字幕在线| 亚洲久久久国产精品| 91麻豆av在线| 国产成人欧美| 叶爱在线成人免费视频播放| 欧美最黄视频在线播放免费| 日韩欧美一区视频在线观看| 婷婷丁香在线五月| 国产一区二区三区视频了| 欧美激情 高清一区二区三区| 久久久久国产一级毛片高清牌| 99久久久亚洲精品蜜臀av| 高清毛片免费观看视频网站| 亚洲一区二区三区不卡视频| 亚洲,欧美精品.| 人妻久久中文字幕网| 免费在线观看影片大全网站| 99在线视频只有这里精品首页| 色老头精品视频在线观看| 听说在线观看完整版免费高清| 久久国产精品人妻蜜桃| 午夜福利在线在线| 日本一区二区免费在线视频| 国产精品99久久99久久久不卡| 一夜夜www| 久久青草综合色| 欧美在线黄色| 99热6这里只有精品| 婷婷亚洲欧美| а√天堂www在线а√下载| 国产成人av激情在线播放| 国产精品久久久av美女十八| 好男人在线观看高清免费视频 | 亚洲九九香蕉| 成人国语在线视频| 亚洲欧美日韩无卡精品| 黄色片一级片一级黄色片| 18禁观看日本| 亚洲成av人片免费观看| 亚洲中文字幕一区二区三区有码在线看 | 黄色成人免费大全| 亚洲一码二码三码区别大吗| 久久久久国产精品人妻aⅴ院| 国产高清有码在线观看视频 | 中文亚洲av片在线观看爽| 在线观看午夜福利视频| 波多野结衣巨乳人妻| 久久久久久国产a免费观看| 成人一区二区视频在线观看| 日韩欧美一区二区三区在线观看| 国内精品久久久久久久电影| 男人舔奶头视频| 午夜免费激情av| 美国免费a级毛片| 欧美日韩亚洲综合一区二区三区_| 九色国产91popny在线| 国产精品免费一区二区三区在线| 国产精品 欧美亚洲| 亚洲精品在线观看二区| 非洲黑人性xxxx精品又粗又长| 亚洲精品在线美女| 亚洲成av人片免费观看| 亚洲av五月六月丁香网| 亚洲久久久国产精品| 国产激情久久老熟女| 身体一侧抽搐| 真人一进一出gif抽搐免费| 手机成人av网站| 男男h啪啪无遮挡| 日韩欧美一区二区三区在线观看| 成人18禁高潮啪啪吃奶动态图| av欧美777| 熟女少妇亚洲综合色aaa.| 久久精品成人免费网站| 婷婷亚洲欧美| 99国产综合亚洲精品| 韩国av一区二区三区四区| 亚洲av美国av| 欧美大码av| 一a级毛片在线观看| 99久久精品国产亚洲精品| 亚洲精品国产区一区二| 久久中文字幕一级| 岛国在线观看网站| 欧美成人性av电影在线观看| 一二三四在线观看免费中文在| 两个人视频免费观看高清| netflix在线观看网站| 视频区欧美日本亚洲| 国产片内射在线| 亚洲国产毛片av蜜桃av| 校园春色视频在线观看| av中文乱码字幕在线| 国产成人欧美在线观看| 2021天堂中文幕一二区在线观 | 一区二区三区激情视频| 国产成人影院久久av| 欧美日韩瑟瑟在线播放| 母亲3免费完整高清在线观看| av免费在线观看网站| 国产一区二区在线av高清观看| 亚洲在线自拍视频| 黑人巨大精品欧美一区二区mp4| 午夜福利在线观看吧| 午夜日韩欧美国产| 亚洲精品中文字幕一二三四区| 国产亚洲精品综合一区在线观看 | 久久久久久大精品| 99精品在免费线老司机午夜| 久久久久久大精品| 看片在线看免费视频| 久久中文字幕一级| 午夜日韩欧美国产| av视频在线观看入口| 国产伦在线观看视频一区| 国产黄色小视频在线观看| 午夜日韩欧美国产| 日本a在线网址| 国产伦一二天堂av在线观看| or卡值多少钱| 亚洲最大成人中文| 俄罗斯特黄特色一大片| 少妇熟女aⅴ在线视频| 久久久久久国产a免费观看| 99国产精品99久久久久| 国产精品久久久av美女十八| 日韩一卡2卡3卡4卡2021年| 岛国在线观看网站| 在线观看午夜福利视频| 欧美国产精品va在线观看不卡| cao死你这个sao货| 天天一区二区日本电影三级| 欧美黑人巨大hd| 两性午夜刺激爽爽歪歪视频在线观看 | 宅男免费午夜| 天堂影院成人在线观看| 久热这里只有精品99| 国产成人av激情在线播放| 人人妻人人澡欧美一区二区| 可以免费在线观看a视频的电影网站| 亚洲国产看品久久| 午夜久久久在线观看| 国产私拍福利视频在线观看| 国产高清videossex| 巨乳人妻的诱惑在线观看| 中文字幕最新亚洲高清| 亚洲av成人不卡在线观看播放网| 免费看日本二区| 午夜福利一区二区在线看| 亚洲第一欧美日韩一区二区三区| 国产蜜桃级精品一区二区三区| 夜夜躁狠狠躁天天躁| 国产精品二区激情视频| 91九色精品人成在线观看| 两个人视频免费观看高清| 中文亚洲av片在线观看爽| 99在线视频只有这里精品首页| 国内少妇人妻偷人精品xxx网站 | 99国产综合亚洲精品| 精品久久久久久久久久免费视频| 免费高清视频大片| svipshipincom国产片| 他把我摸到了高潮在线观看| 九色国产91popny在线| 午夜福利免费观看在线| 精品卡一卡二卡四卡免费| 久久精品国产亚洲av高清一级| 亚洲av成人一区二区三| 久久亚洲真实| 亚洲午夜精品一区,二区,三区| 日韩高清综合在线| 麻豆成人av在线观看| 在线观看免费午夜福利视频| 12—13女人毛片做爰片一| 欧美中文日本在线观看视频| 在线视频色国产色| 亚洲国产精品久久男人天堂| 日韩欧美一区视频在线观看| 中文字幕av电影在线播放| 亚洲欧美一区二区三区黑人| 国产伦一二天堂av在线观看| 免费人成视频x8x8入口观看| 国产精品国产高清国产av| 亚洲精品久久国产高清桃花| 在线十欧美十亚洲十日本专区| 日韩高清综合在线| 女性被躁到高潮视频| a在线观看视频网站| 亚洲熟女毛片儿| 男人的好看免费观看在线视频 | 手机成人av网站| 欧美zozozo另类| 欧美三级亚洲精品| 中文字幕人妻熟女乱码| 一级黄色大片毛片| 一级作爱视频免费观看| 999精品在线视频| 国产精品美女特级片免费视频播放器 | 可以免费在线观看a视频的电影网站| 手机成人av网站| 国产亚洲av嫩草精品影院| 热re99久久国产66热| 搡老妇女老女人老熟妇| 熟妇人妻久久中文字幕3abv| 国产精品久久视频播放| 国产视频一区二区在线看| 后天国语完整版免费观看| 好男人电影高清在线观看| 成人三级做爰电影| 天天躁狠狠躁夜夜躁狠狠躁| 国产亚洲欧美精品永久| 亚洲精品粉嫩美女一区| 欧美乱码精品一区二区三区| 午夜免费激情av| 精品不卡国产一区二区三区| 非洲黑人性xxxx精品又粗又长| 国产亚洲精品久久久久久毛片| 美女 人体艺术 gogo| 99久久久亚洲精品蜜臀av| 免费看日本二区| 99re在线观看精品视频| 久久久久国产一级毛片高清牌| 伊人久久大香线蕉亚洲五| 精品无人区乱码1区二区| 国产成人系列免费观看| 国产单亲对白刺激| 欧美午夜高清在线| 久久热在线av| 制服丝袜大香蕉在线| 亚洲最大成人中文| 亚洲色图 男人天堂 中文字幕| 每晚都被弄得嗷嗷叫到高潮| 91麻豆av在线| 久久精品aⅴ一区二区三区四区| 可以在线观看毛片的网站| 国产伦人伦偷精品视频| 亚洲五月色婷婷综合| 夜夜夜夜夜久久久久| 成人三级黄色视频| 啦啦啦韩国在线观看视频| 免费高清在线观看日韩| 色老头精品视频在线观看| 久久精品成人免费网站| 欧美日本亚洲视频在线播放| 狂野欧美激情性xxxx| 国产精品二区激情视频| 精品欧美一区二区三区在线| www日本黄色视频网| 亚洲免费av在线视频| 十八禁人妻一区二区| 国产激情欧美一区二区| 99re在线观看精品视频| 欧美乱妇无乱码| 在线观看免费视频日本深夜|