• <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)剔除方法研究
    久久精品国产清高在天天线| 亚洲真实伦在线观看| 亚洲天堂国产精品一区在线| 成年女人毛片免费观看观看9| 国产成人欧美在线观看| e午夜精品久久久久久久| 老司机福利观看| 亚洲avbb在线观看| 成人三级做爰电影| av黄色大香蕉| 啦啦啦观看免费观看视频高清| 高潮久久久久久久久久久不卡| 亚洲在线自拍视频| 99热精品在线国产| 老司机深夜福利视频在线观看| av天堂中文字幕网| 国内揄拍国产精品人妻在线| 两个人视频免费观看高清| 97超级碰碰碰精品色视频在线观看| 可以在线观看毛片的网站| 不卡一级毛片| 深夜精品福利| 毛片女人毛片| 黑人欧美特级aaaaaa片| 亚洲国产高清在线一区二区三| 黄色女人牲交| 国产成人精品久久二区二区91| 国产69精品久久久久777片 | 少妇裸体淫交视频免费看高清| 午夜激情福利司机影院| 一区二区三区国产精品乱码| 18禁观看日本| 成人特级av手机在线观看| 亚洲欧美日韩高清专用| 亚洲欧美日韩高清在线视频| 国产亚洲av高清不卡| 热99在线观看视频| 亚洲五月天丁香| 亚洲av电影不卡..在线观看| 成人高潮视频无遮挡免费网站| 亚洲乱码一区二区免费版| 国产一区二区三区在线臀色熟女| 中文字幕精品亚洲无线码一区| 男女那种视频在线观看| 动漫黄色视频在线观看| 免费高清视频大片| 99在线视频只有这里精品首页| 亚洲av成人精品一区久久| 国产亚洲精品久久久久久毛片| 国产伦精品一区二区三区视频9 | 一本一本综合久久| 一级毛片女人18水好多| 久久久久久人人人人人| 欧美性猛交黑人性爽| 三级毛片av免费| 亚洲中文日韩欧美视频| 午夜免费成人在线视频| 叶爱在线成人免费视频播放| 久久精品国产99精品国产亚洲性色| 久久久久免费精品人妻一区二区| 欧美日本亚洲视频在线播放| 亚洲九九香蕉| 国产三级中文精品| 亚洲精品国产精品久久久不卡| 成人av在线播放网站| 99精品在免费线老司机午夜| 国产精品久久久久久久电影 | 国模一区二区三区四区视频 | 国产成人福利小说| 国产精品美女特级片免费视频播放器 | 亚洲在线观看片| 精品一区二区三区视频在线观看免费| 国产欧美日韩精品亚洲av| 1024手机看黄色片| 亚洲av中文字字幕乱码综合| 男女午夜视频在线观看| 男女之事视频高清在线观看| 成年女人毛片免费观看观看9| 一本综合久久免费| 色在线成人网| 亚洲七黄色美女视频| 久久草成人影院| 超碰成人久久| 黄色女人牲交| 午夜精品在线福利| 中文字幕av在线有码专区| 亚洲精品在线观看二区| 91av网站免费观看| 成熟少妇高潮喷水视频| 国语自产精品视频在线第100页| 亚洲中文字幕一区二区三区有码在线看 | 在线观看日韩欧美| 亚洲成av人片在线播放无| 国产av在哪里看| 日本免费a在线| 99久国产av精品| 日韩精品中文字幕看吧| 国产91精品成人一区二区三区| 日韩高清综合在线| 国产激情欧美一区二区| 精品人妻1区二区| 午夜亚洲福利在线播放| www.精华液| 九九在线视频观看精品| 亚洲精品粉嫩美女一区| 99久久精品国产亚洲精品| 熟妇人妻久久中文字幕3abv| 男人舔女人下体高潮全视频| 日本一二三区视频观看| 免费一级毛片在线播放高清视频| 2021天堂中文幕一二区在线观| 免费观看人在逋| 99热这里只有精品一区 | 亚洲精品中文字幕一二三四区| 91麻豆av在线| 日本成人三级电影网站| 国产成人福利小说| 99久久精品国产亚洲精品| 国产伦精品一区二区三区四那| 老司机在亚洲福利影院| 精品国产超薄肉色丝袜足j| av天堂在线播放| 免费av毛片视频| 99在线人妻在线中文字幕| 国产精品av久久久久免费| 真人一进一出gif抽搐免费| 看片在线看免费视频| 免费看a级黄色片| 亚洲七黄色美女视频| 高清在线国产一区| 日韩免费av在线播放| 黄色视频,在线免费观看| 国产v大片淫在线免费观看| 亚洲av五月六月丁香网| 亚洲国产欧洲综合997久久,| 亚洲欧美激情综合另类| 婷婷精品国产亚洲av| 欧美黄色片欧美黄色片| 在线看三级毛片| 色老头精品视频在线观看| 国产91精品成人一区二区三区| 久久久国产成人免费| 99久久99久久久精品蜜桃| 成人三级黄色视频| 日本一二三区视频观看| 91九色精品人成在线观看| 日韩 欧美 亚洲 中文字幕| 欧美成人一区二区免费高清观看 | 久久中文看片网| 91av网站免费观看| 日韩有码中文字幕| 久久伊人香网站| 欧美绝顶高潮抽搐喷水| 国产一级毛片七仙女欲春2| 日本黄色视频三级网站网址| 成人国产一区最新在线观看| 丰满人妻一区二区三区视频av | 午夜a级毛片| 国产精品精品国产色婷婷| 国产精品 国内视频| 毛片女人毛片| 欧美另类亚洲清纯唯美| 一级毛片精品| 成人一区二区视频在线观看| 757午夜福利合集在线观看| 好男人电影高清在线观看| 久久婷婷人人爽人人干人人爱| 国内毛片毛片毛片毛片毛片| 国产黄色小视频在线观看| 麻豆成人午夜福利视频| 欧美丝袜亚洲另类 | 亚洲成a人片在线一区二区| 国产高清videossex| 亚洲一区二区三区不卡视频| 每晚都被弄得嗷嗷叫到高潮| 国产精品久久久人人做人人爽| 亚洲欧美日韩卡通动漫| av片东京热男人的天堂| 波多野结衣巨乳人妻| 欧美日韩黄片免| 18禁国产床啪视频网站| 最近最新中文字幕大全免费视频| 99久久精品国产亚洲精品| 一本精品99久久精品77| 无限看片的www在线观看| 成人三级黄色视频| 嫩草影院精品99| 美女 人体艺术 gogo| 免费看光身美女| 18禁国产床啪视频网站| 性欧美人与动物交配| 女同久久另类99精品国产91| 成年女人永久免费观看视频| 国产精品永久免费网站| av中文乱码字幕在线| 国产精品 国内视频| 国产精品电影一区二区三区| 国内揄拍国产精品人妻在线| 99精品久久久久人妻精品| 亚洲va日本ⅴa欧美va伊人久久| 一进一出好大好爽视频| e午夜精品久久久久久久| 99视频精品全部免费 在线 | 女同久久另类99精品国产91| 免费大片18禁| 亚洲欧美日韩高清在线视频| 国产精品香港三级国产av潘金莲| 热99re8久久精品国产| 成人高潮视频无遮挡免费网站| 亚洲18禁久久av| 狂野欧美激情性xxxx| 欧美三级亚洲精品| 精品免费久久久久久久清纯| 日日摸夜夜添夜夜添小说| 在线观看一区二区三区| 三级毛片av免费| 国产精品99久久99久久久不卡| 久久久久久久午夜电影| 亚洲第一电影网av| 日韩欧美精品v在线| 两性夫妻黄色片| 免费在线观看视频国产中文字幕亚洲| 亚洲精品粉嫩美女一区| 亚洲欧洲精品一区二区精品久久久| 嫩草影视91久久| 欧美黑人巨大hd| 成人精品一区二区免费| 一进一出好大好爽视频| 久久久久国产一级毛片高清牌| 人妻夜夜爽99麻豆av| av视频在线观看入口| 亚洲 欧美一区二区三区| 国内揄拍国产精品人妻在线| 夜夜夜夜夜久久久久| 一区二区三区激情视频| 69av精品久久久久久| 国产激情欧美一区二区| 国产91精品成人一区二区三区| 欧美丝袜亚洲另类 | x7x7x7水蜜桃| 日韩中文字幕欧美一区二区| 欧美zozozo另类| 国产欧美日韩精品一区二区| 国产精品久久久久久精品电影| 97超级碰碰碰精品色视频在线观看| 成人特级黄色片久久久久久久| 一个人免费在线观看电影 | 免费大片18禁| 亚洲av中文字字幕乱码综合| 亚洲精品色激情综合| 欧美一级毛片孕妇| 婷婷精品国产亚洲av在线| www.www免费av| 特级一级黄色大片| 综合色av麻豆| 中文字幕久久专区| 香蕉av资源在线| 国产单亲对白刺激| 97碰自拍视频| 男人和女人高潮做爰伦理| 岛国在线观看网站| 国产精品av久久久久免费| 日韩大尺度精品在线看网址| 国产精品久久久人人做人人爽| 欧美黄色淫秽网站| 此物有八面人人有两片| 禁无遮挡网站| 国产精品98久久久久久宅男小说| 男女之事视频高清在线观看| 韩国av一区二区三区四区| 五月玫瑰六月丁香| 国产综合懂色| 男女做爰动态图高潮gif福利片| 色视频www国产| 身体一侧抽搐| 美女cb高潮喷水在线观看 | 一本一本综合久久| 亚洲人成网站在线播放欧美日韩| 男女视频在线观看网站免费| 亚洲最大成人中文| 91av网站免费观看| 在线观看舔阴道视频| www.999成人在线观看| 蜜桃久久精品国产亚洲av| 免费看a级黄色片| 亚洲国产色片| 亚洲精品在线美女| 午夜免费激情av| 国产毛片a区久久久久| 亚洲av电影在线进入| 国产精品一区二区免费欧美| 精品国产美女av久久久久小说| 国产单亲对白刺激| 国产野战对白在线观看| 熟女人妻精品中文字幕| 国产高潮美女av| 十八禁人妻一区二区| 亚洲国产精品合色在线| 亚洲欧美一区二区三区黑人| 香蕉久久夜色| av在线天堂中文字幕| 精品国内亚洲2022精品成人| 十八禁网站免费在线| 级片在线观看| 国产伦精品一区二区三区视频9 | 午夜精品一区二区三区免费看| 99久久精品热视频| svipshipincom国产片| 欧美av亚洲av综合av国产av| 夜夜夜夜夜久久久久| 日本免费一区二区三区高清不卡| 麻豆国产av国片精品| 国产精品亚洲一级av第二区| 在线观看午夜福利视频| 网址你懂的国产日韩在线| 久久中文字幕一级| 搡老熟女国产l中国老女人| 国产亚洲欧美在线一区二区| 99热精品在线国产| 日韩三级视频一区二区三区| 免费在线观看日本一区| 亚洲国产欧洲综合997久久,| 日韩欧美国产在线观看| 91麻豆精品激情在线观看国产| 精品99又大又爽又粗少妇毛片 | 日本在线视频免费播放| 热99在线观看视频| 国产麻豆成人av免费视频| 久久久久国产精品人妻aⅴ院| 亚洲,欧美精品.| 国产视频内射| 韩国av一区二区三区四区| 九九久久精品国产亚洲av麻豆 | 午夜福利欧美成人| 国产成人aa在线观看| 精品福利观看| 噜噜噜噜噜久久久久久91| 真实男女啪啪啪动态图| 欧美在线黄色| 亚洲人与动物交配视频| 又黄又粗又硬又大视频| 欧美午夜高清在线| 欧美日韩福利视频一区二区| 变态另类成人亚洲欧美熟女| 国产精品亚洲美女久久久| 免费高清视频大片| 级片在线观看| 国产一区二区在线观看日韩 | 1024手机看黄色片| 国内精品一区二区在线观看| 淫秽高清视频在线观看| 真实男女啪啪啪动态图| avwww免费| 欧美3d第一页| 国产高清视频在线播放一区| 麻豆av在线久日| 国内精品一区二区在线观看| 中文字幕熟女人妻在线| 人妻丰满熟妇av一区二区三区| 国产主播在线观看一区二区| 一进一出抽搐动态| 变态另类成人亚洲欧美熟女| 黄色视频,在线免费观看| 亚洲av成人精品一区久久| 男女做爰动态图高潮gif福利片| 日本黄色片子视频| 两个人的视频大全免费| 又大又爽又粗| 曰老女人黄片| 亚洲av电影在线进入| 国产又黄又爽又无遮挡在线| 久久亚洲精品不卡| 国产精华一区二区三区| 两个人看的免费小视频| 亚洲 国产 在线| 视频区欧美日本亚洲| 三级毛片av免费| 午夜福利在线观看免费完整高清在 | 精品国产三级普通话版| 亚洲成人久久爱视频| 国产黄片美女视频| 天堂√8在线中文| 日韩有码中文字幕| 欧美激情久久久久久爽电影| 国产一区二区在线av高清观看| 日本三级黄在线观看| 成人国产一区最新在线观看| 又大又爽又粗| 欧美绝顶高潮抽搐喷水| 国产午夜精品久久久久久| 久久伊人香网站| 最好的美女福利视频网| 亚洲欧美日韩高清专用| 老熟妇乱子伦视频在线观看| 亚洲狠狠婷婷综合久久图片| 男女那种视频在线观看| 啦啦啦免费观看视频1| 国产精品美女特级片免费视频播放器 | av欧美777| 国产97色在线日韩免费| 动漫黄色视频在线观看| 男女视频在线观看网站免费| 欧美日韩国产亚洲二区| www.999成人在线观看| 精品国内亚洲2022精品成人| 露出奶头的视频| 亚洲午夜理论影院| 精品免费久久久久久久清纯| 亚洲欧美精品综合一区二区三区| 久久精品91无色码中文字幕| 不卡一级毛片| 中文字幕精品亚洲无线码一区| 日韩欧美 国产精品| 在线观看午夜福利视频| 色老头精品视频在线观看| 亚洲熟妇熟女久久| 真实男女啪啪啪动态图| 1024香蕉在线观看| 99久久综合精品五月天人人| 国产男靠女视频免费网站| 一本一本综合久久| 在线a可以看的网站| 免费搜索国产男女视频| 天天一区二区日本电影三级| 俄罗斯特黄特色一大片| 国产亚洲精品综合一区在线观看| 黄片小视频在线播放| 熟女人妻精品中文字幕| 欧美日韩精品网址| 精品国产乱码久久久久久男人| 丁香六月欧美| 深夜精品福利| 国产av一区在线观看免费| 日韩成人在线观看一区二区三区| 久99久视频精品免费| 国产精品亚洲一级av第二区| 岛国在线免费视频观看| netflix在线观看网站| 亚洲熟妇中文字幕五十中出| 精品久久久久久久久久久久久| 午夜福利欧美成人| 美女高潮的动态| 日韩欧美国产一区二区入口| 国产欧美日韩一区二区精品| 成人特级黄色片久久久久久久| 久久天躁狠狠躁夜夜2o2o| 麻豆久久精品国产亚洲av| www.www免费av| 人妻夜夜爽99麻豆av| 精品无人区乱码1区二区| 欧美av亚洲av综合av国产av| 亚洲第一电影网av| 19禁男女啪啪无遮挡网站| 欧美乱码精品一区二区三区| 久久亚洲真实| 1000部很黄的大片| 日本撒尿小便嘘嘘汇集6| 丁香欧美五月| 国产又黄又爽又无遮挡在线| 亚洲av日韩精品久久久久久密| 国产成+人综合+亚洲专区| 母亲3免费完整高清在线观看| 一区二区三区激情视频| 精品一区二区三区四区五区乱码| 欧美在线黄色| 亚洲在线自拍视频| 99国产综合亚洲精品| 手机成人av网站| 桃红色精品国产亚洲av| 神马国产精品三级电影在线观看| 欧美日韩亚洲国产一区二区在线观看| 亚洲av美国av| 99久久99久久久精品蜜桃| 欧美zozozo另类| 欧美三级亚洲精品| 国产视频内射| 免费av毛片视频| 国产欧美日韩精品一区二区| а√天堂www在线а√下载| 三级男女做爰猛烈吃奶摸视频| 欧美成狂野欧美在线观看| 97碰自拍视频| 日韩欧美在线二视频| 亚洲精品色激情综合| 亚洲av片天天在线观看| 夜夜爽天天搞| ponron亚洲| 99久国产av精品| 精品乱码久久久久久99久播| 黑人操中国人逼视频| 男女视频在线观看网站免费| 午夜视频精品福利| 精品午夜福利视频在线观看一区| 美女高潮喷水抽搐中文字幕| 亚洲一区二区三区不卡视频| 18禁裸乳无遮挡免费网站照片| 国产乱人视频| 美女扒开内裤让男人捅视频| 长腿黑丝高跟| 特级一级黄色大片| 美女cb高潮喷水在线观看 | 欧美另类亚洲清纯唯美| 99在线人妻在线中文字幕| 国产伦一二天堂av在线观看| 高清毛片免费观看视频网站| 国产av在哪里看| 亚洲第一电影网av| 九九热线精品视视频播放| 老熟妇乱子伦视频在线观看| 成人av在线播放网站| 一级a爱片免费观看的视频| 日本一本二区三区精品| 免费看美女性在线毛片视频| 国产精品亚洲一级av第二区| 国产日本99.免费观看| 草草在线视频免费看| 亚洲 欧美一区二区三区| 精品免费久久久久久久清纯| 不卡一级毛片| 99久久综合精品五月天人人| 亚洲熟妇中文字幕五十中出| 老司机深夜福利视频在线观看| 午夜免费激情av| 最好的美女福利视频网| 国产aⅴ精品一区二区三区波| 国产一区二区三区在线臀色熟女| 麻豆成人av在线观看| 精品一区二区三区四区五区乱码| 在线永久观看黄色视频| 看黄色毛片网站| 一夜夜www| 亚洲电影在线观看av| 亚洲一区高清亚洲精品| 亚洲,欧美精品.| 欧美日韩国产亚洲二区| 国产亚洲精品久久久com| 久久国产乱子伦精品免费另类| 久久热在线av| 性色av乱码一区二区三区2| 亚洲成人免费电影在线观看| xxxwww97欧美| 欧美黄色片欧美黄色片| 老司机深夜福利视频在线观看| 给我免费播放毛片高清在线观看| 欧美av亚洲av综合av国产av| 亚洲欧美日韩高清在线视频| 国产高清有码在线观看视频| 国产高潮美女av| 欧美另类亚洲清纯唯美| 国产精品日韩av在线免费观看| h日本视频在线播放| 美女 人体艺术 gogo| 他把我摸到了高潮在线观看| 亚洲无线观看免费| 欧美乱妇无乱码| 在线视频色国产色| 99热精品在线国产| 男女视频在线观看网站免费| 日韩精品青青久久久久久| 精品久久久久久久久久久久久| 成在线人永久免费视频| 久久精品91蜜桃| 久久久久国产一级毛片高清牌| 久久精品影院6| 亚洲av五月六月丁香网| 三级毛片av免费| 亚洲激情在线av| 欧美乱色亚洲激情| 一本一本综合久久| 在线观看一区二区三区| 久久久久久久久中文| www.www免费av| 精品久久久久久,| 亚洲人成网站高清观看| 欧美成人性av电影在线观看| 大型黄色视频在线免费观看| 国产av一区在线观看免费| 俄罗斯特黄特色一大片| 国内精品美女久久久久久| 欧美极品一区二区三区四区| 亚洲国产欧美人成| 女同久久另类99精品国产91| 欧美中文日本在线观看视频| 国产精品一区二区三区四区久久| 亚洲最大成人中文| 国产精品美女特级片免费视频播放器 | 中文字幕精品亚洲无线码一区| 国产91精品成人一区二区三区| 亚洲自偷自拍图片 自拍| 男人和女人高潮做爰伦理| 午夜福利欧美成人| 亚洲av片天天在线观看| 久久性视频一级片| 欧美在线黄色| 日韩欧美在线乱码| 久久九九热精品免费| www日本黄色视频网| 丁香六月欧美| a在线观看视频网站| 18禁黄网站禁片免费观看直播| 一级毛片精品| 这个男人来自地球电影免费观看| 少妇人妻一区二区三区视频| 国产欧美日韩精品一区二区| 又大又爽又粗| 国产一区在线观看成人免费| 少妇的逼水好多| 国产亚洲av嫩草精品影院| 最新中文字幕久久久久 | 99在线人妻在线中文字幕| 一二三四在线观看免费中文在| 女人被狂操c到高潮| 国产三级黄色录像| 夜夜看夜夜爽夜夜摸| 亚洲avbb在线观看|